Measuring Absolute RNA Copy Numbers at High Temporal Resolution Reveals Transcriptome Kinetics in Development

Soil-dependent responses of bacterial communities, phosphorus and carbon turnover to uranium stress in different soil ecosystems

Uranium (U) can impact microbially driven soil phosphorus (P) and carbon (C) cycling. However, the response of microbial P and C turnover to U in different soils is not well understood. Through the quantitative assay of P pools and soil organic C (SOC) quantitative assay and sequencing of 16S rRNA gene amplicons and metagenomes, we investigated the effect of U on P and C biotransformation in grassland (GL), paddy soil (PY), forest soil (FT). U (60 mg kg-1) impacted the diversity, interaction and stability of soil bacterial communities, leading to a decrease in available P (AP). Under U stress, organophosphate mineralization substantially contributed to the AP in GL and FT, whereas intracellular P metabolism dominated the AP in PY. Also, the reductive citrate cycle (rTCA cycle) promoted the content of SOC in GL, while the rTCA cycle and complex organic C degradation pathways enhanced the SOC in PY and FT, respectively. Notably, functional bacteria carrying organic C degradation genes could decompose SOC to enhance soil AP. Bacteria developed various resistance strategies to cope with U stress. This study reveals soil-dependent response of microbial P and C cycling and its ecological functions under the influence of radioactive contaminants in different soil systems.

The maternal-to-zygotic transition: reprogramming of the cytoplasm and nucleus

A fertilized egg is initially transcriptionally silent and relies on maternally provided factors to initiate development. For embryonic development to proceed, the oocyte-inherited cytoplasm and the nuclear chromatin need to be reprogrammed to create a permissive environment for zygotic genome activation (ZGA). During this maternal-to-zygotic transition (MZT), which is conserved in metazoans, transient totipotency is induced and zygotic transcription is initiated to form the blueprint for future development. Recent technological advances have enhanced our understanding of MZT regulation, revealing common themes across species and leading to new fundamental insights about transcription, mRNA decay and translation.

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 Integrated Database for Exploring Alternative Promoters in Animals

Alternative promoter (AP) events, as a major pre-transcriptional mechanism, can initiate different transcription start sites to generate distinct mRNA isoforms and regulate their expression. At present, hundreds of thousands of APs have been identified across human tissues, and a considerable number of APs have been demonstrated to be associated with complex traits and diseases. Recent researches have also proven important effects of APs on animals. However, the landscape of APs in animals has not been fully recognized. In this study, 102,349 AP profiles from 23,077 samples across 12 species were systematically characterized. We further identified tissue-specific APs and investigated trait-related promoters among various species. In addition, we analyzed the associations between APs and enhancer RNAs (eRNA)/transcription factors (TF) as a means of identifying potential regulatory factors. Integrating these findings, we finally developed Animal-APdb, a database for the searching, browsing, and downloading of information related to Animal APs. Animal-APdb is expected to serve as a valuable resource for exploring the functions and mechanisms of APs in animals.

Maternal and zygotic contributions to H3K4me1 chromatin marking during germ layer formation

An early step in triploblastic embryo differentiation is the formation of the three germ layers. Maternal pioneer transcription factors (TFs) bind to embryonic enhancers before zygotic genome activation, initiating germ layer specification. While maternal TFs' role in establishing epigenetic marks is known, how early pluripotent cells gain spatially restricted epigenetic identities remains unclear. We show that by the early gastrula stage, H3K4me1-marked regions become distinct in each germ layer, with certain chromatin regions forming high density H3K4me1 marked regions (HDRs). Genes associated with these HDRs are more robustly expressed compared to those associated with low density H3K4me1 marked regions (LDRs) in the genome. This process is driven by the sequential actions of maternal and zygotic factors. Knockdown of key maternal endodermal TFs (Otx1, Vegt and Foxh1) leads to a loss of endodermal H3K4me1 marks in endoderm, with a concurrent emergence of ectodermal and mesodermal marks, indicating a shift in chromatin state. This work highlights the importance of coordinated activities of maternal and zygotic TFs in defining the regionally-resolved and dynamic process of chromatin modification conferred by H3K4me1 in the early Xenopus embryo.

Nup107 contributes to the maternal-to-zygotic transition by preventing the premature nuclear export of pri-miR427

Emerging evidence suggests that the nuclear pore complex can have unique compositions and distinct nucleoporin functions in different cells. Here, we show that Nup107, a key component of the NPC scaffold, varies in expression over development: it is expressed at higher levels in the blastula compared to the gastrula, suggesting a crucial role before gastrulation in Xenopus. We find that depletion of Nup107 affects the differentiation of the early germ layers leading to an expansion of the ectoderm at the expense of endoderm and mesoderm. By analyzing an RNA-sequencing time course, we observed that depletion of Nup107 affects the maternal-zygotic transition by delaying the degradation of maternal transcripts that occurs as zygotic transcription begins. The transcripts are enriched in recognition sites for miR427, a conserved microRNA that destabilizes maternal transcripts including REST, which encodes a Kruppel-type zinc-finger transcription factor that we demonstrate is crucial for ectodermal cell fates. Mechanistically, we show that Nup107 is required to prevent the premature export of pri-miR427 transcript before processing. Nup107 depletion leads to the reduced production of mature miR427 and maternal transcript stabilization. We conclude that high levels of Nup107 in the early embryo are crucial for the nuclear retention and subsequent processing of pri-miR427 transcripts that is required for timely maternal RNA clearance to enable gastrulation.

Foxi2 and Sox3 are master regulators controlling ectoderm germ layer specification

In vertebrates, germ layer specification represents a critical transition where pluripotent cells acquire lineage-specific identities. We identify the maternal transcription factors Foxi2 and Sox3 to be pivotal master regulators of ectodermal germ layer specification in Xenopus. Ectopic co-expression of Foxi2 and Sox3 in prospective endodermal tissue induces the expression of ectodermal markers while suppressing mesendodermal markers. Transcriptomics analyses reveal that Foxi2 and Sox3 jointly and independently regulate hundreds of ectodermal target genes. During early cleavage stages, Foxi2 and Sox3 pre-bind to key cis-regulatory modules (CRMs), marking sites that later recruit Ep300 and facilitate H3K27ac deposition, thereby shaping the epigenetic landscape of the ectodermal genome. These CRMs are highly enriched within ectoderm-specific super-enhancers (SEs). Our findings highlight the pivotal role of ectodermal SE-associated CRMs in precise and robust ectodermal gene activation, establishing Foxi2 and Sox3 as central architects of ectodermal lineage specification.

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.

Deletion upstream of MAB21L2 highlights the importance of evolutionarily conserved non-coding sequences for eye development

Anophthalmia, microphthalmia and coloboma (AMC) comprise a spectrum of developmental eye disorders, accounting for approximately 20% of childhood visual impairment. While non-coding regulatory sequences are increasingly recognised as contributing to disease burden, characterising their impact on gene function and phenotype remains challenging. Furthermore, little is known of the nature and extent of their contribution to AMC phenotypes. We report two families with variants in or near MAB21L2, a gene where genetic variants are known to cause AMC in humans and animal models. The first proband, presenting with microphthalmia and coloboma, has a likely pathogenic missense variant (c.338 G > C; p.[Trp113Ser]), segregating within the family. The second individual, presenting with microphthalmia, carries an ~ 113.5 kb homozygous deletion 19.38 kb upstream of MAB21L2. Modelling of the deletion results in transient small lens and coloboma as well as midbrain anomalies in zebrafish, and microphthalmia and coloboma in Xenopus tropicalis. Using conservation analysis, we identify 15 non-coding conserved elements (CEs) within the deleted region, while ChIP-seq data from mouse embryonic stem cells demonstrates that two of these (CE13 and 14) bind Otx2, a protein with an established role in eye development. Targeted disruption of CE14 in Xenopus tropicalis recapitulates an ocular coloboma phenotype, supporting its role in eye development. Together, our data provides insights into regulatory mechanisms underlying eye development and highlights the importance of non-coding sequences as a source of genetic diagnoses in AMC.

Leptin signaling promotes blood vessel formation in the Xenopus tail during the embryo-larval transition

The signals that regulate peripheral blood vessel formation during development are still under investigation. The hormone leptin promotes blood vessel formation, adipose tissue establishment and expansion, tumor growth, and wound healing, but the underlying mechanisms for these actions are currently unknown. We investigated whether leptin promotes angiogenesis in the developing tail fin using embryonic transgenic xflk-1:GFP Xenopus laevis, which express a green fluorescent protein on vascular endothelial cells to mark blood vessels. We found that leptin protein is expressed in endothelial cells of developing blood vessels and that leptin treatment via injection increased phosphorylated STAT3 signaling, which is indicative of leptin activation of its receptor, in blood vessels of the larval tail fin. Leptin administration via media increased vessel length, branching, and reconnection with the cardinal vein, while decreased leptin signaling via immunoneutralization had an opposing effect on vessel development. We also observed disorganization of major vessels and microvessels of the tail fin and muscle when leptin signaling was decreased. Reduced leptin signaling lowered mRNA expression of cenpk, gpx1, and mmp9, markers for cell proliferation, antioxidation, and extracellular matrix remodeling/cell migration, respectively, in the developing tail, providing insight into three possible mechanisms underlying leptin's promotion of angiogenesis. Together these results illustrate that leptin levels are correlated with embryonic angiogenesis and that leptin coordinates multiple aspects of blood vessel growth and development, showing that leptin is an important morphogen during embryonic development.

Zbtb11 interacts with Otx2 and patterns the anterior neuroectoderm in Xenopus

The zinc finger and BTB domain-containing 11 gene (zbtb11) is expressed in the Xenopus anterior neuroectoderm, but the molecular nature of the Zbtb11 protein during embryonic development remains to be elucidated. Here, we show the role of Zbtb11 in anterior patterning of the neuroectoderm and the cooperative action with the transcription factor Otx2. Both overexpression and knockdown of zbtb11 caused similar phenotypes: expanded expression of the posterior gene gbx2 in the neural plate, and later microcephaly with reduced eyes, suggesting that a proper level of zbtb11 expression is necessary for normal patterning of the neuroectoderm, including eye formation. Co-immunoprecipitation assays showed that Zbtb11 formed a complex with itself and with a phosphomimetic and repressive form of Otx2, suggesting that Zbtb11 forms a dimer or oligomer and interacts with Otx2 in a phosphorylation-dependent manner. Reporter analysis further showed that Zbtb11 enhanced the activity of the phosphomimetic Otx2 to repress a silencer element of the posterior gene meis3. These data suggest that Zbtb11 coordinates with phosphorylated Otx2 to specify the anterior neuroectoderm by repressing posterior genes.

The Expansion of Sirtuin Gene Family in Gilthead Sea Bream (Sparus aurata)-Phylogenetic, Syntenic, and Functional Insights across the Vertebrate/Fish Lineage

The Sirtuin (SIRT1-7) family comprises seven evolutionary-conserved enzymes that couple cellular NAD availability with health, nutrition and welfare status in vertebrates. This study re-annotated the sirt3/5 branch in the gilthead sea bream, revealing three paralogues of sirt3 (sirt3.1a/sirt3.1b/sirt3.2) and two of sirt5 (sirt5a/sirt5b) in this Perciform fish. The phylogeny and synteny analyses unveiled that the Sirt3.1/Sirt3.2 dichotomy was retained in teleosts and aquatic-living Sarcopterygian after early vertebrate 2R whole genome duplication (WGD). Additionally, only certain percomorphaceae and gilthead sea bream showed a conserved tandem-duplicated synteny block involving the mammalian-clustered sirt3.1 gene (psmd13-sirt3.1a/b-drd4-cdhr5-ctsd). Conversely, the expansion of the Sirt5 branch was shaped by the teleost-specific 3R WGD. As extensively reviewed in the literature, human-orthologues (sirt3.1/sirt5a) showed a high, conserved expression in skeletal muscle that increased as development advanced. However, recent sirt3.2 and sirt5b suffered an overall muscle transcriptional silencing across life, as well as an enhanced expression on immune-relevant tissues and gills. These findings fill gaps in the ontogeny and differentiation of Sirt genes in the environmentally adaptable gilthead sea bream, becoming a good starting point to advance towards a full understanding of its neo-functionalization. The mechanisms originating from these new paralogs also open new perspectives in the study of cellular energy sensing processes in vertebrates.

The maternal-to-zygotic transition

Rapid cleavage divisions and the transition from maternal to zygotic control of gene expression are the hallmarks of early embryonic development in most species. Early development in insects, fish and amphibians is characterized by several short cell cycles with no gap phases, necessary for the rapid production of cells prior to patterning and morphogenesis. Maternal mRNAs and proteins loaded into the egg during oogenesis are essential to drive these rapid early divisions. Once the function of these maternal inputs is complete, the maternal-to-zygotic transition (MZT) marks the handover of developmental control to the gene products synthesized from the zygotic genome. The MZT requires three major events: the removal of a subset of maternal mRNAs, the initiation of zygotic transcription, and the remodeling of the cell cycle. In each species, the MZT occurs at a highly reproducible time during development due to a series of feedback mechanisms that tightly couple these three processes. Dissecting these feedback mechanisms and their spatiotemporal control will be essential to understanding the control of the MZT. In this primer, we outline the mechanisms that govern the major events of the MZT across species and highlight the role of feedback mechanisms that ensure the MZT is precisely timed and orchestrated.

Changes in repair pathways of radiation-induced DNA double-strand breaks at the midblastula transition in Xenopus embryo

Ionizing radiation (IR) causes DNA damage, particularly DNA double-strand breaks (DSBs), which have significant implications for genome stability. The major pathways of repairing DSBs are homologous recombination (HR) and nonhomologous end joining (NHEJ). However, the repair mechanism of IR-induced DSBs in embryos is not well understood, despite extensive research in somatic cells. The externally developing aquatic organism, Xenopus tropicalis, serves as a valuable model for studying embryo development. A significant increase in zygotic transcription occurs at the midblastula transition (MBT), resulting in a longer cell cycle and asynchronous cell divisions. This study examines the impact of X-ray irradiation on Xenopus embryos before and after the MBT. The findings reveal a heightened X-ray sensitivity in embryos prior to the MBT, indicating a distinct shift in the DNA repair pathway during embryo development. Importantly, we show a transition in the dominant DSB repair pathway from NHEJ to HR before and after the MBT. These results suggest that the MBT plays a crucial role in altering DSB repair mechanisms, thereby influencing the IR sensitivity of developing embryos.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Generation of translucent Xenopus tropicalis through triple knockout of pigmentation genes

Amphibians generally have three types of pigment cells, namely, melanophores (black and brown), xanthophores (yellow and red), and iridophores (iridescent). Single knockout of the tyr, slc2a7, and hps6 genes in Xenopus tropicalis results in the absence of melanophores, xanthophores, and iridophores, respectively. The generation of triple- knockout (3KO) X. tropicalis for these three genes could allow for observation of internal organs without sacrificing the animals, which would be transparent due to the absence of pigments. In this study, we generated 3KO X. tropicalis, which is one of the most widely used model amphibians, through crossing of a slc2a7 single-knockout frog with a tyr and hps6 double-knockout frog, followed by intercrossing of their offspring. The 3KO tadpoles had transparent bodies like the nop mutant and the frogs had translucent bodies. This translucency allowed us to observe the heart, lungs, stomach, liver, and digestive tract through the ventral body skin without surgery. After intravital staining, 3KO X. tropicalis showed much clearer fluorescent signals of mineralized tissues compared with the wild type. These 3KO X. tropicalis provide a useful mutant line for continuous observation of internal organs and fluorescent signals in the body. In particular, such 3KO frogs would revolutionize fluorescence monitoring in transgenic tadpoles and frogs expressing fluorescent proteins.

Deep transcriptome profiling reveals limited conservation of A-to-I RNA editing in Xenopus

BACKGROUND

Xenopus has served as a valuable model system for biomedical research over the past decades. Notably, ADAR was first detected in frog oocytes and embryos as an activity that unwinds RNA duplexes. However, the scope of A-to-I RNA editing by the ADAR enzymes in Xenopus remains underexplored.

RESULTS

Here, we identify millions of editing events in Xenopus with high accuracy and systematically map the editome across developmental stages, adult organs, and species. We report diverse spatiotemporal patterns of editing with deamination activity highest in early embryogenesis before zygotic genome activation and in the ovary. Strikingly, editing events are poorly conserved across different Xenopus species. Even sites that are detected in both X. laevis and X. tropicalis show largely divergent editing levels or developmental profiles. In protein-coding regions, only a small subset of sites that are found mostly in the brain are well conserved between frogs and mammals.

CONCLUSIONS

Collectively, our work provides fresh insights into ADAR activity in vertebrates and suggest that species-specific editing may play a role in each animal's unique physiology or environmental adaptation.

Regulation of gene expression downstream of a novel Fgf/Erk pathway during Xenopus development

Activation of Map kinase/Erk signalling downstream of fibroblast growth factor (Fgf) tyrosine kinase receptors regulates gene expression required for mesoderm induction and patterning of the anteroposterior axis during Xenopus development. We have proposed that a subset of Fgf target genes are activated in the embyo in response to inhibition of a transcriptional repressor. Here we investigate the hypothesis that Cic (Capicua), which was originally identified as a transcriptional repressor negatively regulated by receptor tyrosine kinase/Erk signalling in Drosophila, is involved in regulating Fgf target gene expression in Xenopus. We characterise Xenopus Cic and show that it is widely expressed in the embryo. Fgf overexpression or ectodermal wounding, both of which potently activate Erk, reduce Cic protein levels in embryonic cells. In keeping with our hypothesis, we show that Cic knockdown and Fgf overexpression have overlapping effects on embryo development and gene expression. Transcriptomic analysis identifies a cohort of genes that are up-regulated by Fgf overexpression and Cic knockdown. We investigate two of these genes as putative targets of the proposed Fgf/Erk/Cic axis: fos and rasl11b, which encode a leucine zipper transcription factor and a ras family GTPase, respectively. We identify Cic consensus binding sites in a highly conserved region of intron 1 in the fos gene and Cic sites in the upstream regions of several other Fgf/Cic co-regulated genes, including rasl11b. We show that expression of fos and rasl11b is blocked in the early mesoderm when Fgf and Erk signalling is inhibited. In addition, we show that fos and rasl11b expression is associated with the Fgf independent activation of Erk at the site of ectodermal wounding. Our data support a role for a Fgf/Erk/Cic axis in regulating a subset of Fgf target genes during gastrulation and is suggestive that Erk signalling is involved in regulating Cic target genes at the site of ectodermal wounding.

Hybridization led to a rewired pluripotency network in the allotetraploid Xenopus laevis

After fertilization, maternally contributed factors to the egg initiate the transition to pluripotency to give rise to embryonic stem cells, in large part by activating de novo transcription from the embryonic genome. Diverse mechanisms coordinate this transition across animals, suggesting that pervasive regulatory remodeling has shaped the earliest stages of development. Here, we show that maternal homologs of mammalian pluripotency reprogramming factors OCT4 and SOX2 divergently activate the two subgenomes of Xenopus laevis, an allotetraploid that arose from hybridization of two diploid species ~18 million years ago. Although most genes have been retained as two homeologous copies, we find that a majority of them undergo asymmetric activation in the early embryo. Chromatin accessibility profiling and CUT&RUN for modified histones and transcription factor binding reveal extensive differences in predicted enhancer architecture between the subgenomes, which likely arose through genomic disruptions as a consequence of allotetraploidy. However, comparison with diploid X. tropicalis and zebrafish shows broad conservation of embryonic gene expression levels when divergent homeolog contributions are combined, implying strong selection to maintain dosage in the core vertebrate pluripotency transcriptional program, amid genomic instability following hybridization.

Transcript accumulation rates in the early Caenorhabditis elegans embryo

Dynamic transcriptional changes are widespread in rapidly dividing developing embryos when cell fate decisions are made quickly. The Caenorhabditis elegans embryo overcomes these constraints partly through the rapid production of high levels of transcription factor mRNAs. Transcript accumulation rates for some developmental genes are known at single-cell resolution, but genome-scale measurements are lacking. We estimate zygotic mRNA accumulation rates from single-cell RNA sequencing data calibrated with single-molecule transcript imaging. Rapid transcription is common in the early C. elegans embryo with rates highest soon after zygotic transcription begins. High-rate genes are enriched for recently duplicated cell-fate regulators and share common genomic features. We identify core promoter elements associated with high rate and measure their contributions for two early endomesodermal genes, ceh-51 and sdz-31. Individual motifs modestly affect accumulation rates, suggesting multifactorial control. These results are a step toward estimating absolute transcription kinetics and understanding how transcript dosage drives developmental decisions.

RNA recording in single bacterial cells using reprogrammed tracrRNAs

Capturing an individual cell's transcriptional history is a challenge exacerbated by the functional heterogeneity of cellular communities. Here, we leverage reprogrammed tracrRNAs (Rptrs) to record selected cellular transcripts as stored DNA edits in single living bacterial cells. Rptrs are designed to base pair with sensed transcripts, converting them into guide RNAs. The guide RNAs then direct a Cas9 base editor to target an introduced DNA target. The extent of base editing can then be read in the future by sequencing. We use this approach, called TIGER (transcribed RNAs inferred by genetically encoded records), to record heterologous and endogenous transcripts in individual bacterial cells. TIGER can quantify relative expression, distinguish single-nucleotide differences, record multiple transcripts simultaneously and read out single-cell phenomena. We further apply TIGER to record metabolic bet hedging and antibiotic resistance mobilization in Escherichia coli as well as host cell invasion by Salmonella. Through RNA recording, TIGER connects current cellular states with past transcriptional states to decipher complex cellular responses in single cells.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus. Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities similar to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far particular to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

Activation of P53 pathway contributes to Xenopus hybrid inviability

Hybrid incompatibility as a kind of reproductive isolation contributes to speciation. The nucleocytoplasmic incompatibility between Xenopus tropicalis eggs and Xenopus laevis sperm (te×ls) leads to specific loss of paternal chromosomes 3L and 4L. The hybrids die before gastrulation, of which the lethal causes remain largely unclear. Here, we show that the activation of the tumor suppressor protein P53 at late blastula stage contributes to this early lethality. We find that in stage 9 embryos, P53-binding motif is the most enriched one in the up-regulated Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) peaks between te×ls and wild-type X. tropicalis controls, which correlates with an abrupt stabilization of P53 protein in te×ls hybrids at stage 9. Inhibition of P53 activity via either tp53 knockout or overexpression of a dominant-negative P53 mutant or Murine double minute 2 proto-oncogene (Mdm2), a negative regulator of P53, by mRNA injection can rescue the te×ls early lethality. Our results suggest a causal function of P53 on hybrid lethality prior to gastrulation.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus . Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities analogous to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far unique to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

Xenbase: key features and resources of the Xenopus model organism knowledgebase

Xenbase (https://www.xenbase.org/), the Xenopus model organism knowledgebase, is a web-accessible resource that integrates the diverse genomic and biological data from research on the laboratory frogs Xenopus laevis and Xenopus tropicalis. The goal of Xenbase is to accelerate discovery and empower Xenopus research, to enhance the impact of Xenopus research data, and to facilitate the dissemination of these data. Xenbase also enhances the value of Xenopus data through high-quality curation, data integration, providing bioinformatics tools optimized for Xenopus experiments, and linking Xenopus data to human data, and other model organisms. Xenbase also plays an indispensable role in making Xenopus data interoperable and accessible to the broader biomedical community in accordance with FAIR principles. Xenbase provides annotated data updates to organizations such as NCBI, UniProtKB, Ensembl, the Gene Ontology consortium, and most recently, the Alliance of Genomic Resources, a common clearing house for data from humans and model organisms. This article provides a brief overview of key and recently added features of Xenbase. New features include processing of Xenopus high-throughput sequencing data from the NCBI Gene Expression Omnibus; curation of anatomical, physiological, and expression phenotypes with the newly created Xenopus Phenotype Ontology; Xenopus Gene Ontology annotations; new anatomical drawings of the Normal Table of Xenopus development; and integration of the latest Xenopus laevis v10.1 genome annotations. Finally, we highlight areas for future development at Xenbase as we continue to support the Xenopus research community.

Intron dynamics reveal principles of gene regulation during the maternal-to-zygotic transition

The maternal-to-zygotic transition (MZT) is a conserved embryonic process in animals where developmental control shifts from the maternal to zygotic genome. A key step in this transition is zygotic transcription, and deciphering the MZT requires classifying newly transcribed genes. However, due to current technological limitations, this starting point remains a challenge for studying many species. Here, we present an alternative approach that characterizes transcriptome changes based solely on RNA-seq data. By combining intron-mapping reads and transcript-level quantification, we characterized transcriptome dynamics during the Drosophila melanogaster MZT. Our approach provides an accessible platform to investigate transcriptome dynamics that can be applied to the MZT in nonmodel organisms. In addition to classifying zygotically transcribed genes, our analysis revealed that over 300 genes express different maternal and zygotic transcript isoforms due to alternative splicing, polyadenylation, and promoter usage. The vast majority of these zygotic isoforms have the potential to be subject to different regulatory control, and over two-thirds encode different proteins. Thus, our analysis reveals an additional layer of regulation during the MZT, where new zygotic transcripts can generate additional proteome diversity.

Tissue-specific expression of carbohydrate sulfotransferases drives keratan sulfate biosynthesis in the notochord and otic vesicles of Xenopus embryos

Keratan sulfate (KS) is a glycosaminoglycan that is enriched in vertebrate cornea, cartilage, and brain. During embryonic development, highly sulfated KS (HSKS) is first detected in the developing notochord and then in otic vesicles; therefore, HSKS has been used as a molecular marker of the notochord. However, its biosynthetic pathways and functional roles in organogenesis are little known. Here, I surveyed developmental expression patterns of genes related to HSKS biosynthesis in Xenopus embryos. Of these genes, the KS chain-synthesizing glycosyltransferase genes, beta-1,3-N-acetylglucosaminyltransferase (b3gnt7) and beta-1,4-galactosyltransferase (b4galt4), are strongly expressed in the notochord and otic vesicles, but also in other tissues. In addition, their notochord expression is gradually restricted to the posterior end at the tailbud stage. In contrast, carbohydrate sulfotransferase (Chst) genes, chst2, chst3, and chst5.1, are expressed in both notochord and otic vesicles, whereas chst1, chst4/5-like, and chst7 are confined to otic vesicles. Because the substrate for Chst1 and Chst3 is galactose, while that for others is N-acetylglucosamine, combinatorial, tissue-specific expression patterns of Chst genes should be responsible for tissue-specific HSKS enrichment in embryos. As expected, loss of function of chst1 led to loss of HSKS in otic vesicles and reduction of their size. Loss of chst3 and chst5.1 resulted in HSKS loss in the notochord. These results reveal that Chst genes are critical for HSKS biosynthesis during organogenesis. Being hygroscopic, HSKS forms “water bags” in embryos to physically maintain organ structures. In terms of evolution, in ascidian embryos, b4galt and chst-like genes are also expressed in the notochord and regulate notochord morphogenesis. Furthermore, I found that a chst-like gene is also strongly expressed in the notochord of amphioxus embryos. These conserved expression patterns of Chst genes in the notochord of chordate embryos suggest that Chst is an ancestral component of the chordate notochord.

Dodecaploid Xenopus longipes provides insight into the emergence of size scaling relationships during development

Genome and cell size are strongly correlated across species^1,2,3,4,5,6 and influence physiological traits like developmental rate.^{7,8,9,10,11,12} Although size scaling features such as the nuclear-cytoplasmic (N/C) ratio are precisely maintained in adult tissues,¹³ it is unclear when during embryonic development size scaling relationships are established. Frogs of the genus Xenopus provide a model to investigate this question, since 29 extant Xenopus species vary in ploidy from 2 to 12 copies (N) of the ancestral frog genome, ranging from 20 to 108 chromosomes.^14,15 The most widely studied species, X. laevis (4N = 36) and X. tropicalis (2N = 20), scale at all levels, from body size to cellular and subcellular levels.¹⁶ Paradoxically, the rare, critically endangered dodecaploid (12N = 108) Xenopus longipes (X. longipes) is a small frog.^15,17 We observed that despite some morphological differences, X. longipes and X. laevis embryogenesis occurred with similar timing, with genome to cell size scaling emerging at the swimming tadpole stage. Across the three species, cell size was determined primarily by egg size, whereas nuclear size correlated with genome size during embryogenesis, resulting in different N/C ratios in blastulae prior to gastrulation. At the subcellular level, nuclear size correlated more strongly with genome size, whereas mitotic spindle size scaled with cell size. Our cross-species study indicates that scaling of cell size to ploidy is not due to abrupt changes in cell division timing, that different size scaling regimes occur during embryogenesis, and that the developmental program of Xenopus is remarkably consistent across a wide range of genome and egg sizes.

Developmental exposure to the pesticide malathion enhances expression of Prdm12, a regulator of nociceptor development, in Xenopus laevis

The transcription factor Prdm12 exerts important influences on the development of nociceptors, peripheral touch and pain-sensing neurons, and has been implicated in human pain sensation disorders. We examined the consequences of exposing developing Xenopus laevis embryos to the commonly used pesticide malathion on Prdm12 expression. Using qPCR and western blot analysis we observed that malathion treatment for the first six days of tadpole development significantly increased both prdm12 mRNA levels and Prdm12 protein levels compared to controls. Consequently, early exposure to this pesticide has potential to alter nociceptor development.

Membrane potential drives the exit from pluripotency and cell fate commitment via calcium and mTOR

Transitioning from pluripotency to differentiated cell fates is fundamental to both embryonic development and adult tissue homeostasis. Improving our understanding of this transition would facilitate our ability to manipulate pluripotent cells into tissues for therapeutic use. Here, we show that membrane voltage (Vm) regulates the exit from pluripotency and the onset of germ layer differentiation in the embryo, a process that affects both gastrulation and left-right patterning. By examining candidate genes of congenital heart disease and heterotaxy, we identify KCNH6, a member of the ether-a-go-go class of potassium channels that hyperpolarizes the Vm and thus limits the activation of voltage gated calcium channels, lowering intracellular calcium. In pluripotent embryonic cells, depletion of kcnh6 leads to membrane depolarization, elevation of intracellular calcium levels, and the maintenance of a pluripotent state at the expense of differentiation into ectodermal and myogenic lineages. Using high-resolution temporal transcriptome analysis, we identify the gene regulatory networks downstream of membrane depolarization and calcium signaling and discover that inhibition of the mTOR pathway transitions the pluripotent cell to a differentiated fate. By manipulating Vm using a suite of tools, we establish a bioelectric pathway that regulates pluripotency in vertebrates, including human embryonic stem cells.

Non-coding variants disrupting a tissue-specific regulatory element in HK1 cause congenital hyperinsulinism

Gene expression is tightly regulated, with many genes exhibiting cell-specific silencing when their protein product would disrupt normal cellular function1. This silencing is largely controlled by non-coding elements, and their disruption might cause human disease2. We performed gene-agnostic screening of the non-coding regions to discover new molecular causes of congenital hyperinsulinism. This identified 14 non-coding de novo variants affecting a 42-bp conserved region encompassed by a regulatory element in intron 2 of the hexokinase 1 gene (HK1). HK1 is widely expressed across all tissues except in the liver and pancreatic beta cells and is thus termed a 'disallowed gene' in these specific tissues. We demonstrated that the variants result in a loss of repression of HK1 in pancreatic beta cells, thereby causing insulin secretion and congenital hyperinsulinism. Using epigenomic data accessed from public repositories, we demonstrated that these variants reside within a regulatory region that we determine to be critical for cell-specific silencing. Importantly, this has revealed a disease mechanism for non-coding variants that cause inappropriate expression of a disallowed gene.

Quantitative analysis of transcriptome dynamics provides novel insights into developmental state transitions

BACKGROUND

During embryogenesis, the developmental potential of initially pluripotent cells becomes progressively restricted as they transit to lineage restricted states. The pluripotent cells of Xenopus blastula-stage embryos are an ideal system in which to study cell state transitions during developmental decision-making, as gene expression dynamics can be followed at high temporal resolution.

RESULTS

Here we use transcriptomics to interrogate the process by which pluripotent cells transit to four different lineage-restricted states: neural progenitors, epidermis, endoderm and ventral mesoderm, providing quantitative insights into the dynamics of Waddington's landscape. Our findings provide novel insights into why the neural progenitor state is the default lineage state for pluripotent cells and uncover novel components of lineage-specific gene regulation. These data reveal an unexpected overlap in the transcriptional responses to BMP4/7 and Activin signaling and provide mechanistic insight into how the timing of signaling inputs such as BMP are temporally controlled to ensure correct lineage decisions.

CONCLUSIONS

Together these analyses provide quantitative insights into the logic and dynamics of developmental decision making in early embryos. They also provide valuable lineage-specific time series data following the acquisition of specific lineage states during development.

Differential nuclear import sets the timing of protein access to the embryonic genome

The development of a fertilized egg to an embryo requires the proper temporal control of gene expression. During cell differentiation, timing is often controlled via cascades of transcription factors (TFs). However, in early development, transcription is often inactive, and many TF levels stay constant, suggesting that alternative mechanisms govern the observed rapid and ordered onset of gene expression. Here, we find that in early embryonic development access of maternally deposited nuclear proteins to the genome is temporally ordered via importin affinities, thereby timing the expression of downstream targets. We quantify changes in the nuclear proteome during early development and find that nuclear proteins, such as TFs and RNA polymerases, enter the nucleus sequentially. Moreover, we find that the timing of nuclear proteins' access to the genome corresponds to the timing of downstream gene activation. We show that the affinity of proteins to importin is a major determinant in the timing of protein entry into embryonic nuclei. Thus, we propose a mechanism by which embryos encode the timing of gene expression in early development via biochemical affinities. This process could be critical for embryos to organize themselves before deploying the regulatory cascades that control cell identities.

Dynamic Transcriptional Landscape of Grass Carp (Ctenopharyngodon idella) Reveals Key Transcriptional Features Involved in Fish Development

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

The accuracy of absolute differential abundance analysis from relative count data

Concerns have been raised about the use of relative abundance data derived from next generation sequencing as a proxy for absolute abundances. For example, in the differential abundance setting, compositional effects in relative abundance data may give rise to spurious differences (false positives) when considered from the absolute perspective. In practice however, relative abundances are often transformed by renormalization strategies intended to compensate for these effects and the scope of the practical problem remains unclear. We used simulated data to explore the consistency of differential abundance calling on renormalized relative abundances versus absolute abundances and find that, while overall consistency is high, with a median sensitivity (true positive rates) of 0.91 and specificity (1—false positive rates) of 0.89, consistency can be much lower where there is widespread change in the abundance of features across conditions. We confirm these findings on a large number of real data sets drawn from 16S metabarcoding, expression array, bulk RNA-seq, and single-cell RNA-seq experiments, where data sets with the greatest change between experimental conditions are also those with the highest false positive rates. Finally, we evaluate the predictive utility of summary features of relative abundance data themselves. Estimates of sparsity and the prevalence of feature-level change in relative abundance data give reasonable predictions of discrepancy in differential abundance calling in simulated data and can provide useful bounds for worst-case outcomes in real data.

Molecular sequence counting is a near-ubituiqous method for taking “snapshots” of the state of biological systems at the molecular level and is applied to problems as diverse as profiling gene expression and characterizing bacterial community composition. However, concerns exist about the interpretation of these data, given they are relative counts. In particular some feature-level differences between samples may be technical, not biological, stemming from compositional effects. Here, we quantify the accuracy of estimates of sample-sample differences made from relative versus “absolute” molecular count data, using a comprehensive simulation strategy and published experimental data. We find the accuracy of difference estimation is high in at least 50% of simulated and real data sets but that low accuracy outcomes are far from rare. Further, we observe similar numbers of these low accuracy cases when using any of several popular methods for estimating differences in biological count data. Our results support the use of complementary reference measures of absolute abundance (like RNA spike-ins) for normalizing next-generation sequencing data. We briefly validate the use of these reference quantities and of stringent effect size thresholds as strategies for mitigating interpretational problems with relative count data.

A Comprehensive Analysis of Fibrillar Collagens in Lamprey Suggests a Conserved Role in Vertebrate Musculoskeletal Evolution

Vertebrates have distinct tissues which are not present in invertebrate chordates nor other metazoans. The rise of these tissues also coincided with at least one round of whole-genome duplication as well as a suite of lineage-specific segmental duplications. Understanding whether novel genes lead to the origin and diversification of novel cell types, therefore, is of great importance in vertebrate evolution. Here we were particularly interested in the evolution of the vertebrate musculoskeletal system, the muscles and connective tissues that support a diversity of body plans. A major component of the musculoskeletal extracellular matrix (ECM) is fibrillar collagens, a gene family which has been greatly expanded upon in vertebrates. We thus asked whether the repertoire of fibrillar collagens in vertebrates reflects differences in the musculoskeletal system. To test this, we explored the diversity of fibrillar collagens in lamprey, a jawless vertebrate which diverged from jawed vertebrates (gnathostomes) more than five hundred million years ago and has undergone its own gene duplications. Some of the principal components of vertebrate hyaline cartilage are the fibrillar collagens type II and XI, but their presence in cartilage development across all vertebrate taxa has been disputed. We particularly emphasized the characterization of genes in the lamprey hyaline cartilage, testing if its collagen repertoire was similar to that in gnathostomes. Overall, we discovered thirteen fibrillar collagens from all known gene subfamilies in lamprey and were able to identify several lineage-specific duplications. We found that, while the collagen loci have undergone rearrangement, the Clade A genes have remained linked with the hox clusters, a phenomenon also seen in gnathostomes. While the lamprey muscular tissue was largely similar to that seen in gnathostomes, we saw considerable differences in the larval lamprey skeletal tissue, with distinct collagen combinations pertaining to different cartilage types. Our gene expression analyses were unable to identify type II collagen in the sea lamprey hyaline cartilage nor any other fibrillar collagen during chondrogenesis at the stages observed, meaning that sea lamprey likely no longer require these genes during early cartilage development. Our findings suggest that fibrillar collagens were multifunctional across the musculoskeletal system in the last common ancestor of vertebrates and have been largely conserved, but these genes alone cannot explain the origin of novel cell types.

Modeling Human Genetic Disorders with CRISPR Technologies in Xenopus

Combining the power of Xenopus developmental biology with CRISPR-based technologies promises great discoveries in understanding and treating human genetic disorders. Here we provide a practical pipeline for how to go from known disease gene(s) or risk gene(s) of interest to methods for gaining functional insight into the contribution of these genes to disorder etiology in humans.

Uncovering the mesendoderm gene regulatory network through multi-omic data integration

Mesendodermal specification is one of the earliest events in embryogenesis, where cells first acquire distinct identities. Cell differentiation is a highly regulated process that involves the function of numerous transcription factors (TFs) and signaling molecules, which can be described with gene regulatory networks (GRNs). Cell differentiation GRNs are difficult to build because existing mechanistic methods are low throughput, and high-throughput methods tend to be non-mechanistic. Additionally, integrating highly dimensional data composed of more than two data types is challenging. Here, we use linked self-organizing maps to combine chromatin immunoprecipitation sequencing (ChIP-seq)/ATAC-seq with temporal, spatial, and perturbation RNA sequencing (RNA-seq) data from Xenopus tropicalis mesendoderm development to build a high-resolution genome scale mechanistic GRN. We recover both known and previously unsuspected TF-DNA/TF-TF interactions validated through reporter assays. Our analysis provides insights into transcriptional regulation of early cell fate decisions and provides a general approach to building GRNs using highly dimensional multi-omic datasets.

Temporal and spatial transcriptomic dynamics across brain development in Xenopus laevis tadpoles

Amphibian metamorphosis is a transitional period that involves significant changes in the cell-type populations and biological processes occurring in the brain. Analysis of gene expression dynamics during this process may provide insight into the molecular events underlying these changes. We conducted differential gene expression analyses of the developing Xenopus laevis tadpole brain during this period in two ways: first, over stages of the development in the midbrain and, second, across regions of the brain at a single developmental stage. We found that genes pertaining to positive regulation of neural progenitor cell proliferation as well as known progenitor cell markers were upregulated in the midbrain prior to metamorphic climax; concurrently, expression of cell cycle timing regulators decreased across this period, supporting the notion that cell cycle lengthening contributes to a decrease in proliferation by the end of metamorphosis. We also found that at the start of metamorphosis, neural progenitor populations appeared to be similar across the fore-, mid-, and hindbrain regions. Genes pertaining to negative regulation of differentiation were upregulated in the spinal cord compared to the rest of the brain, however, suggesting that different programs may regulate neurogenesis there. Finally, we found that regulation of biological processes like cell fate commitment and synaptic signaling follow similar trajectories in the brain across early tadpole metamorphosis and mid- to late-embryonic mouse development. By comparing expression across both temporal and spatial conditions, we have been able to illuminate cell-type and biological pathway dynamics in the brain during metamorphosis.

Nanopore Whole Transcriptome Analysis and Pathogen Surveillance by a Novel Solid-Phase Catalysis Approach

The requirement of a large input amount (500 ng) for Nanopore direct RNA-seq presents a major challenge for low input transcriptomic analysis and early pathogen surveillance. The high RNA input requirement is attributed to significant sample loss associated with library preparation using solid-phase reversible immobilization (SPRI) beads. A novel solid-phase catalysis strategy for RNA library preparation to circumvent the need for SPRI bead purification to remove enzymes is reported here. This new approach leverages concurrent processing of non-polyadenylated transcripts with immobilized poly(A) polymerase and T4 DNA ligase, followed by directly loading the prepared library onto a flow cell. Whole transcriptome sequencing, using a human pathogen Listeria monocytogenes as a model, demonstrates this new method displays little sample loss, takes much less time, and generates higher sequencing throughput correlated with reduced nanopore fouling compared to the current library preparation for 500 ng input. Consequently, this approach enables Nanopore low-input direct RNA-seq, improving pathogen detection and transcript identification in a microbial community standard with spike-in transcript controls. Besides, as evident in the bioinformatic analysis, the new method provides accurate RNA consensus with high fidelity and identifies higher numbers of expressed genes for both high and low input RNA amounts.

T-Cell Factors as Transcriptional Inhibitors: Activities and Regulations in Vertebrate Head Development

Since its first discovery in the late 90s, Wnt canonical signaling has been demonstrated to affect a large variety of neural developmental processes, including, but not limited to, embryonic axis formation, neural proliferation, fate determination, and maintenance of neural stem cells. For decades, studies have focused on the mechanisms controlling the activity of β-catenin, the sole mediator of Wnt transcriptional response. More recently, the spotlight of research is directed towards the last cascade component, the T-cell factor (TCF)/Lymphoid-Enhancer binding Factor (LEF), and more specifically, the TCF/LEF-mediated switch from transcriptional activation to repression, which in both embryonic blastomeres and mouse embryonic stem cells pushes the balance from pluri/multipotency towards differentiation. It has been long known that Groucho/Transducin-Like Enhancer of split (Gro/TLE) is the main co-repressor partner of TCF/LEF. More recently, other TCF/LEF-interacting partners have been identified, including the pro-neural BarH-Like 2 (BARHL2), which belongs to the evolutionary highly conserved family of homeodomain-containing transcription factors. This review describes the activities and regulatory modes of TCF/LEF as transcriptional repressors, with a specific focus on the functions of Barhl2 in vertebrate brain development. Specific attention is given to the transcriptional events leading to formation of the Organizer, as well as the roles and regulations of Wnt/β-catenin pathway in growth of the caudal forebrain. We present TCF/LEF activities in both embryonic and neural stem cells and discuss how alterations of this pathway could lead to tumors.

A convergent molecular network underlying autism and congenital heart disease

Patients with neurodevelopmental disorders, including autism, have an elevated incidence of congenital heart disease, but the extent to which these conditions share molecular mechanisms remains unknown. Here, we use network genetics to identify a convergent molecular network underlying autism and congenital heart disease. This network is impacted by damaging genetic variants from both disorders in multiple independent cohorts of patients, pinpointing 101 genes with shared genetic risk. Network analysis also implicates risk genes for each disorder separately, including 27 previously unidentified genes for autism and 46 for congenital heart disease. For 7 genes with shared risk, we create engineered disruptions in Xenopus tropicalis, confirming both heart and brain developmental abnormalities. The network includes a family of ion channels, such as the sodium transporter SCN2A, linking these functions to early heart and brain development. This study provides a road map for identifying risk genes and pathways involved in co-morbid conditions.

cdev: a ground-truth based measure to evaluate RNA-seq normalization performance

Normalization of RNA-seq data has been an active area of research since the problem was first recognized a decade ago. Despite the active development of new normalizers, their performance measures have been given little attention. To evaluate normalizers, researchers have been relying on ad hoc measures, most of which are either qualitative, potentially biased, or easily confounded by parametric choices of downstream analysis. We propose a metric called condition-number based deviation, or cdev, to quantify normalization success. cdev measures how much an expression matrix differs from another. If a ground truth normalization is given, cdev can then be used to evaluate the performance of normalizers. To establish experimental ground truth, we compiled an extensive set of public RNA-seq assays with external spike-ins. This data collection, together with cdev, provides a valuable toolset for benchmarking new and existing normalization methods.

Analysis of the early response to spinal cord injury identified a key role for mTORC1 signaling in the activation of neural stem progenitor cells

Xenopus laevis are able to regenerate the spinal cord during larvae stages through the activation of neural stem progenitor cells (NSPCs). Here we use high-resolution expression profiling to characterize the early transcriptome changes induced after spinal cord injury, aiming to identify the signals that trigger NSPC proliferation. The analysis delineates a pathway that starts with a rapid and transitory activation of immediate early genes, followed by migration processes and immune response genes, the pervasive increase of NSPC-specific ribosome biogenesis factors, and genes involved in stem cell proliferation. Western blot and immunofluorescence analysis showed that mTORC1 is rapidly and transiently activated after SCI, and its pharmacological inhibition impairs spinal cord regeneration and proliferation of NSPC through the downregulation of genes involved in the G1/S transition of cell cycle, with a strong effect on PCNA. We propose that the mTOR signaling pathway is a key player in the activation of NPSCs during the early steps of spinal cord regeneration.

Galloway-Mowat syndrome: New insights from bioinformatics and expression during Xenopus embryogenesis

Galloway-Mowat syndrome (GAMOS) is a rare developmental disease. Patients suffer from congenital brain anomalies combined with renal abnormalities often resulting in an early-onset steroid-resistant nephrotic syndrome. The etiology of GAMOS has a heterogeneous genetic contribution. Mutations in more than 10 different genes have been reported in GAMOS patients. Among these are mutations in four genes encoding members of the human KEOPS (kinase, endopeptidase and other proteins of small size) complex, including OSGEP, TP53RK, TPRKB and LAGE3. Until now, these components have been functionally mainly investigated in bacteria, eukarya and archaea and in humans in the context of the discovery of its role in GAMOS, but the KEOPS complex members' expression and function during embryogenesis in vertebrates is still unknown. In this study, in silico analysis showed that both gene localization and the protein sequences of the three core KEOPS complex members Osgep, Tp53rk and Tprkb are highly conserved across different species including Xenopus laevis. In addition, we examined the spatio-temporal expression pattern of osgep, tp53rk and tprkb using RT-PCR and whole mount in situ hybridization approaches during early Xenopus development. We observed that all three genes were expressed during early embryogenesis and enriched in tissues and organs affected in GAMOS. More precisely, KEOPS complex genes are expressed in the pronephros, but also in neural tissue such as the developing brain, eye and cranial cartilage. These findings suggest that the KEOPS complex plays an important role during vertebrate embryonic development.

The DNA-to-cytoplasm ratio broadly activates zygotic gene expression in Xenopus

In multicellular animals, the first major event after fertilization is the switch from maternal to zygotic control of development. During this transition, zygotic gene transcription is broadly activated in an otherwise quiescent genome in a process known as zygotic genome activation (ZGA). In fast-developing embryos, ZGA often overlaps with the slowing of initially synchronous cell divisions at the mid-blastula transition (MBT). Initial studies of the MBT led to the nuclear-to-cytoplasmic ratio model where MBT timing is regulated by the exponentially increasing amounts of some nuclear component "N" titrated against a fixed cytoplasmic component "C." However, more recent experiments have been interpreted to suggest that ZGA is independent of the N/C ratio. To determine the role of the N/C ratio in ZGA, we generated Xenopus frog embryos with ∼3-fold differences in genomic DNA (i.e., N) by using X. tropicalis sperm to fertilize X. laevis eggs with or without their maternal genome. Resulting embryos have otherwise identical X. tropicalis genome template amounts, embryo sizes, and X. laevis maternal environments. We generated transcriptomic time series across the MBT in both conditions and used X. tropicalis paternally derived mRNA to identify a high-confidence set of exclusively zygotic transcripts. Both ZGA and the increase in cell-cycle duration are delayed in embryos with ∼3-fold less DNA per cell. Thus, DNA is an important component of the N/C ratio, which is a critical regulator of zygotic genome activation in Xenopus embryos.

CRISPR-Cas9 Mutagenesis in Xenopus tropicalis for Phenotypic Analyses in the F0 Generation and Beyond

CRISPR-Cas9 mutagenesis is being widely used to create targeted loss-of-function mutations in the diploid frog Xenopus tropicalis Here we describe a simple mutagenesis protocol using microinjection of Cas9 protein or mRNA, together with synthetic guide RNAs (sgRNAs) targeting specific DNA sequences, into the early embryo. Cas9-catalyzed double-strand breaks undergo error-prone repair, resulting in production of short insertions and/or deletions. Thus, careful selection of target sites can lead to mutations that impair normal function of the protein product. CRISPR-Cas9 can be used to create either mosaic loss-of-function Xenopus embryos that display F₀ generation phenotypes or mutant lines for downstream analysis. In addition to describing how to mutagenize genes using CRISPR-Cas9, we also discuss a simple method to determine the mutagenesis efficiency, some potential problems that can arise, and possible solutions to overcome them. The protocol described here should be applicable to other amphibians and, in principle, many other organisms.

The tweety Gene Family: From Embryo to Disease

The tweety genes encode gated chloride channels that are found in animals, plants, and even simple eukaryotes, signifying their deep evolutionary origin. In vertebrates, the tweety gene family is highly conserved and consists of three members—ttyh1, ttyh2, and ttyh3—that are important for the regulation of cell volume. While research has elucidated potential physiological functions of ttyh1 in neural stem cell maintenance, proliferation, and filopodia formation during neural development, the roles of ttyh2 and ttyh3 are less characterized, though their expression patterns during embryonic and fetal development suggest potential roles in the development of a wide range of tissues including a role in the immune system in response to pathogen-associated molecules. Additionally, members of the tweety gene family have been implicated in various pathologies including cancers, particularly pediatric brain tumors, and neurodegenerative diseases such as Alzheimer’s and Parkinson’s disease. Here, we review the current state of research using information from published articles and open-source databases on the tweety gene family with regard to its structure, evolution, expression during development and adulthood, biochemical and cellular functions, and role in human disease. We also identify promising areas for further research to advance our understanding of this important, yet still understudied, family of genes.

Bayesian Modeling of Spatial Molecular Profiling Data via Gaussian Process

MOTIVATION

The location, timing, and abundance of gene expression (both mRNA and proteins) within a tissue define the molecular mechanisms of cell functions. Recent technology breakthroughs in spatial molecular profiling, including imaging-based technologies and sequencing-based technologies, have enabled the comprehensive molecular characterization of single cells while preserving their spatial and morphological contexts. This new bioinformatics scenario calls for effective and robust computational methods to identify genes with spatial patterns.

RESULTS

We represent a novel Bayesian hierarchical model to analyze spatial transcriptomics data, with several unique characteristics. It models the zero-inflated and over-dispersed counts by deploying a zero-inflated negative binomial model that greatly increases model stability and robustness. Besides, the Bayesian inference framework allows us to borrow strength in parameter estimation in a de novo fashion. As a result, the proposed model shows competitive performances in accuracy and robustness over existing methods in both simulation studies and two real data applications.

AVAILABILITY

The related R/C ++ source code is available at https://github.com/Minzhe/BOOST-GP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Anaplastic lymphoma kinase (alk), a neuroblastoma associated gene, is expressed in neural crest domains during embryonic development of Xenopus

Neuroblastoma is a neural crest-derived paediatric cancer that is the most common and deadly solid extracranial tumour of childhood. It arises when neural crest cells fail to follow their differentiation program to give rise to cells of the sympathoadrenal lineage. These undifferentiated cells can proliferate and migrate, forming tumours mostly found associated with the adrenal glands. Activating mutations in the kinase domain of anaplastic lymphoma kinase (ALK) are linked to high-risk cases, where extensive therapy is ineffective. However, the role of ALK in embryonic development, downstream signal transduction and in metastatic transformation of the neural crest is poorly understood. Here, we demonstrate high conservation of the ALK protein sequences among vertebrates. We then examine alk mRNA expression in the frog models Xenopus laevis and Xenopus tropicalis. Using in situ hybridisation of Xenopus embryos, we show that alk is expressed in neural crest domains throughout development, suggesting a possible role in neuroblastoma initiation. Lastly, RT-qPCR analyses show high levels of alk expression at tadpole stages. Collectively, these data may begin to elucidate how alk functions in neural crest cells and how its deregulation can result in tumorigenesis.