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.

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.

GeneID in "Drosophila"

GeneID is a program to predict genes in anonymous genomic sequences designed with a hierarchical structure. In the first step, splice sites, and start and stop codons are predicted and scored along the sequence using position weight matrices (PWMs). In the second step, exons are built from the sites. Exons are scored as the sum of the scores of the defining sites, plus the log-likelihood ratio of a Markov model for coding DNA. In the last step, from the set of predicted exons, the gene structure is assembled, maximizing the sum of the scores of the assembled exons. In this paper we describe the obtention of PWMs for sites, and the Markov model of coding DNA in Drosophila melanogaster. We also compare other models of coding DNA with the Markov model. Finally, we present and discuss the results obtained when GeneID is used to predict genes in the Adh region. These results show that the accuracy of GeneID predictions compares currently with that of other existing tools but that GeneID is likely to be more efficient in terms of speed and memory usage.

A short motif in Drosophila SECIS Binding Protein 2 provides differential binding affinity to SECIS RNA hairpins

Selenoproteins contain the amino acid selenocysteine which is encoded by a UGA Sec codon. Recoding UGA Sec requires a complex mechanism, comprising the cis-acting SECIS RNA hairpin in the 3′UTR of selenoprotein mRNAs, and trans-acting factors. Among these, the SECIS Binding Protein 2 (SBP2) is central to the mechanism. SBP2 has been so far functionally characterized only in rats and humans. In this work, we report the characterization of the Drosophila melanogaster SBP2 (dSBP2). Despite its shorter length, it retained the same selenoprotein synthesis-promoting capabilities as the mammalian counterpart. However, a major difference resides in the SECIS recognition pattern: while human SBP2 (hSBP2) binds the distinct form 1 and 2 SECIS RNAs with similar affinities, dSBP2 exhibits high affinity toward form 2 only. In addition, we report the identification of a K (lysine)-rich domain in all SBP2s, essential for SECIS and 60S ribosomal subunit binding, differing from the well-characterized L7Ae RNA-binding domain. Swapping only five amino acids between dSBP2 and hSBP2 in the K-rich domain conferred reversed SECIS-binding properties to the proteins, thus unveiling an important sequence for form 1 binding.

Conserved chromosomal clustering of genes governed by chromatin regulators in Drosophila

BACKGROUND: The trithorax group (trxG) and Polycomb group (PcG) proteins are responsible for the maintenance of stable transcriptional patterns of many developmental regulators. They bind to specific regions of DNA and direct the post-translational modifications of histones, playing a role in the dynamics of chromatin structure. RESULTS: We have performed genome-wide expression studies of trx and ash2 mutants in Drosophila melanogaster. Using computational analysis of our microarray data, we have identified 25 clusters of genes potentially regulated by TRX. Most of these clusters consist of genes that encode structural proteins involved in cuticle formation. This organization appears to be a distinctive feature of the regulatory networks of TRX and other chromatin regulators, since we have observed the same arrangement in clusters after experiments performed with ASH2, as well as in experiments performed by others with NURF, dMyc, and ASH1. We have also found many of these clusters to be significantly conserved in D. simulans, D. yakuba, D. pseudoobscura and partially in Anopheles gambiae. CONCLUSION: The analysis of genes governed by chromatin regulators has led to the identification of clusters of functionally related genes conserved in other insect species, suggesting this chromosomal organization is biologically important. Moreover, our results indicate that TRX and other chromatin regulators may act globally on chromatin domains that contain transcriptionally co-regulated genes.

Conserved chromosomal clustering of genes governed by chromatin regulators in Drosophila

Background: The trithorax group (trxG) and Polycomb group (PcG) proteins are responsible for the maintenance of stable transcriptional patterns of many developmental regulators. They bind to specific regions of DNA and direct the post-translational modifications of histones, playing a role in the dynamics of chromatin structure.Results: We have performed genome-wide expression studies of trx and ash2 mutants in Drosophila melanogaster. Using computational analysis of our microarray data, we have identified 25 clusters of genes potentially regulated by TRX. Most of these clusters consist of genes that encode structural proteins involved in cuticle formation. This organization appears to be a distinctive feature of the regulatory networks of TRX and other chromatin regulators, since we have observed the same arrangement in clusters after experiments performed with ASH2, as well as in experiments performed by others with NURF, dMyc, and ASH1. We have also found many of these clusters to be significantly conserved in D. simulans, D. yakuba, D. pseudoobscura and partially in Anopheles gambiae.Conclusion: The analysis of genes governed by chromatin regulators has led to the identification of clusters of functionally related genes conserved in other insect species, suggesting this chromosomal organization is biologically important. Moreover, our results indicate that TRX and other chromatin regulators may act globally on chromatin domains that contain transcriptionally co-regulated genes.

Odorant Receptor (Or) genes: polymorphism and divergence in the D. melanogaster and D. pseudoobscura Lineages

Background: In insects, like in most invertebrates, olfaction is the principal sensory modality, which provides animals with essential information for survival and reproduction. Odorant receptors are involved in this response, mediating interactions between an individual and its environment, as well as between individuals of the same or different species. The adaptive importance of odorant receptors renders them good candidates for having their variation shaped by natural selection. Methodology/Principal Findings: We analyzed nucleotide variation in a subset of eight Or genes located on the 3L chromosomal arm of Drosophila melanogaster in a derived population of this species and also in a population of Drosophila pseudoobscura. Some heterogeneity in the silent polymorphism to divergence ratio was detected in the D. melanogaster/D. simulans comparison, with a single gene (Or67b) contributing ~37% to the test statistic. However, no other signals of a very recent selective event were detected at this gene. In contrast, at the speciation timescale, the MK test uncovered the footprint of positive selection driving the evolution of two of the encoded proteins in both D. melanogaster ¿OR65c and OR67a ¿and D. pseudoobscura ¿OR65b1 and OR67c. Conclusions: The powerful polymorphism/divergence approach provided evidence for adaptive evolution at a rather high proportion of the Or genes studied after relatively recent speciation events. It did not provide, however, clear evidence for very recent selective events in either D. melanogaster or D. pseudoobscura.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

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

Fusion of the human gene for the polyubiquitination co-effector UEV-1 with Kua, a newly identified gene

UEV proteins are enzymatically inactive variants of the E2 ubiquitin-conjugating enzymes that regulate noncanonical elongation of ubiquitin chains. In Saccharomyces cerevisiae, UEV is part of the RAD6-mediated error-free DNA repair pathway. In mammalian cells, UEV proteins can modulate c-FOS transcription and the G2-M transition of the cell cycle. Here we show that the UEV genes from phylogenetically distant organisms present a remarkable conservation in their exon–intron structure. We also show that the human UEV1 gene is fused with the previously unknown gene Kua. In Caenorhabditis elegans and Drosophila melanogaster, Kua and UEV are in separated loci, and are expressed as independent transcripts and proteins. In humans, Kua and UEV1 are adjacent genes, expressed either as separate transcripts encoding independent Kua and UEV1 proteins, or as a hybrid Kua–UEV transcript, encoding a two-domain protein. Kua proteins represent a novel class of conserved proteins with juxtamembrane histidine-rich motifs. Experiments with epitope-tagged proteins show that UEV1A is a nuclear protein, whereas both Kua and Kua–UEV localize to cytoplasmic structures, indicating that the Kua domain determines the cytoplasmic localization of Kua–UEV. Therefore, the addition of a Kua domain to UEV in the fused Kua–UEV protein confers new biological properties to this regulator of variant polyubiquitination.[Kua cDNAs isolated by RT-PCR and described in this paper have been deposited in the GenBank data library under accession nos. AF1155120 (H. sapiens) and AF152361 (D. melanogaster). Genomic clones containing UEV genes: S. cerevisiae, YGL087c (accession no. Z72609); S. pombe, c338 (accession no. AL023781); P. falciparum, MAL3P2 (accession no. AL034558); A. thaliana, F26F24 (accession no. AC005292); C. elegans, F39B2 (accession no. Z92834); D. melanogaster, AC014908; and H. sapiens, 1185N5 (accession no. AL034423). Accession numbers for Kua cDNAs in GenBank dbEST: M. musculus, AA7853; T. cruzi, AI612534. Other Kua-containing sequences: A. thaliana genomic clones F10M23 (accession no. AL035440), F19K23 (accession no. AC000375), and T20K9 (accession no. AC004786).

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

EGASP: the human ENCODE Genome Annotation Assessment Project

Background: We present the results of EGASP, a community experiment to assess the state-ofthe-art in genome annotation within the ENCODE regions, which span 1% of the human genomesequence. The experiment had two major goals: the assessment of the accuracy of computationalmethods to predict protein coding genes; and the overall assessment of the completeness of thecurrent human genome annotations as represented in the ENCODE regions. For thecomputational prediction assessment, eighteen groups contributed gene predictions. Weevaluated these submissions against each other based on a ‘reference set’ of annotationsgenerated as part of the GENCODE project. These annotations were not available to theprediction groups prior to the submission deadline, so that their predictions were blind and anexternal advisory committee could perform a fair assessment.Results: The best methods had at least one gene transcript correctly predicted for close to 70%of the annotated genes. Nevertheless, the multiple transcript accuracy, taking into accountalternative splicing, reached only approximately 40% to 50% accuracy. At the coding nucleotidelevel, the best programs reached an accuracy of 90% in both sensitivity and specificity. Programsrelying on mRNA and protein sequences were the most accurate in reproducing the manuallycurated annotations. Experimental validation shows that only a very small percentage (3.2%) of the selected 221 computationally predicted exons outside of the existing annotation could beverified.Conclusions: This is the first such experiment in human DNA, and we have followed thestandards established in a similar experiment, GASP1, in Drosophila melanogaster. We believe theresults presented here contribute to the value of ongoing large-scale annotation projects and shouldguide further experimental methods when being scaled up to the entire human genome sequence.

Comparative transcriptomics of early dipteran development

BACKGROUND: Modern sequencing technologies have massively increased the amount of data available for comparative genomics. Whole-transcriptome shotgun sequencing (RNA-seq) provides a powerful basis for comparative studies. In particular, this approach holds great promise for emerging model species in fields such as evolutionary developmental biology (evo-devo). RESULTS: We have sequenced early embryonic transcriptomes of two non-drosophilid dipteran species: the moth midge Clogmia albipunctata, and the scuttle fly Megaselia abdita. Our analysis includes a third, published, transcriptome for the hoverfly Episyrphus balteatus. These emerging models for comparative developmental studies close an important phylogenetic gap between Drosophila melanogaster and other insect model systems. In this paper, we provide a comparative analysis of early embryonic transcriptomes across species, and use our data for a phylogenomic re-evaluation of dipteran phylogenetic relationships. CONCLUSIONS: We show how comparative transcriptomics can be used to create useful resources for evo-devo, and to investigate phylogenetic relationships. Our results demonstrate that de novo assembly of short (Illumina) reads yields high-quality, high-coverage transcriptomic data sets. We use these data to investigate deep dipteran phylogenetic relationships. Our results, based on a concatenation of 160 orthologous genes, provide support for the traditional view of Clogmia being the sister group of Brachycera (Megaselia, Episyrphus, Drosophila), rather than that of Culicomorpha (which includes mosquitoes and blackflies).

Toxicity of allyl esters in insect cell lines and in Spodoptera littoralis larvae

We investigated the effects of five allyl esters, two aromatic (allyl cinnamate and allyl 2-furoate) and three aliphatic (allyl hexanoate, allyl heptanoate, and allyl octanoate) in established insect cell lines derived from different species and tissues. We studied embryonic cells of the fruit fly Drosophila melanogaster (S2) (Diptera) and the beet armyworm Spodoptera exigua (Se4) (Lepidoptera), fat body cells of the Colorado potato beetle Leptinotarsa decemlineata (CPB) (Coleoptera), ovarian cells of the silkmoth Bombyx mori (Bm5), and midgut cells of the spruce budworm Choristoneura fumiferana (CF203) (Lepidoptera). Cytotoxicity was determined with use of MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide] and trypan blue. In addition, we tested the entomotoxic action of allyl cinnamate against the cotton leafworm Spodoptera littoralis .The median (50%) cytotoxic concentrations (EC50s) of the five allyl esters in the MTT bioassays ranged between 0.25 and 27 mM with significant differences among allyl esters (P = 0.0012), cell lines (P < 0.0001), and the allyl estercell line interaction (P < 0.0001). Allyl cinnamate was the most active product, and CF203 the most sensitive cell line. In the trypan blue bioassays, cytotoxicity was produced rapidly and followed the same trend observed in the MTT bioassay. In first instars of S. littoralis, allyl cinnamate killed all larvae at 0.25% in the diet after 1 day, while this happened in third instars after 5 days. The LC50 in first instars was 0.08%. In addition, larval weight gain was reduced (P < 0.05) after 1 day of feeding on diet with 0.05%. In conclusion, the data provide evidence of the significant but differential cytotoxicity among allyl esters in insect cells of different species and tissues. Midgut cells show high sensitivity, indicating the insect midgut as a primary target tissue. Allyl cinnamate caused rapid toxic effects in S. littoralis larvae at low concentrations, suggesting further potential for use in pest control.