Genometrics as an essential tool for the assembly of whole genome sequences: the example of the chromosome of Bifidobacterium longum NCC2705.

BACKGROUND: Analysis of the first reported complete genome sequence of Bifidobacterium longum NCC2705, an actinobacterium colonizing the gastrointestinal tract, uncovered its proteomic relatedness to Streptomyces coelicolor and Mycobacterium tuberculosis. However, a rapid scrutiny by genometric methods revealed a genome organization totally different from all so far sequenced high-GC Gram-positive chromosomes. RESULTS: Generally, the cumulative GC- and ORF orientation skew curves of prokaryotic genomes consist of two linear segments of opposite slope: the minimum and the maximum of the curves correspond to the origin and the terminus of chromosome replication, respectively. However, analyses of the B. longum NCC2705 chromosome yielded six, instead of two, linear segments, while its dnaA locus, usually associated with the origin of replication, was not located at the minimum of the curves. Furthermore, the coorientation of gene transcription with replication was very low. Comparison with closely related actinobacteria strongly suggested that the chromosome of B. longum was misassembled, and the identification of two pairs of relatively long homologous DNA sequences offers the possibility for an alternative genome assembly proposed here below. By genometric criteria, this configuration displays all of the characters common to bacteria, in particular to related high-GC Gram-positives. In addition, it is compatible with the partially sequenced genome of DJO10A B. longum strain. Recently, a corrected sequence of B. longum NCC2705, with a configuration similar to the one proposed here below, has been deposited in GenBank, confirming our predictions. CONCLUSION: Genometric analyses, in conjunction with standard bioinformatic tools and knowledge of bacterial chromosome architecture, represent fast and straightforward methods for the evaluation of chromosome assembly.

DNA banking for plant breeding, biotechnology and biodiversity evaluation.

The manipulation of DNA is routine practice in botanical research and has made a huge impact on plant breeding, biotechnology and biodiversity evaluation. DNA is easy to extract from most plant tissues and can be stored for long periods in DNA banks. Curation methods are well developed for other botanical resources such as herbaria, seed banks and botanic gardens, but procedures for the establishment and maintenance of DNA banks have not been well documented. This paper reviews the curation of DNA banks for the characterisation and utilisation of biodiversity and provides guidelines for DNA bank management. It surveys existing DNA banks and outlines their operation. It includes a review of plant DNA collection, preservation, isolation, storage, database management and exchange procedures. We stress that DNA banks require full integration with existing collections such as botanic gardens, herbaria and seed banks, and information retrieval systems that link such facilities, bioinformatic resources and other DNA banks. They also require efficient and well-regulated sample exchange procedures. Only with appropriate curation will maximum utilisation of DNA collections be achieved.

Identification and characterization of the members of the pho1 gene family in moss and rice

Inorganic phosphate (Pi) is one of the main nutrients limiting plant growth anddevelopment in many agro-ecosystems. In plants, phosphate is acquired from the soil by theroots, and is then transferred to the shoot via the xylem. In the model plant Arabidopsisthaliana, PHO1 was previously identified as being involved in loading Pi into the xylem ofroots. AtPHO1, belongs to a multigenic family composed of 10 additional members, namelyAtPHO1;H1 to AtPHO1;10. In this study, we aimed at further investigating the role of thePHO1 gene family in Pi homeostasis in plants, and to this end we isolated and characterizedthe PHO1 members of two main model plants, the moss Physcomitrella patens and the riceOryza sativa.In the bryophyte P. patens, bioinformatic analyses revealed the presence of seven AtPHO1homologues, highly similar to AtPHO1. The seven moss PHO1 genes, namely PpPHO1;1 toPpPHO1;7 appeared to be differentially regulated, both at the tissue level and in response toPi status. However only PpPHO1;1 and PpPHO1;7 were specifically up-regulated upon Pistarvation, suggesting a potential role in Pi homeostasis. We also characterized the responseof P. patens to Pi starvation, showing that higher and lower plants share some commonstrategies to adapt to Pi-deficiency.In the second part, focusing on the monocotyledon rice, we showed the existence of threePHO1 homologues OsPHO1;1 to OsPHO1;3, with the unique particularity of each havingNatural Antisense Transcripts (NATs). Molecular analyses revealed that both the sense andthe antisense OsPHO1;2 transcripts were by far the most abundantly expressed transcripts ofthe family, preferentially expressed in the roots. The stable expression of OsPHO1;2 in allconditions tested, in opposition with the highly induced antisense transcript upon Pistarvation, suggest a putative role for the antisense in regulating the sense transcript.Moreover, mutant analyses revealed that OsPHO1;2 plays a key role in Pi homeostasis, intransferring Pi from the root to the shoot. Finally, complementing the pho1 mutant inArabidopsis, characterized by low Pi in the shoot and reduced growth, with the riceOsPHO1;2 gene revealed a new role for PHO1 in Pi signaling. Indeed, the complementedplants showed normal growth, with however low Pi content.

Secretion of an endogenous subtilisin by Pichia pastoris strains GS115 and KM71.

The methylotrophic yeast Pichia pastoris is widely used for the expression of heterologous enzymes. While the purity of the desired expression product is of major importance for many applications, we found that recombinant enzymes produced in methanol medium were contaminated by a 37-kDa endogenous yeast protease. This enzyme was completely inhibited by phenylmethanesulfonyl fluoride (PMSF) but not by 1,10-phenanthroline, EDTA, and pepstatin A, suggesting the nature of a serine protease. Its secretion was abolished in P. pastoris strains GS115 and KM71 by specific mutagenesis of a subtilisin gene (SUB2) but not by inactivation of the gene encoding vacuolar proteinase B (PRB). Bioinformatic comparisons of Sub2 protein with subtilisins from other fungal genomes and phylogenetic analyses indicated that this enzyme is not an orthologue of the vacuolar protease cerevisin generally present in yeasts but is more closely related to another putative subtilisin found in a small number of yeast genomes. During growth of P. pastoris, Sub2 was produced as a secreted enzyme at a concentration of 10 microg/ml of culture supernatant after overexpression of the full-length SUB2 gene. During fermentative production of recombinant enzymes in methanol medium, 1 ml of P. pastoris culture supernatant was found to contain approximately 3 ng of Sub2, while the enzyme was not detected during growth in a medium containing glycerol as a carbon source. The mutant strain GS115-sub2 was subsequently used as a host for the production of recombinant proteases without endogenous subtilisin contamination.

Comparative gene finding in chicken indicates that we are closing in on the set of multi-exonic widely expressed human genes.

The recent availability of the chicken genome sequence poses the question of whether there are human protein-coding genes conserved in chicken that are currently not included in the human gene catalog. Here, we show, using comparative gene finding followed by experimental verification of exon pairs by RT-PCR, that the addition to the multi-exonic subset of this catalog could be as little as 0.2%, suggesting that we may be closing in on the human gene set. Our protocol, however, has two shortcomings: (i) the bioinformatic screening of the predicted genes, applied to filter out false positives, cannot handle intronless genes; and (ii) the experimental verification could fail to identify expression at a specific developmental time. This highlights the importance of developing methods that could provide a reliable estimate of the number of these two types of genes.

Regulation of 11β-hydroxysteroid dehydrogenase type 2 by microRNA.

The enzyme 11β-hydroxysteroid dehydrogenase type 2 (11β-HSD2) is selectively expressed in aldosterone target tissues, conferring aldosterone selectivity for the mineralocorticoid receptor. A diminished activity causes salt-sensitive hypertension. The mechanism of the variable and distinct 11β-hydroxysteroid dehydrogenase type 2 gene (HSD11B2) expression in the cortical collecting duct is poorly understood. Here, we analyzed for the first time whether the 11β-HSD2 expression is modulated by microRNAs (miRNAs). In silico analysis revealed 53 and 27 miRNAs with potential binding sites on human or rat HSD11B2 3'-untranslated region. A reporter assay demonstrated 3'-untranslated region-dependent regulation of human and rodent HSD11B2. miRNAs were profiled from cortical collecting ducts and proximal convoluted tubules. Bioinformatic analyses showed a distinct clustering for cortical collecting ducts and proximal convoluted tubules with 53 of 375 miRNAs, where 13 were predicted to bind to the rat HSD11B2 3'-untranslated region. To gain insight into potentially relevant miRNAs in vivo, we investigated 2 models with differential 11β-HSD2 activity linked with salt-sensitive hypertension. (1) Comparing Sprague-Dawley with low and Wistar rats with high 11β-HSD2 activity revealed rno-miR-20a-5p, rno-miR-19b-3p, and rno-miR-190a-5p to be differentially expressed. (2) Uninephrectomy lowered 11β-HSD2 activity in the residual kidney with differentially expressed rno-miR-19b-3p, rno-miR-29b-3p, and rno-miR-26-5p. In conclusion, miRNA-dependent mechanisms seem to modulate 11β-HSD2 dosage in health and disease states.

p21 as a transcriptional co-repressor of S-phase and mitotic control genes.

It has been previously described that p21 functions not only as a CDK inhibitor but also as a transcriptional co-repressor in some systems. To investigate the roles of p21 in transcriptional control, we studied the gene expression changes in two human cell systems. Using a human leukemia cell line (K562) with inducible p21 expression and human primary keratinocytes with adenoviral-mediated p21 expression, we carried out microarray-based gene expression profiling. We found that p21 rapidly and strongly repressed the mRNA levels of a number of genes involved in cell cycle and mitosis. One of the most strongly down-regulated genes was CCNE2 (cyclin E2 gene). Mutational analysis in K562 cells showed that the N-terminal region of p21 is required for repression of gene expression of CCNE2 and other genes. Chromatin immunoprecipitation assays indicated that p21 was bound to human CCNE2 and other p21-repressed genes gene in the vicinity of the transcription start site. Moreover, p21 repressed human CCNE2 promoter-luciferase constructs in K562 cells. Bioinformatic analysis revealed that the CDE motif is present in most of the promoters of the p21-regulated genes. Altogether, the results suggest that p21 exerts a repressive effect on a relevant number of genes controlling S phase and mitosis. Thus, p21 activity as inhibitor of cell cycle progression would be mediated not only by the inhibition of CDKs but also by the transcriptional down-regulation of key genes.

Seventy-five genetic loci influencing the human red blood cell.

Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth retardation and impaired physical capacity. To understand further the genetic factors influencing red blood cells, we carried out a genome-wide association study of haemoglobin concentration and related parameters in up to 135,367 individuals. Here we identify 75 independent genetic loci associated with one or more red blood cell phenotypes at P < 10(-8), which together explain 4-9% of the phenotypic variance per trait. Using expression quantitative trait loci and bioinformatic strategies, we identify 121 candidate genes enriched in functions relevant to red blood cell biology. The candidate genes are expressed preferentially in red blood cell precursors, and 43 have haematopoietic phenotypes in Mus musculus or Drosophila melanogaster. Through open-chromatin and coding-variant analyses we identify potential causal genetic variants at 41 loci. Our findings provide extensive new insights into genetic mechanisms and biological pathways controlling red blood cell formation and function.

Plasticity of protein expression during culture of fetal skin cells.

In order to gain insight into the biology of fetal skin during culture, cellular proteins were studied during four culture passages (P00, P01, P04 as well as P10) using high-resolution two-dimensional (2-D) gel electrophoresis and mass spectrometry (MS). Bioinformatic analyses were focused on a region of each gel corresponding to pI between 4 and 8 and M(r) from 8000 to 35 000. In this area, 373 +/- 42 spots were detected (N = 18). Twenty-six spots presented an integrated intensity that increased in the higher passages, whereas five spots showed a progressively lower intensity in subsequent passaging. MS analysis was performed on spots that were unambiguously identified on preparative 2-D gels. Among the 26 spots showing an increased size between P00 and P10, 9 were identified, and corresponded to 3 proteins: (i) peptidyl-prolyl cis-trans isomerase A (P05092; cyclophilin A or cyclosporin A-binding protein), (ii) triosephosphate isomerase (P00938), and (iii) enoyl-CoA hydratase (P30084). Among these nine identified spots, three were absent at P00, but were present at P10. They corresponded to isoforms of peptidyl-prolyl cis-trans isomerase and triosephosphate isomerase, respectively. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) analyses of the acidic isoforms of triosephosphate isomerase showed modifications of cysteine residues to cysteic acid. All these isoforms were clearly present in the skin cells of a 4-year-old child, as well as in skin cells from a 80-year-old man, at P00. These observations probably reflect either an oxidative stress related to cell culture, or, alternatively, maturation, differentiation and the aging of the cells.

Functional genomics of intracellular bacteria.

During the genomic era, a large amount of whole-genome sequences accumulated, which identified many hypothetical proteins of unknown function. Rapidly, functional genomics, which is the research domain that assign a function to a given gene product, has thus been developed. Functional genomics of intracellular pathogenic bacteria exhibit specific peculiarities due to the fastidious growth of most of these intracellular micro-organisms, due to the close interaction with the host cell, due to the risk of contamination of experiments with host cell proteins and, for some strict intracellular bacteria such as Chlamydia, due to the absence of simple genetic system to manipulate the bacterial genome. To identify virulence factors of intracellular pathogenic bacteria, functional genomics often rely on bioinformatic analyses compared with model organisms such as Escherichia coli and Bacillus subtilis. The use of heterologous expression is another common approach. Given the intracellular lifestyle and the many effectors that are used by the intracellular bacteria to corrupt host cell functions, functional genomics is also often targeting the identification of new effectors such as those of the T4SS of Brucella and Legionella.

Evolution and functional characterisation of the IR chemosensory receptors in the Drosophila larval gustatory system

Chemosensory receptor gene families encode divergent proteins capable of detecting a huge diversity of environmental stimuli that are constantly changing over evolutionary time as organisms adapt to distinct ecological niches. While olfaction is dedicated to the detection of volatile compounds, taste is key to assess food quality for nutritional value and presence of toxic substances. The sense of taste also provides initial signals to mediate endocrine regulation of appetite and food metabolism and plays a role in kin recognition. The fruit fly Drosophila melanogaster is a very good model for studying smell and taste because these senses are very important in insects and because a broad variety of genetic tools are available in Drosophila. Recently, a family of 66 chemosensory receptors, the Ionotropic Receptors (IRs) was described in fruit flies. IRs are distantly related to ionotropic glutamate receptors (iGluRs), but their evolutionary origin from these synaptic receptors is unclear. While 16 IRs are expressed in the olfactory system, nothing is known about the other members of this repertoire. In this thesis, I describe bioinformatic, expression and functional analyses of the IRs aimed at understanding how these receptors have evolved, and at characterising the role of the non-olfactory IRs. I show that these have emerged at the basis of the protostome lineage and probably have acquired their sensory function very early. Moreover, although several IRs are conserved across insects, there are rapid and dramatic changes in the size and divergence of IR repertoires across species. I then performed a comprehensive analysis of IR expression in the larva of Drosophila melanogaster, which is a good model to study taste and feeding mechanisms as it spends most of its time eating or foraging. I found that most of the divergent members of the IR repertoire are expressed in both peripheral and internal gustatory neurons, suggesting that these are involved in taste perception. Finally, through the establishment of a new neurophysiological assay in larvae, I identified for the first time subsets of IR neurons that preferentially detect sugars and amino acids, indicating that IRs might be involved in sensing these compounds. Together, my results indicate that IRs are an evolutionarily dynamic and functionally versatile family of receptors. In contrast to the olfactory IRs that are well-conserved, gustatory IRs are rapidly evolving species-specific receptors that are likely to be involved in detecting a wide variety of tastants. - La plupart des animaux possèdent de grandes familles de récepteurs chimiosensoriels dont la fonction est de détecter l'immense diversité de composés chimiques présents dans l'environnement. Ces récepteurs évoluent en même temps que les organismes s'adaptent à leur écosystème. Il existe deux manières de percevoir ces signaux chimiques : l'olfaction et le goût. Alors que le système olfactif perçoit les composés volatiles, le sens du goût permet d'évaluer, par contact, la qualité de la nourriture, de détecter des substances toxiques et de réguler l'appétit et le métabolisme. L'un des organismes modèles les plus pertinents pour étudier le sens du goût est le stade larvaire de la mouche du vinaigre Drosophila melanogaster. En effet, la principale fonction du stade larvaire est de trouver de la nourriture et de manger. De plus, il est possible d'utiliser tous les outils génétiques développés chez la drosophile. Récemment, une nouvelle famille de 66 récepteurs chimiosensoriels appelés Récepteurs Ionotropiques (IRs) a été découverte chez la drosophile. Bien que leur orogine soit peu claire, ces récepteurs sont similaires aux récepteurs ionotropiques glutamatergiques impliqués dans la transmission synaptique. 16 IRs sont exprimés dans le système olfactif de la mouche adulte, mais pour l'instant on ne connaît rien des autres membres de cette famille. Durant ma thèse, j'ai effectué des recherches sur l'évolution de ces récepteurs ainsi que sur l'expression et la fonction des IRs non olfactifs. Je démontre que les IRs sont apparus chez l'ancêtre commun des protostomiens et ont probablement acquis leur fonction sensorielle très rapidement. De plus, bien qu'un certain nombre d'IRs olfactifs soient conservés chez les insectes, d'importantes variations dans la taille et la divergence des répertoires d'IRs entre les espèces ont été constatées. J'ai également découvert qu'un grand nombre d'IRs non olfactifs sont exprimés dans différents organes gustatifs, ce qui leur confère probablement une fonction dans la perception des goûts. Finalement, pour la première fois, des neurones exprimant des IRs ont été identifiés pour leur fonction dans la perception de sucres et d'acides aminés chez la larve. Mes résultats présentent les IRs comme une famille très dynamique, aux fonctions très variées, qui joue un rôle tant dans l'odorat que dans le goût, et dont la fonction est restée importante tout au long de l'évolution. De plus, l'identification de neurones spécialisés dans la perception de certains composés permettra l'étude des circuits neuronaux impliqués dans le traitement de ces informations.

Computational analysis of the genetic and environmental contributions to disease-related human phenotypes: From predicting adverse lipid response of HIV patients receiving antiretroviral therapy to genome-wide association studies of cardiovascular and lipid related disorders

Genetic variants influence the risk to develop certain diseases or give rise to differences in drug response. Recent progresses in cost-effective, high-throughput genome-wide techniques, such as microarrays measuring Single Nucleotide Polymorphisms (SNPs), have facilitated genotyping of large clinical and population cohorts. Combining the massive genotypic data with measurements of phenotypic traits allows for the determination of genetic differences that explain, at least in part, the phenotypic variations within a population. So far, models combining the most significant variants can only explain a small fraction of the variance, indicating the limitations of current models. In particular, researchers have only begun to address the possibility of interactions between genotypes and the environment. Elucidating the contributions of such interactions is a difficult task because of the large number of genetic as well as possible environmental factors.In this thesis, I worked on several projects within this context. My first and main project was the identification of possible SNP-environment interactions, where the phenotypes were serum lipid levels of patients from the Swiss HIV Cohort Study (SHCS) treated with antiretroviral therapy. Here the genotypes consisted of a limited set of SNPs in candidate genes relevant for lipid transport and metabolism. The environmental variables were the specific combinations of drugs given to each patient over the treatment period. My work explored bioinformatic and statistical approaches to relate patients' lipid responses to these SNPs, drugs and, importantly, their interactions. The goal of this project was to improve our understanding and to explore the possibility of predicting dyslipidemia, a well-known adverse drug reaction of antiretroviral therapy. Specifically, I quantified how much of the variance in lipid profiles could be explained by the host genetic variants, the administered drugs and SNP-drug interactions and assessed the predictive power of these features on lipid responses. Using cross-validation stratified by patients, we could not validate our hypothesis that models that select a subset of SNP-drug interactions in a principled way have better predictive power than the control models using "random" subsets. Nevertheless, all models tested containing SNP and/or drug terms, exhibited significant predictive power (as compared to a random predictor) and explained a sizable proportion of variance, in the patient stratified cross-validation context. Importantly, the model containing stepwise selected SNP terms showed higher capacity to predict triglyceride levels than a model containing randomly selected SNPs. Dyslipidemia is a complex trait for which many factors remain to be discovered, thus missing from the data, and possibly explaining the limitations of our analysis. In particular, the interactions of drugs with SNPs selected from the set of candidate genes likely have small effect sizes which we were unable to detect in a sample of the present size (<800 patients).In the second part of my thesis, I performed genome-wide association studies within the Cohorte Lausannoise (CoLaus). I have been involved in several international projects to identify SNPs that are associated with various traits, such as serum calcium, body mass index, two-hour glucose levels, as well as metabolic syndrome and its components. These phenotypes are all related to major human health issues, such as cardiovascular disease. I applied statistical methods to detect new variants associated with these phenotypes, contributing to the identification of new genetic loci that may lead to new insights into the genetic basis of these traits. This kind of research will lead to a better understanding of the mechanisms underlying these pathologies, a better evaluation of disease risk, the identification of new therapeutic leads and may ultimately lead to the realization of "personalized" medicine.

In vitro characterization of PlySK1249, a novel phage lysin, and assessment of its antibacterial activity in a mouse model of Streptococcus agalactiae bacteremia.

Beta-hemolytic Streptococcus agalactiae is the leading cause of bacteremia and invasive infections. These diseases are treated with β-lactams or macrolides, but the emergence of less susceptible and even fully resistant strains is a cause for concern. New bacteriophage lysins could be promising alternatives against such organisms. They hydrolyze the bacterial peptidoglycan at the end of the phage cycle, in order to release the phage progeny. By using a bioinformatic approach to screen several beta-hemolytic streptococci, a gene coding for a lysin was identified on a prophage carried by Streptococcus dysgalactiae subsp. equisimilis SK1249. The gene product, named PlySK1249, harbored an original three-domain structure with a central cell wall-binding domain surrounded by an N-terminal amidase and a C-terminal CHAP domain. Purified PlySK1249 was highly lytic and bactericidal for S. dysgalactiae (2-log10 CFU/ml decrease within 15 min). Moreover, it also efficiently killed S. agalactiae (1.5-log10 CFU/ml decrease within 15 min) but not several streptococcal commensal species. We further investigated the activity of PlySK1249 in a mouse model of S. agalactiae bacteremia. Eighty percent of the animals (n = 10) challenged intraperitoneally with 10(6) CFU of S. agalactiae died within 72 h, whereas repeated injections of PlySK1249 (45 mg/kg 3 times within 24 h) significantly protected the mice (P < 0.01). Thus, PlySK1249, which was isolated from S. dysgalactiae, demonstrated high cross-lytic activity against S. agalactiae both in vitro and in vivo. These encouraging results indicated that PlySK1249 might represent a good candidate to be developed as a new enzybiotic for the treatment of systemic S. agalactiae infections.

An operon of three transcriptional regulators controls horizontal gene transfer of the integrative and conjugative element ICEclc in Pseudomonas knackmussii B13.

The integrative and conjugative element ICEclc is a mobile genetic element in Pseudomonas knackmussii B13, and an experimental model for a widely distributed group of elements in Proteobacteria. ICEclc is transferred from specialized transfer competent cells, which arise at a frequency of 3-5% in a population at stationary phase. Very little is known about the different factors that control the transfer frequency of this ICE family. Here we report the discovery of a three-gene operon encoded by ICEclc, which exerts global control on transfer initiation. The operon consists of three consecutive regulatory genes, encoding a TetR-type repressor MfsR, a MarR-type regulator and a LysR-type activator TciR. We show that MfsR autoregulates expression of the operon, whereas TciR is a global activator of ICEclc gene expression, but no clear role was yet found for MarR. Deletion of mfsR increases expression of tciR and marR, causing the proportion of transfer competent cells to reach almost 100% and transfer frequencies to approach 1 per donor. mfsR deletion also caused a two orders of magnitude loss in population viability, individual cell growth arrest and loss of ICEclc. This indicates that autoregulation is an important feature maintaining ICE transfer but avoiding fitness loss. Bioinformatic analysis showed that the mfsR-marR-tciR operon is unique for ICEclc and a few highly related ICE, whereas tciR orthologues occur more widely in a large variety of suspected ICE among Proteobacteria.

Ionotropic receptors (IRs): chemosensory ionotropic glutamate receptors in Drosophila and beyond.

Ionotropic Receptors (IRs) are a recently characterized family of olfactory receptors in the fruit fly, Drosophila melanogaster. IRs are not related to insect Odorant Receptors (ORs), but rather have evolved from ionotropic glutamate receptors (iGluRs), a conserved family of synaptic ligand-gated ion channels. Here, we review the expression and function of IRs in Drosophila, highlighting similarities and differences with iGluRs. We also briefly describe the organization of the neuronal circuits in which IRs function, comparing and contrasting them with the sensory pathways expressing ORs. Finally, we summarize the bioinformatic identification and initial characterization of IRs in other species, which imply an evolutionarily conserved role for these receptors in chemosensation in insects and other protostomes.