Identification and characterization of causal non-coding variants in human genomes is challenging and requires substantial experimental resources. A new study by Tehranchi et al. describes a cost-effective approach for accurate mapping of molecular quantitative trait loci (QTLs) from pooled samples, a powerful way to link disease-associated changes to molecular functions.

SMAD4 regulates gene expression in response to BMP and TGFβ signal transduction, and is required for diverse morphogenetic processes, but its target genes have remained largely elusive. Here, we identify the SMAD4 target genes in mouse limb buds using an epitope-tagged Smad4 allele for ChIP-seq analysis in combination with transcription profiling. This analysis shows that SMAD4 predominantly mediates BMP signal transduction during early limb bud development. Unexpectedly, the expression of cholesterol biosynthesis enzymes is precociously downregulated and intracellular cholesterol levels are reduced in Smad4-deficient limb bud mesenchymal progenitors. Most importantly, our analysis reveals a predominant function of SMAD4 in upregulating target genes in the anterior limb bud mesenchyme. Analysis of differentially expressed genes shared between Smad4- and Shh-deficient limb buds corroborates this function of SMAD4 and also reveals the repressive effect of SMAD4 on posterior genes that are upregulated in response to SHH signaling. This analysis uncovers opposing trans-regulatory inputs from SHH- and SMAD4-mediated BMP signal transduction on anterior and posterior gene expression during the digit patterning and outgrowth in early limb buds.

ABSTRACT

SMAD4 regulates gene expression in response to BMP and TGFβ signal transduction and is required for diverse morphogenetic processes, but its target genes have remained largely elusive. Here, we use an epitope-tagged Smad4 allele for ChIP-seq analysis together with transcriptome analysis of wild-type and mouse forelimb buds lacking Smad4 in the mesenchyme. This analysis identifies the SMAD4 target genes during establishment of the feedback signaling system and establishes that SMAD4 predominantly mediates BMP signal-transduction during early limb bud development. Unexpectedly, the initial analysis reveals that the expression of cholesterol biosynthesis enzymes is precociously down-regulated and intracellular cholesterol levels reduced in Smad4 -deficient limb bud mesenchymal progenitors. The SMAD4 target GRNs includes genes, whose expression in the anterior limb bud is up-regulated by interactions of SMAD4 complexes with enhancers active in the anterior mesenchyme. This reveals a predominant function of SMAD4 in up-regulating target gene expression in the anterior limb bud mesenchyme. Analysis of differentially expressed genes that are shared between Smad4 - and Shh -deficient limb buds corroborates the positive role of SMAD4 in transcriptional regulation of anterior genes and reveals a repressive effect on posterior genes that are positively regulated by SHH signaling. This analysis uncovers the overall opposing effects of SMAD4-mediated BMP and SHH signalling on transcriptional regulation during early limb bud development. In summary, this analysis indicates that during early digit patterning and limb bud outgrowth, the anterior/proximal and proximo/distal expression dynamics of co-regulated genes are controlled by distinct and contrasting trans-regulatory inputs from SHH and SMAD4-mediated BMP signal transduction.

The combinatorial expression of Hox genes along the body axes is a major determinant of cell fate and plays a pivotal role in generating the animal body plan. Loss of HOXA13 and HOXD13 transcription factors (HOX13) leads to digit agenesis in mice, but how HOX13 proteins regulate transcriptional outcomes and confer identity to the distal-most limb cells has remained elusive. Here, we report on the genome-wide profiling of HOXA13 and HOXD13 in vivo binding and changes of the transcriptome and chromatin state in the transition from the early to the late-distal limb developmental program, as well as in Hoxa13^-/-; Hoxd13^-/- limbs. Our results show that proper termination of the early limb transcriptional program and activation of the late-distal limb program are coordinated by the dual action of HOX13 on cis-regulatory modules.

The HAND2 transcriptional regulator controls cardiac development, and we uncover additional essential functions in the endothelial to mesenchymal transition (EMT) underlying cardiac cushion development in the atrioventricular canal (AVC). In Hand2-deficient mouse embryos, the EMT underlying AVC cardiac cushion formation is disrupted, and we combined ChIP-seq of embryonic hearts with transcriptome analysis of wild-type and mutants AVCs to identify the functionally relevant HAND2 target genes. The HAND2 target gene regulatory network (GRN) includes most genes with known functions in EMT processes and AVC cardiac cushion formation. One of these is Snai1, an EMT master regulator whose expression is lost from Hand2-deficient AVCs. Re-expression of Snai1 in mutant AVC explants partially restores this EMT and mesenchymal cell migration. Furthermore, the HAND2-interacting enhancers in the Snai1 genomic landscape are active in embryonic hearts and other Snai1-expressing tissues. These results show that HAND2 directly regulates the molecular cascades initiating AVC cardiac valve development.

ABSTRACT

During development cell fates are specified by tightly controlled gene expression programs. PBX TALE transcription factors control gene regulatory networks (GRN) that direct vertebrate tissue patterning and organ morphogenesis. How PBX1/2 proteins acquire context-specific functions, despite widespread embryonic expression of Pbx1/2 , remains elusive. In mouse limb buds, mesenchymal-specific loss of PBX1/2 or of the transcriptional regulator HAND2 results in similar phenotypes, suggesting that PBX1/2- and HAND2-dependent programs converge to control limb development. To investigate this scenario, we combined tissue-specific and temporally-controlled mutagenesis with multi-omics approaches using the murine hindlimb bud as a model. We reconstructed a GRN collaboratively directed by PBX1/2 and HAND2, demonstrating that Pbx1 - Hand2 genetically interact in vivo during hindlimb patterning, with PBX1 concomitantly acting as an upstream regulator of Hand2 . At organismal-level resolution the GRN is active within restricted subsets of posterior-proximal hindlimb mesenchymal cells, wherein Pbx1/2 and Hand2 are co-expressed with their target genes. Genome-wide profiling of PBX1 binding across multiple tissues further revealed that HAND2 selects a subset of PBX-bound regions to impart limb patterning functionality. This research elucidates mechanisms underlying limb bud-specific functions by PBX1/2, while informing general principles by which promiscuous transcription factors cooperate with select cofactors to instruct distinct developmental programs.

Epigenomic mapping of enhancer-associated chromatin modifications facilitates the genome-wide discovery of tissue-specific enhancers in vivo. However, reliance on single chromatin marks leads to high rates of false-positive predictions. More sophisticated, integrative methods have been described, but commonly suffer from limited accessibility to the resulting predictions and reduced biological interpretability. Here we present the Limb-Enhancer Genie (LEG), a collection of highly accurate, genome-wide predictions of enhancers in the developing limb, available through a user-friendly online interface. We predict limb enhancers using a combination of >50 published limb-specific datasets and clusters of evolutionarily conserved transcription factor binding sites, taking advantage of the patterns observed at previously in vivo validated elements. By combining different statistical models, our approach outperforms current state-of-the-art methods and provides interpretable measures of feature importance. Our results indicate that including a previously unappreciated score that quantifies tissue-specific nuclease accessibility significantly improves prediction performance. We demonstrate the utility of our approach through in vivo validation of newly predicted elements. Moreover, we describe general features that can guide the type of datasets to include when predicting tissue-specific enhancers genome-wide, while providing an accessible resource to the general biological community and facilitating the functional interpretation of genetic studies of limb malformations.

A lingering question in developmental biology has centered on how transcription factors with widespread distribution in vertebrate embryos can perform tissue-specific functions. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX TALE homeoproteins, viewed primarily as HOX cofactors, attain context-specific developmental roles despite ubiquitous presence in the embryo. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally controlled mutagenesis with multi-omics approaches, we reconstruct a gene regulatory network (GRN) at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 interactions in subsets of posterior hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 interacts with subsets of PBX-bound regions to regulate limb-specific GRNs. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with cofactors that display domain-restricted localization to instruct tissue-specific developmental programs.

A lingering question in developmental biology has centered on how PBX TALE transcription factors, viewed primarily as HOX co-factors, can confer functional specificity to spatially-restricted HOX proteins despite being ubiquitously distributed in vertebrate embryos. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX homeoproteins themselves attain tissue-specific developmental functions. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally-controlled mutagenesis with multi-omics approaches, we then reconstruct a GRN at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 within subsets of posterior-proximal hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 selects a subset of PBX-bound regions to impart limb patterning functionality. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with select cofactors to instruct distinct developmental programs.