Intertumoral heterogeneity impacts oncolytic vesicular stomatitis virus efficacy in mouse pancreatic cancer cells.

Oncolytic virus (OV) therapy is a promising virus-based approach against various malignancies, including pancreatic ductal adenocarcinoma (PDAC). Our previous studies demonstrated that human PDAC cell lines are highly variable in their permissiveness to OVs. Mouse PDAC cell lines, which are widely used for in vivo examination of the adaptive immune responses during OV and other cancer therapies, have never been examined systematically for the impact of intertumoral heterogeneity (the differences observed between tumors in different patients) on OV virus efficacy. Here, we examined phenotypically and genotypically three commonly used allograftable mouse PDAC cell lines (C57BL6 genetic background): Panc02 (derived from chemically induced PDAC; also known as Pan02), and two cell lines originated from PDACs developed in two different KPC (KrasG12D, Trp53R172H, and PDX-1-Cre) mouse models. Our study (i) characterized the ability of a widely used attenuated oncolytic vesicular stomatitis virus VSV-ΔM51-GFP to infect, replicate in, and kill mouse PDAC cells; (ii) examined their innate antiviral responses; (iii) compared their permissiveness to a non-attenuated VSV-Mwt-GFP and chemotherapeutic drugs; and (iv) analyzed their karyotype and exome. Mouse PDAC cell lines showed high divergence in their permissiveness to VSV-ΔM51-GFP, which negatively correlated with their abilities to mount innate antiviral responses, while all three cell lines were highly permissive to VSV-Mwt-GFP. No correlation was found between resistance to VSV-ΔM51-GFP and chemotherapy. Also, mouse PDAC cell lines showed high divergence in their karyotype and exome. The exome analysis demonstrated that more VSV-ΔM51-GFP-permissive mouse PDAC cell lines harbor mutations in multiple important antiviral genes, such as TYK2, JAK2, and JAK3. IMPORTANCE Oncolytic virus (OV) therapy is a promising virus-based approach against various malignancies, including pancreatic ductal adenocarcinoma (PDAC). Our previous studies using various human PDAC cell lines demonstrated that they are highly variable in their permissiveness to OVs. In this study, we examined phenotypically and genotypically three commonly used allograftable mouse PDAC cell lines, which are widely used for in vivo examination of the adaptive immune responses during cancer therapies. Mouse PDAC cell lines showed high divergence in their permissiveness to oncolytic vesicular stomatitis virus (VSV), which negatively correlated with their abilities to mount innate antiviral responses. Also, we discovered that more VSV-permissive mouse PDAC cell lines harbor mutations in multiple important antiviral genes, such as TYK2, JAK2, and JAK3. Our study provides essential information about three model mouse PDAC cell lines and proposes a novel platform to study OV-based therapies against different PDACs in immunocompetent mice.

Disentangling Mechanical and Sensory Modules in the Radiation of Noctilionoid Bats.

AbstractWith diverse mechanical and sensory functions, the vertebrate cranium is a complex anatomical structure whose shifts between modularity and integration, especially in mechanical function, have been implicated in adaptive diversification. Yet how mechanical and sensory systems and their functions coevolve, as well as how their interrelationship contributes to phenotypic disparity, remain largely unexplored. To examine the modularity, integration, and evolutionary rates of sensory and mechanical structures within the head, we analyzed hard and soft tissue scans from ecologically diverse bats in the superfamily Noctilionoidea, a clade that ranges from insectivores and carnivores to frugivores and nectarivores. We identified eight regions that evolved in a coordinated fashion, thus recognizable as evolutionary modules: five associated with bite force and three linked to olfactory, visual, and auditory systems. Interrelationships among these modules differ between Neotropical leaf-nosed bats (family Phyllostomidae) and other noctilionoids. Consistent with the hypothesis that dietary transitions begin with changes in the capacity to detect novel food items followed by adaptations to process them, peak rates of sensory module evolution predate those of some mechanical modules. We propose that the coevolution of structures influencing bite force, olfaction, vision, and hearing constituted a structural opportunity that allowed the phyllostomid ancestor to take advantage of existing ecological opportunities and contributed to the clade's remarkable radiation.

Placing human gene families into their evolutionary context.

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Dietary Diversification and Specialization in Neotropical Bats Facilitated by Early Molecular Evolution.

Dietary adaptation is a major feature of phenotypic and ecological diversification, yet the genetic basis of dietary shifts is poorly understood. Among mammals, Neotropical leaf-nosed bats (family Phyllostomidae) show unmatched diversity in diet; from a putative insectivorous ancestor, phyllostomids have radiated to specialize on diverse food sources including blood, nectar, and fruit. To assess whether dietary diversification in this group was accompanied by molecular adaptations for changing metabolic demands, we sequenced 89 transcriptomes across 58 species and combined these with published data to compare ∼13,000 protein coding genes across 66 species. We tested for positive selection on focal lineages, including those inferred to have undergone dietary shifts. Unexpectedly, we found a broad signature of positive selection in the ancestral phyllostomid branch, spanning genes implicated in the metabolism of all major macronutrients, yet few positively selected genes at the inferred switch to plantivory. Branches corresponding to blood- and nectar-based diets showed selection in loci underpinning nitrogenous waste excretion and glycolysis, respectively. Intriguingly, patterns of selection in metabolism genes were mirrored by those in loci implicated in craniofacial remodeling, a trait previously linked to phyllostomid dietary specialization. Finally, we show that the null model of the widely-used branch-site test is likely to be misspecified, with the implication that the test is too conservative and probably under-reports true cases of positive selection. Our findings point to a complex picture of adaptive radiation, in which the evolution of new dietary specializations has been facilitated by early adaptations combined with the generation of new genetic variation.

Large-scale genome sampling reveals unique immunity and metabolic adaptations in bats.

Comprising more than 1,400 species, bats possess adaptations unique among mammals including powered flight, unexpected longevity, and extraordinary immunity. Some of the molecular mechanisms underlying these unique adaptations includes DNA repair, metabolism and immunity. However, analyses have been limited to a few divergent lineages, reducing the scope of inferences on gene family evolution across the Order Chiroptera. We conducted an exhaustive comparative genomic study of 37 bat species, one generated in this study, encompassing a large number of lineages, with a particular emphasis on multi-gene family evolution across immune and metabolic genes. In agreement with previous analyses, we found lineage-specific expansions of the APOBEC3 and MHC-I gene families, and loss of the proinflammatory PYHIN gene family. We inferred more than 1,000 gene losses unique to bats, including genes involved in the regulation of inflammasome pathways such as epithelial defence receptors, the natural killer gene complex and the interferon-gamma induced pathway. Gene set enrichment analyses revealed genes lost in bats are involved in defence response against pathogen-associated molecular patterns and damage-associated molecular patterns. Gene family evolution and selection analyses indicate bats have evolved fundamental functional differences compared to other mammals in both innate and adaptive immune system, with the potential to enhance antiviral immune response while dampening inflammatory signalling. In addition, metabolic genes have experienced repeated expansions related to convergent shifts to plant-based diets. Our analyses support the hypothesis that, in tandem with flight, ancestral bats had evolved a unique set of immune adaptations whose functional implications remain to be explored.

Foraging shifts and visual preadaptation in ecologically diverse bats.

Changes in behaviour may initiate shifts to new adaptive zones, with physical adaptations for novel environments evolving later. While new mutations are commonly considered engines of adaptive change, sensory evolution enabling access to new resources might also arise from standing genetic diversity, and even gene loss. We examine the relative contribution of molecular adaptations, measured by positive and relaxed selection, acting on eye-expressed genes associated with shifts to new adaptive zones in ecologically diverse bats from the superfamily Noctilionoidea. Collectively, noctilionoids display remarkable ecological breadth, from highly divergent echolocation to flight strategies linked to specialized insectivory, the parallel evolution of diverse plant-based diets (e.g., nectar, pollen and fruit) from ancestral insectivory, and-unusually for echolocating bats-often have large, well-developed eyes. We report contrasting levels of positive selection in genes associated with the development, maintenance and scope of visual function, tracing back to the origins of noctilionoids and Phyllostomidae (the bat family with most dietary diversity), instead of during shifts to novel diets. Generalized plant visiting was not associated with exceptional molecular adaptation, and exploration of these novel niches took place in an ancestral phyllostomid genetic background. In contrast, evidence for positive selection in vision genes was found at subsequent shifts to either nectarivory or frugivory. Thus, neotropical noctilionoids that use visual cues for identifying food and roosts, as well as for orientation, were effectively preadapted, with subsequent molecular adaptations in nectar-feeding lineages and the subfamily Stenodermatinae of fig-eating bats fine-tuning pre-existing visual adaptations for specialized purposes.

Untangling taxonomic confusion and diversification patterns of the Streak-breasted Scimitar Babblers (Timaliidae: Pomatorhinus ruficollis complex) in southern Asia.

Biodiversity patterns in Asia are poorly understood due to inferences drawn from incorrect taxonomy and limited survey effort. The Streak-breasted Scimitar Babblers (Pomatorhinus ruficollis complex) have a wide distribution across southern Asia and exhibit a high degree of plumage variability within and among populations. Continued use of traditional subspecies designations over revised species-limits based on plumage and DNA sequence data furthers confusion in this group and obscures complex biogeographic patterns. In this study, we combined previously published DNA sequences with newly sampled populations to produce a comprehensive dataset for the ruficollis complex. Phylogenetic analysis of these data confirms that traditional subspecies based on plumage alone are paraphyletic and therefore not good descriptors of evolutionary history. With increased sampling, our study supported previous delimitations of phylogenetic species as distinct units, refined the range limits of two taxa - P. reconditus (throughout central China) and P. nigrostellatus (Hainan, Guangxi, N Vietnam), showed two additional clades that may be distinct species, and uncovered a 'suture' zone where populations of multiple species occur in the same localities. Diversification within the ruficollis complex indicates a clade of Sino-Himalayan and SE Asian species sister to a clade distributed in central and southern China species. The 'suture' zone where different ruficollis species are in contact coincides with the meeting of these four major geographic areas in a highly geomorphologically complex region.

Bayesian hierarchical models suggest oldest known plant-visiting bat was omnivorous.

The earliest record of plant visiting in bats dates to the Middle Miocene of La Venta, the world's most diverse tropical palaeocommunity. Palynephyllum antimaster is known from molars that indicate nectarivory. Skull length, an important indicator of key traits such as body size, bite force and trophic specialization, remains unknown. We developed Bayesian models to infer skull length based on dental measurements. These models account for variation within and between species, variation between clades, and phylogenetic error structure. Models relating skull length to trophic level for nectarivorous bats were then used to infer the diet of the fossil. The skull length estimate for Palynephyllum places it among the larger lonchophylline bats. The inferred diet suggests Palynephyllum fed on nectar and insects, similar to its living relatives. Omnivory has persisted since the mid-Miocene. This is the first study to corroborate with fossil data that highly specialized nectarivory in bats requires an omnivorous transition.

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.

The bony cap and its distinction from the distal phalanx in humans, cats, and horses.

It has been recognized as early as the Victorian era that the apex of the distal phalanx has a distinct embryological development from the main shaft of the distal phalanx. Recent studies in regenerative medicine have placed an emphasis on the role of the apex of the distal phalanx in bone regrowth. Despite knowledge about the unique aspects of the distal phalanx, all phalanges are often treated as equivalent. Our morphological study reiterates and highlights the special anatomical and embryological properties of the apex of the distal phalanx, and names the apex "the bony cap" to distinguish it. We posit that the distal phalanx shaft is endochondral, while the bony cap is intramembranous and derived from the ectodermal wall. During development, the bony cap may be a separate structure that will fuse to the endochondral distal phalanx in the adult, as it ossifies well before the distal phalanges across taxa. Our study describes and revives the identity of the bony cap, and we identify it in three mammalian species: humans, cats, and horses (Homo sapiens, Felis catus domestica, and Equus caballus). During the embryonic period, we show the bony cap has a thimble-like shape that surrounds the proximal endochondral distal phalanx. The bony cap may thus play an inductive role in the differentiation of the corresponding nail, claw, or hoof (keratin structures) of the digit. When it is not present or develops erroneously, the corresponding keratin structures are affected, and regeneration is inhibited. By terming the bony cap, we hope to inspire more attention to its distinct identity and role in regeneration.

Cochlea development shapes bat sensory system evolution.

Sensory organs must develop alongside the skull within which they are largely encased, and this relationship can manifest as the skull constraining the organs, organs constraining the skull, or organs constraining one another in relative size. How this interplay between sensory organs and the developing skull plays out during the evolution of sensory diversity; however, remains unknown. Here, we examine the developmental sequence of the cochlea, the organ responsible for hearing and echolocation, in species with distinct diet and echolocation types within the ecologically diverse bat super-family Noctilionoidea. We found the size and shape of the cochlea largely correlates with skull size, with exceptions of Pteronotus parnellii, whose high duty cycle echolocation (nearly constant emission of sound pulses during their echolocation process allowing for detailed information gathering, also called constant frequency echolocation) corresponds to a larger cochlear and basal turn, and Monophyllus redmani, a small-bodied nectarivorous bat, for which interactions with other sensory organs restrict cochlea size. Our findings support the existence of developmental constraints, suggesting that both developmental and anatomical factors may act synergistically during the development of sensory systems in noctilionoid bats.

Aggression is suppressed by acute stress but induced by chronic stress: immobilization effects on aggression, hormones, and cortical 5-HT(1B)/ striatal dopamine D(2) receptor density.

Although it has been established by a number of investigators that a variety of stressors are associated with the induction of aggressive behavior, two specific issues remain unanswered. First, it is unclear whether the contexts surrounding stressors (e.g., stressor length and chance of winning over opponents) change outcomes regarding aggressive behavior. Second, if a relationship exists between stress and aggressive behavior, altered levels of stress-related hormone (e.g., corticosterone [CORT]), as well as aggression-related biomarkers (e.g., testosterone [T], density of prefronto-cortical 5-HT(1B) receptor and striatal dopamine D(2) receptor [D2r]) may contribute to changes in aggressive behavior. Thus, we examined how immobilization (with a 1-, 5-, or 10-day exposure) would impact (1) a longitudinal course of aggression toward different-sized opponents, (2) levels of CORT and T, and (3) densities of 5-HT(1B) receptor (5-HT1Br) in the prefrontal cortex (PFC) and D2r in the striatum. It was found that, regardless of small or large opponents, a single 2-h exposure to immobilization reduced aggressive behavior (stress-suppressed aggression) over time, whereas repeated (10-day) exposure to immobilization escalated aggressive behavior (stress-induced aggression). These stress effects persisted up to 1 week of recovery from immobilization stress. Moreover, immobilized rats demonstrated elevated levels of T, but not CORT, as compared with controls. Finally, acute immobilization altered D2r densities in the shell of the nucleus accumbens, and chronic immobilization changed 5-HT1Br in the PFC, including the downregulation of 5-HT1Br densities in the right prelimbic and orbitolateral cortices. The potential relationships among stress, aggression, and 5-HT1Br/D2r roles are discussed.

Ecological constraints on highly evolvable olfactory receptor genes and morphology in neotropical bats.

Although evolvability of genes and traits may promote specialization during species diversification, how ecology subsequently restricts such variation remains unclear. Chemosensation requires animals to decipher a complex chemical background to locate fitness-related resources, and thus the underlying genomic architecture and morphology must cope with constant exposure to a changing odorant landscape; detecting adaptation amidst extensive chemosensory diversity is an open challenge. In phyllostomid bats, an ecologically diverse clade that evolved plant visiting from a presumed insectivorous ancestor, the evolution of novel food detection mechanisms is suggested to be a key innovation, as plant-visiting species rely strongly on olfaction, supplementarily using echolocation. If this is true, exceptional variation in underlying olfactory genes and phenotypes may have preceded dietary diversification. We compared olfactory receptor (OR) genes sequenced from olfactory epithelium transcriptomes and olfactory epithelium surface area of bats with differing diets. Surprisingly, although OR evolution rates were quite variable and generally high, they are largely independent of diet. Olfactory epithelial surface area, however, is relatively larger in plant-visiting bats and there is an inverse relationship between OR evolution rates and surface area. Relatively larger surface areas suggest greater reliance on olfactory detection and stronger constraint on maintaining an already diverse OR repertoire. Instead of the typical case in which specialization and elaboration are coupled with rapid diversification of associated genes, here the relevant genes are already evolving so quickly that increased reliance on smell has led to stabilizing selection, presumably to maintain the ability to consistently discriminate among specific odorants-a potential ecological constraint on sensory evolution.

Diversity in olfactory receptor repertoires is associated with dietary specialization in a genus of frugivorous bat.

Mammalian olfactory receptor genes (ORs) are a diverse family of genes encoding proteins that directly interact with environmental chemical cues. ORs evolve via gene duplication in a birth-death fashion, neofunctionalizing and pseudogenizing over time. Olfaction is a primary sense used for food detection in plant-visiting bats, but the relationship between dietary specialization and OR repertoire diversity is unclear. Within neotropical Leaf-nosed bats (Phyllostomidae), many lineages are plant specialists, and some have a distinct OR repertoire compared to insectivorous species. Yet, whether specialization on particular plant genera is associated with the evolution of specialized, less diverse OR repertoires has never been tested. Using targeted sequence capture, we sequenced the OR repertoires of three sympatric species of short-tailed fruit bats (Carollia), which vary in their degree of specialization on the fruits of Piper plants. We characterized orthologous vs duplicated receptors among Carollia species, and explored the diversity and redundancy of the receptor gene repertoire. At the species level, the most dedicated Piper specialist, Carollia castanea, had lower OR diversity compared to the two generalists (C. sowelli and C. perspicillata), but we discovered a few unique sets of ORs within C. castanea with high redundancy of similar gene duplicates. These unique receptors potentially enable C. castanea to detect Piper fruit odorants better than its two congeners. Carollia perspicillata, the species with the most generalist diet, had a higher diversity of intact receptors, suggesting the ability to detect a wider range of odorant molecules. Variation among ORs may be a factor in the coexistence of these sympatric species, facilitating the exploitation of different plant resources. Our study sheds light on how gene duplication and changes in OR diversity may play a role in dietary adaptations and underlie ecological interactions between bats and plants.

Find the food first: An omnivorous sensory morphotype predates biomechanical specialization for plant based diets in phyllostomid bats.

The role of mechanical morphologies in the exploitation of novel niche space is well characterized; however, the role of sensory structures in unlocking new niches is less clear. Here, we investigate the relationship between the evolution of sensory structures and diet during the radiation of noctilionoid bats. With a broad range of foraging ecologies and a well-supported phylogeny, noctilionoids constitute an ideal group for studying this relationship. We used diffusible iodine-based contrast enhanced computed tomography scans of 44 noctilionoid species to analyze relationships between the relative volumes of three sensory structures (olfactory bulbs, orbits, and cochleae) and diet. We found a positive relationship between frugivory and both olfactory and orbit size. However, we also found a negative relationship between nectarivory and cochlea size. Ancestral state estimates suggest that larger orbits and olfactory bulbs were present in the common ancestor of family Phyllostomidae, but not in other noctilionoid. This constellation of traits indicates a shift toward omnivory at the base of Phyllostomidae, predating their radiation into an exceptionally broad range of dietary niches. This is consistent with a scenario in which changes in sensory systems associated with foraging and feeding set the stage for subsequent morphological modification and diversification.

Nectar-feeding bats and birds show parallel molecular adaptations in sugar metabolism enzymes.

In most vertebrates, the demand for glucose as the primary substrate for cellular respiration is met by the breakdown of complex carbohydrates, or energy is obtained by protein and lipid catabolism. In contrast, a few bat and bird species have convergently evolved to subsist on nectar, a sugar-rich mixture of glucose, fructose, and sucrose.1-4 How these nectar-feeders have adapted to cope with life-long high sugar intake while avoiding the onset of metabolic syndrome and diabetes5-7 is not understood. We analyzed gene sequences obtained from 127 taxa, including 22 nectar-feeding bat and bird genera that collectively encompass four independent origins of nectarivory. We show these divergent taxa have undergone pervasive molecular adaptation in sugar catabolism pathways, including parallel selection in key glycolytic and fructolytic enzymes. We also uncover convergent amino acid substitutions in the otherwise evolutionarily conserved aldolase B (ALDOB), which catalyzes rate-limiting steps in fructolysis and glycolysis, and the mitochondrial gatekeeper pyruvate dehydrogenase (PDH), which links glycolysis and the tricarboxylic acid cycle. Metabolomic profile and enzyme functional assays are consistent with increased respiratory flux in nectar-feeding bats and help explain how these taxa can both sustain hovering flight and efficiently clear simple sugars. Taken together, our results indicate that nectar-feeding bats and birds have undergone metabolic adaptations that have enabled them to exploit a unique energy-rich dietary niche among vertebrates.

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.

Evaluating the performance of targeted sequence capture, RNA-Seq, and degenerate-primer PCR cloning for sequencing the largest mammalian multigene family.

Multigene families evolve from single-copy ancestral genes via duplication, and typically encode proteins critical to key biological processes. Molecular analyses of these gene families require high-confidence sequences, but the high sequence similarity of the members can create challenges for sequencing and downstream analyses. Focusing on the common vampire bat, Desmodus rotundus, we evaluated how different sequencing approaches performed in recovering the largest mammalian protein-coding multigene family: olfactory receptors (OR). Using the genome as a reference, we determined the proportion of intact protein-coding receptors recovered by: (a) amplicons from degenerate primers sequenced via Sanger technology, (b) RNA-Seq of the main olfactory epithelium, and (c) those genes captured with probes designed from transcriptomes of closely-related species. Our initial re-annotation of the high-quality vampire bat genome resulted in >400 intact OR genes, more than doubling the original estimate. Sanger-sequenced amplicons performed the poorest among the three approaches, detecting <33% of receptors in the genome. In contrast, the transcriptome reliably recovered >50% of the annotated genomic ORs, and targeted sequence capture recovered nearly 75% of annotated genes. Each sequencing approach assembled high-quality sequences, even if it did not recover all receptors in the genome. While some variation may be due to limitations of the study design (e.g., different individuals), variation among approaches was mostly caused by low coverage of some receptors rather than high rates of assembly error. Given this variability, we caution against using the counts of intact receptors per species to model the birth-death process of multigene families. Instead, our results support the use of orthologous sequences to explore and model the evolutionary processes shaping these genes.

Protocols for the Molecular Evolutionary Analysis of Membrane Protein Gene Duplicates.

Gene duplication is an important process in the evolution of gene content in eukaryotic genomes. Understanding when gene duplicates contribute new molecular functions to genomes through molecular adaptation is one important goal in comparative genomics. In large gene families, however, characterizing adaptation and neofunctionalization across species is challenging, as models have traditionally quantified the timing of duplications without considering underlying gene trees. This protocol combines multiple approaches to detect adaptation in protein duplicates at a phylogenetic scale. We include a description of models for gene tree-species tree reconciliation that enable different types of inference, as well as a practical guide to their use. Although simulation-based approaches successfully detect shifts in the rate of duplication/retention, the conflation between the duplication and retention processes, the distinct trajectories of duplicates under non-, sub-, and neofunctionalization, as well as dosage effects offer hitherto unexplored analytical avenues. We introduce mathematical descriptions of these probabilities and offer a road map to computational implementation whose starting point is parsimony reconciliation. Sequence evolution information based on the ratio of nonsynonymous to synonymous nucleotide substitution rates (dN/dS) can be combined with duplicate survival probabilities to better predict the emergence of new molecular functions in retained duplicates. Together, these methods enable characterization of potentially adaptive candidate duplicates whose neofunctionalization may contribute to phenotypic divergence across species.