Global analysis of the transcriptional network controlling Xenopus endoderm formation

Differentiation success of reprogrammed cells is heterogeneous in vivo and modulated by somatic cell identity memory

Nuclear reprogramming can change cellular fates. Yet, reprogramming efficiency is low, and the resulting cell types are often not functional. Here, we used nuclear transfer to eggs to follow single cells during reprogramming in vivo. We show that the differentiation success of reprogrammed cells varies across cell types and depends on the expression of genes specific to the previous cellular identity. We find subsets of reprogramming-resistant cells that fail to form functional cell types, undergo cell death, or disrupt normal body patterning. Reducing expression levels of genes specific to the cell type of origin leads to better reprogramming and improved differentiation trajectories. Thus, our work demonstrates that failing to reprogram in vivo is cell type specific and emphasizes the necessity of minimizing aberrant transcripts of the previous somatic identity for improving reprogramming.

Computational approaches to understand transcription regulation in development

Gene regulatory networks (GRNs) serve as useful abstractions to understand transcriptional dynamics in developmental systems. Computational prediction of GRNs has been successfully applied to genome-wide gene expression measurements with the advent of microarrays and RNA-sequencing. However, these inferred networks are inaccurate and mostly based on correlative rather than causative interactions. In this review, we highlight three approaches that significantly impact GRN inference: (1) moving from one genome-wide functional modality, gene expression, to multi-omics, (2) single cell sequencing, to measure cell type-specific signals and predict context-specific GRNs, and (3) neural networks as flexible models. Together, these experimental and computational developments have the potential to significantly impact the quality of inferred GRNs. Ultimately, accurately modeling the regulatory interactions between transcription factors and their target genes will be essential to understand the role of transcription factors in driving developmental gene expression programs and to derive testable hypotheses for validation.

Nodal and churchill1 position the expression of a notch ligand during Xenopus germ layer segregation

Churchill and Nodal signaling, which participate in vertebrates’ germ layer induction, position a domain of Delta/Notch activity, which refines germ layer boundaries during frog gastrulation.

In vertebrates, Nodal signaling plays a major role in endomesoderm induction, but germ layer delimitation is poorly understood. In avian embryos, the neural/mesoderm boundary is controlled by the transcription factor CHURCHILL1, presumably through the repressor ZEB2, but there is scarce knowledge about its role in other vertebrates. During amphibian gastrulation, Delta/Notch signaling refines germ layer boundaries in the marginal zone, but it is unknown the place this pathway occupies in the network comprising Churchill1 and Nodal. Here, we show that Xenopus churchill1 is expressed in the presumptive neuroectoderm at mid-blastula transition and during gastrulation, upregulates zeb2, prevents dll1 expression in the neuroectoderm, and favors neuroectoderm over endomesoderm development. Nodal signaling prevents dll1 expression in the endoderm but induces it in the presumptive mesoderm, from where it activates Notch1 and its target gene hes4 in the non-involuting marginal zone. We propose a model where Nodal and Churchill1 position Dll1/Notch1/Hes4 domains in the marginal zone, ensuring the delimitation between mesoderm and neuroectoderm.

Nascent transcriptome reveals orchestration of zygotic genome activation in early embryogenesis

Early embryo development requires maternal-to-zygotic transition, during which transcriptionally silent nuclei begin widespread gene expression during zygotic genome activation (ZGA).1-3 ZGA is vital for early cell fating and germ-layer specification,3,4 and ZGA timing is regulated by multiple mechanisms.1-5 However, controversies remain about whether these mechanisms are interrelated and vary among species6-10 and whether the timing of germ-layer-specific gene activation is temporally ordered.11,12 In some embryonic models, widespread ZGA onset is spatiotemporally graded,13,14 yet it is unclear whether the transcriptome follows this pattern. A major challenge in addressing these questions is to accurately measure the timing of each gene activation. Here, we metabolically label and identify the nascent transcriptome using 5-ethynyl uridine (5-EU) in Xenopus blastula embryos. We find that EU-RNA-seq outperforms total RNA-seq in detecting the ZGA transcriptome, which is dominated by transcription from maternal-zygotic genes, enabling improved ZGA timing determination. We uncover discrete spatiotemporal patterns for individual gene activation, a majority following a spatial pattern of ZGA that is correlated with a cell size gradient.14 We further reveal that transcription necessitates a period of developmental progression and that ZGA can be precociously induced by cycloheximide, potentially through elongation of interphase. Finally, most ectodermal genes are activated earlier than endodermal genes, suggesting a temporal orchestration of germ-layer-specific genes, potentially linked to the spatially graded pattern of ZGA. Together, our study provides fundamental new insights into the composition and dynamics of the ZGA transcriptome, mechanisms regulating ZGA timing, and its role in the onset of early cell fating.

Massive Gene Loss and Function Shuffling in Appendicularians Stretch the Boundaries of Chordate Wnt Family Evolution

Gene loss is a pervasive source of genetic variation that influences species evolvability, biodiversity and the innovation of evolutionary adaptations. To better understand the evolutionary patterns and impact of gene loss, here we investigate as a case study the evolution of the wingless (Wnt) family in the appendicularian tunicate Oikopleura dioica, an emergent EvoDevo model characterized by its proneness to lose genes among chordates. Genome survey and phylogenetic analyses reveal that only four of the thirteen Wnt subfamilies have survived in O. dioica—Wnt5, Wnt10, Wnt11, and Wnt16,—representing the minimal Wnt repertoire described in chordates. While the loss of Wnt4 and Wnt8 likely occurred in the last common ancestor of tunicates, representing therefore a synapomorphy of this subphylum, the rest of losses occurred during the evolution of appendicularians. This work provides the first complete Wnt developmental expression atlas in a tunicate and the first insights into the evolution of Wnt developmental functions in appendicularians. Our work highlights three main evolutionary patterns of gene loss: (1) conservation of ancestral Wnt expression domains not affected by gene losses; (2) function shuffling among Wnt paralogs accompanied by gene losses; and (3) extinction of Wnt expression in certain embryonic directly correlated with gene losses. Overall our work reveals that in contrast to “conservative” pattern of evolution of cephalochordates and vertebrates, O. dioica shows an even more radical “liberal” evolutionary pattern than that described ascidian tunicates, stretching the boundaries of the malleability of Wnt family evolution in chordates.

Modeling endoderm development and disease in Xenopus

The endoderm is the innermost germ layer that forms the linings of the respiratory and gastrointestinal tracts, and their associated organs, during embryonic development. Xenopus embryology experiments have provided fundamental insights into how the endoderm develops in vertebrates, including the critical role of TGFβ-signaling in endoderm induction,elucidating the gene regulatory networks controlling germ layer development and the key molecular mechanisms regulating endoderm patterning and morphogenesis. With new genetic, genomic, and imaging approaches, Xenopus is now routinely used to model human disease, discover mechanisms underlying endoderm organogenesis, and inform differentiation protocols for pluripotent stem cell differentiation and regenerative medicine applications. In this chapter, we review historical and current discoveries of endoderm development in Xenopus, then provide examples of modeling human disease and congenital defects of endoderm-derived organs using Xenopus.

A temporally resolved transcriptome for developing "Keller" explants of the Xenopus laevis dorsal marginal zone

BACKGROUND

Explanted tissues from vertebrate embryos reliably develop in culture and have provided essential paradigms for understanding embryogenesis, from early embryological investigations of induction, to the extensive study of Xenopus animal caps, to the current studies of mammalian gastruloids. Cultured explants of the Xenopus dorsal marginal zone ("Keller" explants) serve as a central paradigm for studies of convergent extension cell movements, yet we know little about the global patterns of gene expression in these explants.

RESULTS

In an effort to more thoroughly develop this important model system, we provide here a time-resolved bulk transcriptome for developing Keller explants.

CONCLUSIONS

The dataset reported here provides a useful resource for those using Keller explants for studies of morphogenesis and provide genome-scale insights into the temporal patterns of gene expression in an important tissue when explanted and grown in culture.

Sox17 and β-catenin co-occupy Wnt-responsive enhancers to govern the endoderm gene regulatory network

Lineage specification is governed by gene regulatory networks (GRNs) that integrate the activity of signaling effectors and transcription factors (TFs) on enhancers. Sox17 is a key transcriptional regulator of definitive endoderm development, and yet, its genomic targets remain largely uncharacterized. Here, using genomic approaches and epistasis experiments, we define the Sox17-governed endoderm GRN in Xenopus gastrulae. We show that Sox17 functionally interacts with the canonical Wnt pathway to specify and pattern the endoderm while repressing alternative mesectoderm fates. Sox17 and β-catenin co-occupy hundreds of key enhancers. In some cases, Sox17 and β-catenin synergistically activate transcription apparently independent of Tcfs, whereas on other enhancers, Sox17 represses β-catenin/Tcf-mediated transcription to spatially restrict gene expression domains. Our findings establish Sox17 as a tissue-specific modifier of Wnt responses and point to a novel paradigm where genomic specificity of Wnt/β-catenin transcription is determined through functional interactions between lineage-specific Sox TFs and β-catenin/Tcf transcriptional complexes. Given the ubiquitous nature of Sox TFs and Wnt signaling, this mechanism has important implications across a diverse range of developmental and disease contexts.

Single Amino Acid Change Underlies Distinct Roles of H2A.Z Subtypes in Human Syndrome

The developmental disorder Floating-Harbor syndrome (FHS) is caused by heterozygous truncating mutations in SRCAP, a gene encoding a chromatin remodeler mediating incorporation of histone variant H2A.Z. Here, we demonstrate that FHS-associated mutations result in loss of SRCAP nuclear localization, alter neural crest gene programs in human in vitro models and Xenopus embryos, and cause craniofacial defects. These defects are mediated by one of two H2A.Z subtypes, H2A.Z.2, whose knockdown mimics and whose overexpression rescues the FHS phenotype. Selective rescue by H2A.Z.2 is conferred by one of the three amino acid differences between the H2A.Z subtypes, S38/T38. We further show that H2A.Z.1 and H2A.Z.2 genomic occupancy patterns are qualitatively similar, but quantitatively distinct, and H2A.Z.2 incorporation at AT-rich enhancers and expression of their associated genes are both sensitized to SRCAP truncations. Altogether, our results illuminate the mechanism underlying a human syndrome and uncover selective functions of H2A.Z subtypes during development.

Impaired Gene Expression Due to Iodine Excess in the Development and Differentiation of Endoderm and Thyroid Is Associated with Epigenetic Changes

Background: Thyroid hormone (TH) synthesis is essential for the control of development, growth, and metabolism in vertebrates and depends on a sufficient dietary iodine intake. Importantly, both iodine deficiency and iodine excess (IE) impair TH synthesis, causing serious health problems especially during fetal/neonatal development. While it is known that IE disrupts thyroid function by inhibiting thyroid gene expression, its effects on thyroid development are less clear. Accordingly, this study sought to investigate the effects of IE during the embryonic development/differentiation of endoderm and the thyroid gland. Methods: We used the murine embryonic stem (ES) cell model of in vitro directed differentiation to assess the impact of IE on the generation of endoderm and thyroid cells. Additionally, we subjected endoderm and thyroid explants obtained during early gestation to IE and evaluated gene and protein expression of endodermal markers in both models. Results: ES cells were successfully differentiated into endoderm cells and, subsequently, into thyrocytes expressing the specific thyroid markers Tshr, Slc5a5, Tpo, and Tg. IE exposure decreased the messenger RNA (mRNA) levels of the main endoderm markers Afp, Crcx4, Foxa1, Foxa2, and Sox17 in both ES cell-derived endoderm cells and embryonic explants. Interestingly, IE also decreased the expression of the main thyroid markers in ES cell-derived thyrocytes and thyroid explants. Finally, we demonstrate that DNA methyltransferase expression was increased by exposure to IE, and this was accompanied by hypermethylation and hypoacetylation of histone H3, pointing to an association between the gene repression triggered by IE and the observed epigenetic changes. Conclusions: These data establish that IE treatment is deleterious for embryonic endoderm and thyroid gene expression.

Retinoic Acid-Regulated Target Genes During Development: Integrative Genomics Analysis

Retinoic acid (RA), a major natural active metabolite of vitamin A (VA) is well known to play critical roles in embryonic development. The effects of RA are mediated by nuclear receptors (RARs), which regulate the expression of gene batteries involved in cell growth and differentiation. Since the early 1990s several laboratories have focused on understanding how RA-regulated genes and RAR binding sites operate by studying the differentiation of embryonal carcinoma cells and embryonic stem cells. The development of hybridization-based microarray technology and high performance software analysis programs has allowed the characterization of thousands of RA-regulated genes. During the two last decades, publication of the genome sequence of various organisms has allowed advances in massive parallel sequencing and bioinformatics analysis of genome-wide data sets. These new generation sequencing (NGS) technologies have revolutionized the field by providing a global integrated picture of RA-regulated gene networks and the regulatory programs involved in cell fate decisions during embryonal carcinoma and embryonic stem cells differentiation. Now the challenge is to reconstruct the RA-regulated gene networks at the single cell level during the development of specialized embryonic tissues.

Modeling congenital kidney diseases in Xenopus laevis

Congenital anomalies of the kidney and urinary tract (CAKUT) occur in ∼1/500 live births and are a leading cause of pediatric kidney failure. With an average wait time of 3-5 years for a kidney transplant, the need is high for the development of new strategies aimed at reducing the incidence of CAKUT and preserving renal function. Next-generation sequencing has uncovered a significant number of putative causal genes, but a simple and efficient model system to examine the function of CAKUT genes is needed. Xenopus laevis (frog) embryos are well-suited to model congenital kidney diseases and to explore the mechanisms that cause these developmental defects. Xenopus has many advantages for studying the kidney: the embryos develop externally and are easily manipulated with microinjections, they have a functional kidney in ∼2 days, and 79% of identified human disease genes have a verified ortholog in Xenopus. This facilitates high-throughput screening of candidate CAKUT-causing genes. In this Review, we present the similarities between Xenopus and mammalian kidneys, highlight studies of CAKUT-causing genes in Xenopus and describe how common kidney diseases have been modeled successfully in this model organism. Additionally, we discuss several molecular pathways associated with kidney disease that have been studied in Xenopus and demonstrate why it is a useful model for studying human kidney diseases.

Summary: Understanding how congenital kidney diseases arise is imperative to their treatment. Using Xenopus as a model will aid in elucidating kidney development and congenital kidney diseases.

Encounters across networks: Windows into principles of genomic regulation

Gene regulatory networks account for the ability of the genome to program development in complex multi-cellular organisms. Such networks are based on principles of gene regulation by combinations of transcription factors that bind to specific cis-regulatory DNA sites to activate transcription. These cis-regulatory regions mediate logic processing at each network node, enabling progressive increases in organismal complexity with development. Gene regulatory network explanations of development have been shown to account for patterning and cell type diversification in fly and sea urchin embryonic systems, where networks are characterized by fast coupling between transcriptional inputs and changes in target gene transcription rates, and crucial cis-regulatory elements are concentrated relatively close to the protein coding sequences of the target genes, thus facilitating their identification. Stem cell-based development in post-embryonic mammalian systems also depends on gene networks, but differs from the fly and sea urchin systems. First, the number of regulatory elements per gene and the distances between regulatory elements and the genes they control are considerably larger, forcing searches via genome-wide transcription factor binding surveys rather than functional assays. Second, the intrinsic timing of network state transitions can be slowed considerably by the need to undo stem-cell chromatin configurations, which presumably add stability to stem-cell states but retard responses to transcription factor changes during differentiation. The dispersed, partially redundant cis-regulatory systems controlling gene expression and the slow state transition kinetics in these systems already reveal new insights and opportunities to extend understanding of the repertoire of gene networks and regulatory system logic.

A gene regulatory program controlling early Xenopus mesendoderm formation: Network conservation and motifs

Germ layer formation is among the earliest differentiation events in metazoan embryos. In triploblasts, three germ layers are formed, among which the endoderm gives rise to the epithelial lining of the gut tube and associated organs including the liver, pancreas and lungs. In frogs (Xenopus), where early germ layer formation has been studied extensively, the process of endoderm specification involves the interplay of dozens of transcription factors. Here, we review the interactions between these factors, summarized in a transcriptional gene regulatory network (GRN). We highlight regulatory connections conserved between frog, fish, mouse, and human endodermal lineages. Especially prominent is the conserved role and regulatory targets of the Nodal signaling pathway and the T-box transcription factors, Vegt and Eomes. Additionally, we highlight network topologies and motifs, and speculate on their possible roles in development.

A catalog of Xenopus tropicalis transcription factors and their regional expression in the early gastrula stage embryo

Gene regulatory networks (GRNs) involve highly combinatorial interactions between transcription factors and short sequence motifs in cis-regulatory modules of target genes to control cellular phenotypes. The GRNs specifying most cell types are largely unknown and are the subject of wide interest. A catalog of transcription factors is a valuable tool toward obtaining a deeper understanding of the role of these critical effectors in any biological setting. Here we present a comprehensive catalog of the transcription factors for the diploid frog Xenopus tropicalis. We identify 1235 genes encoding DNA-binding transcription factors, comparable to the numbers found in typical mammalian species. In detail, the repertoire of X. tropicalis transcription factor genes is nearly identical to human and mouse, with the exception of zinc finger family members, and a small number of species/lineage-specific gene duplications and losses relative to the mammalian repertoires. We applied this resource to the identification of transcription factors differentially expressed in the early gastrula stage embryo. We find transcription factor enrichment in Spemann's organizer, the ventral mesoderm, ectoderm and endoderm, and report 218 TFs that show regionalized expression patterns at this stage. Many of these have not been previously reported as expressed in the early embryo, suggesting thus far unappreciated roles for many transcription factors in the GRNs regulating early development. We expect our transcription factor catalog will facilitate myriad studies using Xenopus as a model system to understand basic biology and human disease.

Using Xenopus to study genetic kidney diseases

Modern sequencing technology is revolutionizing our knowledge of inherited kidney disease. However, the molecular role of genes affected by the rapidly rising number of identified mutations is lagging behind. Xenopus is a highly useful, but underutilized model organism with unique properties excellently suited to decipher the molecular mechanisms of kidney development and disease. The embryonic kidney (pronephros) can be manipulated on only one side of the animal and its formation observed directly through the translucent skin. The moderate evolutionary distance between Xenopus and humans is a huge advantage for studying basic principles of kidney development, but still allows us to analyze the function of disease related genes. Optogenetic manipulations and genome editing by CRISPR/Cas are exciting additions to the toolbox for disease modelling and will facilitate the use of Xenopus in translational research. Therefore, the future of Xenopus in kidney research is bright.

Orchestrating liver development

The liver is a central regulator of metabolism, and liver failure thus constitutes a major health burden. Understanding how this complex organ develops during embryogenesis will yield insights into how liver regeneration can be promoted and how functional liver replacement tissue can be engineered. Recent studies of animal models have identified key signaling pathways and complex tissue interactions that progressively generate liver progenitor cells, differentiated lineages and functional tissues. In addition, progress in understanding how these cells interact, and how transcriptional and signaling programs precisely coordinate liver development, has begun to elucidate the molecular mechanisms underlying this complexity. Here, we review the lineage relationships, signaling pathways and transcriptional programs that orchestrate hepatogenesis.

Generating human intestinal tissues from pluripotent stem cells to study development and disease

As one of the largest and most functionally complex organs of the human body, the intestines are primarily responsible for the breakdown and uptake of macromolecules from the lumen and the subsequent excretion of waste from the body. However, the intestine is also an endocrine organ, regulating digestion, metabolism, and feeding behavior. Intricate neuronal, lymphatic, immune, and vascular systems are integrated into the intestine and are required for its digestive and endocrine functions. In addition, the gut houses an extensive population of microbes that play roles in digestion, global metabolism, barrier function, and host-parasite interactions. With such an extensive array of cell types working and performing in one essential organ, derivation of functional intestinal tissues from human pluripotent stem cells (PSCs) represents a significant challenge. Here we will discuss the intricate developmental processes and cell types that are required for assembly of this highly complex organ and how embryonic processes, particularly morphogenesis, have been harnessed to direct differentiation of PSCs into 3-dimensional human intestinal organoids (HIOs) in vitro. We will further describe current uses of HIOs in development and disease research and how additional tissue complexity might be engineered into HIOs for better functionality and disease modeling.

Sox17 regulates insulin secretion in the normal and pathologic mouse β cell

SOX17 is a key transcriptional regulator that can act by regulating other transcription factors including HNF1β and FOXA2, which are known to regulate postnatal β cell function. Given this, we investigated the role of SOX17 in the developing and postnatal pancreas and found a novel role for SOX17 in regulating insulin secretion. Deletion of the Sox17 gene in the pancreas (Sox17-paLOF) had no observable impact on pancreas development. However, Sox17-paLOF mice had higher islet proinsulin protein content, abnormal trafficking of proinsulin, and dilated secretory organelles suggesting that Sox17-paLOF adult mice are prediabetic. Consistant with this, Sox17-paLOF mice were more susceptible to aged-related and high fat diet-induced hyperglycemia and diabetes. Overexpression of Sox17 in mature β cells using Ins2-rtTA driver mice resulted in precocious secretion of proinsulin. Transcriptionally, SOX17 appears to broadly regulate secretory networks since a 24-hour pulse of SOX17 expression resulted in global transcriptional changes in factors that regulate hormone transport and secretion. Lastly, transient SOX17 overexpression was able to reverse the insulin secretory defects observed in MODY4 animals and restored euglycemia. Together, these data demonstrate a critical new role for SOX17 in regulating insulin trafficking and secretion and that modulation of Sox17-regulated pathways might be used therapeutically to improve cell function in the context of diabetes.

High-resolution analysis of gene activity during the Xenopus mid-blastula transition

The Xenopus mid-blastula transition (MBT) marks the onset of large-scale zygotic transcription, as well as an increase in cell cycle length and a loss of synchronous cell divisions. Little is known about what triggers the activation of transcription or how newly expressed genes interact with each other. Here, we use high-resolution expression profiling to identify three waves of gene activity: a post-fertilisation wave involving polyadenylation of maternal transcripts; a broad wave of zygotic transcription detectable as early as the seventh cleavage and extending beyond the MBT at the twelfth cleavage; and a shorter post-MBT wave of transcription that becomes apparent as development proceeds. Our studies have also allowed us to define a set of maternal mRNAs that are deadenylated shortly after fertilisation, and are likely to be degraded thereafter. Experimental analysis indicates that the polyadenylation of maternal transcripts is necessary for the establishment of proper levels of zygotic transcription at the MBT, and that genes activated in the second wave of expression, including Brachyury and Mixer, contribute to the regulation of genes expressed in the third. Together, our high-resolution time series and experimental studies have yielded a deeper understanding of the temporal organisation of gene regulatory networks in the early Xenopus embryo.

How to make an intestine

With the high prevalence of gastrointestinal disorders, there is great interest in establishing in vitro models of human intestinal disease and in developing drug-screening platforms that more accurately represent the complex physiology of the intestine. We will review how recent advances in developmental and stem cell biology have made it possible to generate complex, three-dimensional, human intestinal tissues in vitro through directed differentiation of human pluripotent stem cells. These are currently being used to study human development, genetic forms of disease, intestinal pathogens, metabolic disease and cancer.

Inference of the Xenopus tropicalis embryonic regulatory network and spatial gene expression patterns

BACKGROUND

During embryogenesis, signaling molecules produced by one cell population direct gene regulatory changes in neighboring cells and influence their developmental fates and spatial organization. One of the earliest events in the development of the vertebrate embryo is the establishment of three germ layers, consisting of the ectoderm, mesoderm and endoderm. Attempts to measure gene expression in vivo in different germ layers and cell types are typically complicated by the heterogeneity of cell types within biological samples (i.e., embryos), as the responses of individual cell types are intermingled into an aggregate observation of heterogeneous cell types. Here, we propose a novel method to elucidate gene regulatory circuits from these aggregate measurements in embryos of the frog Xenopus tropicalis using gene network inference algorithms and then test the ability of the inferred networks to predict spatial gene expression patterns.

RESULTS

We use two inference models with different underlying assumptions that incorporate existing network information, an ODE model for steady-state data and a Markov model for time series data, and contrast the performance of the two models. We apply our method to both control and knockdown embryos at multiple time points to reconstruct the core mesoderm and endoderm regulatory circuits. Those inferred networks are then used in combination with known dorsal-ventral spatial expression patterns of a subset of genes to predict spatial expression patterns for other genes. Both models are able to predict spatial expression patterns for some of the core mesoderm and endoderm genes, but interestingly of different gene subsets, suggesting that neither model is sufficient to recapitulate all of the spatial patterns, yet they are complementary for the patterns that they do capture.

CONCLUSION

The presented methodology of gene network inference combined with spatial pattern prediction provides an additional layer of validation to elucidate the regulatory circuits controlling the spatial-temporal dynamics in embryonic development.

Multicolor labeling in developmental gene regulatory network analysis

The sea urchin embryo is an important model system for developmental gene regulatory network (GRN) analysis. This chapter describes the use of multicolor fluorescent in situ hybridization (FISH) as well as a combination of FISH and immunohistochemistry in sea urchin embryonic GRN studies. The methods presented here can be applied to a variety of experimental settings where accurate spatial resolution of multiple gene products is required for constructing a developmental GRN.

A maternally established SoxB1/SoxF axis is a conserved feature of chordate germ layer patterning

Despite deep evolutionary roots in the metazoa, the gene regulatory network driving germ layer specification is surprisingly labile both between and within phyla. In Xenopus laevis, SoxB1- and SoxF-type transcription factors are intimately involved in germ-layer specification, in part through their regulation of Nodal signaling. However, it is unclear if X. laevis is representative of the ancestral vertebrate condition, as the precise roles of SoxF and SoxB1 in germ-layer specification vary among vertebrates, and there is no evidence that SoxF mediates germ-layer specification in any invertebrate. To better understand the evolution of germ-layer specification in the vertebrate lineage, we analyzed the expression of soxB1 and soxF genes in embryos and larvae of the basal vertebrate lamprey, and the basal chordate amphioxus. We find that both species maternally deposit soxB1 mRNA in the animal pole, soxF mRNA in the vegetal hemisphere, and zygotically express soxB1 and soxF throughout nascent ectoderm and mesendoderm, respectively. We also find that soxF is excluded from the vegetalmost blastomeres in lamprey and that, in contrast to vertebrates, amphioxus does not express soxF in the oral epithelium. In the context of recent work, our results suggest that a maternally established animal/vegetal Sox axis is a deeply conserved feature of chordate development that predates the role of Nodal in vertebrate germ-layer specification. Furthermore, exclusion of this axis from the vegetal pole in lamprey is consistent with the presence of an extraembryonic yolk mass, as has been previously proposed. Finally, conserved expression of SoxF in the forming mouth across the vertebrates, but not in amphioxus, lends support to the idea that the larval amphioxus mouth is nonhomologous to the vertebrate mouth.

Prolonged FGF signaling is necessary for lung and liver induction in Xenopus

BACKGROUND

FGF signaling plays numerous roles during organogenesis of the embryonic gut tube. Mouse explant studies suggest that different thresholds of FGF signaling from the cardiogenic mesoderm induce lung, liver, and pancreas lineages from the ventral foregut progenitor cells. The mechanisms that regulate FGF dose in vivo are unknown. Here we use Xenopus embryos to examine the hypothesis that a prolonged duration of FGF signaling from the mesoderm is required to induce foregut organs.

RESULTS

We show that both mesoderm and FGF signaling are required for liver and lung development in Xenopus; formally demonstrating that this important step in organ induction is conserved with other vertebrate species. Prolonged contact with the mesoderm and persistent FGF signaling through both MEK and PI3K over an extended period of time are required for liver and lung specification. Inhibition of FGF signaling results in reduced liver and lung development, with a modest expansion of the pancreas/duodenum progenitor domain. Hyper-activation of FGF signaling has the opposite effect expanding liver and lung gene expression and repressing pancreatic markers. We show that FGF signaling is cell autonomously required in the endoderm and that a dominant negative FGF receptor decreases the ability of ventral foregut progenitor cells to contribute to the lung and liver buds.

CONCLUSIONS

These results suggest that the liver and lungs are specified at progressively later times in development requiring mesoderm contact for different lengths of time. Our data suggest that this is achieved at least in part through prolonged FGF signaling. In addition to providing a foundation for further mechanistic studies on foregut organogenesis using the experimental advantages of the Xenopus system, these data have implications for the directed differentiation of stem cells into foregut lineages.

Molecular pathways controlling pancreas induction

Recent advances in generating pancreatic cell types from human pluripotent stem cells has depended on our knowledge of the developmental processes that regulate pancreas development in vivo. The developmental events between gastrulation and formation of the embryonic pancreatic primordia are both rapid and dynamic and studies in frog, fish, chick, and mouse have identified the molecular basis of how the pancreas develops from multipotent endoderm progenitors. Here, we review the current status of our understanding of molecular mechanisms that control endoderm formation, endoderm patterning, and pancreas specification and highlight how these discoveries have allowed for the development of robust methods to generate pancreatic cells from human pluripotent stem cells.

Xenopus staufen2 is required for anterior endodermal organ formation

Defining the regulatory molecular networks involved in patterning the developing anterior endoderm is essential to understand how the pancreas, liver, stomach, and duodenum are discretely specified from each other. In this study, we analyzed the expression and function of the double-stranded RNA-binding protein Staufen2 in Xenopus laevis endoderm. We found that staufen2 was broadly expressed within the developing endoderm beginning at gastrulation becoming localized to the anterior endoderm at later stages. Through morpholino-mediated knockdown, we demonstrate that Staufen2 function is required for proper formation of the stomach, liver, and pancreas. We define that its function is required during gastrulation for proper patterning of the dorsal-ventral axis and that it acts to regulate expression of BMP signaling components.

The RNA-binding protein XSeb4R regulates maternal Sox3 at the posttranscriptional level during maternal-zygotic transition in Xenopus

The maternal-zygotic transition (MZT) is an embryonic event that overlaps with and plays key roles in primary germ layer specification in vertebrates. During MZT, maternally supplied mRNAs are degraded while zygotic transcripts are synthesized to either reinforce the already specified cell fate or to trigger new cell identity. Here, we show that forced expression of the RNA-binding protein, XSeb4R, in animal pole blastomeres of Xenopus embryos, inappropriately stabilizes transcripts there, including maternal Sox3. This leads to the impaired ability of the ectodermal progenitors to respond to factors regulating brain patterning and their eventual loss by apoptosis. XSeb4R protein binds specifically to the 3'UTR of Sox3 mRNA. XSeb4R gain-of-function in ectodermal explants reveals increased stability of the maternal Sox3 transcripts, associated with a robust Sox3 protein production. Conversely, whereas XSeb4R depletion abolishes VegT expression, the amount of the maternal Sox3 mRNA is rather increased but without augmentation in the amount of Sox3 protein. Moreover, XSeb4R protein knockdown leads to the modification of the ectoderm-mesoderm boundary, marked by expanded/shifted expression of the mesodermal marker genes such as Xbra and Apod, followed by an expression inhibition of Epi. K., an ectodermal marker. Overall, our data suggest XSeb4R as a novel player in gene expression regulation, acting at the posttranscriptional level during ectoderm specification in Xenopus.

Molecular evolution of the transferrin family and associated receptors

BACKGROUND

In vertebrates, serum transferrins are essential iron transporters that have bind and release Fe(III) in response to receptor binding and changes in pH. Some family members such as lactoferrin and melanotransferrin can also bind iron while others have lost this ability and have gained other functions, e.g., inhibitor of carbonic anhydrase (mammals), saxiphilin (frogs) and otolith matrix protein 1 (fish).

SCOPE OF REVIEW

This article provides an overview of the known transferrin family members and their associated receptors and interacting partners.

MAJOR CONCLUSIONS

The number of transferrin genes has proliferated as a result of multiple duplication events, and the resulting paralogs have developed a wide array of new functions. Some homologs in the most primitive metazoan groups resemble both serum and melanotransferrins, but the major yolk proteins show considerable divergence from the rest of the family. Among the transferrin receptors, the lack of TFR2 in birds and reptiles, and the lack of any TFR homologs among the insects draw attention to the differences in iron transport and regulation in those groups.

GENERAL SIGNIFICANCE

The transferrin family members are important because of their clinical significance, interesting biochemical properties, and evolutionary history. More work is needed to better understand the functions and evolution of the non-vertebrate family members. This article is part of a Special Issue entitled Molecular Mechanisms of Iron Transport and Disorders.

Activation of neural and pluripotent stem cell signatures correlates with increased malignancy in human glioma

The presence of stem cell characteristics in glioma cells raises the possibility that mechanisms promoting the maintenance and self-renewal of tissue specific stem cells have a similar function in tumor cells. Here we characterized human gliomas of various malignancy grades for the expression of stem cell regulatory proteins. We show that cells in high grade glioma co-express an array of markers defining neural stem cells (NSCs) and that these proteins can fulfill similar functions in tumor cells as in NSCs. However, in contrast to NSCs glioma cells co-express neural proteins together with pluripotent stem cell markers, including the transcription factors Oct4, Sox2, Nanog and Klf4. In line with this finding, in high grade gliomas mesodermal- and endodermal-specific transcription factors were detected together with neural proteins, a combination of lineage markers not normally present in the central nervous system. Persistent presence of pluripotent stem cell traits could only be detected in solid tumors, and observations based on in vitro studies and xenograft transplantations in mice imply that this presence is dependent on the combined activity of intrinsic and extrinsic regulatory cues. Together these results demonstrate a general deregulated expression of neural and pluripotent stem cell traits in malignant human gliomas, and indicate that stem cell regulatory factors may provide significant targets for therapeutic strategies.

Microarray identification of novel downstream targets of FoxD4L1/D5, a critical component of the neural ectodermal transcriptional network

FoxD4L1/D5 is a forkhead transcription factor that functions as both a transcriptional activator and repressor. FoxD4L1/D5 acts upstream of several other neural transcription factors to maintain neural fate, regulate neural plate patterning, and delay the expression of neural differentiation factors. To identify a more complete list of downstream genes that participate in these earliest steps of neural ectodermal development, we carried out a microarray analysis comparing gene expression in control animal cap ectodermal explants (ACs), which will form epidermis, to that in FoxD4L1/D5-expressing ACs. Forty-four genes were tested for validation by RT-PCR of ACs and/or in situ hybridization assays in embryos; 86% of those genes up-regulated and 100% of those genes down-regulated in the microarray were altered accordingly in one of these independent assays. Eleven of these 44 genes are of unknown function, and we provide herein their developmental expression patterns to begin to reveal their roles in ectodermal development.

Sox17 and chordin are required for formation of Kupffer's vesicle and left-right asymmetry determination in zebrafish

Kupffer's vesicle (KV), a ciliated fluid-filled sphere in the zebrafish embryo with a critical role in laterality determination, is derived from a group of superficial cells in the organizer region of the gastrula named the dorsal forerunner cells (DFC). We have examined the role of the expression of sox17 and chordin (chd) in the DFC in KV formation and laterality determination. Whereas sox17 was known to be expressed in DFC, its function in these cells was not studied before. Further, expression of chd in these cells has not been reported previously. Targeted knockdown of Sox17 and Chd in DFC led to aberrant Left-Right (L-R) asymmetry establishment, as visualized by the expression of southpaw and lefty, and heart and pancreas placement in the embryo. These defects correlated with the formation of small KVs with apparently diminished cilia, consistent with the known requirement for ciliary function in the laterality organ for the establishment of L-R asymmetry.

Klf5 regulates lineage formation in the pre-implantation mouse embryo

Kruppel-like transcription factors (Klfs) are essential for the induction and maintenance of pluripotency of embryonic stem cells (ESCs), yet little is known about their roles in establishing the three lineages of the pre-implantation embryo. Here, we show that Klf5 is required for the formation of the trophectoderm (TE) and the inner cell mass (ICM), and for repressing primitive endoderm (PE) development. Although cell polarity appeared normal, Klf5 mutant embryos arrested at the blastocyst stage and failed to hatch due to defective TE development. Klf5 acted cell-autonomously in the TE, downstream of Fgf4 and upstream of Cdx2, Eomes and Krt8. In the ICM, loss of Klf5 resulted in reduced expression of pluripotency markers Oct4 and Nanog, but led to increased Sox17 expression in the PE, suggesting that Klf5 suppresses the PE lineage. Consistent with this, overexpression of Klf5 in transgenic embryos was sufficient to suppress the Sox17(+) PE lineage in the ICM. Klf5 overexpression led to a dose-dependent decrease in Sox17 promoter activity in reporter assays in cultured cells. Moreover, in chimeric embryos, Klf5(-/-) cells preferentially contributed to the Sox17(+) PE lineage and Cdx2 expression was not rescued in Klf5(-/-) outer cells. Finally, outgrowths from Klf5(-/-) embryos failed to form an ICM/pluripotent colony, had very few Oct4(+) or Cdx2(+) cells, but showed an increase in the percentage of Sox17(+) PE cells. These findings demonstrate that Klf5 is a dynamic regulator of all three lineages in the pre-implantation embryo by promoting the TE and epiblast lineages while suppressing the PE lineage.

Interactions between SOX factors and Wnt/beta-catenin signaling in development and disease

The SOX family of transcription factors have emerged as modulators of canonical Wnt/beta-catenin signaling in diverse development and disease contexts. There are over 20 SOX proteins encoded in the vertebrate genome and recent evidence suggests that many of these can physically interact with beta-catenin and modulate the transcription of Wnt-target genes. The precise mechanisms by which SOX proteins regulate beta-catenin/TCF activity are still being resolved and there is evidence to support a number of models including: protein-protein interactions, the binding of SOX factors to Wnt-target gene promoters, the recruitment of co-repressors or co-activators, modulation of protein stability, and nuclear translocation. In some contexts, Wnt signaling also regulates SOX expression resulting in feedback regulatory loops that fine-tune cellular responses to beta-catenin/TCF activity. In this review, we summarize the examples of Sox-Wnt interactions and examine the underlying mechanisms of this potentially widespread and underappreciated mode of Wnt-regulation.

Multi-scale mechanics from molecules to morphogenesis

Dynamic mechanical processes shape the embryo and organs during development. Little is understood about the basic physics of these processes, what forces are generated, or how tissues resist or guide those forces during morphogenesis. This review offers an outline of some of the basic principles of biomechanics, provides working examples of biomechanical analyses of developing embryos, and reviews the role of structural proteins in establishing and maintaining the mechanical properties of embryonic tissues. Drawing on examples we highlight the importance of investigating mechanics at multiple scales from milliseconds to hours and from individual molecules to whole embryos. Lastly, we pose a series of questions that will need to be addressed if we are to understand the larger integration of molecular and physical mechanical processes during morphogenesis and organogenesis.

Sox17 regulates organ lineage segregation of ventral foregut progenitor cells

The ventral pancreas, biliary system, and liver arise from the posterior ventral foregut, but the cell-intrinsic pathway by which these organ lineages are separated is not known. Here we show that the extrahepatobiliary system shares a common origin with the ventral pancreas and not the liver, as previously thought. These pancreatobiliary progenitor cells coexpress the transcription factors PDX1 and SOX17 at E8.5 and their segregation into a PDX1+ ventral pancreas and a SOX17+ biliary primordium is Sox17-dependent. Deletion of Sox17 at E8.5 results in the loss of biliary structures and ectopic pancreatic tissue in the liver bud and common duct, while Sox17 overexpression suppresses pancreas development and promotes ectopic biliary-like tissue throughout the PDX1+ domain. Restricting SOX17+ biliary progenitor cells to the ventral region of the gut requires the notch effector Hes1. Our results highlight the role of Sox17 and Hes1 in patterning and morphogenetic segregation of ventral foregut lineages.

Sox17 facilitates the differentiation of mouse embryonic stem cells into primitive and definitive endoderm in vitro

The Sox family of HMG (high mobility group)-box transcription factors are highly conserved in vertebrates. Sox members are involved in various developmental processes. Among them Sox17 has been demonstrated to function as an endoderm determinant in zebrafish and Xenopus, respectively. However, little is known about the role of Sox17 in mouse embryonic stem cell (ESC) differentiation. In our research, we investigated the effect of Sox17 on mouse ESC and embryoid body (EB) differentiation. The results demonstrated that Sox17 overexpression upregulated a set of endoderm-specific gene markers, suggesting that Sox17 overexpression induced an ESC differentiation program towards both primitive and definitive endoderm. We believe this finding brings new insights into the understanding of ESC differentiation and the organogenesis of endodermal derivatives.

A distinct H2A.X isoform is enriched in Xenopus laevis eggs and early embryos and is phosphorylated in the absence of a checkpoint

Histone H2A.X is an H2A variant present in multicellular organisms that is specifically phosphorylated on the serine in the C-terminal consensus sequence, canonically "SQEY," in response to DNA damage. We have recently shown the significance of phosphorylation of the penultimate tyrosine for maintenance and processing of the DNA damage response in mammalian cells. Here, we report the identification of distinct H2A.X variants in the eggs and early embryos of the frog Xenopus laevis that contain a C-terminal SQEF, among other changes; we have denoted these proteins as "H2A.X-F." H2A.X-F is present only in late-staged oocytes, eggs, and premidblastula transition embryos and is not present in somatic cells. Similar unannotated isoforms were identified in other rapidly developing aquatic species, such as Xenopus tropicalis, goldfish, and zebrafish, and in Arabidopsis and chickpea. Furthermore, we demonstrate by mass spectrometry and phospho-specific antibodies that H2A.X-F is phosphorylated in the absence of exogenous DNA damage, in both actively dividing, unperturbed embryos and cell-free egg extract in the absence and presence of DNA damage and S-phase checkpoint conditions. We propose that this isoform may be involved in modulating the cellular response to the rapid early cell cycles in externally developing species.

Vertebrate endoderm development and organ formation

The endoderm germ layer contributes to the respiratory and gastrointestinal tracts and to all of their associated organs. Over the past decade, studies in vertebrate model organisms, including frog, fish, chick, and mouse, have greatly enhanced our understanding of the molecular basis of endoderm organ development. We review this progress with a focus on early stages of endoderm organogenesis including endoderm formation, gut tube morphogenesis and patterning, and organ specification. Lastly, we discuss how developmental mechanisms that regulate endoderm organogenesis are used to direct differentiation of embryonic stem cells into specific adult cell types, which function to alleviate disease symptoms in animal models.

The RNA-binding protein XSeb4R: a positive regulator of VegT mRNA stability and translation that is required for germ layer formation in Xenopus

VegT represents a localized maternal determinant essentially required for endoderm formation in Xenopus. Here, we report on the identification of the RNA-binding protein XSeb4R as a positive regulator of VegT. XSeb4R interacts directly with the 3'-untranslated region of VegT mRNA, stabilizes it, and stimulates translation. Ablation of XSeb4R activity results in impairment of endoderm and mesoderm formation, while ectopic expression of XSeb4R in ectodermal cells induces endodermal and mesodermal gene expression. These observations unravel a novel mode of VegT regulation at the post-transcriptional level that is essential for germ layer formation in Xenopus.

Single-species microarrays and comparative transcriptomics

BACKGROUND

Prefabricated expression microarrays are currently available for only a few species but methods have been proposed to extend their application to comparisons between divergent genomes.

METHODOLOGY/PRINCIPAL FINDINGS

Here we demonstrate that the hybridization intensity of genomic DNA is a poor basis on which to select unbiased probes on Affymetrix expression arrays for studies of comparative transcriptomics, and that doing so produces spurious results. We used the Affymetrix Xenopus laevis microarray to evaluate expression divergence between X. laevis, X. borealis, and their F1 hybrids. When data are analyzed with probes that interrogate only sequences with confirmed identity in both species, we recover results that differ substantially analyses that use genomic DNA hybridizations to select probes.

CONCLUSIONS/SIGNIFICANCE

Our findings have implications for the experimental design of comparative expression studies that use single-species microarrays, and for our understanding of divergent expression in hybrid clawed frogs. These findings also highlight important limitations of single-species microarrays for studies of comparative transcriptomics of polyploid species.

The core endodermal gene network of vertebrates: combining developmental precision with evolutionary flexibility

Embryonic development combines paradoxical properties: it has great precision, it is usually conducted at breakneck speed and it is flexible on relatively short evolutionary time scales, particularly at early stages. While these features appear mutually exclusive, we consider how they may be reconciled by the properties of key early regulatory networks. We illustrate these ideas with the network that controls development of endoderm progenitors. We argue that this network enables precision because of its intrinsic stability, self propagation and dependence on signalling. The network enables high developmental speed because it is rapidly established by maternal inputs at multiple points. In turn these properties confer flexibility on an evolutionary time scale because they can be initiated in many ways, while buffering essential progenitor cell populations against changes in their embryonic environment on both evolutionary and developmental time scales. Although stable, these networks must be capable of rapid dissolution as cell differentiation progresses. While we focus on the core early endodermal network of vertebrates, we argue that these properties are likely to be general in early embryonic stem cell populations, such as mammalian ES cells.

Identification of transcripts with enriched expression in the developing and adult pancreas

BACKGROUND

Despite recent advances, the transcriptional hierarchy driving pancreas organogenesis remains largely unknown, in part due to the paucity of comprehensive analyses. To address this deficit we generated ten SAGE libraries from the developing murine pancreas spanning Theiler stages 17-26, making use of available Pdx1 enhanced green fluorescent protein (EGFP) and Neurog3 EGFP reporter strains, as well as tissue from adult islets and ducts.

RESULTS

We used a specificity metric to identify 2,536 tags with pancreas-enriched expression compared to 195 other mouse SAGE libraries. We subsequently grouped co-expressed transcripts with differential expression during pancreas development using K-means clustering. We validated the clusters first using quantitative real time PCR and then by analyzing the Theiler stage 22 pancreas in situ hybridization staining patterns of over 600 of the identified genes using the GenePaint database. These were then categorized into one of the five expression domains within the developing pancreas. Based on these results we identified a cascade of transcriptional regulators expressed in the endocrine pancreas lineage and, from this, we developed a predictive regulatory network describing beta-cell development.

CONCLUSION

Taken together, this work provides evidence that the SAGE libraries generated here are a valuable resource for continuing to elucidate the molecular mechanisms regulating pancreas development. Furthermore, our studies provide a comprehensive analysis of pancreas development, and insights into the regulatory networks driving this process are revealed.

The role of FGF signaling in the establishment and maintenance of mesodermal gene expression in Xenopus

FGF signaling is important for the formation of mesoderm in vertebrates, and when it is perturbed in Xenopus, most trunk and tail mesoderm fails to form. Here we have further dissected the activities of FGF in patterning the embryo by addressing its inductive and maintenance roles. We show that FGF signaling is necessary for the establishment of xbra expression in addition to its well-characterized role in maintaining xbra expression. The role of FGF signaling in organizer formation is not clear in Xenopus. We find that FGF signaling is essential for the initial specification of paraxial mesoderm but not for activation of several pan-mesodermal and most organizer genes; however, early FGF signaling is necessary for the maintenance of organizer gene expression into the neurula stage. Inhibition of FGF signaling prevents VegT activation of specific mesodermal transcripts. These findings illuminate how FGF signaling contributes to the establishment of distinct types of mesoderm.

The efficiency of Xenopus primordial germ cell migration depends on the germplasm mRNA encoding the PDZ domain protein Grip2

A microarray analysis of vegetal pole sequences in the egg and early Xenopus laevis embryo identified Unigene Xl.14891 as a vegetally localized RNA. Analysis of the Xenopus tropicalis genome showed this Unigene to be localized near the 3' end of the Grip2 (glutamate receptor interacting protein 2) transcription unit. RACE showed that the Unigene represented the 3' UTR of Grip2 mRNA. Grip2 mRNA is present in the mitochondrial cloud of late pre-vitellogenic oocytes and then in the germplasm through oogenesis and early development until tailbud tadpole stages. Interference with Grip2 mRNA translation using two antisense morpholino oligos (MOs) impairs primordial germ cell (PGC) migration to the germinal ridges. Both MOs also inhibit swimming movements of the tailbud tadpole, known to involve glutamate receptors. We conclude that Grip2 has several functions in the embryo, including enabling efficient PGC migration.

Duplicate gene evolution and expression in the wake of vertebrate allopolyploidization

Background

The mechanism by which duplicate genes originate – whether by duplication of a whole genome or of a genomic segment – influences their genetic fates. To study events that trigger duplicate gene persistence after whole genome duplication in vertebrates, we have analyzed molecular evolution and expression of hundreds of persistent duplicate gene pairs in allopolyploid clawed frogs (Xenopus and Silurana). We collected comparative data that allowed us to tease apart the molecular events that occurred soon after duplication from those that occurred later on. We also quantified expression profile divergence of hundreds of paralogs during development and in different tissues.

Results

Our analyses indicate that persistent duplicates generated by allopolyploidization are subjected to strong purifying selection soon after duplication. The level of purifying selection is relaxed compared to a singleton ortholog, but not significantly variable over a period spanning about 40 million years. Despite persistent functional constraints, however, analysis of paralogous expression profiles indicates that quantitative aspects of their expression diverged substantially during this period.

Conclusion

These results offer clues into how vertebrate transcriptomes are sculpted in the wake of whole genome duplication (WGD), such as those that occurred in our early ancestors. That functional constraints were relaxed relative to a singleton ortholog but not significantly different in the early compared to the later stage of duplicate gene evolution suggests that the timescale for a return to pre-duplication levels is drawn out over tens of millions of years – beyond the age of these tetraploid species. Quantitative expression divergence can occur soon after WGD and with a magnitude that is not correlated with the rate of protein sequence divergence. On a coarse scale, quantitative expression divergence appears to be more prevalent than spatial and temporal expression divergence, and also faster or more frequent than other processes that operate at the protein level, such as some types of neofunctionalization.