Hoxc8 initiates an ectopic mammary program by regulating Fgf10 and Tbx3 expression and Wnt/β-catenin signaling

Dynamic changes of genomic methylation profiles at different growth stages in Chinese Tan sheep

Background

Tan sheep, an important local sheep breed in China, is famous for their fur quality. One-month-old Tan sheep have white, curly hair with beautiful flower spikes, commonly known as “nine bends”, which has high economic value. However, the “nine bends” characteristic gradually disappears with age; consequently, the economic value of the Tan sheep decreases. Age-related changes in DNA methylation have been reported and may be responsible for age-induced changes in gene expression. Until now, no genome-wide surveys have been conducted to identify potential DNA methylation sites involved in different sheep growth stages. In this study we investigated the dynamic changes of genome-wide DNA methylation profiles in Tan sheep using DNA from skin and deep whole-genome bisulfite sequencing, and compared the DNA methylation levels at three different growth stages: 1, 24, and 48 months old (mon1, mon24, and mon48, respectively).

Results

In this study, 11 skin samples from three growth stages (four for mon1, four for mon24, and three for mon48) were used for DNA methylation analysis and gene expression profiling. There were 52, 288 and 236 differentially methylated genes (DMGs) identified between mon1 and mon24, mon1 and mon48, and mon24 and mon48, respectively. Of the differentially methylated regions, 1.11%, 7.61%, and 7.65% were in the promoter in mon1 vs. mon24, mon24 vs. mon48, and mon1 vs. mon48, respectively. DMGs were enriched in the MAPK and WNT signaling pathways, which are related to age growth and hair follicle morphogenesis processes. There were 51 DMGs associated with age growth and curly fleece formation. Four DMGs between mon1 and mon48 (KRT71, CD44, ROR2 and ZDHHC13) were further validated by bisulfite sequencing.

Conclusions

This study revealed dynamic changes in the genomic methylation profiles of mon1, mon24, and mon48 sheep, and the percentages of methylated cytosines were 3.38%, 2.85% and 4.17%, respectively. Of the DMGs, KRT71 and CD44 were highly methylated in mon1, and ROR2 and ZDHHC13 were highly methylated in mon48. These findings provide foundational information that may be used to develop strategies for potentially retaining the lamb fur and thus improving the economic value of Tan sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00632-9.

The molecular basis of mammary gland development and epithelial differentiation

Our understanding of the molecular events underpinning the development of mammalian organ systems has been increasing rapidly in recent years. With the advent of new and improved next-generation sequencing methods, we are now able to dig deeper than ever before into the genomic and epigenomic events that play critical roles in determining the fates of stem and progenitor cells during the development of an embryo into an adult. In this review, we detail and discuss the genes and pathways that are involved in mammary gland development, from embryogenesis, through maturation into an adult gland, to the role of pregnancy signals in directing the terminal maturation of the mammary gland into a milk producing organ that can nurture the offspring. We also provide an overview of the latest research in the single-cell genomics of mammary gland development, which may help us to understand the lineage commitment of mammary stem cells (MaSCs) into luminal or basal epithelial cells that constitute the mammary gland. Finally, we summarize the use of 3D organoid cultures as a model system to study the molecular events during mammary gland development. Our increased investigation of the molecular requirements for normal mammary gland development will advance the discovery of targets to predict breast cancer risk and the development of new breast cancer therapies.

MicroRNA governs bistable cell differentiation and lineage segregation via a noncanonical feedback

Positive feedback driven by transcriptional regulation has long been considered a key mechanism underlying cell lineage segregation during embryogenesis. Using the developing spinal cord as a paradigm, we found that canonical, transcription-driven feedback cannot explain robust lineage segregation of motor neuron subtypes marked by two cardinal factors, Hoxa5 and Hoxc8. We propose a feedback mechanism involving elementary microRNA-mRNA reaction circuits that differ from known feedback loop-like structures. Strikingly, we show that a wide range of biologically plausible post-transcriptional regulatory parameters are sufficient to generate bistable switches, a hallmark of positive feedback. Through mathematical analysis, we explain intuitively the hidden source of this feedback. Using embryonic stem cell differentiation and mouse genetics, we corroborate that microRNA-mRNA circuits govern tissue boundaries and hysteresis upon motor neuron differentiation with respect to transient morphogen signals. Our findings reveal a previously underappreciated feedback mechanism that may have widespread functions in cell fate decisions and tissue patterning.

Mammary development in the embryo and adult: new insights into the journey of morphogenesis and commitment

The mammary gland is a unique tissue and the defining feature of the class Mammalia. It is a late-evolving epidermal appendage that has the primary function of providing nutrition for the young, although recent studies have highlighted additional benefits of milk including the provision of passive immunity and a microbiome and, in humans, the psychosocial benefits of breastfeeding. In this Review, we outline the various stages of mammary gland development in the mouse, with a particular focus on lineage specification and the new insights that have been gained by the application of recent technological advances in imaging in both real-time and three-dimensions, and in single cell RNA sequencing. These studies have revealed the complexity of subpopulations of cells that contribute to the mammary stem and progenitor cell hierarchy and we suggest a new terminology to distinguish these cells.

Polycomb Repressive Complex 2: a Dimmer Switch of Gene Regulation in Calvarial Bone Development

PURPOSE OF REVIEW

Epigenetic regulation is a distinct mechanism of gene regulation that functions by modulating chromatin structure and accessibility. Polycomb Repressive Complex 2 (PRC2) is a conserved chromatin regulator that is required in the developing embryo to control the expression of key developmental genes. An emerging feature of PRC2 is that it not only allows for binary ON/OFF states of gene expression but can also modulate gene expression in feed-forward loops to change the outcome of gene regulatory networks. This striking feature of epigenetic modulation has improved our understanding of musculoskeletal development.

RECENT FINDINGS

Recent advances in mouse embryos unravel a range of phenotypes that demonstrate the tissue-specific, temporal, and spatial role of PRC2 during organogenesis and cell fate decisions in vivo. Here, we take a detailed view of how PRC2 functions during the development of calvarial bone and skin. Based on the emerging evidence, we propose that PRC2 serves as a "dimmer switch" to modulate gene expression of target genes by altering the expression of activators and inhibitors. This review highlights the findings from contemporary research that allow us to investigate the unique developmental potential of intramembranous calvarial bones.

Fgf10/Fgfr2b Signaling in Mammary Gland Development, Homeostasis, and Cancer

Fibroblast growth factor 10 (Fgf10) is a secreted ligand acting via the Fibroblast growth factor receptor 2b (Fgfr2b). Fgf10/Fgfr2b signaling plays important roles both in the epithelium and in the mesenchyme during mammary gland development. Evidence in mice show that Fgf10 is critical for the induction of four out of five of the mammary placodes and for the formation of the white adipose tissue. Fgfr2b ligands also play important function in the maintenance of the terminal end buds, specialized structures at the tip of the ramified ducts during the postnatal phase of mammary gland development. Finally, in humans, FGF10 has been described to be expressed in 10% of the breast adenocarcinoma and activation of FGFR2b signaling correlates with a worse prognostic. Therefore, Fgf10 plays pleiotropic roles in both mammary gland development, homeostasis and cancer and elucidating its mechanism of action and cellular targets will be crucial to either enhance mammary gland development or to find innovative targets to treat aggressive breast cancer.

HOX genes function in Breast Cancer development

Breast cancer develops in the mammary glands during mammalian adulthood and is considered the second most common type of human carcinoma and the most incident and mortal in the female population. In contrast to other human structures, the female mammary glands continue to develop after birth, undergoing various modifications during pregnancy, lactation and involution under the regulation of hormones and transcription factors, including those encoded by the HOX clusters (A, B, C, and D). Interestingly, HOX gene deregulation is often associated to breast cancer development. Within the HOXB cluster, 8 out of the 10 genes present altered expression levels in breast cancer with an impact in its aggressiveness and resistance to hormone therapy, which highlights the importance of HOXB genes as potential therapeutic targets used to overcome the limitations of tamoxifen-resistant cancer treatments. Here, we review the current state of knowledge on the role of HOX genes in breast cancer, specially focus on HOXB, discussing the causes and consequences of HOXB gene deregulation and their relevance as prognostic factors and therapeutic targets.

Expression profile analysis of dermal papilla cells mRNA in response to WNT10B treatment

Dermal papilla cells (DPCs) are associated with the development of hair follicles (HFs) and the regulation of the hair growth cycle. Previous studies have shown that Wnt family member 10B (WNT10B) plays an important role in the proliferation and survival of DPCs in vitro, and promotes the growth of HFs. However, the underlying mechanisms have not been fully elucidated. The present study evaluated the role of WNT10B in regulating HF morphogenesis by characterizing the differential gene expression profiles between WNT10B-treated DPCs and control DPCs using RNA-sequencing (RNA-seq). A total of 1,073 and 451 genes were upregulated and downregulated, respectively. The RNA-seq data was subsequently validated by reverse-transcription quantitative PCR. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that 442 GO terms and 21 KEGG pathways were significantly enriched. Further functional analysis revealed that WNT10B decreased translation initiation, elongation and termination, and RNA metabolic processes in cultured DPCs compared with controls in vitro. Human signaling networks were compared using pathway analysis, and treatment of DPCs with WNT10B was revealed to downregulate the ribosome biogenesis pathway and decrease protein synthesis in vitro. KEGG pathway analysis showed that WNT10B upregulated the phosphoinositide 3-kinase/protein kinase B signaling pathway. The present study analyzed the expression of mRNA in WNT10B-treated DPCs using next-generation sequencing and uncovered mechanisms regulating the induction of HFs.

A tale of two cities: The genetic mechanisms governing calvarial bone development

The skull bones must grow in a coordinated, three-dimensional manner to coalesce and form the head and face. Mammalian skull bones have a dual embryonic origin from cranial neural crest cells (CNCC) and paraxial mesoderm (PM) and ossify through intramembranous ossification. The calvarial bones, the bones of the cranium which cover the brain, are derived from the supraorbital arch (SOA) region mesenchyme. The SOA is the site of frontal and parietal bone morphogenesis and primary center of ossification. The objective of this review is to frame our current in vivo understanding of the morphogenesis of the calvarial bones and the gene networks regulating calvarial bone initiation in the SOA mesenchyme.

Stage-specific roles of Ezh2 and Retinoic acid signaling ensure calvarial bone lineage commitment

Development of the skull bones requires the coordination of two stem progenitor populations, the cranial neural crest cells (CNCC) and head paraxial mesoderm (PM), to ensure cell fate selection and morphogenesis. The epigenetic methyltransferase, Ezh2, plays a role in skull bone formation, but the spatiotemporal function of Ezh2 between the CNCC- and PM-derived bone formation in vivo remains undefined. Here, using a temporally-inducible conditional deletion of Ezh2 in both the CNCC- and PM- derived cranial mesenchyme between E8.5 and E9.5, we find a reduction of the CNCC-derived calvarial bones and a near complete loss of the PM-derived calvarial bones due to an arrest in calvarial bone fate commitment. In contrast, deletion of Ezh2 after E9.5 permits PM-derived skull bone development, suggesting that Ezh2 is required early to guide calvarial bone progenitor commitment. Furthermore, exposure to all-trans Retinoic acid at E10.0 can mimic the Ezh2 mutant calvarial phenotype, and administration of the pan retinoic acid receptor (RAR) antagonist, BMS-453, to Ezh2 mutants partially restores the commitment to the calvarial bone lineage and PM-derived bone development in vivo. Exogenous RA signaling activation in the Ezh2 mutants leads to synergistic activation of the anti-osteogenic factors in the cranial mesenchyme in vivo. Thus, RA signaling and EZH2 can function in parallel to guide calvarial bone progenitor commitment by balancing the suppression of anti-osteogenic factors.

Improvement Activity of 1-Deoxynojirimycin in the Growth of Dairy Goat Primary Mammary Epithelial Cell through Upregulating LEF-1 Expression

LEF-1/wnt10b is one of the most important signaling pathways regulating mammary gland growth and development and is also a potential target for molecular breeding. In this work, 1-deoxynojirimycin (DNJ), a natural alkaloid extracted from plant mulberry or microorganism, was found to have a positive activity in primary breast epithelial cell growth of dairy goats. The findings showed that, compared to the control, 6 μM DNJ in the DMEM/F12 medium in vitro greatly improved the density of dairy goat breast epithelial cell and significantly increased the LEF-1 mRNA level (P < 0.01) and thus enhanced cell growth. In addition, DNJ displayed a similar function in alleviating the growth suppression of epithelial cell and the decrease of LEF-1 mRNA level resulting from lentiviral-mediated LEF-1 knockdown. Simultaneously, no significant change of the mRNA levels of IGF-1 and Fgf10, the other two key regulators in mammary gland growth and development, could be detected. Furthermore, the mammary duct of DNJ-fed mouse illustrated a better development accompanied with a higher LET-1 mRNA level than that of the control. In conclusion, DNJ could improve breast epithelial cell growth through upregulating LEF-1 expression, which supplied a new means in studying mammary gland growth and development.

HOXC8 promotes proliferation and migration through transcriptional up-regulation of TGFβ1 in non-small cell lung cancer

Homeobox (HOX) genes encode a family of transcription factors, which play crucial roles in numerous processes, and their dysregulation is involved in the carcinogenesis of many human cancers. In the present study, we investigated the roles of HOXC8 in non-small cell lung cancer (NSCLC). We showed that HOXC8 was upregulated in clinical NSCLC specimens compared to normal lung tissues, and the high expression of HOXC8 correlated with tumor node metastasis (TNM) stage, tumor status, lymph nodal status and poor relapse-free survival for lung cancer patients. Functionally, HOXC8 expression significantly promoted the proliferation, anchorage-independent growth and migration of NSCLC, and HOXC8 functioned as a transcription activator to induce the expression of TGFβ1, leading to an increase in the proliferation, anchorage-independent growth and migration of NSCLC. Furthermore, we demonstrated that HOXC8 expression was associated with chemoresistance and anti-apoptosis in NSCLC, suggesting that HOXC8 is a promising therapeutic target for chemosensitization of NSCLC to cisplatin. Altogether, our study defined a critical role of HOXC8 in promoting transcription of TGFβ1 and NSCLC tumorigenesis.

Conditional Loss of Hoxa5 Function Early after Birth Impacts on Expression of Genes with Synaptic Function

Hoxa5 is a member of the Hox gene family that plays critical roles in successive steps of the central nervous system formation during embryonic and fetal development. In the mouse, Hoxa5 was recently shown to be expressed in the medulla oblongata and the pons from fetal stages to adulthood. In these territories, Hoxa5 transcripts are enriched in many precerebellar neurons and several nuclei involved in autonomic functions, while the HOXA5 protein is detected mainly in glutamatergic and GABAergic neurons. However, whether HOXA5 is functionally required in these neurons after birth remains unknown. As a first approach to tackle this question, we aimed at determining the molecular programs downstream of the HOXA5 transcription factor in the context of the postnatal brainstem. A comparative transcriptomic analysis was performed in combination with gene expression localization, using a conditional postnatal Hoxa5 loss-of-function mouse model. After inactivation of Hoxa5 at postnatal days (P)1–P4, we established the transcriptome of the brainstem from P21 Hoxa5 conditional mutants using RNA-Seq analysis. One major finding was the downregulation of several genes associated with synaptic function in Hoxa5 mutant specimens including different actors involved in glutamatergic synapse, calcium signaling pathway, and GABAergic synapse. Data were confirmed and extended by reverse transcription quantitative polymerase chain reaction analysis, and the expression of several HOXA5 candidate targets was shown to co-localize with Hoxa5 transcripts in precerebellar nuclei. Together, these new results revealed that HOXA5, through the regulation of key actors of the glutamatergic/GABAergic synapses and calcium signaling, might be involved in synaptogenesis, synaptic transmission, and synaptic plasticity of the cortico-ponto-cerebellar circuitry in the postnatal brainstem.

HOXC8 regulates self-renewal, differentiation and transformation of breast cancer stem cells

BACKGROUND

Homeobox genes are master regulators of cell fate during embryonic development and their expression is altered in cancer. By regulating the balance between cell proliferation and differentiation, they maintain homeostasis of normal tissues. Here, we screened the expression of homeobox genes in mammary stem cells to establish their role in stem cells transformation in breast cancer.

METHODS

Using a Homeobox Genes PCR array, we screened 83 homeobox genes in normal cancer breast stem/progenitor cells isolated by flow cytometry. The candidate gene HOXC8 epigenetic regulation was studied by DNA methylation and miRNA expression analyses. Self-renewal and differentiation of HOXC8-overexpressing or knockdown cells were assessed by flow cytometry and mammosphere, 3D matrigel and soft agar assays. Clinical relevance of in vitro findings were validated by bioinformatics analysis of patient datasets from TCGA and METABRIC studies.

RESULTS

In this study we demonstrate altered expression of homeobox genes in breast cancer stem/progenitor cells. HOXC8 was consistently downregulated in stem/progenitor cells of all breast molecular subtypes, thus representing an interesting tumour suppressor candidate. We show that downregulated expression of HOXC8 is associated with DNA methylation at the gene promoter and expression of miR196 family members. Functional studies demonstrated that HOXC8 gain of function induces a decrease in the CD44⁺/CD24^-/low cancer stem cell population and proportion of chemoresistant cells, with a concomitant increase in CD24⁺ differentiated cells. Increased HOXC8 levels also decrease the ability of cancer cells to form mammospheres and to grow in anchorage-independent conditions. Furthermore, loss of HOXC8 in non-tumorigenic mammary epithelial cells expands the cancer stem/progenitor cells pool, increases stem cell self-renewal, prevents differentiation induced by retinoic acid and induces a transformed phenotype.

CONCLUSIONS

Taken together, our study points to an important role of homeobox genes in breast cancer stem/progenitor cell function and establishes HOXC8 as a suppressor of stemness and transformation in the mammary gland lineage.

Control of Hoxd gene transcription in the mammary bud by hijacking a preexisting regulatory landscape

Vertebrate Hox genes encode transcription factors operating during the development of multiple organs and structures. However, the evolutionary mechanism underlying this remarkable pleiotropy remains to be fully understood. Here, we show that Hoxd8 and Hoxd9, two genes of the HoxD complex, are transcribed during mammary bud (MB) development. However, unlike in other developmental contexts, their coexpression does not rely on the same regulatory mechanism. Hoxd8 is regulated by the combined activity of closely located sequences and the most distant telomeric gene desert. On the other hand, Hoxd9 is controlled by an enhancer-rich region that is also located within the telomeric gene desert but has no impact on Hoxd8 transcription, thus constituting an exception to the global regulatory logic systematically observed at this locus. The latter DNA region is also involved in Hoxd gene regulation in other contexts and strongly interacts with Hoxd9 in all tissues analyzed thus far, indicating that its regulatory activity was already operational before the appearance of mammary glands. Within this DNA region and neighboring a strong limb enhancer, we identified a short sequence conserved in therian mammals and capable of enhancer activity in the MBs. We propose that Hoxd gene regulation in embryonic MBs evolved by hijacking a preexisting regulatory landscape that was already at work before the emergence of mammals in structures such as the limbs or the intestinal tract.

A Complex Structural Variation on Chromosome 27 Leads to the Ectopic Expression of HOXB8 and the Muffs and Beard Phenotype in Chickens

Muffs and beard (Mb) is a phenotype in chickens where groups of elongated feathers gather from both sides of the face (muffs) and below the beak (beard). It is an autosomal, incomplete dominant phenotype encoded by the Muffs and beard (Mb) locus. Here we use genome-wide association (GWA) analysis, linkage analysis, Identity-by-Descent (IBD) mapping, array-CGH, genome re-sequencing and expression analysis to show that the Mb allele causing the Mb phenotype is a derived allele where a complex structural variation (SV) on GGA27 leads to an altered expression of the gene HOXB8. This Mb allele was shown to be completely associated with the Mb phenotype in nine other independent Mb chicken breeds. The Mb allele differs from the wild-type mb allele by three duplications, one in tandem and two that are translocated to that of the tandem repeat around 1.70 Mb on GGA27. The duplications contain total seven annotated genes and their expression was tested during distinct stages of Mb morphogenesis. A continuous high ectopic expression of HOXB8 was found in the facial skin of Mb chickens, strongly suggesting that HOXB8 directs this regional feather-development. In conclusion, our results provide an interesting example of how genomic structural rearrangements alter the regulation of genes leading to novel phenotypes. Further, it again illustrates the value of utilizing derived phenotypes in domestic animals to dissect the genetic basis of developmental traits, herein providing novel insights into the likely role of HOXB8 in feather development and differentiation.

Genetic variation is a key part for the study of evolution, development and differentiation. In domestic animals, many breeds display striking phenotypes that differentiate them from their wild ancestors. Several of these have been related to structural variations, including Fibromelanosis and Rose-comb in chickens, Double-muscled and Osteopetrosis in cattle, Cone degeneration in dogs, and White coat color in pigs. The feather is a type of skin appendages that exists in multiple variants on different body parts, and the derived feathering phenotypes in domestic birds are perfect resources to decipher the mechanisms regulating feather development and differentiation. Here we study the genetics of the Muffs and beard trait, a variant that alters the feather development in the facial area of chickens. We show that this phenotype is associated with a genomic structural variant that leads to an ectopic expression of HOXB8 in the facial skin during feather development. This is thus another example of how structural variants in the genome lead to novel, derived phenotypic changes in domestic animals and suggests an important role for HOXB8 in feather development.