Pitx1 broadly associates with limb enhancers and is enriched on hindlimb cis-regulatory elements

PITX1 as a grading, prognostic and tumor-infiltrating immune cells marker for chondrosarcoma: a public database-based immunoassay and tissue sample analysis

Background

Chondrosarcoma (CHS) is a rare bone cancer originating from chondrocytes, with high-grade cases associated with high mortality rates. However, the prognostic factors and therapeutic targets for CHS have not been studied.

Methods

Graded gene differential analysis was conducted on 97 CHS tissues to identify genes associated with CHS grading. Additionally, we performed GO and KEGG enrichment analyses of the differentially-expressed genes (DEGs), as well as GSEA analysis, differential expression analysis, survival analysis, and univariable and multifactorial COX analysis of paired-like homology structural domain transcription factor 1 (PITX1). Furthermore, our findings investigated the relationship between tumor-infiltrating immune cells (TICs) in CHS tumors using CIBERSORT to calculate proportions and differences. Our findings also explored the associations among gene expression patterns, survival prognosis, TICs, and immune checkpoints across various cancer types. Finally, immunohistochemical staining was carried out on self-collected clinical samples to assess PITX1 expression levels and correlate them with clinical information.

Results

Gene differential expression analysis revealed a strong correlation between PITX1 expression and tumor grade. GO, KEGG enrichment, and GSEA analysis demonstrated the association of PITX1 with cell proliferation-related processes, such as cell cycle regulation and mitosis, and differentiation-related processes, such as RNA processing. PITX1 expression was associated with tumor stage and survival outcomes. Immunoassay indicated a positive correlation between PITX1 levels and TICs, immune checkpoints, and graded TICs. Pan-cancer analysis confirmed the differential expression of the PITX1 gene across multiple cancers, impacting survival prognosis, TIC patterns, and immune checkpoint regulation. Lastly, our 75 collection of clinical patient tissue samples exhibited varying levels of PITX1 expression across different cancer grades while also demonstrating a significant association with tumor differentiation and metastasis.

Conclusion

PITX1 is a novel biomarker for distinguishing between high-grade and low-grade CHS, serving as a prognostic indicator for patients with this condition and presenting a promising target for immunotherapy. These findings offer innovative insights into the treatment of CHS.

Functional evidence supports the potential role of Tbx4-HLEA in the hindlimb degeneration of cetaceans

The evolution of limb morphology plays an important role in animal adaptation to different ecological niches. To fully adapt to aquatic life, cetaceans underwent hindlimb degeneration and forelimb transformed into flipper; however, the molecular mechanisms underlying the limb changes in cetaceans remain unclear. We previous study had shown that the Tbx4 hindlimb enhancer A (Tbx4-HLEA) in cetaceans exhibited specific deletions and nucleotide substitutions, with significantly reduced regulatory activity. To further investigate whether cetacean HLEA has a potential impact on hindlimb development in vivo, a knock-in mouse model was generated by knocking in the homologous cetacean HLEA in the present study. Phenotypic analysis showed a significant reduction in hindlimb bud development in homozygous knock-in mice at embryonic day (E)10.5; however, the phenotypic difference was rescued after E11.5. Transcriptomic and epigenetic analyses indicated that the cetacean HLEA acts as an enhancer in the mouse embryos and significantly reduces the transcriptional expression levels of Tbx4 at E10.5, supporting that downregulation of cetaceans HLEA regulatory activity reduces the expression of Tbx4. Additionally, both the number of activated non-coding elements and chromatin accessibility near Tbx4 were increased in homozygous knock-in mice at E11.5. The functional redundancy of enhancers compensated for the functional defect of cetacean HLEA, rescuing the expression level of Tbx4, and may account for the phenotype restoration after E11.5. In conclusion, our study suggested that the evolution of cetacean HLEA may be an important link with relevant molecular mechanism for the hindlimb degeneration.

Developmental regulation of conserved non-coding element evolution provides insights into limb loss in squamates

Limb loss shows recurrent phenotypic evolution across squamate lineages. Here, based on three de novo-assembled genomes of limbless lizards from different lineages, we showed that divergence of conserved non-coding elements (CNEs) played an important role in limb development. These CNEs were associated with genes required for limb initiation and outgrowth, and with regulatory signals in the early stage of limb development. Importantly, we identified the extensive existence of insertions and deletions (InDels) in the CNEs, with the numbers ranging from 111 to 756. Most of these CNEs with InDels were lineage-specific in the limbless squamates. Nearby genes of these InDel CNEs were important to early limb formation, such as Tbx4, Fgf10, and Gli3. Based on functional experiments, we found that nucleotide mutations and InDels both affected the regulatory function of the CNEs. Our study provides molecular evidence underlying limb loss in squamate reptiles from a developmental perspective and sheds light on the importance of regulatory element InDels in phenotypic evolution.

Identification of candidate enhancers controlling the transcriptome during the formation of interphalangeal joints

The formation of the synovial joint begins with the visible emergence of a stripe of densely packed mesenchymal cells located between distal ends of the developing skeletal anlagen called the interzone. Recently the transcriptome of the early synovial joint was reported. Knowledge about enhancers would complement these data and lead to a better understanding of the control of gene transcription at the onset of joint development. Using ChIP-sequencing we have mapped the H3-signatures H3K27ac and H3K4me1 to locate regulatory elements specific for the interzone and adjacent phalange, respectively. This one-stage atlas of candidate enhancers (CEs) was used to map the association between these respective joint tissue specific CEs and biological processes. Subsequently, integrative analysis of transcriptomic data and CEs identified new putative regulatory elements of genes expressed in interzone (e.g., GDF5, BMP2 and DACT2) and phalange (e.g., MATN1, HAPLN1 and SNAI1). We also linked such CEs to genes known as crucial in synovial joint hypermobility and osteoarthritis, as well as phalange malformations. These analyses show that the CE atlas can serve as resource for identifying, and as starting point for experimentally validating, putative disease-causing genomic regulatory regions in patients with synovial joint dysfunctions and/or phalange disorders, and enhancer-controlled synovial joint and phalange formation.

Cell-specific alterations in Pitx1 regulatory landscape activation caused by the loss of a single enhancer

Developmental genes are frequently controlled by multiple enhancers sharing similar specificities. As a result, deletions of such regulatory elements have often failed to reveal their full function. Here, we use the Pitx1 testbed locus to characterize in detail the regulatory and cellular identity alterations following the deletion of one of its enhancers (Pen). By combining single cell transcriptomics and an in-embryo cell tracing approach, we observe an increased fraction of Pitx1 non/low-expressing cells and a decreased fraction of Pitx1 high-expressing cells. We find that the over-representation of Pitx1 non/low-expressing cells originates from a failure of the Pitx1 locus to coordinate enhancer activities and 3D chromatin changes. This locus mis-activation induces a localized heterochrony and a concurrent loss of irregular connective tissue, eventually leading to a clubfoot phenotype. This data suggests that, in some cases, redundant enhancers may be used to locally enforce a robust activation of their host regulatory landscapes.

(Non)Parallel developmental mechanisms in vertebrate appendage reduction and loss

Appendages have been reduced or lost hundreds of times during vertebrate evolution. This phenotypic convergence may be underlain by shared or different molecular mechanisms in distantly related vertebrate clades. To investigate, we reviewed the developmental and evolutionary literature of appendage reduction and loss in more than a dozen vertebrate genera from fish to mammals. We found that appendage reduction and loss was nearly always driven by modified gene expression as opposed to changes in coding sequences. Moreover, expression of the same genes was repeatedly modified across vertebrate taxa. However, the specific mechanisms by which expression was modified were rarely shared. The multiple routes to appendage reduction and loss suggest that adaptive loss of function phenotypes might arise routinely through changes in expression of key developmental genes.

Control of osteocyte dendrite formation by Sp7 and its target gene osteocrin

Some osteoblasts embed within bone matrix, change shape, and become dendrite-bearing osteocytes. The circuitry that drives dendrite formation during "osteocytogenesis" is poorly understood. Here we show that deletion of Sp7 in osteoblasts and osteocytes causes defects in osteocyte dendrites. Profiling of Sp7 target genes and binding sites reveals unexpected repurposing of this transcription factor to drive dendrite formation. Osteocrin is a Sp7 target gene that promotes osteocyte dendrite formation and rescues defects in Sp7-deficient mice. Single-cell RNA-sequencing demonstrates defects in osteocyte maturation in the absence of Sp7. Sp7-dependent osteocyte gene networks are associated with human skeletal diseases. Moreover, humans with a SP7R316C mutation show defective osteocyte morphology. Sp7-dependent genes that mark osteocytes are enriched in neurons, highlighting shared features between osteocytic and neuronal connectivity. These findings reveal a role for Sp7 and its target gene Osteocrin in osteocytogenesis, revealing that pathways that control osteocyte development influence human bone diseases.

Tbx4 function during hindlimb development reveals a mechanism that explains the origins of proximal limb defects

We dissect genetically a gene regulatory network that involves the transcription factors Tbx4, Pitx1 and Isl1 acting cooperatively to establish the hindlimb bud, and identify key differences in the pathways that initiate formation of the hindlimb and forelimb. Using live image analysis of murine limb mesenchyme cells undergoing chondrogenesis in micromass culture, we distinguish a series of changes in cellular behaviours and cohesiveness that are required for chondrogenic precursors to undergo differentiation. Furthermore, we provide evidence that the proximal hindlimb defects observed in Tbx4 mutant mice result from a failure in the early differentiation step of chondroprogenitors into chondrocytes, providing an explanation for the origins of proximally biased limb defects.

Understanding and Engineering Chromatin as a Dynamical System across Length and Timescales

Connecting the molecular structure and function of chromatin across length and timescales remains a grand challenge to understanding and engineering cellular behaviors. Across five orders of magnitude, dynamic processes constantly reshape chromatin structures, driving spaciotemporal patterns of gene expression and cell fate. Through the interplay of structure and function, the genome operates as a highly dynamic feedback control system. Recent experimental techniques have provided increasingly detailed data that revise and augment the relatively static, hierarchical view of genomic architecture with an understanding of how dynamic processes drive organization. Here, we review how novel technologies from sequencing, imaging, and synthetic biology refine our understanding of chromatin structure and function and enable chromatin engineering. Finally, we discuss opportunities to use these tools to enhance understanding of the dynamic interrelationship of chromatin structure and function.

The transcription factor PITX1 drives astrocyte differentiation by regulating the SOX9 gene

Astrocytes perform multiple essential functions in the developing and mature brain, including regulation of synapse formation, control of neurotransmitter release and uptake, and maintenance of extracellular ion balance. As a result, astrocytes have been implicated in the progression of neurodegenerative disorders such as Alzheimer's disease, Huntington's disease, and Parkinson's disease. Despite these critical functions, the study of human astrocytes can be difficult because standard differentiation protocols are time-consuming and technically challenging, but a differentiation protocol recently developed in our laboratory enables the efficient derivation of astrocytes from human embryonic stem cells. We used this protocol along with microarrays, luciferase assays, electrophoretic mobility shift assays, and ChIP assays to explore the genes involved in astrocyte differentiation. We demonstrate that paired-like homeodomain transcription factor 1 (PITX1) is critical for astrocyte differentiation. PITX1 overexpression induced early differentiation of astrocytes, and its knockdown blocked astrocyte differentiation. PITX1 overexpression also increased and PITX1 knockdown decreased expression of sex-determining region Y box 9 (SOX9), known initiator of gliogenesis, during early astrocyte differentiation. Moreover, we determined that PITX1 activates the SOX9 promoter through a unique binding motif. Taken together, these findings indicate that PITX1 drives astrocyte differentiation by sustaining activation of the SOX9 promoter.

Network-regulated organ allometry: The developmental regulation of morphological scaling

Morphological scaling relationships, or allometries, describe how traits grow coordinately and covary among individuals in a population. The developmental regulation of scaling is essential to generate correctly proportioned adults across a range of body sizes, while the mis-regulation of scaling may result in congenital birth defects. Research over several decades has identified the developmental mechanisms that regulate the size of individual traits. Nevertheless, we still have poor understanding of how these mechanisms work together to generate correlated size variation among traits in response to environmental and genetic variation. Conceptually, morphological scaling can be generated by size-regulatory factors that act directly on multiple growing traits (trait-autonomous scaling), or indirectly via hormones produced by central endocrine organs (systemically regulated scaling), and there are a number of well-established examples of such mechanisms. There is much less evidence, however, that genetic and environmental variation actually acts on these mechanisms to generate morphological scaling in natural populations. More recent studies indicate that growing organs can themselves regulate the growth of other organs in the body. This suggests that covariation in trait size can be generated by network-regulated scaling mechanisms that respond to changes in the growth of individual traits. Testing this hypothesis, and one of the main challenges of understanding morphological scaling, requires connecting mechanisms elucidated in the laboratory with patterns of scaling observed in the natural world. This article is categorized under: Establishment of Spatial and Temporal Patterns > Regulation of Size, Proportion, and Timing Comparative Development and Evolution > Organ System Comparisons Between Species.

Evolutionary Selection and Constraint on Human Knee Chondrocyte Regulation Impacts Osteoarthritis Risk

During human evolution, the knee adapted to the biomechanical demands of bipedalism by altering chondrocyte developmental programs. This adaptive process was likely not without deleterious consequences to health. Today, osteoarthritis occurs in 250 million people, with risk variants enriched in non-coding sequences near chondrocyte genes, loci that likely became optimized during knee evolution. We explore this relationship by epigenetically profiling joint chondrocytes, revealing ancient selection and recent constraint and drift on knee regulatory elements, which also overlap osteoarthritis variants that contribute to disease heritability by tending to modify constrained functional sequence. We propose a model whereby genetic violations to regulatory constraint, tolerated during knee development, lead to adult pathology. In support, we discover a causal enhancer variant (rs6060369) present in billions of people at a risk locus (GDF5-UQCC1), showing how it impacts mouse knee-shape and osteoarthritis. Overall, our methods link an evolutionarily novel aspect of human anatomy to its pathogenesis.

Developmental and Evolutionary Allometry of the Mammalian Limb Skeleton

The variety of limb skeletal proportions enables a remarkable diversity of behaviors that include powered flight in bats and flipper-propelled swimming in whales using extremes of a range of homologous limb architectures. Even within human limbs, bone lengths span more than an order of magnitude from the short finger and toe bones to the long arm and leg bones. Yet all of this diversity arises from embryonic skeletal elements that are each a very similar size at formation. In this review article, I survey what is and is not yet known of the development and evolution of skeletal proportion at multiple hierarchical levels of biological organization. These include the cellular parameters of skeletal elongation in the cartilage growth plate, genes associated with differential growth, and putative gene regulatory mechanisms that would allow both covariant and independent evolution of the forelimbs and hindlimbs and of individual limb segments. Although the genetic mechanisms that shape skeletal proportion are still largely unknown, and most of what is known is limited to mammals, it is becoming increasingly apparent that the diversity of bone lengths is an emergent property of a complex system that controls elongation of individual skeletal elements using a genetic toolkit shared by all.

Variation in mouse pelvic morphology maps to locations enriched in Sox9 Class II and Pitx1 regulatory features

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high-resolution quantitative trait locus (QTL) mapping in the F₃₄ generation of an advanced intercross experiment (LG,SM-G₃₄ ) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome- or chromosome-wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation-seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

Identifying the homology of the short human pisiform and its lost ossification center

Background

The pisiform and calcaneus are paralogous bones of the wrist and ankle and are the only carpal and tarsal, respectively, to develop from two ossification centers with an associated growth plate in mammals. Human pisiforms and calcanei have undergone drastic evolutionary changes since our last common ancestor with chimpanzees and bonobos. The human pisiform is truncated and has lost an ossification center with the associated growth plate, while the human calcaneus has expanded and retained two ossification centers and a growth plate. Mammalian pisiforms represent a wide range of morphologies but extremely short pisiforms are rare and ossification center loss is even rarer. This raises the question of whether the sole human pisiform ossification center is homologous to the primary center or the secondary center of other species. We performed an ontogenetic study of pisiform and calcaneus ossification patterns and timing in macaques, apes, and humans (n = 907) from museum skeletal collections to address this question.

Results

Human pisiforms ossify irregularly and lack characteristic features of other primates while they develop. Pisiform primary and secondary center ossification timing typically matches that of the calcaneus of non-human primates, while the human pisiform corresponds with calcaneal secondary center ossification. Finally, human pisiforms ossify at the same dental stages as pisiform and calcaneal secondary centers in other hominoids.

Conclusions

These data indicate that the human pisiform is homologous to the pisiform epiphysis of other species, and that humans have lost a primary ossification center and associated growth plate while retaining ossification timing of the secondary center. This represents an exceptional evolutionary event and demonstrates a profound developmental change in the human wrist that is unusual not only among primates, but among mammals.

Genome Sequencing of the Japanese Eel (Anguilla japonica) for Comparative Genomic Studies on tbx4 and a tbx4 Gene Cluster in Teleost Fishes

Limbs originated from paired fish fins are an important innovation in Gnathostomata. Many studies have focused on limb development-related genes, of which the T-box transcription factor 4 gene (tbx4) has been considered as one of the most essential factors in the regulation of the hindlimb development. We previously confirmed pelvic fin loss in tbx4-knockout zebrafish. Here, we report a high-quality genome assembly of the Japanese eel (Anguilla japonica), which is an economically important fish without pelvic fins. The assembled genome is 1.13 Gb in size, with a scaffold N50 of 1.03 Mb. In addition, we collected 24 tbx4 sequences from 22 teleost fishes to explore the correlation between tbx4 and pelvic fin evolution. However, we observed complete exon structures of tbx4 in several pelvic-fin-loss species such as Ocean sunfish (Mola mola) and ricefield eel (Monopterus albus). More interestingly, an inversion of a special tbx4 gene cluster (brip1-tbx4-tbx2b- bcas3) occurred twice independently, which coincides with the presence of fin spines. A nonsynonymous mutation (M82L) was identified in the nuclear localization sequence (NLS) of the Japanese eel tbx4. We also examined variation and loss of hindlimb enhancer B (HLEB), which may account for pelvic fin loss in Tetraodontidae and Diodontidae. In summary, we generated a genome assembly of the Japanese eel, which provides a valuable genomic resource to study the evolution of fish tbx4 and helps elucidate the mechanism of pelvic fin loss in teleost fishes. Our comparative genomic studies, revealed for the first time a potential correlation between the tbx4 gene cluster and the evolutionary development of toxic fin spines. Because fin spines in teleosts are usually venoms, this tbx4 gene cluster may facilitate the genetic engineering of toxin-related marine drugs.

Pigeon foot feathering reveals conserved limb identity networks

The tetrapod limb is a stunning example of evolutionary diversity, with dramatic variation not only among distantly related species, but also between the serially homologous forelimbs (FLs) and hindlimbs (HLs) within species. Despite this variation, highly conserved genetic and developmental programs underlie limb development and identity in all tetrapods, raising the question of how limb diversification is generated from a conserved toolkit. In some breeds of domestic pigeon, shifts in the expression of two conserved limb identity transcription factors, PITX1 and TBX5, are associated with the formation of feathered HLs with partial FL identity. To determine how modulation of PITX1 and TBX5 expression affects downstream gene expression, we compared the transcriptomes of embryonic limb buds from pigeons with scaled and feathered HLs. We identified a set of differentially expressed genes enriched for genes encoding transcription factors, extracellular matrix proteins, and components of developmental signaling pathways with important roles in limb development. A subset of the genes that distinguish scaled and feathered HLs are also differentially expressed between FL and scaled HL buds in pigeons, pinpointing a set of gene expression changes downstream of PITX1 and TBX5 in the partial transformation from HL to FL identity. We extended our analyses by comparing pigeon limb bud transcriptomes to chicken, anole lizard, and mammalian datasets to identify deeply conserved PITX1- and TBX5-responsive components of the limb identity program. Our analyses reveal a suite of predominantly low-level gene expression changes that are conserved across amniotes to regulate the identity of morphologically distinct limbs.

Silenced PITX1 promotes chemotherapeutic resistance to 5-fluorocytosine and cisplatin in gastric cancer cells

Resistance to chemotherapeutic drugs leads to a poor prognosis in gastric cancer (GC). The present study aimed to assess the association between pituitary homeobox paired homeodomain transcription 1 (PITX1) expression and the sensitivity of GC cells to the chemotherapeutic drugs 5-fluorouracil (5-FU) and cisplatin (CDDP). In the present study, the gastric cancer cell lines GES-1, AGS, BGC-823, MCG-803 and SGC-7901 were used. The expression of PITX1 was determined via reverse transcription-quantitative polymerase chain reaction in GC cell lines. AGS and BGC-823 cells, which exhibit a decreased PITX1 expression, were transfected with a PITX1 cDNA construct and its control vector. MCG-803 and SGC-7901 cells, which exhibit an increased PITX1 expression, were transfected with siRNA against PITX1 and its control scramble sequence. A Cell Counting kit-8 assay was performed to determine the impact of PITX1 expression on the sensitivity of GC cells to 5-FU and CDDP. The Cancer Genome Atlas database was used to analyze the expression of PITX1 with GC prognosis in the Asian population and to assess the potential mechanism of PITX1 in 5-FU and CDDP resistance. The results revealed that the overexpression of PIXT1 increased the sensitivity of GC cells to 5-FU/CDDP. The combination of 5-FU/CDDP and PITX1 overexpression also reduced the proliferation of GC cells. Additionally, PIXT1 knockdown decreased the sensitivity of GC cells to 5-FU/CDDP. TCGA data revealed that a lower expression of PITX1 is exhibited in Asian GC patients than in normal individuals. GC patients with a lower expression of PITX1 had a poor prognosis. The expression of PITX1 affected the sensitivity of GC cells to 5-FU/CDDP, indicating that PITX1 may increase the efficacy of treatment in GC patients.

Differential analysis of chromatin accessibility and histone modifications for predicting mouse developmental enhancers

Enhancers are distal cis-regulatory elements that modulate gene expression. They are depleted of nucleosomes and enriched in specific histone modifications; thus, calling DNase-seq and histone mark ChIP-seq peaks can predict enhancers. We evaluated nine peak-calling algorithms for predicting enhancers validated by transgenic mouse assays. DNase and H3K27ac peaks were consistently more predictive than H3K4me1/2/3 and H3K9ac peaks. DFilter and Hotspot2 were the best DNase peak callers, while HOMER, MUSIC, MACS2, DFilter and F-seq were the best H3K27ac peak callers. We observed that the differential DNase or H3K27ac signals between two distant tissues increased the area under the precision-recall curve (PR-AUC) of DNase peaks by 17.5-166.7% and that of H3K27ac peaks by 7.1-22.2%. We further improved this differential signal method using multiple contrast tissues. Evaluated using a blind test, the differential H3K27ac signal method substantially improved PR-AUC from 0.48 to 0.75 for predicting heart enhancers. We further validated our approach using postnatal retina and cerebral cortex enhancers identified by massively parallel reporter assays, and observed improvements for both tissues. In summary, we compared nine peak callers and devised a superior method for predicting tissue-specific mouse developmental enhancers by reranking the called peaks.

Phenotype loss is associated with widespread divergence of the gene regulatory landscape in evolution

Detecting the genomic changes underlying phenotypic changes between species is a main goal of evolutionary biology and genomics. Evolutionary theory predicts that changes in cis-regulatory elements are important for morphological changes. We combined genome sequencing, functional genomics and genome-wide comparative analyses to investigate regulatory elements in lineages that lost morphological traits. We first show that limb loss in snakes is associated with widespread divergence of limb regulatory elements. We next show that eye degeneration in subterranean mammals is associated with widespread divergence of eye regulatory elements. In both cases, sequence divergence results in an extensive loss of transcription factor binding sites. Importantly, diverged regulatory elements are associated with genes required for normal limb patterning or normal eye development and function, suggesting that regulatory divergence contributed to the loss of these phenotypes. Together, our results show that genome-wide decay of the phenotype-specific cis-regulatory landscape is a hallmark of lost morphological traits.

Cis-regulatory elements are important factors for morphological changes. Here, the authors show widespread divergence of limb and eye regulatory elements in limb loss in snakes and eye degeneration in subterranean mammals respectively.

Homeotic transformation of tails into limbs in anurans

Anuran tadpoles can regenerate their tails after amputation. However, they occasionally form ectopic limbs instead of the lost tail part after vitamin A treatment. This is regarded as an example of a homeotic transformation. In this phenomenon, the developmental fate of the tail blastema is apparently altered from that of a tail to that of limbs, indicating a realignment of positional information in the blastema. Morphological observations and analyses of the development of skeletal elements during the process suggest that positional information in the blastema is rewritten from tail to trunk specification under the influence of vitamin A, resulting in limb formation. Despite the extensive information gained from morphological observations, a comprehensive understanding of this phenomenon also requires molecular data. We review previous studies related to anuran homeotic transformation. The findings of these studies provide a basis for evaluating major hypotheses and identifying molecular data that should be prioritized in future studies. Finally, we argue that positional information for the tail blastema changes to that for a part of the trunk, leading to homeotic transformations. To suggest this hypothesis, we present published data that favor the rewriting of positional information.

Comparative Genomics Reveals Accelerated Evolution in Conserved Pathways during the Diversification of Anole Lizards

Squamates include all lizards and snakes, and display some of the most diverse and extreme morphological adaptations among vertebrates. However, compared with birds and mammals, relatively few resources exist for comparative genomic analyses of squamates, hampering efforts to understand the molecular bases of phenotypic diversification in such a speciose clade. In particular, the ∼400 species of anole lizard represent an extensive squamate radiation. Here, we sequence and assemble the draft genomes of three anole species-Anolis frenatus, Anolis auratus, and Anolis apletophallus-for comparison with the available reference genome of Anolis carolinensis. Comparative analyses reveal a rapid background rate of molecular evolution consistent with a model of punctuated equilibrium, and strong purifying selection on functional genomic elements in anoles. We find evidence for accelerated evolution in genes involved in behavior, sensory perception, and reproduction, as well as in genes regulating limb bud development and hindlimb specification. Morphometric analyses of anole fore and hindlimbs corroborated these findings. We detect signatures of positive selection across several genes related to the development and regulation of the forebrain, hormones, and the iguanian lizard dewlap, suggesting molecular changes underlying behavioral adaptations known to reinforce species boundaries were a key component in the diversification of anole lizards.

PITX1 promotes chondrogenesis and myogenesis in mouse hindlimbs through conserved regulatory targets

The PITX1 transcription factor is expressed during hindlimb development, where it plays a critical role in directing hindlimb growth and the specification of hindlimb morphology. While it is known that PITX1 regulates hindlimb formation, in part, through activation of the Tbx4 gene, other transcriptional targets remain to be elucidated. We have used a combination of ChIP-seq and RNA-seq to investigate enhancer regions and target genes that are directly regulated by PITX1 in embryonic mouse hindlimbs. In addition, we have analyzed PITX1 binding sites in hindlimbs of Anolis lizards to identify ancient PITX1 regulatory targets. We find that PITX1-bound regions in both mouse and Anolis hindlimbs are strongly associated with genes implicated in limb and skeletal system development. Gene expression analyses reveal a large number of misexpressed genes in the hindlimbs of Pitx1^-/- mouse embryos. By intersecting misexpressed genes with genes that have neighboring mouse PITX1 binding sites, we identified 440 candidate targets of PITX1. Of these candidates, 68 exhibit ultra-conserved PITX1 binding events that are shared between mouse and Anolis hindlimbs. Among the ancient targets of PITX1 are important regulators of cartilage and skeletal muscle development, including Sox9 and Six1. Our data suggest that PITX1 promotes chondrogenesis and myogenesis in the hindlimb by direct regulation of several key members of the cartilage and muscle transcriptional networks.

Expression of a hindlimb-determining factor Pitx1 in the forelimb of the lizard Pogona vitticeps during morphogenesis

With over 9000 species, squamates, which include lizards and snakes, are the largest group of reptiles and second-largest order of vertebrates, spanning a vast array of appendicular skeletal morphology. As such, they provide a promising system for examining developmental and molecular processes underlying limb morphology. Using the central bearded dragon (Pogona vitticeps) as the primary study model, we examined limb morphometry throughout embryonic development and characterized the expression of three known developmental genes (GHR, Pitx1 and Shh) from early embryonic stage through to hatchling stage via reverse transcription quantitative polymerase chain reaction (RT-qPCR) and immunohistochemistry (IHC). In this study, all genes were found to be transcribed in both the forelimbs and hindlimbs of P. vitticeps. While the highest level of GHR expression occurred at the hatchling stage, Pitx1 and Shh expression was greatest earlier during embryogenesis, which coincides with the onset of the differentiation between forelimb and hindlimb length. We compared our finding of Pitx1 expression—a hindlimb-determining gene—in the forelimbs of P. vitticeps to that in a closely related Australian agamid lizard, Ctenophorus pictus, where we found Pitx1 expression to be more highly expressed in the hindlimb compared with the forelimb during early and late morphogenesis—a result consistent with that found across other tetrapods. Expression of Pitx1 in forelimbs has only rarely been documented, including via in situ hybridization in a chicken and a frog. Our findings from both RT-qPCR and IHC indicate that further research across a wider range of tetrapods is needed to more fully understand evolutionary variation in molecular processes underlying limb morphology.

RNA polymerase II primes Polycomb-repressed developmental genes throughout terminal neuronal differentiation

Polycomb repression in mouse embryonic stem cells (ESCs) is tightly associated with promoter co‐occupancy of RNA polymerase II (RNAPII) which is thought to prime genes for activation during early development. However, it is unknown whether RNAPII poising is a general feature of Polycomb repression, or is lost during differentiation. Here, we map the genome‐wide occupancy of RNAPII and Polycomb from pluripotent ESCs to non‐dividing functional dopaminergic neurons. We find that poised RNAPII complexes are ubiquitously present at Polycomb‐repressed genes at all stages of neuronal differentiation. We observe both loss and acquisition of RNAPII and Polycomb at specific groups of genes reflecting their silencing or activation. Strikingly, RNAPII remains poised at transcription factor genes which are silenced in neurons through Polycomb repression, and have major roles in specifying other, non‐neuronal lineages. We conclude that RNAPII poising is intrinsically associated with Polycomb repression throughout differentiation. Our work suggests that the tight interplay between RNAPII poising and Polycomb repression not only instructs promoter state transitions, but also may enable promoter plasticity in differentiated cells.

Pitx1 directly modulates the core limb development program to implement hindlimb identity

Forelimbs (FLs) and hindlimbs (HLs) develop complex musculoskeletal structures that rely on the deployment of a conserved developmental program. Pitx1, a transcription factor gene with expression restricted to HL and absent from FL, plays an important role in generating HL features. The genomic mechanisms by which Pitx1 effects HL identity remain poorly understood. Here, we use expression profiling and analysis of direct Pitx1 targets to characterize the HL- and FL-restricted genetic programs in mouse and situate the Pitx1-dependent gene network within the context of limb-specific gene regulation. We show that Pitx1 is a crucial component of a narrow network of HL-restricted regulators, acting on a developmental program that is shared between FL and HL. Pitx1 targets sites that are in a similar chromatin state in FL and HL and controls expression of patterning genes as well as the chondrogenic program, consistent with impaired chondrogenesis in Pitx1^-/- HL. These findings support a model in which multifactorial actions of a limited number of HL regulators redirect the generic limb development program in order to generate the unique structural features of the limb.

Limb development: a paradigm of gene regulation

The limb is a commonly used model system for developmental biology. Given the need for precise control of complex signalling pathways to achieve proper patterning, the limb is also becoming a model system for gene regulation studies. Recent developments in genomic technologies have enabled the genome-wide identification of regulatory elements that control limb development, yielding insights into the determination of limb morphology and forelimb versus hindlimb identity. The modulation of regulatory interactions - for example, through the modification of regulatory sequences or chromatin architecture - can lead to morphological evolution, acquired regeneration capacity or limb malformations in diverse species, including humans.

Transcriptional enhancers: Transcription, function and flexibility

Active transcriptional enhancers are often transcribed to eRNAs, whose changing levels mirror those of the target gene mRNA. We discuss some of the reported functions of these eRNAs and their likely diversity to allow utilization of distinct cis regulatory regions to enhance transcription in diverse developmental and cellular contexts.

A unique stylopod patterning mechanism by Shox2-controlled osteogenesis

Vertebrate appendage patterning is programmed by Hox-TALE factor-bound regulatory elements. However, it remains unclear which cell lineages are commissioned by Hox-TALE factors to generate regional specific patterns and whether other Hox-TALE co-factors exist. In this study, we investigated the transcriptional mechanisms controlled by the Shox2 transcriptional regulator in limb patterning. Harnessing an osteogenic lineage-specific Shox2 inactivation approach we show that despite widespread Shox2 expression in multiple cell lineages, lack of the stylopod observed upon Shox2 deficiency is a specific result of Shox2 loss of function in the osteogenic lineage. ChIP-Seq revealed robust interaction of Shox2 with cis-regulatory enhancers clustering around skeletogenic genes that are also bound by Hox-TALE factors, supporting a lineage autonomous function of Shox2 in osteogenic lineage fate determination and skeleton patterning. Pbx ChIP-Seq further allowed the genome-wide identification of cis-regulatory modules exhibiting co-occupancy of Pbx, Meis and Shox2 transcriptional regulators. Integrative analysis of ChIP-Seq and RNA-Seq data and transgenic enhancer assays indicate that Shox2 patterns the stylopod as a repressor via interaction with enhancers active in the proximal limb mesenchyme and antagonizes the repressive function of TALE factors in osteogenesis.

Generation and gene expression profiling of 48 transcription-factor-inducible mouse embryonic stem cell lines

Mouse embryonic stem cells (ESCs) can differentiate into a wide range – and possibly all cell types in vitro, and thus provide an ideal platform to study systematically the action of transcription factors (TFs) in cell differentiation. Previously, we have generated and analyzed 137 TF-inducible mouse ESC lines. As an extension of this “NIA Mouse ESC Bank,” we generated and characterized 48 additional mouse ESC lines, in which single TFs in each line could be induced in a doxycycline-controllable manner. Together, with the previous ESC lines, the bank now comprises 185 TF-manipulable ESC lines (>10% of all mouse TFs). Global gene expression (transcriptome) profiling revealed that the induction of individual TFs in mouse ESCs for 48 hours shifts their transcriptomes toward specific differentiation fates (e.g., neural lineages by Myt1 Isl1, and St18; mesodermal lineages by Pitx1, Pitx2, Barhl2, and Lmx1a; white blood cells by Myb, Etv2, and Tbx6, and ovary by Pitx1, Pitx2, and Dmrtc2). These data also provide and lists of inferred target genes of each TF and possible functions of these TFs. The results demonstrate the utility of mouse ESC lines and their transcriptome data for understanding the mechanism of cell differentiation and the function of TFs.

Molecular shifts in limb identity underlie development of feathered feet in two domestic avian species

Birds display remarkable diversity in the distribution and morphology of scales and feathers on their feet, yet the genetic and developmental mechanisms governing this diversity remain unknown. Domestic pigeons have striking variation in foot feathering within a single species, providing a tractable model to investigate the molecular basis of skin appendage differences. We found that feathered feet in pigeons result from a partial transformation from hindlimb to forelimb identity mediated by cis-regulatory changes in the genes encoding the hindlimb-specific transcription factor Pitx1 and forelimb-specific transcription factor Tbx5. We also found that ectopic expression of Tbx5 is associated with foot feathers in chickens, suggesting similar molecular pathways underlie phenotypic convergence between these two species. These results show how changes in expression of regional patterning genes can generate localized changes in organ fate and morphology, and provide viable molecular mechanisms for diversity in hindlimb scale and feather distribution.

Deletions of 5' HOXC genes are associated with lower extremity malformations, including clubfoot and vertical talus

BACKGROUND

Deletions of the HOXC gene cluster result in variable phenotypes in mice, but have been rarely described in humans.

OBJECTIVE

To report chromosome 12q13.13 microdeletions ranging from 13 to 175 kb and involving the 5' HOXC genes in four families, segregating congenital lower limb malformations, including clubfoot, vertical talus and hip dysplasia.

METHODS

Probands (N=253) with clubfoot or vertical talus were screened for point mutations and copy number variants using multiplexed direct genomic selection, a pooled BAC targeted capture approach. SNP genotyping included 1178 probands with clubfoot or vertical talus and 1775 controls.

RESULTS

The microdeletions share a minimal non-coding region overlap upstream of HOXC13, with variable phenotypes depending upon HOXC13, HOXC12 or the HOTAIR lncRNA inclusion. SNP analysis revealed HOXC11 p.Ser191Phe segregating with clubfoot in a small family and enrichment of HOXC12 p.Asn176Lys in patients with clubfoot or vertical talus (rs189468720, p=0.0057, OR=3.8). Defects in limb morphogenesis include shortened and overlapping toes, as well as peroneus muscle hypoplasia. Finally, HOXC and HOXD gene expression is reduced in fibroblasts from a patient with a 5' HOXC deletion, consistent with previous studies demonstrating that dosage of lncRNAs alters expression of HOXD genes in trans.

CONCLUSIONS

Because HOXD10 has been implicated in the aetiology of congenital vertical talus, variation in its expression may contribute to the lower limb phenotypes occurring with 5' HOXC microdeletions. Identification of 5' HOXC microdeletions highlights the importance of transcriptional regulators in the aetiology of severe lower limb malformations and will improve their diagnosis and management.

Functionally conserved enhancers with divergent sequences in distant vertebrates

BACKGROUND

To examine the contributions of sequence and function conservation in the evolution of enhancers, we systematically identified enhancers whose sequences are not conserved among distant groups of vertebrate species, but have homologous function and are likely to be derived from a common ancestral sequence. Our approach combined comparative genomics and epigenomics to identify potential enhancer sequences in the genomes of three groups of distantly related vertebrate species.

RESULTS

We searched for sequences that were conserved within groups of closely related species but not between groups of more distant species, and were associated with an epigenetic mark of enhancer activity. To facilitate inferring orthology between non-conserved sequences, we limited our search to introns whose orthology could be unambiguously established by mapping the bracketing exons. We show that a subset of these non-conserved but syntenic sequences from the mouse and zebrafish genomes have homologous functions in a zebrafish transgenic enhancer assay. The conserved expression patterns driven by these enhancers are probably associated with short transcription factor-binding motifs present in the divergent sequences.

CONCLUSIONS

We have identified numerous potential enhancers with divergent sequences but a conserved function. These results indicate that selection on function, rather than sequence, may be a common mode of enhancer evolution; evidence for selection at the sequence level is not a necessary criterion to define a gene regulatory element.

Shared Enhancer Activity in the Limbs and Phallus and Functional Divergence of a Limb-Genital cis-Regulatory Element in Snakes

The amniote phallus and limbs differ dramatically in their morphologies but share patterns of signaling and gene expression in early development. Thus far, the extent to which genital and limb transcriptional networks also share cis-regulatory elements has remained unexplored. We show that many limb enhancers are retained in snake genomes, suggesting that these elements may function in non-limb tissues. Consistent with this, our analysis of cis-regulatory activity in mice and Anolis lizards reveals that patterns of enhancer activity in embryonic limbs and genitalia overlap heavily. In mice, deletion of HLEB, an enhancer of Tbx4, produces defects in hindlimbs and genitalia, establishing the importance of this limb-genital enhancer for development of these different appendages. Further analyses demonstrate that the HLEB of snakes has lost hindlimb enhancer function while retaining genital activity. Our findings identify roles for Tbx4 in genital development and highlight deep similarities in cis-regulatory activity between limbs and genitalia.

Strong expression of paired-like homeodomain transcription factor 1 (PITX1) is associated with a favorable outcome in human osteosarcoma

Paired-like homeodomain transcription factor 1 (PITX1) has been implicated as a tumor suppressor in various cancers. However, the biological and clinical significance of PITX1 in osteosarcoma has not been fully elucidated. Here, we studied the expression and clinical significance of PITX1 in 6 normal lower limb bone tissue specimens and 35 osteosarcoma tissue samples by immunohistochemistry. PITX1 was expressed in all normal tissues (6/6, 100 %) and in 85.7 % (30/35) of tumor tissues (P > 0.05). In addition, all normal tissue specimens showed high PITX1 expression (6/6, 100 %) while only 23.3 % (7/30) osteosarcoma tissue specimens had high PITX1 expression (P < 0.05). Patients with median overall survival (OS) >12 months had significantly higher PITX1 levels compared with those whose median OS was less than or equal to 12 months (P < 0.05 or 0.001). Furthermore, patients with lung metastasis had significantly lower PITX1 levels than patients without lung metastasis. In conclusion, PITX1 expression is downregulated in osteosarcoma and correlates with patient survival and lung metastasis.

How the embryo makes a limb: determination, polarity and identity

The vertebrate limb with its complex anatomy develops from a small bud of undifferentiated mesoderm cells encased in ectoderm. The bud has its own intrinsic polarity and can develop autonomously into a limb without reference to the rest of the embryo. In this review, recent advances are integrated with classical embryology, carried out mainly in chick embryos, to present an overview of how the embryo makes a limb bud. We will focus on how mesoderm cells in precise locations in the embryo become determined to form a limb and express the key transcription factors Tbx4 (leg/hindlimb) or Tbx5 (wing/forelimb). These Tbx transcription factors have equivalent functions in the control of bud formation by initiating a signalling cascade involving Wnts and fibroblast growth factors (FGFs) and by regulating recruitment of mesenchymal cells from the coelomic epithelium into the bud. The mesoderm that will form limb buds and the polarity of the buds is determined with respect to both antero-posterior and dorso-ventral axes of the body. The position in which a bud develops along the antero-posterior axis of the body will also determine its identity - wing/forelimb or leg/hindlimb. Hox gene activity, under the influence of retinoic acid signalling, is directly linked with the initiation of Tbx5 gene expression in the region along the antero-posterior axis of the body that will form wings/forelimbs and determines antero-posterior polarity of the buds. In contrast, Tbx4 expression in the regions that will form legs/hindlimbs is regulated by the homeoprotein Pitx1 and there is no evidence that Hox genes determine antero-posterior polarity of the buds. Bone morphogenetic protein (BMP) signalling determines the region along the dorso-ventral axis of the body in which both wings/forelimbs and legs/hindlimbs develop and dorso-ventral polarity of the buds. The polarity of the buds leads to the establishment of signalling regions - the dorsal and ventral ectoderm, producing Wnts and BMPs, respectively, the apical ectodermal ridge producing fibroblast growth factors and the polarizing region, Sonic hedgehog (Shh). These signals are the same in both wings/forelimbs and legs/hindlimbs and control growth and pattern formation by providing the mesoderm cells of the limb bud as it develops with positional information. The precise anatomy of the limb depends on the mesoderm cells in the developing bud interpreting positional information according to their identity - determined by Pitx1 in hindlimbs - and genotype. The competence to form a limb extends along the entire antero-posterior axis of the trunk - with Hox gene activity inhibiting the formation of forelimbs in the interlimb region - and also along the dorso-ventral axis.

Spatiotemporal regulation of GLI target genes in the mammalian limb bud

GLI proteins convert Sonic hedgehog (Shh) signaling into a transcriptional output in a tissue-specific fashion. The Shh pathway has been extensively studied in the limb bud, where it helps regulate growth through a SHH-FGF feedback loop. However, the transcriptional response is still poorly understood. We addressed this by determining the gene expression patterns of approximately 200 candidate GLI-target genes and identified three discrete SHH-responsive expression domains. GLI-target genes expressed in the three domains are predominately regulated by derepression of GLI3 but have different temporal requirements for SHH. The GLI binding regions associated with these genes harbor both distinct and common DNA motifs. Given the potential for interaction between the SHH and FGF pathways, we also measured the response of GLI-target genes to inhibition of FGF signaling and found the majority were either unaffected or upregulated. These results provide the first characterization of the spatiotemporal response of a large group of GLI-target genes and lay the foundation for a systems-level understanding of the gene regulatory networks underlying SHH-mediated limb patterning.

Transcriptome profiling of esophageal squamous cell carcinoma reveals a long noncoding RNA acting as a tumor suppressor

Esophageal Squamous Cell Carcinoma (ESCC) is among the most common malignant cancers worldwide. In the past, extensive efforts have been made to characterize the involvement of protein-coding genes in ESCC tumorigenesis but few for long noncoding RNAs (lncRNAs). To investigate the transcriptome profile and functional relevance of lncRNAs, we performed an integrative analysis of a customized combined lncRNA-mRNA microarray and RNA-seq data on ESCCs and matched normal tissues. We identified numerous lncRNAs that were differentially expressed between the normal and tumor tissues, termed "ESCC-associated lncRNAs (ESCALs)", of which, the majority displayed restricted expression pattern. Also, a subset of ESCALs appeared to be associated with ESCC patient survival. Gene set enrichment analysis (GSEA) further suggested that over half of the ESCALs were positively-or negatively-associated with metastasis. Among these, we identified a novel nuclear-retained lncRNA, named Epist, which is generally highly expressed in esophagus, and which is down-regulated during ESCC progression. Epist over-expression and knockdown studies further suggest that Epist inhibits the metastasis, acting as a tumor suppressor in ESCC. Collectively, our analysis of the ESCC transcriptome identified the potential tumor suppressing lncRNA Epist, and provided a foundation for future efforts to identify functional lncRNAs for cancerous therapeutic targeting.

RNA transcribed from a distal enhancer is required for activating the chromatin at the promoter of the gonadotropin α-subunit gene

Since the discovery that many transcriptional enhancers are transcribed into long noncoding RNAs termed "enhancer RNAs" (eRNAs), their putative role in enhancer function has been debated. Very recent evidence has indicted that some eRNAs play a role in initiating or activating transcription, possibly by helping recruit and/or stabilize binding of the general transcription machinery to the proximal promoter of their target genes. The distal enhancer of the gonadotropin hormone α-subunit gene, chorionic gonadotropin alpha (Cga), is responsible for Cga cell-specific expression in gonadotropes and thyrotropes, and we show here that it encodes two bidirectional nonpolyadenylated RNAs whose levels are increased somewhat by exposure to gonadotropin-releasing hormone but are not necessarily linked to Cga transcriptional activity. Knockdown of the more distal eRNA led to a drop in Cga mRNA levels, initially without effect on the forward eRNA levels. With time, however, the repression on the Cga increased, and the forward eRNA levels were suppressed also. We demonstrate that the interaction of the enhancer with the promoter is lost after eRNA knockdown. Dramatic changes also were seen in the chromatin, with an increase in total histone H3 occupancy throughout this region and a virtual loss of histone H3 Lys 4 trimethylation at the promoter following the eRNA knockdown. Moreover, histone H3 Lys 27 (H3K27) acetylation, which was found at both enhancer and promoter in wild-type cells, appeared to have been replaced by H3K27 trimethylation at the enhancer. Thus, the Cga eRNA mediates the physical interaction between these genomic regions and determines the chromatin structure of the proximal promoter to allow gene expression.

An integrative computational approach for prioritization of genomic variants

An essential step in the discovery of molecular mechanisms contributing to disease phenotypes and efficient experimental planning is the development of weighted hypotheses that estimate the functional effects of sequence variants discovered by high-throughput genomics. With the increasing specialization of the bioinformatics resources, creating analytical workflows that seamlessly integrate data and bioinformatics tools developed by multiple groups becomes inevitable. Here we present a case study of a use of the distributed analytical environment integrating four complementary specialized resources, namely the Lynx platform, VISTA RViewer, the Developmental Brain Disorders Database (DBDB), and the RaptorX server, for the identification of high-confidence candidate genes contributing to pathogenesis of spina bifida. The analysis resulted in prediction and validation of deleterious mutations in the SLC19A placental transporter in mothers of the affected children that causes narrowing of the outlet channel and therefore leads to the reduced folate permeation rate. The described approach also enabled correct identification of several genes, previously shown to contribute to pathogenesis of spina bifida, and suggestion of additional genes for experimental validations. The study demonstrates that the seamless integration of bioinformatics resources enables fast and efficient prioritization and characterization of genomic factors and molecular networks contributing to the phenotypes of interest.

Divergence in a master variator generates distinct phenotypes and transcriptional responses

Genetic basis of phenotypic differences in individuals is an important area in biology and personalized medicine. Analysis of divergent Saccharomyces cerevisiae strains grown under different conditions revealed extensive variation in response to both drugs (e.g., 4-nitroquinoline 1-oxide [4NQO]) and different carbon sources. Differences in 4NQO resistance were due to amino acid variation in the transcription factor Yrr1. Yrr1(YJM789) conferred 4NQO resistance but caused slower growth on glycerol, and vice versa with Yrr1(S96), indicating that alleles of Yrr1 confer distinct phenotypes. The binding targets of Yrr1 alleles from diverse yeast strains varied considerably among different strains grown under the same conditions as well as for the same strain under different conditions, indicating that distinct molecular programs are conferred by the different Yrr1 alleles. Our results demonstrate that genetic variations in one important control gene (YRR1), lead to distinct regulatory programs and phenotypes in individuals. We term these polymorphic control genes "master variators."

Genome-wide data (ChIP-seq) enabled identification of cell wall-related and aquaporin genes as targets of tomato ASR1, a drought stress-responsive transcription factor

BACKGROUND

Identifying the target genes of transcription factors is important for unraveling regulatory networks in all types of organisms. Our interest was precisely to uncover the spectrum of loci regulated by a widespread plant transcription factor involved in physiological adaptation to drought, a type of stress that plants have encountered since the colonization of land habitats 400 MYA. The regulator under study, named ASR1, is exclusive to the plant kingdom (albeit absent in Arabidopsis) and known to alleviate the stress caused by restricted water availability. As its target genes are still unknown despite the original cloning of Asr1 cDNA 20 years ago, we examined the tomato genome for specific loci interacting in vivo with this conspicuous protein.

RESULTS

We performed ChIP followed by high throughput DNA sequencing (ChIP-seq) on leaves from stressed tomato plants, using a high-quality anti-ASR1 antibody. In this way, we unraveled a novel repertoire of target genes, some of which are clearly involved in the response to drought stress. Many of the ASR1-enriched genomic loci we found encode enzymes involved in cell wall synthesis and remodeling as well as channels implicated in water and solute flux, such as aquaporins. In addition, we were able to determine a robust consensus ASR1-binding DNA motif.

CONCLUSIONS

The finding of cell wall synthesis and aquaporin genes as targets of ASR1 is consistent with their suggested role in the physiological adaptation of plants to water loss. The results gain insight into the environmental stress-sensing pathways leading to plant tolerance of drought.