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.

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.

The genetic content of chromosomal inversions across a wide latitudinal gradient

There is increasing evidence regarding the role of chromosomal inversions in relevant biological processes such as local adaptation and speciation. A classic example of the adaptive role of chromosomal polymorphisms is given by the clines of inversion frequencies in Drosophila subobscura, repeatable across continents. Nevertheless, not much is known about the molecular variation associated with these polymorphisms. We characterized the genetic content of ca. 600 individuals from nine European populations following a latitudinal gradient by analysing 19 microsatellite loci from two autosomes (J and U) and the sex chromosome (A), taking into account their chromosomal inversions. Our results clearly demonstrate the molecular genetic uniformity within a given chromosomal inversion across a large latitudinal gradient, particularly from Groningen (Netherlands) in the north to Málaga (Spain) in the south, experiencing highly diverse environmental conditions. This low genetic differentiation within the same gene arrangement across the nine European populations is consistent with the local adaptation hypothesis for th evolutionof chromosomal polymorphisms. We also show the effective role of chromosomal inversions in maintaining different genetic pools within these inverted genomic regions even in the presence of high gene flow. Inversions represent thus an important barrier to gene flux and can help maintain specific allelic combinations with positive effects on fitness. Consistent patterns of microsatellite allele-inversion linkage disequilibrium particularly in loci within inversions were also observed. Finally, we identified areas within inversions presenting clinal variation that might be under selection.

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.

Prototypical Arthropod Gene Content and Genome Organisation in the Centipede Strigamia maritima

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

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

Caracterização molecular e bioquímica de uma esterase macho-específica em Zaprionus indianus

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise de sequências da região intergência ITS-1 do rDNA em espécies de Drosophila do cluster buzzatti, (complexo Buzzatti, subgrupo Mullerri, grupo Repleta). -

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A regulação do elemento transponível Helena em Drosophila

Pós-graduação em Genética - IBILCE

Zaprionus indianus: uma visão da espécie invasora sob o aspecto genético e evolutivo

Pós-graduação em Genética - IBILCE

Investigação de elementos de transposição em populações de Drosophila mojavensis e D. arizonae e seus híbridos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)