Retinoic acid (RA) is a central developmental signal whose perturbation results in teratogenic outcomes. We performed transient physiological RA signaling disturbances during embryogenesis followed by kinetic RNAseq and high-throughput qPCR analysis of the recovery. Unbiased comparative pattern analysis identified the RA network as a key differentially regulated module aimed at achieving RA signaling robustness. We analyzed this module using a principal curve-based approach determining clutch-wise variability, and organized the results into a robustness efficiency matrix comparing the shifts in RA feedback regulation and hox gene expression. We found that feedback autoregulation was sensitive to the direction of the RA perturbation: RA knock-down exhibited an upper response limit, whereas RA addition did not activate a feedback response below a threshold. These results suggest an asymmetric capacity for robust feedback control of RA signal during early development based on genetic polymorphisms, likely a significant contributor to the manifestation of developmental defects.

No abstract

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.

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 (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 46% have lost one of the homoeologs, like the genome average. In contrast, 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 homeolog 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.