Interaction network of human early embryonic transcription factors.

Embryonic genome activation (EGA) occurs during preimplantation development and is characterized by the initiation of de novo transcription from the embryonic genome. Despite its importance, the regulation of EGA and the transcription factors involved in this process are poorly understood. Paired-like homeobox (PRDL) family proteins are implicated as potential transcriptional regulators of EGA, yet the PRDL-mediated gene regulatory networks remain uncharacterized. To investigate the function of PRDL proteins, we are identifying the molecular interactions and the functions of a subset family of the Eutherian Totipotent Cell Homeobox (ETCHbox) proteins, seven PRDL family proteins and six other transcription factors (TFs), all suggested to participate in transcriptional regulation during preimplantation. Using mass spectrometry-based interactomics methods, AP-MS and proximity-dependent biotin labeling, and chromatin immunoprecipitation sequencing we derive the comprehensive regulatory networks of these preimplantation TFs. By these interactomics tools we identify more than a thousand high-confidence interactions for the 21 studied bait proteins with more than 300 interacting proteins. We also establish that TPRX2, currently assigned as pseudogene, is a transcriptional activator.

Embryonic LTR retrotransposons supply promoter modules to somatic tissues.

Long terminal repeat (LTR) retrotransposons are widely distributed across the human genome. They have accumulated through retroviral integration into germline DNA and are latent genetic modules. Active LTR promoters are observed in germline cells; however, little is known about the mechanisms underlying their active transcription in somatic tissues. Here, by integrating our previous transcriptome data set with publicly available data sets, we show that the LTR families MLT2A1 and MLT2A2 are primarily expressed in human four-cell and eight-cell embryos and are also activated in some adult somatic tissues, particularly pineal gland. Three MLT2A elements function as the promoters and first exons of the protein-coding genes ABCE1, COL5A1, and GALNT13 specifically in the pineal gland of humans but not in that of macaques, suggesting that the exaptation of these LTRs as promoters occurred during recent primate evolution. This analysis provides insight into the possible transition from germline insertion to somatic expression of LTR retrotransposons.

The human PRD-like homeobox gene LEUTX has a central role in embryo genome activation.

Leucine twenty homeobox (LEUTX) is a paired (PRD)-like homeobox gene that is expressed almost exclusively in human embryos during preimplantation development. We previously identified a novel transcription start site for the predicted human LEUTX gene based on the transcriptional analysis of human preimplantation embryos. The novel variant encodes a protein with a complete homeodomain. Here, we provide a detailed description of the molecular cloning of the complete homeodomain-containing LEUTX Using a human embryonic stem cell overexpression model we show that the complete homeodomain isoform is functional and sufficient to activate the transcription of a large proportion of the genes that are upregulated in human embryo genome activation (EGA), whereas the previously predicted partial homeodomain isoform is largely inactive. Another PRD-like transcription factor, DPRX, is then upregulated as a powerful repressor of transcription. We propose a two-stage model of human EGA in which LEUTX acts as a transcriptional activator at the 4-cell stage, and DPRX as a balancing repressor at the 8-cell stage. We conclude that LEUTX is a candidate regulator of human EGA.

Characterization of the human RFX transcription factor family by regulatory and target gene analysis.

BACKGROUND: Evolutionarily conserved RFX transcription factors (TFs) regulate their target genes through a DNA sequence motif called the X-box. Thereby they regulate cellular specialization and terminal differentiation. Here, we provide a comprehensive analysis of all the eight human RFX genes (RFX1-8), their spatial and temporal expression profiles, potential upstream regulators and target genes. RESULTS: We extracted all known human RFX1-8 gene expression profiles from the FANTOM5 database derived from transcription start site (TSS) activity as captured by Cap Analysis of Gene Expression (CAGE) technology. RFX genes are broadly (RFX1-3, RFX5, RFX7) and specifically (RFX4, RFX6) expressed in different cell types, with high expression in four organ systems: immune system, gastrointestinal tract, reproductive system and nervous system. Tissue type specific expression profiles link defined RFX family members with the target gene batteries they regulate. We experimentally confirmed novel TSS locations and characterized the previously undescribed RFX8 to be lowly expressed. RFX tissue and cell type specificity arises mainly from differences in TSS architecture. RFX transcript isoforms lacking a DNA binding domain (DBD) open up new possibilities for combinatorial target gene regulation. Our results favor a new grouping of the RFX family based on protein domain composition. We uncovered and experimentally confirmed the TFs SP2 and ESR1 as upstream regulators of specific RFX genes. Using TF binding profiles from the JASPAR database, we determined relevant patterns of X-box motif positioning with respect to gene TSS locations of human RFX target genes. CONCLUSIONS: The wealth of data we provide will serve as the basis for precisely determining the roles RFX TFs play in human development and disease.

An approach to comprehensive genome and proteome expression analyses in Schwann cells and neurons during peripheral nerve myelin formation.

Peripheral nerve myelination is a complex event resulting from spatially and temporally regulated reciprocal interactions between the neuron and myelin-forming Schwann cells. The dynamic process and the protein functional modules and networks that operate throughout the myelination process are poorly understood because of a lack of methodologies suitable for observing specific changes in the Schwann cell/neuron-unit. The identification of the precise roles for the proteins participating in the functional modules and networks that participate in the myelination process is hindered by the cellular and molecular complexity of the nervous tissue itself. We have developed an approach based on a myelinating dorsal root ganglion explant model that allows distinguishing clear, reproducible and predictable differences between the biochemical properties and the genomic and proteomic expression profiles of both cellular components of the Schwann cell/neuron unit at different stages of the myelination process. This model, derived from E13.5 C57BL/6J mouse embryos, is sufficiently robust for use in identifying the protein functional networks and modules related to peripheral nerve myelin formation. The genomic expression profiles of the selected neuronal, Schwann cell and myelin-specific proteins in the cultures reflect in vivo profiles reported in the literature, and the structural and ultrastructural properties of the myelin, as well as the myelination schedule of the cultures, closely resemble those observed in peripheral nerves in situ. The RNA expression data set is available through NCBI gene expression omnibus accession GSE60345. We have developed a reproducible and robust cell culture-based approach, accompanied by a genome-wide expression data set, which allows studying myelination in the peripheral nervous system at the proteomic and transcriptomic levels in Schwann cells and neurons. Myelinating dorsal root explant cultures, prepared from C57BL/6J mouse embryos, present distinct developmental stages comparable to those observed in a peripheral nerve in situ. This model can be used for identifying the protein functional networks and modules related to peripheral nerve myelin formation.

Gene expression analysis of skin grafts and cultured keratinocytes using synthetic RNA normalization reveals insights into differentiation and growth control.

BACKGROUND: Keratinocytes (KCs) are the most frequent cells in the epidermis, and they are often isolated and cultured in vitro to study the molecular biology of the skin. Cultured primary cells and various immortalized cells have been frequently used as skin models but their comparability to intact skin has been questioned. Moreover, when analyzing KC transcriptomes, fluctuation of polyA+ RNA content during the KCs' lifecycle has been omitted. RESULTS: We performed STRT RNA sequencing on 10 ng samples of total RNA from three different sample types: i) epidermal tissue (split-thickness skin grafts), ii) cultured primary KCs, and iii) HaCaT cell line. We observed significant variation in cellular polyA+ RNA content between tissue and cell culture samples of KCs. The use of synthetic RNAs and SAMstrt in normalization enabled comparison of gene expression levels in the highly heterogenous samples and facilitated discovery of differences between the tissue samples and cultured cells. The transcriptome analysis sensitively revealed genes involved in KC differentiation in skin grafts and cell cycle regulation related genes in cultured KCs and emphasized the fluctuation of transcription factors and non-coding RNAs associated to sample types. CONCLUSIONS: The epidermal keratinocytes derived from tissue and cell culture samples showed highly different polyA+ RNA contents. The use of SAMstrt and synthetic RNA based normalization allowed the comparison between tissue and cell culture samples and thus proved to be valuable tools for RNA-seq analysis with translational approach. Transciptomics revealed clear difference both between tissue and cell culture samples and between primary KCs and immortalized HaCaT cells.

Neurofibromatosis type 1 gene mutation analysis using sequence capture and high-throughput sequencing.

Neurofibromatosis type 1 syndrome (NF1) is caused by mutations in the NF1 gene. Availability of new sequencing technology prompted us to search for an alternative method for NF1 mutation analysis. Genomic DNA was isolated from saliva avoiding invasive sampling. The NF1 exons with an additional 50bp of flanking intronic sequences were captured and enriched using the SeqCap EZ Choice Library protocol. The captured DNA was sequenced with the Roche/454 GS Junior system. The mean coverages of the targeted regions were 41x and 74x in 2 separate sets of samples. An NF1 mutation was discovered in 10 out of 16 separate patient samples. Our study provides proof of principle that the sequence capture methodology combined with high-throughput sequencing is applicable to NF1 mutation analysis. Deep intronic mutations may however remain undetectable, and change at the DNA level may not predict the outcome at the mRNA or protein levels.

Comprehensive characterization of the embryonic factor LEUTX.

The paired-like homeobox transcription factor LEUTX is expressed in human preimplantation embryos between the 4- and 8-cell stages, and then silenced in somatic tissues. To characterize the function of LEUTX, we performed a multiomic characterization of LEUTX using two proteomics methods and three genome-wide sequencing approaches. Our results show that LEUTX stably interacts with the EP300 and CBP histone acetyltransferases through its 9 amino acid transactivation domain (9aaTAD), as mutation of this domain abolishes the interactions. LEUTX targets genomic cis-regulatory sequences that overlap with repetitive elements, and through these elements it is suggested to regulate the expression of its downstream genes. We find LEUTX to be a transcriptional activator, upregulating several genes linked to preimplantation development as well as 8-cell-like markers, such as DPPA3 and ZNF280A. Our results support a role for LEUTX in preimplantation development as an enhancer binding protein and as a potent transcriptional activator.

The pathoetiology of neurofibromatosis 1.

Although a mutation in the NF1 gene is the only factor required to initiate the neurocutaneous-skeletal neurofibromatosis 1 (NF1) syndrome, the pathoetiology of the multiple manifestations of this disease in different organ systems seems increasingly complex. The wide spectrum of different clinical phenotypes and their development, severity, and prognosis seem to result from the cross talk between numerous cell types, cell signaling networks, and cell-extracellular matrix interactions. The bi-allelic inactivation of the NF1 gene through a "second hit" seems to be of crucial importance to the development of certain manifestations, such as neurofibromas, café-au-lait macules, and glomus tumors. In each case, the second hit involves only one cell type, which is subsequently clonally expanded in a discrete lesion. Neurofibromas, which are emphasized in this review, and cutaneous neurofibromas in particular, are known to contain a subpopulation of NF1-diploinsufficient Schwann cells and a variety of NF1-haploinsufficient cell types. A recent study identified a multipotent precursor cell population with an NF1(+/-) genotype that resides in human cutaneous neurofibromas and that has been suggested to play a role in their pathogenesis.

The development of cutaneous neurofibromas.

Cutaneous neurofibromas are the hallmarks of neurofibromatosis type 1 (NF1). They are composed of multiple cell types, and traditionally they are believed to arise from small nerve tributaries of the skin. A key finding in the context of this view has been that subpopulations of tumor Schwann cells harbor biallelic inactivation of the NF1 gene (NF1(-/-)). In the present study, our aim was to clarify further the pathogenesis of cutaneous neurofibromas. First, we detected cells expressing multipotency-associated biomarkers in cutaneous neurofibromas. Second, we developed a method for isolating and expanding multipotent neurofibroma-derived precursor cells (NFPs) from dissociated human cutaneous neurofibromas and used it to analyze their growth and differentiation potential. In analogy to solitary cells resident in neurofibromas, NFPs were found to express nestin and had the potential to differentiate to, at least, Schwann cells, neurons, epithelial cells, and adipocytes. Mutation analysis of the NFPs revealed that their genotype was NF1(+/-). The results led us to speculate that the development of cutaneous neurofibromas includes the recruitment of multipotent NF1(+/-) precursor cells. These cells may be derived from the multipotent cells of the hair roots, which often are intimately associated with microscopic neurofibromas.

Oral soft tissue alterations in patients with neurofibromatosis.

Our aim was to characterize the type and frequency of oral soft tissue alterations in neurofibromatosis. A total of 103 patients with neurofibromatosis 1 (NF1) and three patients with neurofibromatosis 2 (NF2) were clinically evaluated for their oral soft tissue alterations. Disturbing growths were removed from nine patients with NF1 and from one patient with NF2. The specimens were analyzed using routine histological methods and with immunohistochemistry using antibodies to S100, type IV collagen, CD34, neurofilament, and neuron-specific tubulin (TUBB3). Alterations including oral tumors, overgrowths of gingival soft tissue, and enlarged papillae of the tongue were discovered in 74% of NF1 patients. The results showed that three tumors clinically classified as plexiform neurofibromas and five out of six discrete mucosal tumors displayed histology and immunohistology consistent with that of neurofibroma. The histology of one palatal lesion resembled that of a scar, and the lesion removed from the patient with NF2 was classified as an amyloid tumor. To conclude, oral soft tissue growths are common findings in NF1, but most lesions do not require treatment and the patients may even not be aware of these alterations. Collagen IV, S100, and CD34 are useful biomarkers in the analysis of NF1-related oral soft tissue tumors. The clinicians should recognize that oral soft tissue alterations are relatively common in NF1. Some of the growths are disturbing, and plexiform neurofibromas may bear a risk of malignant transformation.

CRISPR activation enables high-fidelity reprogramming into human pluripotent stem cells.

Conventional reprogramming methods rely on the ectopic expression of transcription factors to reprogram somatic cells into induced pluripotent stem cells (iPSCs). The forced expression of transcription factors may lead to off-target gene activation and heterogeneous reprogramming, resulting in the emergence of alternative cell types and aberrant iPSCs. Activation of endogenous pluripotency factors by CRISPR activation (CRISPRa) can reduce this heterogeneity. Here, we describe a high-efficiency reprogramming of human somatic cells into iPSCs using optimized CRISPRa. Efficient reprogramming was dependent on the additional targeting of the embryo genome activation-enriched Alu-motif and the miR-302/367 locus. Single-cell transcriptome analysis revealed that the optimized CRISPRa reprogrammed cells more directly and specifically into the pluripotent state when compared to the conventional reprogramming method. These findings support the use of CRISPRa for high-quality pluripotent reprogramming of human cells.

Transient DUX4 expression in human embryonic stem cells induces blastomere-like expression program that is marked by SLC34A2.

Embryonic genome activation (EGA) is critical for embryonic development. However, our understanding of the regulatory mechanisms of human EGA is still incomplete. Human embryonic stem cells (hESCs) are an established model for studying developmental processes, but they resemble epiblast and are sub-optimal for modeling EGA. DUX4 regulates human EGA by inducing cleavage-stage-specific genes, while it also induces cell death. We report here that a short-pulsed expression of DUX4 in primed hESCs activates an EGA-like gene expression program in up to 17% of the cells, retaining cell viability. These DUX4-induced cells resembled eight-cell stage blastomeres and were named induced blastomere-like (iBM) cells. The iBM cells showed marked reduction of POU5F1 protein, as previously observed in mouse two-cell-like cells. Finally, the iBM cells were successfully enriched using an antibody against NaPi2b (SLC34A2), which is expressed in human blastomeres. The iBM cells provide an improved model system to study human EGA transcriptome.

DUX4 is a multifunctional factor priming human embryonic genome activation.

Double homeobox 4 (DUX4) is expressed at the early pre-implantation stage in human embryos. Here we show that induced human DUX4 expression substantially alters the chromatin accessibility of non-coding DNA and activates thousands of newly identified transcribed enhancer-like regions, preferentially located within ERVL-MaLR repeat elements. CRISPR activation of transcribed enhancers by C-terminal DUX4 motifs results in the increased expression of target embryonic genome activation (EGA) genes ZSCAN4 and KHDC1P1. We show that DUX4 is markedly enriched in human zygotes, followed by intense nuclear DUX4 localization preceding and coinciding with minor EGA. DUX4 knockdown in human zygotes led to changes in the EGA transcriptome but did not terminate the embryos. We also show that the DUX4 protein interacts with the Mediator complex via the C-terminal KIX binding motif. Our findings contribute to the understanding of DUX4 as a regulator of the non-coding genome.

Class III beta-tubulin is a component of the mitotic spindle in multiple cell types.

The findings of this study show that Class III beta-tubulin is a component of the mitotic spindle in multiple cell types. Class III beta-tubulin has been widely used as a neuron-specific marker, but it has been detected also in association with breast and pancreatic cancers. In this study, we describe a novel finding of Class III beta-tubulin in a subpopulation of cells in malignant peripheral nerve sheath tumor. The findings of this study also show that Class III beta-tubulin is expressed by normal mesenchymal and epithelial cells (fibroblasts and keratinocytes), two transitional cell carcinoma cell lines, and neurofibroma Schwann cells, as shown by immunolabeling and Western transfer analysis using two different Tuj-1 antibodies that are specific for Class III beta-tubulin. The corresponding mRNA was detected using RT-PCR and whole human genome microarrays. Both antibodies localized Class III beta-tubulin to the mitotic spindle and showed a colocalization with alpha-tubulin. The immunoreaction became visible in early prophase, and the most intense immunoreaction was detected during metaphase and anaphase when microtubules were connected to the kinetochores on chromosomes. Class III beta-tubulin-specific immunoreaction lasted to the point when the midbody of cytokinesis became detectable.

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.

Phylogenetic and mutational analyses of human LEUTX, a homeobox gene implicated in embryogenesis.

Recently, human PAIRED-LIKE homeobox transcription factor (TF) genes were discovered whose expression is limited to the period of embryo genome activation up to the 8-cell stage. One of these TFs is LEUTX, but its importance for human embryogenesis is still subject to debate. We confirmed that human LEUTX acts as a TAATCC-targeting transcriptional activator, like other K50-type PAIRED-LIKE TFs. Phylogenetic comparisons revealed that Leutx proteins are conserved across Placentalia and comprise two conserved domains, the homeodomain, and a Leutx-specific domain containing putative transcriptional activation motifs (9aaTAD). Examination of human genotype resources revealed 116 allelic variants in LEUTX. Twenty-four variants potentially affect function, but they occur only heterozygously at low frequency. One variant affects a DNA-specificity determining residue, mutationally reachable by a one-base transition. In vitro and in silico experiments showed that this LEUTX mutation (alanine to valine at position 54 in the homeodomain) results in a transactivational loss-of-function to a minimal TAATCC-containing promoter and a 36 bp motif enriched in genes involved in embryo genome activation. A compensatory change in residue 47 restores function. The results support the notion that human LEUTX functions as a transcriptional activator important for human embryogenesis.

Characterization and target genes of nine human PRD-like homeobox domain genes expressed exclusively in early embryos.

PAIRED (PRD)-like homeobox genes belong to a class of predicted transcription factor genes. Several of these PRD-like homeobox genes have been predicted in silico from genomic sequence but until recently had no evidence of transcript expression. We found recently that nine PRD-like homeobox genes, ARGFX, CPHX1, CPHX2, DPRX, DUXA, DUXB, NOBOX, TPRX1 and TPRX2, were expressed in human preimplantation embryos. In the current study we characterized these PRD-like homeobox genes in depth and studied their functions as transcription factors. We cloned multiple transcript variants from human embryos and showed that the expression of these genes is specific to embryos and pluripotent stem cells. Overexpression of the genes in human embryonic stem cells confirmed their roles as transcription factors as either activators (CPHX1, CPHX2, ARGFX) or repressors (DPRX, DUXA, TPRX2) with distinct targets that could be explained by the amino acid sequence in homeodomain. Some PRD-like homeodomain transcription factors had high concordance of target genes and showed enrichment for both developmentally important gene sets and a 36 bp DNA recognition motif implicated in Embryo Genome Activation (EGA). Our data implicate a role for these previously uncharacterized PRD-like homeodomain proteins in the regulation of human embryo genome activation and preimplantation embryo development.

Globin mRNA reduction for whole-blood transcriptome sequencing.

The transcriptome analysis of whole-blood RNA by sequencing holds promise for the identification and tracking of biomarkers; however, the high globin mRNA (gmRNA) content of erythrocytes hampers whole-blood and buffy coat analyses. We introduce a novel gmRNA locking assay (GlobinLock, GL) as a robust and simple gmRNA reduction tool to preserve RNA quality, save time and cost. GL consists of a pair of gmRNA-specific oligonucleotides in RNA initial denaturation buffer that is effective immediately after RNA denaturation and adds only ten minutes of incubation to the whole cDNA synthesis procedure when compared to non-blood RNA analysis. We show that GL is fully effective not only for human samples but also for mouse and rat, and so far incompletely studied cow, dog and zebrafish.

Novel PRD-like homeodomain transcription factors and retrotransposon elements in early human development.

Transcriptional program that drives human preimplantation development is largely unknown. Here, by using single-cell RNA sequencing of 348 oocytes, zygotes and single blastomeres from 2- to 3-day-old embryos, we provide a detailed analysis of the human preimplantation transcriptome. By quantifying transcript far 5'-ends (TFEs), we include in our analysis transcripts that derive from alternative promoters. We show that 32 and 129 genes are transcribed during the transition from oocyte to four-cell stage and from four- to eight-cell stage, respectively. A number of identified transcripts originates from previously unannotated genes that include the PRD-like homeobox genes ARGFX, CPHX1, CPHX2, DPRX, DUXA, DUXB and LEUTX. Employing de novo promoter motif extraction on sequences surrounding TFEs, we identify significantly enriched gene regulatory motifs that often overlap with Alu elements. Our high-resolution analysis of the human transcriptome during preimplantation development may have important implications on future studies of human pluripotent stem cells and cell reprograming.