Cellular remodeling and JAK inhibition promote zygotic gene expression in the Ciona germline.

Transcription control is a major determinant of cell fate decisions in somatic tissues. By contrast, early germline fate specification in numerous vertebrate and invertebrate species relies extensively on RNA-level regulation, exerted on asymmetrically inherited maternal supplies, with little-to-no zygotic transcription. However delayed, a maternal-to-zygotic transition is nevertheless poised to complete the deployment of pre-gametic programs in the germline. Here, we focus on early germline specification in the tunicate Ciona to study zygotic genome activation. We first demonstrate that a peculiar cellular remodeling event excludes localized postplasmic Pem-1 mRNA, which encodes the general inhibitor of transcription. Subsequently, zygotic transcription begins in Pem-1-negative primordial germ cells (PGCs), as revealed by histochemical detection of elongating RNA Polymerase II, and nascent Mef2 transcripts. In addition, we uncover a provisional antagonism between JAK and MEK/BMPRI/GSK3 signaling, which controls the onset of zygotic gene expression, following cellular remodeling of PGCs. We propose a 2-step model for the onset of zygotic transcription in the Ciona germline and discuss the significance of germ plasm dislocation and remodeling in the context of developmental fate specification.

Optimizing CRISPR/Cas9 approaches in the polymorphic tunicate Ciona intestinalis.

CRISPR/Cas9 became a powerful tool for genetic engineering and in vivo knockout also in the invertebrate chordate Ciona intestinalis. Ciona (ascidians, tunicates) is an important model organism because it shares developmental features with the vertebrates, considered the sister group of tunicates, and offers outstanding experimental advantages: a compact genome and an invariant developmental cell lineage that, combined with electroporation mediated transgenesis allows for precise and cell type specific targeting in vivo. A high polymorphism and the mosaic expression of electroporated constructs, however, often hamper the efficient CRISPR knockout, and an optimization in Ciona is desirable. Furthermore, seasonality and artificial maintenance settings can profit from in vitro approaches that would save on animals. Here we present improvements for the CRISPR/Cas9 protocol in silico, in vitro and in vivo. Firstly, in designing sgRNAs, prior sequencing of target genomic regions from experimental animals and alignment with reference genomes of C. robusta and C. intestinalis render a correction possible of subspecies polymorphisms. Ideally, the screening for efficient and non-polymorphic sgRNAs will generate a database compatible for worldwide Ciona populations. Secondly, we challenged in vitro assays for sgRNA validation towards reduced in vivo experimentation and report their suitability but also overefficiency concerning mismatch tolerance. Thirdly, when comparing Cas9 with Cas9:Geminin, thought to synchronize editing and homology-direct repair, we could indeed increase the in vivo efficiency and notably the access to an early expressed gene. Finally, for in vivo CRISPR, genotyping by next generation sequencing (NGS) ex vivo streamlined the definition of efficient single guides. Double CRISPR then generates large deletions and reliable phenotypic excision effects. Overall, while these improvements render CRISPR more efficient in Ciona, they are useful when newly establishing the technique and very transferable to CRISPR in other organisms.

Brachyury controls Ciona notochord fate as part of a feed-forward network.

The notochord is a defining feature of the chordates. The transcription factor Brachyury (Bra) is a key regulator of notochord fate but here we show that it is not a unitary master regulator in the model chordate Ciona Ectopic Bra expression only partially reprograms other cell types to a notochord-like transcriptional profile and a subset of notochord-enriched genes is unaffected by CRISPR Bra disruption. We identify Foxa.a and Mnx as potential co-regulators, and find that combinatorial cocktails are more effective at reprogramming other cell types than Bra alone. We reassess the network relationships between Bra, Foxa.a and other components of the notochord gene regulatory network, and find that Foxa.a expression in the notochord is regulated by vegetal FGF signaling. It is a direct activator of Bra expression and has a binding motif that is significantly enriched in the regulatory regions of notochord-enriched genes. These and other results indicate that Bra and Foxa.a act together in a regulatory network dominated by positive feed-forward interactions, with neither being a classically defined master regulator.

The complete cell lineage and MAPK- and Otx-dependent specification of the dopaminergic cells in the Ciona brain.

The Ciona larva has served as a unique model for understanding the development of dopaminergic cells at single-cell resolution owing to the exceptionally small number of neurons in its brain and its fixed cell lineage during embryogenesis. A recent study suggested that the transcription factors Fer2 and Meis directly regulate the dopamine synthesis genes in Ciona, but the dopaminergic cell lineage and the gene regulatory networks that control the development of dopaminergic cells have not been fully elucidated. Here, we reveal that the dopaminergic cells in Ciona are derived from a bilateral pair of cells called a9.37 cells at the center of the neural plate. The a9.37 cells divide along the anterior-posterior axis, and all of the descendants of the posterior daughter cells differentiate into the dopaminergic cells. We show that the MAPK pathway and the transcription factor Otx are required for the expression of Fer2 in the dopaminergic cell lineage. Our findings establish the cellular and molecular framework for fully understanding the commitment to dopaminergic cells in the simple chordate brain.

Differentially expressed chaperone genes reveal a stress response required for unidirectional regeneration in the basal chordate Ciona.

BACKGROUND: Unidirectional regeneration in the basal chordate Ciona intestinalis involves the proliferation of adult stem cells residing in the branchial sac vasculature and the migration of progenitor cells to the site of distal injury. However, after the Ciona body is bisected, regeneration occurs in the proximal but not in the distal fragments, even if the latter include a part of the branchial sac with stem cells. A transcriptome was sequenced and assembled from the isolated branchial sacs of regenerating animals, and the information was used to provide insights into the absence of regeneration in distal body fragments. RESULTS: We identified 1149 differentially expressed genes, which were separated into two major modules by weighted gene correlation network analysis, one consisting of mostly upregulated genes correlated with regeneration and the other consisting of only downregulated genes associated with metabolism and homeostatic processes. The hsp70, dnaJb4, and bag3 genes were among the highest upregulated genes and were predicted to interact in an HSP70 chaperone system. The upregulation of HSP70 chaperone genes was verified and their expression confirmed in BS vasculature cells previously identified as stem and progenitor cells. siRNA-mediated gene knockdown showed that hsp70 and dnaJb4, but not bag3, are required for progenitor cell targeting and distal regeneration. However, neither hsp70 nor dnaJb4 were strongly expressed in the branchial sac vasculature of distal fragments, implying the absence of a stress response. Heat shock treatment of distal body fragments activated hsp70 and dnaJb4 expression indicative of a stress response, induced cell proliferation in branchial sac vasculature cells, and promoted distal regeneration. CONCLUSIONS: The chaperone system genes hsp70, dnaJb4, and bag3 are significantly upregulated in the branchial sac vasculature following distal injury, defining a stress response that is essential for regeneration. The stress response is absent from distal fragments, but can be induced by a heat shock, which activates cell division in the branchial sac vasculature and promotes distal regeneration. This study demonstrates the importance of a stress response for stem cell activation and regeneration in a basal chordate, which may have implications for understanding the limited regenerative activities in other animals, including vertebrates.

Proteomic identification of intracellular vesicle trafficking and protein glycosylation requirements for lumen inflation in Ciona notochord.

Lumen formation and inflation are crucial steps for tubular organ morphogenesis, yet the underling mechanism remains largely unrevealed. Here, we applied 4D proteomics to screen the lumenogenesis-related proteins and revealed the biological pathways potentially that are involved in lumen inflation during notochord lumen formation in the ascidian Ciona savignyi. In total, 910 differentiated expressed proteins (DEPs) were identified before and after notochord lumen formation utilizing Mfuzz analysis. Those DEPs were grouped into four upregulated clusters based on their quantitative expression patterns; the functions of these proteins were enriched in protein metabolic and biosynthetic process, the establishment of localization, and vesicle-mediated transport. We analyzed the vesicle trafficking cluster and focused on several vesicle transport hub proteins. In vivo function-deficient experiments showed that mutation of vesicle transport proteins resulted in an abnormal lumen in notochord development, demonstrating the crucial role of intracellular trafficking for lumen formation. Moreover, abundant extracellular matrix proteins were identified, the majority of which were predicted to be glycosylated proteins. Inhibition of glycosylation markedly reduced the lumen expansion rate in notochord cells, suggesting that protein glycosylation is essential for lumenogenesis. Overall, our study provides an invaluable resource and reveals the crucial mechanisms in lumen formation and expansion.

Ciona embryonic tail bending is driven by asymmetrical notochord contractility and coordinated by epithelial proliferation.

Ventral bending of the embryonic tail within the chorion is an evolutionarily conserved morphogenetic event in both invertebrates and vertebrates. However, the complexity of the anatomical structure of vertebrate embryos makes it difficult to experimentally identify the mechanisms underlying embryonic folding. This study investigated the mechanisms underlying embryonic tail bending in chordates. To further understand the mechanical role of each tissue, we also developed a physical model with experimentally measured parameters to simulate embryonic tail bending. Actomyosin asymmetrically accumulated at the ventral side of the notochord, and cell proliferation of the dorsal tail epidermis was faster than that in the ventral counterpart during embryonic tail bending. Genetic disruption of actomyosin activity and inhibition of cell proliferation dorsally caused abnormal tail bending, indicating that both asymmetrical actomyosin contractility in the notochord and the discrepancy of epidermis cell proliferation are required for tail bending. In addition, asymmetrical notochord contractility was sufficient to drive embryonic tail bending, whereas differential epidermis proliferation was a passive response to mechanical forces. These findings showed that asymmetrical notochord contractility coordinates with differential epidermis proliferation mechanisms to drive embryonic tail bending.This article has an associated 'The people behind the papers' interview.

Different strategies for tissue scaling in dwarf tailbud embryos revealed by single-cell analysis.

The tailbud stage is part of the organogenesis period-an evolutionarily conserved developmental period among chordates that is essential for determining the characteristics of the chordate body plan. When the volume of the egg is artificially decreased by cutting, ascidians produce a normal-looking but miniature (dwarf) tailbud embryo. Although cell lineages during ascidian embryogenesis are invariant, the number of cell divisions in the dwarf embryo is altered by a different mechanism in each tissue (Yamada and Nishida, 1999). Here, to elucidate the size-regulation strategies of the Ciona robusta dwarf tailbud embryo, we compared anatomical structure, developmental speed, and cell number/volume in each tissue between dwarf and wild type (WT) embryos. To do this, we constructed a 3D virtual mid-tailbud embryo (Nakamura et al., 2012). We could make a Ciona dwarf tailbud embryo from eggs with a diameter over 108 â€‹µm (correspond to â€‹> â€‹40% of the wild type egg volume). The timings of cleavage (~St. 12) and subsequent morphogenesis were nearly the same but blastomeres of animal hemisphere slightly delayed the timing of mitosis in the early cleavage period. Intriguingly, the tissue-to-tissue volume ratios of dwarf tailbud embryos were similar to those of wild type embryos suggesting that the ratio of tissue volumes is essential for maintaining the proper shape of the tailbud embryo. The number of cells in the epidermis, nervous system, and mesenchyme was significantly reduced in the dwarf embryos whereas the cell volume distribution of these tissues was similar in the dwarf and wild type. In contrast, the number of cells in the notochord, muscle, heart, and endoderm were maintained in the dwarf embryos; cell volumes were significantly reduced. Neither parameter changed in germline precursors. These results indicate that each tissue uses different scaling strategies to coordinate cell number and cell volume in accordance with the embryo size.

Hox13 is essential for formation of a sensory organ at the terminal end of the sperm duct in Ciona.

Species-specific traits are thought to have been acquired by natural selection. Transcription factors play central roles in the evolution of species-specific traits. Hox genes encode a set of conserved transcription factors essential for establishing the anterior-posterior body axis of animals. Changes in the expression or function of Hox genes can lead to the diversification of animal-body plans. The tunicate ascidian Ciona intestinalis Type A has an orange-colored structure at the sperm duct terminus. This orange-pigmented organ (OPO) is the characteristic that can distinguish this ascidian from other closely related species. The OPO is formed by the accumulation of orange-pigmented cells (OPCs) that are present throughout the adult body. We show that Hox13 is essential for formation of the OPO. Hox13 is expressed in the epithelium of the sperm duct and neurons surrounding the terminal openings for sperm ejection, while OPCs themselves do not express this gene. OPCs are mobile cells that can move through the body vasculature by pseudopodia, suggesting that the OPO is formed by the accumulation of OPCs guided by Hox13-positive cells. Another ascidian species, Ciona savignyi, does not have an OPO. Like Hox13 of C. intestinalis, Hox13 of C. savignyi is expressed at the terminus of its sperm duct; however, its expression domain is limited to the circular area around the openings. The genetic changes responsible for the acquisition or loss of OPO are likely to occur in the expression pattern of Hox13.

Highly dynamic transcriptional reprogramming and shorter isoform shifts under acute stresses during biological invasions.

Phenotypic plasticity has been increasingly recognized for its importance in adaptation to novel environments, and initial rapid plastic response to acute stresses usually serves as the stepping stone for future adaptation. Differential gene expression and alternative splicing have been proposed as two underlying mechanisms for rapid plastic response to environmental stresses. Here, we used an invasive model species, Ciona savignyi, to investigate the temporary plastic changes under temperature stresses on gene expression and alternative splicing. Our results revealed rapid and highly dynamic gene expression reprogramming and alternative splicing switch under acute stresses. Distinct transcriptional response profiles were triggered by two types of temperature stresses, showing resilience recovery and increasing divergence under heat and cold challenges, respectively. Interestingly, alternative exons were more inclined to be skipped under both heat and cold stresses, leading to shorter isoforms but with maintained Open Reading Frames (ORFs). Although similar response patterns were observed between differential gene expression and alternative splicing, low overlap between Differentially Expressed Genes (DEGs) and Differentially Alternative Spliced Genes (DASGs) suggests that distinct gene sets and associated functions should be involved in temperature challenges. Thus, alternative splicing should offer an additional layer of plastic response to environmental challenges. Finally, we identified key plastic genes involved in both gene expression regulation and alternative splicing. The results obtained here shed light on adaptation and accommodation mechanisms during biological invasions, particularly for acute environmental changes at early stages of biological invasions such as transport and introduction.

Transcriptional regulation of the Ciona Gsx gene in the neural plate.

The ascidian neural plate consists of a defined number of identifiable cells organized in a grid of rows and columns, representing a useful model to investigate the molecular mechanisms controlling neural patterning in chordates. Distinct anterior brain lineages are specified via unique combinatorial inputs of signalling pathways with Nodal and Delta-Notch signals patterning along the medial-lateral axis and FGF/MEK/ERK signals patterning along the anterior-posterior axis of the neural plate. The Ciona Gsx gene is specifically expressed in the a9.33 cells in the row III/column 2 position of anterior brain lineages, characterised by a combinatorial input of Nodal-OFF, Notch-ON and FGF-ON. Here, we identify the minimal cis-regulatory element (CRE) of 376 bp, which can recapitulate the early activation of Gsx. We show that this minimal CRE responds in the same way as the endogenous Gsx gene to manipulation of FGF- and Notch-signalling pathways and to overexpression of Snail, a mediator of Nodal signals, and Six3/6, which is required to demarcate the anterior boundary of Gsx expression at the late neurula stage. We reveal that sequences proximal to the transcription start site include a temporal regulatory element required for the precise transcriptional onset of gene expression. We conclude that sufficient spatial and temporal information for Gsx expression is integrated in 376 bp of non-coding cis-regulatory sequences.

Evolution of embryonic cis-regulatory landscapes between divergent Phallusia and Ciona ascidians.

Ascidian species of the Phallusia and Ciona genera are distantly related, their last common ancestor dating several hundred million years ago. Although their genome sequences have extensively diverged since this radiation, Phallusia and Ciona species share almost identical early morphogenesis and stereotyped cell lineages. Here, we explored the evolution of transcriptional control between P. mammillata and C. robusta. We combined genome-wide mapping of open chromatin regions in both species with a comparative analysis of the regulatory sequences of a test set of 10 pairs of orthologous early regulatory genes with conserved expression patterns. We find that ascidian chromatin accessibility landscapes obey similar rules as in other metazoa. Open-chromatin regions are short, highly conserved within each genus and cluster around regulatory genes. The dynamics of chromatin accessibility and closest-gene expression are strongly correlated during early embryogenesis. Open-chromatin regions are highly enriched in cis-regulatory elements: 73% of 49 open chromatin regions around our test genes behaved as either distal enhancers or proximal enhancer/promoters following electroporation in Phallusia eggs. Analysis of this datasets suggests a pervasive use in ascidians of "shadow" enhancers with partially overlapping activities. Cross-species electroporations point to a deep conservation of both the trans-regulatory logic between these distantly-related ascidians and the cis-regulatory activities of individual enhancers. Finally, we found that the relative order and approximate distance to the transcription start site of open chromatin regions can be conserved between Ciona and Phallusia species despite extensive sequence divergence, a property that can be used to identify orthologous enhancers, whose regulatory activity can partially diverge.

Positioning a multifunctional basic helix-loop-helix transcription factor within the Ciona notochord gene regulatory network.

In a multitude of organisms, transcription factors of the basic helix-loop-helix (bHLH) family control the expression of genes required for organ development and tissue differentiation. The functions of different bHLH transcription factors in the specification of nervous system and paraxial mesoderm have been widely investigated in various model systems. Conversely, the knowledge of the role of these regulators in the development of the axial mesoderm, the embryonic territory that gives rise to the notochord, and the identities of their target genes, remain still fragmentary. Here we investigated the transcriptional regulation and target genes of Bhlh-tun1, a bHLH transcription factor expressed in the developing Ciona notochord as well as in additional embryonic territories that contribute to the formation of both larval and adult structures. We describe its possible role in notochord formation, its relationship with the key notochord transcription factor Brachyury, and suggest molecular mechanisms through which Bhlh-tun1 controls the spatial and temporal expression of its effectors.

A complete statistical model for calibration of RNA-seq counts using external spike-ins and maximum likelihood theory.

A fundamental assumption, common to the vast majority of high-throughput transcriptome analyses, is that the expression of most genes is unchanged among samples and that total cellular RNA remains constant. As the number of analyzed experimental systems increases however, different independent studies demonstrate that this assumption is often violated. We present a calibration method using RNA spike-ins that allows for the measurement of absolute cellular abundance of RNA molecules. We apply the method to pooled RNA from cell populations of known sizes. For each transcript, we compute a nominal abundance that can be converted to absolute by dividing by a scale factor determined in separate experiments: the yield coefficient of the transcript relative to that of a reference spike-in measured with the same protocol. The method is derived by maximum likelihood theory in the context of a complete statistical model for sequencing counts contributed by cellular RNA and spike-ins. The counts are based on a sample from a fixed number of cells to which a fixed population of spike-in molecules has been added. We illustrate and evaluate the method with applications to two global expression data sets, one from the model eukaryote Saccharomyces cerevisiae, proliferating at different growth rates, and differentiating cardiopharyngeal cell lineages in the chordate Ciona robusta. We tested the method in a technical replicate dilution study, and in a k-fold validation study.

Anti-virus effects of interferon regulatory factors (IRFs) identified in ascidian Ciona savignyi.

Interferon regulatory factors (IRFs) are key transcription factors that function in the immune system via the interferon (IFN) pathway. In the current study, we identified and characterized three IRFs (CsIRFL1, CsIRFL2, and CsIRFL3) from ascidian Ciona savignyi. Phylogenetic analysis showed that CsIRFL1 was clustered with two IRFs from Ciona robusta and shrimp IRF apart from the vertebrate IRFs, whereas CsIRFL2 and CsIRFL3 were grouped with an unnamed protein from Oikopleura dioica into a sub-branch highly identifying with the vertebrate IRF4, IRF8, and IRF9. Gene expression analysis revealed that CsIRFL1 and CsIRFL2 expressed in all the examined adult tissues (stomach, intestines, eggs, hemocytes, gonad, heart, and pharynx) and predominantly in hemocytes. However, the expression of CsIRFL3 was undetectable in the tested adult tissues. Furthermore, in situ hybridization showed that CsIRFL1 and CsIRFL2 mainly expressed in immunocytes within hemolymph, including phagocytes, macrophage-like cells, morula cells, and amoebocytes, suggesting CsIRFL1 and CsIRFL2 were involved in ascidian immune responses. We then performed LPS and poly(I:C) challenge assay and found that CsIRFL1 highly expressed in the cultured hemocytes following LPS infection for 24 h. After viral analogue poly(I:C) stimulation, the expression of CsIRFL2 was dramatically upregulated from 12 to 24 h. Meanwhile, two critical components of the IFN signaling pathways, STAT and TBK1, showed the increased expression as well after poly(I:C) induction, indicating that CsIRFL2 and IFN pathways genes were activated under the infection of viral analogue. Thus, our findings suggested that CsIRFL2 was a potential transcriptional regulatory factor that participated in regulating the ascidian anti-virus immune response.

Efficient genome editing using CRISPR-Cas-mediated homology directed repair in the ascidian Ciona robusta.

Eliminating or silencing a gene's level of activity is one of the classic approaches developmental biologists employ to determine a gene's function. A recently developed method of gene perturbation called CRISPR-Cas, which was derived from a prokaryotic adaptive immune system, has been adapted for use in eukaryotic cells. This technology has been established in several model organisms as a powerful and efficient tool for knocking out or knocking down the function of a gene of interest. It has been recently shown that CRISPR-Cas functions with fidelity and efficiency in Ciona robusta. Here, we show that in C. robusta CRISPR-Cas mediated genomic knock-ins can be efficiently generated. Electroporating a tissue-specific transgene driving Cas9 and a U6-driven gRNA transgene together with a fluorescent protein-containing homology directed repair (FP-HDR) template results in gene-specific patterns of fluorescence consistent with a targeted genomic insertion. Using the Tyrosinase locus to optimize reagents, we first characterize a new Pol III promoter for expressing gRNAs from the Ciona savignyi H1 gene, and then adapt technology that flanks gRNAs by ribozymes allowing cell-specific expression from Pol II promoters. Next, we examine homology arm-length efficiencies of FP-HDR templates. Reagents were then developed for targeting Brachyury and Pou4 that resulted in expected patterns of fluorescence, and sequenced PCR amplicons derived from single embryos validated predicted genomic insertions. Finally, using two differentially colored FP-HDR templates, we show that biallelic FP-HDR template insertion can be detected in live embryos of the F0 generation.

Enhancer activities of amphioxus Brachyury genes in embryos of the ascidian, Ciona intestinalis.

The notochord and somites are distinctive chordate structures. The T-box transcription factor gene, Brachyury, is expressed in notochord and plays a pivotal role in its formation. In the cephalochordate, Branchiostoma floridae, Brachyury is duplicated into BfBra1 and BfBra2, which are expressed in the somite-formation region as well. In a series of experiments to elucidate the regulatory machinery of chordate Brachyury expression, we carried out a lacZ reporter assay of BfBra in embryos of the urochordate, Ciona intestinalis. Vista analyses suggest the presence of conserved non-coding sequences, not only in the 5'-upstream, but also in the 3'-downstream and in introns of BfBra. We found that: (1) 5'-upstream sequences of both BfBra1 and BfBra2 promote lacZ expression in muscle cells, (2) 3'-downstream sequences have enhancer activity that promotes lacZ expression in notochord cells, and (3) introns of BfBra2 and BfBra1 exhibit lacZ expression preferentially in muscle and notochord cells. These results suggest shared cephalochordate Brachyury enhancer machinery that also works in urochordates. We discuss the results in relation to evolutionary modification of Brachyury expression in formation of chordate-specific organs characteristic of each lineage.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona.

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

Co-expression network analysis of environmental canalization in the ascidian Ciona.

BACKGROUND: Canalization, or buffering, is defined as developmental stability in the face of genetic and/or environmental perturbations. Understanding how canalization works is important in predicting how species survive environmental change, as well as deciphering how development can be altered in the evolutionary process. However, how developmental gene expression is linked to buffering remains unclear. We addressed this by co-expression network analysis, comparing gene expression changes caused by heat stress during development at a whole-embryonic scale in reciprocal hybrid crosses of sibling species of the ascidian Ciona that are adapted to different thermal environments. RESULTS: Since our previous work showed that developmental buffering in this group is maternally inherited, we first identified maternal developmental buffering genes (MDBGs) in which the expression level in embryos is both correlated to the level of environmental canalization and also differentially expressed depending on the species' gender roles in hybrid crosses. We found only 15 MDBGs, all of which showed high correlation coefficient values for expression with a large number of other genes, and 14 of these belonged to a single co-expression module. We then calculated correlation coefficients of expression between MDBGs and transcription factors in the central nervous system (CNS) developmental gene network that had previously been identified experimentally. We found that, compared to the correlation coefficients between MDBGs, which had an average of 0.96, the MDBGs are loosely linked to the CNS developmental genes (average correlation coefficient 0.45). Further, we investigated the correlation of each developmental to MDBGs, showing that only four out of 62 CNS developmental genes showed correlation coefficient > 0.9, comparable to the values between MDBGs, and three of these four genes were signaling molecules: BMP2/4, Wnt7, and Delta-like. CONCLUSIONS: We show that the developmental pathway is not centrally located within the buffering network. We found that out of 62 genes in the developmental gene network, only four genes showed correlation coefficients as high as between MDBGs. We propose that loose links to MDBGs stabilize spatiotemporally dynamic development.

Development of a miRNA Sensor by an Inducible CRISPR-Cas9 Construct in Ciona Embryogenesis.

MicroRNAs (miRNAs) regulate multicellular processes and diverse signaling pathways in organisms. The detection of the spatiotemporal expression of miRNA in vivo is crucial for uncovering the function of miRNA. However, most of the current detecting techniques cannot reflect the dynamics of miRNA sensitively in vivo. Here, we constructed a miRNA-induced CRISPR-Cas9 platform (MICR) used in marine chordate Ciona. The key component of MICR is a pre-single guide RNA (sgRNA) flanked by miRNA-binding sites that can be released by RNA-induced silencing complex (RISC) cleavage to form functional sgRNA in the presence of complementary miRNA. By using the miRNA-inducible CRISPR-on system (MICR-ON), we successfully detected the dynamic expression of a miRNA csa-miR-4018a during development of Ciona embryo. The detected patterns were validated to be consistent with the results by in situ hybridization. It is worth noting that the expression of csa-miR-4018a was examined by MICR-ON to be present in additional tissues, where no obvious signaling was detected by in situ hybridization, suggesting that the MICR-ON might be a more sensitive approach to detect miRNA signal in living animal. Thus, MICR-ON was demonstrated to be a sensitive and highly efficient approach for monitoring the dynamics of expression of miRNA in vivo and will facilitate the exploration of miRNA functions in biological systems.