A function for all posterior Hoxd genes during digit development?

BST2 promotes gastric cancer metastasis under the regulation of HOXD9 and PABPC1

Gastric cancer (GC) constitutes substantial cancer mortality worldwide. Several cancer types aberrantly express bone marrow stromal cell antigen 2 (BST2), yet its functional and underlying mechanisms in GC progression remain unknown. In our study, RNA sequencing data revealed that BST2 was transcriptionally activated by homeobox D9 (HOXD9). BST2 was significantly upregulated in GC tissues and promoted epithelial-mesenchymal transition and metastasis of GC. BST2 knockdown reversed HOXD9's oncogenic effect on GC metastasis. Moreover, BST2 messenger RNA stability could be enhanced by poly(A) binding protein cytoplasmic 1 (PABPC1) through the interaction between BST2 3'-UTR and PABPC1 in GC cells. PABPC1 promoted GC metastasis, which BST2 silencing attenuated in vitro and in vivo. In addition, positive correlations among HOXD9, BST2, and PABPC1 were established in clinical samples. Taken together, increased expression of BST2 induced by HOXD9 synergizing with PABPC1 promoted GC cell migration and invasion capacity.

Dbx2 regulation in limbs suggests interTAD sharing of enhancers

BACKGROUND

During tetrapod limb development, the HOXA13 and HOXD13 transcription factors are critical for the emergence and organization of the autopod, the most distal aspect where digits will develop. Since previous work had suggested that the Dbx2 gene is a target of these factors, we set up to analyze in detail this potential regulatory interaction.

RESULTS

We show that HOX13 proteins bind to mammalian-specific sequences at the vicinity of the Dbx2 locus that have enhancer activity in developing digits. However, the functional inactivation of the DBX2 protein did not elicit any particular phenotype related to Hox genes inactivation in digits, suggesting either redundant or compensatory mechanisms. We report that the neighboring Nell2 and Ano6 genes are also expressed in distal limb buds and are in part controlled by the same Dbx2 enhancers despite being localized into two different topologically associating domains (TADs) flanking the Dbx2 locus.

CONCLUSIONS

We conclude that Hoxa13 and Hoxd genes cooperatively activate Dbx2 expression in developing digits through binding to mammalian specific regulatory sequences in the Dbx2 neighborhood. Furthermore, these enhancers can overcome TAD boundaries in either direction to co-regulate a set of genes located in distinct chromatin domains.

Genomic regions associated with muscularity in beef cattle differ in five contrasting cattle breeds

BACKGROUND

Linear type traits, which reflect the muscular characteristics of an animal, could provide insight into how, in some cases, morphologically very different animals can yield the same carcass weight. Such variability may contribute to differences in the overall value of the carcass since primal cuts vary greatly in price; such variability may also hinder successful genome-based association studies. Therefore, the objective of our study was to identify genomic regions that are associated with five muscularity linear type traits and to determine if these significant regions are common across five different breeds. Analyses were carried out using linear mixed models on imputed whole-genome sequence data in each of the five breeds, separately. Then, the results of the within-breed analyses were used to conduct an across-breed meta-analysis per trait.

RESULTS

We identified many quantitative trait loci (QTL) that are located across the whole genome and associated with each trait in each breed. The only commonality among the breeds and traits was a large-effect pleiotropic QTL on BTA2 that contained the MSTN gene, which was associated with all traits in the Charolais and Limousin breeds. Other plausible candidate genes were identified for muscularity traits including PDE1A, PPP1R1C and multiple collagen and HOXD genes. In addition, associated (gene ontology) GO terms and KEGG pathways tended to differ between breeds and between traits especially in the numerically smaller populations of Angus, Hereford, and Simmental breeds. Most of the SNPs that were associated with any of the traits were intergenic or intronic SNPs located within regulatory regions of the genome.

CONCLUSIONS

The commonality between the Charolais and Limousin breeds indicates that the genetic architecture of the muscularity traits may be similar in these breeds due to their similar origins. Conversely, there were vast differences in the QTL associated with muscularity in Angus, Hereford, and Simmental. Knowledge of these differences in genetic architecture between breeds is useful to develop accurate genomic prediction equations that can operate effectively across breeds. Overall, the associated QTL differed according to trait, which suggests that breeding for a morphologically different (e.g. longer and wider versus shorter and smaller) more efficient animal may become possible in the future.

Heterogeneous combinatorial expression of Hoxd genes in single cells during limb development

BACKGROUND

Global analyses of gene expression during development reveal specific transcription patterns associated with the emergence of various cell types, tissues, and organs. These heterogeneous patterns are instrumental to ensure the proper formation of the different parts of our body, as shown by the phenotypic effects generated by functional genetic approaches. However, variations at the cellular level can be observed within each structure or organ. In the developing mammalian limbs, expression of Hox genes from the HoxD cluster is differentially controlled in space and time, in cells that will pattern the digits and the forearms. While the Hoxd genes broadly share a common regulatory landscape and large-scale analyses have suggested a homogenous Hox gene transcriptional program, it has not previously been clear whether Hoxd genes are expressed together at the same levels in the same cells.

RESULTS

We report a high degree of heterogeneity in the expression of the Hoxd11 and Hoxd13 genes. We analyzed single-limb bud cell transcriptomes and show that Hox genes are expressed in specific combinations that appear to match particular cell types. In cells giving rise to digits, we find that the expression of the five relevant Hoxd genes (Hoxd9 to Hoxd13) is unbalanced, despite their control by known global enhancers. We also report that specific combinatorial expression follows a pseudo-time sequence, which is established based on the transcriptional diversity of limb progenitors.

CONCLUSIONS

Our observations reveal the existence of distinct combinations of Hoxd genes at the single-cell level during limb development. In addition, we document that the increasing combinatorial expression of Hoxd genes in this developing structure is associated with specific transcriptional signatures and that these signatures illustrate a temporal progression in the differentiation of these cells.

Topological Domains, Metagenes, and the Emergence of Pleiotropic Regulations at Hox Loci

Hox gene clusters of jaw vertebrates display a tight genomic organization, which has no equivalent in any other bilateria genomes sequenced thus far. It was previously argued that such a topological consolidation toward a condensed, metagenic structure was due to the accumulation of long-range regulations flanking Hox loci, a phenomenon made possible by the successive genome duplications that occurred at the root of the vertebrate lineage, similar to gene neofunctionalization but applied to a coordinated multigenic system. Here, we propose that the emergence of such large vertebrate regulatory landscapes containing a range of global enhancers was greatly facilitated by the presence of topologically associating domains (TADs). These chromatin domains, mostly constitutive, may have been used as genomic niches where novel regulations could evolve due to both the preexistence of a structural backbone poised to integrate novel regulatory inputs, and a highly adaptive transcriptional readout. We propose a scenario for the coevolution of such TADs and the emergence of pleiotropy at ancestral vertebrate Hox loci.

Reorganisation of Hoxd regulatory landscapes during the evolution of a snake-like body plan

Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.

Selection on different genes with equivalent functions: the convergence story told by Hox genes along the evolution of aquatic mammalian lineages

Background

Convergent evolution has been a challenging topic for decades, being cetaceans, pinnipeds and sirenians textbook examples of three independent origins of equivalent phenotypes. These mammalian lineages acquired similar anatomical features correlated to an aquatic life, and remarkably differ from their terrestrial counterparts. Whether their molecular evolutionary history also involved similar genetic mechanisms underlying such morphological convergence nevertheless remained unknown. To test for the existence of convergent molecular signatures, we studied the molecular evolution of Hox genes in these three aquatic mammalian lineages, comparing their patterns to terrestrial mammals. Hox genes are transcription factors that play a pivotal role in specifying embryonic regional identity of nearly any bilateral animal, and are recognized major agents for diversification of body plans.

Results

We detected few signatures of positive selection on Hox genes across the three aquatic mammalian lineages and verified that purifying selection prevails in these sequences, as expected for pleiotropic genes. Genes found as being positively selected differ across the aquatic mammalian lineages, but we identified a substantial overlap of their developmental functions. Such pattern likely resides on the duplication history of Hox genes, which probably provided different possible evolutionary routes for achieving the same phenotypic solution.

Conclusions

Our results indicate that convergence occurred at a functional level of Hox genes along three independent origins of aquatic mammals. This conclusion reinforces the idea that different changes in developmental genes may lead to similar phenotypes, probably due to the redundancy provided by the participation of Hox paralogous genes in several developmental functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0682-4) contains supplementary material, which is available to authorized users.

The pathogenesis of congenital radial head dislocation/subluxation

The pathogenesis of congenital radial head dislocation/subluxation is unknown and has not been previously investigated. In this review, we explore the pathogenesis and define five different primary insults: collagen abnormalities, abnormal endochondral ossification of the developing growth plate, abnormalities of forearm ossification outside the growth plate, disproportionate growth of the radius and ulna, and altered HOX D expression/activity. Finally, the clinical relevance of our review is discussed.

A genome-wide association study identifies a genomic region for the polycerate phenotype in sheep (Ovis aries)

Horns are a cranial appendage found exclusively in Bovidae, and play important roles in accessing resources and mates. In sheep (Ovies aries), horns vary from polled to six-horned, and human have been selecting polled animals in farming and breeding. Here, we conducted a genome-wide association study on 24 two-horned versus 22 four-horned phenotypes in a native Chinese breed of Sishui Fur sheep. Together with linkage disequilibrium (LD) analyses and haplotype-based association tests, we identified a genomic region comprising 132.0–133.1 Mb on chromosome 2 that contained the top 10 SNPs (including 4 significant SNPs) and 5 most significant haplotypes associated with the polycerate phenotype. In humans and mice, this genomic region contains the HOXD gene cluster and adjacent functional genes EVX2 and KIAA1715, which have a close association with the formation of limbs and genital buds. Our results provide new insights into the genetic basis underlying variable numbers of horns and represent a new resource for use in sheep genetics and breeding.

Genome-wide association reveals the locus responsible for four-horned ruminant

Phenotypic variability in horn characteristics, such as their size, number and shape, offers the opportunity to elucidate the molecular basis of horn development. The objective of this study was to map the genetic determinant controlling the production of four horns in two breeds, Jacob sheep and Navajo-Churro, and examine whether an eyelid abnormality occurring in the same populations is related. Genome-wide association mapping was performed using 125 animals from the two breeds that contain two- and four-horned individuals. A case-control design analysis of 570 712 SNPs genotyped with the ovine HD SNP Beadchip revealed a strong association signal on sheep chromosome 2. The 10 most strongly associated SNPs were all located in a region spanning Mb positions 131.9-132.6, indicating the genetic architecture underpinning the production of four horns is likely to involve a single gene. The closest genes to the most strongly associated marker (OAR2_132568092) were MTX2 and the HOXD cluster, located approximately 93 Kb and 251 Kb upstream respectively. The occurrence of an eyelid malformation across both breeds was restricted to polled animals and those carrying more than two horns. This suggests the eyelid abnormality may be associated with departures from the normal developmental production of two-horned animals and that the two conditions are developmentally linked. This study demonstrated the presence of separate loci responsible for the polled and four-horned phenotypes in sheep and advanced our understanding of the complexity that underpins horn morphology in ruminants.

Key pathways regulated by HoxA9,10,11/HoxD9,10,11 during limb development

Background

The 39 mammalian Hox genes show problematic patterns of functional overlap. In order to more fully define the developmental roles of Hox genes it is necessary to remove multiple combinations of paralogous and flanking genes. In addition, the downstream molecular pathways regulated by Hox genes during limb development remain incompletely delineated.

Results

In this report we examine limb development in mice with frameshift mutations in six Hox genes, Hoxa9,10,11 and Hoxd9,10,11. The mice were made with a novel recombineering method that allows the simultaneous targeting of frameshift mutations into multiple flanking genes. The Hoxa9,10,11^−/−/Hoxd9,10,11^−/− mutant mice show a reduced ulna and radius that is more severe than seen in Hoxa11^−/−/Hoxd11^−/− mice, indicating a minor role for the flanking Hox9,10 genes in zeugopod development, as well as their primary function in stylopod development. The mutant mice also show severe reduction of Shh expression in the zone of polarizing activity, and decreased Fgf8 expression in the apical ectodermal ridge, thereby better defining the roles of these specific Hox genes in the regulation of critical signaling centers during limb development. Importantly, we also used laser capture microdissection coupled with RNA-Seq to characterize the gene expression programs in wild type and mutant limbs. Resting, proliferative and hypertrophic compartments of E15.5 forelimb zeugopods were examined. The results provide an RNA-Seq characterization of the progression of gene expression patterns during normal endochondral bone formation. In addition of the Hox mutants showed strongly altered expression of Pknox2, Zfp467, Gdf5, Bmpr1b, Dkk3, Igf1, Hand2, Shox2, Runx3, Bmp7 and Lef1, all of which have been previously shown to play important roles in bone formation.

Conclusions

The recombineering based frameshift mutation of the six flanking and paralogous Hoxa9,10,11 and Hoxd9,10,11 genes provides a resource for the analysis of their overlapping functions. Analysis of the Hoxa9,10,11^−/−/Hoxd9,10,11^−/− mutant limbs confirms and extends the results of previous studies using mice with Hox mutations in single paralogous groups or with entire Hox cluster deletions. The RNA-Seq analysis of specific compartments of the normal and mutant limbs defines the multiple key perturbed pathways downstream of these Hox genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-015-0078-5) contains supplementary material, which is available to authorized users.

A flexible bayesian model for testing for transmission ratio distortion

Current statistical approaches to investigate the nature and magnitude of transmission ratio distortion (TRD) are scarce and restricted to the most common experimental designs such as F2 populations and backcrosses. In this article, we describe a new Bayesian approach to check TRD within a given biallelic genetic marker in a diploid species, providing a highly flexible framework that can accommodate any kind of population structure. This model relies on the genotype of each offspring and thus integrates all available information from either the parents' genotypes or population-specific allele frequencies and yields TRD estimates that can be corroborated by the calculation of a Bayes factor (BF). This approach has been evaluated on simulated data sets with appealing statistical performance. As a proof of concept, we have also tested TRD in a porcine population with five half-sib families and 352 offspring. All boars and piglets were genotyped with the Porcine SNP60 BeadChip, whereas genotypes from the sows were not available. The SNP-by-SNP screening of the pig genome revealed 84 SNPs with decisive evidences of TRD (BF > 100) after accounting for multiple testing. Many of these regions contained genes related to biological processes (e.g., nucleosome assembly and co-organization, DNA conformation and packaging, and DNA complex assembly) that are critically associated with embryonic viability. The implementation of this method, which overcomes many of the limitations of previous approaches, should contribute to fostering research on TRD in both model and nonmodel organisms.

Unique expression patterns of multiple key genes associated with the evolution of mammalian flight

Bats are the only mammals capable of true flight. Critical adaptations for flight include a pair of dramatically elongated hands with broad wing membranes. To study the molecular mechanisms of bat wing evolution, we perform genomewide mRNA sequencing and in situ hybridization for embryonic bat limbs. We identify seven key genes that display unique expression patterns in embryonic bat wings and feet, compared with mouse fore- and hindlimbs. The expression of all 5'HoxD genes (Hoxd9-13) and Tbx3, six known crucial transcription factors for limb and digit development, is extremely high and prolonged in the elongating wing area. The expression of Fam5c, a tumour suppressor, in bat limbs is bat-specific and significantly high in all short digit regions (the thumb and foot digits). These results suggest multiple genetic changes occurred independently during the evolution of bat wings to elongate the hand digits, promote membrane growth and keep other digits short. Our findings also indicate that the evolution of limb morphology depends on the complex integration of multiple gene regulatory networks and biological processes that control digit formation and identity, chondrogenesis, and interdigital regression or retention.

Thumbs down: a molecular-morphogenetic approach to avian digit homology

Avian forelimb digit homology remains one of the standard themes in comparative biology and EvoDevo research. In order to resolve the apparent contradictions between embryological and paleontological evidence a variety of hypotheses have been presented in recent years. The proposals range from excluding birds from the dinosaur clade, to assignments of homology by different criteria, or even assuming a hexadactyl tetrapod limb ground state. At present two approaches prevail: the frame shift hypothesis and the pyramid reduction hypothesis. While the former postulates a homeotic shift of digit identities, the latter argues for a gradual bilateral reduction of phalanges and digits. Here we present a new model that integrates elements from both hypotheses with the existing experimental and fossil evidence. We start from the main feature common to both earlier concepts, the initiating ontogenetic event: reduction and loss of the anterior-most digit. It is proposed that a concerted mechanism of molecular regulation and developmental mechanics is capable of shifting the boundaries of hoxD expression in embryonic forelimb buds as well as changing the digit phenotypes. Based on a distinction between positional (topological) and compositional (phenotypic) homology criteria, we argue that the identity of the avian digits is II, III, IV, despite a partially altered phenotype. Finally, we introduce an alternative digit reduction scheme that reconciles the current fossil evidence with the presented molecular-morphogenetic model. Our approach identifies specific experiments that allow to test whether gene expression can be shifted and digit phenotypes can be altered by induced digit loss or digit gain.

Regulation of cell polarity in the cartilage growth plate and perichondrium of metacarpal elements by HOXD13 and WNT5A

The morphology of bones is genetically determined, but the molecular mechanisms that control shape, size and the overall gestalt of bones remain unclear. We previously showed that metacarpals in the synpolydactyly homolog (spdh) mouse, which carries a mutation in Hoxd13 similar to the human condition synpolydactyly (SPD), were transformed to carpal-like bones with cuboid shape that lack cortical bone and a perichondrium and are surrounded by a joint surface. Here we provide evidence that spdh metacarpal growth plates have a defect in cell polarization with a random instead of linear orientation. In parallel prospective perichondral cells failed to adopt the characteristic flattened cell shape. We observed a similar cell polarity defect in metacarpals of Wnt5a(-/-) mice. Wnt5a and the closely related Wnt5b were downregulated in spdh handplates, and HOXD13 induced expression of both genes in vitro. Concomitant we observed mislocalization of core planar cell polarity (PCP) components DVL2 and PRICKLE1 in spdh metacarpals indicating a defect in the WNT/PCP pathway. Conversely the WNT/β-CATENIN pathway, a hallmark of joint cells lining carpal bones, was upregulated in the perichondral region. Finally, providing spdh limb explant cultures with cells expressing either HOXD13 or WNT5A led to a non-cell autonomous partial rescue of cell polarity the perichondral region and restored the expression of perichondral markers. This study provides a so far unrecognized link between HOX proteins and cell polarity in the perichondrium and the growth plate, a failure of which leads to transformation of metacarpals to carpal-like structures.

Distinct global shifts in genomic binding profiles of limb malformation-associated HOXD13 mutations

Gene regulation by transcription factors (TFs) determines developmental programs and cell identity. Consequently, mutations in TFs can lead to dramatic phenotypes in humans by disrupting gene regulation. To date, the molecular mechanisms that actually cause these phenotypes have been difficult to address experimentally. ChIP-seq, which couples chromatin immunoprecipitation with high-throughput sequencing, allows TF function to be investigated on a genome-wide scale, enabling new approaches for the investigation of gene regulation. Here, we present the application of ChIP-seq to explore the effect of missense mutations in TFs on their genome-wide binding profile. Using a retroviral expression system in chicken mesenchymal stem cells, we elucidated the mechanism underlying a novel missense mutation in HOXD13 (Q317K) associated with a complex hand and foot malformation phenotype. The mutated glutamine (Q) is conserved in most homeodomains, a notable exception being bicoid-type homeodomains that have lysine (K) at this position. Our results show that the mutation results in a shift in the binding profile of the mutant toward a bicoid/PITX1 motif. Gene expression analysis and functional assays using in vivo overexpression studies confirm that the mutation results in a partial conversion of HOXD13 into a TF with bicoid/PITX1 properties. A similar shift was not observed with another mutation, Q317R, which is associated with brachysyndactyly, suggesting that the bicoid/PITX1-shift observed for Q317K might be related to the severe clinical phenotype. The methodology described can be used to investigate a wide spectrum of TFs and mutations that have not previously been amenable to ChIP-seq experiments.

DNA methylation and differentiation: HOX genes in muscle cells

BACKGROUND

Tight regulation of homeobox genes is essential for vertebrate development. In a study of genome-wide differential methylation, we recently found that homeobox genes, including those in the HOX gene clusters, were highly overrepresented among the genes with hypermethylation in the skeletal muscle lineage. Methylation was analyzed by reduced representation bisulfite sequencing (RRBS) of postnatal myoblasts, myotubes and adult skeletal muscle tissue and 30 types of non-muscle-cell cultures or tissues.

RESULTS

In this study, we found that myogenic hypermethylation was present in specific subregions of all four HOX gene clusters and was associated with various chromatin epigenetic features. Although the 3' half of the HOXD cluster was silenced and enriched in polycomb repression-associated H3 lysine 27 trimethylation in most examined cell types, including myoblasts and myotubes, myogenic samples were unusual in also displaying much DNA methylation in this region. In contrast, both HOXA and HOXC clusters displayed myogenic hypermethylation bordering a central region containing many genes preferentially expressed in myogenic progenitor cells and consisting largely of chromatin with modifications typical of promoters and enhancers in these cells. A particularly interesting example of myogenic hypermethylation was HOTAIR, a HOXC noncoding RNA gene, which can silence HOXD genes in trans via recruitment of polycomb proteins. In myogenic progenitor cells, the preferential expression of HOTAIR was associated with hypermethylation immediately downstream of the gene. Other HOX gene regions also displayed myogenic DNA hypermethylation despite being moderately expressed in myogenic cells. Analysis of representative myogenic hypermethylated sites for 5-hydroxymethylcytosine revealed little or none of this base, except for an intragenic site in HOXB5 which was specifically enriched in this base in skeletal muscle tissue, whereas myoblasts had predominantly 5-methylcytosine at the same CpG site.

CONCLUSIONS

Our results suggest that myogenic hypermethylation of HOX genes helps fine-tune HOX sense and antisense gene expression through effects on 5' promoters, intragenic and intergenic enhancers and internal promoters. Myogenic hypermethylation might also affect the relative abundance of different RNA isoforms, facilitate transcription termination, help stop the spread of activation-associated chromatin domains and stabilize repressive chromatin structures.

Toxicological effects of the different substances in tobacco smoke on human embryonic development by a systems chemo-biology approach

The physiological and molecular effects of tobacco smoke in adult humans and the development of cancer have been well described. In contrast, how tobacco smoke affects embryonic development remains poorly understood. Morphological studies of the fetuses of smoking pregnant women have shown various physical deformities induced by constant fetal exposure to tobacco components, especially nicotine. In addition, nicotine exposure decreases fetal body weight and bone/cartilage growth in addition to decreasing cranial diameter and tibia length. Unfortunately, the molecular pathways leading to these morphological anomalies are not completely understood. In this study, we applied interactome data mining tools and small compound interaction networks to elucidate possible molecular pathways associated with the effects of tobacco smoke components during embryonic development in pregnant female smokers. Our analysis showed a relationship between nicotine and 50 additional harmful substances involved in a variety of biological process that can cause abnormal proliferation, impaired cell differentiation, and increased oxidative stress. We also describe how nicotine can negatively affect retinoic acid signaling and cell differentiation through inhibition of retinoic acid receptors. In addition, nicotine causes a stress reaction and/or a pro-inflammatory response that inhibits the agonistic action of retinoic acid. Moreover, we show that the effect of cigarette smoke on the developing fetus could represent systemic and aggressive impacts in the short term, causing malformations during certain stages of development. Our work provides the first approach describing how different tobacco constituents affect a broad range of biological process in human embryonic development.