Genome-Wide Binding of Posterior HOXA/D Transcription Factors Reveals Subgrouping and Association with CTCF

HOXA10 promotes Gdf5 expression in articular chondrocytes

Growth differentiation factor 5 (GDF5), a BMP family member, is highly expressed in the surface layer of articular cartilage. The GDF5 gene is a key risk locus for osteoarthritis and Gdf5-deficient mice show abnormal joint development, indicating that GDF5 is essential in joint development and homeostasis. In this study, we aimed to identify transcription factors involved in Gdf5 expression by performing two-step screening. We first performed microarray analyses to find transcription factors specifically and highly expressed in the superficial zone (SFZ) cells of articular cartilage, and isolated 11 transcription factors highly expressed in SFZ cells but not in costal chondrocytes. To further proceed with the identification, we generated Gdf5-HiBiT knock-in (Gdf5-HiBiT KI) mice, by which we can easily and reproducibly monitor Gdf5 expression, using CRISPR/Cas9 genome editing. Among the 11 transcription factors, Hoxa10 clearly upregulated HiBiT activity in the SFZ cells isolated from Gdf5-HiBiT KI mice. Hoxa10 overexpression increased Gdf5 expression while Hoxa10 knockdown decreased it in the SFZ cells. Moreover, ChIP and promoter assays proved the direct regulation of Gdf5 expression by HOXA10. Thus, our results indicate the important role played by HOXA10 in Gdf5 regulation and the usefulness of Gdf5-HiBiT KI mice for monitoring Gdf5 expression.

CTCF variant begets to short stature by down-regulation of IGF1

Pathogenic variants in the transcription factor CCCTC-binding factor (CTCF) are associated with mental retardation, autosomal dominant 21 (MRD21, MIM#615502). Current studies supported the strong relationship between CTCF variants and growth, yet the mechanism of CTCF mutation leading to short stature is not known. Clinical information, treatment regimens, and follow-up outcomes of a patient with MRD21 were collected. The possible pathogenic mechanisms of CTCF variants leading to short stature were investigated using immortalized lymphocyte cell lines (LCLs), HEK-293T, and immortalized normal human liver cell lines (LO2). This patient received long-term treatment with recombinant human growth hormone (rhGH) which resulted in an increased height of 1.0 SDS. She had low serum insulin-like growth factor 1 (IGF1) before the treatment and the IGF1 level was not significantly increased during the treatment (-1.38 ± 0.61 SDS). The finding suggested that the CTCF R567W variant could have impaired IGF1 production pathway. We further demonstrated that the mutant CTCF had a reduced ability to bind to the promoter region of IGF1, consequently significantly reducing the transcriptional activation and expression of IGF1. Our novel results demonstrated a direct positive regulation of CTCF on the transcription of the IGF1 promoter. The impaired IGF1 expression due to CTCF mutation may explain the substandard effect of rhGH treatment on MRD21 patients. This study provided novel insights into the molecular basis of CTCF-associated disorder.

Isolation of High-Quality Total RNA from Small Animal Articular Cartilage for Next-Generation Sequencing

Articular cartilage is characterized by a low density of chondrocytes surrounded by an abundant extracellular matrix (ECM) consisting of a dense mixture of collagens, proteoglycans, and glycosaminoglycans. Due to its low cellularity and high proteoglycan content, it is particularly challenging to extract high-quality total RNA suitable for sensitive high-throughput downstream applications such as RNA sequencing (RNA-Seq). Available protocols for high-quality RNA isolation from articular chondrocytes are inconsistent, resulting in suboptimal yield and compromised quality. This poses a significant difficulty in the application of RNA-Seq to study the cartilage transcriptome. Current protocols rely either on dissociation of cartilage ECM by collagenase digestion or pulverizing cartilage using various methods prior to RNA extraction. However, protocols for cartilage processing vary significantly depending on the species and source of cartilage within the body. Protocols for isolating RNA from human or large mammal (e.g., horse or cattle) cartilage samples are available, but this is not the case for chicken cartilage, despite the species being extensively used in cartilage research. Here, we present two improved RNA isolation protocols based on pulverization of fresh articular cartilage using a cryogenic mill or on enzymatic digestion using 1.2% (w/v) collagenase II. Our protocols optimize the collection and tissue processing steps to minimize RNA degradation and enhance RNA purity. Our results show that RNA purified from chicken articular cartilage using these methods has appropriate quality for RNA-Seq experiments. The procedure is applicable for RNA extraction from cartilage from other species such as dog, cat, sheep, and goat. The workflow for RNA-Seq analysis is also described here. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Extraction of total RNA from pulverized chicken articular cartilage Alternate Protocol: Extraction of total RNA from collagen-digested articular cartilage Support Protocol: Dissection of chicken articular cartilage from the knee joint Basic Protocol 2: RNA sequencing of total RNA from chicken articular cartilage.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

Whole genome sequencing of Luxi Black Head sheep for screening selection signatures associated with important traits

Objective

Luxi Black Head sheep (LBH) is the first crossbreed specialized for meat production and was developed by crossbreeding Black Head Dorper sheep (DP) and Small Tailed Han sheep (STH) in the farming areas of northern China. Research on the genomic variations and selection signatures of LBH caused by continuous artificial selection is of great significance for identifying the genetic mechanisms of important traits of sheep and for the continuous breeding of LBH.

Methods

We explored the genetic relationships of LBH, DP, and several Mongolian sheep breeds by constructing phylogenetic tree, principal component analysis and linkage disequilibrium analysis. In addition, we analysed 29 whole genomes of sheep. The genome-wide selection signatures have been scanned with four methods: heterozygosity (H_P), fixation index (F_ST), cross-population extended haplotype homozygosity (XP-EHH) and the nucleotide diversity (θ_π) ratio.

Results

The genetic relationships analysis showed that LBH appeared to be an independent cluster closer to DP. The candidate signatures of positive selection in sheep genome revealed candidate genes for developmental process (HoxA gene cluster, BCL2L11, TSHR), immunity (CXCL6, CXCL1, SKAP2, PTK6, MST1R), growth (PDGFD, FGF18, SRF, SOCS2), and reproduction (BCAS3, TRIM24, ASTL, FNDC3A). Moreover, two signalling pathways closely related to reproduction, the thyroid hormone signalling pathway and the oxytocin signalling pathway, were detected.

Conclusion

The selective sweep analysis of LBH genome revealed candidate genes and signalling pathways associated with developmental process, immunity, growth, and reproduction. Our findings provide a valuable resource for sheep breeding and insight into the mechanisms of artificial selection.

HOXD9‑induced SCNN1A upregulation promotes pancreatic cancer cell proliferation, migration and predicts prognosis by regulating epithelial‑mesenchymal transformation

Pancreatic cancer (PC) is a malignant tumor disease, whose molecular mechanism is not fully understood. Sodium channel epithelial 1α subunit (SCNN1A) serves an important role in tumor progression. The current study explored the role of homeobox D9 (HOXD9) and SCNN1A in the progression of PC. The expression of SCNN1A and HOXD9 in PC samples was predicted on online databases and detected in PC cell lines. The association between SCNN1A expression and PC prognosis was examined by the Gene Expression Profiling Interactive Analysis, The Cancer Genome Atlas and Genotype‑Tissue Expression databases and by a Kaplan‑Meier plotter. Subsequently, the biological effects of SCNN1A on PC cell growth, colony formation, migration and invasion were investigated through RNA interference and cell transfection. Next, the expression of E‑cadherin, N‑cadherin, Vimentin and Snail was detected by western blotting to discover whether HOXD9 dysregulation mediated PC metastasis. Binding sites of HOXD9 and SCNN1A promoters were predicted on JASPAR. Reverse transcription‑quantitative PCR and western blotting were used to detect the expression level of SCNN1A following interference and overexpression of HOXD9. Luciferase assay detected luciferase activity following interference with HOXD9 and the transcriptional activity of SCNN1A following binding site deletion. High expression of SCNN1A and HOXD9 in PC was predicted by online databases, signifying poor prognosis. The present study confirmed the above predictions in PC cell lines. Knockdown of SCNN1A and HOXD9 could effectively inhibit the proliferation, migration, invasion and epithelial‑mesenchymal transition of PC cells. Furthermore, HOXD9 activated SCNN1A transcription, forming a feedback regulatory loop. HOXD9 was demonstrated to activate SCNN1A and promote the malignant biological process of PC.

Mesomelic dysplasias associated with the HOXD locus are caused by regulatory reallocations

Human families with chromosomal rearrangements at 2q31, where the human HOXD locus maps, display mesomelic dysplasia, a severe shortening and bending of the limb. In mice, the dominant Ulnaless inversion of the HoxD cluster produces a similar phenotype suggesting the same origin for these malformations in humans and mice. Here we engineer 1 Mb inversion including the HoxD gene cluster, which positioned Hoxd13 close to proximal limb enhancers. Using this model, we show that these enhancers contact and activate Hoxd13 in proximal cells, inducing the formation of mesomelic dysplasia. We show that a secondary Hoxd13 null mutation in-cis with the inversion completely rescues the alterations, demonstrating that ectopic HOXD13 is directly responsible for this bone anomaly. Single-cell expression analysis and evaluation of HOXD13 binding sites suggests that the phenotype arises primarily by acting through genes normally controlled by HOXD13 in distal limb cells. Altogether, these results provide a conceptual and mechanistic framework to understand and unify the molecular origins of human mesomelic dysplasia associated with 2q31.

Control of mouse limb initiation and antero-posterior patterning by Meis transcription factors

Meis1 and Meis2 are homeodomain transcription factors that regulate organogenesis through cooperation with Hox proteins. Elimination of Meis genes after limb induction has shown their role in limb proximo-distal patterning; however, limb development in the complete absence of Meis function has not been studied. Here, we report that Meis1/2 inactivation in the lateral plate mesoderm of mouse embryos leads to limb agenesis. Meis and Tbx factors converge in this function, extensively co-binding with Tbx to genomic sites and co-regulating enhancers of Fgf10, a critical factor in limb initiation. Limbs with three deleted Meis alleles show proximal-specific skeletal hypoplasia and agenesis of posterior skeletal elements. This failure in posterior specification results from an early role of Meis factors in establishing the limb antero-posterior prepattern required for Shh activation. Our results demonstrate roles for Meis transcription factors in early limb development and identify their involvement in previously undescribed interaction networks that regulate organogenesis.

Androgen and glucocorticoid receptor direct distinct transcriptional programs by receptor-specific and shared DNA binding sites

The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied them in an equivalent cellular context. Analysis of chromatin and sequence suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the result of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.

Folding Keratin Gene Clusters during Skin Regional Specification

Regional specification is critical for skin development, regeneration, and evolution. The contribution of epigenetics in this process remains unknown. Here, using avian epidermis, we find two major strategies regulate β-keratin gene clusters. (1) Over the body, macro-regional specificities (scales, feathers, claws, etc.) established by typical enhancers control five subclusters located within the epidermal differentiation complex on chromosome 25; (2) within a feather, micro-regional specificities are orchestrated by temporospatial chromatin looping of the feather β-keratin gene cluster on chromosome 27. Analyses suggest a three-factor model for regional specification: competence factors (e.g., AP1) make chromatin accessible, regional specifiers (e.g., Zic1) target specific genome regions, and chromatin regulators (e.g., CTCF and SATBs) establish looping configurations. Gene perturbations disrupt morphogenesis and histo-differentiation. This chicken skin paradigm advances our understanding of how regulation of big gene clusters can set up a two-dimensional body surface map.

Differential abilities to engage inaccessible chromatin diversify vertebrate Hox binding patterns

Although Hox genes encode for conserved transcription factors (TFs), they are further divided into anterior, central and posterior groups based on their DNA-binding domain similarity. The posterior Hox group expanded in the deuterostome clade and patterns caudal and distal structures. We aimed to address how similar Hox TFs diverge to induce different positional identities. We studied Hox TF DNA-binding and regulatory activity during an in vitro motor neuron differentiation system that recapitulates embryonic development. We found diversity in the genomic binding profiles of different Hox TFs, even among the posterior group paralogs that share similar DNA-binding domains. These differences in genomic binding were explained by differing abilities to bind to previously inaccessible sites. For example, the posterior group HOXC9 had a greater ability to bind occluded sites than the posterior HOXC10, producing different binding patterns and driving differential gene expression programs. From these results, we propose that the differential abilities of posterior Hox TFs to bind to previously inaccessible chromatin drive patterning diversification.This article has an associated 'The people behind the papers' interview.

Unblending of Transcriptional Condensates in Human Repeat Expansion Disease

Expansions of amino acid repeats occur in >20 inherited human disorders, and many occur in intrinsically disordered regions (IDRs) of transcription factors (TFs). Such diseases are associated with protein aggregation, but the contribution of aggregates to pathology has been controversial. Here, we report that alanine repeat expansions in the HOXD13 TF, which cause hereditary synpolydactyly in humans, alter its phase separation capacity and its capacity to co-condense with transcriptional co-activators. HOXD13 repeat expansions perturb the composition of HOXD13-containing condensates in vitro and in vivo and alter the transcriptional program in a cell-specific manner in a mouse model of synpolydactyly. Disease-associated repeat expansions in other TFs (HOXA13, RUNX2, and TBP) were similarly found to alter their phase separation. These results suggest that unblending of transcriptional condensates may underlie human pathologies. We present a molecular classification of TF IDRs, which provides a framework to dissect TF function in diseases associated with transcriptional dysregulation.

HOX13-dependent chromatin accessibility underlies the transition towards the digit development program

Hox genes encode transcription factors (TFs) that establish morphological diversity in the developing embryo. The similar DNA-binding motifs of the various HOX TFs contrast with the wide-range of HOX-dependent genetic programs. The influence of the chromatin context on HOX binding specificity remains elusive. Here, we used the developing limb as a model system to compare the binding specificity of HOXA13 and HOXD13 (HOX13 hereafter), which are required for digit formation, and HOXA11, involved in forearm/leg development. We find that upon ectopic expression in distal limb buds, HOXA11 binds sites normally HOX13-specific. Importantly, these sites are loci whose chromatin accessibility relies on HOX13. Moreover, we show that chromatin accessibility specific to the distal limb requires HOX13 function. Based on these results, we propose that HOX13 TFs pioneer the distal limb-specific chromatin accessibility landscape for the proper implementation of the distal limb developmental program.

Pioneer factors direct cell fate through switching inaccessible chromatin to an accessible state at specific target enhancers. Here the authors show that HOX13 transcription factors have a pioneer activity which is required for the proper implementation of the distal limb developmental program.

HOXD9 Activates the TGF-β/Smad Signaling Pathway to Promote Gastric Cancer

Background

Gastric cancer (GC) is the most common malignant tumor of the digestive tract and its molecular mechanism is not clear. HOXD9 plays an important role in tumor progression as transcription factor. In the current study, we explored the role of HOXD9 in GC.

Methods

We predicted the expression and potential mechanism of HOXD9 in GC through an online database. The expression of HOXD9 was detected in GC and adjacent tissues, and then we analyzed the relationship between HOXD9 and the prognosis of patients with GC. In vitro, we investigated the effects of HOXD9 on malignant biological behaviors such as proliferation, migration, and invasion of the GC cell line MCG-803. In addition, we have initially studied the underlying mechanism by Western blot.

Results

High expression of HOXD9 in GC was predicted by online database prediction and implied poor prognosis. In the clinical sample, we confirmed the above predictions. In vitro, we found that knockdown of HOXD9 could effectively inhibit the proliferation, migration, and invasion of GC cells. In terms of mechanism, HOXD9 may activate the TGF-β/Smad signaling pathway.

Conclusion

HOXD9 promotes the malignant biological process of GC, which may be a potential therapeutic target for GC.

CGGBP1 regulates CTCF occupancy at repeats

Background

CGGBP1 is a repeat-binding protein with diverse functions in the regulation of gene expression, cytosine methylation, repeat silencing and genomic integrity. CGGBP1 has also been identified as a cooperator of histone-modifying enzymes and as a component of CTCF-containing complexes that regulate the enhancer–promoter looping. CGGBP1–CTCF cross talk in chromatin regulation has been hitherto unknown.

Results

Here, we report that the occupancy of CTCF at repeats depends on CGGBP1. Using ChIP-sequencing for CTCF, we describe its occupancy at repetitive DNA. Our results show that endogenous level of CGGBP1 ensures CTCF occupancy preferentially on repeats over canonical CTCF motifs. By combining CTCF ChIP-sequencing results with ChIP sequencing for three different kinds of histone modifications (H3K4me3, H3K9me3 and H3K27me3), we show that the CGGBP1-dependent repeat-rich CTCF-binding sites regulate histone marks in flanking regions.

Conclusion

CGGBP1 affects the pattern of CTCF occupancy. Our results posit CGGBP1 as a regulator of CTCF and its binding sites in interspersed repeats.

Analysis of the 4q35 chromatin organization reveals distinct long-range interactions in patients affected with Facio-Scapulo-Humeral Dystrophy

Facio-Scapulo Humeral dystrophy (FSHD) is the third most common myopathy, affecting 1 amongst 10,000 individuals (FSHD1, OMIM #158900). This autosomal dominant pathology is associated in 95% of cases with genetic and epigenetic alterations in the subtelomeric region at the extremity of the long arm of chromosome 4 (q arm). A large proportion of the remaining 5% of cases carry a mutation in the SMCHD1 gene (FSHD2, OMIM #158901). Here, we explored the 3D organization of the 4q35 locus by three-dimensions DNA in situ fluorescent hybridization (3D-FISH) in primary fibroblasts isolated from patients and healthy donors. We found that D4Z4 contractions and/or SMCHD1 mutations impact the spatial organization of the 4q35 region and trigger changes in the expression of different genes. Changes in gene expression were corroborated in muscle biopsies suggesting that the modified chromatin landscape impelled a modulation in the level of expression of a number of genes across the 4q35 locus in FSHD. Using induced pluripotent stem cells (hIPSC), we further examined whether chromatin organization is inherited after reprogramming or acquired during differentiation and showed that folding of the 4q35 region is modified upon differentiation. These results together with previous findings highlight the role of the D4Z4 macrosatellite repeat in the topological organization of chromatin and further indicate that the D4Z4-dependent 3D structure induces transcriptional changes of 4q35 genes expression.

Preformed chromatin topology assists transcriptional robustness of Shh during limb development

Long-range gene regulation involves physical proximity between enhancers and promoters to generate precise patterns of gene expression in space and time. However, in some cases, proximity coincides with gene activation, whereas, in others, preformed topologies already exist before activation. In this study, we investigate the preformed configuration underlying the regulation of the Shh gene by its unique limb enhancer, the ZRS, in vivo during mouse development. Abrogating the constitutive transcription covering the ZRS region led to a shift within the Shh-ZRS contacts and a moderate reduction in Shh transcription. Deletion of the CTCF binding sites around the ZRS resulted in the loss of the Shh-ZRS preformed interaction and a 50% decrease in Shh expression but no phenotype, suggesting an additional, CTCF-independent mechanism of promoter-enhancer communication. This residual activity, however, was diminished by combining the loss of CTCF binding with a hypomorphic ZRS allele, resulting in severe Shh loss of function and digit agenesis. Our results indicate that the preformed chromatin structure of the Shh locus is sustained by multiple components and acts to reinforce enhancer-promoter communication for robust transcription.

Direct and Indirect Targeting of HOXA9 Transcription Factor in Acute Myeloid Leukemia

HOXA9 (Homeobox A9) is a homeotic transcription factor known for more than two decades to be associated with leukemia. The expression of HOXA9 homeoprotein is associated with anterior-posterior patterning during embryonic development, and its expression is then abolished in most adult cells, with the exception of hematopoietic progenitor cells. The oncogenic function of HOXA9 was first assessed in human acute myeloid leukemia (AML), particularly in the mixed-phenotype associated lineage leukemia (MPAL) subtype. HOXA9 expression in AML is associated with aggressiveness and a poor prognosis. Since then, HOXA9 has been involved in other hematopoietic malignancies and an increasing number of solid tumors. Despite this, HOXA9 was for a long time not targeted to treat cancer, mainly since, as a transcription factor, it belongs to a class of protein long considered to be an "undruggable" target; however, things have now evolved. The aim of the present review is to focus on the different aspects of HOXA9 targeting that could be achieved through multiple ways: (1) indirectly, through the inhibition of its expression, a strategy acting principally at the epigenetic level; or (2) directly, through the inhibition of its transcription factor function by acting at either the protein/protein interaction or the protein/DNA interaction interfaces.

'Building a perfect body': control of vertebrate organogenesis by PBX-dependent regulatory networks

Pbx genes encode transcription factors that belong to the TALE (three-amino-acid loop extension) superclass of homeodomain proteins. We have witnessed a surge in information about the roles of this gene family as leading actors in the transcriptional control of development. PBX proteins represent a clear example of how transcription factors can regulate developmental processes by combinatorial properties, acting within multimeric complexes to implement activation or repression of transcription depending on their interaction partners. Here, we revisit long-emphasized functions of PBX transcription factors as cofactors for HOX proteins, major architects of the body plan. We further discuss new knowledge on roles of PBX proteins in different developmental contexts as upstream regulators of Hox genes-as factors that interact with non-HOX proteins and can work independently of HOX-as well as potential pioneer factors. Committed to building a perfect body, PBX proteins govern regulatory networks that direct essential morphogenetic processes and organogenesis in vertebrate development. Perturbations of PBX-dependent networks can cause human congenital disease and cancer.

Hoxa13 regulates expression of common Hox target genes involved in cartilage development to coordinate the expansion of the autopodal anlage

To elucidate the role of Hox genes in limb cartilage development, we identified the target genes of HOXA11 and HOXA13 by ChIP‐Seq. The ChIP DNA fragment contained evolutionarily conserved sequences and multiple highly conserved HOX binding sites. A substantial portion of the HOXA11 ChIP fragment overlapped with the HOXA13 ChIP fragment indicating that both factors share common targets. Deletion of the target regions neighboring Bmp2 or Tshz2 reduced their expression in the autopod suggesting that they function as the limb bud‐specific enhancers. We identified the Hox downstream genes as exhibiting expression changes in the Hoxa13 knock out (KO) and Hoxd11‐13 deletion double mutant (Hox13 dKO) autopod by Genechip analysis. The Hox downstream genes neighboring the ChIP fragment were defined as the direct targets of Hox. We analyzed the spatial expression pattern of the Hox target genes that encode two different categories of transcription factors during autopod development and Hox13dKO limb bud. (a) Bcl11a, encoding a repressor of cartilage differentiation, was expressed in the E11.5 autopod and was substantially reduced in the Hox13dKO. (b) The transcription factors Aff3, Bnc2, Nfib and Runx1t1 were expressed in the zeugopodal cartilage but not in the autopod due to the repressive or relatively weak transcriptional activity of Hox13 at E11.5. Interestingly, the expression of these genes was later observed in the autopodal cartilage at E12.5. These results indicate that Hox13 transiently suspends the cartilage differentiation in the autopodal anlage via multiple pathways until establishing the paddle‐shaped structure required to generate five digits.

The Enigmatic HOX Genes: Can We Crack Their Code?

Homeobox genes (HOX) are a large family of transcription factors that direct the formation of many body structures during early embryonic development. There are 39 genes in the subgroup of homeobox genes that constitute the human HOX gene family. Correct embryonic development of flies and vertebrates is, in part, mediated by the unique and highly regulated expression pattern of the HOX genes. Disruptions in these fine-tuned regulatory mechanisms can lead to developmental problems and to human diseases such as cancer. Unfortunately, the molecular mechanisms of action of the HOX family of transcription factors are severely under-studied, likely due to idiosyncratic details of their structure, expression, and function. We suggest that a concerted and collaborative effort to identify interacting protein partners, produce genome-wide binding profiles, and develop HOX network inhibitors in a variety of human cell types will lead to a deeper understanding of human development and disease. Within, we review the technological challenges and possible approaches needed to achieve this goal.

Identification of a limb enhancer that is removed by pathogenic deletions downstream of the SHOX gene

Haploinsufficiency of the human SHOX gene causes Léri-Weill dyschondrosteosis (LWD), characterized by shortening of the middle segments of the limbs and Madelung deformity of the wrist. As many as 35% of LWD cases are caused by deletions of non-coding sequences downstream of SHOX that presumably remove an enhancer or enhancers necessary for SHOX expression in developing limbs. We searched for these active sequences using a transgenic mouse assay and identified a 563 basepair (bp) enhancer with specific activity in the limb regions where SHOX functions. This enhancer has previously escaped notice because of its poor evolutionary conservation, although it does contain 100 bp that are conserved in non-rodent mammals. A primary cell luciferase assay confirmed the enhancer activity of the conserved core sequence and demonstrated that putative HOX binding sites are required for its activity. This enhancer is removed in most non-coding deletions that cause LWD. However, we did not identify any likely pathogenic variants of the enhancer in a screen of 124 LWD individuals for whom no causative mutation had been found, suggesting that only larger deletions in the region commonly cause LWD. We hypothesize that loss of this enhancer contributes to the pathogenicity of deletions downstream of SHOX.

Problems in Fish-to-Tetrapod Transition: Genetic Expeditions Into Old Specimens

The fish-to-tetrapod transition is one of the fundamental problems in evolutionary biology. A significant amount of paleontological data has revealed the morphological trajectories of skeletons, such as those of the skull, vertebrae, and appendages in vertebrate history. Shifts in bone differentiation, from dermal to endochondral bones, are key to explaining skeletal transformations during the transition from water to land. However, the genetic underpinnings underlying the evolution of dermal and endochondral bones are largely missing. Recent genetic approaches utilizing model organisms—zebrafish, frogs, chickens, and mice—reveal the molecular mechanisms underlying vertebrate skeletal development and provide new insights for how the skeletal system has evolved. Currently, our experimental horizons to test evolutionary hypotheses are being expanded to non-model organisms with state-of-the-art techniques in molecular biology and imaging. An integration of functional genomics, developmental genetics, and high-resolution CT scanning into evolutionary inquiries allows us to reevaluate our understanding of old specimens. Here, we summarize the current perspectives in genetic programs underlying the development and evolution of the dermal skull roof, shoulder girdle, and appendages. The ratio shifts of dermal and endochondral bones, and its underlying mechanisms, during the fish-to-tetrapod transition are particularly emphasized. Recent studies have suggested the novel cell origins of dermal bones, and the interchangeability between dermal and endochondral bones, obscuring the ontogenetic distinction of these two types of bones. Assimilation of ontogenetic knowledge of dermal and endochondral bones from different structures demands revisions of the prevalent consensus in the evolutionary mechanisms of vertebrate skeletal shifts.

Functional roles of CTCF in breast cancer

CTCF, Zinc-finger protein, has been identified as a multifunctional transcription factor that regulates gene expression through various mechanisms, including recruitment of other co-activators and binding to promoter regions of target genes. Furthermore, it has been proposed to be an insulator protein that contributes to the establishment of functional three-dimensional chromatin structures. It can disrupt transcription through blocking the connection between an enhancer and a promoter. Previous studies revealed that the onset of various diseases, including breast cancer, could be attributed to the aberrant expression of CTCF itself or one or more of its target genes. In this review, we will describe molecular dysfunction involving CTCF that induces tumorigenesis and summarize the functional roles of CTCF in breast cancer.

Integrating ChIP-seq with other functional genomics data

Transcription is regulated by transcription factor (TF) binding at promoters and distal regulatory elements and histone modifications that control the accessibility of these elements. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) has become the standard assay for identifying genome-wide protein-DNA interactions in vitro and in vivo. As large-scale ChIP-seq data sets have been collected for different TFs and histone modifications, their potential to predict gene expression can be used to test hypotheses about the mechanisms of gene regulation. In addition, complementary functional genomics assays provide a global view of chromatin accessibility and long-range cis-regulatory interactions that are being combined with TF binding and histone remodeling to study the regulation of gene expression. Thus, ChIP-seq analysis is now widely integrated with other functional genomics assays to better understand gene regulatory mechanisms. In this review, we discuss advances and challenges in integrating ChIP-seq data to identify context-specific chromatin states associated with gene activity. We describe the overall computational design of integrating ChIP-seq data with other functional genomics assays. We also discuss the challenges of extending these methods to low-input ChIP-seq assays and related single-cell assays.

Master Regulatory Transcription Factors in Plant Development: A Blooming Perspective

Transcription factors that trigger major developmental decisions in plants and animals are termed "master regulators". Such master regulators are classically seen as acting on the top of a regulatory hierarchy that determines a complete developmental program, and they usually encode transcription factors. Here, we introduce master regulators of flowering time and flower development as examples to show how analysis of molecular interactions and gene-regulatory networks in plants has changed our view on the molecular mechanisms by which these factors control developmental processes. A picture has emerged that emphasizes a complex combinatorial interplay in determining cell-type transcriptional programs, and a high level of feedback control. The expression of master regulators themselves is usually regulated by multiple factors integrating environmental and endogenous spatiotemporal cues. Master regulatory transcription factors regulate gene expression by different mechanisms, including modifications in chromatin status in the bound regions. A poorly understood phenomenon is how developmental master regulators exert functions in different cell- and organ types. This is especially relevant for those factors that have important functions in several developmental processes.

Epigenetic profiling of growth plate chondrocytes sheds insight into regulatory genetic variation influencing height

GWAS have identified hundreds of height-associated loci. However, determining causal mechanisms is challenging, especially since height-relevant tissues (e.g. growth plates) are difficult to study. To uncover mechanisms by which height GWAS variants function, we performed epigenetic profiling of murine femoral growth plates. The profiled open chromatin regions recapitulate known chondrocyte and skeletal biology, are enriched at height GWAS loci, particularly near differentially expressed growth plate genes, and enriched for binding motifs of transcription factors with roles in chondrocyte biology. At specific loci, our analyses identified compelling mechanisms for GWAS variants. For example, at CHSY1, we identified a candidate causal variant (rs9920291) overlapping an open chromatin region. Reporter assays demonstrated that rs9920291 shows allelic regulatory activity, and CRISPR/Cas9 targeting of human chondrocytes demonstrates that the region regulates CHSY1 expression. Thus, integrating biologically relevant epigenetic information (here, from growth plates) with genetic association results can identify biological mechanisms important for human growth.

A Hox complex activates and potentiates the Epidermal Growth Factor signaling pathway to specify Drosophila oenocytes

Hox transcription factors specify distinct cell types along the anterior-posterior axis of metazoans by regulating target genes that modulate signaling pathways. A well-established example is the induction of Epidermal Growth Factor (EGF) signaling by an Abdominal-A (Abd-A) Hox complex during the specification of Drosophila hepatocyte-like cells (oenocytes). Previous studies revealed that Abd-A is non-cell autonomously required to promote oenocyte fate by directly activating a gene (rhomboid) that triggers EGF secretion from sensory organ precursor (SOP) cells. Neighboring cells that receive the EGF signal initiate a largely unknown pathway to promote oenocyte fate. Here, we show that Abd-A also plays a cell autonomous role in inducing oenocyte fate by activating the expression of the Pointed-P1 (PntP1) ETS transcription factor downstream of EGF signaling. Genetic studies demonstrate that both PntP1 and PntP2 are required for oenocyte specification. Moreover, we found that PntP1 contains a conserved enhancer (PntP1OE) that is activated in oenocyte precursor cells by EGF signaling via direct regulation by the Pnt transcription factors as well as a transcription factor complex consisting of Abd-A, Extradenticle, and Homothorax. Our findings demonstrate that the same Abd-A Hox complex required for sending the EGF signal from SOP cells, enhances the competency of receiving cells to select oenocyte cell fate by up-regulating PntP1. Since PntP1 is a downstream effector of EGF signaling, these findings provide insight into how a Hox factor can both trigger and potentiate the EGF signal to promote an essential cell fate along the body plan.

Psip1/p52 regulates posterior Hoxa genes through activation of lncRNA Hottip

Long noncoding RNAs (lncRNAs) have been implicated in various biological functions including the regulation of gene expression, however, the functionality of lncRNAs is not clearly understood and conflicting conclusions have often been reached when comparing different methods to investigate them. Moreover, little is known about the upstream regulation of lncRNAs. Here we show that the short isoform (p52) of a transcriptional co-activator-PC4 and SF2 interacting protein (Psip1), which is known to be involved in linking transcription to RNA processing, specifically regulates the expression of the lncRNA Hottip-located at the 5' end of the Hoxa locus. Using both knockdown and knockout approaches we show that Hottip expression is required for activation of the 5' Hoxa genes (Hoxa13 and Hoxa10/11) and for retaining Mll1 at the 5' end of Hoxa. Moreover, we demonstrate that artificially inducing Hottip expression is sufficient to activate the 5' Hoxa genes and that Hottip RNA binds to the 5' end of Hoxa. By engineering premature transcription termination, we show that it is the Hottip lncRNA molecule itself, not just Hottip transcription that is required to maintains active expression of posterior Hox genes. Our data show a direct role for a lncRNA molecule in regulating the expression of developmentally-regulated mRNA genes in cis.