Genome analysis reveals hepatic transcriptional reprogramming changes mediated by enhancers during chick embryonic development

The liver undergoes a slow process for lipid deposition during chick embryonic period. However, the underlying physiological and molecular mechanisms are still unclear. Therefore, the aim of the current study was to reveal the epigenetic mechanism of hepatic transcriptional reprogramming changes based on the integration analysis of RNA-seq and H3K27ac labeled CUT&Tag. Results showed that lipid contents increased gradually with the embryonic age (E) 11, E15, and E19 based on morphological analysis of Hematoxylin-eosin and Oil Red O staining as well as total triglyceride and cholesterol detection. The hepatic protein level of SREBP-1c was higher in E19 when compared with that in E11 and E15, while H3K27ac and H3K4me2 levels declined from E11 to E19. Differential expression genes (DEGs) among these 3 embryonic ages were determined by transcriptome analysis. A total of 107 and 46 genes were gradually upregulated and downregulated respectively with the embryonic age. Meanwhile, differential H3K27ac occupancy in chromatin was investigated. But the integration analysis of RNA-seq and CUT&Tag data showed that the overlap genes were less between DEGs and target genes of differential peaks in the promoter regions. Further, some KEGG pathways enriched from target genes of typical enhancer were overlapped with those from DEGs in transcriptome analysis such as insulin, FoxO, MAPK signaling pathways which were related to lipid metabolism. DNA motif analysis identify 8 and 10 transcription factors (TFs) based on up and down differential peaks individually among E11, E15, and E19 stages where 7 TFs were overlapped including COUP-TFII, FOXM1, FOXA1, HNF4A, RXR, ERRA, FOXA2. These results indicated that H3K27ac histone modification is involved in the transcriptional reprogramming regulation during embryonic development, which could recruit TFs binding to mediate differential enhancer activation. Differential activated enhancer impels dynamic transcriptional reprogramming towards lipid metabolism to promote the occurrence of special phenotype of hepatic lipid deposition.

Induced pluripotent stem cells of endangered avian species

The number of endangered avian-related species increase in Japan recently. The application of new technologies, such as induced pluripotent stem cells (iPSCs), may contribute to the recovery of the decreasing numbers of endangered animals and conservation of genetic resources. We established novel iPSCs from three endangered avian species (Okinawa rail, Japanese ptarmigan, and Blakiston’s fish owl) with seven reprogramming factors (M3O, Sox2, Klf4, c-Myc, Nanog, Lin28, and Klf2). The iPSCs are pluripotency markers and express pluripotency-related genes and differentiated into three germ layers in vivo and in vitro. These three endangered avian iPSCs displayed different cellular characteristics even though the same reprogramming factors use. Japanese ptarmigan-derived iPSCs have different biological characteristics from those observed in other avian-derived iPSCs. Japanese ptarmigan iPSCs contributed to chimeras part in chicken embryos. To the best of our knowledge, our findings provide the first evidence of the potential value of iPSCs as a resource for endangered avian species conservation.

iPSCs from three endangered avian species (including Okinawa rail, Japanese ptarmigan, and Blakiston’s fish owl) are developed and characterized as a potential resource for their conservation.

Single-Cell DNA Methylation Analysis of Chicken Lampbrush Chromosomes

DNA methylation is an essential epigenetic regulation mechanism implicated in transcription and replication control, developmental reprogramming, retroelements silencing and other genomic processes. During mammalian development, a specific DNA methylation pattern should be established in germ cells to allow embryonic development. Less is known about germ cell DNA methylation in other species. To close this gap, we performed a single-cell methylome analysis of chicken diplotene oocytes. We comprehensively characterized methylation patterns in these cells, obtained methylation-based chicken genome segmentation and identified oocyte-specific methylated gene promoters. Our data show that despite the formation of specific transcriptionally hyperactive genome architecture in chicken diplotene oocytes, methylation patterns in these cells closely resemble genomic distribution observed in somatic tissues.

Genome-wide analysis of alternative splicing differences between oocyte and zygote†

Alternative splicing (AS) of mRNA precursors allows the synthesis of multiple mRNAs from a single primary transcript, significantly expanding the information content and regulatory possibilities of higher eukaryotic genomes. During mammalian development, AS drives certain decisive changes in different physiological processes. As development progresses, the maternal-to-zygotic transition (MZT) will trigger two processes: elimination of a subset of maternal mRNA and transcription of the zygote genome begins. Recent high-throughput technological advancements have facilitated genome-wide AS, whereas its analysis in mouse oocyte transition to the zygote stage has not been reported. We present a high-resolution global analysis of AS transitions and discovered extensive AS transitions between mouse oocyte and zygote. The difference of AS patterns was further confirmed using reverse transcription-polymerase chain reaction analysis. Many genes with specific AS events in mouse oocytes are differentially expressed between oocyte and zygote, but only a few genes with specific AS events in zygote are differentially expressed between oocyte and zygote. We provide a landscape of AS events in mouse oocyte and zygote. Our results advance the understanding of AS transitions during mouse fertilization and its potential functions for MZT and further development.

Cyclin D1 gene expression is essential for cell cycle progression from the maternal-to-zygotic transition during blastoderm development in Japanese quail

Embryogenesis proceeds by a highly regulated series of events. In animals, maternal factors that accumulate in the egg cytoplasm control cell cycle progression at the initial stage of cleavage. However, cell cycle regulation is switched to a system governed by the activated nuclear genome at a specific stage of development, referred to as maternal-to-zygotic transition (MZT). Detailed molecular analyses have been performed on maternal factors and activated zygotic genes in MZT in mammals, fishes and chicken; however, the underlying mechanisms remain unclear in quail. In the present study, we demonstrated that MZT occurred at blastoderm stage V in the Japanese quail using novel gene targeting technology in which the CRISPR/Cas9 and intracytoplasmic sperm injection (ICSI) systems were combined. At blastoderm stage V, we found that maternal retinoblastoma 1 (RB1) protein expression was down-regulated, whereas the gene expression of cyclin D1 (CCND1) was initiated. When a microinjection of sgRNA containing CCND1-targeted sequencing and Cas9 mRNA was administered at the pronuclear stage, blastoderm development stopped at stage V and the down-regulation of RB1 did not occur. This result indicates the most notable difference from mammals in which CCND-knockout embryos are capable of developing beyond MZT. We also showed that CCND1 induced the phosphorylation of the serine/threonine residues of the RB1 protein, which resulted in the degradation of this protein. These results suggest that CCND1 is one of the key factors for RB1 protein degradation at MZT, and the elimination of RB1 may contribute to cell cycle progression after MZT during blastoderm development in the Japanese quail. Our novel technology, which combined the CRISPR/Cas9 system and ICSI, has the potential to become a powerful tool for avian-targeted mutagenesis.

Systematic screening of long intergenic noncoding RNAs expressed during chicken embryogenesis

Long noncoding RNAs (lncRNAs) have emerged as important regulators of many biological processes, including embryogenesis and development. To provide a systematic analysis of lncRNAs expressed during chicken embryogenesis, we used Iso-Seq and RNA-Seq to identify potential lncRNAs at embryonic stages from d 1 to d 8 of incubation: sequential stages covering gastrulation, somitogenesis, and organogenesis. The data characterized an expanded landscape of lncRNAs, yielding 45,410 distinct lncRNAs (31,282 genes). Amongst these, a set of 13,141 filtered intergenic lncRNAs (lincRNAs) transcribed from 9803 lincRNA gene loci, of which, 66.5% were novel, were further analyzed. These lincRNAs were found to share many characteristics with mammalian lincRNAs, including relatively short lengths, fewer exons, lower expression levels, and stage-specific expression patterns. Functional studies motivated by "guilt-by-association" associated individual lincRNAs with specific GO functions, providing an important resource for future studies of lincRNA function. Most importantly, a weighted gene co-expression network analysis suggested that genes of the brown module were specifically associated with the day 2 stage. LincRNAs within this module were co-expressed with proteins involved in hematopoiesis and lipid metabolism. This study presents the systematic identification of lincRNAs in developing chicken embryos and will serve as a powerful resource for the study of lincRNA functions.

Whole-Transcriptome Sequencing-Based Analysis of DAZL and Its Interacting Genes during Germ Cells Specification and Zygotic Genome Activation in Chickens

The deleted in azoospermia like (DAZL) is required for germ cells development and maintenance. In chickens, the mRNA and protein of DAZL, a representative maternally inherited germ plasm factor, are detected in the germ plasm of oocyte, zygote, and all stages of the intrauterine embryos. However, it is still insufficient to explain the origin and specification process of chicken germ cells, because the stage at which the zygotic transcription of DAZL occurs and the stage at which the maternal DAZL RNA/protein clears have not yet been fully identified. Moreover, a comprehensive understanding of the expression of DAZL interacting genes during the germ cells specification and development and zygotic genome activation (ZGA) is lacking in chickens. In this study, we identified a set of DAZL interacting genes in chickens using in silico prediction method. Then, we analyzed the whole-transcriptome sequencing (WTS)-based expression of DAZL and its interacting genes in the chicken oocyte, zygote, and Eyal-Giladi and Kochav (EGK) stage embryos (EGK.I to EGK.X). In the results, DAZL transcripts are increased in the zygote (onset of transcription), maintained the increased level until EGK.VI, and decreased from EGK.VIII (possible clearance of maternal RNAs). Among the DAZL interacting genes, most of them are increased either at 1st ZGA or 2nd ZGA, indicating their involvement in germ cells specification and development.

Zygotic genome activation in the chicken: a comparative review

Maternal RNAs and proteins in the oocyte contribute to early embryonic development. After fertilization, these maternal factors are cleared and embryonic development is determined by an individual's own RNAs and proteins, in a process called the maternal-to-zygotic transition. Zygotic transcription is initially inactive, but is eventually activated by maternal transcription factors. The timing and molecular mechanisms involved in zygotic genome activation (ZGA) have been well-described in many species. Among birds, a transcriptome-based understanding of ZGA has only been explored in chickens by RNA sequencing of intrauterine embryos. RNA sequencing of chicken intrauterine embryos, including oocytes, zygotes, and Eyal-Giladi and Kochav (EGK) stages I-X has enabled the identification of differentially expressed genes between consecutive stages. These studies have revealed that there are two waves of ZGA: a minor wave at the one-cell stage (shortly after fertilization) and a major wave between EGK.III and EGK.VI (during cellularization). In the chicken, the maternal genome is activated during minor ZGA and the paternal genome is quiescent until major ZGA to avoid transcription from supernumerary sperm nuclei. In this review, we provide a detailed overview of events in intrauterine embryonic development in birds (and particularly in chickens), as well as a transcriptome-based analysis of ZGA.

The first whole transcriptomic exploration of pre-oviposited early chicken embryos using single and bulked embryonic RNA-sequencing

Background

The chicken is a valuable model organism, especially in evolutionary and embryology research because its embryonic development occurs in the egg. However, despite its scientific importance, no transcriptome data have been generated for deciphering the early developmental stages of the chicken because of practical and technical constraints in accessing pre-oviposited embryos.

Findings

Here, we determine the entire transcriptome of pre-oviposited avian embryos, including oocyte, zygote, and intrauterine embryos from Eyal-giladi and Kochav stage I (EGK.I) to EGK.X collected using a noninvasive approach for the first time. We also compare RNA-sequencing data obtained using a bulked embryo sequencing and single embryo/cell sequencing technique. The raw sequencing data were preprocessed with two genome builds, Galgal4 and Galgal5, and the expression of 17,108 and 26,102 genes was quantified in the respective builds. There were some differences between the two techniques, as well as between the two genome builds, and these were affected by the emergence of long intergenic noncoding RNA annotations.

Conclusion

The first transcriptome datasets of pre-oviposited early chicken embryos based on bulked and single embryo sequencing techniques will serve as a valuable resource for investigating early avian embryogenesis, for comparative studies among vertebrates, and for novel gene annotation in the chicken genome.

Zygotic gene activation in the chicken occurs in two waves, the first involving only maternally derived genes

The first wave of transcriptional activation occurs after fertilisation in a species-specific pattern. Despite its importance to initial embryonic development, the characteristics of transcription following fertilisation are poorly understood in Aves. Here, we report detailed insights into the onset of genome activation in chickens. We established that two waves of transcriptional activation occurred, one shortly after fertilisation and another at Eyal-Giladi and Kochav Stage V. We found 1544 single nucleotide polymorphisms across 424 transcripts derived from parents that were expressed in offspring during the early embryonic stages. Surprisingly, only the maternal genome was activated in the zygote, and the paternal genome remained silent until the second-wave, regardless of the presence of a paternal pronucleus or supernumerary sperm in the egg. The identified maternal genes involved in cleavage that were replaced by bi-allelic expression. The results demonstrate that only maternal alleles are activated in the chicken zygote upon fertilisation, which could be essential for early embryogenesis and evolutionary outcomes in birds.

The early stages of animal development involve a handover of genetic control. Initially, the egg cell is maintained by genetic information inherited from the mother, but soon after fertilization it starts to depend on its own genes instead. Activating genes inside the fertilized egg cell (zygote) so that they can take control of development is known as zygotic genome activation.

Despite the fact that birds are often used to study how embryos develop, zygotic genome activation in birds is not well understood. Fertilization in birds, including chickens, is different to mammals in that it requires multiple sperm to fertilize an egg cell. As such, zygotic genome activation in birds is likely to differ from that in mammals.

By examining gene expression in embryos from mixed-breed chickens, Hwang, Seo et al. showed that there are two stages of zygotic genome activation in chickens. The genes derived from the mother become active in the first stage, while genes from the father become active in the second stage. Genome activation in birds is therefore very different to the same process in mammals, which involves genome activation of both parents from the first stage. This extra level of control may help to prevent genetic complications resulting from the presence of multiple sperm, each of which carries a different set of genes from the father.