Synchronized Expansion and Contraction of Olfactory, Vomeronasal, and Taste Receptor Gene Families in Hystricomorph Rodents.

Chemical senses, including olfaction, pheromones, and taste, are crucial for the survival of most animals. There has long been a debate about whether different types of senses might influence each other. For instance, primates with a strong sense of vision are thought to have weakened olfactory abilities, although the oversimplified trade-off theory is now being questioned. It is uncertain whether such interactions between different chemical senses occur during evolution. To address this question, we examined four receptor gene families related to olfaction, pheromones, and taste: olfactory receptor (OR), vomeronasal receptor type 1 and type 2 (V1R and V2R), and bitter taste receptor (T2R) genes in Hystricomorpha, which is morphologically and ecologically the most diverse group of rodents. We also sequenced and assembled the genome of the grasscutter, Thryonomys swinderianus. By examining 16 available genome assemblies alongside the grasscutter genome, we identified orthologous gene groups among hystricomorph rodents for these gene families to separate the gene gain and loss events in each phylogenetic branch of the Hystricomorpha evolutionary tree. Our analysis revealed that the expansion or contraction of the four gene families occurred synchronously, indicating that when one chemical sense develops or deteriorates, the others follow suit. The results also showed that V1R/V2R genes underwent the fastest evolution, followed by OR genes, and T2R genes were the most evolutionarily stable. This variation likely reflects the difference in ligands of V1R/V2Rs, ORs, and T2Rs: species-specific pheromones, environment-based scents, and toxic substances common to many animals, respectively.

Genetic variation in the human olfactory receptor OR5AN1 associates with the perception of musks.

Humans have significant individual variations in odor perception, derived from their experience or sometimes from differences in the olfactory receptor (OR) gene repertoire. In several cases, the genetic variation of a single OR affects the perception of its cognate odor ligand. Musks are widely used for fragrance and are known to demonstrate specific anosmia. It, however, remains to be elucidated whether the OR polymorphism contributes to individual variations in musk odor perception. Previous studies reported that responses of the human musk receptor OR5AN1 to a variety of musks in vitro correlated well with perceptual sensitivity to those odors in humans and that the mouse ortholog, Olfr1440 (MOR215-1), plays a critical role in muscone perception. Here, we took advantage of genetic variation in OR5AN1 to examine how changes in receptor sensitivity are associated with human musk perception. We investigated the functional differences between OR5AN1 variants in an in vitro assay and measured both perceived intensity and detection threshold in human subjects with different OR5AN1 genotypes. Human subjects homozygous for the more sensitive L289F allele had a lower detection threshold for muscone and found macrocyclic musks to be more intense than subjects homozygous for the reference allele. These results demonstrate that the genetic variation in OR5AN1 contributes to perceptual differences for some musks. In addition, we found that the more functional variant of OR5A1, a receptor involved in ß-ionone perception, is associated with the less functional variant of OR5AN1, suggesting that the perceived intensities of macrocyclic musks and ß-ionone are inversely correlated.

Structural model for ligand binding and channel opening of an insect gustatory receptor.

Insect gustatory receptors play roles in sensing tastants, such as sugars and bitter substances. We previously demonstrated that the BmGr9 silkworm gustatory receptor is a d-fructose-gated ion channel receptor. However, the molecular mechanism of how d-fructose could initiate channel opening were unclear. Herein, we present a structural model for a channel pore and a d-fructose-binding site in BmGr9. Since the membrane topology and oligomeric state of BmGr9 appeared to be similar to those of an insect odorant receptor coreceptor, Orco, we constructed a structural model of BmGr9 based on the cryo-EM Orco structure. Our site-directed mutagenesis data suggested that the transmembrane region 7 forms channel pore and controls channel gating. This model also suggested that a pocket formed by transmembrane helices 2 to 4 and 6 binds d-fructose. Using mutagenesis experiments in combination with docking simulations, we were able to determine the potent binding mode of d-fructose. Finally, based on these data, we propose a conformational change that leads to channel opening upon d-fructose binding. Taken together, these findings detail the molecular mechanism by which an insect gustatory receptor can be activated by its ligand molecule.

The genomics of ecological flexibility, large brains, and long lives in capuchin monkeys revealed with fecalFACS.

Ecological flexibility, extended lifespans, and large brains have long intrigued evolutionary biologists, and comparative genomics offers an efficient and effective tool for generating new insights into the evolution of such traits. Studies of capuchin monkeys are particularly well situated to shed light on the selective pressures and genetic underpinnings of local adaptation to diverse habitats, longevity, and brain development. Distributed widely across Central and South America, they are inventive and extractive foragers, known for their sensorimotor intelligence. Capuchins have among the largest relative brain size of any monkey and a lifespan that exceeds 50 y, despite their small (3 to 5 kg) body size. We assemble and annotate a de novo reference genome for Cebus imitator Through high-depth sequencing of DNA derived from blood, various tissues, and feces via fluorescence-activated cell sorting (fecalFACS) to isolate monkey epithelial cells, we compared genomes of capuchin populations from tropical dry forests and lowland rainforests and identified population divergence in genes involved in water balance, kidney function, and metabolism. Through a comparative genomics approach spanning a wide diversity of mammals, we identified genes under positive selection associated with longevity and brain development. Additionally, we provide a technological advancement in the use of noninvasive genomics for studies of free-ranging mammals. Our intra- and interspecific comparative study of capuchin genomics provides insights into processes underlying local adaptation to diverse and physiologically challenging environments, as well as the molecular basis of brain evolution and longevity.

Origin and Evolution of the Gene Family of Proteinaceous Pheromones, the Exocrine Gland-Secreting Peptides, in Rodents.

The exocrine-gland secreting peptide (ESP)gene family encodes proteinaceous pheromones that are recognized by the vomeronasal organ in mice. For example, ESP1 is a male pheromone secreted in tear fluid that regulates socio-sexual behavior, and ESP22 is a juvenile pheromone that suppresses adult sexual behavior. The family consists of multiple genes and has been identified only in mouse and rat genomes. The coding region of a mouse ESP gene is separated into two exons, each encoding signal and mature sequences. Here, we report the origin and evolution of the ESP gene family. ESP genes were found only in the Muridea and Cricetidae families of rodents, suggesting a recent origin of ESP genes in the common ancestor of murids and cricetids. ESP genes show a great diversity in number, length, and sequence among different species as well as mouse strains. Some ESPs in rats and golden hamsters are expressed in the lacrimal gland and the salivary gland. We also found that a mature sequence of an ESP gene showed overall sequence similarity to the α-globin gene. The ancestral ESP gene seems to be generated by recombination of a retrotransposed α-globin gene with the signal-encoding exon of the CRISP2 gene located adjacent to the ESP gene cluster. This study provides an intriguing example of molecular tinkering in rapidly evolving species-specific proteinaceous pheromone genes.

Comparative genomic analyses illuminate the distinct evolution of megabats within Chiroptera.

The revision of the sub-order Microchiroptera is one of the most intriguing outcomes in recent mammalian molecular phylogeny. The unexpected sister-taxon relationship between rhinolophoid microbats and megabats, with the exclusion of other microbats, suggests that megabats arose in a relatively short period of time from a microbat-like ancestor. In order to understand the genetic mechanism underlying adaptive evolution in megabats, we determined the whole-genome sequences of two rousette megabats, Leschenault's rousette (Rousettus leschenaultia) and the Egyptian fruit bat (R. aegyptiacus). The sequences were compared with those of 22 other mammals, including nine bats, available in the database. We identified that megabat genomes are distinct in that they have extremely low activity of SINE retrotranspositions, expansion of two chemosensory gene families, including the trace amine receptor (TAAR) and olfactory receptor (OR), and elevation of the dN/dS ratio in genes for immunity and protein catabolism. The adaptive signatures discovered in the genomes of megabats may provide crucial insight into their distinct evolution, including key processes such as virus resistance, loss of echolocation, and frugivorous feeding.

Metabolism of Odorant Molecules in Human Nasal/Oral Cavity Affects the Odorant Perception.

In this study, we examined the mode of metabolism of food odorant molecules in the human nasal/oral cavity in vitro and in vivo. We selected 4 odorants, 2-furfurylthiol (2-FT), hexanal, benzyl acetate, and methyl raspberry ketone, which are potentially important for designing food flavors. In vitro metabolic assays of odorants with saliva/nasal mucus analyzed by gas chromatography mass spectrometry revealed that human saliva and nasal mucus exhibit the following 3 enzymatic activities: (i) methylation of 2-FT into furfuryl methylsulfide (FMS); (ii) reduction of hexanal into hexanol; and (iii) hydrolysis of benzyl acetate into benzyl alcohol. However, (iv) demethylation of methyl raspberry ketone was not observed. Real-time in vivo analysis using proton transfer reaction-mass spectrometry demonstrated that the application of 2-FT and hexanal through 3 different pathways via the nostril or through the mouth generated the metabolites FMS and hexanol within a few seconds. The concentration of FMS and hexanol in the exhaled air was above the perception threshold. A cross-adaptation study based on the activation pattern of human odorant receptors suggested that this metabolism affects odor perception. These results suggest that some odorants in food are metabolized in the human nasal mucus/saliva, and the resulting metabolites are perceived as part of the odor quality of the substrates. Our results help improve the understanding of the mechanism of food odor perception and may enable improved design and development of foods in relation to odor.

Acceleration of Olfactory Receptor Gene Loss in Primate Evolution: Possible Link to Anatomical Change in Sensory Systems and Dietary Transition.

Primates have traditionally been regarded as vision-oriented animals with low olfactory ability, though this "microsmatic primates" view has been challenged recently. To clarify when and how degeneration of the olfactory system occurred and to specify the relevant factors during primate evolution, we here examined the olfactory receptor (OR) genes from 24 phylogenetically and ecologically diverse primate species. The results revealed that strepsirrhines with curved noses had functional OR gene repertoires that were nearly twice as large as those for haplorhines with simple noses. Neither activity pattern (nocturnal/diurnal) nor color vision system showed significant correlation with the number of functional OR genes while phylogeny and nose structure (haplorhine/strepsirrhine) are statistically controlled, but extent of folivory did. We traced the evolutionary fates of individual OR genes by identifying orthologous gene groups, demonstrating that the rates of OR gene losses were accelerated at the ancestral branch of haplorhines, which coincided with the acquisition of acute vision. The highest rate of OR gene loss was observed at the ancestral branch of leaf-eating colobines; this reduction is possibly linked with the dietary transition from frugivory to folivory because odor information is essential for fruit foraging but less so for leaf foraging. Intriguingly, we found accelerations of OR gene losses in an external branch to every hominoid species examined. These findings suggest that the current OR gene repertoire in each species has been shaped by a complex interplay of phylogeny, anatomy, and habitat; therefore, multiple factors may contribute to the olfactory degeneration in primates.

A long-range cis-regulatory element for class I odorant receptor genes.

Individual olfactory sensory neurons express a single odorant receptor gene from either class I genes residing in a single cluster on a single chromosome or class II genes spread over multiple clusters on multiple chromosomes. Here, we identify an enhancer element for mouse class I genes, the J element, that is conserved through mammalian species from the platypus to humans. The J element regulates most class I genes expression by exerting an effect over ~ 3 megabases within the whole cluster. Deletion of the trans J element increases the expression frequencies of class I genes from the intact J allele, indicating that the allelic exclusion of class I genes depends on the activity of the J element. Our data reveal a long-range cis-regulatory element that governs the singular class I gene expression and has been phylogenetically preserved to retain a single cluster organization of class I genes in mammals."Each olfactory sensory neuron expresses a single odorant receptor gene from either class I or class II genes. Here, the authors identify an enhancer for mouse class I genes, that is highly conserved, and regulates most class I genes expression by acting over ~ 3 megabases within the whole cluster."

Comparative genomic analysis of translation initiation mechanisms for genes lacking the Shine-Dalgarno sequence in prokaryotes.

In prokaryotes, translation initiation is believed to occur through an interaction between the 3΄ tail of a 16S rRNA and a corresponding Shine-Dalgarno (SD) sequence in the 5΄ untranslated region (UTR) of an mRNA. However, some genes lack SD sequences (non-SD genes), and the fraction of non-SD genes in a genome varies depending on the prokaryotic species. To elucidate non-SD translation initiation mechanisms in prokaryotes from an evolutionary perspective, we statistically examined the nucleotide frequencies around the initiation codons in non-SD genes from 260 prokaryotes (235 bacteria and 25 archaea). We identified distinct nucleotide frequency biases upstream of the initiation codon in bacteria and archaea, likely because of the presence of leaderless mRNAs lacking a 5΄ UTR. Moreover, we observed overall similarities in the nucleotide patterns between upstream and downstream regions of the initiation codon in all examined phyla. Symmetric nucleotide frequency biases might facilitate translation initiation by preventing the formation of secondary structures around the initiation codon. These features are more prominent in species' genomes that harbor large fractions of non-SD sequences, suggesting that a reduced stability around the initiation codon is important for efficient translation initiation in prokaryotes.

Ligand Specificity and Evolution of Mammalian Musk Odor Receptors: Effect of Single Receptor Deletion on Odor Detection.

Musk odors have been used widely for fragrance and medicine for >2000 years because of their fascinating scent and physiological effects. Therefore, fragrance manufacturers have been eager to develop high-quality musk compounds that are safe and easily synthesized. We recently identified muscone-responsive olfactory receptors (ORs) MOR215-1 and OR5AN1 in mice and humans, respectively (Shirasu et al., 2014). In this study, we identified musk ORs that are evolutionarily closely related to MOR215-1 or OR5AN1 in various primates and investigated structure-activity relationships for various musk odorants and related compounds. We found that each species has one or two functional musk ORs that exhibit specific ligand spectra to musk compounds. Some of them, including the human OR5AN1, responded to nitro musks with chemical properties distinct from muscone. The ligand specificity of OR5AN1 reflects the perception of musk odors in humans. Genetic deletion of MOR215-1 in mice resulted in drastic reduction of sensitivity to muscone, suggesting that MOR215-1 plays a critical role in muscone perception. Therefore, the current study reveals a clear link between the identified OR and muscone perception. Moreover, the strategy established for screening ligands for the muscone OR may facilitate the development of novel and commercially useful musk odors. SIGNIFICANCE STATEMENT: The long-sought musk odor receptor family in mammals was discovered and found to be well conserved and narrowly tuned to musk odors. In mice, deletion of the most sensitive musk receptor resulted in drastic reduction in sensitivity to muscone, demonstrating a strong link between receptor and odor perception. In humans, we found one musk receptor that recognized both macrocyclic and nitro musks that had distinct chemical structures. The structure-activity relationships were in a good agreement with human sensory perception and therefore may be used to develop novel musk aroma in fragrance fields. Finally, identification of a natural ligand(s) for musk receptors in mammals other than musk deer would reveal an evolutionarily pivotal role in each species in the future.

Extreme expansion of the olfactory receptor gene repertoire in African elephants and evolutionary dynamics of orthologous gene groups in 13 placental mammals.

Olfactory receptors (ORs) detect odors in the environment, and OR genes constitute the largest multigene family in mammals. Numbers of OR genes vary greatly among species--reflecting the respective species' lifestyles--and this variation is caused by frequent gene gains and losses during evolution. However, whether the extent of gene gains/losses varies among individual gene lineages and what might generate such variation is unknown. To answer these questions, we used a newly developed phylogeny-based method to classify >10,000 intact OR genes from 13 placental mammal species into 781 orthologous gene groups (OGGs); we then compared the OGGs. Interestingly, African elephants had a surprisingly large repertoire (∼ 2000) of functional OR genes encoded in enlarged gene clusters. Additionally, OR gene lineages that experienced more gene duplication had weaker purifying selection, and Class II OR genes have evolved more dynamically than those in Class I. Some OGGs were highly expanded in a lineage-specific manner, while only three OGGs showed complete one-to-one orthology among the 13 species without any gene gains/losses. These three OGGs also exhibited highly conserved amino acid sequences; therefore, ORs in these OGGs may have physiologically important functions common to every placental mammal. This study provides a basis for inferring OR functions from evolutionary trajectory.

Identification of chemosensory receptor genes from vertebrate genomes.

Chemical senses are essential for the survival of animals. In vertebrates, mainly three different types of receptors, olfactory receptors (ORs), vomeronasal receptors type 1 (V1Rs), and vomeronasal receptors type 2 (V2Rs), are responsible for the detection of chemicals in the environment. Mouse or rat genomes contain >1,000 OR genes, forming the largest multigene family in vertebrates, and have >100 V1R and V2R genes as well. Recent advancement in genome sequencing enabled us to computationally identify nearly complete repertories of OR, V1R, and V2R genes from various organisms, revealing that the numbers of these genes are highly variable among different organisms depending on each species' living environment. Here I would explain bioinformatic methods to identify the entire repertoires of OR, V1R, and V2R genes from vertebrate genome sequences.

The draft genomes of soft-shell turtle and green sea turtle yield insights into the development and evolution of the turtle-specific body plan.

The unique anatomical features of turtles have raised unanswered questions about the origin of their unique body plan. We generated and analyzed draft genomes of the soft-shell turtle (Pelodiscus sinensis) and the green sea turtle (Chelonia mydas); our results indicated the close relationship of the turtles to the bird-crocodilian lineage, from which they split ∼267.9-248.3 million years ago (Upper Permian to Triassic). We also found extensive expansion of olfactory receptor genes in these turtles. Embryonic gene expression analysis identified an hourglass-like divergence of turtle and chicken embryogenesis, with maximal conservation around the vertebrate phylotypic period, rather than at later stages that show the amniote-common pattern. Wnt5a expression was found in the growth zone of the dorsal shell, supporting the possible co-option of limb-associated Wnt signaling in the acquisition of this turtle-specific novelty. Our results suggest that turtle evolution was accompanied by an unexpectedly conservative vertebrate phylotypic period, followed by turtle-specific repatterning of development to yield the novel structure of the shell.

Identification of olfactory receptor genes from mammalian genome sequences.

Olfaction is essential for the survival of mammals. Diverse odorant molecules in the environment are detected by olfactory receptors (ORs) expressed in the olfactory epithelium of the nasal cavity. In general, mammalian genomes harbor ~1,000 OR genes, which form the largest multigene family in mammals. The recent advances in genome sequencing technology have enabled us to computationally identify nearly complete repertoires of OR genes from various organisms. Such studies have revealed that the numbers of OR genes are highly variable among organisms depending on their living environments. Because OR genes are intronless, it is possible to find all OR genes by conducting homology searches against the genome sequences using known OR genes as queries. However, some caution is necessary during the process of extracting intact coding sequences of OR genes and distinguishing among OR and non-OR genes. Presented here is a description of bioinformatics methods to identify the entire OR gene repertoires from mammalian genome sequences.

A methodical microarray design enables surveying of expression of a broader range of genes in Ciona intestinalis.

We provide a new oligo-microarray for Ciona intestinalis, based on the NimbleGen 12-plex×135k format. The array represents 106,285 probes, which is more than double the probe number of the currently available 44k microarray. These probes cover 99.2% of the transcripts in the KyotoHoya (KH) models, published in 2008, and they contain 81.1% of the entries in the UniGene database that are not included in the KH models. In this paper, we show that gene expression levels measured by this new 135k microarray are highly correlated with those obtained by the existing 44k microarray for genes common to both arrays. We also investigated gene expression using samples obtained from the ovary and the neural complex of adult C. intestinalis, showing that the expression of tissue-specific genes is consistent with previous reports. Approximately half of the highly expressed genes identified in the 135k microarray are not included in the previous microarray. The high coverage of gene models by this microarray made it possible to identify splicing variants for a given transcript. The 135k microarray is useful in investigating the functions of genes that are not yet well characterized. Detailed information about this 135k microarray is accessible at no charge from supplemental materials, NCBI Gene Expression Omnibus (GEO), and http://marinegenomics.oist.jp.

Olfactory receptor multigene family in vertebrates: from the viewpoint of evolutionary genomics.

Olfaction is essential for the survival of animals. Diverse odor molecules in the environment are detected by the olfactory receptors (ORs) in the olfactory epithelium of the nasal cavity. There are ~400 and ~1,000 OR genes in the human and mouse genomes, respectively, forming the largest multigene family in mammals. The relationships between ORs and odorants are multiple-to-multiple, which allows for discriminating almost unlimited number of different odorants by a combination of ORs. However, the OR-ligand relationships are still largely unknown, and predicting the quality of odor from its molecular structure is unsuccessful.Extensive bioinformatic analyses using the whole genomes of various organisms revealed a great variation in number of OR genes among species, reflecting the diversity of their living environments. For example, higher primates equipped with a well-developed vision system and dolphins that are secondarily adapted to the aquatic life have considerably smaller numbers of OR genes than most of other mammals do. OR genes are characterized by extremely frequent gene duplications and losses. The OR gene repertories are also diverse among human individuals, explaining the diversity of odor perception such as the specific anosmia.OR genes are present in all vertebrates. The number of OR genes is smaller in teleost fishes than in mammals, while the diversity is higher in the former than the latter. Because the genome of amphioxus, the most basal chordate species, harbors vertebrate-like OR genes, the origin of OR genes can be traced back to the common ancestor of the phylum Chordata.

hzAnalyzer: detection, quantification, and visualization of contiguous homozygosity in high-density genotyping datasets.

The analysis of contiguous homozygosity (runs of homozygous loci) in human genotyping datasets is critical in the search for causal disease variants in monogenic disorders, studies of population history and the identification of targets of natural selection. Here, we report methods for extracting homozygous segments from high-density genotyping datasets, quantifying their local genomic structure, identifying outstanding regions within the genome and visualizing results for comparative analysis between population samples.

Difference in gene duplicability may explain the difference in overall structure of protein-protein interaction networks among eukaryotes.

BACKGROUND: A protein-protein interaction network (PIN) was suggested to be a disassortative network, in which interactions between high- and low-degree nodes are favored while hub-hub interactions are suppressed. It was postulated that a disassortative structure minimizes unfavorable cross-talks between different hub-centric functional modules and was positively selected in evolution. However, by re-examining yeast PIN data, several researchers reported that the disassortative structure observed in a PIN might be an experimental artifact. Therefore, the existence of a disassortative structure and its possible evolutionary mechanism remains unclear. RESULTS: In this study, we investigated PINs from the yeast, worm, fly, human, and malaria parasite including four different yeast PIN datasets. The analyses showed that the yeast, worm, fly, and human PINs are disassortative while the malaria parasite PIN is not. By conducting simulation studies on the basis of a duplication-divergence model, we demonstrated that a preferential duplication of low- and high-degree nodes can generate disassortative and non-disassortative networks, respectively. From this observation, we hypothesized that the difference in degree dependence on gene duplications accounts for the difference in assortativity of PINs among species. Comparison of 55 proteomes in eukaryotes revealed that genes with lower degrees showed higher gene duplicabilities in the yeast, worm, and fly, while high-degree genes tend to have high duplicabilities in the malaria parasite, supporting the above hypothesis. CONCLUSIONS: These results suggest that disassortative structures observed in PINs are merely a byproduct of preferential duplications of low-degree genes, which might be caused by an organism's living environment.