Comparative genomics reveals contraction in olfactory receptor genes in bats

An updated synthesis of and outstanding questions in the olfactory and vomeronasal systems in bats: Genetics asks questions only anatomy can answer

The extensive diversity observed in bat nasal chemosensory systems has been well-documented at the histological level. Understanding how this diversity evolved and developing hypotheses as to why particular patterns exist require a phylogenetic perspective, which was first outlined in the work of anatomist Kunwar Bhatnagar. With the onset of genetics and genomics, it might be assumed that the puzzling patterns observed in the morphological data have been clarified. However, there is still a widespread mismatch of genetic and morphological correlations among bat chemosensory systems. Novel genomic evidence has set up new avenues to explore that demand more evidence from anatomical structures. Here, we outline the progress that has been made in both morphological and molecular studies on the olfactory and vomeronasal systems in bats since the work of Bhatnagar. Genomic data of olfactory and vomeronasal receptors demonstrate the strong need for further morphological sampling, with a particular focus on receiving brain regions, glands, and ducts.

Next Generation Sequencing Revolutionizes Organismal Biology Research in Bats

The advent of next generation sequencing technologies (NGS) has greatly accelerated our understanding of critical aspects of organismal biology from non-model organisms. Bats form a particularly interesting group in this regard, as genomic data have helped unearth a vast spectrum of idiosyncrasies in bat genomes associated with bat biology, physiology, and evolution. Bats are important bioindicators and are keystone species to many eco-systems. They often live in proximity to humans and are frequently associated with emerging infectious diseases, including the COVID-19 pandemic. Nearly four dozen bat genomes have been published to date, ranging from drafts to chromosomal level assemblies. Genomic investigations in bats have also become critical towards our understanding of disease biology and host-pathogen coevolution. In addition to whole genome sequencing, low coverage genomic data like reduced representation libraries, resequencing data, etc. have contributed significantly towards our understanding of the evolution of natural populations, and their responses to climatic and anthropogenic perturbations. In this review, we discuss how genomic data have enhanced our understanding of physiological adaptations in bats (particularly related to ageing, immunity, diet, etc.), pathogen discovery, and host pathogen co-evolution. In comparison, the application of NGS towards population genomics, conservation, biodiversity assessment, and functional genomics has been appreciably slower. We reviewed the current areas of focus, identifying emerging topical research directions and providing a roadmap for future genomic studies in bats.

Comparative analyses of bat genomes identify distinct evolution of immunity in Old World fruit bats

Bats have been identified as natural reservoir hosts of several zoonotic viruses, prompting suggestions that they have unique immunological adaptations. Among bats, Old World fruit bats (Pteropodidae) have been linked to multiple spillovers. To test for lineage-specific molecular adaptations in these bats, we developed a new assembly pipeline to generate a reference-quality genome of the fruit bat Cynopterus sphinx and used this in comparative analyses of 12 bat species, including six pteropodids. Our results reveal that immunity-related genes have higher evolutionary rates in pteropodids than in other bats. Several lineage-specific genetic changes were shared across pteropodids, including the loss of NLRP1, duplications of PGLYRP1 and C5AR2, and amino acid replacements in MyD88. We introduced MyD88 transgenes containing Pteropodidae-specific residues into bat and human cell lines and found evidence of dampened inflammatory responses. By uncovering distinct immune adaptations, our results could help explain why pteropodids are frequently identified as viral hosts.

Molecular evolution and signatures of selective pressures on Bos, focusing on the Nelore breed (Bos indicus)

Evolutionary history leads to genome changes over time, especially for species that have experienced intense selective pressures over a short period. Here, we investigated the genomic evolution of Bos species by searching for potential selection signatures, focusing on Nelore, an economically relevant cattle breed in Brazil. We assessed the genomic processes determining the molecular evolution across Nelore and thirteen other related taxa by evaluating (i) amino acid sequence conservation, (ii) the dN/dS ratio, and (iii) gene families' turnover rate (λ). Low conserved regions potentially associated with fatty acid metabolism seem to reflect differences in meat fat content in taxa with different evolutionary histories. All Bos species presented genes under positive selection, especially B. indicus and Nelore, which include transport protein cobalamin, glycolipid metabolism, and hormone signaling. These findings could be explained by constant selective pressures to obtain higher immune resistance and efficient metabolism. The gene contraction rate across the Nelore + B. indicus branch was almost nine times higher than that in other lineages (λ = 0.01043 vs. 0.00121), indicating gene losses during the domestication process. Amino acid biosynthesis, reproductive and innate immune system-related pathways were associated with genes recognized within the most frequent rapidly evolving gene families and in genes under positive selection, supporting the substantial relevance of such traits from a domestication perspective. Our data provide new insights into how the genome may respond to intense artificial selection in distinct taxa, and reinforces the presence of selective pressures on traits potentially relevant for future animal breeding investments.

Sensory and Cognitive Ecology of Bats

We see stunning morphological diversity across the animal world. Less conspicuous but equally fascinating are the sensory and cognitive adaptations that determine animals’ interactions with their environments and each other. We discuss the development of the fields of sensory and cognitive ecology and the importance of integrating these fields to understand the evolution of adaptive behaviors. Bats, with their extraordinarily high ecological diversity, are ideal animals for this purpose. An explosion in recent research allows for better understanding of the molecular, genetic, neural, and behavioral bases for sensory ecology and cognition in bats. We give examples of studies that illuminate connections between sensory and cognitive features of information filtering, evolutionary trade-offs in sensory and cognitive processing, and multimodal sensing and integration. By investigating the selective pressures underlying information acquisition, processing, and use in bats, we aim to illuminate patterns and processes driving sensory and cognitive evolution.

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.

Novel de Novo Genome of Cynopterus brachyotis Reveals Evolutionarily Abrupt Shifts in Gene Family Composition across Fruit Bats

Major novel physiological or phenotypic adaptations often require accompanying modifications at the genic level. Conversely, the detection of considerable contractions and/or expansions of gene families can be an indicator of fundamental but unrecognized physiological change. We sequenced a novel fruit bat genome (Cynopterus brachyotis) and adopted a comparative approach to reconstruct the evolution of fruit bats, mapping contractions and expansions of gene families along their evolutionary history. Despite a radical change in life history as compared with other bats (e.g., loss of echolocation, large size, and frugivory), fruit bats have undergone surprisingly limited change in their genic composition, perhaps apart from a potentially novel gene family expansion relating to telomere protection and longevity. In sharp contrast, within fruit bats, the new Cynopterus genome bears the signal of unusual gene loss and gene family contraction, despite its similar morphology and lifestyle to two other major fruit bat lineages. Most missing genes are regulatory, immune-related, and olfactory in nature, illustrating the diversity of genomic strategies employed by bats to contend with responses to viral infection and olfactory requirements. Our results underscore that significant fluctuations in gene family composition are not always associated with obvious examples of novel physiological and phenotypic adaptations but may often relate to less-obvious shifts in immune strategies.

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.

Phylotastic: Improving Access to Tree-of-Life Knowledge With Flexible, on-the-Fly Delivery of Trees

A comprehensive phylogeny of species, i.e., a tree of life, has potential uses in a variety of contexts, including research, education, and public policy. Yet, accessing the tree of life typically requires special knowledge, complex software, or long periods of training. The Phylotastic project aims make it as easy to get a phylogeny of species as it is to get driving directions from mapping software. In prior work, we presented a design for an open system to validate and manage taxon names, find phylogeny resources, extract subtrees matching a user's taxon list, scale trees to time, and integrate related resources such as species images. Here, we report the implementation of a set of tools that together represent a robust, accessible system for on-the-fly delivery of phylogenetic knowledge. This set of tools includes a web portal to execute several customizable workflows to obtain species phylogenies (scaled by geologic time and decorated with thumbnail images); more than 30 underlying web services (accessible via a common registry); and code toolkits in R and Python (allowing others to develop custom applications using Phylotastic services). The Phylotastic system, accessible via http://www.phylotastic.org, provides a unique resource to access the current state of phylogenetic knowledge, useful for a variety of cases in which a tree extracted quickly from online resources (as distinct from a tree custom-made from character data) is sufficient, as it is for many casual uses of trees identified here.

Identifying orthologs with OMA: A primer

The Orthologous Matrix (OMA) is a method and database that allows users to identify orthologs among many genomes. OMA provides three different types of orthologs: pairwise orthologs, OMA Groups and Hierarchical Orthologous Groups (HOGs). This Primer is organized in two parts. In the first part, we provide all the necessary background information to understand the concepts of orthology, how we infer them and the different subtypes of orthology in OMA, as well as what types of analyses they should be used for. In the second part, we describe protocols for using the OMA browser to find a specific gene and its various types of orthologs. By the end of the Primer, readers should be able to (i) understand homology and the different types of orthologs reported in OMA, (ii) understand the best type of orthologs to use for a particular analysis; (iii) find particular genes of interest in the OMA browser; and (iv) identify orthologs for a given gene.  The data can be freely accessed from the OMA browser at https://omabrowser.org.

OMA standalone: orthology inference among public and custom genomes and transcriptomes

Genomes and transcriptomes are now typically sequenced by individual laboratories but analyzing them often remains challenging. One essential step in many analyses lies in identifying orthologs—corresponding genes across multiple species—but this is far from trivial. The Orthologous MAtrix (OMA) database is a leading resource for identifying orthologs among publicly available, complete genomes. Here, we describe the OMA pipeline available as a standalone program for Linux and Mac. When run on a cluster, it has native support for the LSF, SGE, PBS Pro, and Slurm job schedulers and can scale up to thousands of parallel processes. Another key feature of OMA standalone is that users can combine their own data with existing public data by exporting genomes and precomputed alignments from the OMA database, which currently contains over 2100 complete genomes. We compare OMA standalone to other methods in the context of phylogenetic tree inference, by inferring a phylogeny of Lophotrochozoa, a challenging clade within the protostomes. We also discuss other potential applications of OMA standalone, including identifying gene families having undergone duplications/losses in specific clades, and identifying potential drug targets in nonmodel organisms. OMA standalone is available under the permissive open source Mozilla Public License Version 2.0.

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.

The Egyptian Rousette Genome Reveals Unexpected Features of Bat Antiviral Immunity

Bats harbor many viruses asymptomatically, including several notorious for causing extreme virulence in humans. To identify differences between antiviral mechanisms in humans and bats, we sequenced, assembled, and analyzed the genome of Rousettus aegyptiacus, a natural reservoir of Marburg virus and the only known reservoir for any filovirus. We found an expanded and diversified KLRC/KLRD family of natural killer cell receptors, MHC class I genes, and type I interferons, which dramatically differ from their functional counterparts in other mammals. Such concerted evolution of key components of bat immunity is strongly suggestive of novel modes of antiviral defense. An evaluation of the theoretical function of these genes suggests that an inhibitory immune state may exist in bats. Based on our findings, we hypothesize that tolerance of viral infection, rather than enhanced potency of antiviral defenses, may be a key mechanism by which bats asymptomatically host viruses that are pathogenic in humans.