Low affinity binding site clusters confer hox specificity and regulatory robustness.

In animals, Hox transcription factors define regional identity in distinct anatomical domains. How Hox genes encode this specificity is a paradox, because different Hox proteins bind with high affinity in vitro to similar DNA sequences. Here, we demonstrate that the Hox protein Ultrabithorax (Ubx) in complex with its cofactor Extradenticle (Exd) bound specifically to clusters of very low affinity sites in enhancers of the shavenbaby gene of Drosophila. These low affinity sites conferred specificity for Ubx binding in vivo, but multiple clustered sites were required for robust expression when embryos developed in variable environments. Although most individual Ubx binding sites are not evolutionarily conserved, the overall enhancer architecture-clusters of low affinity binding sites-is maintained and required for enhancer function. Natural selection therefore works at the level of the enhancer, requiring a particular density of low affinity Ubx sites to confer both specific and robust expression.

Looking across the gap: Understanding the evolution of eyes and vision among insects.

The compound eyes of insects exhibit stunning variation in size, structure, and function, which has allowed these animals to use their vision to adapt to a huge range of different environments and lifestyles, and evolve complex behaviors. Much of our knowledge of eye development has been learned from Drosophila, while visual adaptations and behaviors are often more striking and better understood from studies of other insects. However, recent studies in Drosophila and other insects, including bees, beetles, and butterflies, have begun to address this gap by revealing the genetic and developmental bases of differences in eye morphology and key new aspects of compound eye structure and function. Furthermore, technical advances have facilitated the generation of high-resolution connectomic data from different insect species that enhances our understanding of visual information processing, and the impact of changes in these processes on the evolution of vision and behavior. Here, we review these recent breakthroughs and propose that future integrated research from the development to function of visual systems within and among insect species represents a great opportunity to understand the remarkable diversification of insect eyes and vision.

Evolution of compound eye morphology underlies differences in vision between closely related Drosophila species.

BACKGROUND: Insects have evolved complex visual systems and display an astonishing range of adaptations for diverse ecological niches. Species of Drosophila melanogaster subgroup exhibit extensive intra- and interspecific differences in compound eye size. These differences provide an excellent opportunity to better understand variation in insect eye structure and the impact on vision. Here we further explored the difference in eye size between D. mauritiana and its sibling species D. simulans. RESULTS: We confirmed that D. mauritiana have rapidly evolved larger eyes as a result of more and wider ommatidia than D. simulans since they recently diverged approximately 240,000 years ago. The functional impact of eye size, and specifically ommatidia size, is often only estimated based on the rigid surface morphology of the compound eye. Therefore, we used 3D synchrotron radiation tomography to measure optical parameters in 3D, predict optical capacity, and compare the modelled vision to in vivo optomotor responses. Our optical models predicted higher contrast sensitivity for D. mauritiana, which we verified by presenting sinusoidal gratings to tethered flies in a flight arena. Similarly, we confirmed the higher spatial acuity predicted for Drosophila simulans with smaller ommatidia and found evidence for higher temporal resolution. CONCLUSIONS: Our study demonstrates that even subtle differences in ommatidia size between closely related Drosophila species can impact the vision of these insects. Therefore, further comparative studies of intra- and interspecific variation in eye morphology and the consequences for vision among other Drosophila species, other dipterans and other insects are needed to better understand compound eye structure-function and how the diversification of eye size, shape, and function has helped insects to adapt to the vast range of ecological niches.

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication.

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

The Evolution of Sox Gene Repertoires and Regulation of Segmentation in Arachnids.

The Sox family of transcription factors regulates many processes during metazoan development, including stem cell maintenance and nervous system specification. Characterizing the repertoires and roles of these genes can therefore provide important insights into animal evolution and development. We further characterized the Sox repertoires of several arachnid species with and without an ancestral whole-genome duplication and compared their expression between the spider Parasteatoda tepidariorum and the harvestman Phalangium opilio. We found that most Sox families have been retained as ohnologs after whole-genome duplication and evidence for potential subfunctionalization and/or neofunctionalization events. Our results also suggest that Sox21b-1 likely regulated segmentation ancestrally in arachnids, playing a similar role to the closely related SoxB gene, Dichaete, in insects. We previously showed that Sox21b-1 is required for the simultaneous formation of prosomal segments and sequential addition of opisthosomal segments in P. tepidariorum. We studied the expression and function of Sox21b-1 further in this spider and found that although this gene regulates the generation of both prosomal and opisthosomal segments, it plays different roles in the formation of these tagmata reflecting their contrasting modes of segmentation and deployment of gene regulatory networks with different architectures.

Unraveling the Genetic Basis for the Rapid Diversification of Male Genitalia between Drosophila Species.

In the last 240,000 years, males of the Drosophila simulans species clade have evolved striking differences in the morphology of their epandrial posterior lobes and claspers (surstyli). These appendages are used for grasping the female during mating and so their divergence is most likely driven by sexual selection. Mapping studies indicate a highly polygenic and generally additive genetic basis for these morphological differences. However, we have limited understanding of the gene regulatory networks that control the development of genital structures and how they evolved to result in this rapid phenotypic diversification. Here, we used new D. simulans/D. mauritiana introgression lines on chromosome arm 3L to generate higher resolution maps of posterior lobe and clasper differences between these species. We then carried out RNA-seq on the developing genitalia of both species to identify the expressed genes and those that are differentially expressed between the two species. This allowed us to test the function of expressed positional candidates during genital development in D. melanogaster. We identified several new genes involved in the development and possibly the evolution of these genital structures, including the transcription factors Hairy and Grunge. Furthermore, we discovered that during clasper development Hairy negatively regulates tartan (trn), a gene known to contribute to divergence in clasper morphology. Taken together, our results provide new insights into the regulation of genital development and how this has evolved between species.

tartan underlies the evolution of Drosophila male genital morphology.

Male genital structures are among the most rapidly evolving morphological traits and are often the only features that can distinguish closely related species. This process is thought to be driven by sexual selection and may reinforce species separation. However, while the genetic bases of many phenotypic differences have been identified, we still lack knowledge about the genes underlying evolutionary differences in male genital organs and organ size more generally. The claspers (surstyli) are periphallic structures that play an important role in copulation in insects. Here, we show that divergence in clasper size and bristle number between Drosophila mauritiana and Drosophila simulans is caused by evolutionary changes in tartan (trn), which encodes a transmembrane leucine-rich repeat domain protein that mediates cell-cell interactions and affinity. There are no fixed amino acid differences in trn between D. mauritiana and D. simulans, but differences in the expression of this gene in developing genitalia suggest that cis-regulatory changes in trn underlie the evolution of clasper morphology in these species. Finally, analyses of reciprocal hemizygotes that are genetically identical, except for the species from which the functional allele of trn originates, determined that the trn allele of D. mauritiana specifies larger claspers with more bristles than the allele of D. simulans Therefore, we have identified a gene underlying evolutionary change in the size of a male genital organ, which will help to better understand not only the rapid diversification of these structures, but also the regulation and evolution of organ size more broadly.

Wnt gene regulation and function during maxillary palp development in Drosophila melanogaster.

Wnt genes encode secreted ligands that play many important roles in the development of metazoans. There are thirteen known Wnt gene subfamilies and seven of these are represented in Drosophila melanogaster. While wingless (wg) is the best understood and most widely studied Wnt gene in Drosophila, the functions of many of the other Drosophila Wnt genes are less well understood. For example, relatively little is known about Wnt6, which is an ancient paralog of wg and they form a conserved Wnt cluster together with Wnt9 (Dwnt4) and Wnt10. Wg and Wnt6 encode similar proteins and exhibit overlapping expression in several tissues during development. Both wg and Wnt6 were previously shown to regulate the development of maxillary palps, important olfactory organs in flies, but it remained unclear how these two ligands may combine to carry out specific functions and how this is regulated. Here, we have further analysed Wnt6 function in the context of maxillary palp development. Surprisingly, we found that Wnt6 does not appear to be necessary for development of maxillary palps. While a deletion of the 5' region of Wnt6 results in very small maxillary palps, we show that this effect is more likely to be a consequence of removing cis-regulatory elements that may regulate wg expression in this tissue rather than through the loss of Wnt6 function. Although, we cannot completely exclude the possibility that Wnt6 may subtly regulate maxillary palp development in combination with wg, our analysis of Wnt6 loss of function mutants suggests this ligand plays a more general role in regulating growth during development. Taken together our results provide new insights into maxillary palp formation and Wnt6 functions in Drosophila, and further evidence for a complex cis-regulatory landscape in the Wnt9-wg-Wnt6-Wnt10 cluster, which may help explain its evolutionary conservation.

Gene regulatory network architecture in different developmental contexts influences the genetic basis of morphological evolution.

Convergent phenotypic evolution is often caused by recurrent changes at particular nodes in the underlying gene regulatory networks (GRNs). The genes at such evolutionary 'hotspots' are thought to maximally affect the phenotype with minimal pleiotropic consequences. This has led to the suggestion that if a GRN is understood in sufficient detail, the path of evolution may be predictable. The repeated evolutionary loss of larval trichomes among Drosophila species is caused by the loss of shavenbaby (svb) expression. svb is also required for development of leg trichomes, but the evolutionary gain of trichomes in the 'naked valley' on T2 femurs in Drosophila melanogaster is caused by reduced microRNA-92a (miR-92a) expression rather than changes in svb. We compared the expression and function of components between the larval and leg trichome GRNs to investigate why the genetic basis of trichome pattern evolution differs in these developmental contexts. We found key differences between the two networks in both the genes employed, and in the regulation and function of common genes. These differences in the GRNs reveal why mutations in svb are unlikely to contribute to leg trichome evolution and how instead miR-92a represents the key evolutionary switch in this context. Our work shows that variability in GRNs across different developmental contexts, as well as whether a morphological feature is lost versus gained, influence the nodes at which a GRN evolves to cause morphological change. Therefore, our findings have important implications for understanding the pathways and predictability of evolution.

Homeobox Gene Duplication and Divergence in Arachnids.

Homeobox genes are key toolkit genes that regulate the development of metazoans and changes in their regulation and copy number have contributed to the evolution of phenotypic diversity. We recently identified a whole genome duplication (WGD) event that occurred in an ancestor of spiders and scorpions (Arachnopulmonata), and that many homeobox genes, including two Hox clusters, appear to have been retained in arachnopulmonates. To better understand the consequences of this ancient WGD and the evolution of arachnid homeobox genes, we have characterized and compared the homeobox repertoires in a range of arachnids. We found that many families and clusters of these genes are duplicated in all studied arachnopulmonates (Parasteatoda tepidariorum, Pholcus phalangioides, Centruroides sculpturatus, and Mesobuthus martensii) compared with nonarachnopulmonate arachnids (Phalangium opilio, Neobisium carcinoides, Hesperochernes sp., and Ixodes scapularis). To assess divergence in the roles of homeobox ohnologs, we analyzed the expression of P. tepidariorum homeobox genes during embryogenesis and found pervasive changes in the level and timing of their expression. Furthermore, we compared the spatial expression of a subset of P. tepidariorum ohnologs with their single copy orthologs in P. opilio embryos. We found evidence for likely subfunctionlization and neofunctionalization of these genes in the spider. Overall our results show a high level of retention of homeobox genes in spiders and scorpions post-WGD, which is likely to have made a major contribution to their developmental evolution and diversification through pervasive subfunctionlization and neofunctionalization, and paralleling the outcomes of WGD in vertebrates.

Human eye conditions: insights from the fly eye.

The fruit fly Drosophila melanogaster has served as an excellent model to study and understand the genetics of many human diseases from cancer to neurodegeneration. Studying the regulation of growth, determination and differentiation of the compound eyes of this fly, in particular, have provided key insights into a wide range of diseases. Here we review the regulation of the development of fly eyes in light of shared aspects with human eye development. We also show how understanding conserved regulatory pathways in eye development together with the application of tools for genetic screening and functional analyses makes Drosophila a powerful model to diagnose and characterize the genetics underlying many human eye conditions, such as aniridia and retinitis pigmentosa. This further emphasizes the importance and vast potential of basic research to underpin applied research including identifying and treating the genetic basis of human diseases.

The Wnt and Delta-Notch signalling pathways interact to direct pair-rule gene expression via caudal during segment addition in the spider Parasteatoda tepidariorum.

In short-germ arthropods, posterior segments are added sequentially from a segment addition zone (SAZ) during embryogenesis. Studies in spiders such as Parasteatoda tepidariorum have provided insights into the gene regulatory network (GRN) underlying segment addition, and revealed that Wnt8 is required for dynamic Delta (Dl) expression associated with the formation of new segments. However, it remains unclear how these pathways interact during SAZ formation and segment addition. Here, we show that Delta-Notch signalling is required for Wnt8 expression in posterior SAZ cells, but represses the expression of this Wnt gene in anterior SAZ cells. We also found that these two signalling pathways are required for the expression of the spider orthologues of even-skipped (eve) and runt-1 (run-1), at least in part via caudal (cad). Moreover, it appears that dynamic expression of eve in this spider does not require a feedback loop with run-1, as is found in the pair-rule circuit of the beetle Tribolium Taken together, our results suggest that the development of posterior segments in Parasteatoda is directed by dynamic interactions between Wnt8 and Delta-Notch signalling that are read out by cad, which is necessary but probably not sufficient to regulate the expression of eve and run-1 Our study therefore provides new insights towards better understanding the evolution and developmental regulation of segmentation in other arthropods, including insects.

Duplication and expression of Sox genes in spiders.

BACKGROUND: The Sox family of transcription factors is an important part of the genetic 'toolbox' of all metazoans examined to date and is known to play important developmental roles in vertebrates and insects. However, outside the commonly studied Drosophila model little is known about the repertoire of Sox family transcription factors in other arthropod species. Here we characterise the Sox family in two chelicerate species, the spiders Parasteatoda tepidariorum and Stegodyphus mimosarum, which have experienced a whole genome duplication (WGD) in their evolutionary history. RESULTS: We find that virtually all of the duplicate Sox genes have been retained in these spiders after the WGD. Analysis of the expression of Sox genes in P. tepidariorum embryos suggests that it is likely that some of these genes have neofunctionalised after duplication. Our expression analysis also strengthens the view that an orthologue of vertebrate Group B1 genes, SoxNeuro, is implicated in the earliest events of CNS specification in both vertebrates and invertebrates. In addition, a gene in the Dichaete/Sox21b class is dynamically expressed in the spider segment addition zone, suggestive of an ancient regulatory mechanism controlling arthropod segmentation as recently suggested for flies and beetles. Together with the recent analysis of Sox gene expression in the embryos of other arthropods, our findings support the idea of conserved functions for some of these genes, including a potential role for SoxC and SoxD genes in CNS development and SoxF in limb development. CONCLUSIONS: Our study provides a new chelicerate perspective to understanding the evolution and function of Sox genes and how the retention of duplicates of such important tool-box genes after WGD has contributed to different aspects of spider embryogenesis. Future characterisation of the function of these genes in spiders will help us to better understand the evolution of the regulation of important developmental processes in arthropods and other metazoans including neurogenesis and segmentation.

Molecular control of gut formation in the spider Parasteatoda tepidariorum.

The development of a digestive system is an essential feature of bilaterians. Studies of the molecular control of gut formation in arthropods have been studied in detail in the fruit fly Drosophila melanogaster. However, little is known in other arthropods, especially in noninsect arthropods. To better understand the evolution of arthropod alimentary system, we investigate the molecular control of gut development in the spider Parasteatoda tepidariorum (Pt), the primary chelicerate model species for developmental studies. Orthologs of the ectodermal genes Pt-wingless (Pt-wg) and Pt-hedgehog (Pt-hh), of the endodermal genes, Pt-serpent (Pt-srp) and Pt-hepatocyte-nuclear factor-4 (Pt-hnf4) and of the mesodermal gene Pt-twist (Pt-twi) are expressed in the same germ layers during spider gut development as in D. melanogaster. Thus, our expression data suggest that the downstream molecular components involved in gut development in arthropods are conserved. However, Pt-forkhead (Pt-fkh) expression and function in spiders is considerably different from its D. melanogaster ortholog. Pt-fkh is expressed before gastrulation in a cell population that gives rise to endodermal and mesodermal precursors, suggesting a possible role for this factor in specification of both germ layers. To test this hypothesis, we knocked down Pt-fkh via RNA interference. Pt-fkh RNAi embryos not only fail to develop a proper gut, but also lack the mesodermal Pt-twi expressing cells. Thus, in spiders Pt-fkh specifies endodermal and mesodermal germ layers. We discuss the implications of these findings for the evolution and development of gut formation in Ecdysozoans.

Morphological evolution caused by many subtle-effect substitutions in regulatory DNA.

Morphology evolves often through changes in developmental genes, but the causal mutations, and their effects, remain largely unknown. The evolution of naked cuticle on larvae of Drosophila sechellia resulted from changes in five transcriptional enhancers of shavenbaby (svb), a transcript of the ovo locus that encodes a transcription factor that governs morphogenesis of microtrichiae, hereafter called 'trichomes'. Here we show that the function of one of these enhancers evolved through multiple single-nucleotide substitutions that altered both the timing and level of svb expression. The consequences of these nucleotide substitutions on larval morphology were quantified with a novel functional assay. We found that each substitution had a relatively small phenotypic effect, and that many nucleotide changes account for this large morphological difference. In addition, we observed that the substitutions had non-additive effects. These data provide unprecedented resolution of the phenotypic effects of substitutions and show how individual nucleotide changes in a transcriptional enhancer have caused morphological evolution.

A robust (re-)annotation approach to generate unbiased mapping references for RNA-seq-based analyses of differential expression across closely related species.

BACKGROUND: RNA-seq based on short reads generated by next generation sequencing technologies has become the main approach to study differential gene expression. Until now, the main applications of this technique have been to study the variation of gene expression in a whole organism, tissue or cell type under different conditions or at different developmental stages. However, RNA-seq also has a great potential to be used in evolutionary studies to investigate gene expression divergence in closely related species. RESULTS: We show that the published genomes and annotations of the three closely related Drosophila species D. melanogaster, D. simulans and D. mauritiana have limitations for inter-specific gene expression studies. This is due to missing gene models in at least one of the genome annotations, unclear orthology assignments and significant gene length differences in the different species. A comprehensive evaluation of four statistical frameworks (DESeq2, DESeq2 with length correction, RPKM-limma and RPKM-voom-limma) shows that none of these methods sufficiently accounts for inter-specific gene length differences, which inevitably results in false positive candidate genes. We propose that published reference genomes should be re-annotated before using them as references for RNA-seq experiments to include as many genes as possible and to account for a potential length bias. We present a straight-forward reciprocal re-annotation pipeline that allows to reliably compare the expression for nearly all genes annotated in D. melanogaster. CONCLUSIONS: We conclude that our reciprocal re-annotation of previously published genomes facilitates the analysis of significantly more genes in an inter-specific differential gene expression study. We propose that the established pipeline can easily be applied to re-annotate other genomes of closely related animals and plants to improve comparative expression analyses.

From shavenbaby to the naked valley: trichome formation as a model for evolutionary developmental biology.

Microtrichia or trichomes are non-sensory actin protrusions produced by the epidermal cells of many insects. Studies of trichome formation in Drosophila have over the last 30 years provided key insights towards our understanding of gene regulation, gene regulatory networks (GRNs), development, the genotype to phenotype map, and the evolution of these processes. Here we review classic studies that have used trichome formation as a model to shed light on Drosophila development as well as recent research on the architecture of the GRN underlying trichome formation. This includes the findings that both small peptides and microRNAs play important roles in the regulation and evolution of this network. In addition, we review research on the evolution of trichome patterns that has provided novel insights into the function and architecture of cis-regulatory modules, and into the genetic basis of morphological change. We conclude that further research on these apparently simple and often functionally enigmatic structures will continue to provide new and important knowledge about development and evolution.

Evolutionary crossroads in developmental biology: the spider Parasteatoda tepidariorum.

Spiders belong to the chelicerates, which is an arthropod group that branches basally from myriapods, crustaceans and insects. Spiders are thus useful models with which to investigate whether aspects of development are ancestral or derived with respect to the arthropod common ancestor. Moreover, they serve as an important reference point for comparison with the development of other metazoans. Therefore, studies of spider development have made a major contribution to advancing our understanding of the evolution of development. Much of this knowledge has come from studies of the common house spider, Parasteatoda tepidariorum. Here, we describe how the growing number of experimental tools and resources available to study Parasteatoda development have provided novel insights into the evolution of developmental regulation and have furthered our understanding of metazoan body plan evolution.

Sexual dimorphism and natural variation within and among species in the Drosophila retinal mosaic.

BACKGROUND: Insect compound eyes are composed of ommatidia, which contain photoreceptor cells that are sensitive to different wavelengths of light defined by the specific rhodopsin proteins that they express. The fruit fly Drosophila melanogaster has several different ommatidium types that can be localised to specific retinal regions, such as the dorsal rim area (DRA), or distributed stochastically in a mosaic across the retina, like the 'pale' and 'yellow' types. Variation in these ommatidia patterns very likely has important implications for the vision of insects and could underlie behavioural and environmental adaptations. However, despite the detailed understanding of ommatidia specification in D. melanogaster, the extent to which the frequency and distribution of the different ommatidium types vary between sexes, strains and species of Drosophila is not known. RESULTS: We investigated the frequency and distribution of ommatidium types based on rhodopsin protein expression, and the expression levels of rhodopsin transcripts in the eyes of both sexes of different strains of D. melanogaster, D. simulans and D. mauritiana. We found that while the number of DRA ommatidia was invariant, Rh3 expressing ommatidia were more frequent in the larger eyes of females compared to the males of all species analysed. The frequency and distribution of ommatidium types also differed between strains and species. The D. simulans strain ZOM4 has the highest frequency of Rh3 expressing ommatidia, which is associated with a non-stochastic patch of pale and odd-coupled ommatidia in the dorsal-posterior of their eyes. CONCLUSIONS: Our results show that there is striking variation in the frequency and distribution of ommatidium types between sexes, strains and species of Drosophila. This suggests that evolutionary changes in the underlying regulatory mechanisms can alter the distribution of ommatidium types to promote or restrict their expression in specific regions of the eye within and between species, and that this could cause differences in vision among these flies.

Pleiotropic scaling and QTL data.

Wagner et al. have recently introduced much-needed data to the debate on how complexity of the genotype-phenotype map affects the distribution of mutational effects. They used quantitative trait loci (QTLs) mapping analysis of 70 skeletal characters in mice and regressed the total QTL effect on the number of traits affected (level of pleiotropy). From their results they suggest that mutations with higher pleiotropy have a larger effect, on average, on each of the affected traits-a surprising finding that contradicts previous models. We argue that the possibility of some QTL regions containing multiple mutations, which was not considered by the authors, introduces a bias that can explain the discrepancy between one of the previously suggested models and the new data.