Vanishing chromosomal inversion clines in Drosophila subobscura from Chile: Is behavioral thermoregulation to blame?

Chromosomal inversion clines paralleling the long-standing ones in native Palearctic populations of Drosophila subobscura evolved swiftly after this species invaded the Americas in the late 1970s and early 1980s. However, the new clines did not consistently continue to converge on the Old World baseline. Our recent survey of Chilean populations of D. subobscura shows that inversion clines have faded or even changed sign with latitude. Here, we investigate the hypothesis that this fading of inversion clines might be due to the Bogert effect, namely, that flies' thermoregulatory behavior has eventually compensated for environmental variation in temperature, thus buffering selection on thermal-related traits. We show that latitudinal divergence in thermal preference (T-p) has evolved in Chile for females, with higher-latitude flies having a lower mean T-p. Plastic responses in T-p also lessen latitudinal thermal variation because flies developed at colder temperatures prefer warmer microclimates. Our results are consistent with the idea that active behavioral thermoregulation might buffer environmental variation and reduce the potential effect of thermal selection on other traits as chromosomal arrangements.

Gene flow and gene flux shape evolutionary patterns of variation in Drosophila subobscura

Gene flow (defined as allele exchange between populations) and gene flux (defined as allele exchange during meiosis in heterokaryotypic females) are important factors decreasing genetic differentiation between populations and inversions. Many chromosomal inversions are under strong selection and their role in recombination reduction enhances the maintenance of their genetic distinctness. Here we analyze levels and patterns of nucleotide diversity, selection and demographic history, using 37 individuals of Drosophila subobscura from Mount Parnes (Greece) and Barcelona (Spain). Our sampling focused on two frequent O-chromosome arrangements that differ by two overlapping inversions (OST and O3+4), which are differentially adapted to the environment as observed by their opposing latitudinal clines in inversion frequencies. The six analyzed genes (Pif1A, Abi, Sqd, Yrt, Atpa and Fmr1) were selected for their location across the O-chromosome and their implication in thermal adaptation. Despite the extensive gene flux detected outside the inverted region, significant genetic differentiation between both arrangements was found inside it. However, high levels of gene flow were detected for all six genes when comparing the same arrangement among populations. These results suggest that the adaptive value of inversions is maintained, regardless of the lack of genetic differentiation within arrangements from different populations, and thus favors the Local Adaptation hypothesis over the Coadapted Genome hypothesis as the basis of the selection acting on inversions in these populations.

Inversion polymorphism in two Serbian natural populations of Drosophila subobscura: Analysis of long-term changes

To study whether inversions (or arrangements) by themselves or karyotypes are the main global warming adaptive target of natural selection, two Drosophila subobscura Serbian populations (Apatin and Petnica) were re analyzed using different statistical approaches. Both populations were sampled in an approximately 15 years period: Apatin in 1994 and 2008 + 2009 and Petnica in 1995 and 2010. For all chromosomes, the four collections studied were in Hardy-Weinberg equilibrium. Thus, it seems that inversions (or arrangements) combined at random to constitute populations" karyotypes. However, there were differences in karyotypic fre quencies along the years, although they were significant only for Apatin population. It is possible to conclude that inversions (or arrangements) are likely the target of natural selection, because they presented long term changes, but combine at random to generate the corresponding karyotypic combinations. As a consequence, the frequencies of karyotypes also change along time.

Family Size Evolution in Drosophila Chemosensory Gene Families: A Comparative Analysis with a Critical Appraisal of Methods

Gene turnover rates and the evolution of gene family sizes are important aspects of genome evolution. Here, we use curated sequence data of the major chemosensory gene families from Drosophila-the gustatory receptor, odorant receptor, ionotropic receptor, and odorant-binding protein families-to conduct a comparative analysis among families, exploring different methods to estimate gene birth and death rates, including an ad hoc simulation study. Remarkably, we found that the state-of-the-art methods may produce very different rate estimates, which may lead to disparate conclusions regarding the evolution of chemosensory gene family sizes in Drosophila. Among biological factors, we found that a peculiarity of D. sechellia's gene turnover rates was a major source of bias in global estimates, whereas gene conversion had negligible effects for the families analyzed herein. Turnover rates vary considerably among families, subfamilies, and ortholog groups although all analyzed families were quite dynamic in terms of gene turnover. Computer simulations showed that the methods that use ortholog group information appear to be the most accurate for the Drosophila chemosensory families. Most importantly, these results reveal the potential of rate heterogeneity among lineages to severely bias some turnover rate estimation methods and the need of further evaluating the performance of these methods in a more diverse sampling of gene families and phylogenetic contexts. Using branch-specific codon substitution models, we find further evidence of positive selection in recently duplicated genes, which attests to a nonneutral aspect of the gene birth-and-death process.

Patterns of nucleotide diversity at the regions encompassing the Drosophila insulin-like peptide (dilp) genes: demography vs positive selection in Drosophila melanogaster.

In Drosophila, the insulin-signaling pathway controls some life history traits, such as fertility and lifespan, and it is considered to be the main metabolic pathway involved in establishing adult body size. Several observations concerning variation in body size in the Drosophila genus are suggestive of its adaptive character. Genes encoding proteins in this pathway are, therefore, good candidates to have experienced adaptive changes and to reveal the footprint of positive selection. The Drosophila insulin-like peptides (DILPs) are the ligands that trigger the insulin-signaling cascade. In Drosophila melanogaster, there are several peptides that are structurally similar to the single mammalian insulin peptide. The footprint of recent adaptive changes on nucleotide variation can be unveiled through the analysis of polymorphism and divergence. With this aim, we have surveyed nucleotide sequence variation at the dilp1-7 genes in a natural population of D. melanogaster. The comparison of polymorphism in D. melanogaster and divergence from D. simulans at different functional classes of the dilp genes provided no evidence of adaptive protein evolution after the split of the D. melanogaster and D. simulans lineages. However, our survey of polymorphism at the dilp gene regions of D. melanogaster has provided some evidence for the action of positive selection at or near these genes. The regions encompassing the dilp1-4 genes and the dilp6 gene stand out as likely affected by recent adaptive events.

Gene expression following induction of regeneration in Drosophila wing imaginal discs

Background: Regeneration is the ability of an organism to rebuild a body part that has been damaged or amputated, and can be studied at the molecular level using model organisms. Drosophila imaginal discs, which are the larval primordia of adult cuticular structures, are capable of undergoing regenerative growth after transplantation and in vivo culture into the adult abdomen. Results: Using expression profile analyses, we studied the regenerative behaviour of wing discs at 0, 24 and 72 hours after fragmentation and implantation into adult females. Based on expression level, we generated a catalogue of genes with putative role in wing disc regeneration, identifying four classes: 1) genes with differential expression within the first 24 hours; 2) genes with differential expression between 24 and 72 hours; 3) genes that changed significantly in expression levels between the two time periods; 4) genes with a sustained increase or decrease in their expression levels throughout regeneration. Among these genes, we identified members of the JNK and Notch signalling pathways and chromatin regulators. Through computational analysis, we recognized putative binding sites for transcription factors downstream of these pathways that are conserved in multiple Drosophilids, indicating a potential relationship between members of the different gene classes. Experimental data from genetic mutants provide evidence of a requirement of selected genes in wing disc regeneration. Conclusions: We have been able to distinguish various classes of genes involved in early and late steps of the regeneration process. Our data suggests the integration of signalling pathways in the promoters of regulated genes.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background: The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians. Results: In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern. Conclusions: The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background: The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians. Results: In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern. Conclusions: The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.

Gene expression following induction of regeneration in Drosophila wing imaginal discs

Background: Regeneration is the ability of an organism to rebuild a body part that has been damaged or amputated, and can be studied at the molecular level using model organisms. Drosophila imaginal discs, which are the larval primordia of adult cuticular structures, are capable of undergoing regenerative growth after transplantation and in vivo culture into the adult abdomen. Results: Using expression profile analyses, we studied the regenerative behaviour of wing discs at 0, 24 and 72 hours after fragmentation and implantation into adult females. Based on expression level, we generated a catalogue of genes with putative role in wing disc regeneration, identifying four classes: 1) genes with differential expression within the first 24 hours; 2) genes with differential expression between 24 and 72 hours; 3) genes that changed significantly in expression levels between the two time periods; 4) genes with a sustained increase or decrease in their expression levels throughout regeneration. Among these genes, we identified members of the JNK and Notch signalling pathways and chromatin regulators. Through computational analysis, we recognized putative binding sites for transcription factors downstream of these pathways that are conserved in multiple Drosophilids, indicating a potential relationship between members of the different gene classes. Experimental data from genetic mutants provide evidence of a requirement of selected genes in wing disc regeneration. Conclusions: We have been able to distinguish various classes of genes involved in early and late steps of the regeneration process. Our data suggests the integration of signalling pathways in the promoters of regulated genes.

CBS: an open platform that integrates predictive methods and epigenetics information to characterize conserved regulatory features in multiple Drosophila genomes.

Background: Information about the composition of regulatory regions is of great value for designing experiments to functionally characterize gene expression. The multiplicity of available applications to predict transcription factor binding sites in a particular locus contrasts with the substantial computational expertise that is demanded to manipulate them, which may constitute a potential barrier for the experimental community. Results: CBS (Conserved regulatory Binding Sites, http://compfly.bio.ub.es/CBS) is a public platform of evolutionarily conserved binding sites and enhancers predicted in multiple Drosophila genomes that is furnished with published chromatin signatures associated to transcriptionally active regions and other experimental sources of information. The rapid access to this novel body of knowledge through a user-friendly web interface enables non-expert users to identify the binding sequences available for any particular gene, transcription factor, or genome region. Conclusions: The CBS platform is a powerful resource that provides tools for data mining individual sequences and groups of co-expressed genes with epigenomics information to conduct regulatory screenings in Drosophila.

Genetic markers of the Va/Ba balanced lethal strain of Drosophila subobscura

The Va/Ba strain, constructed by Sperlich et al. (1977), is the only balanced lethal strain in D. subobscura. It allows the production of homozygous O chromosomes and has been a useful tool not only to analyse chromosomal viabilities but also to obtain homokaryotypic lines (Mestres and Serra, 2008). Besides the morphological dominant mutations Va (Varicose) and Ba (Bare), other genetic markers have been characterized in this strain, some of them by our group and not described previously. Here we present a list of these markers.

Identificació de nous gens que participen en la remodelació del sistema traqueal de Drosophila melanogaster

Un dels organismes model més utilitzats en experimentació genètica és la Drosophila melanogaster ja que la facilitat de manipulació genètica i la seva simplicitat permeten estudiar processos biològics amb múltiples aplicabilitats en diferents àmbits d’estudi com el desenvolupament embrionari i la morfogènesis. La morfogènesi es un dels esdeveniments més importants durant el desenvolupament embrionari que permet la formació dels diferent teixits i òrgans, i que depèn de l'expressió genètica i de l'activació i coordinació de diferents vies de senyalització. Entendre com es coordinen aquest processos es fonamental per conèixer com es forma un òrgan. Així, l’objectiu principal d’aquest Treball de Final de Grau és identificar nous gens implicats en la formació del sistema traqueal (el nostre òrgan model) mitjançant un mini-­‐cribratge funcional de gens que s’expressen en la tràquea, a més de generar eines per a l'estudi de la via de senyalització FGF/Bnl durant la remodelació del sistema traqueal mitjançant la tècnica de knock in. Per a dur-­‐ho a terme, amb el suport de la base de dades de Gens i Genomes de Drosophila melanogaster (mod-­‐ENCODE Tissue Expression Data) s’han seleccionat gens candidats expressats a la tràquea en estat larvari. Un cop identificats, s'ha estudiat la seva possible funció en el desenvolupament de les tràquees mitjançant el seu silenciament amb el sistema UAS-­‐Gal4. Així hem vist que Vein (CG10491), CG17098, No Ocelli (CG4491) i Peptidasa (CG4017) presenten diversos fenotips que afecten la formació dels traqueoblasts. També hem vist que Vein, lligand de la via EGF és necessari per a la proliferació i supervivència de les cèl·∙lules traqueals del sac aeri. Finalment s’ha iniciat la generació d'un knock in en el gen branchless (bnl). Per aquest motiu s'han amplificat les regions 5’ i 3’ de l’exó 2 del gen Bnl i s'ha iniciat la seva clonació dirigida al vector de destí pTV-­‐Cherry. Aquesta tècnica generarà eines que permetran entendre la funció del gen bnl durant la remodelació del sistema traqueal.

Study of linkage between miniature and singed genes in Drosophila melanogaster

We have developed a practical exercise for undergraduate students whose main aim is to identify, using genetic crosses, a pair of D. melanogaster mutations (miniature and singed). Each student receives a vial with the problem strain containing two unknown mutations. The first step is to observe and describe both mutations. Then, the students carry out genetic crosses between mutant and normal strains: (P) ♀ mutant strain × ♂ normal strain (P) ♀ normal strain × ♂ mutant strain A different offspring is expected in these crosses: in the first one we will obtain normal females and m sn males, whereas in the second all individuals will present normal phenotype. It is possible to deduce that both are sex linked mutations. With this information and to simplify the amount of work, only F1 individuals from the first cross will be used (m+sn+ / m sn × m sn / Y chrom.) to obtain the F2 generation. By counting the number of miniature (recombinant type), singed (recombinant type), miniature-singed (parental type) and normal (parental type) flies it is possible to estimate the recombination frequency between both genes. Knowing the phenotype, their chromosomal location (X chromosome) and the genetic distance between both mutations, it is possible to identify them by finding all this information in a Drosophila melanogaster genetic map. Additionally, a statistical analysis can be carried out to compare the number of expected F2 individuals with those observed in the experiment. As the distance between both genes is 15.1 m.u., then the expected percentages for each phenotype would be: normal (42.45%), miniature-signed (42.45%), miniature (7.55%) and singed (7.55%). Multiplying the frequency of each class by the total number of individuals obtained in the F2 it is possible to estimate the expected number of flies for each class. Finally, a χ2 test can be computed to ascertain whether there are significant differences between expected and observed number of individuals.

Inbreeding and thermal adaptation in Drosophila subobscura

Using a well-adapted Drosophila subobscura population (Avala, Serbia), a drastic experiment of inbreeding was carried out to assess whether the expected level of homozygosity could be reached or if other evolutionary forces affected the process. In general, no significant changes of inversion (or arrangement) frequencies were detected after 12 brother sister mating generations. Furthermore, no significant differences were obtained between observed and expected (under the inbreeding model) karyotypic frequencies. Thus, these results seemed to indicate that the main evolutionary factor in the experiment was inbreeding. However, in the G12 generation, complete chromosomal fixation was reached only in two out of the eight final inbred lines. In these lines, the chromosomal compositions were difficult to interpret, but they could be likely a consequence of adaptation to particular laboratory conditions (constant 18 °C, food, light period, etc.). Finally, in a second experiment, the inbred lines presented higher fertility at 18 °C than at 13 °C. Also, there was a significant line effect on fertility: inbred line number 6 (A1, J1, U1+2; U1+2+6, E8, and O3+4+7) presented the highest values, which maybe the result of an adaptation to laboratory conditions. Thus, the results obtained in our experiments reflect the adaptive potential of D. subobscura inversions.

A carboxylesterase, Esterase-6, modulates sensory physiological and behavioural response dynamics to pheromone in Drosophila

Background: Insects respond to the spatial and temporal dynamics of a pheromone plume, which implies not only a strong response to"odor on", but also to"odor off". This requires mechanisms geared toward a fast signal termination. Several mechanisms may contribute to signal termination, among which odorant-degrading enzymes. These enzymes putatively play a role in signal dynamics by a rapid inactivation of odorants in the vicinity of the sensory receptors, although direct in vivo experimental evidences are lacking. Here we verified the role of an extracellular carboxylesterase, esterase-6 (Est-6), in the sensory physiological and behavioral dynamics of Drosophila melanogaster response to its pheromone, cis-vaccenyl acetate (cVA). Est-6 was previously linked to post-mating effects in the reproductive system of females. As Est-6 is also known to hydrolyze cVA in vitro and is expressed in the main olfactory organ, the antenna, we tested here its role in olfaction as a putative odorant-degrading enzyme. Results: We first confirm that Est-6 is highly expressed in olfactory sensilla, including cVA-sensitive sensilla, and we show that expression is likely associated with non-neuronal cells. Our electrophysiological approaches show that the dynamics of olfactory receptor neuron (ORN) responses is strongly influenced by Est-6, as in Est-6° null mutants (lacking the Est-6 gene) cVA-sensitive ORN showed increased firing rate and prolonged activity in response to cVA. Est-6° mutant males had a lower threshold of behavioral response to cVA, as revealed by the analysis of two cVAinduced behaviors. In particular, mutant males exhibited a strong decrease of male-male courtship, in association with a delay in courtship initiation. Conclusions: Our study presents evidence that Est-6 plays a role in the physiological and behavioral dynamics of sex pheromone response in Drosophila males and supports a role of Est-6 as an odorant-degrading enzyme (ODE) in male antennae. Our results also expand the role of Est-6 in Drosophila biology, from reproduction to olfaction, and highlight the role of ODEs in insect olfaction. Keywords: carboxylesterase, esterase 6, olfaction, pheromone, signal termination