Exploratory bioinformatics investigation reveals importance of "junk" DNA in early embryo development

Nr5a2 is dispensable for zygotic genome activation but essential for morula development

Early embryogenesis is driven by transcription factors (TFs) that first activate the zygotic genome and then specify the lineages constituting the blastocyst. Although the TFs specifying the blastocyst's lineages are well characterized, those playing earlier roles remain poorly defined. Using mouse models of the TF Nr5a2, we show that Nr5a2-/- embryos arrest at the early morula stage and exhibit altered lineage specification, frequent mitotic failure, and substantial chromosome segregation defects. Although NR5A2 plays a minor but measurable role during zygotic genome activation, it predominantly acts as a master regulator at the eight-cell stage, controlling expression of lineage-specifying TFs and genes involved in mitosis, telomere maintenance, and DNA repair. We conclude that NR5A2 coordinates proliferation, genome stability, and lineage specification to ensure correct morula development.

Rise and SINE: roles of transcription factors and retrotransposons in zygotic genome activation

In sexually reproducing organisms, life begins with the fusion of transcriptionally silent gametes, the oocyte and sperm. Although initiation of transcription in the embryo, known as zygotic genome activation (ZGA), is universally required for development, the transcription factors regulating this process are poorly conserved. In this Perspective, we discuss recent insights into the mechanisms of ZGA in totipotent mammalian embryos, namely ZGA regulation by several transcription factors, including by orphan nuclear receptors (OrphNRs) such as the pioneer transcription factor NR5A2, and by factors of the DUX, TPRX and OBOX families. We performed a meta-analysis and compiled a list of pan-ZGA genes, and found that most of these genes are indeed targets of the above transcription factors. Remarkably, more than a third of these ZGA genes appear to be regulated both by OrphNRs such as NR5A2 and by OBOX proteins, whose motifs co-occur in SINE B1 retrotransposable elements, which are enriched near ZGA genes. We propose that ZGA in mice is activated by recruitment of multiple transcription factors to SINE B1 elements that function as enhancers, and discuss a potential relevance of this mechanism to Alu retrotransposable elements in human ZGA.

Kick-starting the zygotic genome: licensors, specifiers, and beyond

Zygotic genome activation (ZGA), the first transcription event following fertilization, kickstarts the embryonic program that takes over the control of early development from the maternal products. How ZGA occurs, especially in mammals, is poorly understood due to the limited amount of research materials. With the rapid development of single-cell and low-input technologies, remarkable progress made in the past decade has unveiled dramatic transitions of the epigenomes, transcriptomes, proteomes, and metabolomes associated with ZGA. Moreover, functional investigations are yielding insights into the key regulators of ZGA, among which two major classes of players are emerging: licensors and specifiers. Licensors would control the permission of transcription and its timing during ZGA. Accumulating evidence suggests that such licensors of ZGA include regulators of the transcription apparatus and nuclear gatekeepers. Specifiers would instruct the activation of specific genes during ZGA. These specifiers include key transcription factors present at this stage, often facilitated by epigenetic regulators. Based on data primarily from mammals but also results from other species, we discuss in this review how recent research sheds light on the molecular regulation of ZGA and its executors, including the licensors and specifiers.

Obox4 promotes zygotic genome activation upon loss of Dux

Once fertilized, mouse zygotes rapidly proceed to zygotic genome activation (ZGA), during which long terminal repeats (LTRs) of murine endogenous retroviruses with leucine tRNA primer (MERVL) are activated by a conserved homeodomain-containing transcription factor, DUX. However, Dux-knockout embryos produce fertile mice, suggesting that ZGA is redundantly driven by an unknown factor(s). Here, we present multiple lines of evidence that the multicopy homeobox gene, Obox4, encodes a transcription factor that is highly expressed in mouse two-cell embryos and redundantly drives ZGA. Genome-wide profiling revealed that OBOX4 specifically binds and activates MERVL LTRs as well as a subset of murine endogenous retroviruses with lysine tRNA primer (MERVK) LTRs. Depletion of Obox4 is tolerated by embryogenesis, whereas concomitant Obox4/Dux depletion markedly compromises embryonic development. Our study identified OBOX4 as a transcription factor that provides genetic redundancy to preimplantation development.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

NR5A2 connects zygotic genome activation to the first lineage segregation in totipotent embryos

Zygotic genome activation (ZGA) marks the beginning of the embryonic program for a totipotent embryo, which gives rise to the inner cell mass (ICM) where pluripotent epiblast arises, and extraembryonic trophectoderm. However, how ZGA is connected to the first lineage segregation in mammalian embryos remains elusive. Here, we investigated the role of nuclear receptor (NR) transcription factors (TFs), whose motifs are highly enriched and accessible from the 2-cell (2C) to 8-cell (8C) stages in mouse embryos. We found that NR5A2, an NR TF strongly induced upon ZGA, was required for this connection. Upon Nr5a2 knockdown or knockout, embryos developed beyond 2C normally with the zygotic genome largely activated. However, 4-8C-specific gene activation was substantially impaired and Nr5a2-deficient embryos subsequently arrested at the morula stage. Genome-wide chromatin binding analysis showed that NR5A2-bound cis-regulatory elements in both 2C and 8C embryos are strongly enriched for B1 elements where its binding motif is embedded. NR5A2 was not required for the global opening of its binding sites in 2C embryos but was essential to the opening of its 8C-specific binding sites. These 8C-specific, but not 2C-specific, binding sites are enriched near genes involved in blastocyst and stem cell regulation, and are often bound by master pluripotency TFs in blastocysts and embryonic stem cells (ESCs). Importantly, NR5A2 regulated key pluripotency genes Nanog and Pou5f1/Oct4, and primitive endoderm regulatory genes including Gata6 among many early ICM genes, as well as key trophectoderm regulatory genes including Tead4 and Gata3 at the 8C stage. By contrast, master pluripotency TFs NANOG, SOX2, and OCT4 targeted both early and late ICM genes in mouse ESCs. Taken together, these data identify NR5A2 as a key regulator in totipotent embryos that bridges ZGA to the first lineage segregation during mouse early development.

Hypoxia-induced circRTN4IP1 promotes progression and glycolysis of hepatocellular carcinoma cells

Hypoxia is one of the hallmarks of solid tumors, especially in hepatocellular carcinoma (HCC). CircRNAs are reported to be tightly connected to hypoxia and also have essential roles in cancer progression. However, many circRNAs implicated in hypoxia-mediated HCC progression are still unclear and require further exploration. In this study, a hypoxia cell model was structured by exposing cells to hypoxia conditions (1% O2) and normoxia conditions (21% O2) as a control. The effects of hypoxia and normoxia on cell viability, migration, invasion, and glycolysis were examined. The expressions of circRNARTN4IP1 under hypoxia were identified. Finally, molecular mechanisms and biological function of circRTN4IP1 were explored. We confirmed that hypoxia treatment facilitated capacities of proliferation, migration, invasion, and glycolysis in tumor cells. Hypoxia induced a significant increase expression of circRTN4IP1 in cells. Functionally, knockdown of circRTN4IP1 inhibited cell malignant progression and glycolysis under hypoxia HCC cells. Mechanistically, HIF1A targeted the promoter region of circRTN4IP1 and positively regulated the expression of circRTN4IP1. In addition, circRTN4IP1 targeted miR-532-5p/G6PC3 axis. In short, hypoxia induced activation of the HIF1A/circRTN4IP1/miR-532-5p/G6PC3 signaling axis, which promoted proliferation, migration, invasion, and glycolysis of HCC cells. This study may reveal a possible mechanism driving the progression of hypoxia HCC, so as to find potential effective candidates for targeting hypoxia microenvironment therapy.

The composition dynamics of transposable elements in human blastocysts

Transposable elements (TEs) are mobile DNA sequences that can replicate themselves and play significant roles in embryo development and chromosomal structure remodeling. In this study, we investigated the variation of TEs in blastocysts with different parental genetic backgrounds. We analyzed the proportions of 1137 TEs subfamilies from six classes at the DNA level using Bowtie2 and PopoolationTE2 in 196 blastocysts with abnormal parental chromosomal diseases. Our findings revealed that the parental karyotype was the dominant factor influencing TEs frequencies. Out of the 1116 subfamilies, different frequencies were observed in blastocysts with varying parental karyotypes. The development stage of blastocysts was the second most crucial factor influencing TEs proportions. A total of 614 subfamilies exhibited different proportions at distinct blastocyst stages. Notably, subfamily members belonging to the Alu family showed a high proportion at stage 6, while those from the LINE class exhibited a high proportion at stage 3 and a low proportion at stage 6. Moreover, the proportions of some TEs subfamilies also varied depending on blastocyst karyotype, inner cell mass status, and outer trophectoderm status. We found that 48 subfamilies displayed different proportions between balanced and unbalanced blastocysts. Additionally, 19 subfamilies demonstrated varying proportions among different inner cell mass scores, and 43 subfamilies exhibited different proportions among outer trophectoderm scores. This study suggests that the composition of TEs subfamilies may be influenced by various factors and undergoes dynamic modulation during embryo development.

The single-cell expression profile of transposable elements and transcription factors in human early biparental and uniparental embryonic development

Transposable elements (TEs) and transcription factors (TFs) are involved in the precise regulation of gene expression during the preimplantation stage. Activation of TEs is a key event for mammalian embryonic genome activation and preimplantation early embryonic development. TFs are involved in the regulation of drastic changes in gene expression patterns, but an inventory of the interplay between TEs and TFs during normal/abnormal human embryonic development is still lacking. Here we used single-cell RNA sequencing data generated from biparental and uniparental embryos to perform an integrative analysis of TE and TF expression. Our results showed that endogenous retroviruses (ERVs) are mainly expressed during the minor embryonic genome activation (EGA) process of early embryos, while Alu is gradually expressed in the middle and later stages. Some important ERVs (e.g., LTR5_Hs, MLT2A1) and Alu TEs are expressed at significantly lower levels in androgenic embryos. Integrative analysis revealed that the expression of the transcription factors CTCF and POU5F1 is correlated with the differential expression of ERV TEs. Comparative coexpression network analysis further showed distinct expression levels of important TFs (e.g., LEUTX and ZSCAN5A) in dizygotic embryos vs. parthenogenetic and androgenic embryos. This systematic investigation of TE and TF expression in human early embryonic development by single-cell RNA sequencing provides valuable insights into mammalian embryonic development.

Suppression of endogenous retroviral enhancers in mouse embryos derived from somatic cell nuclear transfer

Endogenous retroviruses (ERVs) in the mammalian genome play diverse roles in embryonic development. These developmentally related ERVs are generally repressed in somatic cells and therefore are likely repressed in embryos derived from somatic cell nuclear transfer (SCNT). In this study, we sought to identify ERVs that are repressed in SCNT-derived morulae, which might cause previously unexplained embryonic deaths shortly after implantation. Our transcriptome analysis revealed that, amongst ERV families, ERVK was specifically, and strongly downregulated in SCNT-derived embryos while other transposable elements including LINE and ERVL were unchanged. Among the subfamilies of ERVK, RLTR45-int was most repressed in SCNT-derived embryos despite its highest expression in control fertilized embryos. Interestingly, the nearby genes (within 5–50 kb, n = 18; 50–200 kb, n = 63) of the repressed RLTR45-int loci were also repressed in SCNT-derived embryos, with a significant correlation between them. Furthermore, lysine H3K27 acetylation was enriched around the RLTR45-int loci. These findings indicate that RLTR45-int elements function as enhancers of nearby genes. Indeed, deletion of two sequential RLTR45-int loci on chromosome 4 or 18 resulted in downregulations of nearby genes at the morula stage. We also found that RLTR45-int loci, especially SCNT-low, enhancer-like loci, were strongly enriched with H3K9me3, a repressive histone mark. Importantly, these H3K9me3-enriched regions were not activated by overexpression of H3K9me3 demethylase Kdm4d in SCNT-derived embryos, suggesting the presence of another epigenetic barrier repressing their expressions and enhancer activities in SCNT embryos. Thus, we identified ERVK subfamily RLTR45-int, putative enhancer elements, as a strong reprogramming barrier for SCNT (253 words).

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 Role of Hypoxia-Associated Long Non-Coding RNAs in Breast Cancer

Breast cancer is the leading cause of cancer-related deaths in women worldwide. In the United States, even with earlier diagnosis and treatment improvements, the decline in mortality has stagnated in recent years. More research is needed to provide better diagnostic, prognostic, and therapeutic tools for these patients. Long non-coding RNAs are newly described molecules that have extensive roles in breast cancer. Emerging reports have shown that there is a strong link between these RNAs and the hypoxic response of breast cancer cells, which may be an important factor for enhanced tumoral progression. In this review, we summarize the role of hypoxia-associated lncRNAs in the classic cancer hallmarks, describing their effects on the upstream and downstream hypoxia signaling pathway and the use of them as diagnostic and prognostic tools.

Locus-specific expression of transposable elements in single cells with CELLO-seq

Transposable elements (TEs) regulate diverse biological processes, from early development to cancer. Expression of young TEs is difficult to measure with next-generation, single-cell sequencing technologies because their highly repetitive nature means that short complementary DNA reads cannot be unambiguously mapped to a specific locus. Single CELl LOng-read RNA-sequencing (CELLO-seq) combines long-read single cell RNA-sequencing with computational analyses to measure TE expression at unique loci. We used CELLO-seq to assess the widespread expression of TEs in two-cell mouse blastomeres as well as in human induced pluripotent stem cells. Across both species, old and young TEs showed evidence of locus-specific expression with simulations demonstrating that only a small number of very young elements in the mouse could not be mapped back to the reference with high confidence. Exploring the relationship between the expression of individual elements and putative regulators revealed large heterogeneity, with TEs within a class showing different patterns of correlation and suggesting distinct regulatory mechanisms.

Expression and function of PDGF-C in development and stem cells

Platelet-derived growth factor C (PDGF-C) is a relatively new member of the PDGF family, discovered nearly 20 years after the finding of platelet-derived growth factor A (PDGF-A) and platelet-derived growth factor B (PDGF-B). PDGF-C is generally expressed in most organs and cell types. Studies from the past 20 years have demonstrated critical roles of PDGF-C in numerous biological, physiological and pathological processes, such as development, angiogenesis, tumour growth, tissue remodelling, wound healing, atherosclerosis, fibrosis, stem/progenitor cell regulation and metabolism. Understanding PDGF-C expression and activities thus will be of great importance to various research disciplines. In this review, however, we mainly discuss the expression and functions of PDGF-C and its receptors in development and stem cells.

Dynamics of small non-coding RNAs in bovine scNT embryos through the maternal-to-embryonic transition

The efficiency of somatic cell nuclear transfer (scNT) for production of viable offspring is relatively low as compared to in vitro fertilization (IVF), presumably due to deficiencies in epigenetic reprogramming of the donor cell genome. Such defects may also involve the population of small non-coding RNAs (sncRNAs), which are important during early embryonic development. The objective of this study was to examine dynamic changes in relative abundance of sncRNAs during the maternal-to embryonic transition (MET) in bovine embryos produced by scNT as compared to IVF by using RNA sequencing. When comparing populations of miRNA in scNT versus IVF embryos, only miR-2340, miR-345, and miR34a were differentially expressed in morulae, though many more miRNAs were differentially expressed when comparing across developmental stages. Also of interest, distinct populations of piwi-interacting like RNAs (pilRNAs) were identified in bovine embryos prior to and during embryonic genome activation (EGA) as compared bovine embryos post EGA and differentiated cells. Overall, sncRNA sequencing analysis of preimplantation embryos revealed largely similar profiles of sncRNAs for IVF and scNT embryos at the 2-cell, 8-cell, morula and blastocyst stages of development. However, these sncRNA profiles, including miRNA, piRNA and tRNA fragments, were notably distinct prior to and after completion of the MET.

The methylome of a human polar body reflects that of its sibling oocyte and its aberrance may indicate poor embryo development

STUDY QUESTION

Is it possible to evaluate the methylome of individual oocytes to investigate the DNA methylome alterations in metaphase II (MII) oocytes with reduced embryo developmental potential?

SUMMARY ANSWER

The DNA methylome of each human first polar body (PB1) closely mirrored that of its sibling MII oocyte; hypermethylated long interspersed nuclear element (LINE) and long terminal repeats (LTRs) and methylation aberrations in PB1 promoter regions may indicate poor embryo development.

WHAT IS KNOWN ALREADY

The developmental potential of an embryo is determined by the oocyte's developmental competence, and the PB1 is a good substitute to examine the chromosomal status of the corresponding oocyte. However, DNA methylation, a key epigenetic modification, also regulates gene expression and embryo development.

STUDY DESIGN, SIZE, DURATION

Twelve pairs of PB1s and sibling MII oocytes were biopsied and sequenced to compare their methylomes. To further investigate the methylome of PB1s and the potential epigenetic factors that may affect oocyte quality, MII oocytes (n = 74) were fertilized through ICSI, while PB1s were biopsied and profiled to measure DNA methylation. The corresponding embryos were further cultured to track their development potential. The oocytes and sperm samples used in this study were donated by healthy volunteers with signed informed consent.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Single-cell methylome sequencing was applied to obtain the DNA methylation profiles of PB1s and oocytes. The DNA methylome of PB1s was compared between the respective group of oocytes that progressed to blastocysts and the group of oocytes that failed to develop. DNA methylation levels of corresponding regions and differentially methylated regions were calculated using customized Perl and R scripts. RNA-seq data were downloaded from a previously published paper and reanalysed.

MAIN RESULTS AND THE ROLE OF CHANCE

The results from PB1-MII oocyte pair validated that PB1 contains nearly the same methylome (average Pearson correlation is 0.92) with sibling MII oocyte. LINE and LTR expression increased markedly after fertilization. Moreover, the DNA methylation levels in LINE (including LINE1 and LINE2) and LTR were significantly higher in the PB1s of embryos that could not reach the blastocyst stage (Wilcoxon-Matt-Whitney test, P < 0.05). DNA methylation in PB1 promoters correlated negatively with gene expression of MII oocyte. Regarding the methylation status of the promoter regions, 66 genes were hypermethylated in the developmental arrested group, with their related functions (significantly enriched in several Gene Ontology terms) including transcription, positive regulation of adenylate cyclase activity, mitogen-activated protein kinase (MAPK) cascade and intracellular oestrogen receptor signalling pathway.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

Data analysis performed in this study focused on the competence of human oocytes and compared them with maternal genetic and epigenetic profiles. Therefore, data regarding the potential regulatory roles of paternal genomes in embryo development are lacking.

WIDER IMPLICATIONS OF THE FINDINGS

The results from PB1-oocyte pairs demonstrated that PB1s shared similar methylomes with their sibling oocytes. The selection of the good embryos for transfer should not only rely on morphology but also consider the DNA methylation of the corresponding PB1 and therefore MII oocyte. The application of early-stage analysis of PB1 offers an option for high-quality oocyte and embryo selection, which provides an additional tool for elective single embryo transfer in assisted reproduction.

STUDY FUNDING/COMPETING INTEREST(S)

This study was supported by the National Key Research and Development Program of China (2018YFC1004003, 2017YFA0103801), the National Natural Science Foundation of China (81730038, 3187144, 81521002) and the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA16020703). The authors have no conflicts of interest to declare.

The Regulation and Functions of Endogenous Retrovirus in Embryo Development and Stem Cell Differentiation

Endogenous retroviruses (ERVs) are repetitive sequences in the genome, belonging to the retrotransposon family. During the course of life, ERVs are associated with multiple aspects of chromatin and transcriptional regulation in development and pathological conditions. In mammalian embryos, ERVs are extensively activated in early embryo development, but with a highly restricted spatial-temporal pattern; and they are drastically silenced during differentiation with exceptions in extraembryonic tissue and germlines. The dynamic activation pattern of ERVs raises questions about how ERVs are regulated in the life cycle and whether they are functionally important to cell fate decision during early embryo and somatic cell development. Therefore, in this review, we focus on the pieces of evidence demonstrating regulations and functions of ERVs during stem cell differentiation, which suggests that ERV activation is not a passive result of cell fate transition but the active epigenetic and transcriptional regulation during mammalian development and stem cell differentiation.

Transposable element-derived sequences in vertebrate development

Transposable elements (TEs) are major components of all vertebrate genomes that can cause deleterious insertions and genomic instability. However, depending on the specific genomic context of their insertion site, TE sequences can sometimes get positively selected, leading to what are called "exaptation" events. TE sequence exaptation constitutes an important source of novelties for gene, genome and organism evolution, giving rise to new regulatory sequences, protein-coding exons/genes and non-coding RNAs, which can play various roles beneficial to the host. In this review, we focus on the development of vertebrates, which present many derived traits such as bones, adaptive immunity and a complex brain. We illustrate how TE-derived sequences have given rise to developmental innovations in vertebrates and how they thereby contributed to the evolutionary success of this lineage.

Transcript assembly improves expression quantification of transposable elements in single-cell RNA-seq data

Transposable elements (TEs) are an integral part of the host transcriptome. TE-containing noncoding RNAs (ncRNAs) show considerable tissue specificity and play important roles during development, including stem cell maintenance and cell differentiation. Recent advances in single-cell RNA-seq (scRNA-seq) revolutionized cell type-specific gene expression analysis. However, effective scRNA-seq quantification tools tailored for TEs are lacking, limiting our ability to dissect TE expression dynamics at single-cell resolution. To address this issue, we established a TE expression quantification pipeline that is compatible with scRNA-seq data generated across multiple technology platforms. We constructed TE-containing ncRNA references using bulk RNA-seq data and showed that quantifying TE expression at the transcript level effectively reduces noise. As proof of principle, we applied this strategy to mouse embryonic stem cells and successfully captured the expression profile of endogenous retroviruses in single cells. We further expanded our analysis to scRNA-seq data from early stages of mouse embryogenesis. Our results illustrated the dynamic TE expression at preimplantation stages and revealed 146 TE-containing ncRNA transcripts with substantial tissue specificity during gastrulation and early organogenesis.

Sex-specific epigenetic profile of inner cell mass of mice conceived in vivo or by IVF

The preimplantation stage of development is exquisitely sensitive to environmental stresses, and changes occurring during this developmental phase may have long-term health effects. Animal studies indicate that IVF offspring display metabolic alterations, including hypertension, glucose intolerance and cardiac hypertrophy, often in a sexual dimorphic fashion. The detailed nature of epigenetic changes following in-vitro culture is, however, unknown. This study was performed to evaluate the epigenetic (using whole-genome bisulfite sequencing (WGBS) and assay for transposase-accessible chromatin using sequencing (ATAC-seq)) and transcriptomic changes (using RNA-seq) occurring in the inner cell mass (ICM) of male or female mouse embryos generated in vivo or by IVF. We found that the ICM of IVF embryos, compared to the in-vivo ICM, differed in 3% of differentially methylated regions (DMRs), of which 0.1% were located on CpG islands. ATAC-seq revealed that 293 regions were more accessible and 101 were less accessible in IVF embryos, while RNA-seq revealed that 21 genes were differentially regulated in IVF embryos. Functional enrichment analysis revealed that stress signalling (STAT and NF-kB signalling), developmental processes and cardiac hypertrophy signalling showed consistent changes in WGBS and ATAC-seq platforms. In contrast, male and female embryos showed minimal changes. Male ICM had an increased number of significantly hyper-methylated DMRs, while only 27 regions showed different chromatin accessibility and only one gene was differentially expressed. In summary, this study provides the first comprehensive analysis of DNA methylation, chromatin accessibility and RNA expression changes induced by IVF in male and female ICMs. This dataset can be of value to all researchers interested in the developmental origin of health and disease (DOHaD) hypothesis and might lead to a better understanding of how early embryonic manipulation may affect adult health.

Analyzing a putative enhancer of optic disc morphology

Background

Genome-wide association studies have identified the CDC7-TGFBR3 intergenic region on chromosome 1 to be strongly associated with optic disc area size. The mechanism of its function remained unclear until new data on eQTL markers emerged from the Genotype-Tissue Expression project. The target region was found to contain a strong silencer of the distal (800 kb) Transcription Factor (TF) gene GFI1 (Growth Factor Independent Transcription Repressor 1) specifically in neuroendocrine cells (pituitary gland). GFI1 has also been reported to be involved in the development of sensory neurons and hematopoiesis. Therefore, GFI1, being a developmental gene, is likely to affect optic disc area size by altering the expression of the associated genes via long-range interactions.

Results

Distribution of haplotypes in the putative enhancer region has been assessed using the data on four continental supergroups generated by the 1000 Genomes Project. The East Asian (EAS) populations were shown to manifest a highly homogenous unimodal haplotype distribution pattern within the region with the major haplotype occurring with the frequency of 0.9. Another European specific haplotype was observed with the frequency of 0.21. The major haplotype appears to be involved in silencing GFI1repressor gene expression, which might be the cause of increased optic disc area characteristic of the EAS populations.

The enhancer/eQTL region overlaps AluJo element, which implies that this particular regulatory element is primate-specific and confined to few tissues.

Conclusion

Population specific distribution of GFI1 enhancer alleles may predispose certain ethnic groups to glaucoma.

Duplication and divergence of the retrovirus restriction gene Fv1 in Mus caroli allows protection from multiple retroviruses

Viruses and their hosts are locked in an evolutionary race where resistance to infection is acquired by the hosts while viruses develop strategies to circumvent these host defenses. Forming one arm of the host defense armory are cell autonomous restriction factors like Fv1. Originally described as protecting laboratory mice from infection by murine leukemia virus (MLV), Fv1s from some wild mice have also been found to restrict non-MLV retroviruses, suggesting an important role in the protection against viruses in nature. We surveyed the Fv1 genes of wild mice trapped in Thailand and characterized their restriction activities against a panel of retroviruses. An extra copy of the Fv1 gene, named Fv7, was found on chromosome 6 of three closely related Asian species of mice: Mus caroli, M. cervicolor, and M. cookii. The presence of flanking repeats suggested it arose by LINE-mediated retroduplication within their most recent common ancestor. A high degree of natural variation was observed in both Fv1 and Fv7 and, on top of positive selection at certain residues, insertions and deletions were present that changed the length of the reading frames. These genes exhibited a range of restriction phenotypes, with activities directed against gamma-, spuma-, and lentiviruses. It seems likely, at least in the case of M. caroli, that the observed gene duplication may expand the breadth of restriction beyond the capacity of Fv1 alone and that one or more such viruses have recently driven or continue to drive the evolution of the Fv1 and Fv7 genes.

Investigating the gene expression profiles of cells in seven embryonic stages with machine learning algorithms

The development of embryonic cells involves several continuous stages, and some genes are related to embryogenesis. To date, few studies have systematically investigated changes in gene expression profiles during mammalian embryogenesis. In this study, a computational analysis using machine learning algorithms was performed on the gene expression profiles of mouse embryonic cells at seven stages. First, the profiles were analyzed through a powerful Monte Carlo feature selection method for the generation of a feature list. Second, increment feature selection was applied on the list by incorporating two classification algorithms: support vector machine (SVM) and repeated incremental pruning to produce error reduction (RIPPER). Through SVM, we extracted several latent gene biomarkers, indicating the stages of embryonic cells, and constructed an optimal SVM classifier that produced a nearly perfect classification of embryonic cells. Furthermore, some interesting rules were accessed by the RIPPER algorithm, suggesting different expression patterns for different stages.

Retroelement Insertion in a CRISPR/Cas9 Editing Site in the Early Embryo Intensifies Genetic Mosaicism

Continued CRISPR/Cas9-mediated editing activity that allows differential and asynchronous modification of alleles in successive cell generations expands allelic complexity. To understand the earliest events during CRISPR/Cas9 editing and the allelic selection among the progeny of subsequent cell divisions, we inspected in detail the genotypes of 4- and 8-cell embryos and embryonic stem cells (ESCs) after microinjection of a CRISPR toolkit into the zygotes. We found a higher editing frequency in 8-cell embryos than in 4-cell embryos, indicating that the CRISPR/Cas9 activity persisted through the 8-cell stage. Analysis of a CRISPR/Cas9 transgenic founder mouse revealed that four different alleles were present in its organs in different combinations and that its germline included three different mutant alleles, as shown by the genotypes of the pups. The indel depth, which measured the extent of indels at the sequence level within single embryos, decreased significantly as the embryos advanced to form ESCs, suggesting that exclusion of fatal indels occurred in the subsequent cell generations. Interestingly, we discovered that the CRISPR sites frequently contained introduced retroelement sequences and that this occurred preferentially with certain classes of retroelements. Therefore, in addition to CRISPR/Cas9's innate mechanism of separate, differential enzymatic modifications of alleles, the frequent retroelement insertions that occur in early mouse embryos during CRISPR/Cas9 editing further expand the allelic diversity and mosaicism in the resulting transgenic founders.

Pleomorphic Adenoma Gene 1 Is Needed For Timely Zygotic Genome Activation and Early Embryo Development

Pleomorphic adenoma gene 1 (PLAG1) is a transcription factor involved in cancer and growth. We discovered a de novo DNA motif containing a PLAG1 binding site in the promoters of genes activated during zygotic genome activation (ZGA) in human embryos. This motif was located within an Alu element in a region that was conserved in the murine B1 element. We show that maternally provided Plag1 is needed for timely mouse preimplantation embryo development. Heterozygous mouse embryos lacking maternal Plag1 showed disrupted regulation of 1,089 genes, spent significantly longer time in the 2-cell stage, and started expressing Plag1 ectopically from the paternal allele. The de novo PLAG1 motif was enriched in the promoters of the genes whose activation was delayed in the absence of Plag1. Further, these mouse genes showed a significant overlap with genes upregulated during human ZGA that also contain the motif. By gene ontology, the mouse and human ZGA genes with de novo PLAG1 motifs were involved in ribosome biogenesis and protein synthesis. Collectively, our data suggest that PLAG1 affects embryo development in mice and humans through a conserved DNA motif within Alu/B1 elements located in the promoters of a subset of ZGA genes.

Expression dynamics of repetitive DNA in early human embryonic development

BACKGROUND

The last decade witnessed a number of genome-wide studies on human pre-implantation, which mostly focused on genes and provided only limited information on repeats, excluding the satellites. Considering the fact that repeats constitute a large portion of our genome with reported links to human physiology and disease, a thorough understanding of their spatiotemporal regulation during human embryogenesis will give invaluable clues on chromatin dynamics across time and space. Therefore, we performed a detailed expression analysis of all repetitive DNA elements including the satellites across stages of human pre-implantation and embryonic stem cells.

RESULTS

We uncovered stage-specific expressions of more than a thousand repeat elements whose expressions fluctuated with a mild global decrease at the blastocyst stage. Most satellites were highly expressed at the 4-cell level and expressions of ACRO1 and D20S16 specifically peaked at this point. Whereas all members of the SVA elements were highly upregulated at 8-cell and morula stages, other transposons and small RNA repeats exhibited a high level of variation among their specific subtypes. Our repeat enrichment analysis in gene promoters coupled with expression correlations highlighted potential links between repeat expressions and nearby genes, emphasising mostly 8-cell and morula specific genes together with SVA_D, LTR5_Hs and LTR70 transposons. The DNA methylation analysis further complemented the understanding on the mechanistic aspects of the repeatome's regulation per se and revealed critical stages where DNA methylation levels are negatively correlating with repeat expression.

CONCLUSIONS

Taken together, our study shows that specific expression patterns are not exclusive to genes and long non-coding RNAs but the repeatome also exhibits an intriguingly dynamic pattern at the global scale. Repeats identified in this study; particularly satellites, which were historically associated with heterochromatin, and those with potential links to nearby gene expression provide valuable insights into the understanding of key events in genomic regulation and warrant further research in epigenetics, genomics and developmental biology.

Heterogeneity of transposon expression and activation of the repressive network in human fetal germ cells

Epigenetic resetting in germ cells during development de-represses transposable elements (TEs). piRNAs protect fetal germ cells by targeted mRNA destruction and deposition of repressive epigenetic marks. Here, we provide the first evidence for an active piRNA pathway and TE repression in germ cells of human fetal testis. We identify pre-pachytene piRNAs with features of secondary amplification that map most abundantly to the long interspersed element type 1 (L1) family of TEs. L1-ORF1p expression is heterogeneous in fetal germ cells, peaks at mid-gestation and declines concomitantly with increases in piRNAs, nuclear localization of HIWI2 and an increase in H3K9me3. Surprisingly, the same cells with accumulation of L1-ORF1p display highest levels of HIWI2 and H3K9me3. Conversely, the earliest germ cells with high levels of L1-ORF1p express low levels of the chaperone HSP90α. We propose that a subset of germ cells resists L1 expression, whereas L1-expressing germ cells activate the repression pathway that leads to epigenetic silencing of L1 via H3K9me3.

A Comparison of Pregnancy Outcomes in Patients Undergoing Donor Egg Single Embryo Transfers With and Without Preimplantation Genetic Testing

Two of the many milestone developments in the field of assisted reproduction have been oocyte donation and preimplantation genetic testing for aneuploidy (PGT-A). Because it has been demonstrated that even young women produce a meaningful proportion of aneuploid embryos, screening out such abnormalities could potentially increase the efficacy of donor egg (DE) cycles. In this retrospective cohort study, we investigated the effect of PGT-A on DE cycle outcomes, including implantation rate (IR), spontaneous abortion rate (SABR), and ongoing pregnancy/live birth rate. We used fresh and frozen donor cycles not using PGT-A as comparison groups; all cases involved single embryo transfer. Data analysis revealed that PGT-A did not improve pregnancy outcome metrics in DE cycles, although there was a trend toward decreasing the SABR. There was a significant increase in IR with fresh cycles outperforming all frozen cycles. Overall, these results suggest that the benefits of performing PGT-A on embryos derived from young DEs may be limited and that there is an effect of the freezing process on pregnancy outcomes. These findings may provide useful insights into the science and practice of PGT-A across all of its applications.

iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data

BACKGROUND

RNA-seq is widely used for transcriptomic profiling, but the bioinformatics analysis of resultant data can be time-consuming and challenging, especially for biologists. We aim to streamline the bioinformatic analyses of gene-level data by developing a user-friendly, interactive web application for exploratory data analysis, differential expression, and pathway analysis.

RESULTS

iDEP (integrated Differential Expression and Pathway analysis) seamlessly connects 63 R/Bioconductor packages, 2 web services, and comprehensive annotation and pathway databases for 220 plant and animal species. The workflow can be reproduced by downloading customized R code and related pathway files. As an example, we analyzed an RNA-Seq dataset of lung fibroblasts with Hoxa1 knockdown and revealed the possible roles of SP1 and E2F1 and their target genes, including microRNAs, in blocking G1/S transition. In another example, our analysis shows that in mouse B cells without functional p53, ionizing radiation activates the MYC pathway and its downstream genes involved in cell proliferation, ribosome biogenesis, and non-coding RNA metabolism. In wildtype B cells, radiation induces p53-mediated apoptosis and DNA repair while suppressing the target genes of MYC and E2F1, and leads to growth and cell cycle arrest. iDEP helps unveil the multifaceted functions of p53 and the possible involvement of several microRNAs such as miR-92a, miR-504, and miR-30a. In both examples, we validated known molecular pathways and generated novel, testable hypotheses.

CONCLUSIONS

Combining comprehensive analytic functionalities with massive annotation databases, iDEP ( http://ge-lab.org/idep/ ) enables biologists to easily translate transcriptomic and proteomic data into actionable insights.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

Structure, clustering and functional insights of repeats configurations in the upstream promoter region of the human coding genes

Background

Repetitive DNA sequences (Repeats) are significant regions in the human genome that have a specific genomic distribution, structure, and several binding sites for genome architecture and function. In consequence, the possible configurations of Repeats in specific and dynamic regions like the gene promoters could define footprints for molecular mechanisms, pathways, and cell function beyond their density in the genome. Here we explored the distribution of Repeats in the upstream promoter region of the human coding genes with the aim to identify specific configurations, clusters and functional meaning of those elements. Our method includes structural descriptions, hierarchical clustering, pathway association, and functional enrichment analysis.

Results

We report here several configurations of Repeats in the upstream promoter region (UPR), which define 2729 patterns for the 80% of the human coding genes. There are 47 types of Repeats in these configurations, where the most frequent were Alu, Low_complexity, MIR, Simple_repeat, LINE/L2, LINE/L1, hAT-Charlie, and ERV1. The distribution, length, and the high frequency of Repeats in the UPR defines several patterns and clusters, where the minimum frequency of configuration among Repeats was higher than 0.7. We found those clusters associated with cellular pathways and ontologies; thus, it was plausible to determine groups of Repeats to specific functional insights, for example, pathways for Genetic Information Processing or Metabolism shows particular groups of Repeats with specific configurations.

Conclusion

Based on these findings, we propose that specific configurations of repetitive elements describe frequent patterns in the upstream promoter for sets of human coding genes, which those correlated to specific and essential cell pathways and functions.

Electronic supplementary material

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

Transposable Element Expression in Acute Myeloid Leukemia Transcriptome and Prognosis

Over half of the human genome is comprised of transposable elements (TE). Despite large-scale studies of the transcriptome in cancer, a comprehensive look at TE expression and its relationship to various mutations or prognosis has not been performed. We characterized the expression of TE in 178 adult acute myeloid leukemia (AML) patients using transcriptome data from The Cancer Genome Atlas. We characterized mutation specific dysregulation of TE expression using a multivariate linear model. We identified distinct patterns of TE expression associated with specific mutations and transcriptional networks. Genes regulating methylation was not associated with significant change in TE expression. Using an unpenalized cox regression analysis we identified a TE expression signature that predicted prognosis in AML. We identified 14 candidate prognostic TE transcripts (TEP) that classified AML as high/low-risk and this was independent of mutation-based and coding-gene expression based risk-stratification. TEP was able to predict prognosis in independent cohorts of 284 pediatric AML patients and 19 relapsed adult AML patients. This first comprehensive study of TE expression in AML demonstrates that TE expression can serve as a biomarker for prognosis in AML, and provides novel insights into the biology of AML. Studies characterizing its role in other cancers are warranted.

Human pluripotent reprogramming with CRISPR activators

CRISPR-Cas9-based gene activation (CRISPRa) is an attractive tool for cellular reprogramming applications due to its high multiplexing capacity and direct targeting of endogenous loci. Here we present the reprogramming of primary human skin fibroblasts into induced pluripotent stem cells (iPSCs) using CRISPRa, targeting endogenous OCT4, SOX2, KLF4, MYC, and LIN28A promoters. The low basal reprogramming efficiency can be improved by an order of magnitude by additionally targeting a conserved Alu-motif enriched near genes involved in embryo genome activation (EEA-motif). This effect is mediated in part by more efficient activation of NANOG and REX1. These data demonstrate that human somatic cells can be reprogrammed into iPSCs using only CRISPRa. Furthermore, the results unravel the involvement of EEA-motif-associated mechanisms in cellular reprogramming.

Mouse Obox and Crxos modulate preimplantation transcriptional profiles revealing similarity between paralogous mouse and human homeobox genes

BACKGROUND

ETCHbox genes are eutherian-specific homeobox genes expressed during preimplantation development at a time when the first cell lineage decisions are being made. The mouse has an unusual repertoire of ETCHbox genes with several gene families lost in evolution and the remaining two, Crxos and Obox, greatly divergent in sequence and number. Each has undergone duplication to give a double homeodomain Crxos locus and a large cluster of over 60 Obox loci. The gene content differences between species raise important questions about how evolution can tolerate loss of genes implicated in key developmental events.

RESULTS

We find that Crxos internal duplication occurred in the mouse lineage, while Obox duplication was stepwise, generating subgroups with distinct sequence and expression. Ectopic expression of three Obox genes and a Crxos transcript in primary mouse embryonic cells followed by transcriptome sequencing allowed investigation into their functional roles. We find distinct transcriptomic influences for different Obox subgroups and Crxos, including modulation of genes related to zygotic genome activation and preparation for blastocyst formation. Comparison with similar experiments performed using human homeobox genes reveals striking overlap between genes downstream of mouse Crxos and genes downstream of human ARGFX.

CONCLUSIONS

Mouse Crxos and human ARGFX homeobox genes are paralogous rather than orthologous, yet they have evolved to regulate a common set of genes. This suggests there was compensation of function alongside gene loss through co-option of a different locus. Functional compensation by non-orthologous genes with dissimilar sequences is unusual but may indicate underlying distributed robustness. Compensation may be driven by the strong evolutionary pressure for successful early embryo development.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Suppression of Transposable Elements in Leukemic Stem Cells

Genomic transposable elements (TEs) comprise nearly half of the human genome. The expression of TEs is considered potentially hazardous, as it can lead to insertional mutagenesis and genomic instability. However, recent studies have revealed that TEs are involved in immune-mediated cell clearance. Hypomethylating agents can increase the expression of TEs in cancer cells, inducing ‘viral mimicry’, causing interferon signalling and cancer cell killing. To investigate the role of TEs in the pathogenesis of acute myeloid leukaemia (AML), we studied TE expression in several cell fractions of AML while tracking its development (pre-leukemic haematopoietic stem cells, leukemic stem cells [LSCs], and leukemic blasts). LSCs, which are resistant to chemotherapy and serve as reservoirs for relapse, showed significant suppression of TEs and interferon pathways. Similarly, high-risk cases of myelodysplastic syndrome (MDS) showed far greater suppression of TEs than low-risk cases. We propose TE suppression as a mechanism for immune escape in AML and MDS. Repression of TEs co-occurred with the upregulation of several genes known to modulate TE expression, such as RNA helicases and autophagy genes. Thus, we have identified potential pathways that can be targeted to activate cancer immunogenicity via TEs in AML and MDS.