Uncovering the functional constraints underlying the genomic organization of the odorant-binding protein genes

GALEON: a comprehensive bioinformatic tool to analyse and visualize gene clusters in complete genomes

MOTIVATION

Gene clusters, defined as a set of genes encoding functionally related proteins, are abundant in eukaryotic genomes. Despite the increasing availability of chromosome-level genomes, the comprehensive analysis of gene family evolution remains largely unexplored, particularly for large and highly dynamic gene families or those including very recent family members. These challenges stem from limitations in genome assembly contiguity, particularly in repetitive regions such as large gene clusters. Recent advancements in sequencing technology, such as long reads and chromatin contact mapping, hold promise in addressing these challenges.

RESULTS

To facilitate the identification, analysis, and visualization of physically clustered gene family members within chromosome-level genomes, we introduce GALEON, a user-friendly bioinformatic tool. GALEON identifies gene clusters by studying the spatial distribution of pairwise physical distances among gene family members along with the genome-wide gene density. The pipeline also enables the simultaneous analysis and comparison of two gene families and allows the exploration of the relationship between physical and evolutionary distances. This tool offers a novel approach for studying the origin and evolution of gene families.

AVAILABILITY AND IMPLEMENTATION

GALEON is freely available from https://www.ub.edu/softevol/galeon and https://github.com/molevol-ub/galeon.

Evolution of the odorant-binding protein gene family in Drosophila

Odorant-binding proteins (OBPs) are encoded by a gene family involved in the perception of olfactory signals in insects. This chemosensory gene family has been advocated as a candidate to mediate host preference and host shifts in insects, although it also participates in other physiological processes. Remarkable differences in the OBP gene repertoire have been described across insect groups, suggesting an accelerated gene turnover rate. The genus Drosophila, is a valuable resource for ecological genomics studies since it comprises groups of ecologically diverse species and there are genome data for many of them. Here, we investigate the molecular evolution of this chemosensory gene family across 19 Drosophila genomes, including the melanogaster and repleta species groups, which are mostly associated with rotting fruit and cacti, respectively. We also compared the OBP repertoire among the closely related species of the repleta group, associated with different subfamilies of Cactaceae that represent disparate chemical challenges for the flies. We found that the gene family size varies widely between species, ranging from 39 to 54 candidate OBPs. Indeed, more than 54% of these genes are organized in clusters and located on chromosomes X, 2, and 5, with a distribution conserved throughout the genus. The family sizes in the repleta group and D. virilis (virilis-repleta radiation) were smaller than in the melanogaster group. We tested alternative evolutionary models for OBP family size and turnover rates based on different ecological scenarios. We found heterogeneous gene turnover rates (GR) in comparisons involving columnar cactus specialists, prickly pear specialists, and fruit dwellers lineages, and signals of rapid molecular evolution compatible with positive selection in specific OBP genes. Taking ours and previous results together, we propose that this chemosensory gene family is involved in host adaptation and hypothesize that the adoption of the cactophilic lifestyle in the repleta group accelerated the evolution of members of the family.

Reconstructing Gene Gains and Losses with BadiRate

Estimating gene gain and losses is paramount to understand the molecular mechanisms underlying adaptive evolution. Despite the advent of high-throughput sequencing, such analyses have been so far hampered by the poor contiguity of genome assemblies. The increasing affordability of long-read sequencing technologies will however revolutionize our capacity to identify gene gains and losses at an unprecedented resolution, even in non-model organisms. To thoroughly exploit all such multigene family variation, the software BadiRate implements a collection of birth-and-death stochastic models, aiming at estimating by maximum likelihood the gene turnover rates along the internal and external branches of a given phylogenetic species tree. Its statistical framework also provides versatility for inferring the gene family content at the internal phylogenetic nodes (and to estimate the minimum number of gene gains and losses in each branch), for statistically contrasting competing hypotheses (e.g., accelerations of the gene turnover rates at pre-defined clades), and for pinpointing gene family expansions or contractions likely driven by natural selection. In this chapter we review the theoretical models implemented in BadiRate and illustrate their applicability by analyzing a hypothetical data set of 14 microbial species.

Lipocalins in Arthropod Chemical Communication

Lipocalins represent one of the most successful superfamilies of proteins. Most of them are extracellular carriers for hydrophobic ligands across aqueous media, but other functions have been reported. They are present in most living organisms including bacteria. In animals they have been identified in mammals, molluscs, and arthropods; sequences have also been reported for plants. A subgroup of lipocalins, referred to as odorant-binding proteins (OBPs), mediate chemical communication in mammals by ferrying specific pheromones to the vomeronasal organ. So far, these proteins have not been reported as carriers of semiochemicals in other living organisms; instead chemical communication in arthropods is mediated by other protein families structurally unrelated to lipocalins. A search in the databases has revealed extensive duplication and differentiation of lipocalin genes in some species of insects, crustaceans, and chelicerates. Their large numbers, ranging from a handful to few dozens in the same species, their wide divergence, both within and between species, and their expression in chemosensory organs suggest that such expansion may have occurred under environmental pressure, thus supporting the hypothesis that lipocalins may be involved in chemical communication in arthropods.

Weak Polygenic Selection Drives the Rapid Adaptation of the Chemosensory System: Lessons from the Upstream Regions of the Major Gene Families

The animal chemosensory system is involved in essential biological processes, most of them mediated by proteins encoded in multigene families. These multigene families have been fundamental for the adaptation to new environments, significantly contributing to phenotypic variation. This adaptive potential contrasts, however, with the lack of studies at their upstream regions, especially taking into account the evidence linking their transcriptional changes to certain phenotypic effects. Here, we explicitly characterize the contribution of the upstream sequences of the major chemosensory gene families to rapid adaptive processes. For that, we analyze the genome sequences of 158 lines from a population of Drosophila melanogaster that recently colonized North America, and integrate functional and transcriptional data available for this species. We find that both, strong negative and strong positive selection, shape transcriptional evolution at the genome-wide level. The chemosensory upstream regions, however, exhibit a distinctive adaptive landscape, including multiple mutations of small beneficial effect and a reduced number of cis-regulatory elements. Together, our results suggest that the promiscuous and partially redundant transcription and function of the chemosensory genes provide evolutionarily opportunities for rapid adaptive episodes through weak polygenic selection.

Tissue-specific transcriptomics, chromosomal localization, and phylogeny of chemosensory and odorant binding proteins from the red flour beetle Tribolium castaneum reveal subgroup specificities for olfaction or more general functions

Background

Chemoreception is based on the senses of smell and taste that are crucial for animals to find new food sources, shelter, and mates. The initial step in olfaction involves the translocation of odorants from the periphery through the aqueous lymph of the olfactory sensilla to the odorant receptors most likely by chemosensory proteins (CSPs) or odorant binding proteins (OBPs).

Results

To better understand the roles of CSPs and OBPs in a coleopteran pest species, the red flour beetle Tribolium castaneum (Coleoptera, Tenebrionidae), we performed transcriptome analyses of male and female antennae, heads, mouthparts, legs, and bodies, which revealed that all 20 CSPs and 49 of the 50 previously annotated OBPs are transcribed. Only six of the 20 CSP are significantly transcriptionally enriched in the main chemosensory tissues (antenna and/or mouthparts), whereas of the OBPs all eight members of the antenna binding proteins II (ABPII) subgroup, 18 of the 20 classic OBP subgroup, the C + OBP, and only five of the 21 C-OBPs show increased chemosensory tissue expression. By MALDI-TOF-TOF MS protein fingerprinting, we confirmed three CSPs, four ABPIIs, three classic OBPs, and four C-OBPs in the antennae.

Conclusions

Most of the classic OBPs and all ABPIIs are likely involved in chemoreception. A few are also present in other tissues such as odoriferous glands and testes and may be involved in release or transfer of chemical signals. The majority of the CSPs as well as the C-OBPs are not enriched in antennae or mouthparts, suggesting a more general role in the transport of hydrophobic molecules.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1141) contains supplementary material, which is available to authorized users.