A comparison of the human and chimpanzee olfactory receptor gene repertoires

Fermentation technology as a driver of human brain expansion

Brain tissue is metabolically expensive. Consequently, the evolution of humans' large brains must have occurred via concomitant shifts in energy expenditure and intake. Proposed mechanisms include dietary shifts such as cooking. Importantly, though, any new food source must have been exploitable by hominids with brains a third the size of modern humans'. Here, we propose the initial metabolic trigger of hominid brain expansion was the consumption of externally fermented foods. We define "external fermentation" as occurring outside the body, as opposed to the internal fermentation in the gut. External fermentation could increase the bioavailability of macro- and micronutrients while reducing digestive energy expenditure and is supported by the relative reduction of the human colon. We discuss the explanatory power of our hypothesis and survey external fermentation practices across human cultures to demonstrate its viability across a range of environments and food sources. We close with suggestions for empirical tests.

Comparative physiological anthropogeny: exploring molecular underpinnings of distinctly human phenotypes

Anthropogeny is a classic term encompassing transdisciplinary investigations of the origins of the human species. Comparative anthropogeny is a systematic comparison of humans and other living nonhuman hominids (so-called "great apes"), aiming to identify distinctly human features in health and disease, with the overall goal of explaining human origins. We begin with a historical perspective, briefly describing how the field progressed from the earliest evolutionary insights to the current emphasis on in-depth molecular and genomic investigations of "human-specific" biology and an increased appreciation for cultural impacts on human biology. While many such genetic differences between humans and other hominids have been revealed over the last two decades, this information remains insufficient to explain the most distinctive phenotypic traits distinguishing humans from other living hominids. Here we undertake a complementary approach of "comparative physiological anthropogeny," along the lines of the preclinical medical curriculum, i.e., beginning with anatomy and considering each physiological system and in each case considering genetic and molecular components that are relevant. What is ultimately needed is a systematic comparative approach at all levels from molecular to physiological to sociocultural, building networks of related information, drawing inferences, and generating testable hypotheses. The concluding section will touch on distinctive considerations in the study of human evolution, including the importance of gene-culture interactions.

Genetic and functional odorant receptor variation in the Homo lineage

Humans, Neanderthals, and Denisovans independently adapted to a wide range of geographic environments and their associated food odors. Using ancient DNA sequences, we explored the in vitro function of thirty odorant receptor genes in the genus Homo. Our extinct relatives had highly conserved olfactory receptor sequence, but humans did not. Variations in odorant receptor protein sequence and structure may have produced variation in odor detection and perception. Variants led to minimal changes in specificity but had more influence on functional sensitivity. The few Neanderthal variants disturbed function, whereas Denisovan variants increased sensitivity to sweet and sulfur odors. Geographic adaptations may have produced greater functional variation in our lineage, increasing our olfactory repertoire and expanding our adaptive capacity. Our survey of olfactory genes and odorant receptors suggests that our genus has a shared repertoire with possible local ecological adaptations.

A draft genome of Drung cattle reveals clues to its chromosomal fusion and environmental adaptation

Drung cattle (Bos frontalis) have 58 chromosomes, differing from the Bos taurus 2n = 60 karyotype. To date, its origin and evolution history have not been proven conclusively, and the mechanisms of chromosome fusion and environmental adaptation have not been clearly elucidated. Here, we assembled a high integrity and good contiguity genome of Drung cattle with 13.7-fold contig N50 and 4.1-fold scaffold N50 improvements over the recently published Indian mithun assembly, respectively. Speciation time estimation and phylogenetic analysis showed that Drung cattle diverged from Bos taurus into an independent evolutionary clade. Sequence evidence of centromere regions provides clues to the breakpoints in BTA2 and BTA28 centromere satellites. We furthermore integrated a circulation and contraction-related biological process involving 43 evolutionary genes that participated in pathways associated with the evolution of the cardiovascular system. These findings may have important implications for understanding the molecular mechanisms of chromosome fusion, alpine valleys adaptability and cardiovascular function.

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.

Sensory disruption and sensory inequities in the Anthropocene

Anthropogenic disruptions to animal sensory ecology are as old as our species. But what about the effect on human sensory ecology? Human sensory dysfunction is increasing globally at great economic and health costs (mental, physical, and social). Contemporary sensory problems are directly tied to human behavioral changes and activity as well as anthropogenic pollution. The evolutionary sensory ecology and anthropogenic disruptions to three human senses (vision, audition, olfaction) are examined along with the economic and health costs of functionally reduced senses and demographic risk factors contributing to impairment. The primary goals of the paper are (a) to sew an evolutionary and ecological thread through clinical narratives on sensory dysfunction that highlights the impact of the built environment on the senses, and (b) to highlight structural, demographic, and environmental injustices that create sensory inequities in risk and that promote health disparities.

An integrative understanding of comparative cognition: lessons from human brain evolution

A comprehensive understanding of animal cognition requires the integration of studies on behavior, electrophysiology, neuroanatomy, development, and genomics. Although studies of comparative cognition are receiving increasing attention from organismal biologists, most current studies focus on the comparison of behaviors and anatomical structures to understand their adaptative values. However, to understand the most potentially complex cognitive program of the human brain a greater synthesis of a multitude of disciplines is needed. In this review, we start with extensive neuroanatomic comparisons between humans and other primates. One likely specialization of the human brain is the expansion of neocortex, especially in regions for high-order cognition (e.g., prefrontal cortex). We then discuss how such an expansion can be linked to heterochrony of the brain developmental program, resulting in a greater number of neurons and enhanced computational capacity. Furthermore, alteration of gene expression in the human brain has been associated with positive selection in DNA sequences of gene regulatory regions. These results not only imply that genes associated with brain development are a major factor in the evolution of cognition, but also that high-quality whole-genome sequencing and gene manipulation techniques are needed for an integrative and functional understanding of comparative cognition in non-model organisms.

Therapeutic potential of ectopic olfactory and taste receptors

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

dbCID: a manually curated resource for exploring the driver indels in human cancer

While recent advances in next-generation sequencing technologies have enabled the creation of a multitude of databases in cancer genomic research, there is no comprehensive database focusing on the annotation of driver indels (insertions and deletions) yet. Therefore, we have developed the database of Cancer driver InDels (dbCID), which is a collection of known coding indels that likely to be engaged in cancer development, progression or therapy. dbCID contains experimentally supported and putative driver indels derived from manual curation of literature and is freely available online at http://bioinfo.ahu.edu.cn:8080/dbCID. Using the data deposited in dbCID, we summarized features of driver indels in four levels (gene, DNA, transcript and protein) through comparing with putative neutral indels. We found that most of the genes containing driver indels in dbCID are known cancer genes playing a role in tumorigenesis. Contrary to the expectation, the sequences affected by driver frameshift indels are not larger than those by neutral ones. In addition, the frameshift and inframe driver indels prefer to disrupt high-conservative regions both in DNA sequences and protein domains. Finally, we developed a computational method for discriminating cancer driver from neutral frameshift indels based on the deposited data in dbCID. The proposed method outperformed other widely used non-cancer-specific predictors on an external test set, which demonstrated the usefulness of the data deposited in dbCID. We hope dbCID will be a benchmark for improving and evaluating prediction algorithms, and the characteristics summarized here may assist with investigating the mechanism of indel–cancer association.

Effects of short indels on protein structure and function in human genomes

Insertions and deletions (indels) represent the second most common type of genetic variations in human genomes. Indels can be deleterious and contribute to disease susceptibility as recent genome sequencing projects revealed a large number of indels in various cancer types. In this study, we investigated the possible effects of small coding indels on protein structure and function, and the baseline characteristics of indels in 2504 individuals of 26 populations from the 1000 Genomes Project. We found that each population has a distinct pattern in genes with small indels. Frameshift (FS) indels are enriched in olfactory receptor activity while non-frameshift (NFS) indels are enriched in transcription-related proteins. Structural analysis of NFS indels revealed that they predominantly adopt coil or disordered conformations, especially in proteins with transcription-related NFS indels. These results suggest that the annotated coding indels from the 1000 Genomes Project, while contributing to genetic variations and phenotypic diversity, generally do not affect the core protein structures and have no deleterious effect on essential biological processes. In addition, we found that a number of reference genome annotations might need to be updated due to the high prevalence of annotated homozygous indels in the general population.

ORDB, HORDE, ODORactor and other on-line knowledge resources of olfactory receptor-odorant interactions

We present here an exploration of the evolution of three well-established, web-based resources dedicated to the dissemination of information related to olfactory receptors (ORs) and their functional ligands, odorants. These resources are: the Olfactory Receptor Database (ORDB), the Human Olfactory Data Explorer (HORDE) and ODORactor. ORDB is a repository of genomic and proteomic information related to ORs and other chemosensory receptors, such as taste and pheromone receptors. Three companion databases closely integrated with ORDB are OdorDB, ORModelDB and OdorMapDB; these resources are part of the SenseLab suite of databases (http://senselab.med.yale.edu). HORDE (http://genome.weizmann.ac.il/horde/) is a semi-automatically populated database of the OR repertoires of human and several mammals. ODORactor (http://mdl.shsmu.edu.cn/ODORactor/) provides information related to OR-odorant interactions from the perspective of the odorant. All three resources are connected to each other via web-links.

Database URL: http://senselab.med.yale.edu; http://genome.weizmann.ac.il/horde/; http://mdl.shsmu.edu.cn/ODORactor/

Olfactory Receptors in Non-Chemosensory Organs: The Nervous System in Health and Disease

Olfactory receptors (ORs) and down-stream functional signaling molecules adenylyl cyclase 3 (AC3), olfactory G protein α subunit (Gαolf), OR transporters receptor transporter proteins 1 and 2 (RTP1 and RTP2), receptor expression enhancing protein 1 (REEP1), and UDP-glucuronosyltransferases (UGTs) are expressed in neurons of the human and murine central nervous system (CNS). In vitro studies have shown that these receptors react to external stimuli and therefore are equipped to be functional. However, ORs are not directly related to the detection of odors. Several molecules delivered from the blood, cerebrospinal fluid, neighboring local neurons and glial cells, distant cells through the extracellular space, and the cells’ own self-regulating internal homeostasis can be postulated as possible ligands. Moreover, a single neuron outside the olfactory epithelium expresses more than one receptor, and the mechanism of transcriptional regulation may be different in olfactory epithelia and brain neurons. OR gene expression is altered in several neurodegenerative diseases including Parkinson’s disease (PD), Alzheimer’s disease (AD), progressive supranuclear palsy (PSP) and sporadic Creutzfeldt-Jakob disease (sCJD) subtypes MM1 and VV2 with disease-, region- and subtype-specific patterns. Altered gene expression is also observed in the prefrontal cortex in schizophrenia with a major but not total influence of chlorpromazine treatment. Preliminary parallel observations have also shown the presence of taste receptors (TASRs), mainly of the bitter taste family, in the mammalian brain, whose function is not related to taste. TASRs in brain are also abnormally regulated in neurodegenerative diseases. These seminal observations point to the need for further studies on ORs and TASRs chemoreceptors in the mammalian brain.

Role of ortho-retronasal olfaction in mammalian cortical evolution

Fossils of mammals and their extinct relatives among cynodonts give evidence of correlated transformations affecting olfaction as well as mastication, head movement, and ventilation, and suggest evolutionary coupling of these seemingly separate anatomical regions into a larger integrated system of ortho-retronasal olfaction. Evidence from paleontology and physiology suggests that ortho-retronasal olfaction played a critical role at three stages of mammalian cortical evolution: early mammalian brain development was driven in part by ortho-retronasal olfaction; the bauplan for neocortex had higher-level association functions derived from olfactory cortex; and human cortical evolution was enhanced by ortho-retronasal smell.

The human gene damage index as a gene-level approach to prioritizing exome variants

The protein-coding exome of a patient with a monogenic disease contains about 20,000 variants, only one or two of which are disease causing. We found that 58% of rare variants in the protein-coding exome of the general population are located in only 2% of the genes. Prompted by this observation, we aimed to develop a gene-level approach for predicting whether a given human protein-coding gene is likely to harbor disease-causing mutations. To this end, we derived the gene damage index (GDI): a genome-wide, gene-level metric of the mutational damage that has accumulated in the general population. We found that the GDI was correlated with selective evolutionary pressure, protein complexity, coding sequence length, and the number of paralogs. We compared GDI with the leading gene-level approaches, genic intolerance, and de novo excess, and demonstrated that GDI performed best for the detection of false positives (i.e., removing exome variants in genes irrelevant to disease), whereas genic intolerance and de novo excess performed better for the detection of true positives (i.e., assessing de novo mutations in genes likely to be disease causing). The GDI server, data, and software are freely available to noncommercial users from lab.rockefeller.edu/casanova/GDI.

Genomic approaches to studying human-specific developmental traits

Changes in developmental regulatory programs drive both disease and phenotypic differences among species. Linking human-specific traits to alterations in development is challenging, because we have lacked the tools to assay and manipulate regulatory networks in human and primate embryonic cells. This field was transformed by the sequencing of hundreds of genomes – human and non-human – that can be compared to discover the regulatory machinery of genes involved in human development. This approach has identified thousands of human-specific genome alterations in developmental genes and their regulatory regions. With recent advances in stem cell techniques, genome engineering, and genomics, we can now test these sequences for effects on developmental gene regulation and downstream phenotypes in human cells and tissues.

Led by the nose: Olfaction in primate feeding ecology

Olfaction, the sense of smell, was a latecomer to the systematic investigation of primate sensory ecology after long years in which it was considered to be of minor importance. This view shifted with the growing understanding of its role in social behavior and the accumulation of physiological studies demonstrating that the olfactory abilities of some primates are on a par with those of olfactory-dependent mammals such as dogs and rodents. Recent years have seen a proliferation of physiological, behavioral, anatomical, and genetic investigations of primate olfaction. These investigations have begun to shed light on the importance of olfaction in the process of food acquisition. However, integration of these works has been limited. It is therefore still difficult to pinpoint large-scale evolutionary scenarios, namely the functions that the sense of smell fulfills in primates' feeding ecology and the ecological niches that favor heavier reliance on olfaction. Here, we review available behavioral and physiological studies of primates in the field or captivity and try to elucidate how and when the sense of smell can help them acquire food.

Global Survey of Variation in a Human Olfactory Receptor Gene Reveals Signatures of Non-Neutral Evolution

Allelic variation at 4 loci in the human olfactory receptor gene OR7D4 is associated with perceptual variation in the sex steroid-derived odorants, androstenone, and androstadienone. Androstadienone has been linked with chemosensory identification whereas androstenone makes pork from uncastrated pigs distasteful ("boar taint"). In a sample of 2224 individuals from 43 populations, we identified 45 OR7D4 single nucleotide polymorphisms. Coalescent modeling of frequency-site-spectrum-based statistics identified significant deviation from neutrality in human OR7D4; individual populations with statistically significant deviations from neutrality include Gujarati, Beijing Han, Great Britain, Iberia, and Puerto Rico. Analysis of molecular variation values indicated statistically significant population differentiation driven mainly by the 4 alleles associated with androstenone perception variation; however, fixation values were low suggesting that genetic structure may not have played a strong role in creating these group divisions. We also studied OR7D4 in the genomes of extinct members of the human lineage: Altai Neandertal and Denisovan. No variants were identified in Altai but 2 were in Denisova, one of which is shared by modern humans and one of which is novel. A functional test of modern human and a synthesized mutant Denisova OR7D4 indicated no statistically significant difference in responses to androstenone between the 2 species. Our results suggest non-neutral evolution for an olfactory receptor gene.

Inference of Evolutionary Forces Acting on Human Biological Pathways

Because natural selection is likely to act on multiple genes underlying a given phenotypic trait, we study here the potential effect of ongoing and past selection on the genetic diversity of human biological pathways. We first show that genes included in gene sets are generally under stronger selective constraints than other genes and that their evolutionary response is correlated. We then introduce a new procedure to detect selection at the pathway level based on a decomposition of the classical McDonald-Kreitman test extended to multiple genes. This new test, called 2DNS, detects outlier gene sets and takes into account past demographic effects and evolutionary constraints specific to gene sets. Selective forces acting on gene sets can be easily identified by a mere visual inspection of the position of the gene sets relative to their two-dimensional null distribution. We thus find several outlier gene sets that show signals of positive, balancing, or purifying selection but also others showing an ancient relaxation of selective constraints. The principle of the 2DNS test can also be applied to other genomic contrasts. For instance, the comparison of patterns of polymorphisms private to African and non-African populations reveals that most pathways show a higher proportion of nonsynonymous mutations in non-Africans than in Africans, potentially due to different demographic histories and selective pressures.

Population diversity and adaptive evolution in keratinization genes: impact of environment in shaping skin phenotypes

Several studies have demonstrated the role of climatic factors in shaping skin phenotypes, particularly pigmentation. Keratinization is another well-designed feature of human skin, which is involved in modulating transepidermal water loss (TEWL). Although this physiological process is closely linked to climate, presently it is not clear whether genetic diversity is observed in keratinization and whether this process also responds to the environmental pressure. To address this, we adopted a multipronged approach, which involved analysis of 1) copy number variations in diverse Indian and HapMap populations from varied geographical regions; 2) genetic association with geoclimatic parameters in 61 populations of dbCLINE database in a set of 549 genes from four processes namely keratinization, pigmentation, epidermal differentiation, and housekeeping functions; 3) sequence divergence in 4,316 orthologous promoters and corresponding exonic regions of human and chimpanzee with macaque as outgroup, and 4) protein sequence divergence (Ka/Ks) across nine vertebrate classes, which differ in their extent of TEWL. Our analyses demonstrate that keratinization and epidermal differentiation genes are under accelerated evolution in the human lineage, relative to pigmentation and housekeeping genes. We show that this entire pathway may have been driven by environmental selection pressure through concordant functional polymorphisms across several genes involved in skin keratinization. Remarkably, this underappreciated function of skin may be a crucial determinant of adaptation to diverse environmental pressures across world populations.

Mammalian olfactory receptors: molecular mechanisms of odorant detection, 3D-modeling, and structure-activity relationships

This chapter describes the main characteristics of olfactory receptor (OR) genes of vertebrates, including generation of this large multigenic family and pseudogenization. OR genes are compared in relation to evolution and among species. OR gene structure and selection of a given gene for expression in an olfactory sensory neuron (OSN) are tackled. The specificities of OR proteins, their expression, and their function are presented. The expression of OR proteins in locations other than the nasal cavity is regulated by different mechanisms, and ORs display various additional functions. A conventional olfactory signal transduction cascade is observed in OSNs, but individual ORs can also mediate different signaling pathways, through the involvement of other molecular partners and depending on the odorant ligand encountered. ORs are engaged in constitutive dimers. Ligand binding induces conformational changes in the ORs that regulate their level of activity depending on odorant dose. When present, odorant binding proteins induce an allosteric modulation of OR activity. Since no 3D structure of an OR has been yet resolved, modeling has to be performed using the closest G-protein-coupled receptor 3D structures available, to facilitate virtual ligand screening using the models. The study of odorant binding modes and affinities may infer best-bet OR ligands, to be subsequently checked experimentally. The relationship between spatial and steric features of odorants and their activity in terms of perceived odor quality are also fields of research that development of computing tools may enhance.

Positive Darwinian selection in the singularly large taste receptor gene family of an 'ancient' fish, Latimeria chalumnae

Background

Chemical senses are one of the foremost means by which organisms make sense of their environment, among them the olfactory and gustatory sense of vertebrates and arthropods. Both senses use large repertoires of receptors to achieve perception of complex chemosensory stimuli. High evolutionary dynamics of some olfactory and gustatory receptor gene families result in considerable variance of chemosensory perception between species. Interestingly, both ora/v1r genes and the closely related t2r genes constitute small and rather conserved families in teleost fish, but show rapid evolution and large species differences in tetrapods. To understand this transition, chemosensory gene repertoires of earlier diverging members of the tetrapod lineage, i.e. lobe-finned fish such as Latimeria would be of high interest.

Results

We report here the complete T2R repertoire of Latimeria chalumnae, using thorough data mining and extensive phylogenetic analysis. Eighty t2r genes were identified, by far the largest family reported for any species so far. The genomic neighborhood of t2r genes is enriched in repeat elements, which may have facilitated the extensive gene duplication events resulting in such a large family. Examination of non-synonymous vs. synonymous substitution rates (dN/dS) suggests pronounced positive Darwinian selection in Latimeria T2Rs, conceivably ensuring efficient neo-functionalization of newly born t2r genes. Notably, both traits, positive selection and enrichment of repeat elements in the genomic neighborhood, are absent in the twenty v1r genes of Latimeria. Sequence divergence in Latimeria T2Rs and V1Rs is high, reminescent of the corresponding teleost families. Some conserved sequence motifs of Latimeria T2Rs and V1Rs are shared with the respective teleost but not tetrapod genes, consistent with a potential role of such motifs in detection of aquatic chemosensory stimuli.

Conclusions

The singularly large T2R repertoire of Latimeria may have been generated by facilitating local gene duplication via increased density of repeat elements, and efficient neofunctionalization via positive Darwinian selection.

The high evolutionary dynamics of tetrapod t2r gene families precedes the emergence of tetrapods, i.e. the water-to-land transition, and thus constitutes a basal feature of the lobe-finned lineage of vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-650) contains supplementary material, which is available to authorized users.

Databases in SenseLab for the genomics, proteomics, and function of olfactory receptors

We present here, the salient aspects of three databases: Olfactory Receptor Database (ORDB) is a repository of genomics and proteomics information of ORs; OdorDB stores information related to odorous compounds, specifically identifying those that have been shown to interact with olfactory rectors; and OdorModelDB disseminates information related to computational models of olfactory receptors (ORs). The data stored among these databases is integrated. Presented in this chapter are descriptions of these resources, which are part of the SenseLab suite of databases, a discussion of the computational infrastructure that enhances the efficacy of information storage, retrieval, dissemination, and automated data population from external sources.

Olfactory receptors in non-chemosensory tissues

Olfactory receptors (ORs) detect volatile chemicals that lead to the initial perception of smell in the brain. The olfactory receptor (OR) is the first protein that recognizes odorants in the olfactory signal pathway and it is present in over 1,000 genes in mice. It is also the largest member of the G protein-coupled receptors (GPCRs). Most ORs are extensively expressed in the nasal olfactory epithelium where they perform the appropriate physiological functions that fit their location. However, recent whole-genome sequencing shows that ORs have been found outside of the olfactory system, suggesting that ORs may play an important role in the ectopic expression of non-chemosensory tissues. The ectopic expressions of ORs and their physiological functions have attracted more attention recently since MOR23 and testicular hOR17-4 have been found to be involved in skeletal muscle development, regeneration, and human sperm chemotaxis, respectively. When identifying additional expression profiles and functions of ORs in non-olfactory tissues, there are limitations posed by the small number of antibodies available for similar OR genes. This review presents the results of a research series that identifies ectopic expressions and functions of ORs in non-chemosensory tissues to provide insight into future research directions. [BMB Reports 2012; 45(11): 612-622]

Membrane environment imposes unique selection pressures on transmembrane domains of G protein-coupled receptors

We have investigated the influence of the plasma membrane environment on the molecular evolution of G protein-coupled receptors (GPCRs), the largest receptor family in Metazoa. In particular, we have analyzed the site-specific rate variation across the two primary structural partitions, transmembrane (TM) and extramembrane (EM), of these membrane proteins. We find that TM domains evolve more slowly than do EM domains, though TM domains display increased rate heterogeneity relative to their EM counterparts. Although the majority of residues across GPCRs experience strong to weak purifying selection, many GPCRs experience positive selection at both TM and EM residues, albeit with a slight bias towards the EM. Further, a subset of GPCRs, chemosensory receptors (including olfactory and taste receptors), exhibit increased rates of evolution relative to other GPCRs, an effect which is more pronounced in their TM spans. Although it has been previously suggested that the TM's low evolutionary rate is caused by their high percentage of buried residues, we show that their attenuated rate seems to stem from the strong biophysical constraints of the membrane itself, or by functional requirements. In spite of the strong evolutionary constraints acting on the TM spans of GPCRs, positive selection and high levels of evolutionary rate variability are common. Thus, biophysical constraints should not be presumed to preclude a protein's ability to evolve.

HORDE: comprehensive resource for olfactory receptor genomics

Olfactory receptors (ORs) constitute the largest gene family in the mammalian genome. The existence of these proteins underlies the nature of, and variability in, odorant perception. The Human Olfactory Receptor Data Explorer (HORDE, http://genome.weizmann.ac.il/horde/ ) is a free online resource, which presents a complete compendium of all OR genes and pseudogenes in the genome of human and four other vertebrates. HORDE includes three parts: (1) an automated pipeline, which mines OR gene and pseudogene sequences out of complete genomes, and generates gene symbols based on sequence similarity; (2) a card generator that obtains and displays annotative information on individual ORs retrieved from external databases and relevant studies; and (3) a search engine that allows user retrieval of OR information. For human ORs, HORDE specifically addresses the universe of interindividual variation, as obtained from several sources, including whole genome sequences made possible by next-generation sequencing. This encompasses single nucleotide polymorphisms (SNP) and copy number variation (CNV), including deleterious mutational events. HORDE also hosts a number of tools designed specifically to assist in the study of OR evolution and function. In this chapter, we describe the status of HORDE (build #43). We also discuss plans for future enhancements and a road map for HORDE to become a better community-based bioinformatics tool. We highlight HORDE's role as a major research tool in the study of an expanding cohort of OR repertoires.

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

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

Functional evolution of mammalian odorant receptors

The mammalian odorant receptor (OR) repertoire is an attractive model to study evolution, because ORs have been subjected to rapid evolution between species, presumably caused by changes of the olfactory system to adapt to the environment. However, functional assessment of ORs in related species remains largely untested. Here we investigated the functional properties of primate and rodent ORs to determine how well evolutionary distance predicts functional characteristics. Using human and mouse ORs with previously identified ligands, we cloned 18 OR orthologs from chimpanzee and rhesus macaque and 17 mouse-rat orthologous pairs that are broadly representative of the OR repertoire. We functionally characterized the in vitro responses of ORs to a wide panel of odors and found similar ligand selectivity but dramatic differences in response magnitude. 87% of human-primate orthologs and 94% of mouse-rat orthologs showed differences in receptor potency (EC50) and/or efficacy (dynamic range) to an individual ligand. Notably dN/dS ratio, an indication of selective pressure during evolution, does not predict functional similarities between orthologs. Additionally, we found that orthologs responded to a common ligand 82% of the time, while human OR paralogs of the same subfamily responded to the common ligand only 33% of the time. Our results suggest that, while OR orthologs tend to show conserved ligand selectivity, their potency and/or efficacy dynamically change during evolution, even in closely related species. These functional changes in orthologs provide a platform for examining how the evolution of ORs can meet species-specific demands.

Predicting the effects of frameshifting indels

Each human has approximately 50 to 280 frameshifting indels, yet their implications are unknown. We created SIFT Indel, a prediction method for frameshifting indels that has 84% accuracy. The percentage of human frameshifting indels predicted to be gene-damaging is negatively correlated with allele frequency. We also show that although the first frameshifting indel in a gene causes loss of function, there is a tendency for the second frameshifting indel to compensate and restore protein function. SIFT Indel is available at http://sift-dna.org/www/SIFT_indels2.html

The evolutionary history of the SAL1 gene family in eutherian mammals

Background

SAL1 (salivary lipocalin) is a member of the OBP (Odorant Binding Protein) family and is involved in chemical sexual communication in pig. SAL1 and its relatives may be involved in pheromone and olfactory receptor binding and in pre-mating behaviour. The evolutionary history and the selective pressures acting on SAL1 and its orthologous genes have not yet been exhaustively described. The aim of the present work was to study the evolution of these genes, to elucidate the role of selective pressures in their evolution and the consequences for their functions.

Results

Here, we present the evolutionary history of SAL1 gene and its orthologous genes in mammals. We found that (1) SAL1 and its related genes arose in eutherian mammals with lineage-specific duplications in rodents, horse and cow and are lost in human, mouse lemur, bushbaby and orangutan, (2) the evolution of duplicated genes of horse, rat, mouse and guinea pig is driven by concerted evolution with extensive gene conversion events in mouse and guinea pig and by positive selection mainly acting on paralogous genes in horse and guinea pig, (3) positive selection was detected for amino acids involved in pheromone binding and amino acids putatively involved in olfactory receptor binding, (4) positive selection was also found for lineage, indicating a species-specific strategy for amino acid selection.

Conclusions

This work provides new insights into the evolutionary history of SAL1 and its orthologs. On one hand, some genes are subject to concerted evolution and to an increase in dosage, suggesting the need for homogeneity of sequence and function in certain species. On the other hand, positive selection plays a role in the diversification of the functions of the family and in lineage, suggesting adaptive evolution, with possible consequences for speciation and for the reinforcement of prezygotic barriers.

Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA

BACKGROUND

Artificial selection has caused rapid evolution in domesticated species. The identification of selection footprints across domesticated genomes can contribute to uncover the genetic basis of phenotypic diversity.

METHODOLOGY/MAIN FINDINGS

Genome wide footprints of pig domestication and selection were identified using massive parallel sequencing of pooled reduced representation libraries (RRL) representing ∼2% of the genome from wild boar and four domestic pig breeds (Large White, Landrace, Duroc and Pietrain) which have been under strong selection for muscle development, growth, behavior and coat color. Using specifically developed statistical methods that account for DNA pooling, low mean sequencing depth, and sequencing errors, we provide genome-wide estimates of nucleotide diversity and genetic differentiation in pig. Widespread signals suggestive of positive and balancing selection were found and the strongest signals were observed in Pietrain, one of the breeds most intensively selected for muscle development. Most signals were population-specific but affected genomic regions which harbored genes for common biological categories including coat color, brain development, muscle development, growth, metabolism, olfaction and immunity. Genetic differentiation in regions harboring genes related to muscle development and growth was higher between breeds than between a given breed and the wild boar.

CONCLUSIONS/SIGNIFICANCE

These results, suggest that although domesticated breeds have experienced similar selective pressures, selection has acted upon different genes. This might reflect the multiple domestication events of European breeds or could be the result of subsequent introgression of Asian alleles. Overall, it was estimated that approximately 7% of the porcine genome has been affected by selection events. This study illustrates that the massive parallel sequencing of genomic pools is a cost-effective approach to identify footprints of selection.

Comparative genetic approaches to the evolution of human brain and behavior

With advances in genomic technologies, the amount of genetic data available to scientists today is vast. Genomes are now available or planned for 14 different primate species and complete resequencing of numerous human individuals from numerous populations is underway. Moreover, high-throughput deep sequencing is quickly making whole genome efforts within the reach of single laboratories allowing for unprecedented studies. Comparative genetic approaches to the identification of the underlying basis of human brain, behavior, and cognitive ability are moving to the forefront. Two approaches predominate: inter-species divergence comparisons and intra-species polymorphism studies. These methodological differences are useful for different time scales of evolution and necessarily focus on different evolutionary events in the history of primate and hominin evolution. Inter-species divergence is more useful in studying large scale primate, or hominoid, evolution whereas intra-species polymorphism can be more illuminating of recent hominin evolution. These differences in methodological utility also extend to studies of differing genetic substrates; current divergence studies focus primarily on protein evolution whereas polymorphism studies are substrate ambivalent. Some of the issues inherent in these studies can be ameliorated by current sequencing capabilities whereas others remain intractable. New avenues are also being opened that allow for the incorporation of novel substrates and approaches. In the post-genomic era, the study of human evolution, specifically as it relates to the brain, is becoming more complete focusing increasingly on the totality of the system and better conceptualizing the entirety of the genetic changes that have lead to the human phenotype today.

Study of bioengineered zebra fish olfactory receptor 131-2: receptor purification and secondary structure analysis

How fishes are able to detect trace molecules in large bodies of water is not understood. It is plausible that they use olfactory receptors to detect water-soluble compounds. How the zebra fish Danio Rerio, an organism with only 98 functional olfactory receptors, is able to selectively detect and recognize numerous compounds in water remains a puzzling phenomenon. We are interested in studying the biochemical and molecular mechanisms of olfaction in fish. Here, we report on the study of a bioengineered zebra fish olfactory receptor OR131-2, affinity-purified from a HEK293S tetracycline-inducible system. This receptor was expressed and translocated to the cell plasma membrane as revealed by confocal microscopy. Circular dichroism spectroscopy showed that the purified zebra fish receptor folded into an α-helical structure, as observed for other G-protein coupled receptors (GPCRs). Our study shows that it is possible to produce viable quantities of the zebra fish olfactory receptor. This will not only enable detailed structural and functional analyses, but also aid in the design of biosensor devices in order to detect water-soluble metabolites or its intermediates, which are associated with human health.

Global features of neural activity in the olfactory system form a parallel code that predicts olfactory behavior and perception

Odor identity is coded in spatiotemporal patterns of neural activity in the olfactory bulb. Here we asked whether meaningful olfactory information could also be read from the global olfactory neural population response. We applied standard statistical methods of dimensionality-reduction to neural activity from 12 previously published studies using seven different species. Four studies reported olfactory receptor activity, seven reported glomerulus activity, and one reported the activity of projection-neurons. We found two linear axes of neural population activity that accounted for more than half of the variance in neural response across species. The first axis was correlated with the total sum of odor-induced neural activity, and reflected the behavior of approach or withdrawal in animals, and odorant pleasantness in humans. The second and orthogonal axis reflected odorant toxicity across species. We conclude that in parallel with spatiotemporal pattern coding, the olfactory system can use simple global computations to read vital olfactory information from the neural population response.

Olfactory receptor gene family evolution in stickleback and medaka fishes

Interaction of olfactory receptor (OR) genes with environmental odors is regarded as the first step of olfaction. In this study, OR genes of two fish, medaka (Oryzias latipes) and stickleback (Gasterosteus aculeatus), were identified and an evolutional analysis was conducted. The selection pressure of different TM regions and complete coding region were compared. Three TM regions (TM4, TM5 and TM6) were found to have higher average Ka/Ks values, which might be partly caused by positive selection as suggested by subsequent positive selection analysis. Further analysis showed that many PTSs overlap, or are adjacent to previously deduced binding sites in mammals. These results support the hypothesis that binding sites of fish OR genes may evolved under positive selection.

Evolutionary dynamics of olfactory receptor genes in chordates: interaction between environments and genomic contents

Olfaction is essential for the survival of animals. Versatile odour molecules in the environment are received by olfactory receptors (ORs), which form the largest multigene family in vertebrates. Identification of the entire repertories of OR genes using bioinformatics methods from the whole-genome sequences of diverse organisms revealed that the numbers of OR genes vary enormously, ranging from ~1,200 in rats and ~400 in humans to ~150 in zebrafish and ~15 in pufferfish. Most species have a considerable fraction of pseudogenes. Extensive phylogenetic analyses have suggested that the numbers of gene gains and losses are extremely large in the OR gene family, which is a striking example of the birth-and-death evolution. It appears that OR gene repertoires change dynamically, depending on each organism's living environment. For example, higher primates equipped with a well-developed vision system have lost a large number of OR genes. Moreover, two groups of OR genes for detecting airborne odorants greatly expanded after the time of terrestrial adaption in the tetrapod lineage, whereas fishes retain diverse repertoires of genes that were present in aquatic ancestral species. The origin of vertebrate OR genes can be traced back to the common ancestor of all chordate species, but insects, nematodes and echinoderms utilise distinctive families of chemoreceptors, suggesting that chemoreceptor genes have evolved many times independently in animal evolution.

Degeneration of olfactory receptor gene repertories in primates: no direct link to full trichromatic vision

Odor molecules in the environment are detected by olfactory receptors (ORs), being encoded by a large multigene family in mammalian genomes. It is generally thought that primates are vision oriented and dependent weakly on olfaction. Previous studies suggested that Old World monkeys (OWMs) and hominoids lost many functional OR genes after the divergence from New World monkeys (NWMs) due to the acquisition of well-developed trichromatic vision. To examine this hypothesis, here we analyzed OR gene repertoires of five primate species including NWMs, OWMs, and hominoids for which high-coverage genome sequences are available, together with two prosimians and tree shrews with low-coverage genomes. The results showed no significant differences in the number of functional OR genes between NWMs (marmosets) and OWMs/hominoids. Two independent analyses, identification of orthologous genes among the five primates and estimation of the numbers of ancestral genes by the reconciled tree method, did not support a sudden loss of OR genes at the branch of the OWMs/hominoids ancestor but suggested a gradual loss in every lineage. Moreover, we found that humans retain larger numbers of ancestral OR genes that were in the common ancestor of NWMs/OWMs/hominoids than orangutans and macaques and that the OR gene repertoire in humans is more similar to that of marmosets than those of orangutans and macaques. These results suggest that the degeneration of OR genes in primates cannot simply be explained by the acquisition of trichromatic vision, and our sense of smell may not be inferior to other primate species.

Dynamic functional evolution of an odorant receptor for sex-steroid-derived odors in primates

Odorant receptors are among the fastest evolving genes in animals. However, little is known about the functional changes of individual odorant receptors during evolution. We have recently demonstrated a link between the in vitro function of a human odorant receptor, OR7D4, and in vivo olfactory perception of 2 steroidal ligands--androstenone and androstadienone--chemicals that are shown to affect physiological responses in humans. In this study, we analyzed the in vitro function of OR7D4 in primate evolution. Orthologs of OR7D4 were cloned from different primate species. Ancestral reconstruction allowed us to reconstitute additional putative OR7D4 orthologs in hypothetical ancestral species. Functional analysis of these orthologs showed an extremely diverse range of OR7D4 responses to the ligands in various primate species. Functional analysis of the nonsynonymous changes in the Old World Monkey and Great Ape lineages revealed a number of sites causing increases or decreases in sensitivity. We found that the majority of the functionally important residues in OR7D4 were not predicted by the maximum likelihood analysis detecting positive Darwinian selection.

Phylogenomic analyses reveal convergent patterns of adaptive evolution in elephant and human ancestries

Specific sets of brain-expressed genes, such as aerobic energy metabolism genes, evolved adaptively in the ancestry of humans and may have evolved adaptively in the ancestry of other large-brained mammals. The recent addition of genomes from two afrotherians (elephant and tenrec) to the expanding set of publically available sequenced mammalian genomes provided an opportunity to test this hypothesis. Elephants resemble humans by having large brains and long life spans; tenrecs, in contrast, have small brains and short life spans. Thus, we investigated whether the phylogenomic patterns of adaptive evolution are more similar between elephant and human than between either elephant and tenrec lineages or human and mouse lineages, and whether aerobic energy metabolism genes are especially well represented in the elephant and human patterns. Our analyses encompassed approximately 6,000 genes in each of these lineages with each gene yielding extensive coding sequence matches in interordinal comparisons. Each gene's nonsynonymous and synonymous nucleotide substitution rates and dN/dS ratios were determined. Then, from gene ontology information on genes with the higher dN/dS ratios, we identified the more prevalent sets of genes that belong to specific functional categories and that evolved adaptively. Elephant and human lineages showed much slower nucleotide substitution rates than tenrec and mouse lineages but more adaptively evolved genes. In correlation with absolute brain size and brain oxygen consumption being largest in elephants and next largest in humans, adaptively evolved aerobic energy metabolism genes were most evident in the elephant lineage and next most evident in the human lineage.

Inferring the evolutionary history of gene clusters from phylogenetic and gene order data

Gene duplication is frequent within gene clusters and plays a fundamental role in evolution by providing a source of new genetic material upon which natural selection can act. Although classical phylogenetic inference methods provide some insight into the evolutionary history of a gene cluster, they are not sufficient alone to differentiate single- from multiple gene duplication events and to answer other questions regarding the nature and size of evolutionary events. In this paper, we present an algorithm allowing to infer a set of optimal evolutionary histories for a gene cluster in a single species, according to a general cost model involving variable length duplications (in tandem or inverted), deletions, and inversions. We applied our algorithm to the human olfactory receptor and protocadherin gene clusters, showing that the duplication size distribution differs significantly between the two gene families. The algorithm is available through a web interface at http://www-lbit.iro.umontreal.ca/DILTAG/.

Heterogeneous patterns of gene-expression diversification in mammalian gene duplicates

Gene duplication is a major mechanism for molecular evolutionary innovation. Young gene duplicates typically exhibit elevated rates of protein evolution and, according to a number of recent studies, increased expression divergence. However, the nature of these changes is still poorly understood. To gain novel insights into the functional consequences of gene duplication, we have undertaken an in-depth analysis of a large data set of gene families containing primate- and/or rodent-specific gene duplicates. We have found a clear tendency toward an increase in protein, promoter, and expression divergence with increasing number of duplication events undergone by each gene since the human-mouse split. In addition, gene duplication is significantly associated with a reduction in expression breadth and intensity. Interestingly, it is possible to identify three main groups regarding the evolution of gene expression following gene duplication. The first group, which comprises around 25% of the families, shows patterns compatible with tissue-expression partitioning. The second and largest group, comprising 33-53% of the families, shows broad expression of one of the gene copies and reduced, overlapping, expression of the other copy or copies. This can be attributed, in most cases, to loss of expression in several tissues of one or more gene copies. Finally, a substantial number of families, 19-35%, maintain a very high level of tissue-expression overlap (>0.8) after tens of millions of years of evolution. These families may have been subject to selection for increased gene dosage.

Genomics of olfactory receptors

In many species, the sense of smell plays important roles in locating food, detecting predators, navigating, and communicating social information. The olfactory system has evolved complex repertoires of odor receptors (ORs) to fulfill these functions. Through computational data mining, OR repertoires of multiple species were identified, revealing a surprisingly large OR gene family in rodents and evolutionary fluctuation among different organisms. Characteristics of OR genes were explored through computational and experimental methods, showing a complicated gene structure and special genomic distribution. Utilizing high-throughput OR microarrays, expression profiles of the mouse and human OR repertoire were examined, their olfactory functions verified, and their zonal, ectopic and developmental expression determined. Variation in human smelling abilities results from different functional OR repertoires, variable expressional levels and polymorphisms in the copy number of the OR genes. These genomic approaches have both provided new data and generated new questions.

Evolution of olfactory receptor genes in primates dominated by birth-and-death process

Olfactory receptor (OR) is a large family of G protein–coupled receptors that can detect odorant in order to generate the sense of smell. They constitute one of the largest multiple gene families in animals including primates. To better understand the variation in odor perception and evolution of OR genes among primates, we computationally identified OR gene repertoires in orangutans, marmosets, and mouse lemurs and investigated the birth-and-death process of OR genes in the primate lineage. The results showed that 1) all the primate species studied have no more than 400 intact OR genes, fewer than rodents and canine; 2) Despite the similar number of OR genes in the genome, the makeup of the OR gene repertoires between different primate species is quite different as they had undergone dramatic birth-and-death evolution with extensive gene losses in the lineages leading to current species; 3) Apes and Old World monkey (OWM) have similar fraction of pseudogenes, whereas New World monkey (NWM) have fewer pseudogenes. To measure the selective pressure that had affected the OR gene repertoires in primates, we compared the ratio of nonsynonymous with synonymous substitution rates by using 70 one-to-one orthologous quintets among five primate species. We found that OR genes showed relaxed selective constraints in apes (humans, chimpanzees, and orangutans) than in OWMs (macaques) and NWMs (marmosets). We concluded that OR gene repertoires in primates have evolved in such a way to adapt to their respective living environments. Differential selective constraints might play important role in the primate OR gene evolution in each primate species.

Temporal pattern of loss/persistence of duplicate genes involved in signal transduction and metabolic pathways after teleost-specific genome duplication

Background

Recent genomic studies have revealed a teleost-specific third-round whole genome duplication (3R-WGD) event occurred in a common ancestor of teleost fishes. However, it is unclear how the genes duplicated in this event were lost or persisted during the diversification of teleosts, and therefore, how many of the duplicated genes contribute to the genetic differences among teleosts. This subject is also important for understanding the process of vertebrate evolution through WGD events. We applied a comparative evolutionary approach to this question by focusing on the genes involved in long-term potentiation, taste and olfactory transduction, and the tricarboxylic acid cycle, based on the whole genome sequences of four teleosts; zebrafish, medaka, stickleback, and green spotted puffer fish.

Results

We applied a state-of-the-art method of maximum-likelihood phylogenetic inference and conserved synteny analyses to each of 130 genes involved in the above biological systems of human. These analyses identified 116 orthologous gene groups between teleosts and tetrapods, and 45 pairs of 3R-WGD-derived duplicate genes among them. This suggests that more than half [(45×2)/(116+45)] = 56.5%) of the loci, probably more than ten thousand genes, present in a common ancestor of the four teleosts were still duplicated after the 3R-WGD. The estimated temporal pattern of gene loss suggested that, after the 3R-WGD, many (71/116) of the duplicated genes were rapidly lost during the initial 75 million years (MY), whereas on average more than half (27.3/45) of the duplicated genes remaining in the ancestor of the four teleosts (45/116) have persisted for about 275 MY. The 3R-WGD-derived duplicates that have persisted for a long evolutionary periods of time had significantly larger number of interacting partners and longer length of protein coding sequence, implying that they tend to be more multifunctional than the singletons after the 3R-WGD.

Conclusion

We have shown firstly the temporal pattern of gene loss process after 3R-WGD on the basis of teleost phylogeny and divergence time frameworks. The 3R-WGD-derived duplicates have not undergone constant exponential decay, suggesting that selection favoured the long-term persistence of a subset of duplicates that tend to be multi-functional. On the basis of these results obtained from the analysis of 116 orthologous gene groups, we propose that more than ten thousand of 3R-WGD-derived duplicates have experienced lineage-specific evolution, that is, the differential sub-/neo-functionalization or secondary loss between lineages, and contributed to teleost diversity.

Human olfaction: from genomic variation to phenotypic diversity

The sense of smell is a complex molecular device, encompassing several hundred olfactory receptor proteins (ORs). These receptors, encoded by the largest human gene superfamily, integrate odorant signals into an accurate 'odor image' in the brain. Widespread phenotypic diversity in human olfaction is, in part, attributable to prevalent genetic variation in OR genes, owing to copy number variation, deletion alleles and deleterious single nucleotide polymorphisms. The development of new genomic tools, including next generation sequencing and CNV assays, provides opportunities to characterize the genetic variations of this system. The advent of large-scale functional screens of expressed ORs, combined with genetic association studies, has the potential to link variations in ORs to human chemosensory phenotypes. This promises to provide a genome-wide view of human olfaction, resulting in a deeper understanding of personalized odor coding, with the potential to decipher flavor and fragrance preferences.

Hotspots of biased nucleotide substitutions in human genes

Genes that have experienced accelerated evolutionary rates on the human lineage during recent evolution are candidates for involvement in human-specific adaptations. To determine the forces that cause increased evolutionary rates in certain genes, we analyzed alignments of 10,238 human genes to their orthologues in chimpanzee and macaque. Using a likelihood ratio test, we identified protein-coding sequences with an accelerated rate of base substitutions along the human lineage. Exons evolving at a fast rate in humans have a significant tendency to contain clusters of AT-to-GC (weak-to-strong) biased substitutions. This pattern is also observed in noncoding sequence flanking rapidly evolving exons. Accelerated exons occur in regions with elevated male recombination rates and exhibit an excess of nonsynonymous substitutions relative to the genomic average. We next analyzed genes with significantly elevated ratios of nonsynonymous to synonymous rates of base substitution (d_N/d_S) along the human lineage, and those with an excess of amino acid replacement substitutions relative to human polymorphism. These genes also show evidence of clusters of weak-to-strong biased substitutions. These findings indicate that a recombination-associated process, such as biased gene conversion (BGC), is driving fixation of GC alleles in the human genome. This process can lead to accelerated evolution in coding sequences and excess amino acid replacement substitutions, thereby generating significant results for tests of positive selection.

Regions of the human genome that appear to evolve rapidly may have been under strong positive selection and could contain the genetic changes responsible for the uniqueness of our species. However, neutral (nonadaptive) evolutionary processes can give rise to signals that can be mistaken as signs of selection. In this article, we identify coding sequences that have undergone accelerated rates of change in humans, affecting the divergence of the proteins they encode. By analyzing patterns of molecular evolution in these genes and their distribution in the genome, we show that many protein-coding changes in the fastest-changing genes are not a result of selection operating on the genes, but instead result from biased fixation of AT-to-GC mutations. Our findings are consistent with a model of recombination-driven biased gene conversion. This leads to the provocative hypothesis that many of the genetic changes leading to human-specific characters may have been prompted by fixation of deleterious mutations.

Natural selection is commonly believed to be the main engine of functional genetic change, but a separate neutral evolutionary process linked to recombination may have contributed significantly to the divergence of human proteins.

Human pseudogenes of the ABO family show a complex evolutionary dynamics and loss of function

The GT6 glycosyltransferases gene family, that includes the ABO blood group, shows a complex evolution pattern, with multiple events of gain and loss in different mammal species. In humans the ABO gene is considered the sole functional member although the O allele is null and is fixed in certain populations. Here, we analyze the human GT6 pseudogene sequences (Forssman, IGB3, GGTA1, GT6m5, GT6m6, and GT6m7) from an evolutionary perspective, by the study of (i) their diversity levels in populations through the resequencing analysis of European and African individuals; (ii) the interpopulation differentiation, with genotyping data from a survey of populations covering most of human genetic diversity; and (iii) the interespecific divergence, by the comparison of the human and some other primate species sequences. Since pseudogenes are expected to evolve under neutrality, they should show an evolutionary pattern different to that of functional sequences, with higher levels of diversity as well as a ratio of nonsynonymous to synonymous changes close to 1. We describe some departures from these expectations, including selection for inactivation in IGB3, GGTA1, and the interesting case of FS (Forssman) with a probable shift of its initial function in the primate lineage, which put it apart from a pure neutral pseudogene. These results suggest that some of these GT6 human pseudogenes may still be functional and retain some valuable unknown function in humans, in some case even at the protein level. The evolutionary analysis of all members of the GT6 family in humans allows an insight into their functional history, a process likely due to the interaction of the host glycans that they synthesize with pathogens; the past process that can be unraveled through the footprints left by natural selection in the extant genome variation.

Genome-wide analysis of chimpanzee genes with premature termination codons

Background

Premature termination codons (PTCs) cause mRNA degradation or a truncated protein and thereby contribute to the transcriptome and proteome divergence between species. Here we present the first genome-wide study of PTCs in the chimpanzee. By comparing the human and chimpanzee genome sequences we identify and characterize genes with PTCs, in order to understand the contribution of these mutations to the transcriptome diversity between the species.

Results

We have studied a total of 13,487 human-chimpanzee gene pairs and found that ~8% were affected by PTCs in the chimpanzee. A majority (764/1,109) of PTCs were caused by insertions or deletions and the remaining part was caused by substitutions. The distribution of PTC genes varied between chromosomes, with Y having the highest proportion. Furthermore, the density of PTC genes varied on a megabasepair scale within chromosomes and we found the density to be correlated both with indel divergence and proximity to the telomere. Within genes, PTCs were more common close to the 5' and 3' ends of the amino acid sequence. Gene Ontology classification revealed that olfactory receptor genes were over represented among the PTC genes.

Conclusion

Our results showed that the density of PTC genes fluctuated across the genome depending on the local genomic context. PTCs were preferentially located in the terminal parts of the transcript, which generally have a lower frequency of functional domains, indicating that selection was operating against PTCs at sites central to protein function. The enrichment of GO terms associated with olfaction suggests that PTCs may have influenced the difference in the repertoire of olfactory genes between humans and chimpanzees. In summary, 8% of the chimpanzee genes were affected by PTCs and this type of variation is likely to have an important effect on the transcript and proteomic divergence between humans and chimpanzees.

Modelling the evolution of multi-gene families

A number of biological processes can lead to genes being copied within the genome of some given species. Duplicate genes of this form are called paralogs and such genes share a high degree sequence similarity as well as often having closely related functions. Some genes have become widely duplicated to form multigene families in which the copies are distributed both within the genomes of individual species and across different species. Statistical modelling of gene duplication and the evolution of multi-gene families currently lags behind well-established models of DNA sequence evolution despite an increasing volume of available data, but the analysis of multi-gene families is important as part of a wider effort to understand evolution at the genomic level. This article reviews existing approaches to modelling multi-gene families and presents various challenges and possibilities for this exciting area of research.