Analysis of embryonic development in the unsequenced axolotl: Waves of transcriptomic upheaval and stability

Pre-Embryonic Period Observation Shows a Unique Reproductive Strategy of the Critically Endangered Anji Salamander (Hynobius amjiensis)

Hynobius amjiensis, also known as the Anji salamander, is an amphibian species currently categorized as endangered due to its limited geographical distribution, primarily in China. To address the critical conservation status of this species, artificial breeding is essential for population expansion. However, progress in artificial breeding efforts has been hindered by the scarcity of research on the reproductive biology of the Anji salamander. In this study, we identified 25 distinct early stages of embryo development. Additionally, we observed that Anji salamander embryos contain a lesser amount of yolk compared to other salamanders or frogs. We further discovered that the Anji salamander employs a highly competitive reproductive strategy, producing a smaller number of high-quality offspring. This strategy aims to generate adaptive individuals through intense intraspecific competition, driven by three factors: extremely confined breeding habitats, a substantial number of eggs, and a reduced yolk content. We introduce the term "mass escape" to describe this strategy, which provides a novel perspective on cannibalization, focusing on the consumption of specific body parts rather than a single-individual activity. This study offers valuable insights into artificial breeding techniques designed to mitigate inherent intraspecific competitive pressure, thereby improving metamorphosis and survival rates. Additionally, it provides a comprehensive table detailing the pre-embryonic developmental stages of the Anji salamander.

RNA localization during early development of the axolotl

The asymmetric localization of biomolecules is critical for body plan development. One of the most popular model organisms for early embryogenesis studies is Xenopus laevis but there is a lack of information in other animal species. Here, we compared the early development of two amphibian species-the frog X. laevis and the axolotl Ambystoma mexicanum. This study aimed to identify asymmetrically localized RNAs along the animal-vegetal axis during the early development of A. mexicanum. For that purpose, we performed spatial transcriptome-wide analysis at low resolution, which revealed dynamic changes along the animal-vegetal axis classified into the following categories: profile alteration, de novo synthesis and degradation. Surprisingly, our results showed that many of the vegetally localized genes, which are important for germ cell development, are degraded during early development. Furthermore, we assessed the motif presence in UTRs of degraded mRNAs and revealed the enrichment of several motifs in RNAs of germ cell markers. Our results suggest novel reorganization of the transcriptome during embryogenesis of A. mexicanum to converge to the similar developmental pattern as the X. laevis.

NANOG is required to establish the competence for germ-layer differentiation in the basal tetrapod axolotl

Pluripotency defines the unlimited potential of individual cells of vertebrate embryos, from which all adult somatic cells and germ cells are derived. Understanding how the programming of pluripotency evolved has been obscured in part by a lack of data from lower vertebrates; in model systems such as frogs and zebrafish, the function of the pluripotency genes NANOG and POU5F1 have diverged. Here, we investigated how the axolotl ortholog of NANOG programs pluripotency during development. Axolotl NANOG is absolutely required for gastrulation and germ-layer commitment. We show that in axolotl primitive ectoderm (animal caps; ACs) NANOG and NODAL activity, as well as the epigenetic modifying enzyme DPY30, are required for the mass deposition of H3K4me3 in pluripotent chromatin. We also demonstrate that all 3 protein activities are required for ACs to establish the competency to differentiate toward mesoderm. Our results suggest the ancient function of NANOG may be establishing the competence for lineage differentiation in early cells. These observations provide insights into embryonic development in the tetrapod ancestor from which terrestrial vertebrates evolved.

A Nile rat transcriptomic landscape across 22 organs by ultra-deep sequencing and comparative RNA-seq pipeline (CRSP)

RNA sequencing (RNA-seq) has been a widely used high-throughput method to characterize transcriptomic dynamics spatiotemporally. However, RNA-seq data analysis pipelines typically depend on either a sequenced genome and/or corresponding reference transcripts. This limitation is a challenge for species lacking sequenced genomes and corresponding reference transcripts. The Nile rat (Arvicanthis niloticus) has two key features - it is daytime active, and it is prone to diet-induced diabetes, which makes it more similar to humans than regular laboratory rodents. However, at the time of this study, neither a Nile rat genome nor a reference transcript set were available, making it technically challenging to perform large-scale RNA-seq based transcriptomic studies. This genome-independent work progressed concurrently with our generation of a Nile rat genome. A well-annotated genome requires several iterations of manually reviewing curated transcripts and takes years to achieve. Here, we developed a Comparative RNA-Seq Pipeline (CRSP), integrating a comparative species strategy independent of a specific sequenced genome or species-matched reference transcripts. We performed benchmarking to validate that our CRSP tool can accurately quantify gene expression levels. In this study, we generated the first ultra-deep (2.3 billion × 2 paired-end) Nile rat RNA-seq data from 59 biopsy samples representing 22 major organs, providing a unique resource and spatial gene expression reference for Nile rat researchers. Importantly, CRSP is not limited to the Nile rat species and can be applied to any species without prior genomic knowledge. To facilitate a general use of CRSP, we also characterized the number of RNA-seq reads required for accurate estimation via simulation studies. CRSP and documents are available at: https://github.com/pjiang1105/CRSP.

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.

The amazing and anomalous axolotls as scientific models

Ambystoma mexicanum (axolotl) embryos and juveniles have been used as model organisms for developmental and regenerative research for many years. This neotenic aquatic species maintains the unique capability to regenerate most, if not all, of its tissues well into adulthood. With large externally developing embryos, axolotls were one of the original model species for developmental biology. However, increased access to, and use of, organisms with sequenced and annotated genomes, such as Xenopus laevis and tropicalis and Danio rerio, reduced the prevalence of axolotls as models in embryogenesis studies. Recent sequencing of the large axolotl genome opens up new possibilities for defining the recipes that drive the formation and regeneration of tissues like the limbs and spinal cord. However, to decode the large Ambystoma mexicanum genome will take a herculean effort, community resources, and the development of novel techniques. Here, we provide an updated axolotl-staging chart ranging from 1-cell stage to immature adult paired with a perspective on both historical and current axolotl research that spans from their use in early studies of development to the recent cutting-edge research, employment of transgenesis, high resolution imaging, and study of mechanisms deployed in regeneration. This article is protected by copyright. All rights reserved.

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.

TimeMeter assesses temporal gene expression similarity and identifies differentially progressing genes

Comparative time series transcriptome analysis is a powerful tool to study development, evolution, aging, disease progression and cancer prognosis. We develop TimeMeter, a statistical method and tool to assess temporal gene expression similarity, and identify differentially progressing genes where one pattern is more temporally advanced than the other. We apply TimeMeter to several datasets, and show that TimeMeter is capable of characterizing complicated temporal gene expression associations. Interestingly, we find: (i) the measurement of differential progression provides a novel feature in addition to pattern similarity that can characterize early developmental divergence between two species; (ii) genes exhibiting similar temporal patterns between human and mouse during neural differentiation are under strong negative (purifying) selection during evolution; (iii) analysis of genes with similar temporal patterns in mouse digit regeneration and axolotl blastema differentiation reveals common gene groups for appendage regeneration with potential implications in regenerative medicine.

Endodermal Maternal Transcription Factors Establish Super-Enhancers during Zygotic Genome Activation

Elucidation of the sequence of events underlying the dynamic interaction between transcription factors and chromatin states is essential. Maternal transcription factors function at the top of the regulatory hierarchy to specify the primary germ layers at the onset of zygotic genome activation (ZGA). We focus on the formation of endoderm progenitor cells and examine the interactions between maternal transcription factors and chromatin state changes underlying the cell specification process. Endoderm-specific factors Otx1 and Vegt together with Foxh1 orchestrate endoderm formation by coordinated binding to select regulatory regions. These interactions occur before the deposition of enhancer histone marks around the regulatory regions, and these TFs recruit RNA polymerase II, regulate enhancer activity, and establish super-enhancers associated with important endodermal genes. Therefore, maternal transcription factors Otx1, Vegt, and Foxh1 combinatorially regulate the activity of super-enhancers, which in turn activate key lineage-specifying genes during ZGA.

How do maternal transcription factors interact with chromatin regions to coordinate the endodermal gene regulatory program? Paraiso et al. demonstrate that combinatorial binding of maternal Otx1, Vegt, and Foxh1 to select enhancers and super-enhancers in the genome controls endodermal cell fate specification during zygotic gene activation.

Transcriptome analysis of axolotl oropharyngeal explants during taste bud differentiation stages

The Mexican salamander, Ambystoma mexicanum (Axolotl), is an excellent vertebrate model system to understand development and regeneration. Studies in axolotl embryos have provided important insights into taste bud development. Taste bud specification and determination occur in the oropharyngeal endoderm of axolotl embryos during gastrulation and neurulation, respectively, whereas taste bud innervation and taste cell differentiation occur later in development. Axolotl embryos are amenable to microsurgery, and tissue explants develop readily in vitro. We performed RNA-seq analysis to investigate the differential expression of genes in oropharyngeal explants at several stages of taste cell differentiation. Since the axolotl genome has only recently been sequenced, we used a Trinity pipeline to perform de novo assembly of sequencing reads. Linear models for RNA-seq data were used to identify differentially expressed genes. We found 1234 unique genes differentially expressed during taste cell differentiation stages. We validated four of these genes using RTqPCR and performed GO functional analysis. The differential expression of these genes suggests that they may play a role in taste cell differentiation in axolotls.

Integrative Analysis of Axolotl Gene Expression Data from Regenerative and Wound Healing Limb Tissues

Axolotl (Ambystoma mexicanum) is a urodele amphibian endowed with remarkable regenerative capacities manifested in scarless wound healing and restoration of amputated limbs, which makes it a powerful experimental model for regenerative biology and medicine. Previous studies have utilized microarrays and RNA-Seq technologies for detecting differentially expressed (DE) genes in different phases of the axolotl limb regeneration. However, sufficient consistency may be lacking due to statistical limitations arising from intra-laboratory analyses. This study aims to bridge such gaps by performing an integrative analysis of publicly available microarray and RNA-Seq data from axolotl limb samples having comparable study designs using the "merging" method. A total of 351 genes were found DE in regenerative samples compared to the control in data of both technologies, showing an adjusted p-value < 0.01 and log fold change magnitudes >1. Downstream analyses illustrated consistent correlations of the directionality of DE genes within and between data of both technologies, as well as concordance with the literature on regeneration related biological processes. qRT-PCR analysis validated the observed expression level differences of five of the top DE genes. Future studies may benefit from the utilized concept and approach for enhanced statistical power and robust discovery of biomarkers of regeneration.

A comprehensive reference transcriptome resource for the Iberian ribbed newt Pleurodeles waltl, an emerging model for developmental and regeneration biology

Urodele newts have unique biological properties, notably including prominent regeneration ability. The Iberian ribbed newt, Pleurodeles waltl, is a promising model amphibian distinguished by ease of breeding and efficient transgenic and genome editing methods. However, limited genetic information is available for P. waltl. We conducted an intensive transcriptome analysis of P. waltl using RNA-sequencing to build and annotate gene models. We generated 1.2 billion Illumina reads from a wide variety of samples across 12 different tissues/organs, unfertilized egg, and embryos at eight different developmental stages. These reads were assembled into 1,395,387 contigs, from which 202,788 non-redundant ORF models were constructed. The set is expected to cover a large fraction of P. waltl protein-coding genes, as confirmed by BUSCO analysis, where 98% of universal single-copy orthologs were identified. Ortholog analyses revealed the gene repertoire evolution of urodele amphibians. Using the gene set as a reference, gene network analysis identified regeneration-, developmental-stage-, and tissue-specific co-expressed gene modules. Our transcriptome resource is expected to enhance future research employing this emerging model animal for regeneration research as well as for investigations in other areas including developmental biology, stem cell biology, and cancer research. These data are available via our portal website, iNewt (http://www.nibb.ac.jp/imori/main/).

Longitudinal 16S rRNA data derived from limb regenerative tissue samples of axolotl Ambystoma mexicanum

The Mexican axolotl (Ambystoma mexicanum) is a critically endangered species and a fruitful amphibian model for regenerative biology. Despite growing body of research on the cellular and molecular biology of axolotl limb regeneration, microbiological aspects of this process remain poorly understood. Here, we describe bacterial 16S rRNA amplicon dataset derived from axolotl limb tissue samples in the course of limb regeneration. The raw data was obtained by sequencing V3–V4 region of 16S rRNA gene and comprised 14,569,756 paired-end raw reads generated from 21 samples. Initial data analysis using DADA2 pipeline resulted in amplicon sequence variant (ASV) table containing a total of ca. 5.9 million chimera-removed, high-quality reads and a median of 296,971 reads per sample. The data constitute a useful resource for the research on the microbiological aspects of axolotl limb regeneration and will also broadly facilitate comparative studies in the developmental and conservation biology of this critically endangered species.

Design Type(s)	time series design • organism development design
Measurement Type(s)	rRNA 16S
Technology Type(s)	DNA sequencing
Factor Type(s)	biological replicate • temporal_interval
Sample Characteristic(s)	Ambystoma mexicanum • limb • laboratory environment

Machine-accessible metadata file describing the reported data (ISA-Tab format)

Effect of carotenoid class and dose on the larval growth and development of the critically endangered southern corroboree frog

Dietary carotenoids are expected to improve vertebrate growth and development, though evidence for beneficial effects remains limited. One reason for this might be that few studies have directly compared the effects of carotenoids from different classes (carotenes versus xanthophylls) at more than one dose. Here, we tested the effect of two doses of dietary β-carotene and lutein (representing two different carotenoid classes) on the growth and development of larval southern corroboree frogs (Pseudophryne corroboree). Individuals were supplemented with either β-carotene or lutein at one of two doses (0.1 mg g^-1, 1 mg g^-1), or given a diet without carotenoids (control). Each dietary treatment included 36 replicate individuals, and individuals remained on the same diet until metamorphosis (25-39 weeks). We measured larval survival, larval growth (body length), time to metamorphosis, metamorphic body size (mass and SVL), and body condition. Lutein had no detectable effect on larval growth and development. However, larvae receiving a high dose (1 mg g^-1) of β-carotene metamorphosed significantly faster than all other dietary treatments, despite no significant differences in growth rate. This result indicates that β-carotene supplementation in P. corroboree has positive effects on development independent of growth effects. Our study provides new evidence for differential effects of carotenoid class and dose on vertebrate development. From a conservation perspective, our findings are expected to assist with the recovery of P. corroboree by expediting the generation of frogs required for the maintenance of captive insurance colonies, or the provision of frogs for release. More broadly, our study highlights the potential for dietary manipulation to assist with the ex situ management of threatened amphibian species worldwide.

AP-1cFos/JunB/miR-200a regulate the pro-regenerative glial cell response during axolotl spinal cord regeneration

Salamanders have the remarkable ability to functionally regenerate after spinal cord transection. In response to injury, GFAP⁺ glial cells in the axolotl spinal cord proliferate and migrate to replace the missing neural tube and create a permissive environment for axon regeneration. Molecular pathways that regulate the pro-regenerative axolotl glial cell response are poorly understood. Here we show axolotl glial cells up-regulate AP-1^cFos/JunB after injury, which promotes a pro-regenerative glial cell response. Injury induced upregulation of miR-200a in glial cells supresses c-Jun expression in these cells. Inhibition of miR-200a during regeneration causes defects in axonal regrowth and transcriptomic analysis revealed that miR-200a inhibition leads to differential regulation of genes involved with reactive gliosis, the glial scar, extracellular matrix remodeling and axon guidance. This work identifies a unique role for miR-200a in inhibiting reactive gliosis in axolotl glial cells during spinal cord regeneration.

Keith Sabin et al. showed that upregulation of the AP-1 complex, composed of c-Fos and JunB, in the axolotl spinal cord promotes a pro-regenerative glial cell response. This response is impaired by inhibition of miR-200a; suggesting an important role for this microRNA in axolotl spinal cord regeneration.

Trendy: segmented regression analysis of expression dynamics in high-throughput ordered profiling experiments

BACKGROUND

High-throughput expression profiling experiments with ordered conditions (e.g. time-course or spatial-course) are becoming more common for studying detailed differentiation processes or spatial patterns. Identifying dynamic changes at both the individual gene and whole transcriptome level can provide important insights about genes, pathways, and critical time points.

RESULTS

We present an R package, Trendy, which utilizes segmented regression models to simultaneously characterize each gene's expression pattern and summarize overall dynamic activity in ordered condition experiments. For each gene, Trendy finds the optimal segmented regression model and provides the location and direction of dynamic changes in expression. We demonstrate the utility of Trendy to provide biologically relevant results on both microarray and RNA-sequencing (RNA-seq) datasets.

CONCLUSIONS

Trendy is a flexible R package which characterizes gene-specific expression patterns and summarizes changes of global dynamics over ordered conditions. Trendy is freely available on Bioconductor with a full vignette at https://bioconductor.org/packages/release/bioc/html/Trendy.html .

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.

A Novel Approach to Single Cell RNA-Sequence Analysis Facilitates In Silico Gene Reporting of Human Pluripotent Stem Cell-Derived Retinal Cell Types

Cell type-specific investigations commonly use gene reporters or single-cell analytical techniques. However, reporter line development is arduous and generally limited to a single gene of interest, while single-cell RNA (scRNA)-sequencing (seq) frequently yields equivocal results that preclude definitive cell identification. To examine gene expression profiles of multiple retinal cell types derived from human pluripotent stem cells (hPSCs), we performed scRNA-seq on optic vesicle (OV)-like structures cultured under cGMP-compatible conditions. However, efforts to apply traditional scRNA-seq analytical methods based on unbiased algorithms were unrevealing. Therefore, we developed a simple, versatile, and universally applicable approach that generates gene expression data akin to those obtained from reporter lines. This method ranks single cells by expression level of a bait gene and searches the transcriptome for genes whose cell-to-cell rank order expression most closely matches that of the bait. Moreover, multiple bait genes can be combined to refine datasets. Using this approach, we provide further evidence for the authenticity of hPSC-derived retinal cell types. Stem Cells 2018;36:313-324.

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 .

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.

Advances in Decoding Axolotl Limb Regeneration

Humans and other mammals are limited in their natural abilities to regenerate lost body parts. By contrast, many salamanders are highly regenerative and can spontaneously replace lost limbs even as adults. Because salamander limbs are anatomically similar to human limbs, knowing how they regenerate should provide important clues for regenerative medicine. Although interest in understanding the mechanics of this process has never wavered, until recently researchers have been vexed by seemingly impenetrable logistics of working with these creatures at a molecular level. Chief among the problems has been the very large size of salamander genomes, and not a single salamander genome has been fully sequenced to date. Recently the enormous gap in sequence information has been bridged by approaches that leverage mRNA as the starting point. Together with functional experimentation, these data are rapidly enabling researchers to finally uncover the molecular mechanisms underpinning the astonishing biological process of limb regeneration.