Cross-Species Interaction between Rapidly Evolving Telomere-Specific Drosophila Proteins

Telomere integrity in Drosophila melanogaster is maintained by a putative multisubunit complex called terminin that is believed to act in analogy to the mammalian shelterin complex in protecting chromosome ends from being recognized as sites of DNA damage. The five proteins supposed to form the terminin complex are HP1-ORC associated protein, HP1-HOAP interacting protein, Verrocchio, Drosophila Telomere Loss/Modigliani and Heterochromatic Protein 1. Four of these proteins evolve rapidly within the Drosophila genus. The accelerated evolution of terminin components may indicate the involvement of these proteins in the process by which new species arise, as the resulting divergence of terminin proteins might prevent hybrid formation, thus driving speciation. However, terminin is not an experimentally proven entity, and no biochemical studies have been performed to investigate its assembly and action in detail. Motivated by these facts in order to initiate biochemical studies on terminin function, we attempted to reconstitute terminin by co-expressing its subunits in bacteria and investigated the possible role of the fast-evolving parts of terminin components in complex assembly. Our results suggest formation of stable subcomplexes of terminin, but not of the whole complex in vitro. We found that the accelerated evolution is restricted to definable regions of terminin components, and that the divergence of D. melanogaster Drosophila Telomere Loss and D. yakuba Verrocchio proteins does not preclude their stable interaction.     

Meiotic,genomic and evolutionary properties of crossover distribution in Drosophila yakuba

The number and location of crossovers across genomes are highly regulated during meiosis, yet the key components controlling them are fast evolving, hindering our understanding of the mechanistic causes and evolutionary consequences of changes in crossover rates. Drosophila melanogaster has been a model species to study meiosis for more than a century, with an available high-resolution crossover map that is, nonetheless, missing for closely related species, thus preventing evolutionary context. Here, we applied a novel and highly efficient approach to generate whole-genome high-resolution crossover maps in D. yakuba to tackle multiple questions that benefit from being addressed collectively within an appropriate phylogenetic framework, in our case the D. melanogaster species subgroup. The genotyping of more than 1,600 individual meiotic events allowed us to identify several key distinct properties relative to D. melanogaster. We show that D. yakuba, in addition to higher crossover rates than D. melanogaster, has a stronger centromere effect and crossover assurance than any Drosophila species analyzed to date. We also report the presence of an active crossover-associated meiotic drive mechanism for the X chromosome that results in the preferential inclusion in oocytes of chromatids with crossovers. Our evolutionary and genomic analyses suggest that the genome-wide landscape of crossover rates in D. yakuba has been fairly stable and captures a significant signal of the ancestral crossover landscape for the whole D. melanogaster subgroup, even informative for the D. melanogaster lineage. Contemporary crossover rates in D. melanogaster, on the other hand, do not recapitulate ancestral crossovers landscapes. As a result, the temporal stability of crossover landscapes observed in D. yakuba makes this species an ideal system for applying population genetic models of selection and linkage, given that these models assume temporal constancy in linkage effects. Our studies emphasize the importance of generating multiple high-resolution crossover rate maps within a coherent phylogenetic context to broaden our understanding of crossover control during meiosis and to improve studies on the evolutionary consequences of variable crossover rates across genomes and time.     

Human Impacts Flatten Rainforest-Savanna Gradient and Reduce Adaptive Diversity in a Rainforest Bird

Ecological gradients have long been recognized as important regions for diversification and speciation. However, little attention has been paid to the evolutionary consequences or conservation implications of human activities that fundamentally change the environmental features of such gradients. Here we show that recent deforestation in West Africa has homogenized the rainforest-savanna gradient, causing a loss of adaptive phenotypic diversity in a common rainforest bird, the little greenbul (Andropadus virens). Previously, this species was shown to exhibit morphological and song divergence along this gradient in Central Africa. Using satellite-based estimates of forest cover, recent morphological data, and historical data from museum specimens collected prior to widespread deforestation, we show that the gradient has become shallower in West Africa and that A. virens populations there have lost morphological variation in traits important to fitness. In contrast, we find no loss of morphological variation in Central Africa where there has been less deforestation and gradients have remained more intact. While rainforest deforestation is a leading cause of species extinction, the potential of deforestation to flatten gradients and inhibit rainforest diversification has not been previously recognized. More deforestation will likely lead to further flattening of the gradient and loss of diversity, and may limit the ability of species to persist under future environmental conditions.     

Current problems of environmental gradients and species response curves in relation to continuum theory

Abstract. Comparisons of the positions of species on Grimes'C-S-R triangular ordination model with their responses to individual environmental gradients indicates that the C-S-R model does not necessarily predict species ecological behaviour. The importance of the stress, productivity and disturbance gradients relative to other environmental gradients needs to be determined. In studies of species behaviour along a biomass/productivity gradient the collective vegetation property, biomass, has been confused with the environmental factor, fertility. Patterns of responses to biomass gradients e.g. Keddy's centrifugal model, should be examined in a two-dimensional environmental space to avoid such confounding effects. Assumptions regarding the shapes of species responses to environmental gradients remain untested. A recent model of species response functions to environmental gradients suggested that skewed responses curves show a pattern in the direction of the skew, always with the tail towards the presumed most mesic position on the gradient. Further evidence is presented to support this model for a temperature gradient in eucalypt forest in south-eastern Australia. 21 out of 24 species tested conform to the model.     

Novel Drosophila Viruses Encode Host-Specific Suppressors of RNAi

The ongoing conflict between viruses and their hosts can drive the co-evolution between host immune genes and viral suppressors of immunity. It has been suggested that an evolutionary ‘arms race’ may occur between rapidly evolving components of the antiviral RNAi pathway of Drosophila and viral genes that antagonize it. We have recently shown that viral protein 1 (VP1) of Drosophila melanogaster Nora virus (DmelNV) suppresses Argonaute-2 (AGO2)-mediated target RNA cleavage (slicer activity) to antagonize antiviral RNAi. Here we show that viral AGO2 antagonists of divergent Nora-like viruses can have host specific activities. We have identified novel Nora-like viruses in wild-caught populations of D. immigrans (DimmNV) and D. subobscura (DsubNV) that are 36% and 26% divergent from DmelNV at the amino acid level. We show that DimmNV and DsubNV VP1 are unable to suppress RNAi in D. melanogaster S2 cells, whereas DmelNV VP1 potently suppresses RNAi in this host species. Moreover, we show that the RNAi suppressor activity of DimmNV VP1 is restricted to its natural host species, D. immigrans. Specifically, we find that DimmNV VP1 interacts with D. immigrans AGO2, but not with D. melanogaster AGO2, and that it suppresses slicer activity in embryo lysates from D. immigrans, but not in lysates from D. melanogaster. This species-specific interaction is reflected in the ability of DimmNV VP1 to enhance RNA production by a recombinant Sindbis virus in a host-specific manner. Our results emphasize the importance of analyzing viral RNAi suppressor activity in the relevant host species. We suggest that rapid co-evolution between RNA viruses and their hosts may result in host species-specific activities of RNAi suppressor proteins, and therefore that viral RNAi suppressors could be host-specificity factors.     

A Staging Scheme for the Development of the Moth Midge Clogmia albipunctata

Model organisms, such as Drosophila melanogaster, allow us to address a wide range of biological questions with experimental rigour. However, studies in model species need to be complemented by comparative studies if we are to fully understand the functional properties and evolutionary history of developmental processes. The establishment of new model organisms is crucial for this purpose. One of the first essential steps to establish a species as an experimental model is to carefully describe its life cycle and development. The resulting staging scheme serves as a framework for molecular studies, and allows us to homologise developmental processes between species. In this paper, we have characterised the life cycle and development of an emerging non-drosophilid dipteran model system: the moth midge Clogmia albipunctata. In particular, we focus on early embryogenesis (cleavage and blastoderm cycles before gastrulation), on formation and retraction of extraembryonic tissues, and on formation of the germ line. Considering the large evolutionary distance between the two species (approximately 250 million years), we find that the development of C. albipunctata is remarkably conserved compared to D. melanogaster. On the other hand, we detect significant differences in morphology and timing affecting the development of extraembryonic tissues and the germ line. Moreover, C. albipunctata shows several heterochronic shifts, and lacks head involution and associated processes during late stages of development.     

Improving the understanding of cytoneme-mediated morphogen gradients by in silico modeling

Morphogen gradients are crucial for the development of organisms. The biochemical properties of many morphogens prevent their extracellular free diffusion, indicating the need of an active mechanism for transport. The involvement of filopodial structures (cytonemes) has been proposed for morphogen signaling. Here, we describe an in silico model based on the main general features of cytoneme-meditated gradient formation and its implementation into Cytomorph, an open software tool. We have tested the spatial and temporal adaptability of our model quantifying Hedgehog (Hh) gradient formation in two Drosophila tissues. Cytomorph is able to reproduce the gradient and explain the different scaling between the two epithelia. After experimental validation, we studied the predicted impact of a range of features such as length, size, density, dynamics and contact behavior of cytonemes on Hh morphogen distribution. Our results illustrate Cytomorph as an adaptive tool to test different morphogen gradients and to generate hypotheses that are difficult to study experimentally.     

PRD-Class Homeobox Genes in Bovine Early Embryos: Function,Evolution, and Overlapping Roles

Eutherian Totipotent Cell Homeobox (ETCHbox) genes are mammalian-specific PRD-class homeobox genes with conserved expression in the preimplantation embryo but fast-evolving and highly divergent sequences. Here, we exploit an ectopic expression approach to examine the role of bovine ETCHbox genes and show that ARGFX and LEUTX homeodomain proteins upregulate genes normally expressed in the blastocyst; the identities of the regulated genes suggest that, in vivo, the ETCHbox genes play a role in coordinating the physical formation of the blastocyst structure. Both genes also downregulate genes expressed earlier during development and genes associated with an undifferentiated cell state, possibly via the JAK/STAT pathway. We find evidence that bovine ARGFX and LEUTX have overlapping functions, in contrast to their antagonistic roles in humans. Finally, we characterize a mutant bovine ARGFX allele which eliminates the homeodomain and show that homozygous mutants are viable. These data support the hypothesis of functional overlap between ETCHbox genes within a species, roles for ETCHbox genes in blastocyst formation and the change of their functions over evolutionary time.     

Signalogs: orthology-based identification of novel signaling pathway components in three metazoans

Background

Uncovering novel components of signal transduction pathways and their interactions within species is a central task in current biological research. Orthology alignment and functional genomics approaches allow the effective identification of signaling proteins by cross-species data integration. Recently, functional annotation of orthologs was transferred across organisms to predict novel roles for proteins. Despite the wide use of these methods, annotation of complete signaling pathways has not yet been transferred systematically between species.

Principal Findings

Here we introduce the concept of ‘signalog’ to describe potential novel signaling function of a protein on the basis of the known signaling role(s) of its ortholog(s). To identify signalogs on genomic scale, we systematically transferred signaling pathway annotations among three animal species, the nematode Caenorhabditis elegans, the fruit fly Drosophila melanogaster, and humans. Using orthology data from InParanoid and signaling pathway information from the SignaLink database, we predict 88 worm, 92 fly, and 73 human novel signaling components. Furthermore, we developed an on-line tool and an interactive orthology network viewer to allow users to predict and visualize components of orthologous pathways. We verified the novelty of the predicted signalogs by literature search and comparison to known pathway annotations. In C. elegans, 6 out of the predicted novel Notch pathway members were validated experimentally. Our approach predicts signaling roles for 19 human orthodisease proteins and 5 known drug targets, and suggests 14 novel drug target candidates.

Conclusions

Orthology-based pathway membership prediction between species enables the identification of novel signaling pathway components that we referred to as signalogs. Signalogs can be used to build a comprehensive signaling network in a given species. Such networks may increase the biomedical utilization of C. elegans and D. melanogaster. In humans, signalogs may identify novel drug targets and new signaling mechanisms for approved drugs.     

The Non-Proliferative Nature of Ascidian Folliculogenesis as a Model of Highly Ordered Cellular Topology Distinct from Proliferative Epithelia

Previous studies have addressed why and how mono‐stratified epithelia adopt a polygonal topology. One major additional, and yet unanswered question is how the frequency of different cell shapes is achieved and whether the same distribution applies between non-proliferative and proliferative epithelia. We compared different proliferative and non-proliferative epithelia from a range of organisms as well as Drosophila melanogaster mutants, deficient for apoptosis or hyperproliferative. We show that the distribution of cell shapes in non‐proliferative epithelia (follicular cells of five species of tunicates) is distinctly, and more stringently organized than proliferative ones (cultured epithelial cells and Drosophila melanogaster imaginal discs). The discrepancy is not supported by geometrical constraints (spherical versus flat monolayers), number of cells, or apoptosis events. We have developed a theoretical model of epithelial morphogenesis, based on the physics of divided media, that takes into account biological parameters such as cell‐cell contact adhesions and tensions, cell and tissue growth, and which reproduces the effects of proliferation by increasing the topological heterogeneity observed experimentally. We therefore present a model for the morphogenesis of epithelia where, in a proliferative context, an extended distribution of cell shapes (range of 4 to 10 neighbors per cell) contrasts with the narrower range of 5-7 neighbors per cell that characterizes non proliferative epithelia.     

BMP morphogen gradients in flies

Bone morphogenetic proteins (BMPs) act as morphogens to control patterning and growth in a variety of developing tissues in different species. How BMP morphogen gradients are established and interpreted in the target tissues has been extensively studied in Drosophila melanogaster. In Drosophila, Decapentaplegic (Dpp), a homologue of vertebrate BMP2/4, acts as a morphogen to control dorsal–ventral patterning of the early embryo and anterior–posterior patterning and growth of the wing imaginal disc. Despite intensive efforts over the last twenty years, how the Dpp morphogen gradient in the wing imaginal disc forms remains controversial, while gradient formation in the early embryo is well understood. In this review, we first focus on the current models of Dpp morphogen gradient formation in these two tissues, and then discuss new strategies using genome engineering and nanobodies to tackle open questions.     

Eco-Evolutionary Community Dynamics: Covariation between Diversity and Invasibility across Temperature Gradients *

Abstract Understanding biodiversity gradients is a long-standing challenge, and progress requires theory unifying ecology and evolution. Here, we unify concepts related to the speed of evolution, the influence of species richness on diversification, and niche-based coexistence. We focus on the dynamics, through evolutionary time, of community invasibility and species richness across a broad thermal gradient. In our framework, the evolution of body size influences the ecological structure and dynamics of a trophic network, and organismal metabolism ties temperature to eco-evolutionary processes. The framework distinguishes ecological invasibility (governed by ecological interactions) from evolutionary invasibility (governed by local ecology and constraints imposed by small phenotypic effects of mutation). The model yields four primary predictions: (1) ecological invasibility declines through time and with increasing temperature; (2) average evolutionary invasibility across communities increases and then decreases through time as the richness-temperature gradient flattens; (3) in the early stages of diversification, richness and evolutionary invasibility both increase with increasing temperature; and (4) at equilibrium, richness does not vary with temperature, yet evolutionary invasibility decreases with increasing temperature. These predictions emerge from the "evolutionary-speed" hypothesis, which attempts to account for latitudinal species richness gradients by invoking faster biological rates in warmer, tropical regions. The model contrasts with predictions from other richness-gradient hypotheses, such as "niche conservatism" and "species energy." Empirically testing our model's predictions should help distinguish among these hypotheses.     

Drosophila Interspecific Hybrids Phenocopy piRNA-Pathway Mutants

The Piwi-interacting RNA (piRNA) pathway defends the germline of animals from the deleterious activity of selfish transposable elements (TEs) through small-RNA mediated silencing. Adaptation to novel invasive TEs is proposed to occur by incorporating their sequences into the piRNA pool that females produce and deposit into their eggs, which then propagates immunity against specific TEs to future generations. In support of this model, the F1 offspring of crosses between strains of the same Drosophila species sometimes suffer from germline derepression of paternally inherited TE families, caused by a failure of the maternal strain to produce the piRNAs necessary for their regulation. However, many protein components of the Drosophila piRNA pathway exhibit signatures of positive selection, suggesting that they also contribute to the evolution of host genome defense. Here we investigate piRNA pathway function and TE regulation in the F1 hybrids of interspecific crosses between D. melanogaster and D. simulans and compare them with intraspecific control crosses of D. melanogaster. We confirm previous reports showing that intraspecific crosses are characterized by derepression of paternally inherited TE families that are rare or absent from the maternal genome and piRNA pool, consistent with the role of maternally deposited piRNAs in shaping TE silencing. In contrast to the intraspecific cross, we discover that interspecific hybrids are characterized by widespread derepression of both maternally and paternally inherited TE families. Furthermore, the pattern of derepression of TE families in interspecific hybrids cannot be attributed to their paucity or absence from the piRNA pool of the maternal species. Rather, we demonstrate that interspecific hybrids closely resemble piRNA effector-protein mutants in both TE misregulation and aberrant piRNA production. We suggest that TE derepression in interspecific hybrids largely reflects adaptive divergence of piRNA pathway genes rather than species-specific differences in TE-derived piRNAs.     

Comparative copulation anatomy of the Drosophila melanogaster species complex (Diptera: Drosophilidae)

The Drosophila melanogaster species complex consists of four species: D. melanogaster, D. simulans, D. sechellia and D. mauritiana. To identify these closely related species, researchers often examine the male genitalia, especially species‐specific shapes of the posterior process, as the most reliable and easily observable character. However, compared to genetic aspects, the evolutionary significance of the posterior process and other genital parts remains largely unexplained. By comparing genital coupling among these species, we revealed that the posterior processes, which are hidden under the female abdominal tergite VII when genital coupling is established, mesh with different parts of the intersegmental membrane between the tergite VIII and the oviscapts and that this membrane region broadens in a species‐specific manner. Furthermore, in D. simulans and D. sechellia, this membrane region is likely to incur wounds from the sharply pointed tip of the posterior process. On the basis of the use and functions of these and other genital parts, we discuss possible evolutionary forces underlying the diversification of genitalia in this group.     

Effects of body-size variation on flight-related traits in latitudinal populations of Drosophila melanogaster

In the present study, we tested the hypothesis whether flight-related traits such as wing area, flight-muscle ratio, wing loading and dispersal yield evidence of geographical variation in nine wild-collected as well as laboratory-reared (at 21°C) latitudinal populations of Drosophila melanogaster from the Indian subcontinent. We observed positive clinal variation in the wing–thorax ratio, wing aspect ratio and wing area, along a latitudinal gradient for both the sexes. In contrast, geographical changes in three parameters of flight ability, i.e. flight-muscle ratio, wing loading and dispersal, showed negative correlation withlatitude. On the basis of isofemale line variability, we observed positive correlation of wing loading with flight-muscle ratio as well as dispersal behaviour in both the sexes. We also found positive correlation between duration of development and wing area. Interestingly, southern populations of D. melanogaster from warm and humid habitats exhibited higher flight-muscle ratio as well as the higher wing loading than northern populations which occur in cooler and drier climatic conditions. Laboratory tests for dispersal-related walking behaviour showed significantly higher values for southern populations compared with northern populations of D. melanogaster. Multiple regression analysis of geographical changes in flight-muscle ratio, wing loading as well as walking behaviour as a function of average temperature and relative humidity of the origin of populations in wild-collected flies have suggested adaptive changes in flight-related traits in response to steeper gradients of climatic factors in the Indian subcontinent. Finally, adaptive latitudinal variations in flight-related traits in D. melanogaster are consistent with results of other studies from different continents despite differences due to specific climatic conditions in the Indian subontinent.     

Seasonal variation in life history traits in two Drosophila species

Seasonal environmental heterogeneity is cyclic, persistent and geographically widespread. In species that reproduce multiple times annually, environmental changes across seasonal time may create different selection regimes that may shape the population ecology and life history adaptation in these species. Here, we investigate how two closely related species of Drosophila in a temperate orchard respond to environmental changes across seasonal time. Natural populations of Drosophila melanogaster and Drosophila simulans were sampled at four timepoints from June through November to assess seasonal change in fundamental aspects of population dynamics as well as life history traits. D. melanogaster exhibit pronounced change across seasonal time: early in the season, the population is inferred to be uniformly young and potentially represents the early generation following overwintering survivorship. D. melanogaster isofemale lines derived from the early population and reared in a common garden are characterized by high tolerance to a variety of stressors as well as a fast rate of development in the laboratory environment that declines across seasonal time. In contrast, wild D. simulans populations were inferred to be consistently heterogeneous in age distribution across seasonal collections; only starvation tolerance changed predictably over seasonal time in a parallel manner as in D. melanogaster. These results suggest fundamental differences in population and evolutionary dynamics between these two taxa associated with seasonal heterogeneity in environmental parameters and associated selection pressures.     

Gene Circuit Analysis of the Terminal Gap Gene huckebein

The early embryo of Drosophila melanogaster provides a powerful model system to study the role of genes in pattern formation. The gap gene network constitutes the first zygotic regulatory tier in the hierarchy of the segmentation genes involved in specifying the position of body segments. Here, we use an integrative, systems-level approach to investigate the regulatory effect of the terminal gap gene huckebein (hkb) on gap gene expression. We present quantitative expression data for the Hkb protein, which enable us to include hkb in gap gene circuit models. Gap gene circuits are mathematical models of gene networks used as computational tools to extract regulatory information from spatial expression data. This is achieved by fitting the model to gap gene expression patterns, in order to obtain estimates for regulatory parameters which predict a specific network topology. We show how considering variability in the data combined with analysis of parameter determinability significantly improves the biological relevance and consistency of the approach. Our models are in agreement with earlier results, which they extend in two important respects: First, we show that Hkb is involved in the regulation of the posterior hunchback (hb) domain, but does not have any other essential function. Specifically, Hkb is required for the anterior shift in the posterior border of this domain, which is now reproduced correctly in our models. Second, gap gene circuits presented here are able to reproduce mutants of terminal gap genes, while previously published models were unable to reproduce any null mutants correctly. As a consequence, our models now capture the expression dynamics of all posterior gap genes and some variational properties of the system correctly. This is an important step towards a better, quantitative understanding of the developmental and evolutionary dynamics of the gap gene network.