The Homeodomain Resource: a comprehensive collection of sequence, structure, interaction, genomic and functional information on the homeodomain protein family

Effect of polycyclic aromatic hydrocarbons on homeobox gene expression during embryonic development of cuttlefish, Sepia officinalis

Polycyclic aromatic hydrocarbons (PAHs) are persistent organic pollutants (POPs) commonly found in marine environments. Their bioaccumulation can cause harm to aquatic organisms, including invertebrates, particularly during the early stages of embryonic development. In this study, we evaluated, for the first time, the patterns of PAH accumulation in both capsule and embryo of common cuttlefish (Sepia officinalis). In addition, we explored the effects of PAHs by analysing the expression profiles of seven homeobox genes [i.e., gastrulation brain homeobox (GBX), paralogy group labial/Hox1 (HOX1), paralogy group Hox3 (HOX3), dorsal root ganglia homeobox (DRGX), visual system homeobox (VSX), aristaless-like homeobox (ARX) and LIM-homeodomain transcription factor (LHX3/4)]. We found that PAH levels in egg capsules were higher than those observed in chorion membranes (35.1 ± 13.3 ng/g vs 16.4 ± 5.9 ng/g). Furthermore, PAHs were also found in perivitellin fluid (11.5 ± 5.0 ng/ml). Naphthalene and acenaphthene were the congeners present at highest concentrations in each analysed egg component suggesting higher bioaccumulation rates. Embryos with high concentrations of PAHs also showed a significant increase in mRNA expression for each of the analysed homeobox genes. In particular, we observed a 15-fold increase in the ARX expression levels. Additionally, the statistically significant variation in homeobox gene expression patterns was accompanied by a concomitant increase in mRNA levels of both aryl hydrocarbon receptor (AhR) and estrogen receptor (ER). These findings suggest that bioaccumulation of PAHs may modulate developmental processes of cuttlefish embryos by targeting homeobox gene-mediated transcriptional outcomes. Mechanisms underlying the upregulation of homeobox genes could be related to the ability of PAHs to directly activate AhR- or ER-related signaling pathways.

Virtual screening of potential anticancer drugs based on microbial products

The development of microbial products for cancer treatment has been in the spotlight in recent years. In order to accelerate the lengthy and expensive drug development process, in silico screening tools are systematically employed, especially during the initial discovery phase. Moreover, considering the steadily increasing number of molecules approved by authorities for commercial use, there is a demand for faster methods to repurpose such drugs. Here we present a review on virtual screening web tools, such as publicly available databases of molecular targets and libraries of ligands, with the aim to facilitate the discovery of potential anticancer drugs based on microbial products. We provide an entry-level step-by-step description of the workflow for virtual screening of microbial metabolites with known protein targets, as well as two practical examples using freely available web tools. The first case presents a virtual screening study of drugs developed from microbial products using Caver Web, a web tool that performs docking along a tunnel. The second case comprises a comparative analysis between a wild type isocitrate dehydrogenase 1 and a mutant that results in cancer, using the recently developed web tool PredictSNP^Onco. In summary, this review provides the basic and essential background information necessary for virtual screening experiments, which may accelerate the discovery of novel anticancer drugs.

Twist-1 is upregulated by NSD2 and contributes to tumour dissemination and an epithelial-mesenchymal transition-like gene expression signature in t(4;14)-positive multiple myeloma

Approximately 15% of patients with multiple myeloma (MM) harbour the t(4;14) chromosomal translocation, leading to the overexpression of the histone methyltransferase NSD2. Patients with this translocation display increased tumour dissemination, accelerated disease progression and rapid relapse. Using publicly available gene expression profile data from NSD2^high (n = 135) and NSD2^low (n = 878) MM patients, we identified 39 epithelial-mesenchymal transition (EMT)-associated genes which are overexpressed in NSD2^high MM plasma cells. In addition, our analyses identified Twist-1 as a key transcription factor upregulated in NSD2^high MM patients and t(4;14)-positive cell lines. Overexpression and knockdown studies confirmed that Twist-1 is involved in driving the expression of EMT-associated genes in the human MM cell line KMS11 and promoted the migration of myeloma cell lines in vitro. Notably, Twist-1 overexpression in the mouse MM cell line 5TGM1 significantly increased tumour dissemination in an intratibial tumour model. These findings demonstrate that Twist-1, downstream of NSD2, contributes to the induction of an EMT-like signature in t(4;14)-positive MM and enhances the dissemination of MM plasma cells in vivo, which may, in part, explain the aggressive disease features associated with t(4;14)-positive MM.

Characterisation of transcription factor profiles in polycystic kidney disease (PKD): identification and validation of STAT3 and RUNX1 in the injury/repair response and PKD progression

Abstract

Autosomal dominant polycystic kidney disease (ADPKD) is the most common genetic renal disease, caused in the majority of the cases by a mutation in either the PKD1 or the PKD2 gene. ADPKD is characterised by a progressive increase in the number and size of cysts, together with fibrosis and distortion of the renal architecture, over the years. This is accompanied by alterations in a complex network of signalling pathways. However, the underlying molecular mechanisms are not well characterised. Previously, we defined the PKD Signature, a set of genes typically dysregulated in PKD across different disease models from a meta-analysis of expression profiles. Given the importance of transcription factors (TFs) in modulating disease, we focused in this paper on characterising TFs from the PKD Signature. Our results revealed that out of the 1515 genes in the PKD Signature, 92 were TFs with altered expression in PKD, and 32 of those were also implicated in tissue injury/repair mechanisms. Validating the dysregulation of these TFs by qPCR in independent PKD and injury models largely confirmed these findings. STAT3 and RUNX1 displayed the strongest activation in cystic kidneys, as demonstrated by chromatin immunoprecipitation (ChIP) followed by qPCR. Using immunohistochemistry, we showed a dramatic increase of expression after renal injury in mice and cystic renal tissue of mice and humans. Our results suggest a role for STAT3 and RUNX1 and their downstream targets in the aetiology of ADPKD and indicate that the meta-analysis approach is a viable strategy for new target discovery in PKD.

Key messages

We identified a list of transcription factors (TFs) commonly dysregulated in ADPKD.

Out of the 92 TFs identified in the PKD Signature, 35% are also involved in injury/repair processes.

STAT3 and RUNX1 are the most significantly dysregulated TFs after injury and during PKD progression.

STAT3 and RUNX1 activity is increased in cystic compared to non-cystic mouse kidneys.

Increased expression of STAT3 and RUNX1 is observed in the nuclei of renal epithelial cells, also in human ADPKD samples.

Electronic supplementary material

The online version of this article (10.1007/s00109-019-01852-3) contains supplementary material, which is available to authorized users.

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

New Insights into Cooperative Binding of Homeodomain Transcription Factors PREP1 and PBX1 to DNA

PREP1 and PBX1 are homeodomain (HD) transcription factors that play crucial roles in embryonic development. Here, we present the first biophysical characterization of a PREP1 HD, and the NMR spectroscopic study of its DNA binding pocket. The data show that residues flanking the HD participate in DNA binding. The kinetic parameters for DNA binding of individual PREP1 and PBX1 HDs, and of their combination, show that isolated PREP1 and PBX1 HDs bind to DNA in a cooperative manner. A novel PREP1 motif, flanking the HD at the C-terminus, is required for cooperativity.

Extending the functions of the homeotic transcription factor Cdx2 in the digestive system through nontranscriptional activities

The homeoprotein encoded by the intestinal-specific Cdx2 gene is a major regulator of gut development and homeostasis, also involved in colon cancer as well as in intestinal-type metaplasias when it is abnormally expressed outside the gut. At the molecular level, structure/function studies have demonstrated that the Cdx2 protein is a transcription factor containing a conserved homeotic DNA-binding domain made of three alpha helixes arranged in a helix-turn-helix motif, preceded by a transcriptional domain and followed by a regulatory domain. The protein interacts with several thousand sites on the chromatin and widely regulates intestinal functions in stem/progenitor cells as well as in mature differentiated cells. Yet, this transcription factor also acts trough original nontranscriptional mechanisms. Indeed, the identification of novel protein partners of Cdx2 and also of a splicing variant revealed unexpected functions in the control of signaling pathways like the Wnt and NF-κB pathways, in double-strand break DNA repair and in premessenger RNA splicing. These novel functions of Cdx2 must be considered to fully understand the complexity of the role of Cdx2 in the healthy intestine and in diseases.

Establishment and analysis of a reference transcriptome for Spodoptera frugiperda

Background

Spodoptera frugiperda (Noctuidae) is a major agricultural pest throughout the American continent. The highly polyphagous larvae are frequently devastating crops of importance such as corn, sorghum, cotton and grass. In addition, the Sf9 cell line, widely used in biochemistry for in vitro protein production, is derived from S. frugiperda tissues. Many research groups are using S. frugiperda as a model organism to investigate questions such as plant adaptation, pest behavior or resistance to pesticides.

Results

In this study, we constructed a reference transcriptome assembly (Sf_TR2012b) of RNA sequences obtained from more than 35 S. frugiperda developmental time-points and tissue samples. We assessed the quality of this reference transcriptome by annotating a ubiquitous gene family - ribosomal proteins - as well as gene families that have a more constrained spatio-temporal expression and are involved in development, immunity and olfaction. We also provide a time-course of expression that we used to characterize the transcriptional regulation of the gene families studied.

Conclusion

We conclude that the Sf_TR2012b transcriptome is a valid reference transcriptome. While its reliability decreases for the detection and annotation of genes under strong transcriptional constraint we still recover a fair percentage of tissue-specific transcripts. That allowed us to explore the spatial and temporal expression of genes and to observe that some olfactory receptors are expressed in antennae and palps but also in other non related tissues such as fat bodies. Similarly, we observed an interesting interplay of gene families involved in immunity between fat bodies and antennae.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-704) contains supplementary material, which is available to authorized users.

Annotation enrichment analysis: an alternative method for evaluating the functional properties of gene sets

Gene annotation databases (compendiums maintained by the scientific community that describe the biological functions performed by individual genes) are commonly used to evaluate the functional properties of experimentally derived gene sets. Overlap statistics, such as Fishers Exact test (FET), are often employed to assess these associations, but don't account for non-uniformity in the number of genes annotated to individual functions or the number of functions associated with individual genes. We find FET is strongly biased toward over-estimating overlap significance if a gene set has an unusually high number of annotations. To correct for these biases, we develop Annotation Enrichment Analysis (AEA), which properly accounts for the non-uniformity of annotations. We show that AEA is able to identify biologically meaningful functional enrichments that are obscured by numerous false-positive enrichment scores in FET, and we therefore suggest it be used to more accurately assess the biological properties of gene sets.

Disruption of aryl hydrocarbon receptor homeostatic levels during embryonic stem cell differentiation alters expression of homeobox transcription factors that control cardiomyogenesis

BACKGROUND

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor that regulates the expression of xenobiotic detoxification genes and is a critical mediator of gene-environment interactions. Many AHR target genes identified by genome-wide gene expression profiling have morphogenetic functions, suggesting that AHR may play a role in embryonic development.

OBJECTIVES

To characterize the developmental functions of the AHR, we studied the consequences of AHR activation by the agonist 2,3,7,8-tetrachlorodibenzo-p-doxin (TCDD), and the result of its repression by the antagonists 6,2,4-trimethoxyflavone and CH 223191 or by short-hairpin RNA (shRNA)-mediated Ahr knockdown during spontaneous differentiation of embryonic stem (ES) cells into cardiomyocytes.

METHODS

We generated an AHR-positive cardiomyocyte lineage differentiated from mouse ES cells that expresses puromycin resistance and enhanced green fluorescent protein (eGFP) under the control of the Cyp1a1 (cytochrome P450 1a1) promoter. We used RNA sequencing (RNA.Seq) to analyze temporal trajectories of TCDD-dependent global gene expression in these cells during differentiation.

RESULTS

Activation, inhibition, and knockdown of Ahr significantly inhibited the formation of contractile cardiomyocyte nodes. Global expression analysis of AHR-positive cells showed that activation of the AHR/TCDD axis disrupted the concerted expression of genes that regulate multiple signaling pathways involved in cardiac and neural morphogenesis and differentiation, including dozens of genes encoding homeobox transcription factors and Polycomb and trithorax group proteins.

CONCLUSIONS

Disruption of AHR expression levels resulted in gene expression changes that perturbed cardiomyocyte differentiation. The main function of the AHR during development appears to be the coordination of a complex regulatory network responsible for attainment and maintenance of cardiovascular homeostasis.

Searching for likeness in a database of macromolecular complexes

A software tool and workflow based on distance geometry is presented that can be used to search for local similarity in substructures in a comprehensive database of experimentally derived macromolecular structure. The method does not rely on fold annotation, specific secondary structure assignments, or sequence homology and may be used to locate compound substructures of multiple segments spanning different macromolecules that share a queried backbone geometry. This generalized substructure searching capability is intended to allow users to play an active part in exploring the role specific substructures play in larger protein domains, quaternary assemblies of proteins, and macromolecular complexes of proteins and polynucleotides. The user may select any portion or portions of an existing structure or complex to serve as a template for searching, and other structures that share the same structural features are identified, retrieved and overlaid to emphasize substructural likeness. Matching structures may be compared using a variety of integrated tools including molecular graphics for structure visualization and matching substructure sequence logos. A number of examples are provided that illustrate how generalized substructure searching may be used to understand both the similarity, and individuality of specific macromolecular structures. Web-based access to our substructure searching services is freely available at https://drugsite.msi.umn.edu.

p53-Dependent transcriptional responses to interleukin-3 signaling

p53 is critical in the normal response to a variety of cellular stresses including DNA damage and loss of p53 function is a common feature of many cancers. In hematological malignancies, p53 deletion is less common than in solid malignancies but is associated with poor prognosis and resistance to chemotherapy. Compared to their wild-type (WT) counterparts, hematopoietic progenitor cells lacking p53 have a greater propensity to survive cytokine loss, in part, due to the failure to transcribe Puma, a proapoptotic Bcl-2 family member. Using expression arrays, we have further characterized the differences that distinguish p53^−/− cells from WT myeloid cells in the presence of Interleukin-3 (IL-3) to determine if such differences contribute to the increased clonogenicity and survival responses observed in p53^−/− cells. We show that p53^−/− cells have a deregulated intracellular signaling environment and display a more rapid and sustained response to IL-3. This was accompanied by an increase in active ERK1/2 and a dependence on an intact MAP kinase signaling pathway. Contrastingly, we find that p53^−/− cells are independent on AKT for their survival. Thus, loss of p53 in myeloid cells results in an altered transcriptional and kinase signaling environment that favors enhanced cytokine signaling.

DNA search efficiency is modulated by charge composition and distribution in the intrinsically disordered tail

Intrinsically disordered tails are common in DNA-binding proteins and can affect their search efficiency on nonspecific DNA by promoting the brachiation dynamics of intersegment transfer. During brachiation, the protein jumps between distant DNA regions via an intermediate state in which the tail and globular moieties are bound to different DNA segments. While the disordered tail must be long and positively charged to facilitate DNA search, the effect of its residue sequence on brachiation is unknown. We explored this issue using the NK-2 and Antp homeodomain transcription factors. We designed 566 NK-2 tail-variants and 55 Antp tail-variants having different net charges and positive charge distributions and studied their dynamics and DNA search efficiencies using coarse-grained molecular dynamics simulations. More intersegment transfers occur when the tail is moderately positively charged and the positive charges are clustered together in the middle of the tail or towards its N terminus. The presence of a negatively charged residue does not significantly affect protein brachiation, although it is likely that the presence of many negatively charged residues will complicate the DNA search mechanism. A bioinformatic analysis of 1,384 wild-type homeodomains illustrates that the charge composition and distribution in their N-tail sequences are consistent with an optimal charge pattern to promote intersegment transfer. Our study thus indicates that the residue sequence of the disordered tails of DNA-binding proteins has unique characteristics that were evolutionarily selected to achieve optimized function and suggests that the sequence-structure-function paradigm known for structured proteins is valid for intrinsically disordered proteins as well.

The homeodomain complement of the ctenophore Mnemiopsis leidyi suggests that Ctenophora and Porifera diverged prior to the ParaHoxozoa

Background

The much-debated phylogenetic relationships of the five early branching metazoan lineages (Bilateria, Cnidaria, Ctenophora, Placozoa and Porifera) are of fundamental importance in piecing together events that occurred early in animal evolution. Comparisons of gene content between organismal lineages have been identified as a potentially useful methodology for phylogenetic reconstruction. However, these comparisons require complete genomes that, until now, did not exist for the ctenophore lineage. The homeobox superfamily of genes is particularly suited for these kinds of gene content comparisons, since it is large, diverse, and features a highly conserved domain.

Results

We have used a next-generation sequencing approach to generate a high-quality rough draft of the genome of the ctenophore Mnemiopsis leidyi and subsequently identified a set of 76 homeobox-containing genes from this draft. We phylogenetically categorized this set into established gene families and classes and then compared this set to the homeodomain repertoire of species from the other four early branching metazoan lineages. We have identified several important classes and subclasses of homeodomains that appear to be absent from Mnemiopsis and from the poriferan Amphimedon queenslandica. We have also determined that, based on lineage-specific paralog retention and average branch lengths, it is unlikely that these missing classes and subclasses are due to extensive gene loss or unusually high rates of evolution in Mnemiopsis.

Conclusions

This paper provides a first glimpse of the first sequenced ctenophore genome. We have characterized the full complement of Mnemiopsis homeodomains from this species and have compared them to species from other early branching lineages. Our results suggest that Porifera and Ctenophora were the first two extant lineages to diverge from the rest of animals. Based on this analysis, we also propose a new name - ParaHoxozoa - for the remaining group that includes Placozoa, Cnidaria and Bilateria.