Genome evolution in the allotetraploid frog Xenopus laevis

Xenopus laevis as an infection model for human pathogenic bacteria

Animal infection models are essential for understanding bacterial pathogenicity and corresponding host immune responses. In this study, we investigated whether juvenile Xenopus laevis could be used as an infection model for human pathogenic bacteria. Xenopus frogs succumbed to intraperitoneal injection containing the human pathogenic bacteria Staphylococcus aureus, Pseudomonas aeruginosa, and Listeria monocytogenes. In contrast, non-pathogenic bacteria Bacillus subtilis and Escherichia coli did not induce mortality in Xenopus frogs. The administration of appropriate antibiotics suppressed mortality caused by S. aureus and P. aeruginosa. Strains lacking the agr locus, cvfA (rny) gene, or hemolysin genes in S. aureus, LIPI-1-deleted mutant of L. monocytogenes, which attenuate virulence within mammals, exhibited reduced virulence in Xenopus frogs compared with their respective wild-type counterparts. Bacterial distribution analysis revealed that S. aureus persisted in the blood, liver, heart, and muscles of Xenopus frogs until death. These results suggested that intraperitoneal injection of human pathogenic bacteria induces sepsis-like symptoms in Xenopus frogs, supporting their use as a valuable animal model for evaluating antimicrobial efficacy and identifying virulence genes in various human pathogenic bacteria.

Why some hearts heal and others don't: The phylogenetic landscape of cardiac regenerative capacity

The limited ability of adult humans to replenish lost heart muscle cells after a heart attack has attracted scientists to explore natural heart regeneration capabilities in the animal kingdom. In particular, research has accelerated since the landmark discovery more than twenty years ago that zebrafish can completely regrow myocardial tissue. In this review, we survey heart regeneration studies in diverse model and non-model animals, aiming to gain insights into both the evolutionary trends in cardiac regenerative potential and the variations among closely related species. Differences in cardiomyogenesis, vasculature formation, and the communication between cardiovascular cells and other players have been investigated to understand the cellular basis, although the precise molecular and genetic causes underlying the stark differences in cardiac regenerative potential among certain close cousins remain largely unknown. By studying cardiovascular regeneration and repair in diverse organisms, we may uncover distinct mechanisms, offering new perspectives for advancing regenerative medicine.

Exploring genome architecture as a source of phylogenetic characters for resolving the apulmonate arachnid polytomy

Chromosome-level genome assemblies are powerful tools for identifying the presence of rare genomic changes that can overcome phylogenetically intractable problems. Chelicerata, the sister group to the remaining arthropods, harbors a soft polytomy at the base of an internal node named Euchelicerata, which is variably resolved across phylogenomic studies. As a result, seven orders, comprising horseshoe crabs and six apulmonate arachnid lineages, exhibit highly unstable placements from one study to the next, typically with maximal nodal support. Here, we analyzed recently released chromosome-level genomes of two of these orders, Opiliones (harvestmen) and Solifugae (camel spiders). We show that both Opiliones and Solifugae exhibit an unduplicated genome condition, as inferred from analysis of gene clusters, microRNAs, and macrosynteny. These results are congruent with phylogenomic studies that have refuted traditional morphological placements of Opiliones and Solifugae as close relatives of orders within Arachnopulmonata, a subset of six arachnid orders that are united by a shared whole genome duplication. Additionally, we examine irreversible chromosome fusion-with-mixing events as potential sources of phylogenetic data. We show that while fusion and mixing events are common in apulmonate arachnids, multiple mixing events support incompatible unrooted tree topologies. These results suggest that fusion and mixing events have evolved convergently in the chelicerate tree of life, particularly for extant lineages with a small number of chromosomes. Overall, our findings demonstrate that broader sampling of chelicerate genomes and establishment of genomic resources for key missing orders are essential to unlocking the potential of rare genomic changes as phylogenetic data sources.

Probiotic colonization of Xenopus laevis skin causes short-term changes in skin microbiomes and gene expression

Probiotic therapies have been suggested for amelioration efforts of wildlife disease such as chytridiomycosis caused by Batrachochytrium spp. in amphibians. However, there is a lack of information on how probiotic application affects resident microbial communities and immune responses. To better understand these interactions, we hypothesized that probiotic application would alter microbial community composition and host immune expression in Xenopus laevis. Accordingly, we applied three amphibian-derived and anti-Batrachochytrium bacteria strains (two Pseudomonas spp. and one Stenotrophomonas sp.) to X. laevis in monoculture and also as a cocktail. We quantified microbial community structure using 16S rRNA gene sequencing. We also quantified genes involved in X. laevis immune responses using quantitative reverse transcriptase polymerase chain reaction (RT-qPCR) and skin transcriptomics over 1 and 3-week periods. All probiotic treatments successfully colonized X. laevis skin for 3 weeks, but with differential amplicon sequence variant (ASV) sequence counts over time. Bacterial community and immune gene effects were most pronounced at week 1 post-probiotic exposure and decreased thereafter. All probiotic treatments caused initial changes to bacterial community alpha and beta diversity, including reduction in diversity from pre-exposure anti-Batrachochytrium bacterial ASV relative abundance. Probiotic colonization by Pseudomonas probiotic strain RSB5.4 reduced expression of regulatory T cell marker (FOXP3, measured with RT-qPCR) and caused the greatest gene expression changes detected by transcriptomics. Single bacterial strains and mixed cultures, therefore, altered amphibian microbiome-immune interactions. This work will help to improve our understanding of the role of the microbiome-immune interface underlying both disease dynamics and emergent eco-evolutionary processes.IMPORTANCEAmphibian skin microbial communities have an important role in determining disease outcomes, in part through complex yet poorly understood interactions with host immune systems. Here we report that probiotic-induced changes to the Xenopus laevis frog skin microbial communities also result in significant alterations to these animals' immune gene expression. These findings underscore the interdependence of amphibian skin immune-microbiome interactions.

Genomic Complexity of ccdc40 in Xenopus : Implications for CRISPR Targeting and Disease Modeling

Mutations in CCDC40 cause primary ciliary dyskinesia in humans. To evaluate the pathogenicity of variants in CCDC40 , we examined the genomic structure of this gene in Xenopus tropicalis , a diploid frog suitable as a model for genetic studies. We identified inconsistencies in the current ccdc40 gene model and discovered two distinct ccdc40 genes near the previously annotated locus. Surprisingly, Xenopus laevis , an allotetraploid species that typically has two homoeologs, contains only one homoeolog ( ccdc40.S ), making it a more suitable genetic model for studying ccdc40 function and potentially expediting the functional characterization of CCDC40 variants linked to primary ciliary dyskinesia.

Large tandem repeats of grass frog (Rana temporaria) in silico and in situ

BACKGROUND

Genomes of higher eukaryotes contain a large fraction of non-coding repetitive DNA, including tandem repeats (TRs) and transposable elements (TEs). The impact of TRs on genome structure and function and the importance of TR transcripts have been described for several model species. Amphibians have one of the most diverse genome sizes among vertebrates, attributed to the abundance of repetitive non-coding DNA. Consequently, amphibians are good models for the analysis of repetitive sequences, including TRs. However, few studies have focused on amphibian genomes.

RESULTS

Bioinformatic analyses were performed to characterise the content and localisation of TRs in the sequenced grass frog Rana temporaria genome. By applying different bioinformatic approaches, 76 TR families and 314 single TR arrays (not grouped into families) were identified. Each TR was characterised on the basis of chromosomal position, monomer length and variability and GC content. Bioinformatic analysis revealed a great diversity of TRs, with a clear predominance of TRs with short monomers (< 100 bp), although TRs with long monomers (> 1000 bp) also exist. The six most abundant TRs were successfully mapped by fluorescence in situ hybridization (FISH), which highlighted the presence of specific TR sequences in strategic chromosomal regions, i.e., the pericentromeric regions. A comparison of the results of in situ and in silico TR mapping revealed some inaccuracies in the assembly of heterochromatic regions. A putative new non-autonomous TE called "FEDoR" (Frog Element Dispersed organised Repeat) is also described. FEDoR is ∼ 3.5 kb in length, has no significant similarity to any known TE family, contains multiple internal TR motifs, and is flanked on both sides by pairs of inverted repeat sequences (IRSs) and target site duplications (TSDs).

CONCLUSION

Characterisation of TRs in this frog species has provided some insights regarding TR biology in Anuran amphibians.

Establishing genome sequencing and assembly for non-model and emerging model organisms: a brief guide

Reference genome assemblies are the basis for comprehensive genomic analyses and comparisons. Due to declining sequencing costs and growing computational power, genome projects are now feasible in smaller labs. De novo genome sequencing for non-model or emerging model organisms requires knowledge about genome size and techniques for extracting high molecular weight DNA. Next to quality, the amount of DNA obtained from single individuals is crucial, especially, when dealing with small organisms. While long-read sequencing technologies are the methods of choice for creating high quality genome assemblies, pure short-read assemblies might bear most of the coding parts of a genome but are usually much more fragmented and do not well resolve repeat elements or structural variants. Several genome initiatives produce more and more non-model organism genomes and provide rules for standards in genome sequencing and assembly. However, sometimes the organism of choice is not part of such an initiative or does not meet its standards. Therefore, if the scientific question can be answered with a genome of low contiguity in intergenic parts, missing the high standards of chromosome scale assembly should not prevent publication. This review describes how to set up an animal genome sequencing project in the lab, how to estimate costs and resources, and how to deal with suboptimal conditions. Thus, we aim to suggest optimal strategies for genome sequencing that fulfil the needs according to specific research questions, e.g. "How are species related to each other based on whole genomes?" (phylogenomics), "How do genomes of populations within a species differ?" (population genomics), "Are differences between populations relevant for conservation?" (conservation genomics), "Which selection pressure is acting on certain genes?" (identification of genes under selection), "Did repeats expand or contract recently?" (repeat dynamics).

Impact of a horizontally transferred Helitron family on genome evolution in Xenopus laevis

BACKGROUND

Within eukaryotes, most horizontal transfer of genetic material involves mobile DNA sequences and such events are called horizontal transposable element transfer (HTT). Although thousands of HTT examples have been reported, the transfer mechanisms and their impacts on host genomes remain elusive.

RESULTS

In this work, we carefully annotated three Helitron families within several Xenopus frog genomes. One of the Helitron family, Heli1Xen1, is recurrently involved in capturing and shuffling Xenopus laevis genes required in early embryonic development. Remarkably, we found that Heli1Xen1 is seemingly expressed in X. laevis and has produced multiple genomic polymorphisms within the X. laevis population. To identify the origin of Heli1Xen1, we searched its consensus sequence against available genome assemblies. We found highly similar copies in the genomes of another 13 vertebrate species from divergent vertebrate lineages, including reptiles, ray-finned fishes and amphibians. Further phylogenetic analysis provides evidence showing that Heli1Xen1 invaded these lineages via HTT quite recently, around 0.58-10.74 million years ago.

CONCLUSIONS

The frequently Heli1Xen1-involved HTT events among reptiles, fishes and amphibians could provide insights into possible vectors for transfer, such as shared viruses across lineages. Furthermore, we propose that the Heli1Xen1 sequence could be an ideal candidate for studying the mechanism and genomic impact of Helitron transposition.

Transposable element dynamics in Xenopus laevis embryogenesis: a tale of two coexisting subgenomes

The African clawed frog Xenopus laevis has an allotetraploid genome consisting of two subgenomes referred as L relating to the Long chromosomes and S relating to the Short chromosomes. While the L subgenome presents conserved synteny with X. tropicalis chromosomes, the S subgenome has undergone rearrangements and deletions leading to differences in gene and transposable element (TE) content between the two subgenomes. The asymmetry in the evolution of the two subgenomes is also detectable in gene expression levels and TE mobility. TEs, also known as "jumping genes", are mobile genetic elements having a key role in genome evolution and gene regulation. However, due to their potential deleterious effects, TEs are controlled by host defense mechanisms such as the nucleosome remodeling and deacetylase (NuRD) complex and the Argonaute proteins that mainly modify the heterochromatin environment. In embryogenesis, TEs can escape the silencing mechanisms during the maternal-to-zygotic transition when a transcriptionally permissive environment is created. Moreover, further evidence highlighted that the reactivation of TEs during early developmental stages is not the result of this genome-wide reorganization of chromatin but it is class and stage-specific, suggesting a precise regulation. In line with these premises, we explored the impact of TE transcriptional contribution in six developmental stages of X. laevis. Overall, the expression pattern referred to the entire set of transcribed TEs was constant across the six developmental stages and in line with their abundance in the genome. However, focusing on subgenome-specific TEs, our analyses revealed a distinctive transcriptional pattern dominated by LTR retroelements in the L subgenome and LINE retroelements in the S subgenome attributable to young copies. Interestingly, genes encoding proteins involved in maintaining the repressive chromatin environment were active in both subgenomes highlighting that TE controlling systems were active in X. laevis embryogenesis and evolved symmetrically.

Sperm derived H2AK119ub1 is required for embryonic development in Xenopus laevis

Ubiquitylation of H2A (H2AK119ub1) by the polycomb repressive complexe-1 plays a key role in the initiation of facultative heterochromatin formation in somatic cells. Here we evaluate the contribution of sperm derived H2AK119ub1 to embryo development. In Xenopus laevis we found that H2AK119ub1 is present during spermiogenesis and into early embryonic development, highlighting its credential for a role in the transmission of epigenetic information from the sperm to the embryo. In vitro treatment of sperm with USP21, a H2AK119ub1 deubiquitylase, just prior to injection to egg, results in developmental defects associated with gene upregulation. Sperm H2AK119ub1 editing disrupts egg factor mediated paternal chromatin remodelling processes. It leads to post-replication accumulation of H2AK119ub1 on repeat element of the genome instead of CpG islands. This shift in post-replication H2AK119ub1 distribution triggered by sperm epigenome editing entails a loss of H2AK119ub1 from genes misregulated in embryos derived from USP21 treated sperm. We conclude that sperm derived H2AK119ub1 instructs egg factor mediated epigenetic remodelling of paternal chromatin and is required for embryonic development.

Genomic Insights into Fertilization: Tracing PLCZ1 Orthologs Across Amphibian Lineages

Fertilization triggers a cascade of events, including a rise in egg cytosolic calcium that marks the onset of embryonic development. In mammals and birds, this critical process is mediated by the sperm-derived phospholipase C zeta (PLCζ), which is pivotal in releasing calcium from the endoplasmic reticulum in the egg and initiating embryonic activation. Intriguingly, Xenopus laevis, a key model organism in reproductive biology, lacks an annotated PLCZ1 gene, prompting questions about its calcium release mechanism during fertilization. Using bioinformatics and RNA sequencing of adult X. laevis testes, we investigated the presence of a PLCZ1 ortholog in amphibians. While we identified PLCZ1 homologs in 25 amphibian species, including 14 previously uncharacterized orthologs, we found none in X. laevis or its close relative, Xenopus tropicalis. Additionally, we found no compensatory expression of other PLC isoforms in these species. Synteny analysis revealed a PLCZ1 deletion in species within the Pipidae family and another intriguing deletion of potential sperm factor PLCD4 in the mountain slow frog, Nanorana parkeri. Our findings indicate that the calcium release mechanism in frog eggs involves a signaling pathway distinct from the PLCζ-mediated process observed in mammals.

FISH mapping in Xenopus pygmaeus refines understanding of genomic rearrangements and reveals jumping NORs in African clawed frogs

Chromosomal rearrangements are fundamental evolutionary drivers leading to genomic diversification. African clawed frogs (genus Xenopus, subgenera Silurana and Xenopus) represent an allopolyploid model system with conserved chromosome numbers in species with the same ploidy within each subgenus. Two significant interchromosomal rearrangements have been identified: a translocation between chromosomes 9 and 2, found in subgenus Silurana, and a fusion between chromosomes 9 and 10, probably widespread in subgenus Xenopus. Here, we study the allotetraploid Xenopus pygmaeus (subgenus Xenopus) based on in-depth karyotype analysis using chromosome measurements and fluorescent in situ hybridization (FISH). We designed FISH probes for genes associated with translocation and fusion to test for the presence of the two main types of rearrangements. We also examined the locations of 5S and 28S ribosomal tandem repeats, with the former often associated with telomeric regions and the latter with nucleolus organizer regions (NORs). The translocation-associated gene mapping did not detect the translocation in X. pygmaeus, supporting the hypothesis that the translocation is restricted to Silurana, but instead identified a pericentromeric inversion on chromosome 2S. The fusion-associated gene mapping confirmed the fusion of chromosomes 9 and 10, supporting this fusion as an ancestral state in subgenus Xenopus. As expected, the 5S repeats were found predominantly in telomere regions on almost all chromosomes. The nucleolar 28S repeats were localized on chromosome 6S, a position previously found only in the closely related species X. parafraseri, whereas other, phylogenetically more distant species have NORs located on different chromosomes. We therefore hypothesize that a jumping mechanism could explain the relatively frequent changes in the location of NORs during Xenopus evolution.

Evolution and function of galectins in Xenopus laevis: Comparison with mammals and new perspectives

Galectins are metal-independent sugar-binding proteins that recognize galactose (the β-galactoside structure) and regulate the cross-linking of sugar chains between cells and the extracellular matrix. Their specificity for galactose is attributed to their highly conserved carbohydrate recognition domain. Galectins participate in biological processes across species, including development, differentiation, morphogenesis, tumor progression, metastasis, immunity, and apoptosis. However, the relationship between the binding of galectin to sugar chains (glycans) and their biological functions remains unclear. Thus, a comprehensive functional analysis of galectins is required to better characterize their evolutionarily conserved and unique functions. We have previously identified and characterized 12 Xenopus laevis galectins (xgalectins), the only non-mammalian vertebrate species in which galectins have been comprehensively characterized to date. In this review, we present the latest findings on the types and functions of xgalectins and discuss prospects for elucidating their diverse functions from an evolutionary perspective through comparisons with mammalian galectins.

The Osgin Gene Family: Underexplored Yet Essential Mediators of Oxidative Stress

The Osgin gene family consists of two members, Osgin1 and Osgin2, involved in the cellular oxidative stress response. While many members of this essential cellular pathway have been extensively characterized, the Osgin gene family, despite its broad phylogenetic distribution, has received far less attention. Here, we review published articles and open-source databases to synthesize the current research on the evolutionary history, structure, biochemical and physiological functions, expression patterns, and role in disease of the Osgin gene family. Although Osgin displays broad spatiotemporal expression during development and adulthood, there is ambiguity regarding the cellular functions of the OSGIN proteins. A recent study identified OSGIN-1 as a flavin-dependent monooxygenase, but the biochemical role of OSGIN-2 has not yet been defined. Moreover, while the Osgin genes are implicated as mediators of cell proliferation, apoptosis, and autophagy, these functions have not been connected to the enzymatic classification of OSGIN. Misregulation of Osgin expression has long been associated with various disease states, yet recent analyses highlight the mechanistic role of OSGIN in pathogenesis and disease progression, underscoring the therapeutic potential of targeting OSGIN. In light of these findings, we suggest further avenues of research to advance our understanding of this essential, yet underexplored, gene family.

An Extensive Survey of Vertebrate-specific, Nonvisual Opsins Identifies a Novel Subfamily, Q113-Bistable Opsin

A group of nonvisual opsins specific to vertebrates is essential to understand evolution of lateral eyes, one of the most prominent innovations in this lineage. Nevertheless, our knowledge of their evolutionary history remains limited. To develop an integrated view of their evolution, we surveyed these non-visual opsins (VA opsin, pinopsin, parapinopsin, parietopsin, and parapinopsin-like) in 451 vertebrate genomes. Through extensive manual curation, we completed a high-quality catalog. We could not find them in 202 mammals, supporting previous reports of their loss. VA opsins are highly conserved among nonmammals. In contrast, other opsin subfamilies experienced more dynamic molecular evolution with many secondary losses. In addition, we found a previously unreported opsin subfamily that we named Q113-Bistable (QB) opsin. We found its orthologs only in several lizards and the tuatara. Nevertheless, QB opsin pseudogenes were discovered in diverse taxa, including ray-finned fishes, indicating its ancient origin. QB opsin, parapinopsin, and parietopsin are extremely prone to be lost in the course of evolution, and loss events involving these opsins seem to occur concomitantly. Furthermore, we demonstrated the spectral properties of QB opsin as a UV-sensitive, bistable photo-pigment. This study provides the first integrated view of the entire evolutionary history of this group of opsins.

k-mer approaches for biodiversity genomics

The wide array of currently available genomes displays a wonderful diversity in size, composition, and structure and is quickly expanding thanks to several global biodiversity genomics initiatives. However, sequencing of genomes, even with the latest technologies, can still be challenging for both technical (e.g., small physical size, contaminated samples, or access to appropriate sequencing platforms) and biological reasons (e.g., germline-restricted DNA, variable ploidy levels, sex chromosomes, or very large genomes). In recent years, k-mer-based techniques have become popular to overcome some of these challenges. They are based on the simple process of dividing the analyzed sequences (e.g., raw reads or genomes) into a set of subsequences of length k, called k-mers, and then analyzing the frequency or sequences of those k-mers. Analyses based on k-mers allow for a rapid and intuitive assessment of complex sequencing data sets. Here, we provide a comprehensive review to the theoretical properties and practical applications of k-mers in biodiversity genomics with a special focus on genome modeling.

Fertilization and the fast block to polyspermy in the African Clawed Frog, Xenopus laevis: A historical perspective

The African clawed frog, Xenopus laevis, has long been a model organism for studying fertilization due to its large and abundant eggs that are easily manipulated and rapidly undergo embryonic development. Research on this model organism has provided significant insights into the mechanisms that ensure successful fertilization, including the prevention of polyspermy. Polyspermy, the fertilization of an egg by multiple sperm, poses a significant threat to successful embryonic development in most sexually reproducing animals. To counter this, eggs have evolved mechanisms known as polyspermy blocks, which prevent additional sperm from entering once fertilization has occurred. This review focuses on fertilization research in general, and specifically on studies of the fast block to polyspermy in X. laevis. We trace key discoveries and experimental advancements that have shaped our current understanding. Indeed, studies on X. laevis have revealed that fertilization triggers a depolarization of the egg membrane mediated by an efflux of Cl- through the Ca2+-activated Cl- channel TMEM16A, effectively preventing polyspermy. Despite these advances, several questions remain regarding the precise molecular interactions and signaling pathways involved. Continued research on X. laevis promises to uncover further details about the earliest events in embryogenesis and the voltage-dependent mechanisms of fertilization, offering broader insights into reproductive biology across species.

Every Gain Comes With Loss: Ecological and Physiological Shifts Associated With Polyploidization in a Pygmy Frog

Polyploidization plays a pivotal role in vertebrate evolution and diversification. However, the effects of polyploidization on animals across various biological levels, and how these differences drive ecological shifts, remain unclear. Through karyotype analysis and whole-genome sequencing, we identified an autotetraploid Microhyla fissipes from Hainan Island, which shows reproductive isolation and geographic differentiation from its diploid counterpart. Tetraploids exhibited larger cell size, improved tadpole growth rates, and greater whole-body size, along with reduced cell cycle activity. Rather than being simple scaled-up diploids, tetraploids showed shifts in physiological performance, organ allometry, gene expression profiles, and metabolic patterns. Tetraploid adults demonstrated superior jumping ability and increased reproductive investment (e.g. larger gonads and steeper slopes in the relationship between gonadal weight and body weight), suggesting a potential competitive advantage over diploids. However, tetraploids exhibited higher energy expenditure at elevated temperatures, reduced hepatic energy storage, and altered pulmonary regulatory metabolites at 25 °C. Males had smaller relative heart sizes, and females showed flatter slopes in the relationship between heart and lung weight and body weight, indicating reduced investment in cardiopulmonary system. These variations suggest an increased risk of metabolic constraints under heat stress, putting tetraploids at a disadvantage in warmer regions. Importantly, the physiological tradeoffs associated with polyploidization help explain the geographical differentiation between diploids and tetraploids, which reflects a climatic boundary, with tetraploids occupying cooler northeastern areas. Our findings identify an autotetraploid frog, report the first autotetraploid genome in amphibians, and demonstrate how vertebrate polyploids physiologically and ecologically diverge from their diploid counterparts.

Mechanistic basis of temperature adaptation in microtubule dynamics across frog species

Cellular processes are remarkably effective across diverse temperature ranges, even with highly conserved proteins. In the context of the microtubule cytoskeleton, which is critically involved in a wide range of cellular activities, this is particularly striking, as tubulin is one of the most conserved proteins while microtubule dynamic instability is highly temperature sensitive. Here, we leverage the diversity of natural tubulin variants from three closely related frog species that live at different temperatures. We determine the microtubule structure across all three species at between 3.0 and 3.6 Å resolution by cryo-electron microscopy and find small differences at the β-tubulin lateral interactions. Using in vitro reconstitution assays and quantitative biochemistry, we show that tubulin's free energy scales inversely with temperature. The observed weakening of lateral contacts and the low apparent activation energy for tubulin incorporation provide an explanation for the overall stability and higher growth rates of microtubules in cold-adapted frog species. This study thus broadens our conceptual framework for understanding microtubule dynamics and provides insights into how conserved cellular processes are tailored to different ecological niches.

Self-organization from organs to embryoids by activin in early amphibian development

Embryonic development is a complex self-organizing process orchestrated by a series of regulatory events at the molecular and cellular levels, resulting in the formation of a fully functional organism. This review focuses on activin protein as a mesoderm-inducing factor and the self-organizing properties it confers. Activin has been detected in both unfertilized eggs and embryos, suggesting its involvement in early developmental processes. To explore its effects, animal cap cells-pluripotent cells from the animal pole of amphibian blastula-stage embryos-were treated with varying concentrations of activin. The results showed that activin induced mesodermal tissues, including blood, muscle, and notochord, in a dose-dependent manner. Co-treatment with activin and retinoic acid further promoted the development of kidney and pancreatic tissues, while activin alone stimulated the formation of beating cardiac tissue. In subsequent experiments, high concentrations of activin conferred an organizer-like activity on animal cap cells. The pretreatment duration affected outcomes: longer exposure induced anterior structures, such as eyes, while shorter exposure resulted in posterior structures, like tails. These findings reflect moderate self-assembly, where cells become increasingly organized. In another experiment, activin was used to create an artificial gradient. Explants cultured on this gradient developed into embryoids with well-defined anteroposterior, dorsoventral, and left-right axes, exemplifying higher-order self-organization. These results demonstrate that controlled activin gradients can drive the formation of nearly complete tadpole-like larvae, effectively recapitulating the processes of early embryogenesis. This system offers valuable insights into the mechanisms underlying axis formation and organogenesis, providing a promising platform for future research in developmental biology.

Population Genomics Reveals Local Adaptation Related to Temperature Variation in Two Stream Frog Species: Implications for Vulnerability to Climate Warming

Identifying populations at highest risk from climate change is a critical component of conservation efforts. However, vulnerability assessments are usually applied at the species level, even though intraspecific variation in exposure, sensitivity and adaptive capacity play a crucial role in determining vulnerability. Genomic data can inform intraspecific vulnerability by identifying signatures of local adaptation that reflect population-level variation in sensitivity and adaptive capacity. Here, we address the question of local adaptation to temperature and the genetic basis of thermal tolerance in two stream frogs (Ascaphus truei and A. montanus). Building on previous physiological and temperature data, we used whole-genome resequencing of tadpoles from four sites spanning temperature gradients in each species to test for signatures of local adaptation. To support these analyses, we developed the first annotated reference genome for A. truei. We then expanded the geographic scope of our analysis using targeted capture at an additional 11 sites per species. We found evidence of local adaptation to temperature based on physiological and genomic data in A. montanus and genomic data in A. truei, suggesting similar levels of sensitivity (i.e., susceptibility) among populations regardless of stream temperature. However, invariant thermal tolerances across temperatures in A. truei suggest that populations occupying warmer streams may be most sensitive. We identified high levels of evolutionary potential in both species based on genomic and physiological data. While further integration of these data is needed to comprehensively evaluate spatial variation in vulnerability, this work illustrates the value of genomics in identifying spatial patterns of climate change vulnerability.

Vesicular glutamate release is necessary for neural tube formation

The brain and spinal cord originate from a neural tube that is preceded by a flat structure known as the neural plate during early embryogenesis. In humans, failure of the neural plate to convert into a tube by the fourth week of pregnancy leads to neural tube defects (NTDs), birth defects with serious neurological consequences. The signaling mechanisms governing the process of neural tube morphogenesis are unclear. Here we show that in Xenopus laevis embryos, neural plate cells release glutamate during neural plate folding in a Ca 2+ and vesicular glutamate transporter-1 (VGluT1)-dependent manner. Vesicular release of glutamate elicits Ca 2+ transients in neural plate cells that correlate with activation of Erk1/2. Knocking down or out VGluT1 leads to NTDs through increased expression of Sox2, neural stem cell transcription factor, and neural plate cell proliferation. Exposure during early pregnancy to neuroactive drugs that disrupt these signaling mechanisms might increase the risk of NTDs in offspring.

Interplay of Light, Melatonin, and Circadian Genes in Skin Pigmentation Regulation

Circadian regulation of skin pigmentation is essential for thermoregulation, ultraviolet (UV) protection, and synchronization of skin cell renewal. This regulation involves both cell-autonomous photic responses and non-cell-autonomous hormonal control, particularly through melatonin produced in a light-sensitive manner. Photosensitive opsins, cryptochromes, and melatonin regulate circadian rhythms in skin pigment cells. We studied light/dark cycles and melatonin coordination in melanin synthesis and cell proliferation of Xenopus laevis melanophores. In vivo, tadpole pigmentation shows robust circadian regulation mainly hormone-driven, in that isolated melanophores respond strongly to melatonin but only slightly to light. Melanophore proliferation is faster in the dark and slower with melatonin as compared to a 12/12 light/dark cycle. Expression of circadian core genes (clock, bmal1, per1, per2, per3, cry1, cry2, and cry4) in melatonin-treated cells during the light phase mimics dark phase expression. Overexpression of individual Crys did not affect melanization or cell proliferation, likely due to their cooperative actions. Melanin synthesis was inhibited by circadian cycle deregulation through (a) pharmacological inhibition of Cry1 and Cry2 degradation with KL001, (b) continuous light or dark conditions, and (c) melatonin treatment. Our findings suggest that circadian cycle regulation, rather than proliferative capacity, alters melanization of melanophores.

Xenopus as a model system for studying pigmentation and pigmentary disorders

Human pigmentary disorders encompass a broad spectrum of phenotypic changes arising from disruptions in various stages of melanocyte formation, the melanogenesis process, or the transfer of pigment from melanocytes to keratinocytes. A large number of pigmentation genes associated with pigmentary disorders have been identified, many of them awaiting in vivo confirmation. A more comprehensive understanding of the molecular basis of pigmentary disorders requires a vertebrate animal model where changes in pigmentation are easily observable in vivo and can be combined to genomic modifications and gain/loss-of-function tools. Here we present the amphibian Xenopus with its unique features that fulfill these requirements. Changes in pigmentation are particularly easy to score in Xenopus embryos, allowing whole-organism based phenotypic screening. The development and behavior of Xenopus melanocytes closely mimic those observed in mammals. Interestingly, both Xenopus and mammalian skins exhibit comparable reactions to ultraviolet radiation. This review highlights how Xenopus constitutes an alternative and complementary model to the more commonly used mouse and zebrafish, contributing to the advancement of knowledge in melanocyte cell biology and related diseases.

Rapid expansion and specialization of the TAS2R bitter taste receptor family in amphibians

TAS2Rs are a family of G protein-coupled receptors that function as bitter taste receptors in vertebrates. Mammalian TAS2Rs have historically garnered the most attention, leading to our understanding of their roles in taste perception relevant to human physiology and behaviors. However, the evolution and functional implications of TAS2Rs in other vertebrate lineages remain less explored. Here, we identify 9,291 TAS2Rs from 661 vertebrate genomes. Large-scale phylogenomic analyses reveal that frogs and salamanders contain unusually high TAS2R gene content, in stark contrast to other vertebrate lineages. In most species, TAS2R genes are found in clusters; compared to other vertebrates, amphibians have additional clusters and more genes per cluster. We find that vertebrate TAS2Rs have few one-to-one orthologs between closely related species, although total TAS2R count is stable in most lineages. Interestingly, TAS2R count is proportional to the receptors expressed solely in extra-oral tissues. In vitro receptor activity assays uncover that many amphibian TAS2Rs function as tissue-specific chemosensors to detect ecologically important xenobiotics.

Histological and gene-expression analyses of pyloric sphincter formation during stomach metamorphosis in Xenopus laevis

During anuran metamorphosis from herbivorous tadpoles to carnivorous frogs, the gastrointestinal (GI) tract undergoes drastic remodeling, such as the formation of the stomach-intestine boundary and the development of the pyloric sphincter at the posterior end of the stomach. However, the morphogenetic process and molecular mechanisms of how the pyloric sphincter is formed during metamorphosis, instead of during embryogenesis as in amniotes, are largely uninvestigated. Using the African clawed frog Xenopus laevis, we histologically examined the development of the pylorus region from embryonic to froglet stages and performed spatiotemporal gene expression analyses. We found that the pyloric sphincter is formed at a flexure within the pyloric region during metamorphic climax, and that the pyloric and duodenal epithelia, which are morphologically indistinguishable before sphincter formation, become clearly demarcated by the sphincter at the end of metamorphosis. Consistent with these morphological changes, expression domains of a stomach marker barx1 and an intestine marker cdx2 overlapped until late metamorphic climax, but became separated after metamorphosis. Despite the absence of the sphincter before metamorphosis, various genes crucial for sphincter formation in amniotes were already expressed in the pylorus region of Xenopus embryos. RNA-sequencing analysis at pre-metamorphic and metamorphic-climax stages suggest unappreciated roles of genes, such as those for retinoic acid signaling and various transcription factors, in suppressing or promoting sphincter formation. These data provide histological and molecular insights into the heterochrony of the pyloric sphincter formation in amniotes and anurans.

The Synthesis of GABA during the Tailbud Stage Is Required for Axial Elongation in Xenopus laevis embryos

In Xenopus laevis , axial elongation beyond the tailbud stage requires gamma-aminobutyric acid (GABA). However, the role of GABA synthesized during early development in this process remains unclear. In this study, by treating embryos with allylglycine (AG), an inhibitor of GABA synthesis, we observed a significant reduction in axial elongation. This inhibition was rescued by exogenous GABA, demonstrating that GABA synthesis via glutamate decarboxylase (GAD) is essential for axial elongation after the tailbud stage. Our findings suggest that GABA-dependent elongation functions independently of mechanisms like convergent extension, which are crucial during early development.

Phylogeny and Polyploidy Evolution of the Suckers (Teleostei: Catostomidae)

Fishes in the cypriniform family Catostomidae (suckers) are evolutionary tetraploids. The use of nuclear markers in the phylogenetic study of this important group has been greatly hindered by the challenge of identifying paralogous copies of genes. In the present study, we used two different methods to separate the gene copies of five single-copy nuclear genes (i.e., RAG1, EGR2B, EGR3, IRBP2, and RAG2). For each gene, all sequences of Copy I formed a clade that was sister to the clade formed by all sequences of Copy II in the phylogenetic trees. The maternal and paternal progenitor of the tetraploid ancestor of the Catostomidae could not be determined. We also constructed a mitochondrial tree to reflect the maternal relationships among major catostomid lineages. Our data appear to support a sister relationship between Catostominae and a monophyletic group composed of Myxocyprininae, Cycleptinae, and Ictiobinae. However, within Catostominae, there is significant conflict between mitochondrial and nuclear data regarding the relationships among Erimyzonini, Catostomini, and Moxostomatini/Thoburnini. Many indels, unexpected stop codons, and possible gene loss were identified in one gene copy of RAG1, RAG2, and IRBP2. We believe that additional nuclear genome data are needed to better resolve the phylogenetic relationships within the family Catostomidae.

Single-cell atlas comparison across vertebrates reveals auditory cell evolution and mechanisms for hair cell regeneration

Mammals suffer permanent hearing impairment from the loss of auditory hair cells due to their inability to regenerate. In contrast, lower vertebrates exhibit extraordinary capacity for hair cell regeneration and hearing restoration, but the mechanisms remain unclear. Here we characterize the single-cell atlas of Xenopus laevis inner ear and perform a comprehensive comparison with mouse model. An exceptionally conserved inner ear neuronal cell type is discovered. The results reveal that the outer hair cells (OHCs) exist exclusively in mammals. Importantly, our analyses reveal an orchestrated gene expression program in Xenopus, characterized by upregulation of hair cell regeneration-related genes, coupled with downregulation of proliferation inhibitory genes. These findings unveil a natural feature of regenerative capacity in Xenopus, and provide molecular and evolutionary evidences for differential regenerative capacities across vertebrates. This work offers insights from amphibians into developing strategies to solve the challenges of hair cell regeneration in humans.

Marcks and Marcks-like 1 proteins promote spinal cord development and regeneration in Xenopus

Marcks and Marcksl1 are abundant proteins that shuttle between the cytoplasm and membrane to modulate multiple cellular processes, including cytoskeletal dynamics, proliferation, and secretion. Here, we performed loss- and gain-of-function experiments in Xenopus laevis to reveal the novel roles of these proteins in spinal cord development and regeneration. We show that Marcks and Marcksl1 have partly redundant functions and are required for normal neurite formation and proliferation of neuro-glial progenitors during embryonic spinal cord development and for its regeneration during tadpole stages. Rescue experiments in Marcks and Marcksl1 loss-of-function animals further suggested that some of the functions of Marcks and Marcksl1 in the spinal cord are mediated by phospholipid signaling. Taken together, these findings identify Marcks and Marcksl1 as critical new players in spinal cord development and regeneration and suggest new pathways to be targeted for therapeutic stimulation of spinal cord regeneration in human patients.

Rapid Sex Chromosome Turnover in African Clawed Frogs (Xenopus) and the Origins of New Sex Chromosomes

Sex chromosomes of some closely related species are not homologous, and sex chromosome turnover is often attributed to mechanisms that involve linkage to or recombination arrest around sex-determining loci. We examined sex chromosome turnover and recombination landscapes in African clawed frogs (genus Xenopus) with reduced representation genome sequences from 929 individuals from 19 species. We recovered extensive variation in sex chromosomes, including at least eight nonhomologous sex-associated regions-five newly reported here, with most maintaining female heterogamety, but two independent origins of Y chromosomes. Seven of these regions are found in allopolyploid species in the subgenus Xenopus, and all of these reside in one of their two subgenomes, which highlights functional asymmetry between subgenomes. In three species with chromosome-scale genome assemblies (Xenopus borealis, Xenopus laevis, and Xenopus tropicalis), sex-specific recombination landscapes have similar patterns of sex differences in rates and locations of recombination. Across these Xenopus species, sex-associated regions are significantly nearer chromosome ends than expected by chance, even though this is where the ancestral recombination rate is highest in both sexes before the regions became sex associated. As well, expansions of sex-associated recombination arrest occurred multiple times. New information on sex linkage along with among-species variation in female specificity of the sex-determining gene dm-w argues against a "jumping gene" model, where dm-w moves around the genome. The diversity of sex chromosomes in Xenopus raises questions about the roles of natural and sexual selection, polyploidy, the recombination landscape, and neutral processes in driving sex chromosome turnover in animal groups with mostly heterogametic females.

Rethinking retinoic acid self-regulation: A signaling robustness network approach

All-trans retinoic acid (ATRA) signaling is a major pathway regulating numerous differentiation, proliferation, and patterning processes throughout life. ATRA biosynthesis depends on the nutritional availability of vitamin A and other retinoids and carotenoids, while it is sensitive to dietary and environmental toxicants. This nutritional and environmental influence requires a robustness response that constantly fine-tunes the ATRA metabolism to maintain a context-specific, physiological range of signaling levels. The ATRA metabolic and signaling network is characterized by the existence of multiple enzymes, transcription factors, and binding proteins capable of performing the same activity. The partial spatiotemporal expression overlap of these enzymes and proteins yields different network compositions in the cells and tissues where this pathway is active. Genetic polymorphisms affecting the activity of individual network components further impact the network composition variability and the self-regulatory feedback response to ATRA fluctuations. Experiments directly challenging the robustness response uncovered a Pareto optimality in the ATRA network, such that some genetic backgrounds efficiently deal with excess ATRA but are very limited in their robustness response to reduced ATRA and vice versa. We discuss a network-focused framework to describe the robustness response and the Pareto optimality of the ATRA metabolic and signaling network.

Characterisation and automated quantification of induced seizure-related behaviours in Xenopus laevis tadpoles

Epilepsy, a clinical diagnosis characterised by paroxysmal episodes known as seizures, affects 1% of people worldwide. Safe and patient-specific treatment is vital and can be achieved by the development of rapid pre-clinical models of for identified epilepsy genes. Epilepsy can result from either brain injury or gene mutations, and can also be induced chemically. Xenopus laevis tadpoles could be a useful model for confirmation of variants of unknown significance found in epilepsy patients, and for drug re-purposing screens that could eventually lead to benefits for patients. Here, we characterise and quantify seizure-related behaviours in X. laevis tadpoles arrayed in 24-well plates. To provoke acute seizure behaviours, tadpoles were chemically induced with either pentylenetetrazole (PTZ) or 4-aminopyridine (4-AP). To test the capacity to adapt this method for drug testing, we also exposed induced tadpoles to the anti-seizure drug valproate (VPA). Four induced seizure-like behaviours were described and manually quantified, and two of these (darting, circling) could be accurately detected automatically, using the video analysis software TopScan. Additionally, we recorded swimming trajectories and mean swimming velocity. Automatic detection showed that either PTZ or 4-AP induced darting behaviour and increased mean swimming velocity compared to untreated controls. Both parameters were significantly reduced in the presence of VPA. In particular, darting behaviour was a shown to be a sensitive measure of epileptic seizure activity. While we could not automatically detect the full range of seizure behaviours, this method shows promise for future studies since X. laevis is a well-characterised and genetically tractable model organism.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Unveiling axolotl transcriptome for tissue regeneration with high-resolution annotation via long-read sequencing

Axolotls are known for their remarkable regeneration ability. Exploring their transcriptome provides insight into regenerative mechanisms. However, the current annotation of the axolotl transcriptome is limited, leaving the role of unannotated transcripts in regeneration unknown. To discourse this challenge, we exploited long-read sequencing technology, which enables direct observation of full-length RNA transcripts, greatly enhancing the coverage and accuracy of axolotl transcriptome annotation. By utilizing this method, we identified 222 novel gene loci and 4775 novel transcripts, which were quantified using short-read sequencing data. Through the inclusive analysis, we discovered novel homologs, potential functional proteins, noncoding RNAs, and alternative splicing events in key regeneration pathways. In particular, we identified novel transcripts with high protein-coding potential implicated in cell cycle regulation and musculoskeletal development, and regeneration were identified. Interestingly, alternative splice variants were also detected across diverse pathways critical to regeneration. This specifies that these novel transcripts potentially play vital roles underpinning the robust regenerative capacities of axolotls. Single-cell transcriptomic analysis further revealed these isoforms to predominantly exist in axolotl limb chondrocytes and mature tissue cell populations. Overall, the findings significantly advanced consideration of the axolotl transcriptome and provided a new perspective for understanding the mechanisms of regenerative abilities of axolotls.

Repeat-Rich Regions Cause False-Positive Detection of NUMTs: A Case Study in Amphibians Using an Improved Cane Toad Reference Genome

Mitochondrial DNA (mtDNA) has been widely used in genetics research for decades. Contamination from nuclear DNA of mitochondrial origin (NUMTs) can confound studies of phylogenetic relationships and mtDNA heteroplasmy. Homology searches with mtDNA are widely used to detect NUMTs in the nuclear genome. Nevertheless, false-positive detection of NUMTs is common when handling repeat-rich sequences, while fragmented genomes might result in missing true NUMTs. In this study, we investigated different NUMT detection methods and how the quality of the genome assembly affects them. We presented an improved nuclear genome assembly (aRhiMar1.3) of the invasive cane toad (Rhinella marina) with additional long-read Nanopore and 10× linked-read sequencing. The final assembly was 3.47 Gb in length with 91.3% of tetrapod universal single-copy orthologs (n = 5,310), indicating the gene-containing regions were well assembled. We used 3 complementary methods (NUMTFinder, dinumt, and PALMER) to study the NUMT landscape of the cane toad genome. All 3 methods yielded consistent results, showing very few NUMTs in the cane toad genome. Furthermore, we expanded NUMT detection analyses to other amphibians and confirmed a weak relationship between genome size and the number of NUMTs present in the nuclear genome. Amphibians are repeat-rich, and we show that the number of NUMTs found in highly repetitive genomes is prone to inflation when using homology-based detection without filters. Together, this study provides an exemplar of how to robustly identify NUMTs in complex genomes when confounding effects on mtDNA analyses are a concern.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, anti-predator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium ) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

Carboxy-terminal polyglutamylation regulates signaling and phase separation of the Dishevelled protein

Polyglutamylation is a reversible posttranslational modification that is catalyzed by enzymes of the tubulin tyrosine ligase-like (TTLL) family. Here, we found that TTLL11 generates a previously unknown type of polyglutamylation that is initiated by the addition of a glutamate residue to the free C-terminal carboxyl group of a substrate protein. TTLL11 efficiently polyglutamylates the Wnt signaling protein Dishevelled 3 (DVL3), thereby changing the interactome of DVL3. Polyglutamylation increases the capacity of DVL3 to get phosphorylated, to undergo phase separation, and to act in the noncanonical Wnt pathway. Both carboxy-terminal polyglutamylation and the resulting reduction in phase separation capacity of DVL3 can be reverted by the deglutamylating enzyme CCP6, demonstrating a causal relationship between TTLL11-mediated polyglutamylation and phase separation. Thus, C-terminal polyglutamylation represents a new type of posttranslational modification, broadening the range of proteins that can be modified by polyglutamylation and providing the first evidence that polyglutamylation can modulate protein phase separation.

Functional and regulatory diversification of Period genes responsible for circadian rhythm in vertebrates

The Period genes (Per) play essential roles in modulating the molecular circadian clock timing in a broad range of species, which regulates the physiological and cellular rhythms through the transcription-translation feedback loop. While the Period gene paralogs are widely observed among vertebrates, the evolutionary history and the functional diversification of Per genes across vertebrates are not well known. In this study, we comprehensively investigated the evolution of Per genes at the copy number and sequence levels, including de novo binding motif discovery by comparative genomics. We also determined the lineage-specific transcriptome landscape across tissues and developmental stages and phenotypic effects in public RNA-seq data sets of model species. We observed multiple lineage-specific gain and loss events Per genes, though no simple association was observed between ecological factors and Per gene numbers in each species. Among salmonid fish species, the per3 gene has been lost in the majority, whereas those retaining the per3 gene exhibit not a signature of relaxed selective constraint but rather a signature of intensified selection. We also determined the signature of adaptive diversification of the CRY-binding region in Per1 and Per3, which modulates the circadian rhythm. We also discovered putative regulatory sequences, which are lineage-specific, suggesting that these cis-regulatory elements may have evolved rapidly and divergently across different lineages. Collectively, our findings revealed the evolution of Per genes and their fine-tuned contribution to the plastic and precise regulation of circadian rhythms in various vertebrate taxa.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, antipredator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

Plant chromosome polytenization contributes to suppression of root growth in high polyploids

Autopolyploidization, which refers to a polyploidization via genome duplication without hybridization, promotes growth in autotetraploids, but suppresses growth in high polyploids (autohexaploids or auto-octoploids). The mechanism underlying this growth suppression (i.e. 'high-ploidy syndrome') has not been comprehensively characterized. In this study, we conducted a kinematic analysis of the root apical meristem cells in Arabidopsis thaliana autopolyploids (diploid, tetraploid, hexaploid, and octoploid) to determine the effects of the progression of genome duplication on root growth. The results of the root growth analysis showed that tetraploidization increases the cell volume, but decreases cell proliferation. However, cell proliferation and volume growth are suppressed in high polyploids. Whole-mount fluorescence in situ hybridization analysis revealed extensive chromosome polytenization in the region where cell proliferation does not usually occur in the roots of high polyploids, which is likely to be at least partly correlated with the suppression of endoreduplication. The study findings indicate that chromosome polytenization is important for the suppressed growth of high polyploids.

Discovery of numerous novel Helitron-like elements in eukaryote genomes using HELIANO

Helitron-like elements (HLEs) are widespread eukaryotic DNA transposons employing a rolling-circle transposition mechanism. Despite their prevalence in fungi, animals, and plant genomes, identifying Helitrons remains a formidable challenge. We introduce HELIANO, a software for annotating and classifying autonomous and non-autonomous HLE sequences from whole genomes. HELIANO overcomes several limitations of existing tools in speed and accuracy, demonstrated through benchmarking and its application to the complex genomes of frogs (Xenopus tropicalis and Xenopus laevis) and rice (Oryza sativa), where it uncovered numerous previously unidentified HLEs. In an extensive analysis of 404 eukaryote genomes, we found HLEs widely distributed across phyla, with exceptions in specific taxa. HELIANO's application led to the discovery of numerous new HLEs in land plants and identified 20 protein domains captured by certain autonomous HLE families. A comprehensive phylogenetic analysis further classified HLEs into two primary clades, HLE1 and HLE2, and revealed nine subgroups, some of which are enriched within specific taxa. The future use of HELIANO promises to improve the global analysis of HLEs across genomes, significantly advancing our understanding of this fascinating transposon superfamily.

The evolution of tenascins

BACKGROUND

The evolution of extracellular matrix is tightly linked to the evolution of organogenesis in metazoans. Tenascins are extracellular matrix glycoproteins of chordates that participate in integrin-signaling and morphogenetic events. Single tenascins are encoded by invertebrate chordates, and multiple tenascin paralogs are found in vertebrates (designated tenascin-C, tenascin-R, tenascin-W and tenascin-X) yet, overall, the evolution of this family has remained unclear.

RESULTS

This study examines the genomes of hemichordates, cephalochordates, tunicates, agnathans, cartilaginous fishes, lobe-finned fishes, ray-finned fishes and representative tetrapods to identify predicted tenascin proteins. We comprehensively assess their evolutionary relationships by sequence conservation, molecular phylogeny and examination of conservation of synteny of the encoding genes. The resulting new evolutionary model posits the origin of tenascin in an ancestral chordate, with tenascin-C-like and tenascin-R-like paralogs emerging after a whole genome duplication event in an ancestral vertebrate. Tenascin-X appeared following a second round of whole genome duplication in an ancestral gnathostome, most likely from duplication of the gene encoding the tenascin-R homolog. The fourth gene, encoding tenascin-W (also known as tenascin-N), apparently arose from a local duplication of tenascin-R.

CONCLUSIONS

The diversity of tenascin paralogs observed in agnathans and gnathostomes has evolved through selective retention of novel genes that arose from a combination of whole genome and local duplication events. The evolutionary appearance of specific tenascin paralogs coincides with the appearance of vertebrate-specific cell and tissue types where the paralogs are abundantly expressed, such as the endocranium and facial skeleton (tenascin-C), an expanded central nervous system (tenascin-R), and bone (tenascin-W).

Redesigning the Drosophila histone gene cluster: an improved genetic platform for spatiotemporal manipulation of histone function

Mutating replication-dependent (RD) histone genes is an important tool for understanding chromatin-based epigenetic regulation. Deploying this tool in metazoans is particularly challenging because RD histones in these organisms are typically encoded by many genes, often located at multiple loci. Such gene arrangements make the ability to generate homogenous histone mutant genotypes by site-specific gene editing quite difficult. Drosophila melanogaster provides a solution to this problem because the RD histone genes are organized into a single large tandem array that can be deleted and replaced with transgenes containing mutant histone genes. In the last ∼15 years several different RD histone gene replacement platforms were developed using this simple strategy. However, each platform contains weaknesses that preclude full use of the powerful developmental genetic capabilities available to Drosophila researchers. Here we describe the development of a newly engineered platform that rectifies many of these weaknesses. We used CRISPR to precisely delete the RD histone gene array (HisC), replacing it with a multifunctional cassette that permits site-specific insertion of either one or two synthetic gene arrays using selectable markers. We designed this cassette with the ability to selectively delete each of the integrated gene arrays in specific tissues using site-specific recombinases. We also present a method for rapidly synthesizing histone gene arrays of any genotype using Golden Gate cloning technologies. These improvements facilitate the generation of histone mutant cells in various tissues at different stages of Drosophila development and provide an opportunity to apply forward genetic strategies to interrogate chromatin structure and gene regulation.

Time of day dependent changes in embryonic heart rate are detectable after maturation of rhythmic circadian gene expression in the eye, but before the heart in Xenopus laevis tadpoles cultured in LD

We systematically characterized onset of expression of circadian genes in the embryonic eye and heart of Xenopus laevis tadpoles. We found that period1 (per1) and nr1d1 ( rev-erbα) were the first circadian genes to display significant 24-hour rhythms in the developing eye and heart in a 12-hour light-dark cycle (LD). Rhythmic expression of both oscillator and output genes were present in the eye by 2.75 days post fertilization (dpf), but not in 15 dpf hearts. Surprisingly, rhythmic oscillation of heart rate occurred after 3.2 dpf suggesting that heart rate may be controlled directly by light or indirectly by the pineal in LD.

Correlation Between Subgenome-biased DNA Loss and DNA Transposon Activation Following Hybridization in the Allotetraploid Xenopus Frogs

In certain tetraploid species resulting from interspecific hybridization, one parent's subgenome is known to selectively undergo DNA loss. The molecular mechanisms behind this remain unclear. In our study, we compared the genomes of a standard diploid species with two allotetraploid species from the Xenopus genus, both possessing L (longer) and S (shorter) homoeologous subgenomes. We observed substantial gene losses and intergenic DNA deletions in both the S and L subgenomes of the tetraploid species. Gene losses were around 1,000 to 3,000 for L and 4,000 to 6,000 for S, with especially prominent losses in the S subgenome. Many of these losses likely occurred shortly after interspecific hybridization in both L/S subgenomes. We also deduced frequent large inversions in the S subgenome. Upon reassessing transposon dynamics using updated genome databases, we reaffirmed heightened DNA transposon activity during the hybridization, as previously reported. We next investigated whether S subgenome-biased DNA loss could be correlated with the activation of DNA transposons following hybridization. Notably, distinct patterns were observed in the dynamics of DNA transposons between the L and S subgenomes. Several DNA transposon subfamilies correlated positively with DNA deletions in the S subgenome and negatively in the L subgenome. Based on these results, we propose a model that, upon and after hybridization between two related diploid Xenopus species, the mixture of their genomes resulted in the derepression of DNA transposons, especially in the S subgenome, leading to selective DNA loss in the S subgenome.

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

Ionizing radiation induces neurotoxicity in Xenopus laevis embryos through neuroactive ligand-receptor interaction pathway

Ionizing radiation (IR) poses a significant threat to both the natural environment and biological health. Exposure to specific doses of ionizing radiation early in an organism's development can lead to developmental toxicity, particularly neurotoxicity. Through experimentation with Xenopus laevis (X. laevis), we examined the effects of radiation on early developmental stage. Our findings revealed that radiation led to developmental abnormalities and mortality in X. laevis embryos in a dose-dependent manner, disrupting redox homeostasis and inducing cell apoptosis. Additionally, radiation caused neurotoxic effects, resulting in abnormal behavior and neuron damage in the embryos. Further investigation into the underlying mechanisms of radiation-induced neurotoxicity indicated the potential involvement of the neuroactive ligand-receptor interaction pathway, which was supported by RNA-Seq analysis. Validation of gene expression associated with this pathway and analysis of neurotransmitter levels confirmed our hypothesis. In addition, we further validated the important role of this signaling pathway in radiation-induced neurotoxicity through edaravone rescue experiments. This research establishes a valuable model for radiation damage studying and provides some insight into radiation-induced neurotoxicity mechanisms.

Temporal Transcriptomic Profiling of the Developing Xenopus laevis Eye

Retinal progenitor cells (RPCs) are a multipotent and highly proliferative population that give rise to all retinal cell types during organogenesis. Defining their molecular signature is a key step towards identifying suitable approaches to treat visual impairments. Here, we performed RNA sequencing of whole eyes from Xenopus at three embryonic stages and used differential expression analysis to define the transcriptomic profiles of optic tissues containing proliferating and differentiating RPCs during retinogenesis. Gene Ontology and KEGG pathway analyses showed that genes associated with developmental pathways (including Wnt and Hedgehog signaling) were upregulated during the period of active RPC proliferation in early retinal development (Nieuwkoop Faber st. 24 and 27). Developing eyes had dynamic expression profiles and shifted to enrichment for metabolic processes and phototransduction during RPC progeny specification and differentiation (st. 35). Furthermore, conserved adult eye regeneration genes were also expressed during early retinal development, including sox2, pax6, nrl, and Notch signaling components. The eye transcriptomic profiles presented here span RPC proliferation to retinogenesis and include regrowth-competent stages. Thus, our dataset provides a rich resource to uncover molecular regulators of RPC activity and will allow future studies to address regulators of RPC proliferation during eye repair and regrowth.