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

EVOLUTIONARY CO-OPTION OF AN ANCESTRAL CLOACAL REGULATORY LANDSCAPE DURING THE EMERGENCE OF DIGITS AND GENITALS

The transition from fins to limbs has been a rich source of discussion for more than a century. One open and important issue is understanding how the mechanisms that pattern digits arose during vertebrate evolution. In this context, the analysis of Hox gene expression and functions to infer evolutionary scenarios has been a productive approach to explain the changes in organ formation, particularly in limbs. In tetrapods, the transcription of Hoxd genes in developing digits depends on a well-characterized set of enhancers forming a large regulatory landscape1,2. This control system has a syntenic counterpart in zebrafish, even though they lack bona fide digits, suggestive of deep homology3 between distal fin and limb developmental mechanisms. We tested the global function of this landscape to assess ancestry and source of limb and fin variation. In contrast to results in mice, we show here that the deletion of the homologous control region in zebrafish has a limited effect on the transcription of hoxd genes during fin development. However, it fully abrogates hoxd expression within the developing cloaca, an ancestral structure related to the mammalian urogenital sinus. We show that similar to the limb, Hoxd gene function in the urogenital sinus of the mouse also depends on enhancers located in this same genomic domain. Thus, we conclude that the current regulation underlying Hoxd gene expression in distal limbs was co-opted in tetrapods from a preexisting cloacal program. The orthologous chromatin domain in fishes may illustrate a rudimentary or partial step in this evolutionary co-option.

Developmental and evolutionary comparative analysis of a regulatory landscape in mouse and chicken

Modifications in gene regulation are driving forces in the evolution of organisms. Part of these changes involve cis-regulatory elements (CREs), which contact their target genes through higher-order chromatin structures. However, how such architectures and variations in CREs contribute to transcriptional evolvability remains elusive. We use Hoxd genes as a paradigm for the emergence of regulatory innovations, as many relevant enhancers are located in a regulatory landscape highly conserved in amniotes. Here, we analysed their regulation in murine vibrissae and chicken feather primordia, two skin appendages expressing different Hoxd gene subsets, and compared the regulation of these genes in these appendages with that in the elongation of the posterior trunk. In the two former structures, distinct subsets of Hoxd genes are contacted by different lineage-specific enhancers, probably as a result of using an ancestral chromatin topology as an evolutionary playground, whereas the gene regulation that occurs in the mouse and chicken embryonic trunk partially relies on conserved CREs. A high proportion of these non-coding sequences active in the trunk have functionally diverged between species, suggesting that transcriptional robustness is maintained, despite considerable divergence in enhancer sequences.

Summary: Analyses of the relationships between chromatin architecture and regulatory activities at the HoxD locus show that ancestral transcription patterns can be maintained while new regulations evolve.

HOXA1, a breast cancer oncogene

More than 25 years ago, the first literature records mentioned HOXA1 expression in human breast cancer. A few years later, HOXA1 was confirmed as a proper oncogene in mammary tissue. In the following two decades, molecular data about the mode of action of the HOXA1 protein, the factors contributing to activate and maintain HOXA1 gene expression and the identity of its target genes have accumulated and provide a wider view on the association of this transcription factor to breast oncogenesis. Large-scale transcriptomic data gathered from wide cohorts of patients further allowed refining the relationship between breast cancer type and HOXA1 expression. Several recent reports have reviewed the connection between cancer hallmarks and the biology of HOX genes in general. Here we take HOXA1 as a paradigm and propose an extensive overview of the molecular data centered on this oncoprotein, from what its expression modulators, to the interactors contributing to its oncogenic activities, and to the pathways and genes it controls. The data converge to an intricate picture that answers questions on the multi-modality of its oncogene activities, point towards better understanding of breast cancer aetiology and thereby provides an appraisal for treatment opportunities.

Intestinal Dysbiosis in the Infant and the Future of Lacto-Engineering to Shape the Developing Intestinal Microbiome

PURPOSE

The goal of this study was to review the role of human milk in shaping the infant intestinal microbiota and the potential of human milk bioactive molecules to reverse trends of increasing intestinal dysbiosis and dysbiosis-associated diseases.

METHODS

This narrative review was based on recent and historic literature.

FINDINGS

Human milk immunoglobulins, oligosaccharides, lactoferrin, lysozyme, milk fat globule membranes, and bile salt-stimulating lipase are complex multifunctional bioactive molecules that, among other important functions, shape the composition of the infant intestinal microbiota.

IMPLICATIONS

The co-evolution of human milk components and human milk-consuming commensal anaerobes many thousands of years ago resulted in a stable low-diversity infant microbiota. Over the past century, the introduction of antibiotics and modern hygiene practices plus changes in the care of newborns have led to significant alterations in the intestinal microbiota, with associated increases in risk of dysbiosis-associated disease. A better understanding of mechanisms by which human milk shapes the intestinal microbiota of the infant during a vulnerable period of development of the immune system is needed to alter the current trajectory and decrease intestinal dysbiosis and associated diseases.

Sequential in cis mutagenesis in vivo reveals various functions for CTCF sites at the mouse HoxD cluster

Mammalian Hox gene clusters contain a range of CTCF binding sites. In addition to their importance in organizing a TAD border, which isolates the most posterior genes from the rest of the cluster, the positions and orientations of these sites suggest that CTCF may be instrumental in the selection of various subsets of contiguous genes, which are targets of distinct remote enhancers located in the flanking regulatory landscapes. We examined this possibility by producing an allelic series of cumulative in cis mutations in these sites, up to the abrogation of CTCF binding in the five sites located on one side of the TAD border. In the most impactful alleles, the global chromatin architecture of the locus was modified, yet not drastically, illustrating that CTCF sites located on one side of a strong TAD border are sufficient to organize at least part of this insulation. Spatial colinearity in the expression of these genes along the major body axis was nevertheless maintained, despite abnormal expression boundaries. In contrast, strong effects were scored in the selection of target genes responding to particular enhancers, leading to the misregulation of Hoxd genes in specific structures. Altogether, while most enhancer-promoter interactions can occur in the absence of this series of CTCF sites, the binding of CTCF in the Hox cluster is required to properly transform a rather unprecise process into a highly discriminative mechanism of interactions, which is translated into various patterns of transcription accompanied by the distinctive chromatin topology found at this locus. Our allelic series also allowed us to reveal the distinct functional contributions for CTCF sites within this Hox cluster, some acting as insulator elements, others being necessary to anchor or stabilize enhancer-promoter interactions, and some doing both, whereas they all together contribute to the formation of a TAD border. This variety of tasks may explain the amazing evolutionary conservation in the distribution of these sites among paralogous Hox clusters or between various vertebrates.

Analysis of Polycerate Mutants Reveals the Evolutionary Co-option of HOXD1 for Horn Patterning in Bovidae

In the course of evolution, pecorans (i.e., higher ruminants) developed a remarkable diversity of osseous cranial appendages, collectively referred to as “headgear,” which likely share the same origin and genetic basis. However, the nature and function of the genetic determinants underlying their number and position remain elusive. Jacob and other rare populations of sheep and goats are characterized by polyceraty, the presence of more than two horns. Here, we characterize distinct POLYCERATE alleles in each species, both associated with defective HOXD1 function. We show that haploinsufficiency at this locus results in the splitting of horn bud primordia, likely following the abnormal extension of an initial morphogenetic field. These results highlight the key role played by this gene in headgear patterning and illustrate the evolutionary co-option of a gene involved in the early development of bilateria to properly fix the position and number of these distinctive organs of Bovidae.

Induction of a chromatin boundary in vivo upon insertion of a TAD border

Mammalian genomes are partitioned into sub-megabase to megabase-sized units of preferential interactions called topologically associating domains or TADs, which are likely important for the proper implementation of gene regulatory processes. These domains provide structural scaffolds for distant cis regulatory elements to interact with their target genes within the three-dimensional nuclear space and architectural proteins such as CTCF as well as the cohesin complex participate in the formation of the boundaries between them. However, the importance of the genomic context in providing a given DNA sequence the capacity to act as a boundary element remains to be fully investigated. To address this question, we randomly relocated a topological boundary functionally associated with the mouse HoxD gene cluster and show that it can indeed act similarly outside its initial genomic context. In particular, the relocated DNA segment recruited the required architectural proteins and induced a significant depletion of contacts between genomic regions located across the integration site. The host chromatin landscape was re-organized, with the splitting of the TAD wherein the boundary had integrated. These results provide evidence that topological boundaries can function independently of their site of origin, under physiological conditions during mouse development.

During development, enhancer sequences tightly regulate the spatio-temporal expression of target genes often located hundreds of kilobases away. This complex process is made possible by the folding of chromatin into domains, which are separated from one another by specific genomic regions referred to as boundaries. In order to understand whether such boundary sequences require their particular genomic contexts to achieve their isolating effect, we analyzed the impact of introducing one such boundary, taken from the HoxD locus, into a distinct topological domain. We show that this ectopic boundary splits the host domain into two sub-domains and affects the expression levels of a neighboring gene. We conclude that this sequence can work independently from its genomic context and thus carries all the information necessary to act as a boundary element.

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.

Chromatin topology and the timing of enhancer function at the HoxD locus

The HoxD gene cluster is critical for proper limb formation in tetrapods. In the emerging limb buds, different subgroups of Hoxd genes respond first to a proximal regulatory signal, then to a distal signal that organizes digits. These two regulations are exclusive from one another and emanate from two distinct topologically associating domains (TADs) flanking HoxD, both containing a range of appropriate enhancer sequences. The telomeric TAD (T-DOM) contains several enhancers active in presumptive forearm cells and is divided into two sub-TADs separated by a CTCF-rich boundary, which defines two regulatory submodules. To understand the importance of this particular regulatory topology to control Hoxd gene transcription in time and space, we either deleted or inverted this sub-TAD boundary, eliminated the CTCF binding sites, or inverted the entire T-DOM to exchange the respective positions of the two sub-TADs. The effects of such perturbations on the transcriptional regulation of Hoxd genes illustrate the requirement of this regulatory topology for the precise timing of gene activation. However, the spatial distribution of transcripts was eventually resumed, showing that the presence of enhancer sequences, rather than either their exact topology or a particular chromatin architecture, is the key factor. We also show that the affinity of enhancers to find their natural target genes can overcome the presence of both a strong TAD border and an unfavorable orientation of CTCF sites.

Mammalian-specific ectodermal enhancers control the expression of Hoxc genes in developing nails and hair follicles

Vertebrate Hox genes are critical for the establishment of structures during the development of the main body axis. Subsequently, they play important roles either in organizing secondary axial structures such as the appendages, or during homeostasis in postnatal stages and adulthood. Here, we set up to analyze their elusive function in the ectodermal compartment, using the mouse limb bud as a model. We report that the HoxC gene cluster was co-opted to be transcribed in the distal limb ectoderm, where it is activated following the rule of temporal colinearity. These ectodermal cells subsequently produce various keratinized organs such as nails or claws. Accordingly, deletion of the HoxC cluster led to mice lacking nails (anonychia), a condition stronger than the previously reported loss of function of Hoxc13, which is the causative gene of the ectodermal dysplasia 9 (ECTD9) in human patients. We further identified two mammalian-specific ectodermal enhancers located upstream of the HoxC gene cluster, which together regulate Hoxc gene expression in the hair and nail ectodermal organs. Deletion of these regulatory elements alone or in combination revealed a strong quantitative component in the regulation of Hoxc genes in the ectoderm, suggesting that these two enhancers may have evolved along with the mammalian taxon to provide the level of HOXC proteins necessary for the full development of hair and nail.

Three distinct mechanisms, Notch instructive, permissive, and independent, regulate the expression of two different pericardial genes to specify cardiac cell subtypes

The development of a complex organ involves the specification and differentiation of diverse cell types constituting that organ. Two major cell subtypes, contractile cardial cells (CCs) and nephrocytic pericardial cells (PCs), comprise the Drosophila heart. Binding sites for Suppressor of Hairless [Su(H)], an integral transcription factor in the Notch signaling pathway, are enriched in the enhancers of PC-specific genes. Here we show three distinct mechanisms regulating the expression of two different PC-specific genes, Holes in muscle (Him), and Zn finger homeodomain 1 (zfh1). Him transcription is activated in PCs in a permissive manner by Notch signaling: in the absence of Notch signaling, Su(H) forms a repressor complex with co-repressors and binds to the Him enhancer, repressing its transcription; upon alleviation of this repression by Notch signaling, Him transcription is activated. In contrast, zfh1 is transcribed by a Notch-instructive mechanism in most PCs, where mere alleviation of repression by preventing the binding of Su(H)-co-repressor complex is not sufficient to activate transcription. Our results suggest that upon activation of Notch signaling, the Notch intracellular domain associates with Su(H) to form an activator complex that binds to the zfh1 enhancer, and that this activator complex is necessary for bringing about zfh1 transcription in these PCs. Finally, a third, Notch-independent mechanism activates zfh1 transcription in the remaining, even skipped-expressing, PCs. Collectively, our data show how the same feature, enrichment of Su(H) binding sites in PC-specific gene enhancers, is utilized by two very distinct mechanisms, one permissive, the other instructive, to contribute to the same overall goal: the specification and differentiation of a cardiac cell subtype by activation of the pericardial gene program. Furthermore, our results demonstrate that the zfh1 enhancer drives expression in two different domains using distinct Notch-instructive and Notch-independent mechanisms.

Polycomb Repressive Complexes in Hox Gene Regulation: Silencing and Beyond: The Functional Dynamics of Polycomb Repressive Complexes in Hox Gene Regulation

The coordinated expression of the Hox gene family encoding transcription factors is critical for proper embryonic development and patterning. Major efforts have thus been dedicated to understanding mechanisms controlling Hox expression. In addition to the temporal and spatial sequential activation of Hox genes, proper embryonic development requires that Hox genes get differentially silenced in a cell-type specific manner as development proceeds. Factors contributing to Hox silencing include the polycomb repressive complexes (PRCs), which control gene expression through epigenetic modifications. This review focuses on PRC-dependent regulation of the Hox genes and is aimed at integrating the growing complexity of PRC functional properties in the context of Hox regulation. In particular, mechanisms underlying PRC binding dynamics as well as a series of studies that have revealed the impact of PRC on the 3D organization of the genome is discussed, which has a significant role on Hox regulation during development.

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.

Comparative Genomic Characterization of the Multimammate Mouse Mastomys coucha

Mastomys are the most widespread African rodent and carriers of various diseases such as the plague or Lassa virus. In addition, mastomys have rapidly gained a large number of mammary glands. Here, we generated a genome, variome, and transcriptomes for Mastomys coucha. As mastomys diverged at similar times from mouse and rat, we demonstrate their utility as a comparative genomic tool for these commonly used animal models. Furthermore, we identified over 500 mastomys accelerated regions, often residing near important mammary developmental genes or within their exons leading to protein sequence changes. Functional characterization of a noncoding mastomys accelerated region, located in the HoxD locus, showed enhancer activity in mouse developing mammary glands. Combined, our results provide genomic resources for mastomys and highlight their potential both as a comparative genomic tool and for the identification of mammary gland number determining factors.

Impact of genome architecture on the functional activation and repression of Hox regulatory landscapes

BACKGROUND

The spatial organization of the mammalian genome relies upon the formation of chromatin domains of various scales. At the level of gene regulation in cis, collections of enhancer sequences define large regulatory landscapes that usually match with the presence of topologically associating domains (TADs). These domains often contain ranges of enhancers displaying similar or related tissue specificity, suggesting that in some cases, such domains may act as coherent regulatory units, with a global on or off state. By using the HoxD gene cluster, which specifies the topology of the developing limbs via highly orchestrated regulation of gene expression, as a paradigm, we investigated how the arrangement of regulatory domains determines their activity and function.

RESULTS

Proximal and distal cells in the developing limb express different levels of Hoxd genes, regulated by flanking 3' and 5' TADs, respectively. We characterized the effect of large genomic rearrangements affecting these two TADs, including their fusion into a single chromatin domain. We show that, within a single hybrid TAD, the activation of both proximal and distal limb enhancers globally occurred as when both TADs are intact. However, the activity of the 3' TAD in distal cells is generally increased in the fused TAD, when compared to wild type where it is silenced. Also, target gene activity in distal cells depends on whether or not these genes had previously responded to proximal enhancers, which determines the presence or absence of H3K27me3 marks. We also show that the polycomb repressive complex 2 is mainly recruited at the Hox gene cluster and can extend its coverage to far-cis regulatory sequences as long as confined to the neighboring TAD structure.

CONCLUSIONS

We conclude that antagonistic limb proximal and distal enhancers can exert their specific effects when positioned into the same TAD and in the absence of their genuine target genes. We also conclude that removing these target genes reduced the coverage of a regulatory landscape by chromatin marks associated with silencing, which correlates with its prolonged activity in time.

Embryonic timing, axial stem cells, chromatin dynamics, and the Hox clock

Collinear regulation of Hox genes in space and time has been an outstanding question ever since the initial work of Ed Lewis in 1978. Here we discuss recent advances in our understanding of this phenomenon in relation to novel concepts associated with large-scale regulation and chromatin structure during the development of both axial and limb patterns. We further discuss how this sequential transcriptional activation marks embryonic stem cell-like axial progenitors in mammals and, consequently, how a temporal genetic system is further translated into spatial coordinates via the fate of these progenitors. In this context, we argue the benefit and necessity of implementing this unique mechanism as well as the difficulty in evolving an alternative strategy to deliver this critical positional information.

The HoxD cluster is a dynamic and resilient TAD boundary controlling the segregation of antagonistic regulatory landscapes

The mammalian HoxD cluster lies between two topologically associating domains (TADs) matching distinct enhancer-rich regulatory landscapes. During limb development, the telomeric TAD controls the early transcription of Hoxd genes in forearm cells, whereas the centromeric TAD subsequently regulates more posterior Hoxd genes in digit cells. Therefore, the TAD boundary prevents the terminal Hoxd13 gene from responding to forearm enhancers, thereby allowing proper limb patterning. To assess the nature and function of this CTCF-rich DNA region in embryos, we compared chromatin interaction profiles between proximal and distal limb bud cells isolated from mutant stocks where various parts of this boundary region were removed. The resulting progressive release in boundary effect triggered inter-TAD contacts, favored by the activity of the newly accessed enhancers. However, the boundary was highly resilient, and only a 400-kb deletion, including the whole-gene cluster, was eventually able to merge the neighboring TADs into a single structure. In this unified TAD, both proximal and distal limb enhancers nevertheless continued to work independently over a targeted transgenic reporter construct. We propose that the whole HoxD cluster is a dynamic TAD border and that the exact boundary position varies depending on both the transcriptional status and the developmental context.

Control of growth and gut maturation by HoxD genes and the associated lncRNA Haglr

During embryonic development, Hox genes participate in the building of a functional digestive system in metazoans, and genetic conditions involving these genes lead to important, sometimes lethal, growth retardation. Recently, this phenotype was obtained after deletion of Haglr, the Hoxd antisense growth-associated long noncoding RNA (lncRNA) located between Hoxd1 and Hoxd3 In this study, we have analyzed the function of Hoxd genes in delayed growth trajectories by looking at several nested targeted deficiencies of the mouse HoxD cluster. Mutant pups were severely stunted during the suckling period, but many recovered after weaning. After comparing seven distinct HoxD alleles, including CRISPR/Cas9 deletions involving Haglr, we identified Hoxd3 as the critical component for the gut to maintain milk-digestive competence. This essential function could be abrogated by the dominant-negative effect of HOXD10 as shown by a genetic rescue approach, thus further illustrating the importance of posterior prevalence in Hox gene function. A role for the lncRNA Haglr in the control of postnatal growth could not be corroborated.