Hox genes, evo-devo, and the case of the ftz gene

Expression of posterior Hox genes and opisthosomal appendage development in a mygalomorph spider

Spiders represent an evolutionary successful group of chelicerate arthropods. The body of spiders is subdivided into two regions (tagmata). The anterior tagma, the prosoma, bears the head appendages and four pairs of walking legs. The segments of the posterior tagma, the opisthosoma, either lost their appendages during the course of evolution or their appendages were substantially modified to fulfill new tasks such as reproduction, gas exchange, and silk production. Previous work has shown that the homeotic Hox genes are involved in shaping the posterior appendages of spiders. In this paper, we investigate the expression of the posterior Hox genes in a tarantula that possesses some key differences of posterior appendages compared to true spiders, such as the lack of the anterior pair of spinnerets and a second set of book lungs instead of trachea. Based on the observed differences in posterior Hox gene expression in true spiders and tarantulas, we argue that subtle changes in the Hox gene expression of the Hox genes abdA and AbdB are possibly responsible for at least some of the morphological differences seen in true spiders versus tarantulas.

Comparative genomics of Hox and ParaHox genes among major lineages of Branchiopoda with emphasis on tadpole shrimps

Hox and ParaHox genes (HPHGs) are key developmental genes that pattern regional identity along the anterior–posterior body axis of most animals. Here, we identified HPHGs in tadpole shrimps (Pancrustacea, Branchiopoda, Notostraca), an iconic example of the so-called “living fossils” and performed a comparative genomics analysis of HPHGs and the Hox cluster among major branchiopod lineages. Notostraca possess the entire Hox complement, and the Hox cluster seems to be split into two different subclusters, although we were not able to support this finding with chromosome-level assemblies. However, the genomic structure of Hox genes in Notostraca appears more derived than that of Daphnia spp., which instead retains the plesiomorphic condition of a single compact cluster. Spinicaudata and Artemia franciscana show instead a Hox cluster subdivided across two or more genomic scaffolds with some orthologs either duplicated or missing. Yet, branchiopod HPHGs are similar among the various clades in terms of both intron length and number, as well as in their pattern of molecular evolution. Sequence substitution rates are in fact generally similar for most of the branchiopod Hox genes and the few differences we found cannot be traced back to natural selection, as they are not associated with any signals of diversifying selection or substantial switches in selective modes. Altogether, these findings do not support a significant stasis in the Notostraca Hox cluster and further confirm how morphological evolution is not tightly associated with genome dynamics.

Hox genes in spiders: Their significance for development and evolution

Hox genes are known for their role in the specification of typical body plan features in animals. Evolutionary changes in Hox gene function are believed to be involved in the emergence of the diverse body plans we observe in animals today. Spiders share many body plan features with other arthropods, but also have numerous unique traits of their own. Studies of spider Hox genes have already provided insights into evolutionarily conserved and derived features of the spider body plan and their genetic basis. However, many aspects of Hox gene biology have been insufficiently studied in spiders so far. In this review, we highlight previous comparative studies of Hox genes in spiders and their significance for our understanding of the evolution of the spider body plan. We also identify aspects of Hox gene biology that need to be studied in greater detail. Many spider Hox genes have not been investigated beyond their mRNA expression patterns, and the role of Hox genes with apparently plesiomorphic or dual functions, like ftz and Hox3 is still unclear. Spiders have a duplicated Hox gene cluster, but possible sub- or neofunctionalisation of duplicates have not yet been studied systematically. Future research should therefore focus on these issues, in addition to the role of Polycomb and trithorax-mediated regulation, the identification of regulatory regions, cofactors or spider-specific target genes, and the significance of non-coding RNAs transcribed from within the Hox cluster and even from the antisense strand of particular Hox genes.

Patterning, From Conifers to Consciousness: Turing’s Theory and Order From Fluctuations

This is a brief account of Turing’s ideas on biological pattern and the events that led to their wider acceptance by biologists as a valid way to investigate developmental pattern, and of the value of theory more generally in biology. Periodic patterns have played a key role in this process, especially 2D arrays of oriented stripes, which proved a disappointment in theoretical terms in the case of Drosophila segmentation, but a boost to theory as applied to skin patterns in fish and model chemical reactions. The concept of “order from fluctuations” is a key component of Turing’s theory, wherein pattern arises by selective amplification of spatial components concealed in the random disorder of molecular and/or cellular processes. For biological examples, a crucial point from an analytical standpoint is knowing the nature of the fluctuations, where the amplifier resides, and the timescale over which selective amplification occurs. The answer clarifies the difference between “inelegant” examples such as Drosophila segmentation, which is perhaps better understood as a programmatic assembly process, and “elegant” ones expressible in equations like Turing’s: that the fluctuations and selection process occur predominantly in evolutionary time for the former, but in real time for the latter, and likewise for error suppression, which for Drosophila is historical, in being lodged firmly in past evolutionary events. The prospects for a further extension of Turing’s ideas to the complexities of brain development and consciousness is discussed, where a case can be made that it could well be in neuroscience that his ideas find their most important application.

Annotation of Hox cluster and Hox cofactor genes in the Asian citrus psyllid, Diaphorina citri, reveals novel features

Hox genes and their cofactors are essential developmental genes specifying regional identity in animals. Hox genes have a conserved arrangement in clusters in the same order in which they specify identity along the anterior-posterior axis. A few insect species have breaks in the cluster, but these are exceptions. We annotated the 10 Hox genes of the Asian citrus psyllid Diaphorina citri, and found a split in its Hox cluster between the Deformed and Sex combs reduced genes - the first time a break at this position has been observed in an insect Hox cluster. We also annotated D. citri orthologs of the Hox cofactor genes homothorax, PKNOX and extradenticle and found an additional copy of extradenticle in D. citri that appears to be a retrogene. Expression data and sequence conservation suggest that the extradenticle retrogene may have retained the original extradenticle function and allowed divergence of the parental extradenticle gene.

Panarthropod tiptop/teashirt and spalt orthologs and their potential role as "trunk"-selector genes

BACKGROUND

In the vinegar fly Drosophila melanogaster, the homeodomain containing transcription factor Teashirt (Tsh) appears to specify trunk identity in concert with the function of the Hox genes. While in Drosophila there is a second gene closely related to tsh, called tiptop (tio), in other arthropods species only one copy exists (called tio/tsh). The expression of tsh and tio/tsh, respectively, is surprisingly similar among arthropods suggesting that its function as trunk selector gene may be conserved. Other research, for example on the beetle Tribolium castaneum, questions even conservation of Tsh function among insects. The zinc-finger transcription factor Spalt (Sal) is involved in the regulation of Drosophila tsh, but this regulatory interaction does not appear to be conserved in Tribolium either. Whether the function and interaction of tsh and sal as potential trunk-specifiers, however, is conserved is still unclear because comparative studies on sal expression (except for Tribolium) are lacking, and functional data are (if at all existing) restricted to Insecta.

RESULTS

Here, we provide additional data on arthropod tsh expression, show the first data on onychophoran tio/tsh expression, and provide a comprehensive investigation on sal expression patterns in arthropods and an onychophoran.

CONCLUSIONS

Our data support the idea that tio/tsh genes are involved in the development of "trunk" segments by regulating limb development. Our data suggest further that the function of Sal is indeed unlikely to be conserved in trunk vs head development like in Drosophila, but early expression of sal is in line with a potential homeotic function, at least in Arthropoda.

speck, First Identified in Drosophila melanogaster in 1910, Is Encoded by the Arylalkalamine N-Acetyltransferase (AANAT1) Gene

The pigmentation mutation speck is a commonly used recombination marker characterized by a darkly pigmented region at the wing hinge. Identified in 1910 by Thomas Hunt Morgan, speck was characterized by Sturtevant as the most "workable" mutant in the rightmost region of the second chromosome and eventually localized to 2-107.0 and 60C1-2. Though the first speck mutation was isolated over 110 years ago, speck is still not associated with any gene. Here, as part of an undergraduate-led research effort, we show that speck is encoded by the Arylalkylamine N-acetyltransferase 1 (AANAT1) gene. Both alleles from the Morgan lab contain a retrotransposon in exon 1 of the RB transcript of the AANAT1 gene. We have also identified a new insertion allele and generated multiple deletion alleles in AANAT1 that all give a strong speck phenotype. In addition, expression of AANAT1 RNAi constructs either ubiquitously or in the dorsal portion of the developing wing generates a similar speck phenotype. We find that speck alleles have additional phenotypes, including ectopic pigmentation in the posterior pupal case, leg joints, cuticular sutures and overall body color. We propose that the acetylated dopamine generated by AANAT1 decreases the dopamine pool available for melanin production. When AANAT1 function is decreased, the excess dopamine enters the melanin pathway to generate the speck phenotype.

RNAi-Based Functional Genomics in Hemiptera

RNA interference (RNAi) is a powerful approach for sequence-specific gene silencing, displaying tremendous potential for functional genomics studies in hemipteran insects. Exploiting RNAi allows the biological roles of critical genes to be defined and aids the development of RNAi-based biopesticides. In this review, we provide context to the rapidly expanding field of RNAi-based functional genomics studies in hemipteran insects. We highlight the most widely used RNAi delivery strategies, including microinjection, oral ingestion and topical application. Additionally, we discuss the key variables affecting RNAi efficacy in hemipteran insects, including insect life-stage, gene selection, the presence of nucleases, and the role of core RNAi machinery. In conclusion, we summarise the application of RNAi in functional genomics studies in Hemiptera, focusing on genes involved in reproduction, behaviour, metabolism, immunity and chemical resistance across 33 species belonging to 14 families.

Bridging the gap between theory and practice in elucidating modular gene regulatory sequence organisation within genomes

Changes to promoter regions probably have been responsible for many morphological evolutionary transitions, especially in animals. This idea is becoming testable, as data from genome projects amass and enable bioinformaticians to conduct comparative sequence analyses and test for correlations between genotypic similarities or differences and phenotypic likeness or disparity. Although such practical pursuits have initiated some theoretical considerations, a conceptual framework for understanding promoter region evolution, potentially effecting morphological evolution, is only starting to emerge, predominantly resulting from computational research. We contribute to this framework by specifying three big problems for promoter region research; reviewing computational research on promoter region evolution; and exemplifying a topic for future promoter region research - module evolution.

Genomic insights into mite phylogeny, fitness, development, and reproduction

BACKGROUND

Predatory mites (Acari: Phytoseiidae) are the most important beneficial arthropods used in augmentative biological pest control of protected crops around the world. However, the genomes of mites are far less well understood than those of insects and the evolutionary relationships among mite and other chelicerate orders are contested, with the enigmatic origin of mites at one of the centres in discussion of the evolution of Arachnida.

RESULTS

We here report the 173 Mb nuclear genome (from 51.75 Gb pairs of Illumina reads) of the predatory mite, Neoseiulus cucumeris, a biocontrol agent against pests such as mites and thrips worldwide. We identified nearly 20.6 Mb (~ 11.93% of this genome) of repetitive sequences and annotated 18,735 protein-coding genes (a typical gene 2888 bp in size); the total length of protein-coding genes was about 50.55 Mb (29.2% of this assembly). About 37% (6981) of the genes are unique to N. cucumeris based on comparison with other arachnid genomes. Our phylogenomic analysis supported the monophyly of Acari, therefore rejecting the biphyletic origin of mites advocated by other studies based on limited gene fragments or few taxa in recent years. Our transcriptomic analyses of different life stages of N. cucumeris provide new insights into genes involved in its development. Putative genes involved in vitellogenesis, regulation of oviposition, sex determination, development of legs, signal perception, detoxification and stress-resistance, and innate immune systems are identified.

CONCLUSIONS

Our genomics and developmental transcriptomics analyses of N. cucumeris provide invaluable resources for further research on the development, reproduction, and fitness of this economically important mite in particular and Arachnida in general.

Arthropod segmentation

There is now compelling evidence that many arthropods pattern their segments using a clock-and-wavefront mechanism, analogous to that operating during vertebrate somitogenesis. In this Review, we discuss how the arthropod segmentation clock generates a repeating sequence of pair-rule gene expression, and how this is converted into a segment-polarity pattern by ‘timing factor’ wavefronts associated with axial extension. We argue that the gene regulatory network that patterns segments may be relatively conserved, although the timing of segmentation varies widely, and double-segment periodicity appears to have evolved at least twice. Finally, we describe how the repeated evolution of a simultaneous (Drosophila-like) mode of segmentation within holometabolan insects can be explained by heterochronic shifts in timing factor expression plus extensive pre-patterning of the pair-rule genes.

Integrated transcriptomics reveals master regulators of lung adenocarcinoma and novel repositioning of drug candidates

BACKGROUND

Lung adenocarcinoma is the major cause of cancer-related deaths in the world. Given this, the importance of research on its pathophysiology and therapy remains a key health issue. To assist in this endeavor, recent oncology studies are adopting Systems Biology approaches and bioinformatics to analyze and understand omics data, bringing new insights about this disease and its treatment.

METHODS

We used reverse engineering of transcriptomic data to reconstruct nontumorous lung reference networks, focusing on transcription factors (TFs) and their inferred target genes, referred as regulatory units or regulons. Afterwards, we used 13 case-control studies to identify TFs acting as master regulators of the disease and their regulatory units. Furthermore, the inferred activation patterns of regulons were used to evaluate patient survival and search drug candidates for repositioning.

RESULTS

The regulatory units under the influence of ATOH8, DACH1, EPAS1, ETV5, FOXA2, FOXM1, HOXA4, SMAD6, and UHRF1 transcription factors were consistently associated with the pathological phenotype, suggesting that they may be master regulators of lung adenocarcinoma. We also observed that the inferred activity of FOXA2, FOXM1, and UHRF1 was significantly associated with risk of death in patients. Finally, we obtained deptropine, promazine, valproic acid, azacyclonol, methotrexate, and ChemBridge ID compound 5109870 as potential candidates to revert the molecular profile leading to decreased survival.

CONCLUSION

Using an integrated transcriptomics approach, we identified master regulator candidates involved with the development and prognostic of lung adenocarcinoma, as well as potential drugs for repurposing.

Seeing is believing: the Bicoid protein reveals its path

In this commentary, I will review the latest findings on the Bicoid (Bcd) morphogen in Drosophila, a paradigm for gradient formation taught to biology students for more than two decades. “Seeing is believing” also summarizes the erroneous steps that were needed to elucidate the mechanisms of gradient formation and the path of movement of Bcd. Initially proclaimed as a dogma in 1988 and later incorporated into the SDD model where the broad diffusion of Bcd throughout the embryo was the predominant step leading to gradient formation, the SDD model was irrefutable for more than two decades until first doubts were raised in 2007 regarding the diffusion properties of Bcd associated with the SDD model. This led to re-thinking of the issue and the definition of a new model, termed the ARTS model which could explain most of the physical constraints that were inherently associated with the SDD model. In the ARTS model, gradient formation is mediated by the mRNA which is redistributed along cortical microtubules to form a mRNA gradient which is translated to form the protein gradient. Contrary to the SDD model, there is no Bcd diffusion from the tip. The ARTS model is also compatible with the observed cortical movement of Bcd. I will critically compare the SDD and the ARTS models as well as other models, analyze the major differences, and highlight the path where Bcd is localized during early nuclear cycles.

Characterization of NlHox3, an essential gene for embryonic development in Nilaparvata lugens

Hox genes encode transcriptional regulatory proteins that control axial patterning in all bilaterians. The brown planthopper (BPH), Nilaparvata lugens (Hemiptera: Delphacidae), is a destructive insect pest of rice plants in Asian countries. During analysis of the N. lugens transcriptome, we identified a Hox3-like gene (NlHox3) that was highly and specifically expressed in the embryonic stage. We performed functional analysis on the gene to identify its roles in embryonic development and its potential use as a target in RNA interference (RNAi) based pest control. The sequence analysis showed that NlHox3 was homologous to the Hox3 gene and was most closely related with zen of Drosophila. There were no significant differences in oviposition between the treated and control females after injecting double-stranded RNA of NlHox3 (dsNlHox3) into newly emerged female adult BPHs; however, there was a significant difference in the hatchability of those eggs laid, which no egg from the treated group hatched normally. Injecting female adult BPHs with dsNlHox3 led to necrosis of these offspring embryos, with eye reversal and undeveloped organs, suggesting that NlHox3 was an essential gene for embryonic development and might be a potential target for RNAi-based control of this insect pest.

Identification and characterization of homeobox (Hox) genes and conservation of the single Hox cluster (324.6 kb) in the water flea Daphnia magna

We report the complete sequence analysis of the entire complement of eight typical homeobox (Hox) genes (Lab, Pb, Dfd, Scr, Antp, Ubx, Abd-A, and Abd-B) and two other genes (Hox3 and Ftz) in a 324.6-kb region in the water flea Daphnia magna. In the cluster of D. magna Hox genes, we found one long interspersed nuclear element (LINE)/R2-NeSL between Ubx and Abd-A that was not present in Daphnia pulex Hox genes. In basal expression of Hox genes at different developmental stages, biothorax complex genes (Ubx, Abd-A, and Abd-B) and some antennapedia complex genes (Lab, Scr, Antp) were moderately expressed, but the Hox3 gene was barely expressed. Three homeobox genes (Antp, Ubx, Abd-A) were highly expressed at 6-7 days after release from the brood chamber and/or in the adult stage. The structural array and transcribed orientation of Dm-Hox genes were identical to those of the sister species D. pulex (∼340 kb), indicating that the Hox gene structure in daphnids is highly conserved. However, Dm- and Dp-Hox3, -deformed (Dfd), and -fushi tarazu (Ftz) genes varied from orthologous genes in pancrustacean species.

Conservation and variation in pair-rule gene expression and function in the intermediate-germ beetle Dermestes maculatus

A set of pair-rule (PR) segmentation genes (PRGs) promotes the formation of alternate body segments in Drosophila melanogaster Whereas Drosophila embryos are long-germ, with segments specified more or less simultaneously, most insects add segments sequentially as the germband elongates. The hide beetle Dermestes maculatus represents an intermediate between short- and long-germ development, ideal for comparative study of PRGs. We show that eight of nine Drosophila PRG orthologs are expressed in stripes in Dermestes Functional results parse these genes into three groups: Dmac-eve, -odd and -run play roles in both germband elongation and PR patterning; Dmac-slp and -prd function exclusively as complementary, classic PRGs, supporting functional decoupling of elongation and segment formation; and orthologs of ftz, ftz-f1, h and opa show more variable function in Dermestes and other species. While extensive cell death generally prefigured Dermestes PRG RNAi-mediated cuticle defects, an organized region with high mitotic activity near the margin of the segment addition zone is likely to have contributed to truncation of eve^RNAi embryos. Our results suggest general conservation of clock-like regulation of PR stripe addition in sequentially segmenting species while highlighting regulatory rewiring involving a subset of PRG orthologs.

A molecular view of onychophoran segmentation

This paper summarizes our current knowledge on the expression and assumed function of Drosophila and (other) arthropod segmentation gene orthologs in Onychophora, a closely related outgroup to Arthropoda. This includes orthologs of the so-called Drosophila segmentation gene cascade including the Hox genes, as well as other genetic factors and pathways involved in non-drosophilid arthropods. Open questions about and around the topic are addressed, such as the definition of segments in onychophorans, the unclear regulation of conserved expression patterns downstream of non-conserved factors, and the potential role of mesodermal patterning in onychophoran segmentation.

Activation of Ftz-F1-Responsive Genes through Ftz/Ftz-F1 Dependent Enhancers

The orphan nuclear receptor Ftz-F1 is expressed in all somatic nuclei in Drosophila embryos, but mutations result in a pair-rule phenotype. This was explained by the interaction of Ftz-F1 with the homeodomain protein Ftz that is expressed in stripes in the primordia of segments missing in either ftz-f1 or ftz mutants. Ftz-F1 and Ftz were shown to physically interact and coordinately activate the expression of ftz itself and engrailed by synergistic binding to composite Ftz-F1/Ftz binding sites. However, attempts to identify additional target genes on the basis of Ftz-F1/ Ftz binding alone has met with only limited success. To discern rules for Ftz-F1 target site selection in vivo and to identify additional target genes, a microarray analysis was performed comparing wildtype and ftz-f1 mutant embryos. Ftz-F1-responsive genes most highly regulated included engrailed and nine additional genes expressed in patterns dependent on both ftz and ftz-f1. Candidate enhancers for these genes were identified by combining BDTNP Ftz ChIP-chip data with a computational search for Ftz-F1 binding sites. Of eight enhancer reporter genes tested in transgenic embryos, six generated expression patterns similar to the corresponding endogenous gene and expression was lost in ftz mutants. These studies identified a new set of Ftz-F1 targets, all of which are co-regulated by Ftz. Comparative analysis of enhancers containing Ftz/Ftz-F1 binding sites that were or were not bona fide targets in vivo suggested that GAF negatively regulates enhancers that contain Ftz/Ftz-F1 binding sites but are not actually utilized. These targets include other regulatory factors as well as genes involved directly in morphogenesis, providing insight into how pair-rule genes establish the body pattern.

Polycomb and trithorax opposition in development and disease

Early discoveries in chromatin biology and epigenetics heralded new insights into organismal development. From these studies, two mediators of cellular differentiation were discovered: the Polycomb group (PcG) of transcriptional repressors, and the trithorax group (trxG) of transcriptional activators. These protein families, while opposed in function, work together to coordinate the appropriate cellular developmental programs that allow for both embryonic stem cell self-renewal and differentiation. Recently, both the PcG and trxG chromatin modulators have been observed to be deregulated in a wide spectrum diseases including developmental disorders and cancer. To understand the impact of these findings we outline the past, present, and future. WIREs Dev Biol 2016, 5:659-688. doi: 10.1002/wdev.244 For further resources related to this article, please visit the WIREs website.

Homeodomain proteins: an update

Here, we provide an update of our review on homeobox genes that we wrote together with Walter Gehring in 1994. Since then, comprehensive surveys of homeobox genes have become possible due to genome sequencing projects. Using the 103 Drosophila homeobox genes as example, we present an updated classification. In animals, there are 16 major classes, ANTP, PRD, PRD-LIKE, POU, HNF, CUT (with four subclasses: ONECUT, CUX, SATB, and CMP), LIM, ZF, CERS, PROS, SIX/SO, plus the TALE superclass with the classes IRO, MKX, TGIF, PBC, and MEIS. In plants, there are 11 major classes, i.e., HD-ZIP (with four subclasses: I to IV), WOX, NDX, PHD, PLINC, LD, DDT, SAWADEE, PINTOX, and the two TALE classes KNOX and BEL. Most of these classes encode additional domains apart from the homeodomain. Numerous insights have been obtained in the last two decades into how homeodomain proteins bind to DNA and increase their specificity by interacting with other proteins to regulate cell- and tissue-specific gene expression. Not only protein-DNA base pair contacts are important for proper target selection; recent experiments also reveal that the shape of the DNA plays a role in specificity. Using selected examples, we highlight different mechanisms of homeodomain protein-DNA interaction. The PRD class of homeobox genes was of special interest to Walter Gehring in the last two decades. The PRD class comprises six families in Bilateria, and tinkers with four different motifs, i.e., the PAIRED domain, the Groucho-interacting motif EH1 (aka Octapeptide or TN), the homeodomain, and the OAR motif. Homologs of the co-repressor protein Groucho are also present in plants (TOPLESS), where they have been shown to interact with small amphipathic motives (EAR), and in yeast (TUP1), where we find an EH1-like motif in MATα2.