Molecular evolution of HoxA13 and the multiple origins of limbless morphologies in amphibians and reptiles

Developmental processes and their results, morphological characters, are inherited through transmission of genes regulating development. While there is ample evidence that cis-regulatory elements tend to be modular, with sequence segments dedicated to different roles, the situation for proteins is less clear, being particularly complex for transcription factors with multiple functions. Some motifs mediating protein-protein interactions may be exclusive to particular developmental roles, but it is also possible that motifs are mostly shared among different processes. Here we focus on HoxA13, a protein essential for limb development. We asked whether the HoxA13 amino acid sequence evolved similarly in three limbless clades: Gymnophiona, Amphisbaenia and Serpentes. We explored variation in ω (dN/dS) using a maximum-likelihood framework and HoxA13sequences from 47 species. Comparisons of evolutionary models provided low ω global values and no evidence that HoxA13 experienced relaxed selection in limbless clades. Branch-site models failed to detect evidence for positive selection acting on any site along branches of Amphisbaena and Gymnophiona, while three sites were identified in Serpentes. Examination of alignments did not reveal consistent sequence differences between limbed and limbless species. We conclude that HoxA13 has no modules exclusive to limb development, which may be explained by its involvement in multiple developmental processes

Enhanced genome assembly and a new official gene set for Tribolium castaneum

Background. The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality. Results. Here, we present an improved genome assembly (Tcas5.2) and an enhanced genome annotation resulting in a new official gene set (OGS3) for Tribolium castaneum, which significantly increase the quality of the genomic resources. By adding large-distance jumping library DNA sequencing to join scaffolds and fill small gaps, the gaps in the genome assembly were reduced and the N50 increased to 4753kbp. The precision of the gene models was enhanced by the use of a large body of RNA-Seq reads of different life history stages and tissue types, leading to the discovery of 1452 novel gene sequences. We also added new features such as alternative splicing, well defined UTRs and microRNA target predictions. For quality control, 399 gene models were evaluated by manual inspection. The current gene set was submitted to Genbank and accepted as a RefSeq genome by NCBI. Conclusions. The new genome assembly (Tcas5.2) and the official gene set (OGS3) provide enhanced genomic resources for genetic work in Tribolium castaneum. The much improved information on transcription start sites supports transgenic and gene editing approaches. Further, novel types of information such as splice variants and microRNA target genes open additional possibilities for analysis

Ubx Regulates Differential Enlargement and Diversification of Insect Hind Legs

Differential enlargement of hind (T3) legs represents one of the hallmarks of insect evolution. However, the actual mechanism(s) responsible are yet to be determined. To address this issue, we have now studied the molecular basis of T3 leg enlargement in Oncopeltus fasciatus (milkweed bug) and Acheta domesticus (house cricket). In Oncopeltus, the T3 tibia displays a moderate increase in size, whereas in Acheta, the T3 femur, tibia, and tarsus are all greatly enlarged. Here, we show that the hox gene Ultrabithorax (Ubx) is expressed in the enlarged segments of hind legs. Furthermore, we demonstrate that depletion of Ubx during embryogenesis has a primary effect in T3 legs and causes shortening of leg segments that are enlarged in a wild type. This result shows that Ubx is regulating the differential growth and enlargement of T3 legs in both Oncopeltus and Acheta. The emerging view suggests that Ubx was co-opted for a novel role in regulating leg growth and that the transcriptional modification of its expression may be a universal mechanism for the evolutionary diversification of insect hind legs

Measuring the functional sequence complexity of proteins

<p>Abstract</p> <p>Background</p> <p>Abel and Trevors have delineated three aspects of sequence complexity, Random Sequence Complexity (RSC), Ordered Sequence Complexity (OSC) and Functional Sequence Complexity (FSC) observed in biosequences such as proteins. In this paper, we provide a method to measure functional sequence complexity.</p> <p>Methods and Results</p> <p>We have extended Shannon uncertainty by incorporating the data variable with a functionality variable. The resulting measured unit, which we call Functional bit (Fit), is calculated from the sequence data jointly with the defined functionality variable. To demonstrate the relevance to functional bioinformatics, a method to measure functional sequence complexity was developed and applied to 35 protein families. Considerations were made in determining how the measure can be used to correlate functionality when relating to the whole molecule and sub-molecule. In the experiment, we show that when the proposed measure is applied to the aligned protein sequences of ubiquitin, 6 of the 7 highest value sites correlate with the binding domain.</p> <p>Conclusion</p> <p>For future extensions, measures of functional bioinformatics may provide a means to evaluate potential evolving pathways from effects such as mutations, as well as analyzing the internal structural and functional relationships within the 3-D structure of proteins.</p

Sex-Biased Expression of MicroRNAs in Schistosoma mansoni

Schistosomiasis is an important neglected tropical disease caused by digenean helminth parasites of the genus Schistosoma. Schistosomes are unusual in that they are dioecious and the adult worms live in the blood system. MicroRNAs play crucial roles during gene regulation and are likely to be important in sex differentiation in dioecious species. Here we characterize 112 microRNAs from adult Schistosoma mansoni individuals, including 84 novel microRNA families, and investigate the expression pattern in different sexes. By deep sequencing, we measured the relative expression levels of conserved and newly identified microRNAs between male and female samples. We observed that 13 microRNAs exhibited sex-biased expression, 10 of which are more abundant in females than in males. Sex chromosomes showed a paucity of female-biased genes, as predicted by theoretical evolutionary models. We propose that the recent emergence of separate sexes in Schistosoma had an effect on the chromosomal distribution and evolution of microRNAs, and that microRNAs are likely to participate in the sex differentiation/maintenance process

MiR-10 Represses HoxB1a and HoxB3a in Zebrafish

BACKGROUND: The Hox genes are involved in patterning the anterior-posterior axis. In addition to the protein coding Hox genes, the miR-10, miR-196 and miR-615 families of microRNA genes are conserved within the vertebrate Hox clusters. The members of the miR-10 family are located at positions associated with Hox-4 paralogues. No function is yet known for this microRNA family but the genomic positions of its members suggest a role in anterior-posterior patterning. METHODOLOGY/PRINCIPAL FINDINGS: Using sensor constructs, overexpression and morpholino knockdown, we show in Zebrafish that miR-10 targets HoxB1a and HoxB3a and synergizes with HoxB4 in the repression of these target genes. Overexpression of miR-10 also induces specific phenotypes related to the loss of function of these targets. HoxB1a and HoxB3a have a dominant hindbrain expression domain anterior to that of miR-10 but overlap in a weaker expression domain in the spinal cord. In this latter domain, miR-10 knockdown results in upregulation of the target genes. In the case of a HoxB3a splice variant that includes miR-10c within its primary transcript, we show that the microRNA acts in an autoregulatory fashion. CONCLUSIONS/SIGNIFICANCE: We find that miR-10 acts to repress HoxB1a and HoxB3a within the spinal cord and show that this repression works cooperatively with HoxB4. As with the previously described interactions between miR-196 and HoxA7 and Hox-8 paralogues, the target genes are located in close proximity to the microRNA. We present a model in which we postulate a link between the clustering of Hox genes and post-transcriptional gene regulation. We speculate that the high density of transcription units and enhancers within the Hox clusters places constraints on the precision of the transcriptional control that can be achieved within these clusters and requires the involvement of post-transcriptional gene silencing to define functional domains of genes appropriately

Insights into Hox Protein Function from a Large Scale Combinatorial Analysis of Protein Domains

Protein function is encoded within protein sequence and protein domains. However, how protein domains cooperate within a protein to modulate overall activity and how this impacts functional diversification at the molecular and organism levels remains largely unaddressed. Focusing on three domains of the central class Drosophila Hox transcription factor AbdominalA (AbdA), we used combinatorial domain mutations and most known AbdA developmental functions as biological readouts to investigate how protein domains collectively shape protein activity. The results uncover redundancy, interactivity, and multifunctionality of protein domains as salient features underlying overall AbdA protein activity, providing means to apprehend functional diversity and accounting for the robustness of Hox-controlled developmental programs. Importantly, the results highlight context-dependency in protein domain usage and interaction, allowing major modifications in domains to be tolerated without general functional loss. The non-pleoitropic effect of domain mutation suggests that protein modification may contribute more broadly to molecular changes underlying morphological diversification during evolution, so far thought to rely largely on modification in gene cis-regulatory sequences

Three Drosophila Hox Complex microRNAs Do Not Have Major Effects on Expression of Evolutionarily Conserved Hox Gene Targets during Embryogenesis

The discovery of microRNAs has resulted in a major expansion of the number of molecules known to be involved in gene regulation. Elucidating the functions of animal microRNAs has posed a significant challenge as their target interactions with messenger RNAs do not adhere to simple rules. Of the thousands of known animal microRNAs, relatively few microRNA:messenger RNA regulatory interactions have been biologically validated in an normal organismal context. Here we present evidence that three microRNAs from the Hox complex in Drosophila (miR-10-5p, miR-10-3p, miR-iab-4-5p) do not have significant effects during embryogenesis on the expression of Hox genes that contain high confidence microRNAs target sites in the 3′ untranslated regions of their messenger RNAs. This is significant, in that it suggests that many predicted microRNA-target interactions may not be biologically relevant, or that the outcomes of these interactions may be so subtle that mutants may only show phenotypes in specific contexts, such as in environmental stress conditions, or in combinations with other microRNA mutations

Isolation of Hox Cluster Genes from Insects Reveals an Accelerated Sequence Evolution Rate

Among gene families it is the Hox genes and among metazoan animals it is the insects (Hexapoda) that have attracted particular attention for studying the evolution of development. Surprisingly though, no Hox genes have been isolated from 26 out of 35 insect orders yet, and the existing sequences derive mainly from only two orders (61% from Hymenoptera and 22% from Diptera). We have designed insect specific primers and isolated 37 new partial homeobox sequences of Hox cluster genes (lab, pb, Hox3, ftz, Antp, Scr, abd-a, Abd-B, Dfd, and Ubx) from six insect orders, which are crucial to insect phylogenetics. These new gene sequences provide a first step towards comparative Hox gene studies in insects. Furthermore, comparative distance analyses of homeobox sequences reveal a correlation between gene divergence rate and species radiation success with insects showing the highest rate of homeobox sequence evolution

Are Hox Genes Ancestrally Involved in Axial Patterning? Evidence from the Hydrozoan Clytia hemisphaerica (Cnidaria)

Background: The early evolution and diversification of Hox-related genes in eumetazoans has been the subject of conflicting hypotheses concerning the evolutionary conservation of their role in axial patterning and the pre-bilaterian origin of the Hox and ParaHox clusters. The diversification of Hox/ParaHox genes clearly predates the origin of bilaterians. However, the existence of a "Hox code' predating the cnidarian-bilaterian ancestor and supporting the deep homology of axes is more controversial. This assumption was mainly based on the interpretation of Hox expression data from the sea anemone, but growing evidence from other cnidarian taxa puts into question this hypothesis. Methodology/Principal Findings: Hox, ParaHox and Hox-related genes have been investigated here by phylogenetic analysis and in situ hybridisation in Clytia hemisphaerica, an hydrozoan species with medusa and polyp stages alternating in the life cycle. Our phylogenetic analyses do not support an origin of ParaHox and Hox genes by duplication of an ancestral ProtoHox cluster, and reveal a diversification of the cnidarian HOX9-14 genes into three groups called A, B, C. Among the 7 examined genes, only those belonging to the HOX9-14 and the CDX groups exhibit a restricted expression along the oralaboral axis during development and in the planula larva, while the others are expressed in very specialised areas at the medusa stage. Conclusions/Significance: Cross species comparison reveals a strong variability of gene expression along the oral-aboral axis and during the life cycle among cnidarian lineages. The most parsimonious interpretation is that the Hox code, collinearity and conservative role along the antero-posterior axis are bilaterian innovations