When population and evolutionary genetics met behaviour

In this review, we analyse the impact of a population and evolutionary genetics approach on the study of insect behaviour. Our attention is focused on the model organism Drosophila melanogaster and several other insect species. In particular, we explore the relationship between rhythmic behaviours and the molecular evolution of clock and ion channel genes.

Molecular and physiological characterisation of lonotropic receptors in the sacculus in drosophila

Chemosensation is the detection of chemical signals in the environment that enable an animal to make informed decisions about food choice, mate preference or predator detection. Dissecting the molecular and neural mechanisms by which animals detect chemical cues is an important goal towards understanding how they interact with the environment. An attractive system to dissect the mechanisms of chemosensation is the olfactory system. One of the most-investigated olfactory systems is that of Drosophila melanogaster, a model organism that is amenable to a powerful combination of genetic and physiological analyses. Embedded within the antennal olfactory organ of Drosophila is an unusual sensory structure called the sacculus. The sacculus is comprised of three distinct chambers, each lined with several sensilla housing two to three neurons. Previous morphological, anatomical and surgical studies of sacculus neurons have implicated sacculus neurons in chemosensation, hygrosensation and/or thermosensation. While a subset of sacculus neurons have been physiologically characterised as temperature sensors, the role of this organ has remained largely mysterious, due to its inaccessibility to peripheral electrophysiological analysis. Recently a new family of olfactory receptors, the lonotropic Receptors (IRs), was identified. Five IRs are expressed in sacculus neurons providing the first selective molecular markers for these cells. In this thesis I describe the molecular, physiological and anatomical characterisation of these neurons. Genetic labelling of specific populations of sacculus neurons with anatomical (CD8:GFP) reporters has identified neurons in sacculus chambers I and II express IR40a+IR93a together with their co- receptor IR25a, while neurons in chamber III express IR64a with its co-receptor IR8a. Both these sets of neurons project to two distinct glomeruli in the antennal lobe; IR40a neurons project to the column and arm, IR64a neurons project to DC4 and DP1m. Through a live optical imaging screen I showed that these neurons are indeed olfactory and IR64a neurons recognise acidic ligands, while IR40a neurons recognise amine ligands. IR40a and IR64a neurons are in fact composed of anatomically and physiologically distinct subpopulations, strongly implying the existence of other factors that define their functional properties. My thesis identifies the sacculus as a specialised olfactory organ capable of detecting acids and bases, which are of widespread importance to insects. The data from my thesis along with data from other labs show the sacculus is composed of different populations of olfactory sensory neurons and thermosensory neurons. Comparative genomic analysis of sacculus IRs across insects reveals them to be among the most conserved of this receptor repertoire, suggesting that the sacculus represents an evolutionarily ancient insect olfactory acid-base sensor. - La détection des produits chimiques se trouvant dans l'environnement (perception chimiosensorielle) permet à un animal de choisir sa nourriture, son partenaire ou encore d'identifier ses prédateurs. Décortiquer les mécanismes moléculaires et neuronaux grâce auxquels les animaux détectent ces signaux chimiques permet de comprendre comment ces animaux interagissent avec leur environnement. Un système intéressant pour décortiquer ces mécanismes de perception chimiosensorielle est le système olfactif, de la drosophile (Drosophila melanogaster), aussi appelée mouche du vinaigre. C'est un animal modèle très utile grâce à la combinaison d'outils génétiques puissants et d'analyses physiologiques facilement réalisables. Dans l'antenne de la drosophile, qui est l'organe olfactif principal de cet animal, se trouve une structure appelée sacculus. Celui-ci est composé de trois chambres distinctes, chacune comprenant plusieurs sensilles à l'intérieur desquelles se trouvent deux à trois neurones. De précédentes études morphologiques et anatomiques des ces neurones ont déterminé qu'ils sont impliqués dans la perception des odeurs, de l'humidité et de la température. Malgré ceci, la fonction principale de cet organe reste largement inconnue, principalement car il est inaccessible aux analyses électrophysiologiques. Récemment, une nouvelle famille de soixante-six récepteurs olfactifs, nommés Récepteurs lonotropiques (IRs), a été découverte chez la drosophile. Cinq IRs sont exprimés dans les neurones du sacculus. Pour la première fois, une sélection de marqueurs moléculaires est disponible pour l'étude de ces cellules. Dans cette thèse, les caractéristiques moléculaires, physiologiques et anatomiques des neurones du sacculus sont décrites. Ces populations de neurones situés dans le sacculus ont été marquées avec des gènes rapporteurs (CD8:GFP). Ceci a montré que les récepteurs IR40a et IR93a sont exprimés ensemble avec le co-récepteur IR25a dans les chambres I et II, tandis que les neurones de la chambre III expriment IR64a avec son co-récepteur IR8a. Ces deux groupes de neurones projettent vers deux glomérules distincts du lobe antennaire : les neurones IR40a projettent vers la column et le arm, alors que les neurones IR64a projettent vers DC4 et DP1m. Un screen d'imagerie optique a démontré que ces neurones sont en effet des neurones olfactifs, et que les neurones IR64a reconnaissent des ligands acides, tandis que les neurones IR40a reconnaissent des ligands aminés. De plus, les neurones IR40a et IR64a sont séparés en sous-populations distinctes anatomiquement et physiologiquement, et d'autres facteurs permettant de définir leurs propriétés fonctionnelles sont probablement impliqués. Cette thèse identifie ainsi le sacculus comme un organe olfactif spécialisé capable de détecter des acides et amines, lesquels sont très importants pour les insectes. Toutes les données collectées durant cette thèse, combinées aux données d'autres laboratoires, montrent que le sacculus est composé de différentes populations de neurones olfactifs et thermosenseurs. Ces IRs sont très conservés parmi les insectes, suggérant que le sacculus représente révolution d'un ancien détecteur olfactif d'acides et de bases chez l'insecte. - Tous les animaux sont capables de percevoir les signaux chimiques dans leur environnement, comme les odeurs ou le goût, via différents organes. L'odorat est le sens qui permet de percevoir les odeurs, et il est implique des neurones olfactifs qui se trouvent dans le nez des mammifères ou les antennes des insectes. La capacité d'un neurone olfactif à détecter une molécule odorante dépend des types de récepteurs olfactifs qu'il exprime. Il existe deux grandes familles de récepteurs qui perçoivent les odeurs : les Récepteurs Olfactifs, ORs, et Récepteurs lonotropiques IRs, qui détectent différents types d'odeurs avec différents mécanismes. Lorsqu'un récepteur reconnaît une molécule odorante, il convertit ce signal en un signal électrique qui est ensuite transmis au centre olfactif dans le cerveau. La drosophile (Drosophila melanogaster), aussi appelée mouche du vinaigre, est utilisée comme animal modèle pour étudier l'odorat, parce que son génome entier a été séquencé et que ses gènes sont facilement manipulables. De plus, l'anatomie du système olfactif de la mouche est similaire à celui des mammifères, malgré qu'il possède moins de neurones, ce qui le rend moins complexe. Ma thèse a pour objectif d'étudier les Récepteurs lonotropiques dans un organe spécifique, appelé le sacculus, situé dans les antennes. Les neurones du sacculus exprimant des IRs envoient leurs projections au centre olfactif du cerveau, suggérant que ces neurones perçoivent les odeurs. Une technique d'imagerie optique a été utilisée sur le cerveau de mouches vivantes afin de mesurer la réponse des neurones du le sacculus à différentes odeurs. J'ai démontré que ces récepteurs détectent des acides et des amines, qui sont très importants pour les insectes. Par exemple, les acides se retrouvent dans les fruits mûrs sur lesquels les mouches vont se nourrir, s'accoupler et poser leurs oeufs, et les amines sont souvent produites par des bactéries pouvant être nuisible pour la mouche. La principale découverte de ma thèse est donc l'identification du sacculus comme un organe capable de détecter deux des principales odeurs importantes pour la mouche. Ces récepteurs sont aussi présents dans d'autres insectes où ils jouent peut-être des rôles différents. Les acides et les amines se retrouvent aussi dans les excrétions (comme la sueur ou l'urine) de beaucoup de mammifères, qui pourraient potentiellement être dangereux pour la mouche, mais qui attirent les moustiques se nourrissant de leur sang.

New gene functions in megakaryopoiesis and platelet formation.

Platelets are the second most abundant cell type in blood and are essential for maintaining haemostasis. Their count and volume are tightly controlled within narrow physiological ranges, but there is only limited understanding of the molecular processes controlling both traits. Here we carried out a high-powered meta-analysis of genome-wide association studies (GWAS) in up to 66,867 individuals of European ancestry, followed by extensive biological and functional assessment. We identified 68 genomic loci reliably associated with platelet count and volume mapping to established and putative novel regulators of megakaryopoiesis and platelet formation. These genes show megakaryocyte-specific gene expression patterns and extensive network connectivity. Using gene silencing in Danio rerio and Drosophila melanogaster, we identified 11 of the genes as novel regulators of blood cell formation. Taken together, our findings advance understanding of novel gene functions controlling fate-determining events during megakaryopoiesis and platelet formation, providing a new example of successful translation of GWAS to function.

Life-history consequences of adaptation to larval nutritional stress in Drosophila.

Many animal species face periods of chronic nutritional stress during which the individuals must continue to develop, grow, and/or reproduce despite low quantity or quality of food. Here, we use experimental evolution to study adaptation to such chronic nutritional stress in six replicate Drosophila melanogaster populations selected for the ability to survive and develop within a limited time on a very poor larval food. In unselected control populations, this poor food resulted in 20% lower egg-to-adult viability, 70% longer egg-to-adult development, and 50% lower adult body weight (compared to the standard food on which the flies were normally maintained). The evolutionary changes associated with adaptation to the poor food were assayed by comparing the selected and control lines in a common environment for different traits after 29-64 generations of selection. The selected populations evolved improved egg-to-adult viability and faster development on poor food. Even though the adult dry weight of selected flies when raised on the poor food was lower than that of controls, their average larval growth rate was higher. No differences in proportional pupal lipid content were observed. When raised on the standard food, the selected flies showed the same egg-to-adult viability and the same resistance to larval heat and cold shock as the controls and a slightly shorter developmental time. However, despite only 4% shorter development time, the adults of selected populations raised on the standard food were 13% smaller and showed 20% lower early-life fecundity than the controls, with no differences in life span. The selected flies also turned out less tolerant to adult malnutrition. Thus, fruit flies have the genetic potential to adapt to poor larval food, with no detectable loss of larval performance on the standard food. However, adaptation to larval nutritional stress is associated with trade-offs with adult fitness components, including adult tolerance to nutritional stress.

FTZ-F1-related orphan receptors in Xenopus laevis: transcriptional regulators differentially expressed during early embryogenesis.

Orphan receptors of the FTZ-F1-related group of nuclear receptors (xFF1r) were identified in Xenopus laevis by isolation of cDNAs from a neurula stage library. Two cDNAs were found, which encode full length, highly related receptor proteins, xFF1rA and B, whose closet relative known so far is the murine LRH-1 orphan receptor. xFF1rA protein expressed by a recombinant vaccinia virus system specifically binds to FTZ-F1 response elements (FRE; PyCAAGGPyCPu). In cotransfection studies, xFF1rA constitutively activates transcription, in a manner dependent on the number of FREs. The amounts of at least four mRNAs encoding full-length receptors greatly increase between gastrula and early tailbud stages and decrease at later stages. At early tailbud stages, xFTZ-F1-related antigens are found in all nuclei of the embryo.

An olfactory receptor for food-derived odours promotes male courtship in Drosophila.

Many animals attract mating partners through the release of volatile sex pheromones, which can convey information on the species, gender and receptivity of the sender to induce innate courtship and mating behaviours by the receiver. Male Drosophila melanogaster fruitflies display stereotyped reproductive behaviours towards females, and these behaviours are controlled by the neural circuitry expressing male-specific isoforms of the transcription factor Fruitless (FRU(M)). However, the volatile pheromone ligands, receptors and olfactory sensory neurons (OSNs) that promote male courtship have not been identified in this important model organism. Here we describe a novel courtship function of Ionotropic receptor 84a (IR84a), which is a member of the chemosensory ionotropic glutamate receptor family, in a previously uncharacterized population of FRU(M)-positive OSNs. IR84a-expressing neurons are activated not by fly-derived chemicals but by the aromatic odours phenylacetic acid and phenylacetaldehyde, which are widely found in fruit and other plant tissues that serve as food sources and oviposition sites for drosophilid flies. Mutation of Ir84a abolishes both odour-evoked and spontaneous electrophysiological activity in these neurons and markedly reduces male courtship behaviour. Conversely, male courtship is increased--in an IR84a-dependent manner--in the presence of phenylacetic acid but not in the presence of another fruit odour that does not activate IR84a. Interneurons downstream of IR84a-expressing OSNs innervate a pheromone-processing centre in the brain. Whereas IR84a orthologues and phenylacetic-acid-responsive neurons are present in diverse drosophilid species, IR84a is absent from insects that rely on long-range sex pheromones. Our results suggest a model in which IR84a couples food presence to the activation of the fru(M) courtship circuitry in fruitflies. These findings reveal an unusual but effective evolutionary solution to coordinate feeding and oviposition site selection with reproductive behaviours through a specific sensory pathway.

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).

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 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.

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.

Assessment of transcript reconstruction methods for RNA-seq.

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.

Disruption of gene expression in hybrids of the fire ants Solenopsis invicta and Solenopsis richteri.

Transcriptome analysis is a powerful tool for unveiling the distribution and magnitude of genetic incompatibilities between hybridizing taxa. The nature of such incompatibilities is closely associated with the evolutionary histories of the parental species and may differ across tissues and between the sexes. In eusocial insects, the presence of castes that experience divergent selection regimes may result in additional distinct patterns of caste-specific hybrid incompatibilities. We analysed levels of expression of >14 000 genes in two life stages of each caste in the fire ants Solenopsis invicta and Solenopsis richteri and in their hybrids. We found strong contributions of both developmental stage and caste to gene expression patterns. In contrast, variability in expression was only weakly associated with taxonomic identity, with hybrid scores falling between those of the two parental species. Hybrid incompatibilities were surprisingly modest, with only 32 genes being mis-expressed, indicating low levels of disruption in gene regulation in hybrids; males and workers each mis-expressed at least seven times as many genes as queens. Interestingly, homologues of many of the mis-expressed genes have been implicated in behavioural variation in Drosophila melanogaster. General expression profiles of hybrids consistently were more similar to those of S. richteri than S. invicta, presumably because S. richteri trans-regulatory elements tend to be dominant and/or because there is an overall bias in the genetic composition of the hybrids towards S. richteri. Altogether, our results suggest that selection acting on each caste may contribute differently to interspecific divergence and speciation in this group of ants.

No evidence that within-group male relatedness reduces harm to females in Drosophila.

Conflict between males and females over whether, when, and how often to mate often leads to the evolution of sexually antagonistic interactions that reduce female reproductive success. Because the offspring of relatives contribute to inclusive fitness, high relatedness between rival males might be expected to reduce competition and result in the evolution of reduced harm to females. A recent study investigated this possibility in Drosophila melanogaster and concluded that groups of brothers cause less harm to females than groups of unrelated males, attributing the effect to kin selection. That study did not control for the rearing environment of males, rendering the results impossible to interpret in the context of kin selection. Here, we conducted a similar experiment while manipulating whether males developed with kin prior to being placed with females. We found no difference between related and unrelated males in the harm caused to females when males were reared separately. In contrast, when related males developed and emerged together before the experiment, female reproductive output was higher. Our results show that relatedness among males is insufficient to reduce harm to females, while a shared rearing environment - resulting in males similar to or familiar with one another - is necessary to generate this pattern.

Functional interactions between the estrogen receptor and the transcription activator Sp1 regulate the estrogen-dependent transcriptional activity of the vitellogenin A1 io promoter.

Two distinct, TATA box-containing promoters regulate the transcriptional activity of the Xenopus vitellogenin A1 gene. These two promoters are of different strength and are separated by 1.8 kilobase pairs of untranslated sequence. Estrogen receptor (ER) and its ligand, 17beta-estradiol, induce the activity of both promoters. The estrogen response elements (EREs) are located proximal to the downstream i promoter while no ERE-like sequences have been identified in the vicinity of the upstream io promoter. We show here, that transcriptional activity of the upstream io promoter is Sp1-dependent. Moreover, we demonstrate that estrogen inducibility of the io promoter results from functional interactions between the io bound Sp1 and the ER bound at the proximity of i. Functional interactions between Sp1 and ER do not require the presence of a TATA box for transcriptional activation, as is demonstrated using the acyl-CoA oxidase promoter. The relative positions that ER and Sp1 occupy with respect to the initiation site determines whether these two transcription activators can synergize for transcription initiation.