Life history and chromosome organization determine chemoreceptor gene expression in rattlesnakes

Predatory species who hunt for their prey rely on a suite of integrated characters, including sensory traits that are also used for nonpredatory behaviors. Linking the evolution of sensory traits to specific selection pressures therefore requires a deep understanding of the underlying genetics and molecular mechanisms producing these complex phenotypes. However, this relationship remains poorly understood for complex sensory systems that consist of proteins encoded by large gene families. The chemosensory repertoire of rattlesnakes includes hundreds of type-2 vomeronasal receptors and olfactory receptors, representing the two largest gene families found in the genome. To investigate the biological importance of this chemoreceptor diversity, we assessed gene expression in the eastern diamondback rattlesnake (Crotalus adamanteus) and identified sex- and age-biased genes. We found type-2 vomeronasal receptor expression in the vomeronasal epithelium was limited to juvenile snakes, suggesting the sensory programming of this tissue may be correlated with early life development. In the olfactory epithelium, we found subtle expression biases that were more indicative of life history rather than development. We also found transcriptional evidence for dosage compensation of sex-linked genes and trait integration in the expression of transcription factors. We overlay our molecular characterizations in Crotalus adamanteus onto updated olfactory receptor and type-2 vomeronasal receptor phylogenies, providing a genetic road map for future research on these receptors. Finally, we investigated the deeper macroevolutionary context of the most highly expressed type-2 vomeronasal receptor gene spanning the rise of tetrapods and estimated the strength of positive selection for individual amino acid residues in the predicted protein structure. We hypothesize that this gene may have evolved as a conserved signaling subunit to ensure consistent G-protein coupled receptor functionality, potentially relaxing signaling constraints on other type-2 vomeronasal receptor paralogs and promoting ligand binding specificity.

The Importance of Olfaction for Mixed Paternity in Birds

Olfaction can aid individuals in finding genetically compatible mates in many animals, while high levels of mixed paternity may result from a limited ability to evaluate their mate's genetic profile against their own before mating. To test this suggestion and explore if olfaction may indeed influence mating patterns in birds, we combined published measures of olfactory ability with data on genetic mating pattern in the same species, across a phylogenetically broad range of species. We used three measures of olfaction: (1) olfactory bulb diameter, (2) olfactory bulb volume and (3) number of olfactory receptor genes (148, 134 and 48 species, respectively). These measures were then matched to species-specific estimates of mating pattern, measured as percentage of broods with mixed paternity (> 1 male siring offspring in the same brood). Limited overlaps between the datasets resulted in 30 matched species for olfactory bulb diameter, 31 for olfactory bulb volume and 15 for olfactory receptor genes. Controlling for brain size (telencephalon), we then correlated olfaction to mating pattern, and found that the bigger the relative olfactory bulb diameter, the lower the proportion of mixed paternity. In contrast, there was no significant correlation between olfactory bulb volume or number of receptor genes and paternity. This study thus indicates that mating patterns in birds may be influenced by olfactory ability, measured as olfactory bulb diameter. Next, we suggest expanding the datasets by collecting olfactory-focused measures, targeting species for which paternity measures already exist, to allow a full phylogenetic analysis.

The scent of habitat shift: Olfactory receptor evolution is associated with environmental transitions in turtles

The transition between aquatic and terrestrial habitats leads to extreme structural changes in sensorial systems. Olfactory receptors (OR) are involved in the detection of odorant molecules both in water and on land. Therefore, ORs are affected by evolutionary habitat transitions experienced by organisms. In this study, we used turtles, a group of vertebrates which inhabit many distinct environments, to explore whether functional olfactory gene receptor repertoires are correlated to habitat. We found that the proportion of class I vs class II functional olfactory receptor genes (used for waterborne odorant detection and volatile odorant detection, respectively) was closely linked to habitat. Fully terrestrial turtles had the largest proportion of class II functional receptor genes while marine turtles had a larger proportion of class I receptor genes. Freshwater turtles had more balanced numbers of class I and class II functional receptor genes, but showed a gradient of OR type proportions likely reflecting species-specific amphibious preferences. Interestingly, freshwater turtles had by far the largest number of functional OR genes compared to those in other habitats, challenging the hypothesis that secondary adaptions to water may have reduced OR repertoires in amniotes. Our study provides novel results which shed new light on the relationship between chemical communication and habitat.

Evolution of olfactory receptor superfamily in bats based on high throughput molecular modelling

The origin of flight and laryngeal echolocation in bats is likely to have been accompanied by evolutionary changes in other aspects of their sensory biology. Of all sensory modalities in bats, olfaction is perhaps the least well understood. Olfactory receptors (ORs) function in recognizing odour molecules, with crucial roles in evaluating food, as well as in processing social information. Here we compare OR repertoire sizes across taxa and apply a new pipeline that integrates comparative genome data with protein structure modelling and then we employ molecular docking techniques with small molecules to analyse OR functionality based on binding energies. Our results suggest a sharp contraction in odorant recognition of the functional OR repertoire during the origin of bats, consistent with a reduced dependence on olfaction. We also compared bat lineages with contrasting different ecological characteristics and found evidence of differences in OR gene expansion and contraction, and in the composition of ORs with different tuning breadths. The strongest binding energies of ORs in non-echolocating fruit-eating bats were seen to correspond to ester odorants, although we did not detect a quantitative advantage of functional OR repertoires in these bats compared with echolocating insectivorous species. Overall, our findings based on molecular modelling and computational docking suggest that bats have undergone olfactory evolution linked to dietary adaptation. Our results from extant and ancestral bats help to lay the groundwork for targeted experimental functional tests in the future.

A family of olfactory receptors uniquely expanded in marsupial and monotreme genomes are expressed by a T cell subset also unique to marsupials and monotremes

Olfactory receptors (OR), expressed on olfactory neurons, mediate the sense of smell. Recently, OR have also been shown to be expressed in non-olfactory tissues, including cells of the immune system. An analysis of single-cell transcriptomes of splenocytes of the grey short-tailed opossum (Monodelphis domestica) found OR are expressed on a subset of T cells, the γμ T cells, that are unique to marsupials and monotremes. A majority of opossum γμ T cells transcriptomes contain OR family 14 transcripts, specifically, from the OR14C subfamily. Amongst the mammals, the OR14 gene family is expanded in the genomes of marsupials and monotremes, and rarer or absent in placental mammals. In summary, here we demonstrate the intriguing correlation that a family of OR genes, abundant in the genomes of marsupials and monotremes, are ectopically expressed in a particular subset of T cells unique to the marsupials and monotremes.

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.

Phylogeny, Ecology, and Gene Families Covariation Shaped the Olfactory Subgenome of Rodents

Olfactory receptor (OR) genes represent the largest multigenic family in mammalian genomes and encode proteins that bind environmental odorant molecules. The OR repertoire is extremely variable among species and is subject to many gene duplications and losses, which have been linked to ecological adaptations in mammals. Although they have been studied on a broad taxonomic scale (i.e., placental), finer sampling has rarely been explored in order to better capture the mechanisms that drove the evolution of the OR repertoire. Among placental mammals, rodents are well-suited for this task, as they exhibit diverse life history traits, and genomic data are available for most major families and a diverse array of lifestyles. In this study, 53 rodent published genomes were mined for their OR subgenomes. We retrieved more than 85,000 functional and pseudogene OR sequences that were subsequently classified into phylogenetic clusters. Copy number variation among rodents is similar to that of other mammals. Using our OR counts along with comparative phylogenetic approaches, we demonstrated that ecological niches such as diet, period of activity, and a fossorial lifestyle strongly impacted the proportion of OR pseudogenes. Within the OR subgenome, phylogenetic inertia was the main factor explaining the relative variations of the 13 OR gene families. However, a striking exception was a convergent 10-fold expansion of the OR family 14 among the phylogenetically divergent subterranean mole-rat lineages belonging to Bathyergidae and Spalacidae families. This study illustrates how the diversity of the OR repertoire has evolved among rodents, both shaped by selective forces stemming from species life history traits and neutral evolution along the rodent phylogeny.

Parallel evolution of fish bi-modal breathing and expansion of olfactory receptor (OR) genes: toward a universal ORs nomenclature

Olfactory receptors (ORs) play a key role in the prime sensorial perception, being highly relevant for intra/interspecific interactions. ORs are a subgroup of G-protein coupled receptors that exhibit highly complex subgenomes in vertebrates. However, OR repertoires remain poorly studied in fish lineages, precluding finely retracing their origin, evolution, and diversification, especially in the most basal groups. Here, we conduct an exhaustive gene screening upon 43 high-quality fish genomes exhibiting varied gene repertoires (2-583 genes). While the early vertebrates performed gas exchange through gills, we hypothesize that the emergence of new breathing structures (swim bladder and paired lungs) in early osteichthyans may be associated with expansions in the ORs gene families sensitive to airborne molecules. Additionally, we verify that the OR repertoire of moderns actinopterygians has not increased as expected following a whole genome duplication, likely due to regulatory mechanisms compensating the gene load excess. Finally, we identify 25 distinct OR families, allowing us to propose an updated universal nomenclature for the fish ORs.

Phylogenetically under-dispersed gut microbiomes are not correlated with host genomic heterozygosity in a genetically diverse reptile community

While key elements of fitness in vertebrate animals are impacted by their microbiomes, the host genetic characteristics that factor into microbiome composition are not fully understood. Here, we correlate host genomic heterozygosity and gut microbiome phylogenetic diversity across a community of reptiles in southwestern New Mexico to test hypotheses about the behaviour of host genes that drive microbiome assembly. We find that microbiome communities are phylogenetically under-dispersed relative to random expectations, and that host heterozygosity is not correlated with microbiome diversity. Our analyses reinforce results from functional genomic work that identify conserved host immune and nonimmune genes as key players in microbiome assembly, rather than gene families that rely on heterozygosity for their function.

Amplification of olfactory signals by Anoctamin 9 is important for mammalian olfaction

Sensing smells of foods, prey, or predators determines animal survival. Olfactory sensory neurons in the olfactory epithelium (OE) detect odorants, where cAMP and Ca2+ play a significant role in transducing odorant inputs to electrical activity. Here we show Anoctamin 9, a cation channel activated by cAMP/PKA pathway, is expressed in the OE and amplifies olfactory signals. Ano9-deficient mice had reduced olfactory behavioral sensitivity, electro-olfactogram signals, and neural activity in the olfactory bulb. In line with the difference in olfaction between birds and other vertebrates, chick ANO9 failed to respond to odorants, whereas chick CNGA2, a major transduction channel, showed greater responses to cAMP. Thus, we concluded that the signal amplification by ANO9 is important for mammalian olfactory transduction.

Genomic insight into the nocturnal adaptation of the black-crowned night heron (Nycticorax nycticorax)

BACKGROUND

The black-crowned night heron (Nycticorax nycticorax) is an ardeid bird successfully adapted to the nocturnal environment. Previous studies had indicated that the eyes of the night herons have evolved several specialized morphological traits favoring nocturnal vision. However, the molecular mechanisms of the nocturnal vision adaptation of night herons remained inattentions. In this study, the whole genome of N. nycticorax was sequenced and comparative analyses were performed on the vision-related and olfactory receptor (OR) genes to understand the molecular mechanisms of the visual and olfactory adaptation of night herons.

RESULTS

The results indicated that a number of vision genes were under positive or relaxed selection in N. nycticorax, whereas a number of other vision genes were under relaxed or intensified selection in the boat-billed heron (Cochlearius cochlearius), which suggested that the two species adapt to nocturnality with different genetic mechanisms. The different selections acting on vision genes are probably associated with the enlargement of eye size and the enhancement of visual sensitivity in night herons. The analyses on olfactory receptor (OR) genes indicated that the total number of OR genes in the genomes of N. nycticorax and C. cochlearius were about half those in the little egret (Egretta garzetta), whereas the diversity of their OR genes was not remarkably different. Additionally, the number of expressed OR genes in the transcriptomes of N. nycticorax was also fewer than that in E. garzetta. These results suggest a reduced olfactory capability in night herons compared with E. garzetta.

CONCLUSIONS

Our results provided evidence that several vision genes of the night herons were subjected to different natural selections, which can contribute to a better understanding of the genetic mechanisms of visual adaptions of the night heron. In addition, the finding of the reduced number of total and expressed OR genes in night herons may reflect a trade-off between olfaction and vision.

Comparative Genomics Provides Insights into Adaptive Evolution in Tactile-Foraging Birds

Tactile-foraging birds have evolved an enlarged principal sensory nucleus (PrV) but smaller brain regions related to the visual system, which reflects the difference in sensory dependence. The “trade-off” may exist between different senses in tactile foragers, as well as between corresponding sensory-processing areas in the brain. We explored the mechanism underlying the adaptive evolution of sensory systems in three tactile foragers (kiwi, mallard, and crested ibis). The results showed that olfaction-related genes in kiwi and mallard and hearing-related genes in crested ibis were expanded, indicating they may also have sensitive olfaction or hearing, respectively. However, some genes required for visual development were positively selected or had convergent amino acid substitutions in all three tactile branches, and it seems to show the possibility of visual degradation. In addition, we may provide a new visual-degradation candidate gene PDLIM1 who suffered dense convergent amino acid substitutions within the ZM domain. At last, two genes responsible for regulating the proliferation and differentiation of neuronal progenitor cells may play roles in determining the relative sizes of sensory areas in brain. This exploration offers insight into the relationship between specialized tactile-forging behavior and the evolution of sensory abilities and brain structures.

A high-quality assembly reveals genomic characteristics, phylogenetic status, and causal genes for leucism plumage of Indian peafowl

BACKGROUND

The dazzling phenotypic characteristics of male Indian peafowl (Pavo cristatus) are attractive both to the female of the species and to humans. However, little is known about the evolution of the phenotype and phylogeny of these birds at the whole-genome level. So far, there are no reports regarding the genetic mechanism of the formation of leucism plumage in this variant of Indian peafowl.

RESULTS

A draft genome of Indian peafowl was assembled, with a genome size of 1.05 Gb (the sequencing depth is 362×), and contig and scaffold N50 were up to 6.2 and 11.4 Mb, respectively. Compared with other birds, Indian peafowl showed changes in terms of metabolism, immunity, and skeletal and feather development, which provided a novel insight into the phenotypic evolution of peafowl, such as the large body size and feather morphologies. Moreover, we determined that the phylogeny of Indian peafowl was more closely linked to turkey than chicken. Specifically, we first identified that PMEL was a potential causal gene leading to the formation of the leucism plumage variant in Indian peafowl.

CONCLUSIONS

This study provides an Indian peafowl genome of high quality, as well as a novel understanding of phenotypic evolution and phylogeny of Indian peafowl. These results provide a valuable reference for the study of avian genome evolution. Furthermore, the discovery of the genetic mechanism for the development of leucism plumage is both a breakthrough in the exploration of peafowl plumage and also offers clues and directions for further investigations of the avian plumage coloration and artificial breeding in peafowl.

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.

Olfactory receptor subgenome and expression in a highly olfactory procellariiform seabird

Procellariiform seabirds rely on their sense of smell for foraging and homing. Both genomes and transcriptomes yield important clues about how olfactory receptor (OR) subgenomes are shaped by natural and sexual selection, yet no transcriptomes have been made of any olfactory epithelium of any bird species thus far. Here we assembled a high-quality genome and nasal epithelium transcriptome of the Leach's storm-petrel (Oceanodroma leucorhoa) to extensively characterize their OR repertoire. Using a depth-of-coverage-assisted counting method, we estimated over 160 intact OR genes (∼500 including OR fragments). This method reveals the highest number of intact OR genes and the lowest proportion of pseudogenes compared to other waterbirds studied, and suggests that rates of OR gene duplication vary between major clades of birds, with particularly high rates in passerines. OR expression patterns reveal two OR genes (OR6-6 and OR5-11) highly expressed in adults, and four OR genes (OR14-14, OR14-12, OR10-2, and OR14-9) differentially expressed between age classes of storm-petrels. All four genes differentially expressed between age classes were more highly expressed in chicks compared to adults, suggesting that ORs genes may exhibit ontogenetic specializations. Three highly differentially expressed OR genes also had high copy number ratios, suggesting that expression variation may be linked to copy number in the genome. We provide better estimates of OR gene number by using a copy number-assisted counting method, and document ontogenetic changes in OR gene expression that may be linked to olfactory specialization. These results provide valuable insight into the expression, development, and macroevolution of olfaction in seabirds.

Transposable elements mark a repeat-rich region associated with migratory phenotypes of willow warblers (Phylloscopus trochilus)

The genetic basis of bird migration has been the focus of several studies. Two willow warbler subspecies (Phylloscopus trochilus trochilus and Phylloscopus trochilus acredula) follow different migratory routes to wintering grounds in Africa. Their breeding populations overlap in contact areas or "migratory divides" located in central Scandinavia and in eastern Poland. Earlier analyses demonstrated that the genetic differences between these two migratory phenotypes are few and cluster on chromosomes 1 and 5. In addition, an amplified fragment length polymorphism-derived biallelic marker (known as WW2) presents steep clines across both migratory divides but failed to be mapped in the genome. Here, we characterize the WW2 marker and describe its two variants (WW2 ancestral and WW2 derived) as portions of long terminal repeat retrotransposons originating from an ancient infection by an endogenous retrovirus. We used quantitative polymerase chain reaction techniques to quantify copy numbers of the WW2 derived variant in the two subspecies and their hybrids. This, together with genome analyses revealed that WW2 derived variants are much more abundant in P. t. acredula and appear embedded in a large repeat-rich region (>12 Mbp), not associated with the divergent regions of chromosomes 1 or 5. However, it might interact with genetic elements controlling migration direction. Testing this hypothesis further will require knowing the exact location of this region, such as by obtaining more complete genome assemblies preferably in combination with techniques like fluorescence in situ hybridization applied to a willow warbler karyotype, and finally to investigate the copy number of this marker in hybrids with known migratory tracks.

What Have We Learned from the First 500 Avian Genomes?

The increased capacity of DNA sequencing has significantly advanced our understanding of the phylogeny of birds and the proximate and ultimate mechanisms molding their genomic diversity. In less than a decade, the number of available avian reference genomes has increased to over 500—approximately 5% of bird diversity—placing birds in a privileged position to advance the fields of phylogenomics and comparative, functional, and population genomics. Whole-genome sequence data, as well as indels and rare genomic changes, are further resolving the avian tree of life. The accumulation of bird genomes, increasingly with long-read sequence data, greatly improves the resolution of genomic features such as germline-restricted chromosomes and the W chromosome, and is facilitating the comparative integration of genotypes and phenotypes. Community-based initiatives such as the Bird 10,000 Genomes Project and Vertebrate Genome Project are playing a fundamental role in amplifying and coalescing a vibrant international program in avian comparative genomics.

Olfactory receptor gene evolution is unusually rapid across Tetrapoda and outpaces chemosensory phenotypic change

Chemosensation is the most ubiquitous sense in animals, enacted by the products of complex gene families that detect environmental chemical cues and larger-scale sensory structures that process these cues. While there is a general conception that olfactory receptor (OR) genes evolve rapidly, the universality of this phenomenon across vertebrates, and its magnitude, are unclear. The supposed correlation between molecular rates of chemosensory evolution and phenotypic diversity of chemosensory systems is largely untested. We combine comparative genomics and sensory morphology to test whether OR genes and olfactory phenotypic traits evolve at faster rates than other genes or traits. Using published genomes, we identified ORs in 21 tetrapods, including amphibians, reptiles, birds, and mammals and compared their rates of evolution to those of orthologous non-OR protein-coding genes. We found that, for all clades investigated, most OR genes evolve nearly an order of magnitude faster than other protein-coding genes, with many OR genes showing signatures of diversifying selection across nearly all taxa in this study. This rapid rate of evolution suggests that chemoreceptor genes are in "evolutionary overdrive," perhaps evolving in response to the ever-changing chemical space of the environment. To obtain complementary morphological data, we stained whole fixed specimens with iodine, µCT-scanned the specimens, and digitally segmented chemosensory and nonchemosensory brain regions. We then estimated phenotypic variation within traits and among tetrapods. While we found considerable variation in chemosensory structures, they were no more diverse than nonchemosensory regions. We suggest chemoreceptor genes evolve quickly in reflection of an ever-changing chemical space, whereas chemosensory phenotypes and processing regions are more conserved because they use a standardized or constrained architecture to receive and process a range of chemical cues.

Adaptive selection signatures in river buffalo with emphasis on immune and major histocompatibility complex genes

River buffalo is an agriculturally important species with many traits, such as disease tolerance, which promote its use worldwide. Highly contiguous genome assemblies of the river buffalo, goat, pig, human and two cattle subspecies were aligned to study gene gains and losses and signs of positive selection. The gene families that have changed significantly in river buffalo since divergence from cattle play important roles in protein degradation, the olfactory receptor system, detoxification and the immune system. We used the branch site model in PAML to analyse single-copy orthologs to identify positively selected genes that may be involved in skin differentiation, mammary development and bone formation in the river buffalo branch. The high contiguity of the genomes enabled evaluation of differences among species in the major histocompatibility complex. We identified a Babesia-like L1 LINE insertion in the DRB1-like gene in the river buffalo and discuss the implication of this finding.

Highly contiguous genomes improve the understanding of avian olfactory receptor repertoires

Third generation (long read-based) sequencing technologies are reshaping our understanding of genome structure and function. One of the most persistent challenges in genome biology has been confidently reconstructing radiations of complex gene families. Olfactory receptors (ORs) represent just such a gene family with upwards of 1000s of receptors in some mammalian taxa. Whereas in birds olfaction was historically an overlooked sensory modality, new studies have revealed an important role for smell. Chromosome-level assemblies for birds allow a new opportunity to characterize patterns of OR diversity among major bird lineages. Previous studies of short read (second-generation) genome assemblies have associated OR gene family size with avian ecology, but such conclusions could be premature if new assembly methods reshape our understanding of avian OR evolution. Here we provide a fundamental characterization of OR repertoires in five recent genome assemblies, including the most recent assembly of golden-collared manakin (Manacus vitellinus). We find that short read-based assemblies systematically undercount the avian-specific gamma-c OR subfamily, a subfamily that comprises over 65% of avian OR diversity. Therefore, in contrast to previous studies we find a high diversity of gamma-c ORs across the avian tree of life. Building on these findings, ongoing sequencing efforts and improved genome assemblies will clarify the relationship between OR diversity and avian ecology.

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.

Evolutionary dynamics of the OR gene repertoire in teleost fishes: evidence of an association with changes in olfactory epithelium shape

Teleost fishes perceive their environment through a range of sensory modalities, among which olfaction often plays an important role. Richness of the olfactory repertoire depends on the diversity of receptors coded by homologous genes classified into four families: OR, TAAR, VR1, and VR2. Herein, we focus on the OR gene repertoire. While independent large contractions of the OR gene repertoire associated with ecological transitions have been found in mammals, little is known about the diversity of the OR gene repertoire and its evolution in teleost fishes, a group that includes more than 34,000 living species. We analyzed genomes of 163 species representing diversity in this large group. We found a large range of variation in the number of functional OR genes, from 15 in the Broad-nose Pipefish Syngnathus typhle and the Ocean Sunfish Mola mola, to 429 in the Zig-zag Eel Mastacembelus armatus. The number of OR genes was higher in species when a multilamellar olfactory rosette was present. Moreover, the number of lamellae was correlated with the richness of the OR gene repertoire. While a slow and balanced birth-and-death process generally drives the evolution of the OR gene repertoire, we inferred several episodes of high rates of gene loss, sometimes followed by large gains in the number of OR genes. These gains coincide with morphological changes of the olfactory organ and suggest a strong functional association between changes in the morphology and the evolution of the OR gene repertoire.

The Chemosensory Repertoire of the Eastern Diamondback Rattlesnake (Crotalus adamanteus) Reveals Complementary Genetics of Olfactory and Vomeronasal-Type Receptors

Pitviper sensory perception incorporates diverse stimuli through the integration of trichromatic color vision, bifocal heat-sensing, and dual-system chemoperception. Chemoperception, or olfaction, is mediated by chemoreceptors in the olfactory bulb and the vomeronasal organ, but the true genomic complexity of the gene families and their relative contributions is unknown. A full genomic accounting of pitviper chemoperception directly complements our current understanding of their venoms by generating a more complete polyphenic representation of their predatory arsenal. To characterize the genetic repertoire of pitviper chemoperception, we analyzed a full-genome assembly for Crotalus adamanteus, the eastern diamondback rattlesnake. We identified hundreds of genes encoding both olfactory receptors (ORs; 362 full-length genes) and type-2 vomeronasal receptors (V2Rs; 430 full-length genes). Many chemoreceptor genes are organized into large tandem repeat arrays. Comparative analysis of V2R orthologs across squamates demonstrates how gene array expansion and contraction underlies the evolution of the chemoreceptor repertoire, which likely reflects shifts in life history traits. Chromosomal assignments of chemosensory genes identified sex chromosome specific chemoreceptor genes, providing gene candidates underlying observed sex-specific chemosensory-based behaviors. We detected widespread episodic evolution in the extracellular, ligand-binding domains of both ORs and V2Rs, suggesting the diversification of chemoreceptors is driven by transient periods of positive selection. We provide a robust genetic framework for studying pitviper chemosensory ecology and evolution.

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.

Putative Antimicrobial Peptides of the Posterior Salivary Glands from the Cephalopod Octopus vulgaris Revealed by Exploring a Composite Protein Database

Cephalopods, successful predators, can use a mixture of substances to subdue their prey, becoming interesting sources of bioactive compounds. In addition to neurotoxins and enzymes, the presence of antimicrobial compounds has been reported. Recently, the transcriptome and the whole proteome of the Octopus vulgaris salivary apparatus were released, but the role of some compounds—e.g., histones, antimicrobial peptides (AMPs), and toxins—remains unclear. Herein, we profiled the proteome of the posterior salivary glands (PSGs) of O. vulgaris using two sample preparation protocols combined with a shotgun-proteomics approach. Protein identification was performed against a composite database comprising data from the UniProtKB, all transcriptomes available from the cephalopods’ PSGs, and a comprehensive non-redundant AMPs database. Out of the 10,075 proteins clustered in 1868 protein groups, 90 clusters corresponded to venom protein toxin families. Additionally, we detected putative AMPs clustered with histones previously found as abundant proteins in the saliva of O. vulgaris. Some of these histones, such as H2A and H2B, are involved in systemic inflammatory responses and their antimicrobial effects have been demonstrated. These results not only confirm the production of enzymes and toxins by the O. vulgaris PSGs but also suggest their involvement in the first line of defense against microbes.

Functional evolution of vertebrate sensory receptors

Sensory receptors enable animals to perceive their external world, and functional properties of receptors evolve to detect the specific cues relevant for an organism's survival. Changes in sensory receptor function or tuning can directly impact an organism's behavior. Functional tests of receptors from multiple species and the generation of chimeric receptors between orthologs with different properties allow for the dissection of the molecular basis of receptor function and identification of the key residues that impart functional changes in different species. Knowledge of these functionally important sites facilitates investigation into questions regarding the role of epistasis and the extent of convergence, as well as the timing of sensory shifts relative to other phenotypic changes. However, as receptors can also play roles in non-sensory tissues, and receptor responses can be modulated by numerous other factors including varying expression levels, alternative splicing, and morphological features of the sensory cell, behavioral validation can be instrumental in confirming that responses observed in heterologous systems play a sensory role. Expression profiling of sensory cells and comparative genomics approaches can shed light on cell-type specific modifications and identify other proteins that may affect receptor function and can provide insight into the correlated evolution of complex suites of traits. Here we review the evolutionary history and diversity of functional responses of the major classes of sensory receptors in vertebrates, including opsins, chemosensory receptors, and ion channels involved in temperature-sensing, mechanosensation and electroreception.

Measuring phylogenetic signal between categorical traits and phylogenies

MOTIVATION

Determining whether a trait and phylogeny share some degree of phylogenetic signal is a flagship goal in evolutionary biology. Signatures of phylogenetic signal can assist the resolution of a broad range of evolutionary questions regarding the tempo and mode of phenotypic evolution. However, despite the considerable number of strategies to measure it, few and limited approaches exist for categorical traits. Here, we used the concept of Shannon entropy and propose the δ statistic for evaluating the degree of phylogenetic signal between a phylogeny and categorical traits.

RESULTS

We validated δ as a measure of phylogenetic signal: the higher the δ-value the higher the degree of phylogenetic signal between a given tree and a trait. Based on simulated data we proposed a threshold-based classification test to pinpoint cases of phylogenetic signal. The assessment of the test's specificity and sensitivity suggested that the δ approach should only be applied to 20 or more species. We have further tested the performance of δ in scenarios of branch length and topology uncertainty, unbiased and biased trait evolution and trait saturation. Our results showed that δ may be applied in a wide range of phylogenetic contexts. Finally, we investigated our method in 14 360 mammalian gene trees and found that olfactory receptor genes are significantly associated with the mammalian activity patterns, a result that is congruent with expectations and experiments from the literature. Our application shows that δ can successfully detect molecular signatures of phenotypic evolution. We conclude that δ represents a useful measure of phylogenetic signal since many phenotypes can only be measured in categories.

AVAILABILITY AND IMPLEMENTATION

https://github.com/mrborges23/delta_statistic.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

How Natural Genetic Variation Shapes Behavior

Nervous systems allow animals to acutely respond and behaviorally adapt to changes and recurring patterns in their environment at multiple timescales-from milliseconds to years. Behavior is further shaped at intergenerational timescales by genetic variation, drift, and selection. This sophistication and flexibility of behavior makes it challenging to measure behavior consistently in individual subjects and to compare it across individuals. In spite of these challenges, careful behavioral observations in nature and controlled measurements in the laboratory, combined with modern technologies and powerful genetic approaches, have led to important discoveries about the way genetic variation shapes behavior. A critical mass of genes whose variation is known to modulate behavior in nature is finally accumulating, allowing us to recognize emerging patterns. In this review, we first discuss genetic mapping approaches useful for studying behavior. We then survey how variation acts at different levels-in environmental sensation, in internal neuronal circuits, and outside the nervous system altogether-and then discuss the sources and types of molecular variation linked to behavior and the mechanisms that shape such variation. We end by discussing remaining questions in the field.

Avian Binocularity and Adaptation to Nocturnal Environments: Genomic Insights from a Highly Derived Visual Phenotype

Typical avian eyes are phenotypically engineered for photopic vision (daylight). In contrast, the highly derived eyes of the barn owl (Tyto alba) are adapted for scotopic vision (dim light). The dramatic modifications distinguishing barn owl eyes from other birds include: 1) shifts in frontal orientation to improve binocularity, 2) rod-dominated retina, and 3) enlarged corneas and lenses. Some of these features parallel mammalian eye patterns, which are hypothesized to have initially evolved in nocturnal environments. Here, we used an integrative approach combining phylogenomics and functional phenotypes of 211 eye-development genes across 48 avian genomes representing most avian orders, including the stem lineage of the scotopic-adapted barn owl. Overall, we identified 25 eye-development genes that coevolved under intensified or relaxed selection in the retina, lens, cornea, and optic nerves of the barn owl. The agtpbp1 gene, which is associated with the survival of photoreceptor populations, was pseudogenized in the barn owl genome. Our results further revealed that barn owl retinal genes responsible for the maintenance, proliferation, and differentiation of photoreceptors experienced an evolutionary relaxation. Signatures of relaxed selection were also observed in the lens and cornea morphology-associated genes, suggesting that adaptive evolution in these structures was essentially structural. Four eye-development genes (ephb1, phactr4, prph2, and rs1) evolved in positive association with the orbit convergence in birds and under relaxed selection in the barn owl lineage, likely contributing to an increased reliance on binocular vision in the barn owl. Moreover, we found evidence of coevolutionary interactions among genes that are expressed in the retina, lens, and optic nerve, suggesting synergetic adaptive events. Our study disentangles the genomic changes governing the binocularity and low-light perception adaptations of barn owls to nocturnal environments while revealing the molecular mechanisms contributing to the shift from the typical avian photopic vision to the more-novel scotopic-adapted eye.

Signature of adaptive evolution in olfactory receptor genes in Cory's Shearwater supports molecular basis for smell in procellariiform seabirds

Olfactory receptors (ORs), encoded by the largest vertebrate multigene family, enable the detection of thousands of unique odorants in the environment and consequently play a critical role in species survival. Here, we advance our knowledge of OR gene evolution in procellariiform seabirds, an avian group which relies on the sense of olfaction for critical ecological functions. We built a cosmid library of Cory's Shearwater (Calonectris borealis) genomic DNA, a model species for the study of olfaction-based navigation, and sequence OR gene-positive cosmid clones with a combination of sequencing technologies. We identified 220 OR open reading frames, 20 of which are full length, intact OR genes, and found a large ratio of partial and pseudogenes to intact OR genes (2:1), suggestive of a dynamic mode of evolution. Phylogenetic analyses revealed that while a few genes cluster with those of other sauropsid species in a γ (gamma) clade that predates the divergence of different avian lineages, most genes belong to an avian-specific γ-c clade, within which sequences cluster by species, suggesting frequent duplication and/or gene conversion events. We identified evidence of positive selection on full length γ-c clade genes. These patterns are consistent with a key role of adaptation in the functional diversification of olfactory receptor genes in a bird lineage that relies extensively on olfaction.

LMAP_S: Lightweight Multigene Alignment and Phylogeny eStimation

Background

Recent advances in genome sequencing technologies and the cost drop in high-throughput sequencing continue to give rise to a deluge of data available for downstream analyses. Among others, evolutionary biologists often make use of genomic data to uncover phenotypic diversity and adaptive evolution in protein-coding genes. Therefore, multiple sequence alignments (MSA) and phylogenetic trees (PT) need to be estimated with optimal results. However, the preparation of an initial dataset of multiple sequence file(s) (MSF) and the steps involved can be challenging when considering extensive amount of data. Thus, it becomes necessary the development of a tool that removes the potential source of error and automates the time-consuming steps of a typical workflow with high-throughput and optimal MSA and PT estimations.

Results

We introduce LMAP_S (Lightweight Multigene Alignment and Phylogeny eStimation), a user-friendly command-line and interactive package, designed to handle an improved alignment and phylogeny estimation workflow: MSF preparation, MSA estimation, outlier detection, refinement, consensus, phylogeny estimation, comparison and editing, among which file and directory organization, execution, manipulation of information are automated, with minimal manual user intervention. LMAP_S was developed for the workstation multi-core environment and provides a unique advantage for processing multiple datasets. Our software, proved to be efficient throughout the workflow, including, the (unlimited) handling of more than 20 datasets.

Conclusions

We have developed a simple and versatile LMAP_S package enabling researchers to effectively estimate multiple datasets MSAs and PTs in a high-throughput fashion. LMAP_S integrates more than 25 software providing overall more than 65 algorithm choices distributed in five stages. At minimum, one FASTA file is required within a single input directory. To our knowledge, no other software combines MSA and phylogeny estimation with as many alternatives and provides means to find optimal MSAs and phylogenies. Moreover, we used a case study comparing methodologies that highlighted the usefulness of our software. LMAP_S has been developed as an open-source package, allowing its integration into more complex open-source bioinformatics pipelines. LMAP_S package is released under GPLv3 license and is freely available at https://lmap-s.sourceforge.io/.

Massive Loss of Olfactory Receptors But Not Trace Amine-Associated Receptors in the World's Deepest-Living Fish (Pseudoliparis swirei)

Olfactory receptor repertoires show highly dynamic evolution associated with ecological adaptations in different species. The Mariana snailfish (Pseudoliparis swirei) living below a depth of 6000 m in the Mariana Trench evolved degraded vision and occupies a specific feeding habitat in a dark, low-food environment. However, whether such adaptations involve adaptive changes in the chemosensory receptor repertoire is not known. Here, we conducted a comparative analysis of the olfactory receptor (OR) and trace amine-associated receptor (TAAR) gene repertoires in nine teleosts with a focus on the evolutionary divergence between the Mariana snailfish and its shallow-sea relative, Tanaka's snailfish (Liparis tanakae). We found many fewer functional OR genes and a significantly higher fraction of pseudogenes in the Mariana snailfish, but the numbers of functional TAAR genes in the two species were comparable. Phylogenetic analysis showed that the expansion patterns of the gene families were shared by the two species, but that Mariana snailfish underwent massive gene losses in its OR repertoire. Despite an overall decreased size in OR subfamilies and a reduced number of TAAR subfamilies in the Mariana snailfish, expansion of certain subfamilies was observed. Selective pressure analysis indicated greatly relaxed selective strength in ORs but a slightly enhanced selective strength in TAARs of Mariana snailfish. Overall, our study reveals simplified but specific OR and TAAR repertoires in the Mariana snailfish shaped by natural selection with respect to ecological adaptations in the hadal environment. This is the first study on the chemosensation evolution in vertebrates living in the hadal zone, which could provide new insights into evolutionary adaptation to the hadal environment.

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.

The Vertebrate TLR Supergene Family Evolved Dynamically by Gene Gain/Loss and Positive Selection Revealing a Host–Pathogen Arms Race in Birds

The vertebrate toll-like receptor (TLRs) supergene family is a first-line immune defense against viral and non-viral pathogens. Here, comparative evolutionary-genomics of 79 vertebrate species (8 mammals, 48 birds, 11 reptiles, 1 amphibian, and 11 fishes) revealed differential gain/loss of 26 TLRs, including 6 (TLR3, TLR7, TLR8, TLR14, TLR21, and TLR22) that originated early in vertebrate evolution before the diversification of Agnatha and Gnathostomata. Subsequent dynamic gene gain/loss led to lineage-specific diversification with TLR repertoires ranging from 8 subfamilies in birds to 20 in fishes. Lineage-specific loss of TLR8-9 and TLR13 in birds and gains of TLR6 and TLR10-12 in mammals and TLR19-20 and TLR23-27 in fishes. Among avian species, 5–10% of the sites were under positive selection (PS) (omega 1.5–2.5) with radical amino-acid changes likely affecting TLR structure/functionality. In non-viral TLR4 the 20 PS sites (posterior probability PP > 0.99) likely increased ability to cope with diversified ligands (e.g., lipopolysaccharide and lipoteichoic). For viral TLR7, 23 PS sites (PP > 0.99) possibly improved recognition of highly variable viral ssRNAs. Rapid evolution of the TLR supergene family reflects the host–pathogen arms race and the coevolution of ligands/receptors, which follows the premise that birds have been important vectors of zoonotic pathogens and reservoirs for viruses.

Exploring the molecular basis of neuronal excitability in a vocal learner

BACKGROUND

Vocal learning, the ability to learn to produce vocalizations through imitation, relies on specialized brain circuitry known in songbirds as the song system. While the connectivity and various physiological properties of this system have been characterized, the molecular genetic basis of neuronal excitability in song nuclei remains understudied. We have focused our efforts on examining voltage-gated ion channels to gain insight into electrophysiological and functional features of vocal nuclei. A previous investigation of potassium channel genes in zebra finches (Taeniopygia guttata) revealed evolutionary modifications unique to songbirds, as well as transcriptional specializations in the song system [Lovell PV, Carleton JB, Mello CV. BMC Genomics 14:470 2013]. Here, we expand this approach to sodium, calcium, and chloride channels along with their modulatory subunits using comparative genomics and gene expression analysis encompassing microarrays and in situ hybridization.

RESULTS

We found 23 sodium, 38 calcium, and 33 chloride channel genes (HGNC-based classification) in the zebra finch genome, several of which were previously unannotated. We determined 15 genes are missing relative to mammals, including several genes (CLCAs, BEST2) linked to olfactory transduction. The majority of sodium and calcium but few chloride channels showed differential expression in the song system, among them SCN8A and CACNA1E in the direct motor pathway, and CACNG4 and RYR2 in the anterior forebrain pathway. In several cases, we noted a seemingly coordinated pattern across multiple nuclei (SCN1B, SCN3B, SCN4B, CACNB4) or sparse expression (SCN1A, CACNG5, CACNA1B).

CONCLUSION

The gene families examined are highly conserved between avian and mammalian lineages. Several cases of differential expression likely support high-frequency and burst firing in specific song nuclei, whereas cases of sparse patterns of expression may contribute to the unique electrophysiological signatures of distinct cell populations. These observations lay the groundwork for manipulations to determine how ion channels contribute to the neuronal excitability properties of vocal learning systems.

Expansion of vomeronasal receptor genes (OlfC) in the evolution of fright reaction in Ostariophysan fishes

Ostariophysans are the most diverse group of freshwater fishes and feature a pheromone-elicited fright reaction. However, the genetic basis of fright reaction is unclear. Here, we compared vomeronasal type 2 receptor-like (OlfC) genes from fishes having and lacking fright reaction, to provide insight into evolution of pheromonal olfaction in fishes. We found OlfC genes expanded remarkably in ostariophysans having fright reaction compared with fishes lacking fright reaction. Phylogenetic analysis indicates OlfC subfamily 9 expanded specifically in ostariophysans having fright reaction. Principle component and phylogenetic logistic regression analysis partitioned fishes by ecotype (having or lacking fright reaction) and identified OlfC subfamily 9 as being an important factor for fright reaction. Expression levels of expanded OlfC subfamily genes after fright reaction in zebrafish changed more than did genes that had not expanded. Furthermore, evidence of positive selection was found in the expanded OlfC proteins in ostariophysan fishes having fright reaction. These results provide new insight into the genetic basis of fright reaction in ostariophysan fish and will enable future research into the mechanism of action of OlfC proteins.

Olfactory receptor repertoire size in dinosaurs

The olfactory bulb (OB) ratio is the size of the OB relative to the cerebral hemisphere, and is used to estimate the proportion of the forebrain devoted to smell. In birds, OB ratio correlates with the number of olfactory receptor (OR) genes and therefore has been used as a proxy for olfactory acuity. By coupling OB ratios with known OR gene repertoires in birds, we infer minimum repertoire sizes for extinct taxa, including non-avian dinosaurs, using phylogenetic modelling, ancestral state reconstruction and comparative genomics. We highlight a shift in the scaling of OB ratio to body size along the lineage leading to modern birds, demonstrating variable OR repertoires present in different dinosaur and crown-bird lineages, with varying factors potentially influencing sensory evolution in theropods. We investigate the ancestral sensory space available to extinct taxa, highlighting potential adaptations to ecological niches. Through combining morphological and genomic data, we show that, while genetic information for extinct taxa is forever lost, it is potentially feasible to investigate evolutionary trajectories in extinct genomes.

Comparative analyses identify genomic features potentially involved in the evolution of birds-of-paradise

The diverse array of phenotypes and courtship displays exhibited by birds-of-paradise have long fascinated scientists and nonscientists alike. Remarkably, almost nothing is known about the genomics of this iconic radiation. There are 41 species in 16 genera currently recognized within the birds-of-paradise family (Paradisaeidae), most of which are endemic to the island of New Guinea. In this study, we sequenced genomes of representatives from all five major clades within this family to characterize genomic changes that may have played a role in the evolution of the group's extensive phenotypic diversity. We found genes important for coloration, morphology, and feather and eye development to be under positive selection. In birds-of-paradise with complex lekking systems and strong sexual dimorphism, the core birds-of-paradise, we found Gene Ontology categories for "startle response" and "olfactory receptor activity" to be enriched among the gene families expanding significantly faster compared to the other birds in our study. Furthermore, we found novel families of retrovirus-like retrotransposons active in all three de novo genomes since the early diversification of the birds-of-paradise group, which might have played a role in the evolution of this fascinating group of birds.

OR14I1 is a receptor for the human cytomegalovirus pentameric complex and defines viral epithelial cell tropism

A human cytomegalovirus (HCMV) pentameric glycoprotein complex (PC), gH-gL-UL128-UL130-UL131A, is necessary for viral infection of clinically relevant cell types, including epithelial cells, which are important for interhost transmission and disease. We performed genome-wide CRISPR/Cas9 screens of different cell types in parallel to identify host genes specifically required for HCMV infection of epithelial cells. This effort identified a multipass membrane protein, OR14I1, as a receptor for HCMV infection. This olfactory receptor family member is required for HCMV attachment, entry, and infection of epithelial cells and is dependent on the presence of viral PC. OR14I1 is required for AKT activation and mediates endocytosis entry of HCMV. We further found that HCMV infection of epithelial cells is blocked by a synthetic OR14I1 peptide and inhibitors of adenylate cyclase and protein kinase A (PKA) signaling. Identification of OR14I1 as a PC-dependent HCMV host receptor associated with epithelial tropism and the role of the adenylate cyclase/PKA/AKT-mediated signaling pathway in HCMV infection reveal previously unappreciated targets for the development of vaccines and antiviral therapies.

Genome-Wide Identification and Characterization of Olfactory Receptor Genes in Chinese Perch, Siniperca chuatsi

Olfaction, which is mediated by olfactory receptor (OR) genes, is essential in the daily life of fish, especially in foraging. However, Chinese perch (Siniperca chuatsi) is believed to prey with reliance on vision and lateral sensation, but not on olfaction. Therefore, understanding the evolutionary dynamics of the Chinese perch OR repertoire could provide insights into genetic evidence for adapting to a decreasing reliance on olfaction. Here, we reported a whole-genome analysis of the Chinese perch OR repertoire. Our analysis identified a total of 152 OR genes, including 123 functional genes and 29 pseudogenes, and showed their genomic organization. A phylogenetic tree was constructed, and the phylogenetic relationships of teleosts ORs was illustrated. The dN/dS (global ratios of non-synonymous to synonymous) analysis demonstrated that OR groups all appeared to be under purifying selection. Among the five Percomorpha fishes, Chinese perch only had 22 subfamilies, suggesting a decrease in OR diversities. The species-specific loss of subfamily 56 and 66 in Chinese perch, of which the genes belonged to subfamily 66, were orthologs of OR51E2, which recognized the plant odorant β-ionone, indicating that extremely piscivorous fish which might lose those receptors responded to plant-related odors. Finally, the expression profiles of OR genes in the olfactory epithelium at different developmental stages were investigated using RNA-seq data. From the aforementioned results, the evolution of the OR repertoire may be shaped by the adaption of vision-dependent specializations for foraging in Chinese perch. The first systematic study of OR genes in Chinese perch could provide valuable genomic resources for the further investigation of olfactory function in teleosts.

Molecular evolution of umami/sweet taste receptor genes in reptiles

Sensory systems play an important role in animal survival. Changes to these systems may be critical in evolution of species in new environments. Previous studies exploring the correlation between feeding ecology and Tas1r evolution mainly focused on mammals and birds, and found that the relationship was complex. However, in reptiles, the correlation between Tas1r evolution and dietary preferences is still unclear. Here, we attempted to explore this relationship in representative species of the major groups of reptiles (turtles, snakes, lizards, crocodilians), for which the genome information is known. We first predicted the functionality (intact, partial, or defective) of Tas1r, and then related it to the feeding preferences. As a result, we identified 11 Tas1r1, 12 Tas1r2, and 12 Tas1r3 genes to be partial or intact and another 22 Tas1r genes to be absent or pseudogenized in the 19 reptiles. We found that, as it was revealed in some other vertebrate groups, no correlation existed between feeding ecology and Tas1r evolution in reptiles: genomic prediction indicated that the Tas1r genes possibly have been lost or pseudogenized in snakes, but in crocodylia and testudines Tas1r genes are either intact or partial, regardless of their feeding habits. Thus, we suggest that the driving force of Tas1r evolution in reptiles is complex, and the feeding habit of swallowing food whole without chewing or the absence of taste buds in certain species may account for the possible umami/sweet perception loss. In addition, we propose that caution should be taken when predicting gene functionality from the publicly available genome database.

The Birth and Death of Olfactory Receptor Gene Families in Mammalian Niche Adaptation

The olfactory receptor (OR) gene families, which govern mammalian olfaction, have undergone extensive expansion and contraction through duplication and pseudogenization. Previous studies have shown that broadly defined environmental adaptations (e.g., terrestrial vs. aquatic) are correlated with the number of functional and non-functional OR genes retained. However, to date, no study has examined species-specific gene duplications in multiple phylogenetically divergent mammals to elucidate OR evolution and adaptation. Here, we identify the OR gene families driving adaptation to different ecological niches by mapping the fate of species-specific gene duplications in the OR repertoire of 94 diverse mammalian taxa, using molecular phylogenomic methods. We analyze >70,000 OR gene sequences mined from whole genomes, generated from novel amplicon sequencing data, and collated with data from previous studies, comprising one of the largest OR studies to date. For the first time, we demonstrate statistically significant patterns of OR species-specific gene duplications associated with the presence of a functioning vomeronasal organ. With respect to dietary niche, we uncover a novel link between a large number of duplications in OR family 5/8/9 and herbivory. Our results also highlight differences between social and solitary niches, indicating that a greater OR repertoire expansion may be associated with a solitary lifestyle. This study demonstrates the utility of species-specific duplications in elucidating gene family evolution, revealing how the OR repertoire has undergone expansion and contraction with respect to a number of ecological adaptations in mammals.

Avian genomics lends insights into endocrine function in birds

The genomics era has brought along the completed sequencing of a large number of bird genomes that cover a broad range of the avian phylogenetic tree (>30 orders), leading to major novel insights into avian biology and evolution. Among recent findings, the discovery that birds lack a large number of protein coding genes that are organized in highly conserved syntenic clusters in other vertebrates is very intriguing, given the physiological importance of many of these genes. A considerable number of them play prominent endocrine roles, suggesting that birds evolved compensatory genetic or physiological mechanisms that allowed them to survive and thrive in spite of these losses. While further studies are needed to establish the exact extent of avian gene losses, these findings point to birds as potentially highly relevant model organisms for exploring the genetic basis and possible therapeutic approaches for a wide range of endocrine functions and disorders.

Identifying human diamine sensors for death related putrescine and cadaverine molecules

Pungent chemical compounds originating from decaying tissue are strong drivers of animal behavior. Two of the best-characterized death smell components are putrescine (PUT) and cadaverine (CAD), foul-smelling molecules produced by decarboxylation of amino acids during decomposition. These volatile polyamines act as ‘necromones’, triggering avoidance or attractive responses, which are fundamental for the survival of a wide range of species. The few studies that have attempted to identify the cognate receptors for these molecules have suggested the involvement of the seven-helix trace amine-associated receptors (TAARs), localized in the olfactory epithelium. However, very little is known about the precise chemosensory receptors that sense these compounds in the majority of organisms and the molecular basis of their interactions. In this work, we have used computational strategies to characterize the binding between PUT and CAD with the TAAR6 and TAAR8 human receptors. Sequence analysis, homology modeling, docking and molecular dynamics studies suggest a tandem of negatively charged aspartates in the binding pocket of these receptors which are likely to be involved in the recognition of these small biogenic diamines.

The distinctive dead smell comes largely from molecules like cadaverine and putrescine that are produced during decomposition of organic tissues. These volatile compounds act as powerful chemical signals important for the survival of a wide range of species. Previous studies have identified the trace amine-associated receptor 13c (or TAAR13c) in zebrafish as the cognate receptor of cadaverine in bony fishes. In this work, we employed computational strategies to disclose the human TAAR6 and TAAR8 receptors as sensors of the putrescine and cadaverine molecules. Our results indicate that several negatively charged residues in the ligand binding pocket of these receptors constitute the molecular basis for recognition of these necromones in humans.

Datamonkey 2.0: A Modern Web Application for Characterizing Selective and Other Evolutionary Processes

Inference of how evolutionary forces have shaped extant genetic diversity is a cornerstone of modern comparative sequence analysis. Advances in sequence generation and increased statistical sophistication of relevant methods now allow researchers to extract ever more evolutionary signal from the data, albeit at an increased computational cost. Here, we announce the release of Datamonkey 2.0, a completely re-engineered version of the Datamonkey web-server for analyzing evolutionary signatures in sequence data. For this endeavor, we leveraged recent developments in open-source libraries that facilitate interactive, robust, and scalable web application development. Datamonkey 2.0 provides a carefully curated collection of methods for interrogating coding-sequence alignments for imprints of natural selection, packaged as a responsive (i.e. can be viewed on tablet and mobile devices), fully interactive, and API-enabled web application. To complement Datamonkey 2.0, we additionally release HyPhy Vision, an accompanying JavaScript application for visualizing analysis results. HyPhy Vision can also be used separately from Datamonkey 2.0 to visualize locally executed HyPhy analyses. Together, Datamonkey 2.0 and HyPhy Vision showcase how scientific software development can benefit from general-purpose open-source frameworks. Datamonkey 2.0 is freely and publicly available at http://www.datamonkey.org, and the underlying codebase is available from https://github.com/veg/datamonkey-js.

The Evolving Neural and Genetic Architecture of Vertebrate Olfaction

Evolution sculpts the olfactory nervous system in response to the unique sensory challenges facing each species. In vertebrates, dramatic and diverse adaptations to the chemical environment are possible because of the hierarchical structure of the olfactory receptor (OR) gene superfamily: expansion or contraction of OR subfamilies accompanies major changes in habitat and lifestyle; independent selection on OR subfamilies can permit local adaptation or conserved chemical communication; and genetic variation in single OR genes can alter odor percepts and behaviors driven by precise chemical cues. However, this genetic flexibility contrasts with the relatively fixed neural architecture of the vertebrate olfactory system, which requires that new olfactory receptors integrate into segregated and functionally distinct neural pathways. This organization allows evolution to couple critical chemical signals with selectively advantageous responses, but also constrains relationships between olfactory receptors and behavior. The coevolution of the OR repertoire and the olfactory system therefore reveals general principles of how the brain solves specific sensory problems and how it adapts to new ones.

Visualizing the GPCR Network: Classification and Evolution

In this study, we delineate an unsupervised clustering algorithm, minimum span clustering (MSC), and apply it to detect G-protein coupled receptor (GPCR) sequences and to study the GPCR network using a base dataset of 2770 GPCR and 652 non-GPCR sequences. High detection accuracy can be achieved with a proper dataset. The clustering results of GPCRs derived from MSC show a strong correlation between their sequences and functions. By comparing our level 1 MSC results with the GPCRdb classification, the consistency is 87.9% for the fourth level of GPCRdb, 89.2% for the third level, 98.4% for the second level, and 100% for the top level (the lowest resolution level of GPCRdb). The MSC results of GPCRs can be well explained by estimating the selective pressure of GPCRs, as exemplified by investigating the largest two subfamilies, peptide receptors (PRs) and olfactory receptors (ORs), in class A GPCRs. PRs are decomposed into three groups due to a positive selective pressure, whilst ORs remain as a single group due to a negative selective pressure. Finally, we construct and compare phylogenetic trees using distance-based and character-based methods, a combination of which could convey more comprehensive information about the evolution of GPCRs.