Transcriptomic insights into genetic diversity of protein-coding genes in X. laevis

Enhanced Loss of Retinoic Acid Network Genes in Xenopus laevis Achieves a Tighter Signal Regulation

Retinoic acid (RA) is a major regulatory signal during embryogenesis produced from vitamin A (retinol) by an extensive, autoregulating metabolic and signaling network to prevent fluctuations that result in developmental malformations. Xenopus laevis is an allotetraploid hybrid frog species whose genome includes L (long) and S (short) chromosomes from the originating species. Evolutionarily, the X. laevis subgenomes have been losing either L or S homoeologs in about 43% of genes to generate singletons. In the RA network, out of the 47 genes, about 47% have lost one of the homoeologs, like the genome average. Interestingly, RA metabolism genes from storage (retinyl esters) to retinaldehyde production exhibit enhanced gene loss with 75% singletons out of 28 genes. The effect of this gene loss on RA signaling autoregulation was studied. Employing transient RA manipulations, homoeolog gene pairs were identified in which one homoeolog exhibits enhanced responses or looser regulation than the other, while in other pairs both homoeologs exhibit similar RA responses. CRISPR/Cas9 targeting of individual homoeologs to reduce their activity supports the hypothesis where the RA metabolic network gene loss results in tighter network regulation and more efficient RA robustness responses to overcome complex regulation conditions.

Retinoic Acid Fluctuation Activates an Uneven, Direction-Dependent Network-Wide Robustness Response in Early Embryogenesis

Robustness is a feature of regulatory pathways to ensure signal consistency in light of environmental changes or genetic polymorphisms. The retinoic acid (RA) pathway, is a central developmental and tissue homeostasis regulatory signal, strongly dependent on nutritional sources of retinoids and affected by environmental chemicals. This pathway is characterized by multiple proteins or enzymes capable of performing each step and their integration into a self-regulating network. We studied RA network robustness by transient physiological RA signaling disturbances followed by kinetic transcriptomic analysis of the recovery during embryogenesis. The RA metabolic network was identified as the main regulated module to achieve signaling robustness using an unbiased pattern analysis. We describe the network-wide responses to RA signal manipulation and found the feedback autoregulation to be sensitive to the direction of the RA perturbation: RA knockdown exhibited an upper response limit, whereas RA addition had a minimal feedback-activation threshold. Surprisingly, our robustness response analysis suggests that the RA metabolic network regulation exhibits a multi-objective optimization, known as Pareto optimization, characterized by trade-offs between competing functionalities. We observe that efficient robustness to increasing RA is accompanied by worsening robustness to reduced RA levels and vice versa. This direction-dependent trade-off in the network-wide feedback response, results in an uneven robustness capacity of the RA network during early embryogenesis, likely a significant contributor to the manifestation of developmental defects.

Natural size variation among embryos leads to the corresponding scaling in gene expression

Xenopus laevis frogs from laboratory stocks normally lay eggs exhibiting extensive size variability. We find that these initial size differences subsequently affect the size of the embryos prior to the onset of growth, and the size of tadpoles during the growth period. Even though these tadpoles differ in size, their tissues, organs, and structures always seem to be properly proportioned, i.e. they display static allometry. Initial axial patterning events in Xenopus occur in a spherical embryo, allowing easy documentation of their size-dependent features. We examined the size distribution of early Xenopus laevis embryos and measured diameters that differed by about 38% with a median of about 1.43 ​mm. This range of embryo sizes corresponds to about a 1.9-fold difference in surface area and a 2.6-fold difference in volume. We examined the relationship between embryo size and gene expression and observed a significant correlation between diameter and RNA content during gastrula stages. In addition, we investigated the expression levels of genes that pattern the mesoderm, induce the nervous system and mediate the progression of ectodermal cells to neural precursors in large and small embryos. We found that most of these factors were expressed at levels that scaled with the different embryo sizes and total embryo RNA content. In agreement with the changes in transcript levels, the expression domains in larger embryos increased proportionally with the increase in surface area, maintaining their relative expression domain size in relation to the total size of the embryo. Thus, our study identified a mechanism for adapting gene expression domains to embryo size by adjusting the transcript levels of the genes regulating mesoderm induction and patterning. In the neural plate, besides the scaling of the expression domains, we observed similar cell sizes and cell densities in small and large embryos suggesting that additional cell divisions took place in large embryos to compensate for the increased size. Our results show in detail the size variability among Xenopus laevis embryos and the transcriptional adaptation to scale gene expression with size. The observations further support the involvement of BMP/ADMP signaling in the scaling process.

Genomics Methods for Xenopus Embryos and Tissues

High-throughput sequencing methods have created exciting opportunities to explore the regulatory landscape of the entire genome. Here we introduce methods to characterize the genomic locations of bound proteins, open chromatin, and sites of DNA-DNA contact in Xenopus embryos. These methods include chromatin immunoprecipitation followed by sequencing (ChIP-seq), a combination of DNase I digestion and sequencing (DNase-seq), the assay for transposase-accessible chromatin and sequencing (ATAC-seq), and the use of proximity-based DNA ligation followed by sequencing (Hi-C).

The dynamics of gene expression in vertebrate embryogenesis at single-cell resolution

Time series of single-cell transcriptome measurements can reveal dynamic features of cell differentiation pathways. From measurements of whole frog embryos spanning zygotic genome activation through early organogenesis, we derived a detailed catalog of cell states in vertebrate development and a map of differentiation across all lineages over time. The inferred map recapitulates most if not all developmental relationships and associates new regulators and marker genes with each cell state. We find that many embryonic cell states appear earlier than previously appreciated. We also assess conflicting models of neural crest development. Incorporating a matched time series of zebrafish development from a companion paper, we reveal conserved and divergent features of vertebrate early developmental gene expression programs.