hoxc12/c13 as key regulators for rebooting the developmental program in Xenopus limb regeneration

During organ regeneration, after the initial responses to injury, gene expression patterns similar to those in normal development are reestablished during subsequent morphogenesis phases. This supports the idea that regeneration recapitulates development and predicts the existence of genes that reboot the developmental program after the initial responses. However, such rebooting mechanisms are largely unknown. Here, we explore core rebooting factors that operate during Xenopus limb regeneration. Transcriptomic analysis of larval limb blastema reveals that hoxc12/c13 show the highest regeneration specificity in expression. Knocking out each of them through genome editing inhibits cell proliferation and expression of a group of genes that are essential for development, resulting in autopod regeneration failure, while limb development and initial blastema formation are not affected. Furthermore, the induction of hoxc12/c13 expression partially restores froglet regenerative capacity which is normally very limited compared to larval regeneration. Thus, we demonstrate the existence of genes that have a profound impact alone on rebooting of the developmental program in a regeneration-specific manner.

A CRISPR-Cas9-mediated versatile method for targeted integration of a fluorescent protein gene to visualize endogenous gene expression in Xenopus laevis

Xenopus laevis is a widely used model organism in developmental and regeneration studies. Despite several reports regarding targeted integration techniques in Xenopus, there is still room for improvement of them, especially in creating reporter lines that rely on endogenous regulatory enhancers/promoters. We developed a CRISPR-Cas9-based simple method to efficiently introduce a fluorescent protein gene into 5' untranslated regions (5'UTRs) of target genes in Xenopus laevis. A donor plasmid DNA encoding an enhanced green fluorescent protein (eGFP) flanked by a genomic fragment ranging from 66 bp to 878 bp including target 5'UTR was co-injected into fertilized eggs with a single guide RNA and Cas9 protein. Injections for krt12.2.L, myod1.S, sox2.L or brevican.S resulted in embryos expressing eGFP fluorescence in a tissue-specific manner, recapitulating endogenous expression of target genes. Integrations of the donor DNA into the target regions were examined by genotyping PCR for the eGFP-expressing embryos. The rate of embryos expressing the specific eGFP varied from 2.1% to 13.2% depending on the target locus and length of the genomic fragment in the donor plasmids. Germline transmission of an integrated DNA was observed. This simple method provides a powerful tool for exploring gene expression and function in developmental and regeneration research in X. laevis.

Protocols for transgenesis at a safe harbor site in the Xenopus laevis genome using CRISPR-Cas9

We have established a new transgenesis protocol based on CRISPR-Cas9, "New and Easy XenopusTransgenesis (NEXTrans)," and identified a novel safe harbor site in African clawed frogs, Xenopus laevis. We describe steps in detail for the construction of NEXTrans plasmid and guide RNA, CRISPR-Cas9-mediated NEXTrans plasmid integration into the locus, and its validation by genomic PCR. This improved strategy allows us to simply generate transgenic animals that stably express the transgene. For complete details on the use and execution of this protocol, please refer to Shibata et al. (2022).1.

CRISPR-Cas9-Based Functional Analysis in Amphibians: Xenopus laevis, Xenopus tropicalis, and Pleurodeles waltl

Amphibians have made many fundamental contributions to our knowledge, from basic biology to biomedical research on human diseases. Current genome editing tools based on the CRISPR-Cas system enable us to perform gene functional analysis in vivo, even in non-model organisms. We introduce here a highly efficient and easy protocol for gene knockout, which can be used in three different amphibians seamlessly: Xenopus laevis, Xenopus tropicalis, and Pleurodeles waltl. As it utilizes Cas9 ribonucleoprotein complex (RNP) for injection, this cloning-free method enables researchers to obtain founder embryos with a nearly complete knockout phenotype within a week. To evaluate somatic mutation rate and its correlation to the phenotype of a Cas9 RNP-injected embryo (crispant), we also present accurate and cost-effective genotyping methods using pooled amplicon-sequencing and a user-friendly web-based tool.

CRISPR/Cas9-based simple transgenesis in Xenopus laevis

Transgenic techniques have greatly increased our understanding of the transcriptional regulation of target genes through live reporter imaging, as well as the spatiotemporal function of a gene using loss- and gain-of-function constructs. In Xenopus species, two well-established transgenic methods, restriction enzyme-mediated integration and I-SceI meganuclease-mediated transgenesis, have been used to generate transgenic animals. However, donor plasmids are randomly integrated into the Xenopus genome in both methods. Here, we established a new and simple targeted transgenesis technique based on CRISPR/Cas9 in Xenopus laevis. In this method, Cas9 ribonucleoprotein (RNP) targeting a putative harbor site (the transforming growth factor beta receptor 2-like (tgfbr2l) locus) and a preset donor plasmid DNA were co-injected into the one-cell stage embryos of X. laevis. Approximately 10% of faithful reporter expression was detected in F0 crispants in a promoter/enhancer-specific manner. Importantly, efficient germline transmission and stable transgene expression were observed in the F1 offspring. The simplicity of this method only required preparation of a donor vector containing the tgfbr2l genome fragment and Cas9 RNP targeting this site, which are common experimental procedures used in Xenopus laboratories. Our improved technique allows the simple generation of transgenic X. laevis, so is expected to become a powerful tool for reporter assay and gene function analysis.

Amiodarone bioconcentration and suppression of metamorphosis in Xenopus

Trace concentrations of a number of pharmaceutically active compounds have been detected in the aquatic environment in many countries, where they are thought to have the potential to exert adverse effects on non-target organisms. Amiodarone (AMD) is one such high-risk compound commonly used in general hospitals. AMD is known to alter normal thyroid hormone (TH) function, although little information is available regarding the specific mechanism by which this disruption occurs. Anuran tadpole metamorphosis is a TH-controlled developmental process and has proven to be useful as a screening tool for environmental pollutants suspected of disrupting TH functions. In the present study, our objective was to clarify the effects of AMD on Xenopus metamorphosis as well as to assess the bioconcentration of this pharmaceutical in the liver. We found that AMD suppressed spontaneous metamorphosis, including tail regression and hindlimb elongation in pro-metamorphic stage tadpoles, which is controlled by endogenous circulating TH, indicating that AMD is a TH antagonist. In transgenic X. laevis tadpoles carrying plasmid DNA containing TH-responsive element (TRE) and a 5'-upstream promoter region of the TH receptor (TR) βA1 gene linked to a green fluorescent protein (EGFP) gene, triiodothyronine (T₃) exposure induced a strong EGFP expression in the hind limbs, whereas the addition of AMD to T₃ suppressed EGFP expression, suggesting that this drug interferes with the binding of T₃ to TR, leading to the inhibition of TR-mediated gene expression. We also found AMD to be highly bioconcentrated in the liver of pro-metamorphic X. tropicalis tadpoles, and we monitored hepatic accumulation of this drug using mass spectrometry imaging (MSI). Our findings suggest that AMD imposes potential risk to aquatic wildlife by disrupting TH homeostasis, with further possibility of accumulating in organisms higher up in the food chain.

A simple and practical workflow for genotyping of CRISPR-Cas9-based knockout phenotypes using multiplexed amplicon sequencing

Founder animals carrying high proportions of somatic mutation induced by CRISPR-Cas9 enable a rapid and scalable strategy for the functional screening of numerous target genes in vivo. In this functional screening, genotyping using pooled amplicons with next-generation sequencing is the most suitable approach for large-scale management of multiple samples and accurate evaluation of the efficiency of Cas9-induced somatic mutations at target sites. Here, we present a simple workflow for genotyping of multiple CRISPR-Cas9-based knockout founders by pooled amplicon sequencing. Using custom barcoded primers, pooled amplicons from multiple individuals can be run in a single-indexed library on the Illumina MiSeq platform. Additionally, a user-friendly web tool, CLiCKAR, is available to simultaneously perform demultiplexing of pooled sequence data and evaluation of somatic mutation in each phenotype. CLiCKAR provides users with practical reports regarding the positions of insertions/deletions, as well as the frameshift ratio and tables containing mutation sequences, and read counts of each phenotype, with just a few clicks by the implementation of demultiplexing for pooled sample data and calculation of the frameshift ratio. This genotyping workflow can be harnessed to evaluate genotype-phenotype correlations in CRISPR-Cas9-based loss-of-function screening of numerous target genes in various organisms.

In vitro assessment of effects of persistent organic pollutants on the transactivation of estrogen receptor α and β (ERα and ERβ) from the Baikal seal (Pusa sibirica)

To assess the effect of exposure to persistent organic pollutants (POPs) on the estrogen receptor (ER) signaling pathway in Baikal seals (Pusa sibirica), we investigated the molecular characterizations and functions of two Baikal seal ER (bsER) isoforms, bsERα and bsERβ. The bsERα and bsERβ cDNA clones isolated have an open reading frame of 595 and 530 amino acid residues, respectively. The tissue distribution analyses of bsER mRNAs showed that bsERα transcripts were primarily found in the ovary and uterus, and bsERβ in the muscle in wild Baikal seals. The immunofluorescence staining assay showed that 17β-estradiol (E₂) treatment promoted the nuclear translocation of in vitro-expressed bsERα. Transient transfection of bsERα in U2OS cells enhanced the transcription of pS2, an ER target gene of E₂. We then measured bsER-mediated transactivation potencies of POPs in an in vitro reporter gene assay system, in which a bsERα or bsERβ expression vector was transfected into COS-1 cells. For comparison, transactivation potencies of POPs on mouse ERs (mERα and mERβ) were also evaluated in the same manner. Results showed significant dose-dependent responses of bsERs and mERs when treated with p,p'-dichlorodiphenyltrichloroethane (p,p'-DDT), and p,p'-dichlorodiphenyldichloroethylene (p,p'-DDE). bsERs and mERs showed no response when exposed to polychlorinated biphenyls (PCBs) or 2,3,7,8-tetrachlorodibenzo-p-dioxin. Comparison of the dose-response curves of DDTs across species (bsERs vs. mERs) showed that bsERα had a response similar to mERα, but bsERβ was less sensitive than mERβ. Comparing the lowest observable effective concentrations of p,p'-DDT (2.8 μM) and p,p'-DDE (10 μM) for in vitro bsERα-mediated transactivation with their hepatic concentrations in wild Baikal seals indicated that some individuals accumulated these compounds at levels comparable to the effective concentrations, suggesting the potential disruption of the bsERα signaling pathway in the wild population by these compounds. Co-transfection experiments with bsER and the aryl hydrocarbon receptor (AHR) suggested that high accumulation of estrogenic compounds exerts a synergistic effect with dioxin-like congeners on ER signaling through AHR activation in the wild seal population.

Functional analysis of thyroid hormone receptor beta in Xenopus tropicalis founders using CRISPR-Cas

Amphibians provide an ideal model to study the actions of thyroid hormone (TH) in animal development because TH signaling via two TH receptors, TRα and TRβ, is indispensable for amphibian metamorphosis. However, specific roles for the TRβ isoform in metamorphosis are poorly understood. To address this issue, we generated trβ-disrupted Xenopus tropicalis tadpoles using the CRISPR-Cas system. We first established a highly efficient and rapid workflow for gene disruption in the founder generation (F0) by injecting sgRNA and Cas9 ribonucleoprotein. Most embryos showed severe mutant phenotypes carrying high somatic mutation rates. Utilizing this founder analysis system, we examined the role of trβ in metamorphosis. trβ-disrupted pre-metamorphic tadpoles exhibited mixed responsiveness to exogenous TH. Specifically, gill resorption and activation of several TH-response genes, including trβ itself and two protease genes, were impaired. However, hind limb outgrowth and induction of the TH-response genes, klf9 and fra-2, were not affected by loss of trβ Surprisingly, trβ-disrupted tadpoles were able to undergo spontaneous metamorphosis normally, except for a slight delay in tail resorption. These results indicate TRβ is not required but contributes to the timing of resorptive events of metamorphosis.

Reactivation of larval keratin gene (krt62.L) in blastema epithelium during Xenopus froglet limb regeneration

Limb regeneration is considered a form of limb redevelopment because of the molecular and morphological similarities. Forming a regeneration blastema is, in essence, creating a developing limb bud in an adult body. This reactivation of a developmental process in a mature body is worth studying. Xenopus laevis has a biphasic life cycle that involves distinct larval and adult stages. These distinct developmental stages are useful for investigating the reactivation of developmental processes in post-metamorphic frogs (froglets). In this study, we focused on the re-expression of a larval gene (krt62.L) during Xenopus froglet limb regeneration. Recently renamed krt62.L, this gene was known as the larval keratin (xlk) gene, which is specific to larval-tadpole stages. During limb regeneration in a froglet, krt62.L was re-expressed in a basal layer of blastema epithelium, where adult-specific keratin (Krt12.6.S) expression was also observable. Nerves produce important regulatory factors for amphibian limb regeneration, and also play a role in blastema formation and maintenance. The effect of nerve function on krt62.L expression could be seen in the maintenance of krt62.L expression, but not in its induction. When an epidermis-stripped limb bud was grafted in a froglet blastema, the grafted limb bud could reach the digit-forming stage. This suggests that krt62.L-positive froglet blastema epithelium is able to support the limb development process. These findings imply that the developmental process is locally reactivated in an postmetamorphic body during limb regeneration.

Dynamic changes in the interchromosomal interaction of early histone gene loci during development of sea urchin

The nuclear positioning and chromatin dynamics of eukaryotic genes are closely related to the regulation of gene expression, but they have not been well examined during early development, which is accompanied by rapid cell cycle progression and dynamic changes in nuclear organization, such as nuclear size and chromatin constitution. In this study, we focused on the early development of the sea urchin Hemicentrotus pulcherrimus and performed three-dimensional fluorescence in situ hybridization of gene loci encoding early histones (one of the types of histone in sea urchin). There are two non-allelic early histone gene loci per sea urchin genome. We found that during the morula stage, when the early histone gene expression levels are at their maximum, interchromosomal interactions were often formed between the early histone gene loci on separate chromosomes and that the gene loci were directed to locate to more interior positions. Furthermore, these interactions were associated with the active transcription of the early histone genes. Thus, such dynamic interchromosomal interactions may contribute to the efficient synthesis of early histone mRNA during the morula stage of sea urchin development.

Developmental changes in drug-metabolizing enzyme expression during metamorphosis of Xenopus tropicalis

A large number of chemicals are routinely detected in aquatic environments, and these chemicals may adversely affect aquatic organisms. Accurate risk assessment requires understanding drug-metabolizing systems in aquatic organisms because metabolism of these chemicals is a critical determinant of chemical bioaccumulation and related toxicity. In this study, we evaluated mRNA expression levels of nuclear receptors and drug-metabolizing enzymes as well as cytochrome P450 (CYP) activities in pro-metamorphic tadpoles, froglets, and adult frogs to determine how drug-metabolizing systems are altered at different life stages. We found that drug-metabolizing systems in tadpoles were entirely immature, and therefore, tadpoles appeared to be more susceptible to chemicals compared with metamorphosed frogs. On the other hand, cyp1a mRNA expression and CYP1A-like activity were higher in tadpoles. We found that thyroid hormone (TH), which increases during metamorphosis, induced CYP1A-like activity. Because endogenous TH concentration is significantly increased during metamorphosis, endogenous TH would induce CYP1A-like activity in tadpoles.

Clustered Xenopus keratin genes: A genomic, transcriptomic, and proteomic analysis

Keratin genes belong to the intermediate filament superfamily and their expression is altered following morphological and physiological changes in vertebrate epithelial cells. Keratin genes are divided into two groups, type I and II, and are clustered on vertebrate genomes, including those of Xenopus species. Various keratin genes have been identified and characterized by their unique expression patterns throughout ontogeny in Xenopus laevis; however, compilation of previously reported and newly identified keratin genes in two Xenopus species is required for our further understanding of keratin gene evolution, not only in amphibians but also in all terrestrial vertebrates. In this study, 120 putative type I and II keratin genes in total were identified based on the genome data from two Xenopus species. We revealed that most of these genes are highly clustered on two homeologous chromosomes, XLA9_10 and XLA2 in X. laevis, and XTR10 and XTR2 in X. tropicalis, which are orthologous to those of human, showing conserved synteny among tetrapods. RNA-Seq data from various embryonic stages and adult tissues highlighted the unique expression profiles of orthologous and homeologous keratin genes in developmental stage- and tissue-specific manners. Moreover, we identified dozens of epidermal keratin proteins from the whole embryo, larval skin, tail, and adult skin using shotgun proteomics. In light of our results, we discuss the radiation, diversification, and unique expression of the clustered keratin genes, which are closely related to epidermal development and terrestrial adaptation during amphibian evolution, including Xenopus speciation.

Unliganded thyroid hormone receptor α regulates developmental timing via gene repression in Xenopus tropicalis

Thyroid hormone (TH) receptor (TR) expression begins early in development in all vertebrates when circulating TH levels are absent or minimal, yet few developmental roles for unliganded TRs have been established. Unliganded TRs are expected to repress TH-response genes, increase tissue responsivity to TH, and regulate the timing of developmental events. Here we examined the role of unliganded TRα in gene repression and development in Xenopus tropicalis. We used transcription activator-like effector nuclease gene disruption technology to generate founder animals with mutations in the TRα gene and bred them to produce F1 offspring with a normal phenotype and a mutant phenotype, characterized by precocious hind limb development. Offspring with a normal phenotype had zero or one disrupted TRα alleles, and tadpoles with the mutant hind limb phenotype had two truncated TRα alleles with frame shift mutations between the two zinc fingers followed by 40-50 mutant amino acids and then an out-of-frame stop codon. We examined TH-response gene expression and early larval development with and without exogenous TH in F1 offspring. As hypothesized, mutant phenotype tadpoles had increased expression of TH-response genes in the absence of TH and impaired induction of these same genes after exogenous TH treatment, compared with normal phenotype animals. Also, mutant hind limb phenotype animals had reduced hind limb and gill responsivity to exogenous TH. Similar results in methimazole-treated tadpoles showed that increased TH-response gene expression and precocious development were not due to early production of TH. These results indicate that unliganded TRα delays developmental progression by repressing TH-response genes.

Targeted mutagenesis of multiple and paralogous genes in Xenopus laevis using two pairs of transcription activator-like effector nucleases

Transcription activator-like effector nucleases (TALENs) have been extensively used in genome editing in various organisms. In some cases, however, it is difficult to efficiently disrupt both paralogous genes using a single pair of TALENs in Xenopus laevis because of its polyploidy. Here, we report targeted mutagenesis of multiple and paralogous genes using two pairs of TALENs in X. laevis. First, we show simultaneous targeted mutagenesis of three genes, tyrosinase paralogues (tyra and tyrb) and enhanced green fluorescent protein (egfp) by injection of two TALENs pairs in transgenic embryos carrying egfp. Consistent with the high frequency of both severe phenotypic traits, albinism and loss of GFP fluorescence, frameshift mutation rates of tyr paralogues and egfp reached 40-80%. Next, we show early introduction of TALEN-mediated mutagenesis of these target loci during embryogenesis. Finally, we also demonstrate that two different pairs of TALENs can simultaneously introduce mutations to both paralogues encoding histone chaperone with high efficiency. Our results suggest that targeted mutagenesis of multiple genes using TALENs can be applied to analyze the functions of paralogous genes with redundancy in X. laevis.

High efficiency TALENs enable F0 functional analysis by targeted gene disruption in Xenopus laevis embryos

Recently, gene editing with transcription activator-like effector nucleases (TALENs) has been used in the life sciences. TALENs can be easily customized to recognize a specific DNA sequence and efficiently introduce double-strand breaks at the targeted genomic locus. Subsequent non-homologous end-joining repair leads to targeted gene disruption by base insertion, deletion, or both. Here, to readily evaluate the efficacy of TALENs in Xenopus laevis embryos, we performed the targeted gene disruption of tyrosinase (tyr) and pax6 genes that are involved in pigmentation and eye formation, respectively. We constructed TALENs targeting tyr and pax6 and injected their mRNAs into fertilized eggs at the one-cell stage. Expectedly, introduction of tyr TALEN mRNA resulted in drastic loss of pigmentation with high efficiency. Similarly, for pax6, TALENs led to deformed eyes in the injected embryos. We confirmed mutations of the target alleles by restriction enzyme digestion and sequence analyses of genomic PCR products. Surprisingly, not only biallelic but also paralogous, gene disruption was observed. Our results demonstrate that targeted gene disruption by TALENs provides a method comparable to antisense morpholinos in analyzing gene function in Xenopus F0 embryos, but also applies beyond embryogenesis to any life stage.

Spatiotemporal expression profile of no29/nucleophosmin3 in the intestine of Xenopus laevis during metamorphosis

A Xenopus laevis homolog of nucleophosmin/nucleoplasmin3 (NPM3), no29, has been previously identified as a thyroid hormone (TH)-response gene during TH-induced metamorphosis. X. laevis has another NPM3 homolog (npm3) in the pseudo-tetraploid genome, whereas X. tropicalis possesses one ortholog in the diploid genome. To assess the possible roles of these NPM3 homologs in amphibian metamorphosis, we have analyzed their expression profiles in X. laevis tadpoles. Levels of no29 and npm3 mRNA are rapidly up-regulated by exogenous TH in various organs of the premetamorphic tadpoles. Notably, in the small intestine, no29 and npm3 mRNA levels are transiently up-regulated during metamorphic climax, when progenitor/stem cells of the adult epithelium appear and actively proliferate. In situ hybridization analysis has revealed that the no29 transcript is specifically localized in adult epithelial progenitor/stem cells of the intestine during natural and TH-induced metamorphosis. Double-staining for in situ hybridization and immunohistochemistry has shown co-expression of no29 mRNA and no38 protein (an ortholog of NPM1), which is known to interact with NPM3 and to regulate cell proliferation in mammals. Thus, no29/npm3 might serve as a stem cell marker in the intestine during metamorphosis.