Transcriptional landscapes of Axolotl (Ambystoma mexicanum)

Gene functions of the Ambystoma altamirani skin microbiome vary across space and time but potential antifungal genes are widespread and prevalent

Amphibian skin microbiomes can play a critical role in host survival against emerging diseases by protecting their host against pathogens. While a plethora of biotic and abiotic factors have been shown to influence the taxonomic diversity of amphibian skin microbiomes it remains unclear whether functional genomic diversity varies in response to temporal and environmental factors. Here we applied a metagenomic approach to evaluate whether seasonality, distinct elevations/sites, and pathogen presence influenced the functional genomic diversity of the A. altamirani skin microbiome. We obtained a gene catalogue of 92 107 nonredundant annotated genes and a set of 50 unique metagenome assembled genomes (MAGs). Our analysis showed that genes linked to general and potential antifungal traits significantly differed across seasons and sampling locations at different elevations. Moreover, we found that the functional genomic diversity of A. altamirani skin microbiome differed between B. dendrobatidis infected and not infected axolotls only during winter, suggesting an interaction between seasonality and pathogen infection. In addition, we identified the presence of genes and biosynthetic gene clusters (BGCs) linked to potential antifungal functions such as biofilm formation, quorum sensing, secretion systems, secondary metabolite biosynthesis, and chitin degradation. Interestingly genes linked to these potential antifungal traits were mainly identified in Burkholderiales and Chitinophagales MAGs. Overall, our results identified functional traits linked to potential antifungal functions in the A. altamirani skin microbiome regardless of variation in the functional diversity across seasons, elevations/sites, and pathogen presence. Our findings suggest that potential antifungal traits found in Burkholderiales and Chitinophagales taxa could be related to the capacity of A. altamirani to survive in the presence of Bd, although further experimental analyses are required to test this hypothesis.

Applying a Knock-In Strategy to Create Reporter-Tagged Knockout Alleles in Axolotls

Tetrapod species axolotls exhibit the powerful capacity to fully regenerate their tail and limbs upon injury, hence serving as an excellent model organism in regenerative biology research. Developing proper molecular and genetic tools in axolotls is an absolute necessity for deep dissection of tissue regeneration mechanisms. Previously, CRISPR-/Cas9-based knockout and targeted gene knock-in approaches have been established in axolotls, allowing genetically deciphering gene function, labeling, and tracing particular types of cells. Here, we further extend the CRISPR/Cas9 technology application and describe a method to create reporter-tagged knockout allele in axolotls. This method combines gene knockout and knock-in and achieves loss of function of a given gene and simultaneous labeling of cells expressing this particular gene, that allows identification, tracking of the "knocking out" cells. Our method offers a useful gene function analysis tool to the field of axolotl developmental and regenerative research.

W, Chen H, Yu C, Yang L, Fu Y, Li J, Fu S, Sun Z, Fei L, Guo Q, Wang J, Xiao Y, Wang X, Zhang P, Ma L, Ge D, Xu S, Caballero-Pérez J, Cruz-Ramírez A, Zhou Y, Chen M, Fei JF, Han X, Guo G. Construction of the axolotl cell landscape using combinatorial hybridization sequencing at single-cell resolution

The Mexican axolotl (Ambystoma mexicanum) is a well-established tetrapod model for regeneration and developmental studies. Remarkably, neotenic axolotls may undergo metamorphosis, a process that triggers many dramatic changes in diverse organs, accompanied by gradually decline of their regeneration capacity and lifespan. However, the molecular regulation and cellular changes in neotenic and metamorphosed axolotls are still poorly investigated. Here, we develop a single-cell sequencing method based on combinatorial hybridization to generate a tissue-based transcriptomic landscape of the neotenic and metamorphosed axolotls. We perform gene expression profiling of over 1 million single cells across 19 tissues to construct the first adult axolotl cell landscape. Comparison of single-cell transcriptomes between the tissues of neotenic and metamorphosed axolotls reveal the heterogeneity of non-immune parenchymal cells in different tissues and established their regulatory network. Furthermore, we describe dynamic gene expression patterns during limb development in neotenic axolotls. This system-level single-cell analysis of molecular characteristics in neotenic and metamorphosed axolotls, serves as a resource to explore the molecular identity of the axolotl and facilitates better understanding of metamorphosis.

The Mexican axolotl is a well-established tetrapod model for regeneration and development. Here the authors report a scRNA-seq method to profile neotenic, metamorphic and limb development stages, highlighting unique perturbation patterns of cell type-related gene expression throughout metamorphosis.

Microtubule Dynamics Following Central and Peripheral Nervous System Axotomy

Disturbance in the neuronal network leads to instability in the microtubule (MT) railroad of axons, causing hindrance in the intra-axonal transport and making it difficult to re-establish the broken network. Peripheral nervous system (PNS) neurons can stabilize their MTs, leading to the formation of regeneration-promoting structures called "growth cones". However, central nervous system (CNS) neurons lack this intrinsic reparative capability and, instead, form growth-incompetent structures called "retraction bulbs", which have a disarrayed MT network. It is evident from various studies that although axonal regeneration depends on both cell-extrinsic and cell-intrinsic factors, any therapy that aims at axonal regeneration ultimately converges onto MTs. Understanding the neuronal MT dynamics will help develop effective therapeutic strategies in diseases where the MT network gets disrupted, such as spinal cord injury, traumatic brain injury, multiple sclerosis, and amyotrophic lateral sclerosis. It is also essential to know the factors that aid or inhibit MT stabilization. In this review, we have discussed the MT dynamics postaxotomy in the CNS and PNS, and factors that can directly influence MT stability in various diseases.

Evidence of requirement for homologous-mediated DNA repair during Ambystoma mexicanum limb regeneration

BACKGROUND

Limb regeneration in the axolotl is achieved by epimorphosis, thus depending on the blastema formation, a mass of progenitor cells capable of proliferating and differentiating to recover all lost structures functionally. During regeneration, the blastema cells accelerate the cell cycle and duplicate its genome, which is inherently difficult to replicate because of its length and composition, thus being prone to suffer double-strand breaks.

RESULTS

We identified and characterized two remarkable components of the homologous recombination repair pathway (Amex.RAD51 and Amex.MRE11), which were heterologously expressed, biochemically characterized, and inhibited by specific chemicals. These same inhibitors were applied at different time points after amputation to study their effects during limb regeneration. We observed an increase in cellular senescent accompanied by a slight delay in regeneration at 28 days post-amputation regenerated tissues; moreover, inhibitors caused a rise in the double-strand break signaling as a response to the inhibition of the repair mechanisms.

CONCLUSIONS

We confirmed the participation and importance of homologous recombination during limb regeneration. Where the chemical inhibition induces double-strand breaks that lead to DNA damage associated senescence, or in an alternatively way, this damage could be possibly repaired by a different DNA repair pathway, permitting proper regeneration and avoiding senescence. This article is protected by copyright. All rights reserved.

Gene Conversion Facilitates the Adaptive Evolution of Self-Resistance in Highly Toxic Newts

Reconstructing the histories of complex adaptations and identifying the evolutionary mechanisms underlying their origins are two of the primary goals of evolutionary biology. Taricha newts, which contain high concentrations of the deadly toxin tetrodotoxin (TTX) as an antipredator defense, have evolved resistance to self-intoxication, which is a complex adaptation requiring changes in six paralogs of the voltage-gated sodium channel (Na_v) gene family, the physiological target of TTX. Here, we reconstruct the origins of TTX self-resistance by sequencing the entire Na_v gene family in newts and related salamanders. We show that moderate TTX resistance evolved early in the salamander lineage in three of the six Na_v paralogs, preceding the proposed appearance of tetrodotoxic newts by ∼100 My. TTX-bearing newts possess additional unique substitutions across the entire Na_v gene family that provide physiological TTX resistance. These substitutions coincide with signatures of positive selection and relaxed purifying selection, as well as gene conversion events, that together likely facilitated their evolution. We also identify a novel exon duplication within Na_v1.4 encoding an expressed TTX-binding site. Two resistance-conferring changes within newts appear to have spread via nonallelic gene conversion: in one case, one codon was copied between paralogs, and in the second, multiple substitutions were homogenized between the duplicate exons of Na_v1.4. Our results demonstrate that gene conversion can accelerate the coordinated evolution of gene families in response to a common selection pressure.

Towards comparative analyses of salamander limb regeneration

Among tetrapods, only salamanders can regenerate their limbs and tails throughout life. This amazing regenerative ability has attracted the attention of scientists for hundreds of years. Now that large, salamander genomes are beginning to be sequenced for the first time, omics tools and approaches can be used to integrate new perspectives into the study of tissue regeneration. Here we argue the need to move beyond the primary salamander models to investigate regeneration in other species. Salamanders at first glance come across as a phylogenetically conservative group that has not diverged greatly from their ancestors. While salamanders do present ancestral characteristics of basal tetrapods, including the ability to regenerate limbs, data from fossils and data from studies that have tested for species differences suggest there may be considerable variation in how salamanders develop and regenerate their limbs. We review the case for expanded studies of salamander tissue regeneration and identify questions and approaches that are most likely to reveal commonalities and differences in regeneration among species. We also address challenges that confront such an initiative, some of which are regulatory and not scientific. The time is right to gain evolutionary perspective about mechanisms of tissue regeneration from comparative studies of salamander species.

Functional Characterization of the Lin28/let-7 Circuit During Forelimb Regeneration in Ambystoma mexicanum and Its Influence on Metabolic Reprogramming

The axolotl (Ambystoma mexicanum) is a caudate amphibian, which has an extraordinary ability to restore a wide variety of damaged structures by a process denominated epimorphosis. While the origin and potentiality of progenitor cells that take part during epimorphic regeneration are known to some extent, the metabolic changes experienced and their associated implications, remain unexplored. However, a circuit with a potential role as a modulator of cellular metabolism along regeneration is that formed by Lin28/let-7. In this study, we report two Lin28 paralogs and eight mature let-7 microRNAs encoded in the axolotl genome. Particularly, in the proliferative blastema stage amxLin28B is more abundant in the nuclei of blastemal cells, while the microRNAs amx-let-7c and amx-let-7a are most downregulated. Functional inhibition of Lin28 factors increase the levels of most mature let-7 microRNAs, consistent with an increment of intermediary metabolites of the Krebs cycle, and phenotypic alterations in the outgrowth of the blastema. In summary, we describe the primary components of the Lin28/let-7 circuit and their function during axolotl regeneration, acting upstream of metabolic reprogramming events.

Conservation analysis of core cell cycle regulators and their transcriptional behavior during limb regeneration in Ambystoma mexicanum

Ambystoma mexicanum (axolotl) has been one of the major experimental models for the study of regeneration during the past 100 years. Axolotl limb regeneration takes place through a multi-stage and complex developmental process called epimorphosis that involves diverse events of cell reprogramming. Such events start with dedifferentiation of somatic cells and the proliferation of quiescent stem cells to generate a population of proliferative cells called blastema. Once the blastema reaches a mature stage, cells undergo progressive differentiation into the diverse cell lineages that will form the new limb. Such pivotal cell reprogramming phenomena depend on the fine-tuned regulation of the cell cycle in each regeneration stage, where cell populations display specific proliferative capacities and differentiation status. The axolotl genome has been fully sequenced and released recently, and diverse RNA-seq approaches have also been generated, enabling the identification and conservatory analysis of core cell cycle regulators in this species. We report here our results from such analyses and present the transcriptional behavior of key regulatory factors during axolotl limb regeneration. We also found conserved protein interactions between axolotl Cyclin Dependent Kinases 2, 4 and 6 and Cyclins type D and E. Canonical CYC-CDK interactions that play major roles in modulating cell cycle progression in eukaryotes.

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.

Exploring the role of microRNAs in axolotl regeneration

The axolotl, Ambystoma mexicanum, is used extensively for research in developmental biology, particularly for its ability to regenerate and restore lost organs, including in the nervous system, to full functionality. Regeneration in mammals typically depends on the healing process and scar formation with limited replacement of lost tissue. Other organisms, such as spiny mice (Acomys cahirinus), salamanders, and zebrafish, are able to regenerate some damaged body components. Blastema is a tissue that is formed after tissue injury in such organisms and is composed of progenitor cells or dedifferentiated cells that differentiate into various cell types during regeneration. Thus, identifying the molecules responsible for initiation of blastema formation is an important aspect for understanding regeneration. Introns, a major source of noncoding RNAs (ncRNAs), have characteristic sizes in the axolotl, particularly in genes associated with development. These ncRNAs, particularly microRNAs (miRNAs), exhibit dynamic regulation during regeneration. These miRNAs play an essential role in timing and control of gene expression to order and organize processes necessary for blastema creation. Master keys or molecules that underlie the remarkable regenerative abilities of the axolotl remain to be fully explored and exploited. Further and ongoing research on regeneration promises new knowledge that may allow improved repair and renewal of human tissues.

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.

Response of transgenic tobacco overexpressing the CchGLP gene to cadmium and aluminium: phenotypic and microRNAs expression changes

Transgenic tobacco (N. tabacum cv. Xanthi nc) expressing Capsicum chinense CchGLP gene that encodes an Mn-SOD, constitutively produces hydrogen peroxide that increase endogenous ROS levels. Previous studies using these plants against geminivirus infections as well as drought stress confirmed that CchGLP expression conferred resistance against biotic and abiotic stresses. Cadmium (Cd) and Aluminium (Al) contamination in soils are a major ecological concern since they are two of the most widespread toxic elements in terrestrial environments. Trying to explore additional possible tolerance to another stresses in these plants, the aim of this work was to analyse the response to cadmium and aluminium salts during germination and early stages of plantlet development and a differential transcriptome of microRNAs (miRNAs) expression in expressing CchGLP transgenic lines and an azygote non-CchGLP expressing line. Plants were grown in vitro with addition of CdCl₂ and AlCl₃ at three different concentrations: 100, 300 and 500 μM and 50, 150 and 300 μM, respectively. The results showed higher tolerance to Cd and Al salts evaluated in two CchGLP-expressing transgenic lines L8 and L26 in comparison with the azygous non-CchGLP expressing line L1. Interestingly, L8 under Al stress presented vigorous roots and development of radicular hairs in comparison with azygous control (L1). Differentially expressed miRNAs in the comparison between L8 and L1 were associated with up and down-regulation of target genes related with structural molecule activity and ribosome constituents, as well as down-regulation in proton-transporting V-type ATPase (Vacuolar ATPase or V-ATPase). Moreover, KEGG analysis of the target genes for the differentially expressed miRNAs, led to identification of genes related with metabolic pathways and biosynthesis of secondary metabolites. One possible explanation of the tolerance to Cd and Al displayed in the transgenic tobaccos evaluated, might involve the fact that several down-regulated miRNAs, were found associated with target genes expressing V-ATPase. Specifically, miR7904-5p was down regulated and related with the up-regulation of one V-ATPase. The expression levels of these genes was confirmed by qRT-PCR assays, thus suggesting that a cation transport activity driven by the V-ATPases-dependent proton motive force, might significantly contribute as one mechanism for Cd and Al detoxification by vacuolar compartmentation in these transgenic tobacco plants.

Too Many False Targets for MicroRNAs: Challenges and Pitfalls in Prediction of miRNA Targets and Their Gene Ontology in Model and Non-model Organisms

Short ("seed") or extended base pairing between microRNAs (miRNAs) and their target RNAs enables post-transcriptional silencing in many organisms. These interactions allow the computational prediction of potential targets. In model organisms, predicted targets are frequently validated experimentally; hence meaningful miRNA-regulated processes are reported. However, in non-models, these reports mostly rely on computational prediction alone. Many times, further bioinformatic analyses such as Gene Ontology (GO) enrichment are based on these in silico projections. Here such approaches are reviewed, their caveats are highlighted and the ease of picking false targets from predicted lists is demonstrated. Discoveries that shed new light on how miRNAs evolved to regulate targets in various phyletic groups are discussed, in addition to the pitfalls of target identification in non-model organisms. The goal is to prevent the misuse of bioinformatic tools, as they cannot bypass the biological understanding of miRNA-target regulation.

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.

Using transcriptomics to enable a plethodontid salamander (Bolitoglossa ramosi) for limb regeneration research

Background

Tissue regeneration is widely distributed across the tree of life. Among vertebrates, salamanders possess an exceptional ability to regenerate amputated limbs and other complex structures. Thus far, molecular insights about limb regeneration have come from a relatively limited number of species from two closely related salamander families. To gain a broader perspective on the molecular basis of limb regeneration and enhance the molecular toolkit of an emerging plethodontid salamander (Bolitoglossa ramosi), we used RNA-Seq to generate a de novo reference transcriptome and identify differentially expressed genes during limb regeneration.

Results

Using paired-end Illumina sequencing technology and Trinity assembly, a total of 433,809 transcripts were recovered and we obtained functional annotation for 142,926 non-redundant transcripts of the B. ramosi de novo reference transcriptome. Among the annotated transcripts, 602 genes were identified as differentially expressed during limb regeneration. This list was further processed to identify a core set of genes that exhibit conserved expression changes between B. ramosi and the Mexican axolotl (Ambystoma mexicanum), and presumably their common ancestor from approximately 180 million years ago.

Conclusions

We identified genes from B. ramosi that are differentially expressed during limb regeneration, including multiple conserved protein-coding genes and possible putative species-specific genes. Comparative analyses reveal a subset of genes that show similar patterns of expression with ambystomatid species, which highlights the importance of developing comparative gene expression data for studies of limb regeneration among salamanders.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5076-0) contains supplementary material, which is available to authorized users.

Forever young: Endocrinology of paedomorphosis in the Mexican axolotl (Ambystoma mexicanum)

The Mexican axolotl (Ambystoma mexicanum) is a salamander species that does not undergo metamorphosis, resulting in the retention of juvenile characteristics in the mature breeding stage (paedomorphosis). Here we review the endocrinological studies investigating the proximate cause of axolotl paedomorphosis with a focus on the hypothalamo-pituitary-thyroid (HPT) axis. It is well established that axolotl paedomorphosis is a consequence of low activity of the HPT axis. The pituitary hormone thyrotropin (TSH) is capable of inducing metamorphosis in the axolotl, which indicates that all processes and interactions in the HPT axis below the pituitary level are functional, but that TSH release is impaired. In metamorphosing species, TSH secretion is largely controlled by the hypothalamic neuropeptide corticotropin-releasing hormone (CRH), which seems to have lost its thyrotropic activity in the axolotl. However, preliminary experiments have not yet confirmed a role for faulty CRH signalling in axolotl paedomorphosis. Other hypothalamic factors and potential pituitary inhibitors need to be investigated to identify their roles in amphibian metamorphosis and axolotl paedomorphosis.

Comparative transcriptomics of limb regeneration: Identification of conserved expression changes among three species of Ambystoma

Transcriptome studies are revealing the complex gene expression basis of limb regeneration in the primary salamander model - Ambystoma mexicanum (axolotl). To better understand this complexity, there is need to extend analyses to additional salamander species. Using microarray and RNA-Seq, we performed a comparative transcriptomic study using A. mexicanum and two other ambystomatid salamanders: A. andersoni, and A. maculatum. Salamanders were administered forelimb amputations and RNA was isolated and analyzed to identify 405 non-redundant genes that were commonly, differentially expressed 24 h post amputation. Many of the upregulated genes are predicted to function in wound healing and developmental processes, while many of the downregulated genes are typically expressed in muscle. The conserved transcriptional changes identified in this study provide a high-confidence dataset for identifying factors that simultaneous orchestrate wound healing and regeneration processes in response to injury, and more generally for identifying genes that are essential for salamander limb regeneration.