Development of a heat-stable alkaline phosphatase reporter system for cis-regulatory analysis and its application to 3D digital imaging of Xenopus embryonic tissues.

Xenopus is one of the essential model systems for studying vertebrate development. However, one drawback of this system is that, because of the opacity of Xenopus embryos, 3D imaging analysis is limited to surface structures, explant cultures, and post-embryonic tadpoles. To develop a technique for 3D tissue/organ imaging in whole Xenopus embryos, we identified optimal conditions for using placental alkaline phosphatase (PLAP) as a transgenic reporter and applied it to the correlative light microscopy and block-face imaging (CoMBI) method for visualization of PLAP-expressing tissues/organs. In embryos whose endogenous alkaline phosphatase activities were heat-inactivated, PLAP staining visualized various tissue-specific enhancer/promoter activities in a manner consistent with green fluorescent protein (GFP) fluorescence. Furthermore, PLAP staining appeared to be more sensitive than GFP fluorescence as a reporter, and the resulting expression patterns were not mosaic, in striking contrast to the mosaic staining pattern of ß-galactosidase expressed from the lacZ gene that was introduced by the same transgenesis method. Owing to efficient penetration of alkaline phosphatase substrates, PLAP activity was detected in deep tissues, such as the developing brain, spinal cord, heart, and somites, by whole-mount staining. The stained embryos were analyzed by the CoMBI method, resulting in the digital reconstruction of 3D images of the PLAP-expressing tissues. These results demonstrate the efficacy of the PLAP reporter system for detecting enhancer/promoter activities driving deep tissue expression and its combination with the CoMBI method as a powerful approach for 3D digital imaging analysis of specific tissue/organ structures in Xenopus embryos.

Evolutionary Tuning of Transient Receptor Potential Ankyrin 1 Underlies the Variation in Heat Avoidance Behaviors among Frog Species Inhabiting Diverse Thermal Niches.

Environmental temperature is a critical factor for all forms of life, and thermal tolerance defines the habitats utilized by a species. Moreover, the evolutionary tuning of thermal perception can also play a key role in habitat selection. Yet, the relative importance of thermal tolerance and perception in environmental adaptation remains poorly understood. Thermal conditions experienced by anuran tadpoles differ among species due to the variation in breeding seasons and water environments selected by parental frogs. In the present study, heat tolerance and avoidance temperatures were compared in tadpoles from five anuran species that spatially and temporally inhabit different thermal niches. These two parameters were positively correlated with each other and were consistent with the thermal conditions of habitats. The species difference in avoidance temperature was 2.6 times larger than that in heat tolerance, suggesting the importance of heat avoidance responses in habitat selection. In addition, the avoidance temperature increased after warm acclimation, especially in the species frequently exposed to heat in their habitats. Characterization of the heat-sensing transient receptor potential ankyrin 1 (TRPA1) ion channel revealed an amphibian-specific alternatively spliced variant containing a single valine insertion relative to the canonical alternative spliced variant of TRPA1, and this novel variant altered the response to thermal stimuli. The two alternatively spliced variants of TRPA1 exhibited different thermal responses in a species-specific manner, which are likely to be associated with a difference in avoidance temperatures among species. Together, our findings suggest that the functional change in TRPA1 plays a crucial role in thermal adaptation processes.

Phenotype-genotype relationships in Xenopus sox9 crispants provide insights into campomelic dysplasia and vertebrate jaw evolution.

Since CRISPR-based genome editing technology works effectively in the diploid frog Xenopus tropicalis, a growing number of studies have successfully modeled human genetic diseases in this species. However, most of their targets were limited to non-syndromic diseases that exhibit abnormalities in a small fraction of tissues or organs in the body. This is likely because of the complexity of interpreting the phenotypic variations resulting from somatic mosaic mutations generated in the founder animals (crispants). In this study, we attempted to model the syndromic disease campomelic dysplasia (CD) by generating sox9 crispants in X. tropicalis. The resulting crispants failed to form neural crest cells at neurula stages and exhibited various combinations of jaw, gill, ear, heart, and gut defects at tadpole stages, recapitulating part of the syndromic phenotype of CD patients. Genotyping of the crispants with a variety of allelic series of mutations suggested that the heart and gut defects depend primarily on frame-shift mutations expected to be null, whereas the jaw, gill, and ear defects could be induced not only by such mutations but also by in-frame deletion mutations expected to delete part of the jawed vertebrate-specific domain from the encoded Sox9 protein. These results demonstrate that Xenopus crispants are useful for investigating the phenotype-genotype relationships behind syndromic diseases and examining the tissue-specific role of each functional domain within a single protein, providing novel insights into vertebrate jaw evolution.

Cryo-injury procedure-induced cardiac regeneration shows unique gene expression profiles in the newt Pleurodeles waltl.

BACKGROUND: Cardiac regeneration in the adult mouse is not substantial. Some vertebrates, such as newts and zebrafish, regenerate the heart throughout their lives. To understand how regenerative abilities differ among animal species, comparative research has been conducted in animals like mouse, zebrafish, and newt. For those purposes, cryo-injury is suitable as an experimental model for the pathological condition of human myocardial infarction. In fact, cryo-injury procedures are common in mouse and zebrafish. RESULTS: In the present study, we induced cryo-damage on the ventricle in Iberian ribbed newts using a liquid nitrogen-chilled probe. We observed that the injured area recovered within 8 weeks, with remodeling of scar tissue and proliferation of cardiomyocytes. We investigated the subsequent recovery of cryo-injured and amputated tissues by comparative analysis of the gene expression profiles following these two procedures. CONCLUSIONS: Notably, we established a cryo-injury procedure for the newt and confirmed that regeneration of the cryo-damaged myocardial tissue is achieved by changes in gene expression that are milder than those observed in the amputation model. Our results suggest that the cryo-injury method is suitable for comparing the process of cardiac regeneration in the newt with that in other animal models.

Optimization of CRISPR/Cas9-mediated gene disruption in Xenopus laevis using a phenotypic image analysis technique.

The CRISPR/Cas9 method has become popular for gene disruption experiments in Xenopus laevis. However, the experimental conditions that influence the efficiency of CRISPR/Cas9 remain unclear. To that end, we developed an image analysis technique for the semi-quantitative evaluation of the pigment phenotype resulting from the disruption of tyrosinase genes in X. laevis using a CRISPR/Cas9 approach, and then examined the effects of varying five experimental parameters (timing of the CRISPR reagent injection into developing embryos; amount of Cas9 mRNA in the injection reagent; total injection volume per embryo; number of injection sites per embryo; and the culture temperature of the injected embryos) on the gene disruption efficiency. The results of this systematic analysis suggest that the highest possible efficiency of target gene disruption can be achieved by injecting a total of 20 nL of the CRISPR reagent containing 1500 pg of Cas9 mRNA or 4 ng of Cas9 protein into two separate locations (10 nL each) of one-cell stage embryos cultured at 22°C. This study also highlights the importance of balancing the experimental parameters for increasing gene disruption efficiency and provides valuable insights into the optimal conditions for applying the CRISPR/Cas9 system to new experimental organisms.

Heterologous expression and biochemical comparison of two homologous SoxX proteins of endosymbiotic Candidatus Vesicomyosocius okutanii and free-living Hydrogenovibrio crunogenus from deep-sea environments.

Candidatus Vesicomyosocius okutanii is a currently uncultured endosymbiotic bacterium of Phreagena okutanii, a clam that inhabits deep-sea vent environments. The genome of Ca. V. okutanii encodes a sulfur-oxidizing (Sox) enzyme complex, presumably generating biological energy for the host from inorganic sulfur compounds. Here, Ca. V. okutanii SoxX (VoSoxX), a mono-heme cytochrome c component of the Sox complex, was shown to be phylogenetically related to its homologous counterpart (HcSoxX) from a free-living deep-sea bacterium, Hydrogenovibrio crunogenus. Both proteins were heterologously expressed in Escherichia coli co-expressing cytochrome c maturation genes for comparative biochemical analysis. The VoSoxX recombinant had significantly lower thermal stability than HcSoxX, reflecting the difference in growth conditions of the source bacteria. The endosymbiont inhabits a mild intracellular environment, whereas the free-living bacterium dwells in a harsh environment. This study represents the first successful case of heterologous expression of genes from Ca. V. okutanii, allowing further biochemical studies of the molecular mechanism of sulfur oxidation in deep-sea environments.

Transcriptome analyses of immune tissues from three Japanese frogs (genus Rana ) reveals their utility in characterizing major histocompatibility complex class II.

BACKGROUND: In Japan and East Asia, endemic frogs appear to be tolerant or not susceptible to chytridiomycosis, a deadly amphibian disease caused by the chytrid fungus Batrachochytridium dendrobatidis (Bd). Japanese frogs may have evolved mechanisms of immune resistance to pathogens such as Bd. This study characterizes immune genes expressed in various tissues of healthy Japanese Rana frogs. RESULTS: We generated transcriptome data sets of skin, spleen and blood from three adult Japanese Ranidae frogs (Japanese brown frog Rana japonica, the montane brown frog Rana ornativentris, and Tago's brown frog Rana tagoi tagoi) as well as whole body of R. japonica and R. ornativentris tadpoles. From this, we identified tissue- and stage-specific differentially expressed genes; in particular, the spleen was most enriched for immune-related genes. A specific immune gene, major histocompatibility complex class IIB (MHC-IIB), was further characterized due to its role in pathogen recognition. We identified a total of 33 MHC-IIB variants from the three focal species (n = 7 individuals each), which displayed evolutionary signatures related to increased MHC variation, including balancing selection. Our supertyping analyses of MHC-IIB variants from Japanese frogs and previously studied frog species identified potential physiochemical properties of MHC-II that may be important for recognizing and binding chytrid-related antigens. CONCLUSIONS: This is one of the first studies to generate transcriptomic resources for Japanese frogs, and contributes to further understanding the immunogenetic factors associated with resistance to infectious diseases in amphibians such as chytridiomycosis. Notably, MHC-IIB supertyping analyses identified unique functional properties of specific MHC-IIB alleles that may partially contribute to Bd resistance, and such properties provide a springboard for future experimental validation.

Characterisation of major histocompatibility complex class I genes in Japanese Ranidae frogs.

The major histocompatibility complex (MHC) is a key component of adaptive immunity in all jawed vertebrates, and understanding the evolutionary mechanisms that have shaped these genes in amphibians, one of the earliest terrestrial tetrapods, is important. We characterised MHC class I variation in three common Japanese Rana species (Rana japonica, Rana ornativentris and Rana tagoi tagoi) and identified a total of 60 variants from 21 individuals. We also found evolutionary signatures of gene duplication, recombination and balancing selection (including trans-species polymorphism), all of which drive increased MHC diversity. A unique feature of MHC class I from these three Ranidae species includes low synonymous differences per site (d S) within species, which we attribute to a more recent diversification of these sequences or recent gene duplication. The resulting higher d N/d S ratio relative to other anurans studied could be related to stronger selection pressure at peptide binding sites. This is one of the first studies to investigate MHC in Japanese amphibians and permits further exploration of the polygenetic factors associated with resistance to infectious diseases.

Extraordinary Diversity in the Origins of Sex Chromosomes in Anurans Inferred from Comparative Gene Mapping.

Sex determination in frogs (anurans) is genetic and includes both male and female heterogamety. However, the origins of the sex chromosomes and their differentiation processes are poorly known. To investigate diversity in the origins of anuran sex chromosomes, we compared the chromosomal locations of sex-linked genes in 4 species: the African clawed frog (Xenopus laevis), the Western clawed frog (Silurana/X. tropicalis), the Japanese bell-ring frog (Buergeria buergeri), and the Japanese wrinkled frog (Rana rugosa). Comparative mapping data revealed that the sex chromosomes of X. laevis, X. tropicalis and R. rugosa are different chromosome pairs; however, the sex chromosomes of X. tropicalis and B. buergeri are homologous, although this may represent distinct evolutionary origins. We also examined the status of sex chromosomal differentiation in B. buergeri, which possesses heteromorphic ZW sex chromosomes, using comparative genomic hybridization and chromosome painting with DNA probes from the microdissected W chromosome. At least 3 rearrangement events have occurred in the proto-W chromosome: deletion of the nucleolus organizer region and a paracentric inversion followed by amplification of non-W-specific repetitive sequences.

Microsatellite marker development by multiplex ion torrent PGM sequencing: a case study of the endangered Odorrana narina complex of frogs.

The endangered Ryukyu tip-nosed frog Odorrana narina and its related species, Odorrana amamiensis, Odorrana supranarina, and Odorrana utsunomiyaorum, belong to the family Ranidae and are endemically distributed in Okinawa (O. narina), Amami and Tokunoshima (O. amamiensis), and Ishigaki and Iriomote (O. supranarina and O. utsunomiyaorum) Islands. Because of varying distribution patterns, this species complex is an intrinsic model for speciation and adaptation. For effective conservation and molecular ecological studies, further genetic information is needed. For rapid, cost-effective development of several microsatellite markers for these and 2 other species, we used next-generation sequencing technology of Ion Torrent PGM™. Distribution patterns of repeat motifs of microsatellite loci in these modern frog species (Neobatrachia) were similarly skewed. We isolated and characterized 20 new microsatellite loci of O. narina and validated cross-amplification in the three-related species. Seventeen, 16, and 13 loci were cross-amplified in O. amamiensis, O. supranarina, and O. utsunomiyaorum, respectively, reflecting close genetic relationships between them. Mean number of alleles and expected heterozygosity of newly isolated loci varied depending on the size of each inhabited island. Our findings suggested the suitability of Ion Torrent PGM™ for microsatellite marker development. The new markers developed for the O. narina complex will be applicable in conservation genetics and molecular ecological studies.

Sex chromosome cycle as a mechanism of stable sex determination.

Recent advances in DNA sequencing technology have enabled the precise decoding of genomes in non-model organisms, providing a basis for unraveling the patterns and mechanisms of sex chromosome evolution. Studies of different species have yielded conflicting results regarding the traditional theory that sex chromosomes evolve from autosomes via the accumulation of deleterious mutations and degeneration of the Y (or W) chromosome. The concept of the 'sex chromosome cycle,' emerging from this context, posits that at any stage of the cycle (i.e., differentiation, degeneration, or loss), sex chromosome turnover can occur while maintaining stable sex determination. Thus, understanding the mechanisms that drive both the persistence and turnover of sex chromosomes at each stage of the cycle is crucial. In this review, we integrate recent findings on the mechanisms underlying maintenance and turnover, with a special focus on several organisms having unique sex chromosomes. Our review suggests that the diversity of sex chromosomes in the maintenance of stable sex determination is underappreciated and emphasizes the need for more research on the sex chromosome cycle.

Stable aquaculture of the Japanese lancelet Branchiostoma japonicum for 7 years.

Despite advances in the study on animal evolution in the last two decades, paucity of experimental data on cephalochordates comparable to those on the other chordates hinders an integrative understanding of chordate evolutionary history. To obtain lancelet data under well-controlled experiments, laboratory cultures of lancelets have been performed at several institutions. In a mass culture started in 2005, we have obtained up to three consecutive generations of Branchiostoma japonicum. Using sand substratum, survival rates of laboratory lancelets until maturation have improved to greater than 30%, much higher than compared to previously, and for adults the annual average survival rate was 82.3%. The high survival rate allows maintaining animals at least 6 years and potentially longer. Water temperatures lower than 23°C obviously reduced the frequency of spawning even after the onset of spawning period, and 1-2 days after changing the temperature at 25°C animals became spawned well. We also observed obvious sex reversal from male to female in individuals that had been cultured for 3 years or more. Our continuous culture has provided sufficient materials for vital experiments on early development and for studying metamorphosis, as well as for the conservation of wild populations. The subculture of successive laboratory generations will provide a valuable resource for genetic studies.

Factors associated with spinal instability in low back lumbar diseases with leg pain: Analysis of sagittal translation and segmental angulation.

BACKGROUND: Determining the association between radiographic spinal instability assessment and lower back lumbar diseases with lower limb symptoms can contribute to evidence-based assessment and treatment in clinical practice and rehabilitation. Therefore, radiological evidence of lumbar spine instability assessment, such as sagittal translation (ST) and segmental angulation (SA), is clinically important. OBJECTIVE: To identify factors associated with the assessment of spinal instability in lumbar disc herniation with leg pain and discogenic low back pain using ST and SA. METHODS: We examined 112 patients with lumbar disc herniation with leg pain and 116 with discogenic low back pain at our clinic from 2016 to 2021. Data on age, gender, sports activities, and occupation were collected from medical records. Additionally, ST and SA of L4 and L5 during maximum trunk flexion and extension were measured using radiography. Simple and multiple logistic regression analyses were used for statistical analysis. RESULTS: Simple logistic regression analysis showed that ST and SA (odds ratio [OR]: 1.11; 95% confidence interval [CI]: 1.03-1.19) were associated with lumbar disc herniation. Multiple logistic regression analysis showed that only ST was associated with lumbar disc herniation (OR: 2.29; 95% CI: 1.78-3.00). CONCLUSION: Multiple logistic regression analysis showed that ST was associated with lumbar disc herniation with leg pain and had a stronger association than SA.

Rapid and collective determination of the complete "hot-spring frog" mitochondrial genome containing long repeat regions using Nanopore sequencing.

The mitochondrial genome (mt-genome) is one of the promising molecular markers for phylogenetics and population genetics. Recently, various mt-genomes have been determined rapidly by using massively parallel sequencers. However, the control region (CR, also called D-loop) in mt-genomes remain difficult to precisely determine due to the presence of repeat regions. Here, using Nanopore sequencing, we succeeded in rapid and collective determination of complete mt-genome of the hot-spring frog, Buergeria japonica, and found that its mt-genome size was 22,274 bp including CR (6,929 bp) with two types of tandem repeat motifs forming repeat regions. Comparison of assembly strategies revealed that the long- and short-read data combined together enabled efficient determination of the CR, but the short-read data alone did not. The B. japonica CR was longer than that of a congenic species inhabiting cooler climate areas, Buergeria buergeri, because of the long repeat regions in the former. During the thermal adaptation of B. japonica, the longer repeat regions in its CR may have accumulated within a period after divergence from B. buergeri.

Identification of tumor-related genes via RNA sequencing of tumor tissues in Xenopus tropicalis.

Cancer treatment is still challenging because the disease is often caused by multiple mutations. Although genomic studies have identified many oncogenes and tumor suppressor genes, gene sets involved in tumorigenesis remain poorly understood. Xenopus, a genus of aquatic frogs, is a useful model to identify gene sets because it can be genetically and experimentally analyzed. Here, we analyzed gene expression in tumor tissues of three individuals in Xenopus tropicalis and identified 55 differentially expressed genes (DEGs). Gene ontology (GO) analysis showed that the upregulated genes in the tumor tissues were enriched in GO terms related to the extracellular matrix and collagen fibril organization. Hierarchical clustering showed that the gene expression patterns of tumor tissues in X. tropicalis were comparable to those of human connective, soft, and subcutaneous tissue-derived cancers. Additionally, pathway analysis revealed that these DEGs were associated with multiple pathways, including the extracellular matrix, collagen fibril organization, MET signaling, and keratan sulfate. We also found that the expression tendency of some DEGs that have not been well analyzed in the cancer field clearly determines the prognosis of human cancer patients. This study provides a remarkable reference for future experimental work on X. tropicalis to identify gene sets involved in human cancer.

Conserved Evolution of MHC Supertypes among Japanese Frogs Suggests Selection for Bd Resistance.

The chytrid fungus Batrachochytrium dendrobatidis (Bd) is a major threat to amphibians, yet there are no reports of major disease impacts in East Asian frogs. Genetic variation of the major histocompatibility complex (MHC) has been associated with resistance to Bd in frogs from East Asia and worldwide. Using transcriptomic data collated from 11 Japanese frog species (one individual per species), we isolated MHC class I and IIb sequences and validated using molecular cloning. We then compared MHC from Japanese frogs and other species worldwide, with varying Bd susceptibility. Supertyping analysis, which groups MHC alleles based on physicochemical properties of peptide binding sites, identified that all examined East Asian frogs contained at least one MHC-IIb allele belonging to supertype ST-1. This indicates that, despite the large divergence times between some Japanese frogs (up to 145 million years), particular functional properties in the peptide binding sites of MHC-II are conserved among East Asian frogs. Furthermore, preliminary analysis using NetMHCIIpan-4.0, which predicts potential Bd-peptide binding ability, suggests that MHC-IIb ST-1 and ST-2 have higher overall peptide binding ability than other supertypes, irrespective of whether the peptides are derived from Bd, other fungi, or bacteria. Our findings suggest that MHC-IIb among East Asian frogs may have co-evolved under the same selective pressure. Given that Bd originated in this region, it may be a major driver of MHC evolution in East Asian frogs.