E2f6 and Bmi1 cooperate in axial skeletal development

The osteoporosis susceptibility SNP rs188303909 at 2q14.2 regulates EN1 expression by modulating DNA methylation and E2F6 binding

EN1 encodes a homeodomain-containing transcription factor and is a determinant of bone density and fracture. Previous powerful genome-wide association studies (GWASs) have identified multiple single-nucleotide polymorphisms (SNPs) near EN1 at 2q14.2 locus for osteoporosis, but the causal SNPs and functional mechanisms underlying these associations are poorly understood. The target genes regulated by the transcription factor EN1 are also unclear. In this study, we identified rs188303909, a functional CpG-SNP, as a causal SNP for osteoporosis at 2q14.2 through the integration of functional and epigenomic analyses. Functional experiments demonstrated that unmethylated rs188303909 acted as a strong allele-specific distal enhancer to regulate EN1 expression by modifying the binding of transcription factor E2F6, but rs188303909 methylation attenuated the active effect of E2F6 on EN1 expression. Importantly, transcription factor EN1 could differentially bind osteoporosis GWAS lead SNPs rs4869739-T and rs4355801-G to upregulate CCDC170 and COLEC10 expression, thus promoting bone formation. Our study provided a mechanistic insight into expression regulation of the osteoporosis susceptibility gene EN1, which could be a potential therapeutic target for osteoporosis precision medicine. KEY MESSAGES: CpG-SNP rs188303909 is a causal SNP at the osteoporosis susceptibility locus 2q14.2. Rs188303909 distally regulates EN1 expression by modulating DNA methylation and E2F6 binding. EN1 upregulates CCDC170 and COLEC10 expression through osteoporosis GWAS lead SNPs rs4869739 and rs4355801.

Angular limb deformity associated with TSPAN18, NRG3 and NOVA2 in Rambouillet rams

Angular limb deformity (ALD) affects many species of livestock and companion animals. The mechanisms of ALD development are not well understood, but previous research suggests the involvement of genetic risk factors. A case-control genome-wide association study (GWAS) was conducted with 40 ALD-affected and 302 unaffected Rambouillet rams and 40,945 single nucleotide polymorphisms (SNPs). Forelimbs of 6 ALD-affected rams were examined and diagnosed with osteochondrosis. Genome-wide or chromosome-wide significant SNPs were positioned exonic, intronic or within the 3'UTR of genes TSPAN18, NRG3 and NOVA2, respectively. These genes have previously described roles related to angiogenesis and osteoblast, osteoclast and chondrocyte proliferation and differentiation, which suggests the possibility for their involvement in the pathogenesis of osteochondrosis. Functional consequences of SNPs were evaluated through transcription factor binding site analysis, which predicted binding sites for transcription factors of known importance to bone growth, including SOX6, SOX9 and RUNX2. The identification of genetic risk factors for ALD may help to improve animal welfare and production in Rambouillet, a breed known to be at risk for ALD development. This study proposes genes TSPAN18, NRG3 and NOVA2 as targets for further research towards understanding the etiology of ALD in Rambouillet sheep.

An association study on imputed whole-genome resequencing from high-throughput sequencing data for body traits in crossbred pigs

Body traits are important economic factors in the pig industry. Genome-wide association studies (GWASs) have been widely applied using high-density genotype data to detect QTL in pigs. The aim of the present study was to detect the genetic variants significantly associated with body traits in crossbred pigs using the Illumina Porcine SNP50 BeadChip and imputed whole-genome sequence data. A set of seven body traits - body length, body height, chest circumference, cannon bone circumference, leg buttock circumference, back fat thickness and loin muscle depth - were measured. Moderate to high heritabilities were obtained for most traits (from 0.14 to 0.46), and significant genetic and phenotypic correlations among them were observed. GWAS identified 714 significantly associated SNPs located at 39 regions on all autosomes for body traits, and a total of seven functionally related candidate genes: PIK3CD, HOXA, PCGF2, CHST11, COL2A1, BMI1 and OSR2. Functional enrichment analysis revealed that candidate genes were enriched in the estrogen signaling pathway, embryonic skeletal system morphogenesis and embryonic skeletal system development. These results aim to uncover the genetic mechanisms underlying body development and marker-assisted selection programs focusing on body traits in pigs.

E2F6 initiates stable epigenetic silencing of germline genes during embryonic development

In mouse development, long-term silencing by CpG island DNA methylation is specifically targeted to germline genes; however, the molecular mechanisms of this specificity remain unclear. Here, we demonstrate that the transcription factor E2F6, a member of the polycomb repressive complex 1.6 (PRC1.6), is critical to target and initiate epigenetic silencing at germline genes in early embryogenesis. Genome-wide, E2F6 binds preferentially to CpG islands in embryonic cells. E2F6 cooperates with MGA to silence a subgroup of germline genes in mouse embryonic stem cells and in embryos, a function that critically depends on the E2F6 marked box domain. Inactivation of E2f6 leads to a failure to deposit CpG island DNA methylation at these genes during implantation. Furthermore, E2F6 is required to initiate epigenetic silencing in early embryonic cells but becomes dispensable for the maintenance in differentiated cells. Our findings elucidate the mechanisms of epigenetic targeting of germline genes and provide a paradigm for how transient repression signals by DNA-binding factors in early embryonic cells are translated into long-term epigenetic silencing during mouse development.

DNA methylation targets CpG island promoters of germline genes to repress their expression in mouse somatic cells. Here the authors show that a transcription factor E2F6 is required to target CpG island DNA methylation and epigenetic silencing to germline genes during early mouse development.

Lin28a/let-7 pathway modulates the Hox code via Polycomb regulation during axial patterning in vertebrates

The body plan along the anteroposterior axis and regional identities are specified by the spatiotemporal expression of Hox genes. Multistep controls are required for their unique expression patterns; however, the molecular mechanisms behind the tight control of Hox genes are not fully understood. In this study, we demonstrated that the Lin28a/let-7 pathway is critical for axial elongation. Lin28a^–/– mice exhibited axial shortening with mild skeletal transformations of vertebrae, which were consistent with results in mice with tail bud-specific mutants of Lin28a. The accumulation of let-7 in Lin28a^–/– mice resulted in the reduction of PRC1 occupancy at the Hox cluster loci by targeting Cbx2. Consistently, Lin28a loss in embryonic stem-like cells led to aberrant induction of posterior Hox genes, which was rescued by the knockdown of let-7. These results suggest that the Lin28/let-7 pathway is involved in the modulation of the ‘Hox code’ via Polycomb regulation during axial patterning.

Variation in zygotic CRISPR/Cas9 gene editing outcomes generates novel reporter and deletion alleles at the Gdf11 locus

Recent advances in CRISPR/Cas gene editing technology have significantly expanded the possibilities and accelerated the pace of creating genetically engineered animal models. However, CRISPR/Cas-based strategies designed to precisely edit the genome can often yield unintended outcomes. Here, we report the use of zygotic CRISPR/Cas9 injections to generate a knock-in GFP reporter mouse at the Gdf11 locus. Phenotypic and genomic characterization of founder animals from these injections revealed a subset that contained the correct targeting event and exhibited GFP expression that, within the hematopoietic system, was restricted predominantly to lymphoid cells. Yet, in another subset of founder mice, we detected aberrant integration events at the target site that dramatically and inaccurately shifted hematopoietic GFP expression from the lymphoid to the myeloid lineage. Additionally, we recovered multiple Gdf11 deletion alleles that modified the C-terminus of the GDF11 protein. When bred to homozygosity, most of these alleles recapitulated skeletal phenotypes reported previously for Gdf11 knockout mice, suggesting that these represent null alleles. However, we also recovered one Gdf11 deletion allele that encodes a novel GDF11 variant protein ("GDF11-WE") predicted to contain two additional amino acids (tryptophan (W) and glutamic acid (E)) at the C-terminus of the mature ligand. Unlike the other Gdf11 deletion alleles recovered in this study, homozygosity for the Gdf11WE allele did not phenocopy Gdf11 knockout skeletal phenotypes. Further investigation using in vivo and in vitro approaches demonstrated that GDF11-WE retains substantial physiological function, indicating that GDF11 can tolerate at least some modifications of its C-terminus and providing unexpected insights into its biochemical activities. Altogether, our study confirms that one-step zygotic injections of CRISPR/Cas gene editing complexes provide a quick and powerful tool to generate gene-modified mouse models. Moreover, our findings underscore the critical importance of thorough characterization and validation of any modified alleles generated by CRISPR, as unintended on-target effects that fail to be detected by simple PCR screening can produce substantially altered phenotypic readouts.

Bioinformatics analysis of novel transcription factors and related differentially regulated modules in non-union skeletal fractures

OBJECTIVE

This study aimed to further clarify the underlying pathomechanism of non-union skeletal fractures.

METHODS

Gene expression profile dataset GSE494 obtained from six non-union skeletal fracture and six normal samples was downloaded from the Gene Expression Omnibus database. Overlapping genes in at least two platforms were analyzed, and differentially expressed genes (DEGs) between normal and disease groups were screened. Transcriptional regulatory relationships and differentially regulated modules of various transcription factors (TFs) were determined. Differentially regulated modules with unknown functions were subjected to functional enrichment analysis.

RESULTS

Overall, 4,252 overlapping genes in at least two platforms and 77 DEGs, including 31 up and 46 downregulated genes, were obtained. Overall, 64,623 transcriptional regulatory relationships, including 49 TFs and 3,900 target genes, and 9 significant modules for differential regulation were identified. Three modules with unknown functions regulated by TFs, including zinc finger, ZZ-type containing 3 (ZZZ3), nuclear TF Y, alpha (NFYA), and POU class 2 homeobox 2 (POU2F2), were identified. Enriched GO-BP terms of NFYA and POU2F2 modules included cell adhesion and related terms and those of ZZ3 included cell cycle, cell proliferation, and associated terms.

CONCLUSION

Three TFs, including ZZZ3, POU2F2, and NFYA, and their regulated modules may have important effects on non-union skeletal fractures. Cell proliferation may be related with ZZZ3; cell adhesion and its similar process may be related with POU2F2 and NFYA.

Evolutionary changes of Hox genes and relevant regulatory factors provide novel insights into mammalian morphological modifications

The diversity of body plans of mammals accelerates the innovation of lifestyles and the extensive adaptation to different habitats, including terrestrial, aerial and aquatic habitats. However, the genetic basis of those phenotypic modifications, which have occurred during mammalian evolution, remains poorly explored. In the present study, we synthetically surveyed the evolutionary pattern of Hox clusters that played a powerful role in the morphogenesis along the head–tail axis of animal embryos and the main regulatory factors (Mll, Bmi1 and E2f6) that control the expression of Hox genes. A deflected density of repetitive elements and lineage‐specific radical mutations of Mll have been determined in marine mammals with morphological changes, suggesting that evolutionary changes may alter Hox gene expression in these lineages, leading to the morphological modification of these lineages. Although no positive selection was detected at certain ancestor nodes of lineages, the increased ω values of Hox genes implied the relaxation of functional constraints of these genes during the mammalian evolutionary process. More importantly, 49 positively‐selected sites were identified in mammalian lineages with phenotypic modifications, indicating adaptive evolution acting on Hox genes and regulatory factors. In addition, 3 parallel amino acid substitutions in some Hox genes were examined in marine mammals, which might be responsible for their streamlined body.

The Snail repressor recruits EZH2 to specific genomic sites through the enrollment of the lncRNA HOTAIR in epithelial-to-mesenchymal transition

The transcription factor Snail is a master regulator of cellular identity and epithelial-to-mesenchymal transition (EMT) directly repressing a broad repertoire of epithelial genes. How chromatin modifiers instrumental to its activity are recruited to Snail-specific binding sites is unclear. Here we report that the long non-coding RNA (lncRNA) HOTAIR (for HOX Transcript Antisense Intergenic RNA) mediates a physical interaction between Snail and enhancer of zeste homolog 2 (EZH2), an enzymatic subunit of the polycomb-repressive complex 2 and the main writer of chromatin-repressive marks. The Snail-repressive activity, here monitored on genes with a pivotal function in epithelial and hepatic morphogenesis, differentiation and cell-type identity, depends on the formation of a tripartite Snail/HOTAIR/EZH2 complex. These results demonstrate an lncRNA-mediated mechanism by which a transcriptional factor conveys a general chromatin modifier to specific genes, thereby allowing the execution of hepatocyte transdifferentiation; moreover, they highlight HOTAIR as a crucial player in the Snail-mediated EMT.

Disruption of histone methylation in developing sperm impairs offspring health transgenerationally

A father's lifetime experiences can be transmitted to his offspring to affect health and development. However, the mechanisms underlying paternal epigenetic transmission are unclear. Unlike in somatic cells, there are few nucleosomes in sperm, and their function in epigenetic inheritance is unknown. We generated transgenic mice in which overexpression of the histone H3 lysine 4 (H3K4) demethylase KDM1A (also known as LSD1) during spermatogenesis reduced H3K4 dimethylation in sperm. KDM1A overexpression in one generation severely impaired development and survivability of offspring. These defects persisted transgenerationally in the absence of KDM1A germline expression and were associated with altered RNA profiles in sperm and offspring. We show that epigenetic inheritance of aberrant development can be initiated by histone demethylase activity in developing sperm, without changes to DNA methylation at CpG-rich regions.

Targeted Deletion of Btg1 and Btg2 Results in Homeotic Transformation of the Axial Skeleton

Btg1 and Btg2 encode highly homologous proteins that are broadly expressed in different cell lineages, and have been implicated in different types of cancer. Btg1 and Btg2 have been shown to modulate the function of different transcriptional regulators, including Hox and Smad transcription factors. In this study, we examined the in vivo role of the mouse Btg1 and Btg2 genes in specifying the regional identity of the axial skeleton. Therefore, we examined the phenotype of Btg1 and Btg2 single knockout mice, as well as novel generated Btg1^-/-;Btg2^-/- double knockout mice, which were viable, but displayed a non-mendelian inheritance and smaller litter size. We observed both unique and overlapping phenotypes reminiscent of homeotic transformation along the anterior-posterior axis in the single and combined Btg1 and Btg2 knockout animals. Both Btg1^-/- and Btg2^-/- mice displayed partial posterior transformation of the seventh cervical vertebra, which was more pronounced in Btg1^-/-;Btg2^-/- mice, demonstrating that Btg1 and Btg2 act in synergy. Loss of Btg2, but not Btg1, was sufficient for complete posterior transformation of the thirteenth thoracic vertebra to the first lumbar vertebra. Moreover, Btg2^-/- animals displayed complete posterior transformation of the sixth lumbar vertebra to the first sacral vertebra, which was only partially present at a low frequency in Btg1^-/- mice. The Btg1^-/-;Btg2^-/- animals showed an even stronger phenotype, with L5 to S1 transformation. Together, these data show that both Btg1 and Btg2 are required for normal vertebral patterning of the axial skeleton, but each gene contributes differently in specifying the identity along the anterior-posterior axis of the skeleton.

E2f6-mediated repression of the meiotic Stag3 and Smc1β genes during early embryonic development requires Ezh2 and not the de novo methyltransferase Dnmt3b

The E2f6 transcriptional repressor is an E2F-family member essential for the silencing of a group of meiosis-specific genes in somatic tissues. Although E2f6 has been shown to associate with both polycomb repressive complexes (PRC) and the methyltransferase Dnmt3b, the cross-talk between these repressive machineries during E2f6-mediated gene silencing has not been clearly demonstrated yet. In particular, it remains largely undetermined when and how E2f6 establishes repression of meiotic genes during embryonic development. We demonstrate here that the inactivation of a group of E2f6 targeted genes, including Stag3 and Smc1β, first occurs at the transition from mouse embryonic stem cells (ESCs) to epiblast stem cells (EpiSCs), which represent pre- and post-implantation stages, respectively. This process was accompanied by de novo methylation of their promoters. Of interest, despite a clear difference in DNA methylation status, E2f6 was similarly bound to the proximal promoter regions both in ESCs and EpiSCs. Neither E2f6 nor Dnmt3b overexpression in ESCs decreased meiotic gene expression or increased DNA methylation, indicating that additional factors are required for E2f6-mediated repression during the transition. When the SET domain of Ezh2, a core subunit of the PRC2 complex, was deleted, however, repression of Stag3 and Smc1β during embryoid body differentiation was largely impaired, indicating that the event required the enzymatic activity of Ezh2. In addition, repression of Stag3 and Smc1β occurred in the absence of Dnmt3b. The data presented here suggest a primary role of PRC2 in E2f6-mediated gene silencing of the meiotic genes.

In vivo regulation of E2F1 by Polycomb group genes in Drosophila

The E2F transcription factors are important regulators of the cell cycle whose function is commonly misregulated in cancer. To identify novel regulators of E2F1 activity in vivo, we used Drosophila to conduct genetic screens. For this, we generated transgenic lines that allow the tissue-specific depletion of dE2F1 by RNAi. Expression of these transgenes using Gal4 drivers in the eyes and wings generated reliable and modifiable phenotypes. We then conducted genetic screens testing the capacity of Exelixis deficiencies to modify these E2F1-RNAi phenotypes. From these screens, we identified mutant alleles of Suppressor of zeste 2 [Su(z)2] and multiple Polycomb group genes as strong suppressors of the E2F1-RNA interference phenotypes. In validation of our genetic data, we find that depleting Su(z)2 in cultured Drosophila cells restores the cell-proliferation defects caused by reduction of dE2F1 by elevating the level of dE2f1. Furthermore, analyses of methylation status of histone H3 lysine 27 (H3K27me) from the published modENCODE data sets suggest that the genomic regions harboring dE2f1 gene and certain dE2f1 target genes display H3K27me during development and in several Drosophila cell lines. These in vivo observations suggest that the Polycomb group may regulate cell proliferation by repressing the transcription of dE2f1 and certain dE2F1 target genes. This mechanism may play an important role in coordinating cellular differentiation and proliferation during Drosophila development.

BMI1 as a novel target for drug discovery in cancer

Growing evidence has demonstrated that clonogenic cancer stem (initiating) cells are responsible for tumor regrowth and disease relapse. Bmi-1 plays a critical role in the self-renewal of adult stem cells. The Bmi-1 protein is elevated in many types of cancers, and experimental reduction of Bmi-1 protein levels by small interfering RNA (siRNA) causes apoptosis and/or senescence in tumor cells in vitro and increases susceptibility to cytotoxic agents. The Bmi-1 protein has no known enzymatic activity, but serves as the key regulatory component of the PRC1 complex (polycomb repressive complex-1). This complex influences chromatin structure and regulates transcriptional activity of a number of important loci including the Ink4a locus which encodes the tumor suppressor proteins p16(Ink4a) and p14(Arf) . In this prospective study, we will discuss the implication of BMI1 in cancers, the biology of BMI1, and the regulatory control of BMI1 expression. The target validation and the future prospects of targeting BMI1 in cancer therapy are also discussed.

G1 arrest and differentiation can occur independently of Rb family function

Repression of E2F target genes is required for cell cycle arrest in Rb family (Rb, p107, and p130)-deficient cells.

The ability of progenitor cells to exit the cell cycle is essential for proper embryonic development and homeostasis, but the mechanisms governing cell cycle exit are still not fully understood. Here, we tested the requirement for the retinoblastoma (Rb) protein and its family members p107 and p130 in G0/G1 arrest and differentiation in mammalian cells. We found that Rb family triple knockout (TKO) mouse embryos survive until days 9–11 of gestation. Strikingly, some TKO cells, including in epithelial and neural lineages, are able to exit the cell cycle in G0/G1 and differentiate in teratomas and in culture. This ability of TKO cells to arrest in G0/G1 is associated with the repression of key E2F target genes. Thus, G1 arrest is not always dependent on Rb family members, which illustrates the robustness of cell cycle regulatory networks during differentiation and allows for the identification of candidate pathways to inhibit the expansion of cancer cells with mutations in the Rb pathway.

Dnmt3b recruitment through E2F6 transcriptional repressor mediates germ-line gene silencing in murine somatic tissues

Methylation of cytosine residues within the CpG dinucleotide in mammalian cells is an important mediator of gene expression, genome stability, X-chromosome inactivation, genomic imprinting, chromatin structure, and embryonic development. The majority of CpG sites in mammalian cells is methylated in a nonrandom fashion, raising the question of how DNA methylation is distributed along the genome. Here, we focused on the functions of DNA methyltransferase-3b (Dnmt3b), of which deregulated activity is linked to several human pathologies. We generated Dnmt3b hypomorphic mutant mice with reduced catalytic activity, which first revealed a deregulation of Hox genes expression, consistent with the observed homeotic transformations of the posterior axis. In addition, analysis of deregulated expression programs in Dnmt3b mutant embryos, using DNA microarrays, highlighted illegitimate activation of several germ-line genes in somatic tissues that appeared to be linked directly to their hypomethylation in mutant embryos. We provide evidence that these genes are direct targets of Dnmt3b. Moreover, the recruitment of Dnmt3b to their proximal promoter is dependant on the binding of the E2F6 transcriptional repressor, which emerges as a common hallmark in the promoters of genes found to be up-regulated as a consequence of impaired Dnmt3b activity. Therefore, our results unraveled a coordinated regulation of genes involved in meiosis, through E2F6-dependant methylation and transcriptional silencing in somatic tissues.

Cell-context specific role of the E2F/Rb pathway in development and disease

Development of the central nervous system (CNS) requires the generation of neuronal and glial cell subtypes in appropriate numbers, and this demands the careful coordination of cell-cycle exit, survival, and differentiation. The E2F/Rb pathway is critical for cell-cycle regulation and also modulates survival and differentiation of distinct cell types in the developing and adult CNS. In this review, we first present the specific temporal patterns of expression of the E2F and Rb family members during CNS development and then discuss the genetic ablation of single or multiple members of these two families. Overall, the available data suggest a time-dependent and cell-context specific role of E2F and Rb family members in the developing and adult CNS.

Polycomb group complexes--many combinations, many functions

Polycomb Group (PcG) proteins are transcription regulatory proteins that control the expression of a variety of genes from early embryogenesis through birth to adulthood. PcG proteins form several complexes that are thought to collaborate to repress gene transcription. Individual PcG proteins have unique characteristics, and mutations in genes encoding different PcG proteins cause distinct phenotypes. Histone modifications have important roles in some PcG protein functions, but they are not universally required. The mechanisms of gene-specific recruitment, transcription repression, and selective derepression of genes by vertebrate PcG proteins are incompletely understood. Future studies of this enigmatic group of developmental regulators are certain to produce unanticipated discoveries.

E2f3a and E2f3b make overlapping but different contributions to total E2f3 activity

The E2f transcription factors are key downstream targets of the retinoblastoma protein tumor suppressor that control cell proliferation. E2F3 has garnered particular attention because it is amplified in various human tumors. E2f3 mutant mice typically die around birth and E2f3-deficient cells have a proliferation defect that correlates with impaired E2f target gene activation and also induction of p19(Arf) and p53. The E2f3 locus encodes two isoforms, E2f3a and E2f3b, which differ in their N-termini. However, it is unclear how E2f3a versus E2f3b contributes to E2f3's requirement in either proliferation or development. To address this, we use E2f3a- and E2f3b-specific knockouts. We show that inactivation of E2f3a results in a low penetrance proliferation defect in vitro whereas loss of E2f3b has no effect. This proliferation defect appears insufficient to disrupt normal development as E2f3a and E2f3b mutant mice are both fully viable and have no detectable defects. However, when combined with E2f1 mutation, inactivation of E2f3a, but not E2f3b, causes significant proliferation defects in vitro, neonatal lethality and also a striking cartilage defect. Thus, we conclude that E2f3a and E2f3b have largely overlapping functions in vivo and that E2f3a can fully substitute for E2f1 and E2f3 in most murine tissues.