Gene expression patterns that support novel developmental stress buffering in embryos of the annual killifish Austrofundulus limnaeus

Axis formation in annual killifish: Nodal and β-catenin regulate morphogenesis without Huluwa prepatterning

Axis formation in fish and amphibians typically begins with a prepattern of maternal gene products. Annual killifish embryogenesis, however, challenges prepatterning models as blastomeres disperse and then aggregate to form the germ layers and body axes. We show that huluwa, a prepatterning factor thought to break symmetry by stabilizing β-catenin, is truncated and inactive in Nothobranchius furzeri. Nuclear β-catenin is not selectively stabilized on one side of the blastula but accumulates in cells forming the aggregate. Blocking β-catenin activity or Nodal signaling disrupts aggregate formation and germ layer specification. Nodal signaling coordinates cell migration, establishing an early role for this signaling pathway. These results reveal a surprising departure from established mechanisms of axis formation: Huluwa-mediated prepatterning is dispensable, and β-catenin and Nodal regulate morphogenesis.

Origin, form and function of extraembryonic structures in teleost fishes

Teleost eggs have evolved a highly derived early developmental pattern within vertebrates as a result of the meroblastic cleavage pattern, giving rise to a polar stratified architecture containing a large acellular yolk and a small cellular blastoderm on top. Besides the acellular yolk, the teleost-specific yolk syncytial layer (YSL) and the superficial epithelial enveloping layer are recognized as extraembryonic structures that play critical roles throughout embryonic development. They provide enriched microenvironments in which molecular feedback loops, cellular interactions and mechanical signals emerge to sculpt, among other things, embryonic patterning along the dorsoventral and left-right axes, mesendodermal specification and the execution of morphogenetic movements in the early embryo and during organogenesis. An emerging concept points to a critical role of extraembryonic structures in reinforcing early genetic and morphogenetic programmes in reciprocal coordination with the embryonic blastoderm, providing the necessary boundary conditions for development to proceed. In addition, the role of the enveloping cell layer in providing mechanical, osmotic and immunological protection during early stages of development, and the autonomous nutritional support provided by the yolk and YSL, have probably been key aspects that have enabled the massive radiation of teleosts to colonize every ecological niche on the Earth. This article is part of the theme issue 'Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom'.

A pan-metazoan concept for adult stem cells: the wobbling Penrose landscape

Adult stem cells (ASCs) in vertebrates and model invertebrates (e.g. Drosophila melanogaster) are typically long‐lived, lineage‐restricted, clonogenic and quiescent cells with somatic descendants and tissue/organ‐restricted activities. Such ASCs are mostly rare, morphologically undifferentiated, and undergo asymmetric cell division. Characterized by ‘stemness’ gene expression, they can regulate tissue/organ homeostasis, repair and regeneration. By contrast, analysis of other animal phyla shows that ASCs emerge at different life stages, present both differentiated and undifferentiated phenotypes, and may possess amoeboid movement. Usually pluri/totipotent, they may express germ‐cell markers, but often lack germ‐line sequestering, and typically do not reside in discrete niches. ASCs may constitute up to 40% of animal cells, and participate in a range of biological phenomena, from whole‐body regeneration, dormancy, and agametic asexual reproduction, to indeterminate growth. They are considered legitimate units of selection. Conceptualizing this divergence, we present an alternative stemness metaphor to the Waddington landscape: the ‘wobbling Penrose’ landscape. Here, totipotent ASCs adopt ascending/descending courses of an ‘Escherian stairwell’, in a lifelong totipotency pathway. ASCs may also travel along lower stemness echelons to reach fully differentiated states. However, from any starting state, cells can change their stemness status, underscoring their dynamic cellular potencies. Thus, vertebrate ASCs may reflect just one metazoan ASC archetype.

Heat stress-induced alterations in the expression of genes associated with gonadal integrity of the teleost Puntius sophore

Temperature plays an important role on reproductive physiology of vertebrates including mammals, fish, and birds. It has varying effects on fish reproduction depending on the species; higher temperatures favor the spring-spawning species, while lower temperatures stimulate reproduction in autumn spawners. To evaluate the impact of high temperature on the reproductive physiology of minnow Puntius sophore, we carried out expression analysis of selected genes associated with gamete quality (hsp60, hsp70, hsp90, hsf1, vtg), pleuripotency (sox2, oct4, nanog), and sex determination (dmrt1) in gonads (ovary and testis) of P. sophore, heat stressed for different time periods (36 °C/7 days or 60 days) using real-time quantitative polymerase chain reaction (RT-qPCR). Expression of most of the hsp, vtg, and pleuripotency marker genes sox-2, oct-4, and nanog genes was downregulated in both ovary and testis of heat-stressed fish. The expression of dmrt-1 was upregulated in testis but downregulated in ovary of the heat-stressed fish which could be a male favoring effect of high temperature in P. sophore. This study suggests that the reproductive physiology and health of the nutrient dense P. sophore would be negatively affected by high temperature stress.

Small noncoding RNA profiles along alternative developmental trajectories in an annual killifish

Embryonic development of Austrofundulus limnaeus can occur along two phenotypic trajectories that are physiologically and biochemically distinct. Phenotype appears to be influenced by maternal provisioning based on the observation that young females produce predominately non-diapausing embryos and older females produce mostly diapausing embryos. Embryonic incubation temperature can override this pattern and alter trajectory. We hypothesized that temperature-induced phenotypic plasticity may be regulated by post-transcriptional modification via noncoding RNAs. As a first step to exploring this possibility, RNA-seq was used to generate transcriptomic profiles of small noncoding RNAs in embryos developing along the two alternative trajectories. We find distinct profiles of mature sequences belonging to the miR-10 family expressed in increasing abundance during development and mature sequences of miR-430 that follow the opposite pattern. Furthermore, miR-430 sequences are enriched in escape trajectory embryos. MiR-430 family members are known to target maternally provisioned mRNAs in zebrafish and may operate similarly in A. limnaeus in the context of normal development, and also by targeting trajectory-specific mRNAs. This expression pattern and function for miR-430 presents a potentially novel model for maternal-embryonic conflict in gene regulation that provides the embryo the ability to override maternal programming in the face of altered environmental conditions.

The genome of Austrofundulus limnaeus offers insights into extreme vertebrate stress tolerance and embryonic development

Background

The annual killifish Austrofundulus limnaeus inhabits ephemeral ponds in northern Venezuela, South America, and is an emerging extremophile model for vertebrate diapause, stress tolerance, and evolution. Embryos of A. limnaeus regularly experience extended periods of desiccation and anoxia as a part of their natural history and have unique metabolic and developmental adaptations. Currently, there are limited genomic resources available for gene expression and evolutionary studies that can take advantage of A. limnaeus as a unique model system.

Results

We describe the first draft genome sequence of A. limnaeus. The genome was assembled de novo using a merged assembly strategy and was annotated using the NCBI Eukaryotic Annotation Pipeline. We show that the assembled genome has a high degree of completeness in genic regions that is on par with several other teleost genomes. Using RNA-seq and phylogenetic-based approaches, we identify several candidate genes that may be important for embryonic stress tolerance and post-diapause development in A. limnaeus. Several of these genes include heat shock proteins that have unique expression patterns in A. limnaeus embryos and at least one of these may be under positive selection.

Conclusion

The A. limnaeus genome is the first South American annual killifish genome made publicly available. This genome will be a valuable resource for comparative genomics to determine the genetic and evolutionary mechanisms that support the unique biology of annual killifishes. In a broader context, this genome will be a valuable tool for exploring genome-environment interactions and their impacts on vertebrate physiology and evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4539-7) contains supplementary material, which is available to authorized users.

Small noncoding RNA expression during extreme anoxia tolerance of annual killifish (Austrofundulus limnaeus) embryos

Small noncoding RNAs (sncRNA) have recently emerged as specific and rapid regulators of gene expression, involved in a myriad of cellular and organismal processes. MicroRNAs, a class of sncRNAs, are differentially expressed in diverse taxa in response to environmental stress, including anoxia. In most vertebrates, a brief period of oxygen deprivation results in severe tissue damage or death. Studies on sncRNA and anoxia have focused on these anoxia-sensitive species. Studying sncRNAs in anoxia-tolerant organisms may provide insight into adaptive mechanisms supporting anoxia tolerance. Embryos of the annual killifish Austrofundulus limnaeus are the most anoxia-tolerant vertebrates known, surviving over 100 days at their peak tolerance at 25°C. Their anoxia tolerance and physiology vary over development, such that both anoxia-tolerant and anoxia-sensitive phenotypes comprise the species. This allows for a robust comparison to identify sncRNAs essential to anoxia-tolerance. For this study, RNA sequencing was used to identify and quantify expression of sncRNAs in four embryonic stages of A. limnaeus in response to an exposure to anoxia and subsequent aerobic recovery. Unique stage-specific patterns of expression were identified that correlate with anoxia tolerance. In addition, embryos of A. limnaeus appear to constitutively express stress-responsive miRNAs. Most differentially expressed sncRNAs were expressed at higher levels during recovery. Many novel groups of sncRNAs with expression profiles suggesting a key role in anoxia tolerance were identified, including sncRNAs derived from mitochondrial tRNAs. This global analysis has revealed groups of candidate sncRNAs that we hypothesize support anoxia tolerance.

Extra-embryonic tissue spreading directs early embryo morphogenesis in killifish

The spreading of mesenchymal-like cell layers is critical for embryo morphogenesis and tissue repair, yet we know little of this process in vivo. Here we take advantage of unique developmental features of the non-conventional annual killifish embryo to study the principles underlying tissue spreading in a simple cellular environment, devoid of patterning signals and major morphogenetic cell movements. Using in vivo experimentation and physical modelling we reveal that the extra-embryonic epithelial enveloping cell layer, thought mainly to provide protection to the embryo, directs cell migration and the spreading of embryonic tissue during early development. This function relies on the ability of embryonic cells to couple their autonomous random motility to non-autonomous signals arising from the expansion of the extra-embryonic epithelium, mediated by cell membrane adhesion and tension. Thus, we present a mechanism of extra-embryonic control of embryo morphogenesis that couples the mechanical properties of adjacent tissues in the early killifish embryo.

Transcriptomic analysis of maternally provisioned cues for phenotypic plasticity in the annual killifish, Austrofundulus limnaeus

Background

Genotype and environment can interact during development to produce novel adaptive traits that support life in extreme conditions. The development of the annual killifish Austrofundulus limnaeus is unique among vertebrates because the embryos have distinct cell movements that separate epiboly from axis formation during early development, can enter into a state of metabolic dormancy known as diapause and can survive extreme environmental conditions. The ability to enter into diapause can be maternally programmed, with young females producing embryos that do not enter into diapause. Alternately, embryos can be programmed to “escape” from diapause and develop directly by both maternal factors and embryonic incubation conditions. Thus, maternally packaged gene products are hypothesized to regulate developmental trajectory and perhaps the other unique developmental characters in this species.

Results

Using high-throughput RNA sequencing, we generated transcriptomic profiles of mRNAs, long non-coding RNAs and small non-coding RNAs (sncRNAs) in 1–2 cell stage embryos of A. limnaeus. Transcriptomic analyses suggest maternal programming of embryos through alternatively spliced mRNAs and antisense sncRNAs. Comparison of these results to those of comparable studies on zebrafish and other fishes reveals a surprisingly high abundance of transcripts involved in the cellular response to stress and a relatively lower expression of genes required for rapid transition through the cell cycle.

Conclusions

Maternal programming of developmental trajectory is unlikely accomplished by differential expression of diapause-specific genes. Rather, evidence suggests a role for trajectory-specific splice variants of genes expressed in both phenotypes. In addition, based on comparative studies with zebrafish, the A. limnaeus 1–2 cell stage transcriptome is unique in ways that are consistent with their unique life history. These results not only impact our understanding of the genetic mechanisms that regulate entrance into diapause, but also provide insight into the epigenetic regulation of gene expression during development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-017-0069-7) contains supplementary material, which is available to authorized users.

Evolutionary Conservation of pou5f3 Genomic Organization and Its Dynamic Distribution during Embryogenesis and in Adult Gonads in Japanese Flounder Paralichthys olivaceus

Octamer-binding transcription factor 4 (Oct4) is a member of POU (Pit-Oct-Unc) transcription factor family Class V that plays a crucial role in maintaining the pluripotency and self-renewal of stem cells. Though it has been deeply investigated in mammals, its lower vertebrate homologue, especially in the marine fish, is poorly studied. In this study, we isolated the full-length sequence of Paralichthys olivaceus pou5f3 (Popou5f3), and we found that it is homologous to mammalian Oct4. We identified two transcript variants with different lengths of 3′-untranslated regions (UTRs) generated by alternative polyadenylation (APA). Quantitative real-time RT-PCR (qRT-PCR), in situ hybridization (ISH) and immunohistochemistry (IHC) were implemented to characterize the spatial and temporal expression pattern of Popou5f3 during early development and in adult tissues. Our results show that Popou5f3 is maternally inherited, abundantly expressed at the blastula and early gastrula stages, then greatly diminishes at the end of gastrulation. It is hardly detectable from the heart-beating stage onward. We found that Popou5f3 expression is restricted to the adult gonads, and continuously expresses during oogenesis while its dynamics are downregulated during spermatogenesis. Additionally, numerous cis-regulatory elements (CRE) on both sides of the flanking regions show potential roles in regulating the expression of Popou5f3. Taken together, these findings could further our understanding of the functions and evolution of pou5f3 in lower vertebrates, and also provides fundamental information for stem cell tracing and genetic manipulation in Paralichthys olivaceus.

Mitochondrial DNA Sequence and Lack of Response to Anoxia in the Annual Killifish Austrofundulus limnaeus

The annual killifish Austrofundulus limnaeus inhabits ephemeral ponds in regions of Venezuela, South America. Permanent populations of A. limnaeus are maintained by production of stress-tolerant embryos that are able to persist in the desiccated sediment. Previous work has demonstrated that A. limnaeus have a remarkable ability to tolerate extended periods of anoxia and desiccating conditions. After considering temperature, A. limnaeus embryos have the highest known tolerance to anoxia when compared to any other vertebrate yet studied. Oxygen is completely essential for the process of oxidative phosphorylation by mitochondria, the intracellular organelle responsible for the majority of adenosine triphosphate production. Thus, understanding the unique properties of A. limnaeus mitochondria is of great interest. In this work, we describe the first complete mitochondrial genome (mtgenome) sequence of a single adult A. limnaeus individual and compare both coding and non-coding regions to several other closely related fish mtgenomes. Mitochondrial features were predicted using MitoAnnotator and polyadenylation sites were predicted using RNAseq mapping. To estimate the responsiveness of A. limnaeus mitochondria to anoxia treatment, we measure relative mitochondrial DNA copy number and total citrate synthase activity in both relatively anoxia-tolerant and anoxia-sensitive embryonic stages. Our cross-species comparative approach identifies unique features of ND1, ND5, ND6, and ATPase-6 that may facilitate the unique phenotype of A. limnaeus embryos. Additionally, we do not find evidence for mitochondrial degradation or biogenesis during anoxia/reoxygenation treatment in A. limnaeus embryos, suggesting that anoxia-tolerant mitochondria do not respond to anoxia in a manner similar to anoxia-sensitive mitochondria.

Mechanisms of animal diapause: recent developments from nematodes, crustaceans, insects, and fish

Life cycle delays are beneficial for opportunistic species encountering suboptimal environments. Many animals display a programmed arrest of development (diapause) at some stage(s) of their development, and the diapause state may or may not be associated with some degree of metabolic depression. In this review, we will evaluate current advancements in our understanding of the mechanisms responsible for the remarkable phenotype, as well as environmental cues that signal entry and termination of the state. The developmental stage at which diapause occurs dictates and constrains the mechanisms governing diapause. Considerable progress has been made in clarifying proximal mechanisms of metabolic arrest and the signaling pathways like insulin/Foxo that control gene expression patterns. Overlapping themes are also seen in mechanisms that control cell cycle arrest. Evidence is emerging for epigenetic contributions to diapause regulation via small RNAs in nematodes, crustaceans, insects, and fish. Knockdown of circadian clock genes in selected insect species supports the importance of clock genes in the photoperiodic response that cues diapause. A large suite of chaperone-like proteins, expressed during diapause, protects biological structures during long periods of energy-limited stasis. More information is needed to paint a complete picture of how environmental cues are coupled to the signal transduction that initiates the complex diapause phenotype, as well as molecular explanations for how the state is terminated. Excellent examples of molecular memory in post-dauer animals have been documented in Caenorhabditis elegans It is clear that a single suite of mechanisms does not regulate diapause across all species and developmental stages.

Evolution and functions of Oct4 homologs in non-mammalian vertebrates

PouV class transcription factor Oct4/Pou5f1 is a central regulator of indefinite pluripotency in mammalian embryonic stem cells (ESCs) but also participates in cell lineage specification in mouse embryos and in differentiating cell cultures. The molecular basis for this versatility, which is shared between Oct4 and its non-mammalian homologs Pou5f1 and Pou5f3, is not yet completely understood. Here, I review the current understanding of the evolution of PouV class transcription factors and discuss equivalent and diverse roles of Oct4 homologs in pluripotency, differentiation, and cell behavior in different vertebrate embryos. This article is part of a Special Issue entitled: The Oct Transcription Factor Family, edited by Dr. Dean Tantin.

Zygotic Genome Activators, Developmental Timing, and Pluripotency

The transcription factors Pou5f1, Sox2, and Nanog are central regulators of pluripotency in mammalian ES and iPS cells. In vertebrate embryos, Pou5f1/3, SoxB1, and Nanog control zygotic genome activation and participate in lineage decisions. We review the current knowledge of the roles of these genes in developing vertebrate embryos from fish to mammals and suggest a model for pluripotency gene regulatory network functions in early development.

Molecular Cloning, Promoter Analysis and Expression Profiles of the sox3 Gene in Japanese Flounder, Paralichthys olivaceus

Sox3, which belongs to the SoxB1 subgroup, plays major roles in neural and gonadal development. In the present study, Japanese flounder Paralichthys olivaceus sox3 gene (Posox3) and its promoter sequence were isolated and characterized. The deduced PoSox3 protein contained 298 amino acids with a characteristic HMG-box domain. Alignment and phylogenetic analyses indicated that PoSox3 shares highly identical sequence with Sox3 homologues from different species. The promoter region of Posox3 has many potential transcription factor (TF) binding sites. The expression profiles of Posox3 in different developmental stages and diverse adult tissues were analyzed by quantitative real-time RT-PCR (qRT-PCR). Posox3 mRNA was maternally inherited, and maintained at a considerably high expression level between the blastula stage and the hatching stage during embryonic development. Posox3 was abundantly expressed in the adult brain and showed sexually dimorphic expression pattern. In situ hybridization (ISH) was carried out to investigate the cellular distribution of Posox3 in the ovary, and results showed the uniform distribution of Posox3 throughout the cytoplasm of oogonia and stage I–III oocytes. These results indicate that Posox3 has potentially vital roles in embryonic and neural development and may be involved in the oogenesis process. Our work provides a fundamental understanding of the structure and potential functions of Sox3 in Paralichthys olivaceus.