Evolution of olfactory receptor genes in the human genome

DNA damage drives antigen diversification through mosaic Variant Surface Glycoprotein (VSG) formation in Trypanosoma brucei

Antigenic variation, using large genomic repertoires of antigen-encoding genes, allows pathogens to evade host antibody. Many pathogens, including the African trypanosome Trypanosoma brucei, extend their antigenic repertoire through genomic diversification. While evidence suggests that T. brucei depends on the generation of new variant surface glycoprotein (VSG) genes to maintain a chronic infection, a lack of experimentally tractable tools for studying this process has obscured its underlying mechanisms. Here, we present a highly sensitive targeted sequencing approach for measuring VSG diversification. Using this method, we demonstrate that a Cas9-induced DNA double-strand break within the VSG coding sequence can induce VSG recombination with patterns identical to those observed during infection. These newly generated VSGs are antigenically distinct from parental clones and thus capable of facilitating immune evasion. Together, these results provide insight into the mechanisms of VSG diversification and an experimental framework for studying the evolution of antigen repertoires in pathogenic microbes.

NAD activates olfactory receptor 1386 to regulate type I interferon responses in Plasmodium yoelii YM infection

Olfactory receptors (Olfr) are G protein-coupled receptors that are normally expressed on olfactory sensory neurons to detect volatile chemicals or odorants. Interestingly, many Olfrs are also expressed in diverse tissues and function in cell-cell recognition, migration, and proliferation as well as immune responses and disease processes. Here, we showed that many Olfr genes were expressed in the mouse spleen, linked to Plasmodium yoelii genetic loci significantly, and/or had genome-wide patterns of LOD scores (GPLSs) similar to those of host Toll-like receptor genes. Expression of specific Olfr genes such as Olfr1386 in HEK293T cells significantly increased luciferase signals driven by IFN-β and NF-κB promoters, with elevated levels of phosphorylated TBK1, IRF3, P38, and JNK. Mice without Olfr1386 were generated using the CRISPR/Cas9 method, and the Olfr1386-/- mice showed significantly lower IFN-α/β levels and longer survival than wild-type (WT) littermates after infection with P. yoelii YM parasites. Inhibition of G protein signaling and P38 activity could affect cyclic AMP-responsive element promoter-driven luciferase signals and IFN-β mRNA levels in HEK293T cells expressing the Olfr1386 gene, respectively. Screening of malaria parasite metabolites identified nicotinamide adenine dinucleotide (NAD) as a potential ligand for Olfr1386, and NAD could stimulate IFN-β responses and phosphorylation of TBK1 and STAT1/2 in RAW264.7 cells. Additionally, parasite RNA (pRNA) could significantly increase Olfr1386 mRNA levels. This study links multiple Olfrs to host immune response pathways, identifies a candidate ligand for Olfr1386, and demonstrates the important roles of Olfr1386 in regulating type I interferon (IFN-I) responses during malaria parasite infections.

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.

Unique transcriptomes of sensory and non-sensory neurons: insights from Splicing Regulatory States

Alternative Splicing (AS) programs serve as instructive signals of cell type specificity, particularly within the brain, which comprises dozens of molecularly and functionally distinct cell types. Among them, retinal photoreceptors stand out due to their unique transcriptome, making them a particularly well-suited system for studying how AS shapes cell type-specific molecular functions. Here, we use the Splicing Regulatory State (SRS) as a novel framework to discuss the splicing factors governing the unique AS pattern of photoreceptors, and how this pattern may aid in the specification of their highly specialized sensory cilia. In addition, we discuss how other sensory cells with ciliated structures, for which data is much scarcer, also rely on specific SRSs to implement a proteome specialized in the detection of sensory stimuli. By reviewing the general rules of cell type- and tissue-specific AS programs, firstly in the brain and subsequently in specialized sensory neurons, we propose a novel paradigm on how SRSs are established and how they can diversify. Finally, we illustrate how SRSs shape the outcome of mutations in splicing factors to produce cell type-specific phenotypes that can lead to various human diseases.

Localized Expression of Olfactory Receptor Genes in the Olfactory Organ of Common Minke Whales

Baleen whales (Mysticeti) possess the necessary anatomical structures and genetic elements for olfaction. Nevertheless, the olfactory receptor gene (OR) repertoire has undergone substantial degeneration in the cetacean lineage following the divergence of the Artiodactyla and Cetacea. The functionality of highly degenerated mysticete ORs within their olfactory epithelium remains unknown. In this study, we extracted total RNA from the nasal mucosae of common minke whales (Balaenoptera acutorostrata) to investigate ORs' localized expression. All three sections of the mucosae examined in the nasal chamber displayed comparable histological structure. However, the posterior portion of the frontoturbinal region exhibited notably high OR expression. Neither the olfactory bulb nor the external skin exhibited the expression of these genes. Although this species possesses four intact non-class-2 ORs, all the ORs expressed in the nasal mucosae belong to class-2, implying the loss of aversion to specific odorants. These anatomical and genomic analyses suggest that ORs are still responsible for olfaction within the nasal region of baleen whales, enabling them to detect desirable scents such as prey and potential mating partners.

M2OR: a database of olfactory receptor-odorant pairs for understanding the molecular mechanisms of olfaction

Mammalian sense of smell is triggered by interaction between odorant molecules and a class of proteins, called olfactory receptors (ORs). These receptors, expressed at the surface of olfactory sensory neurons, encode myriad of distinct odors via a sophisticated activation pattern. However, determining the molecular recognition spectrum of ORs remains a major challenge. The Molecule to Olfactory Receptor database (M2OR, https://m2or.chemsensim.fr/) provides curated data that allows an easy exploration of the current state of the research on OR-molecule interaction. We have gathered a database of 75,050 bioassay experiments for 51 395 distinct OR-molecule pairs. Drawn from published literature and public databases, M2OR contains information about OR responses to molecules and their mixtures, receptor sequences and experimental details. Users can obtain information on the activity of a chosen molecule or a group of molecules, or search for agonists for a specific OR or a group of ORs. Advanced search allows for fine-grained queries using various metadata such as species or experimental assay system, and the database can be queried by multiple inputs via a batch search. Finally, for a given search query, users can access and download a curated aggregation of the experimental data into a binarized combinatorial code of olfaction.

Tracking the Diversity and Chromosomal Distribution of the Olfactory Receptor Gene Repertoires of Three Anurans Species

Olfaction is a crucial capability for most vertebrates and is realized through olfactory receptors in the nasal cavity. The enormous diversity of olfactory receptors has been created by gene duplication, following a birth-and-death model of evolution. The olfactory receptor genes of the amphibians have received relatively little attention up to now, although recent studies have increased the number of species for which data are available. This study analyzed the diversity and chromosomal distribution of the OR genes of three anuran species (Engystomops pustulosus, Bufo bufo and Hymenochirus boettgeri). The OR genes were identified through searches for homologies, and sequence filtering and alignment using bioinformatic tools and scripts. A high diversity of OR genes was found in all three species, ranging from 917 in B. bufo to 1194 in H. boettgeri, and a total of 2076 OR genes in E. pustulosus. Six OR groups were recognized using an evolutionary gene tree analysis. While E. pustulosus has one of the highest numbers of genes of the gamma group (which detect airborne odorants) yet recorded in an anuran, B. bufo presented the smallest number of pseudogene sequences ever identified, with no pseudogenes in either the beta or epsilon groups. Although H. boettgeri shares many morphological adaptations for an aquatic lifestyle with Xenopus, and presented a similar number of genes related to the detection of water-soluble odorants, it had comparatively far fewer genes related to the detection of airborne odorants. This study is the first to describe the complete OR repertoire of the three study species and represents an important contribution to the understanding of the evolution and function of the sense of smell in vertebrates.

Recent positive selection signatures reveal phenotypic evolution in the Han Chinese population

Characterizing natural selection signatures and relationships with phenotype spectra is important for understanding human evolution and both biological and pathological mechanisms. Here, we identified 24 genetic loci under recent selection by analyzing rare singletons in 3946 high-depth whole-genome sequencing data of Han Chinese. The loci include immune-related gene regions (MHC cluster, IGH cluster, STING1, and PSG), alcohol metabolism-related gene regions (ADH1B, ALDH2, and ALDH3B2), and the olfactory perception gene OR4C16, in which the MHC cluster, ADH1B, and ALDH2 were also identified by TOPMed and WestLake Biobank. Among the signals, the IGH cluster is particularly interesting, in which the favored allele of variant 14_105737776_C_T (rs117518546, IgG1-G396R) promotes immune response, but also increases the risk of an autoimmune disease systemic lupus erythematosus (SLE). It is also surprising that our newly discovered ALDH3B2 evolved in the opposite direction to ALDH2 for alcohol metabolism. Besides monogenic traits, we found that multiple complex traits experienced polygenic adaptation. Particularly, multi-methods consistently revealed that lower blood pressure was favored in natural selection. Finally, we built a database named RePoS (recent positive selection, http://bigdata.ibp.ac.cn/RePoS/) to integrate and display multi-population selection signals. Our study extended our understanding of natural evolution and phenotype adaptation in Han Chinese as well as other populations.

Insights from a chum salmon (Oncorhynchus keta) genome assembly regarding whole-genome duplication and nucleotide variation influencing gene function

Chum salmon are ecologically important to Pacific Ocean ecosystems and commercially important to fisheries. To improve the genetic resources available for this species, we sequenced and assembled the genome of a male chum salmon using Oxford Nanopore read technology and the Flye genome assembly software (contig N50: ∼2 Mbp, complete BUSCOs: ∼98.1%). We also resequenced the genomes of 59 chum salmon from hatchery sources to better characterize the genome assembly and the diversity of nucleotide variants impacting phenotype variation. With genomic sequences from a doubled haploid individual, we were able to identify regions of the genome assembly that have been collapsed due to high sequence similarity between homeologous (duplicated) chromosomes. The homeologous chromosomes are relics of an ancient salmonid-specific genome duplication. These regions were enriched with genes whose functions are related to the immune system and responses to toxins. From analyzing nucleotide variant annotations of the resequenced genomes, we were also able to identify genes that have increased levels of variants thought to moderately impact gene function. Genes related to the immune system and the detection of chemical stimuli (olfaction) had increased levels of these variants based on a gene ontology enrichment analysis. The tandem organization of many of the enriched genes raises the question of why they have this organization.

Specific olfactory deficit patterns observed in seniors and associated with cognitive decline

Olfactory dysfunction is a common symptom in neurodegenerative disorders and is regarded as a potential early predictor of impending cognitive decline. This study was undertaken in order to determine if olfactory dysfunction observed in the elderly is due to a general loss of smell or the inability to detect specific odours, and if misidentification of odours correlates with cognitive scores. Seniors for the Olfactory Response and Cognition in Aging (ORCA) sub-study were recruited from the Quebec Nutrition and Successful Aging (NuAge) cohort. The University of Pennsylvania smell identification test (UPSIT) was performed to measure olfactory function and the telephone Mini Mental State Examination (t-MMSE) and the French version of the Telephone Interview for Cognitive Status Modified (F-TICS-m) for cognitive status. The results demonstrate that seniors exhibit specific olfactory loss and had severe difficulty in particular in identifying lemon, pizza, fruit punch, cheddar cheese and lime. Furthermore, there was a significant difference in the ability to detect certain odours between the sexes. Results also showed that misidentification of certain scents was associated with cognitive scores, and when the sexes were assessed separately sex-specific misidentification of cognitive-associated odours was observed. The relationship between the cognitive scores and scent misidentification suggests that impending cognitive decline may be highlighted by the inability to smell specific odours. Our study provides additional support for the testing of olfactory function in the elderly and suggests that loss of smell for particular scents may become a useful diagnostic tool.

Genome-Wide Identification and Characterization of Olfactory Receptor Genes in Silver Sillago (Sillago sihama)

Olfactory receptor (OR) genes are essential in the specific recognition of diverse stimuli in fish. In this study, a total of 141 OR genes were identified in silver sillago (Sillago sihama), a marine fish sensitive to environmental stimuli, including 112 intact genes, 26 truncated genes, and three pseudogenes. A phylogenetic tree analysis elucidated that the OR genes of S. sihama were classified into six groups, of which β, γ, δ, ε, and ζ groups belonged to type I, and the η group belonged to type II. The type I OR genes contained almost all conserved motifs (n = 62), while type II OR genes mainly retained conserved motifs 7(3), 1, 10, 4, and 2 (n = 39). OR genes were mainly distributed on LG1, LG9, LG11, and LG12. Of all OR genes, 36.23% (50 genes) showed significant expansion in S. sihama. Ka/Ks analysis demonstrated that 227 sites were under purifying selection, while 12 sites were under positive selection, including eight genes in the OR2A12 gene subfamily. Sixty-one genes (44.20%) displayed differential expression under hypoxic stress. The identified OR genes explored the mechanism of environmental stress and ecological adaptation of S. sihama, and provided valuable genomic resources for further research on the olfaction of teleosts.

Repetitive Direct Comparison Method for Odor Sensing

Olfactory sensors are one of the most anticipated applications of gas sensors. To distinguish odors-complex mixtures of gas species, it is necessary to extract sensor responses originating from the target odors. However, the responses of gas sensors tend to be affected by interfering gases with much higher concentrations than target odor molecules. To realize practical applications of olfactory sensors, extracting minute sensor responses of odors from major interfering gases is required. In this study, we propose a repetitive direct comparison (rDC) method, which can highlight the difference in odors by alternately injecting the two target odors into a gas sensor. We verified the feasibility of the rDC method on chocolates with two different flavors by using a sensor system based on membrane-type surface stress sensors (MSS). The odors of the chocolates were measured by the rDC method, and the signal-to-noise ratios (S/N) of the measurements were evaluated. The results showed that the rDC method achieved improved S/N compared to a typical measurement. The result also indicates that sensing signals could be enhanced for a specific combination of receptor materials of MSS and target odors.

Comparative Analysis of Olfactory Receptor Repertoires Sheds Light on the Diet Adaptation of the Bamboo-Eating Giant Panda Based on the Chromosome-Level Genome

The giant panda (Ailuropoda melanoleuca) is the epitome of a flagship species for wildlife conservation and also an ideal model of adaptive evolution. As an obligate bamboo feeder, the giant panda relies on the olfaction for food recognition. The number of olfactory receptor (OR) genes and the rate of pseudogenes are the main factors affecting the olfactory ability of animals. In this study, we used the chromosome-level genome of the giant panda to identify OR genes and compared the genome sequences of OR genes with five other Ursidae species (spectacled bear (Tremarctos ornatus), American black bear (Ursus americanus), brown bear (Ursus arctos), polar bear (Ursus maritimus) and Asian black bear (Ursus thibetanus)). The giant panda had 639 OR genes, including 408 functional genes, 94 partial OR genes and 137 pseudogenes. Among them, 222 OR genes were detected and distributed on 18 chromosomes, and chromosome 8 had the most OR genes. A total of 448, 617, 582, 521 and 792 OR genes were identified in the spectacled bear, American black bear, brown bear, polar bear and Asian black bear, respectively. Clustering analysis based on the OR protein sequences of the six species showed that the OR genes distributed in 69 families and 438 subfamilies based on sequence similarity, and the six mammals shared 72 OR gene subfamilies, while the giant panda had 31 unique OR gene subfamilies (containing 35 genes). Among the 35 genes, there are 10 genes clustered into 8 clusters with 10 known human OR genes (OR8J3, OR51I1, OR10AC1, OR1S2, OR1S1, OR51S1, OR4M1, OR4M2, OR51T1 and OR5W2). However, the kind of odor molecules can be recognized by the 10 known human OR genes separately, which needs further research. The phylogenetic tree showed that 345 (about 84.56%) functional OR genes were clustered as Class-II, while only 63 (about 15.44%) functional OR genes were clustered as Class-I, which required further and more in-depth research. The potential odor specificity of some giant panda OR genes was identified through the similarity to human protein sequences. Sequences similar to OR2B1, OR10G3, OR11H6 and OR11H7P were giant panda-specific lacking, which may be related to the transformation and specialization from carnivore to herbivore of the giant panda. Since our reference to flavoring agents comes from human research, the possible flavoring agents from giant panda-specific OR genes need further investigation. Moreover, the conserved motifs of OR genes were highly conserved in Ursidae species. This systematic study of OR genes in the giant panda will provide a solid foundation for further research on the olfactory function and variation of the giant panda.

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.

Building the Chordata Olfactory Receptor Database using more than 400,000 receptors annotated by Genome2OR

Olfactory receptors are poorly annotated for most genome-sequenced chordates. To address this deficiency, we developed a nhmmer-based olfactory receptor annotation tool Genome2OR ( https://github.com/ToHanwei/Genome2OR.git ), and used it to process 1,695 sequenced chordate genomes in the NCBI Assembly database as of January, 2021. In total, 765,248 olfactory receptor genes were annotated, with 404,426 functional genes and 360,822 pseudogenes, which represents a four-fold increase in the number of annotated olfactory receptors. Based on the annotation data, we built a database called Chordata Olfactory Receptor Database (CORD, https://cord.ihuman.shanghaitech.edu.cn ) for archiving, analysing and disseminating the data. Beyond the primary data, we offer derivative information, including pictures of species, cross references to public databases, structural models, sequence similarity networks and sequence profiles in the CORD. Furthermore, we did brief analyses on these receptors, including building a huge protein sequence similarity network covering all receptors in the database, and clustering them into 20 communities, classifying the 20 communities into three categories based on their presences/absences in ray-finned fish and/or lobe-finned fish. We infer that olfactory receptors should have unique activation and desensitization mechanisms by analysing their sequences and structural models. We believe the CORD can benefit the researchers and the general public who are interested in olfaction.

Functional Divergence in a Multi-gene Family Is a Key Evolutionary Innovation for Anaerobic Growth in Saccharomyces cerevisiae

The amplification and diversification of genes into large multi-gene families often mark key evolutionary innovations, but this process often creates genetic redundancy that hinders functional investigations. When the model budding yeast Saccharomyces cerevisiae transitions to anaerobic growth conditions, the cell massively induces the expression of seven serine/threonine-rich anaerobically-induced cell wall mannoproteins (anCWMPs): TIP1, TIR1, TIR2, TIR3, TIR4, DAN1, and DAN4. Here, we show that these genes likely derive evolutionarily from a single ancestral anCWMP locus, which was duplicated and translocated to new genomic contexts several times both prior to and following the budding yeast whole genome duplication (WGD) event. Based on synteny and their phylogeny, we separate the anCWMPs into four gene subfamilies. To resolve prior inconclusive genetic investigations of these genes, we constructed a set of combinatorial deletion mutants to determine their contributions toward anaerobic growth in S. cerevisiae. We found that two genes, TIR1 and TIR3, were together necessary and sufficient for the anCWMP contribution to anaerobic growth. Overexpressing either gene alone was insufficient for anaerobic growth, implying that they encode non-overlapping functional roles in the cell during anaerobic growth. We infer from the phylogeny of the anCWMP genes that these two important genes derive from an ancient duplication that predates the WGD event, whereas the TIR1 subfamily experienced gene family amplification after the WGD event. Taken together, the genetic and molecular evidence suggests that one key anCWMP gene duplication event, several auxiliary gene duplication events, and functional divergence underpin the evolution of anaerobic growth in budding yeasts.

Identification and characterization of olfactory receptor genes and olfactory perception in rapa whelk Rapana venosa (Valenciennes, 1846) during larval settlement and metamorphosis

The rapa whelk Rapana venosa, an economically important marine fishery resource in China but a major invader all over the world, changes from a phytophagous to a carnivorous form following settlement and metamorphosis. However, the low settlement and metamorphosis rates (<1%) of larvae limit the abundance of R. venosa. This critical step (settlement and metamorphosis) remains poorly characterized but may be related to how larvae perceive the presence of shellfish, their new source of food. Here, we report that larvae may use olfactory perception to sense shellfish. Olfactory receptor (OR) genes are involved in odor sensing in animals. We identified a total of 463 OR genes, which could be grouped into nine clades based on phylogenetic analysis. When assessing the attraction of larvae at different developmental stages to oyster odor, R. venosa showed active settlement and metamorphosis behavior only at the J4 stage (competent larva, 1000-1500 μm shell length) and in the presence of shellfish odor at the same time. Expression of OR gene family members differed between stage 2 (four-spiral whorl stage) and stage 1 (single- to three-spiral whorl stage), indicating significant changes in the olfactory system during larval development. These findings broaden our understanding of olfactory perception, settlement, and metamorphosis in gastropods and can be used to improve R. venosa harvesting, as well as the sustainable development and utilization of this resource.

Characterization of olfactory receptor repertoires provides insights into the high-altitude adaptation of the yak based on the chromosome-level genome

Olfaction in vertebrates plays pivotal parts in many aspects, such as localizing prey or food, mating behavior, avoiding predators, and social communication. Yak (Bos grunniens) is the only Bos species that can thrive in high-altitude areas. In view of the critical role of olfactory receptors (ORs) in the specific recognition of diverse stimuli, investigating the evolutionary dynamics of ORs in the yak means a lot. In this study, we used the chromosome-level genome of the yak to identify the ORs genes and discussed the effects of high altitude on the yak's olfaction by comparing the yak with other low-altitude living Bos species (Bos frontalis (gayal), Bos gaurus (gaur), Bos indicus (zebu) and Bos taurus (cattle)). The yak had 400 OR genes, including 264 functional genes, 16 partial genes and 120 OR pseudo genes. There were 387 OR genes mapped to yak 31 chromosomes, and chromosomes 13 and 8 had the most OR genes and functional OR genes. Among these five Bos species, yak had the least number of OR gene subfamilies, OR genes and functional OR genes, while the total number of OR genes in gayal (n = 784) was almost twice as many as that of yak, indicating that the olfaction of yak may be less developed. In addition, the phylogenetic relationships of the functional Bos OR genes were illustrated, which comprised 79 families and 466 subfamilies distributed in two classes (Class I and Class II). There were 76 OR gene subfamilies shared by these five Bos species and 17 OR gene subfamilies were unique to the yak. The potential odor specificity of 44 yak OR genes was identified through the similarity to human OR protein sequences. Remarkably, yak lacks β-ionone and Isovaleric acid(IVA)-related ORs, which may be related to the decline of high-altitude herbaceous plant diversity and underdeveloped yak sweat glands. The conserved motifs of OR genes were highly conserved in Bos species. These results provided a solid foundation for further studies on the molecular mechanisms of the yak's adaptation to the high-altitude environment in olfaction.

Shared Mechanisms for Mutually Exclusive Expression and Antigenic Variation by Protozoan Parasites

Cellular decision-making at the level of gene expression is a key process in the development and evolution of every organism. Variations in gene expression can lead to phenotypic diversity and the development of subpopulations with adaptive advantages. A prime example is the mutually exclusive activation of a single gene from within a multicopy gene family. In mammals, this ranges from the activation of one of the two immunoglobulin (Ig) alleles to the choice in olfactory sensory neurons of a single odorant receptor (OR) gene from a family of more than 1,000. Similarly, in parasites like Trypanosoma brucei, Giardia lamblia or Plasmodium falciparum, the process of antigenic variation required to escape recognition by the host immune system involves the monoallelic expression of vsg, vsp or var genes, respectively. Despite the importance of this process, understanding how this choice is made remains an enigma. The development of powerful techniques such as single cell RNA-seq and Hi-C has provided new insights into the mechanisms these different systems employ to achieve monoallelic gene expression. Studies utilizing these techniques have shown how the complex interplay between nuclear architecture, physical interactions between chromosomes and different chromatin states lead to single allele expression. Additionally, in several instances it has been observed that high-level expression of a single gene is preceded by a transient state where multiple genes are expressed at a low level. In this review, we will describe and compare the different strategies that organisms have evolved to choose one gene from within a large family and how parasites employ this strategy to ensure survival within their hosts.

Olfactory receptor coding sequences cause silencing of episomal constructs in multiple cell lines

The mammalian olfactory system consists of sensory neurons with specialized odorant-binding capability accomplished by mutually exclusive odorant receptor (OR) expression. Mutually exclusive OR expression is a complex multi-step process regulated by a number of cis and trans factors, including pan-silencing of all OR genes preceding the robust and stable expression of the one OR selected in each sensory neuron. We transfected two olfactory-placode-derived cell lines modeling immature odorant sensory neurons, as well as the GD25 fibroblast cell line, with episomes containing CMV-driven GFP and TK-driven hygromycin reporter genes. We inserted various coding sequences, along with an IRES, immediately upstream of the GFP gene to produce bicistronic mRNAs driven from the local CMV promoter. We found that the presence of several OR coding sequences resulted in significantly diminished episomal expression of GFP in all three cell lines. These findings suggest that OR coding sequences have intrinsic self-silencing capability that might facilitate mutually exclusive OR expression in olfactory sensory neurons by making it less likely that multiple ORs acquire an above-threshold level of expression at once.

Odorant metabolism of the olfactory cleft mucus in idiopathic olfactory impairment patients and healthy volunteers

BACKGROUND

It remains unclear whether the metabolic activity of nasal mucus in the olfactory and respiratory areas is different. Moreover, age- and olfaction-related changes may affect metabolism.

METHODS

Hexanal, octanal, and 2-methylbutanal were selected for in vitro metabolism analysis and compared between the olfactory cleft and respiratory mucus of participants < 50-year-old with normal olfaction using gas chromatography mass spectrometry. The metabolic activity of hexanal in the olfactory cleft mucus was further compared between three groups, (1) normal olfaction, age < 50 years old, (2) normal olfaction, age ≥50 years old, and (3) idiopathic olfactory impairment. To characterize the enzyme(s) responsible for aldehyde reduction, we also tested if epalr22897estat and 3,5-dichlorosalicylic acid, types of reductase inhibitors, affect metabolism.

RESULTS

Conversion of aldehydes to their corresponding alcohols was observed in the olfactory cleft and respiratory mucus. The metabolic production of hexanol, octanol, and 2-methybutanol was significantly higher in the olfactory cleft mucus than in the respiratory mucus (p < 0.01). The metabolic conversion of hexanal to hexanol in the mucus of the idiopathic olfactory impairment group was significantly lower than that in the age-matched normal olfaction group. Excluding the nicotinamide adenine dinucleotide phosphate (NADPH) regenerating system from the reaction mixture inhibited metabolism. The addition of either epalr22897estat or 3,5-dichlorosalicylic acid did not inhibit this metabolic conversion.

CONCLUSIONS

The enzymatic metabolism of odorants in the olfactory cleft mucus is markedly higher than in the respiratory mucus and decreases in patients with idiopathic olfactory impairment.

OdoriFy: A conglomerate of artificial intelligence-driven prediction engines for olfactory decoding

The molecular mechanisms of olfaction, or the sense of smell, are relatively underexplored compared with other sensory systems, primarily because of its underlying molecular complexity and the limited availability of dedicated predictive computational tools. Odorant receptors (ORs) allow the detection and discrimination of a myriad of odorant molecules and therefore mediate the first step of the olfactory signaling cascade. To date, odorant (or agonist) information for the majority of these receptors is still unknown, limiting our understanding of their functional relevance in odor-induced behavioral responses. In this study, we introduce OdoriFy, a Web server featuring powerful deep neural network-based prediction engines. OdoriFy enables (1) identification of odorant molecules for wildtype or mutant human ORs (Odor Finder); (2) classification of user-provided chemicals as odorants/nonodorants (Odorant Predictor); (3) identification of responsive ORs for a query odorant (OR Finder); and (4) interaction validation using Odorant-OR Pair Analysis. In addition, OdoriFy provides the rationale behind every prediction it makes by leveraging explainable artificial intelligence. This module highlights the basis of the prediction of odorants/nonodorants at atomic resolution and for the ORs at amino acid levels. A key distinguishing feature of OdoriFy is that it is built on a comprehensive repertoire of manually curated information of human ORs with their known agonists and nonagonists, making it a highly interactive and resource-enriched Web server. Moreover, comparative analysis of OdoriFy predictions with an alternative structure-based ligand interaction method revealed comparable results. OdoriFy is available freely as a web service at https://odorify.ahujalab.iiitd.edu.in/olfy/.

Improved inference of tandem domain duplications

MOTIVATION

Protein domain duplications are a major contributor to the functional diversification of protein families. These duplications can occur one at a time through single domain duplications, or as tandem duplications where several consecutive domains are duplicated together as part of a single evolutionary event. Existing methods for inferring domain-level evolutionary events are based on reconciling domain trees with gene trees. While some formulations consider multiple domain duplications, they do not explicitly model tandem duplications; this leads to inaccurate inference of which domains duplicated together over the course of evolution.

RESULTS

Here, we introduce a reconciliation-based framework that considers the relative positions of domains within extant sequences. We use this information to uncover tandem domain duplications within the evolutionary history of these genes. We devise an integer linear programming approach that solves our problem exactly, and a heuristic approach that works well in practice. We perform extensive simulation studies to demonstrate that our approaches can accurately uncover single and tandem domain duplications, and additionally test our approach on a well-studied orthogroup where lineage-specific domain expansions exhibit varying and complex domain duplication patterns.

AVAILABILITY AND IMPLEMENTATION

Code is available on github at https://github.com/Singh-Lab/TandemDuplications.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Limbic Expression of mRNA Coding for Chemoreceptors in Human Brain-Lessons from Brain Atlases

Animals strongly rely on chemical senses to uncover the outside world and adjust their behaviour. Chemical signals are perceived by facial sensitive chemosensors that can be clustered into three families, namely the gustatory (TASR), olfactory (OR, TAAR) and pheromonal (VNR, FPR) receptors. Over recent decades, chemoreceptors were identified in non-facial parts of the body, including the brain. In order to map chemoreceptors within the encephalon, we performed a study based on four brain atlases. The transcript expression of selected members of the three chemoreceptor families and their canonical partners was analysed in major areas of healthy and demented human brains. Genes encoding all studied chemoreceptors are transcribed in the central nervous system, particularly in the limbic system. RNA of their canonical transduction partners (G proteins, ion channels) are also observed in all studied brain areas, reinforcing the suggestion that cerebral chemoreceptors are functional. In addition, we noticed that: (i) bitterness-associated receptors display an enriched expression, (ii) the brain is equipped to sense trace amines and pheromonal cues and (iii) chemoreceptor RNA expression varies with age, but not dementia or brain trauma. Extensive studies are now required to further understand how the brain makes sense of endogenous chemicals.

A chromosome-level assembly of blunt snout bream (Megalobrama amblycephala) reveals an expansion of olfactory receptor genes in freshwater fish

The number of olfactory receptor genes (ORs), which are responsible for detecting diverse odor molecules varies extensively among mammals as a result of frequent gene gains and losses that contribute to olfactory specialization. However, how OR expansions/contractions in fish are influenced by habitat and feeding habit and which OR subfamilies are important in each ecological niche is unknown. Here, we report a major OR expansion in a freshwater herbivorous fish, Megalobrama amblycephala, using a highly contiguous, chromosome-level assembly. We evaluate the possible contribution of OR expansion to habitat and feeding specialization by comparing the OR repertoire in 28 phylogenetically and ecologically diverse teleosts. In total, we analyzed > 4,000 ORs including 3,253 intact, 122 truncated, and 913 pseudogenes. The number of intact ORs is highly variable ranging from 20 to 279. We estimate that the most recent common ancestor of Osteichthyes had 62 intact ORs, which declined in most lineages except the freshwater Otophysa clade that has a substantial expansion in subfamily β and ε ORs. Across teleosts, we found a strong association between duplications of β and ε ORs and freshwater habitat. Nearly, all ORs were expressed in the olfactory epithelium (OE) in three tested fish species. Specifically, all the expanded β and ε ORs were highly expressed in OE of M. amblycephala. Together, we provide molecular and functional evidence for how OR repertoires in fish have undergone gain and loss with respect to ecological factors and highlight the role of β and ε OR in freshwater adaptation.

The Transposable Element Environment of Human Genes Differs According to Their Duplication Status and Essentiality

Transposable elements (TEs) are major components of eukaryotic genomes and represent approximately 45% of the human genome. TEs can be important sources of novelty in genomes and there is increasing evidence that TEs contribute to the evolution of gene regulation in mammals. Gene duplication is an evolutionary mechanism that also provides new genetic material and opportunities to acquire new functions. To investigate how duplicated genes are maintained in genomes, here, we explored the TE environment of duplicated and singleton genes. We found that singleton genes have more short-interspersed nuclear elements and DNA transposons in their vicinity than duplicated genes, whereas long-interspersed nuclear elements and long-terminal repeat retrotransposons have accumulated more near duplicated genes. We also discovered that this result is highly associated with the degree of essentiality of the genes with an unexpected accumulation of short-interspersed nuclear elements and DNA transposons around the more-essential genes. Our results underline the importance of taking into account the TE environment of genes to better understand how duplicated genes are maintained in genomes.

Structural characterization and duplication modes of pseudogenes in plants

We identified and characterized the pseudogene complements of five plant species: four dicots (Arabidopsis thaliana, Vitis vinifera, Populus trichocarpa and Phaseolus vulgaris) and one monocot (Oryza sativa). Retroposition was considered of modest importance for pseudogene formation in all investigated species except V. vinifera, which showed an unusually high number of retro-pseudogenes in non coding genic regions. By using a pipeline for the classification of sequence duplicates in plant genomes, we compared the relative importance of whole genome, tandem, proximal, transposed and dispersed duplication modes in the pseudo and functional gene complements. Pseudogenes showed higher tendencies than functional genes to genomic dispersion. Dispersed pseudogenes were prevalently fragmented and showed high sequence divergence at flanking regions. On the contrary, those deriving from whole genome duplication were proportionally less than expected based on observations on functional loci and showed higher levels of flanking sequence conservation than dispersed pseudogenes. Pseudogenes deriving from tandem and proximal duplications were in excess compared to functional loci, probably reflecting the high evolutionary rate associated with these duplication modes in plant genomes. These data are compatible with high rates of sequence turnover at neutral sites and double strand break repairs mediated duplication mechanisms.

Characterization of Olfactory Receptor Repertoires in the Endangered Snow Leopard Based on the Chromosome-Level Genome

Olfaction is a complicated process that begins with the specific binding of volatile odorant molecules to dedicated olfactory receptors (ORs) in the olfactory epithelium and plays a pivotal role in the survival of mammals. The OR subgenome of the snow leopard has remained largely unexplored, and thus, investigation of the OR system would shed light on the evolutionary dynamics of the snow leopard OR repertoires and genetic evidence for environmental adaptation. In this study, we conducted genome-wide identification and characterization of OR genes in the snow leopard and compared them to all other Panthera species. A total of 213, 294, 624, 305, and 253 functional OR genes were identified in the snow leopard, lion, jaguar, leopard, and tiger, respectively. The phylogenetic relationships of functional Panthera OR genes were illustrated, which comprised 69 families and 350 subfamilies distributed in two classes (Class I and Class II). Comparative analysis of the five Panthera species indicated 115 shared and 5 snow leopard-specific clusters. The potential odorant specificity of certain snow leopard OR genes was identified by similarities to human protein sequences and we identified odorants such as eugenol methyl ether that had the most OR genes. Since our references for odorants were from human studies, possible odorants from snow leopard-specific OR genes need further investigation. The lowest number of OR genes for the snow leopard among Panthera species possibly revealed the association between OR gene family contraction and high-altitude adaptation, which needed further and deeper investigation. This systematic study of OR genes in the snow leopard will provide a solid foundation for further study of olfactory function and variation in the snow leopard.

Genome-wide identification and characterization of olfactory receptor genes in common carp (Cyprinus carpio)

The environment contains a large extent of chemical information, which could be detected as olfactory sense. Olfactory in vertebrates plays important roles on many aspects during life time, including localizing prey or food, avoiding predators, mating behavior and social communication. Considering the essential role of olfactory receptors in the specific recognition of diverse stimuli, understanding the evolutionary dynamics of olfactory receptors in teleost means a lot, especially in the allotetraploid common carp, who has undergone the fourth whole-genome duplication event. Here, we identified the whole set of olfactory receptor genes in representative teleosts and found a significant contraction in common carp when compared with other teleosts. Odorant receptor genes (OR) occupy the most among four groups of olfactory receptors, including 33 functional genes and 16 pseudogenes. Furthermore, 6 trace amine-associated receptor (TAAR) genes (including 1 pseudogene), 7 odorant-related-A receptor genes, and 10 olfactory C family receptor genes (including 3 pseudogenes) were identified in common carp. Phylogenetic and motif analysis were performed to illustrate the phylogenetic relationship and structural conservation of teleost olfactory receptors. Selection pressure analysis suggested that olfactory receptor groups in common carp were all under relaxed purifying-selection. Additionally, gene expression divergences for olfactory receptor genes were investigated during embryonic development stages of common carp. We aim to determine the abundance of common carp olfactory receptor genes, explore the evolutionary fate and expression dynamics, and provide some genomic clues for the evolution of polyploid olfactory after whole-genome duplication and for future studies of teleost olfactory.

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.

Analysis of single-cell transcriptomes links enrichment of olfactory receptors with cancer cell differentiation status and prognosis

Ectopically expressed olfactory receptors (ORs) have been linked with multiple clinically-relevant physiological processes. Previously used tissue-level expression estimation largely shadowed the potential role of ORs due to their overall low expression levels. Even after the introduction of the single-cell transcriptomics, a comprehensive delineation of expression dynamics of ORs in tumors remained unexplored. Our targeted investigation into single malignant cells revealed a complex landscape of combinatorial OR expression events. We observed differentiation-dependent decline in expressed OR counts per cell as well as their expression intensities in malignant cells. Further, we constructed expression signatures based on a large spectrum of ORs and tracked their enrichment in bulk expression profiles of tumor samples from The Cancer Genome Atlas (TCGA). TCGA tumor samples stratified based on OR-centric signatures exhibited divergent survival probabilities. In summary, our comprehensive analysis positions ORs at the cross-road of tumor cell differentiation status and cancer prognosis.

A unified nomenclature for vertebrate olfactory receptors

BACKGROUND

Olfactory receptors (ORs) are G protein-coupled receptors with a crucial role in odor detection. A typical mammalian genome harbors ~ 1000 OR genes and pseudogenes; however, different gene duplication/deletion events have occurred in each species, resulting in complex orthology relationships. While the human OR nomenclature is widely accepted and based on phylogenetic classification into 18 families and further into subfamilies, for other mammals different and multiple nomenclature systems are currently in use, thus concealing important evolutionary and functional insights.

RESULTS

Here, we describe the Mutual Maximum Similarity (MMS) algorithm, a systematic classifier for assigning a human-centric nomenclature to any OR gene based on inter-species hierarchical pairwise similarities. MMS was applied to the OR repertoires of seven mammals and zebrafish. Altogether, we assigned symbols to 10,249 ORs. This nomenclature is supported by both phylogenetic and synteny analyses. The availability of a unified nomenclature provides a framework for diverse studies, where textual symbol comparison allows immediate identification of potential ortholog groups as well as species-specific expansions/deletions; for example, Or52e5 and Or52e5b represent a rat-specific duplication of OR52E5. Another example is the complete absence of OR subfamily OR6Z among primate OR symbols. In other mammals, OR6Z members are located in one genomic cluster, suggesting a large deletion in the great ape lineage. An additional 14 mammalian OR subfamilies are missing from the primate genomes. While in chimpanzee 87% of the symbols were identical to human symbols, this number decreased to ~ 50% in dog and cow and to ~ 30% in rodents, reflecting the adaptive changes of the OR gene superfamily across diverse ecological niches. Application of the proposed nomenclature to zebrafish revealed similarity to mammalian ORs that could not be detected from the current zebrafish olfactory receptor gene nomenclature.

CONCLUSIONS

We have consolidated a unified standard nomenclature system for the vertebrate OR superfamily. The new nomenclature system will be applied to cow, horse, dog and chimpanzee by the Vertebrate Gene Nomenclature Committee and its implementation is currently under consideration by other relevant species-specific nomenclature committees.

The giant pouched rat (Cricetomys ansorgei) olfactory receptor repertoire

For rodents, olfaction is essential for locating food, recognizing mates and competitors, avoiding predators, and navigating their environment. It is thought that rodents may have expanded olfactory receptor repertoires in order to specialize in olfactory behavior. Despite being the largest clade of mammals and depending on olfaction relatively little work has documented olfactory repertoires outside of conventional laboratory species. Here we report the olfactory receptor repertoire of the African giant pouched rat (Cricetomys ansorgei), a Muroid rodent distantly related to mice and rats. The African giant pouched rat is notable for its large cortex and olfactory bulbs relative to its body size compared to other sympatric rodents, which suggests anatomical elaboration of olfactory capabilities. We hypothesized that in addition to anatomical elaboration for olfaction, these pouched rats might also have an expanded olfactory receptor repertoire to enable their olfactory behavior. We examined the composition of the olfactory receptor repertoire to better understand how their sensory capabilities have evolved. We identified 1145 intact olfactory genes, and 260 additional pseudogenes within 301 subfamilies from the African giant pouched rat genome. This repertoire is similar to mice and rats in terms of size, pseudogene percentage and number of subfamilies. Analyses of olfactory receptor gene trees revealed that the pouched rat has 6 expansions in different subfamilies compared to mice, rats and squirrels. We identified 81 orthologous genes conserved among 4 rodent species and an additional 147 conserved genes within the Muroid rodents. The orthologous genes shared within Muroidea suggests that there may be a conserved Muroid-specific olfactory receptor repertoire. We also note that the description of this repertoire can serve as a complement to other studies of rodent olfaction, as the pouched rat is an outgroup within Muroidea. Thus, our data suggest that African giant pouched rats are capable of both natural and trained olfactory behaviors with a typical Muriod olfactory receptor repertoire.

Holothurians have a reduced GPCR and odorant receptor-like repertoire compared to other echinoderms

Sea cucumbers lack vision and rely on chemical sensing to reproduce and survive. However, how they recognize and respond to environmental cues remains unknown. Possible candidates are the odorant receptors (ORs), a diverse family of G protein-coupled receptors (GPCRs) involved in olfaction. The present study aimed at characterizing the chemosensory GPCRs in sea cucumbers. At least 246 distinct GPCRs, of which ca. 20% putative ORs, were found in a transcriptome assembly of putative chemosensory (tentacles, oral cavity, calcareous ring, and papillae/tegument) and reproductive (ovary and testis) tissues from Holothuria arguinensis (57 ORs) and in the Apostichopus japonicus genome (79 ORs). The sea cucumber ORs clustered with those of sea urchin and starfish into four main clades of gene expansions sharing a common ancestor and evolving under purifying selection. However, the sea cucumber ORs repertoire was the smallest among the echinoderms and the olfactory receptor signature motif LxxPxYxxxxxLxxxDxxxxxxxxP was better conserved in cluster OR-l1 which also had more members. ORs were expressed in tentacles, oral cavity, calcareous ring, and papillae/tegument, supporting their potential role in chemosensing. This study is the first comprehensive survey of chemosensory GPCRs in sea cucumbers, and provides the molecular basis to understand how they communicate.

Convergent degeneration of olfactory receptor gene repertoires in marine mammals

Background

Olfactory receptors (ORs) can bind odor molecules and play a crucial role in odor sensation. Due to the frequent gains and losses of genes during evolution, the number of OR members varies greatly among different species. However, whether the extent of gene gains/losses varies between marine mammals and related terrestrial mammals has not been clarified, and the factors that might underlie these variations are unknown.

Results

To address these questions, we identified more than 10,000 members of the OR family in 23 mammals and classified them into 830 orthologous gene groups (OGGs) and 281 singletons. Significant differences occurred in the number of OR repertoires and OGGs among different species. We found that all marine mammals had fewer OR genes than their related terrestrial lineages, with the fewest OR genes found in cetaceans, which may be closely related to olfactory degradation. ORs with more gene duplications or loss events tended to be under weaker purifying selection. The average gain and loss rates of OR genes in terrestrial mammals were higher than those of mammalian gene families, while the average gain and loss rates of OR genes in marine mammals were significantly lower and much higher than those of mammalian gene families, respectively. Additionally, we failed to detect any one-to-one orthologous genes in the focal species, suggesting that OR genes are not well conserved among marine mammals.

Conclusions

Marine mammals have experienced large numbers of OR gene losses compared with their related terrestrial lineages, which may result from the frequent birth-and-death evolution under varied functional constrains. Due to their independent degeneration, OR genes present in each lineage are not well conserved among marine mammals. Our study provides a basis for future research on the olfactory receptor function in mammals from the perspective of evolutionary trajectories.

Aquatic Adaptation and Depleted Diversity: A Deep Dive into the Genomes of the Sea Otter and Giant Otter

Despite its recent invasion into the marine realm, the sea otter (Enhydra lutris) has evolved a suite of adaptations for life in cold coastal waters, including limb modifications and dense insulating fur. This uniquely dense coat led to the near-extinction of sea otters during the 18th-20th century fur trade and an extreme population bottleneck. We used the de novo genome of the southern sea otter (E. l. nereis) to reconstruct its evolutionary history, identify genes influencing aquatic adaptation, and detect signals of population bottlenecks. We compared the genome of the southern sea otter with the tropical freshwater-living giant otter (Pteronura brasiliensis) to assess common and divergent genomic trends between otter species, and with the closely related northern sea otter (E. l. kenyoni) to uncover population-level trends. We found signals of positive selection in genes related to aquatic adaptations, particularly limb development and polygenic selection on genes related to hair follicle development. We found extensive pseudogenization of olfactory receptor genes in both the sea otter and giant otter lineages, consistent with patterns of sensory gene loss in other aquatic mammals. At the population level, the southern sea otter and the northern sea otter showed extremely low genomic diversity, signals of recent inbreeding, and demographic histories marked by population declines. These declines may predate the fur trade and appear to have resulted in an increase in putatively deleterious variants that could impact the future recovery of the sea otter.

Therapeutic potential of ectopic olfactory and taste receptors

Olfactory and taste receptors are expressed primarily in the nasal olfactory epithelium and gustatory taste bud cells, where they transmit real-time sensory signals to the brain. However, they are also expressed in multiple extra-nasal and extra-oral tissues, being implicated in diverse biological processes including sperm chemotaxis, muscle regeneration, bronchoconstriction and bronchodilatation, inflammation, appetite regulation and energy metabolism. Elucidation of the physiological roles of these ectopic receptors is revealing potential therapeutic and diagnostic applications in conditions including wounds, hair loss, asthma, obesity and cancers. This Review outlines current understanding of the diverse functions of ectopic olfactory and taste receptors and assesses their potential to be therapeutically exploited.

ClusterScan: simple and generalistic identification of genomic clusters

Summary

Studies on gene clusters proved to be an excellent source of information to understand genomes evolution and identifying specific metabolic pathways or gene families. Improvements in sequencing methods have resulted in a large increase of sequenced genomes for which cluster annotation could be performed and standardized. Currently available programs are developed to search for specific cluster types and none of them is suitable for a broad range of user-based choices. We have developed ClusterScan which allows identifying clusters of any kind of feature simply based on their genomic coordinates and user-defined categorical annotations.

Availability and implementation

The tool is written in Python, distributed under the GNU General Public License (GPL) and available on Github at http://bit.ly/ClusterScan or as Docker image at sangeslab/clusterscan: latest. It is supported through a mailing-list on http://bit.ly/ClusterScanSupport.

Supplementary information

Supplementary data are available at Bioinformatics online.

Raptor genomes reveal evolutionary signatures of predatory and nocturnal lifestyles

Background

Birds of prey (raptors) are dominant apex predators in terrestrial communities, with hawks (Accipitriformes) and falcons (Falconiformes) hunting by day and owls (Strigiformes) hunting by night.

Results

Here, we report new genomes and transcriptomes for 20 species of birds, including 16 species of birds of prey, and high-quality reference genomes for the Eurasian eagle-owl (Bubo bubo), oriental scops owl (Otus sunia), eastern buzzard (Buteo japonicus), and common kestrel (Falco tinnunculus). Our extensive genomic analysis and comparisons with non-raptor genomes identify common molecular signatures that underpin anatomical structure and sensory, muscle, circulatory, and respiratory systems related to a predatory lifestyle. Compared with diurnal birds, owls exhibit striking adaptations to the nocturnal environment, including functional trade-offs in the sensory systems, such as loss of color vision genes and selection for enhancement of nocturnal vision and other sensory systems that are convergent with other nocturnal avian orders. Additionally, we find that a suite of genes associated with vision and circadian rhythm are differentially expressed in blood tissue between nocturnal and diurnal raptors, possibly indicating adaptive expression change during the transition to nocturnality.

Conclusions

Overall, raptor genomes show genomic signatures associated with the origin and maintenance of several specialized physiological and morphological features essential to be apex predators.

Electronic supplementary material

The online version of this article (10.1186/s13059-019-1793-1) contains supplementary material, which is available to authorized users.

Characterizing novel olfactory receptors expressed in the murine renal cortex

The kidney uses specialized G protein-coupled receptors, including olfactory receptors (ORs), to act as sensors of molecules and metabolites. In the present study, we cloned and studied seven renal ORs, which we previously found to be expressed in the murine renal cortex. As most ORs are orphan receptors, our goal was to identify ligands for these ORs in the hope that this will guide future research into their functional roles. We identified novel ligands for two ORs: Olfr558 and Olfr90. For Olfr558, we confirmed activation by previously reported ligands and identified 16 additional carboxylic acids that activated this OR. The strongest activation of Olfr558 was produced by butyric, cyclobutanecarboxylic, isovaleric, 2-methylvaleric, 3-methylvaleric, 4-methylvaleric, and valeric acids. The primary in vivo source of both butyric and isovaleric acids is gut microbial metabolism. We also identified 14 novel ligands that activated Olfr90, the strongest of which were 2-methyl-4-propyl-1,3-oxathiane, 1-octen-3-ol, 2-octanol, and 3-octanol. Interestingly, 8 of these 14 ligands are of fungal origin. We also investigated the tissue distribution of these receptors and found that they are each found in a subset of "nonsensory" tissues. Finally, we examined the putative human orthologs of Olfr558 and Olfr90 and found that the human ortholog of Olfr558 (OR51E1) has a similar ligand profile, indicating that the role of this OR is likely evolutionarily conserved. In summary, we examined seven novel renal ORs and identified new ligands for Olfr558 and Olfr90, which imply that both of these receptors serve to detect metabolites produced by microorganisms.

Mammalian class I odorant receptors exhibit a conserved vestibular-binding pocket

Odorant receptors represent the largest family of mammalian G protein-coupled receptors. Phylogenetically, they are split into two classes (I and II). By analyzing the entire subclass I odorant receptors sequences, we identified two class I-specific and highly conserved motifs. These are predicted to face each other at the extra-cellular portion of the transmembrane domain, forming a vestibular site at the entrance to the orthosteric-binding cavity. Molecular dynamics simulation combined with site-directed mutagenesis and in vitro functional assays confirm the functional role of this vestibular site in ligand-driven activation. Mutations at this part of the receptor differentially affect the receptor response to four agonists. Since this vestibular site is involved in ligand recognition, it could serve ligand design that targets specifically this sub-genome of mammalian odorant receptors.

Landscape of copy number variations in Bos taurus: individual - and inter-breed variability

Background

The number of studies of Copy Number Variation in cattle has increased in recent years. This has been prompted by the increased availability of data on polymorphisms and their relationship with phenotypes. In addition, livestock species are good models for some human phenotypes. In the present study, we described the landscape of CNV driven genetic variation in a large population of 146 individuals representing 13 cattle breeds, using whole genome DNA sequence.

Results

A highly significant variation among all individuals and within each breed was observed in the number of duplications (P < 10⁻¹⁵) and in the number of deletions (P < 10⁻¹⁵). We also observed significant differences between breeds for duplication (P = 0.01932) and deletion (P = 0.01006) counts. The same variation CNV length - inter-individual and inter-breed differences were significant for duplications (P < 10⁻¹⁵) and deletions (P < 10⁻¹⁵). Moreover, breed-specific variants were identified, with the largest proportion of breed-specific duplications (9.57%) found for Fleckvieh and breed-specific deletions found for Brown Swiss (5.00%). Such breed-specific CNVs were predominantly located in intragenic regions, however in Simmental, one deletion present in five individuals was found in the coding sequence of a novel gene ENSBTAG00000000688 on chromosome 18. In Brown Swiss, Norwegian Red and Simmental breed-specific deletions were located within KIT and MC1R genes, which are responsible for a coat colour. The functional annotation of coding regions underlying the breed-specific CNVs showed that in Norwegian Red, Guernsey, and Simmental significantly under- and overrepresented GO terms were related to chemical stimulus involved in sensory perception of smell and the KEGG pathways for olfactory transduction. In addition, specifically for the Norwegian Red breed, the dopaminergic synapse KEGG pathway was significantly enriched within deleted parts of the genome.

Conclusions

The CNV landscape in Bos taurus genome revealed by this study was highly complex, with inter-breed differences, but also a significant variation within breeds. The former, may explain some of the phenotypic differences among analysed breeds, and the latter contributes to within-breed variation available for selection.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4815-6) contains supplementary material, which is available to authorized users.

Olfactory receptors are expressed in pancreatic β-cells and promote glucose-stimulated insulin secretion

Olfactory receptors (ORs) mediate olfactory chemo-sensation in OR neurons. Herein, we have demonstrated that the OR chemo-sensing machinery functions in pancreatic β-cells and modulates insulin secretion. First, we found several OR isoforms, including OLFR15 and OLFR821, to be expressed in pancreatic islets and a β-cell line, MIN6. Immunostaining revealed OLFR15 and OLFR821 to be uniformly expressed in pancreatic β-cells. In addition, mRNAs of Olfr15 and Olfr821 were detected in single MIN6 cells. These results indicate that multiple ORs are simultaneously expressed in individual β-cells. Octanoic acid, which is a medium-chain fatty acid contained in food and reportedly interacts with OLFR15, potentiated glucose-stimulated insulin secretion (GSIS), thereby improving glucose tolerance in vivo. GSIS potentiation by octanoic acid was confirmed in isolated pancreatic islets and MIN6 cells and was blocked by OLFR15 knockdown. While Gα_olf expression was not detectable in β-cells, experiments using inhibitors and siRNA revealed that the pathway dependent on phospholipase C-inositol triphosphate, rather than cAMP-protein kinase A, mediates GSIS potentiation via OLFR15. These findings suggest that the OR system in pancreatic β-cells has a chemo-sensor function allowing recognition of environmental substances obtained from food, and potentiates insulin secretion in a cell-autonomous manner, thereby modulating systemic glucose metabolism.

Evolution of Trichocyte Keratin Associated Proteins

The major components of hair are keratins and keratin associated proteins (KRTAPs). KRTAPs form the interfilamentous matrix between intermediate filament bundles through extensive disulfide bond cross-linking with the numerous cysteine residues in hair keratins. A variable number of approximately100-180 genes compose the KRTAP gene family in mammals. KRTAP gene family members present a typical pattern of concerted evolution, and its evolutionary features are consistent with the evolution of mammalian hair. KRATP genes might be more important in determining the structure of cashmere fibers in domestic mammals like sheep and goats. KRTAP gene variants thus should provide information for improved wool by sheep and goat breeding.

Annotated Draft Genome Assemblies for the Northern Bobwhite (Colinus virginianus) and the Scaled Quail (Callipepla squamata) Reveal Disparate Estimates of Modern Genome Diversity and Historic Effective Population Size

Northern bobwhite (Colinus virginianus; hereafter bobwhite) and scaled quail (Callipepla squamata) populations have suffered precipitous declines across most of their US ranges. Illumina-based first- (v1.0) and second- (v2.0) generation draft genome assemblies for the scaled quail and the bobwhite produced N50 scaffold sizes of 1.035 and 2.042 Mb, thereby producing a 45-fold improvement in contiguity over the existing bobwhite assembly, and ≥90% of the assembled genomes were captured within 1313 and 8990 scaffolds, respectively. The scaled quail assembly (v1.0 = 1.045 Gb) was ∼20% smaller than the bobwhite (v2.0 = 1.254 Gb), which was supported by kmer-based estimates of genome size. Nevertheless, estimates of GC content (41.72%; 42.66%), genome-wide repetitive content (10.40%; 10.43%), and MAKER-predicted protein coding genes (17,131; 17,165) were similar for the scaled quail (v1.0) and bobwhite (v2.0) assemblies, respectively. BUSCO analyses utilizing 3023 single-copy orthologs revealed a high level of assembly completeness for the scaled quail (v1.0; 84.8%) and the bobwhite (v2.0; 82.5%), as verified by comparison with well-established avian genomes. We also detected 273 putative segmental duplications in the scaled quail genome (v1.0), and 711 in the bobwhite genome (v2.0), including some that were shared among both species. Autosomal variant prediction revealed ∼2.48 and 4.17 heterozygous variants per kilobase within the scaled quail (v1.0) and bobwhite (v2.0) genomes, respectively, and estimates of historic effective population size were uniformly higher for the bobwhite across all time points in a coalescent model. However, large-scale declines were predicted for both species beginning ∼15-20 KYA.

β2-microglobulin gene duplication in cetartiodactyla remains intact only in pigs and possibly confers selective advantage to the species

Several β2-microglobulin (B2M) -bound protein complexes undertake key roles in various immune system pathways, including the neonatal Fc receptor (FcRn), cluster of differentiation 1 (CD1) protein, non-classical major histocompatibility complex (MHC), and well-known MHC class I molecules. Therefore, the duplication of B2M may lead to an increase in the biological competence of organisms to the environment. Based on the pig genome assembly SSC10.2, a segmental duplication of ~45.5 kb, encoding the entire B2M protein, was identified in pig chromosome 1. Through experimental validation, we confirmed the functional duplication of the B2M gene with a completely identical coding sequence between two copies in pigs. Considering the importance of B2M in the immune system, we performed the phylogenetic analysis of B2M duplication in ten mammalian species, confirming the presence of B2M duplication in cetartioldactyls, like cattle, sheep, goats, pigs and whales, but non-cetartiodactyl species, like mice, cats, dogs, horses, and humans. The density of long interspersed nuclear element (LINE) at the edges of duplicated blocks (39 to 66%) was found to be 2 to 3-fold higher than the average (20.12%) of the pig genome, suggesting its role in the duplication event. The B2M mRNA expression level in pigs was 12.71 and 7.57 times (2-ΔΔCt values) higher than humans and mice, respectively. However, we were unable to experimentally demonstrate the difference in the level of B2M protein because species specific anti-B2M antibodies are not available. We reported, for the first time, the functional duplication of the B2M gene in animals. The identification of partially remaining duplicated B2M sequences in the genomes of only cetartiodactyls indicates that the event was lineage specific. B2M duplication could be beneficial to the immune system of pigs by increasing the availability of MHC class I light chain protein, B2M, to complex with the proteins encoded by the relatively large number of MHC class I heavy chain genes in pigs. Further studies are necessary to address the biological meaning of increased expression of B2M.

Evolutionary History of Chemosensory-Related Gene Families across the Arthropoda

Chemosensory-related gene (CRG) families have been studied extensively in insects, but their evolutionary history across the Arthropoda had remained relatively unexplored. Here, we address current hypotheses and prior conclusions on CRG family evolution using a more comprehensive data set. In particular, odorant receptors were hypothesized to have proliferated during terrestrial colonization by insects (hexapods), but their association with other pancrustacean clades and with independent terrestrial colonizations in other arthropod subphyla have been unclear. We also examine hypotheses on which arthropod CRG family is most ancient. Thus, we reconstructed phylogenies of CRGs, including those from new arthropod genomes and transcriptomes, and mapped CRG gains and losses across arthropod lineages. Our analysis was strengthened by including crustaceans, especially copepods, which reside outside the hexapod/branchiopod clade within the subphylum Pancrustacea. We generated the first high-resolution genome sequence of the copepod Eurytemora affinis and annotated its CRGs. We found odorant receptors and odorant binding proteins present only in hexapods (insects) and absent from all other arthropod lineages, indicating that they are not universal adaptations to land. Gustatory receptors likely represent the oldest chemosensory receptors among CRGs, dating back to the Placozoa. We also clarified and confirmed the evolutionary history of antennal ionotropic receptors across the Arthropoda. All antennal ionotropic receptors in E. affinis were expressed more highly in males than in females, suggestive of an association with male mate-recognition behavior. This study is the most comprehensive comparative analysis to date of CRG family evolution across the largest and most speciose metazoan phylum Arthropoda.