Extensive gains and losses of olfactory receptor genes in mammalian evolution

Structural differences and differential expression among rhabdomeric opsins reveal functional change after gene duplication in the bay scallop, Argopecten irradians (Pectinidae)

Background

Opsins are the only class of proteins used for light perception in image-forming eyes. Gene duplication and subsequent functional divergence of opsins have played an important role in expanding photoreceptive capabilities of organisms by altering what wavelengths of light are absorbed by photoreceptors (spectral tuning). However, new opsin copies may also acquire novel function or subdivide ancestral functions through changes to temporal, spatial or the level of gene expression. Here, we test how opsin gene copies diversify in function and evolutionary fate by characterizing four rhabdomeric (G_q-protein coupled) opsins in the scallop, Argopecten irradians, identified from tissue-specific transcriptomes.

Results

Under a phylogenetic analysis, we recovered a pattern consistent with two rounds of duplication that generated the genetic diversity of scallop G_q-opsins. We found strong support for differential expression of paralogous G_q-opsins across ocular and extra-ocular photosensitive tissues, suggesting that scallop G_q-opsins are used in different biological contexts due to molecular alternations outside and within the protein-coding regions. Finally, we used available protein models to predict which amino acid residues interact with the light-absorbing chromophore. Variation in these residues suggests that the four G_q-opsin paralogs absorb different wavelengths of light.

Conclusions

Our results uncover novel genetic and functional diversity in the light-sensing structures of the scallop, demonstrating the complicated nature of G_q-opsin diversification after gene duplication. Our results highlight a change in the nearly ubiquitous shadow response in molluscs to a narrowed functional specificity for visual processes in the eyed scallop. Our findings provide a starting point to study how gene duplication may coincide with eye evolution, and more specifically, different ways neofunctionalization of G_q-opsins may occur.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0823-9) contains supplementary material, which is available to authorized users.

The Genomes of Two Bat Species with Long Constant Frequency Echolocation Calls

Bats can perceive the world by using a wide range of sensory systems, and some of the systems have become highly specialized, such as auditory sensory perception. Among bat species, the Old World leaf-nosed bats and horseshoe bats (rhinolophoid bats) possess the most sophisticated echolocation systems. Here, we reported the whole-genome sequencing and de novo assembles of two rhinolophoid bats – the great leaf-nosed bat (Hipposideros armiger) and the Chinese rufous horseshoe bat (Rhinolophus sinicus). Comparative genomic analyses revealed the adaptation of auditory sensory perception in the rhinolophoid bat lineages, probably resulting from the extreme selectivity used in the auditory processing by these bats. Pseudogenization of some vision-related genes in rhinolophoid bats was observed, suggesting that these genes have undergone relaxed natural selection. An extensive contraction of olfactory receptor gene repertoires was observed in the lineage leading to the common ancestor of bats. Further extensive gene contractions can be observed in the branch leading to the rhinolophoid bats. Such concordance suggested that molecular changes at one sensory gene might have direct consequences for genes controlling for other sensory modalities. To characterize the population genetic structure and patterns of evolution, we re-sequenced the genome of 20 great leaf-nosed bats from four different geographical locations of China. The result showed similar sequence diversity values and little differentiation among populations. Moreover, evidence of genetic adaptations to high altitudes in the great leaf-nosed bats was observed. Taken together, our work provided a useful resource for future research on the evolution of bats.

Identification of Complete Repertoire of Apis florea Odorant Receptors Reveals Complex Orthologous Relationships with Apis mellifera

We developed a computational pipeline for homology based identification of the complete repertoire of olfactory receptor (OR) genes in the Asian honey bee species, Apis florea. Apis florea is phylogenetically the most basal honey bee species and also the most distant sister species to the Western honey bee Apis mellifera, for which all OR genes had been identified before. Using our pipeline, we identified 180 OR genes in A. florea, which is very similar to the number of ORs identified in A. mellifera (177 ORs). Many characteristics of the ORs including gene structure, synteny of tandemly repeated ORs and basic phylogenetic clustering are highly conserved. The composite phylogenetic tree of A. florea and A. mellifera ORs could be divided into 21 clades which are in harmony with the existing Hymenopteran tree. However, we found a few nonorthologous OR relationships between both species as well as independent pseudogenization of ORs suggesting separate evolutionary changes. Particularly, a subgroup of the OR gene clade XI, which had been hypothesized to code cuticular hydrocarbon receptors showed a high number of species-specific ORs. RNAseq analysis detected a total number of 145 OR transcripts in male and 162 in female antennae. Most of the OR genes were highly expressed on the female antennae. However, we detected five distinct male-biased OR genes, out of which three genes (AfOr11, AfOr18, AfOr170P) were shown to be male-biased in A. mellifera, too, thus corroborating a behavioral function in sex-pheromone communication.

Evolution of brain elaboration

Large, complex brains have evolved independently in several lineages of protostomes and deuterostomes. Sensory centres in the brain increase in size and complexity in proportion to the importance of a particular sensory modality, yet often share circuit architecture because of constraints in processing sensory inputs. The selective pressures driving enlargement of higher, integrative brain centres has been more difficult to determine, and may differ across taxa. The capacity for flexible, innovative behaviours, including learning and memory and other cognitive abilities, is commonly observed in animals with large higher brain centres. Other factors, such as social grouping and interaction, appear to be important in a more limited range of taxa, while the importance of spatial learning may be a common feature in insects with large higher brain centres. Despite differences in the exact behaviours under selection, evolutionary increases in brain size tend to derive from common modifications in development and generate common architectural features, even when comparing widely divergent groups such as vertebrates and insects. These similarities may in part be influenced by the deep homology of the brains of all Bilateria, in which shared patterns of developmental gene expression give rise to positionally, and perhaps functionally, homologous domains. Other shared modifications of development appear to be the result of homoplasy, such as the repeated, independent expansion of neuroblast numbers through changes in genes regulating cell division. The common features of large brains in so many groups of animals suggest that given their common ancestry, a limited set of mechanisms exist for increasing structural and functional diversity, resulting in many instances of homoplasy in bilaterian nervous systems.

Penguins reduced olfactory receptor genes common to other waterbirds

The sense of smell, or olfaction, is fundamental in the life of animals. However, penguins (Aves: Sphenisciformes) possess relatively small olfactory bulbs compared with most other waterbirds such as Procellariiformes and Gaviiformes. To test whether penguins have a reduced reliance on olfaction, we analyzed the draft genome sequences of the two penguins, which diverged at the origin of the order Sphenisciformes; we also examined six closely related species with available genomes, and identified 29 one-to-one orthologous olfactory receptor genes (i.e. ORs) that are putatively functionally conserved and important across the eight birds. To survey the 29 one-to-one orthologous ORs in penguins and their relatives, we newly generated 34 sequences that are missing from the draft genomes. Through the analysis of totaling 378 OR sequences, we found that, of these functionally important ORs common to other waterbirds, penguins have a significantly greater percentage of OR pseudogenes than other waterbirds, suggesting a reduction of olfactory capability. The penguin-specific reduction of olfactory capability arose in the common ancestor of penguins between 23 and 60 Ma, which may have resulted from the aquatic specializations for underwater vision. Our study provides genetic evidence for a possible reduction of reliance on olfaction in penguins.

Combinatorial effects of odorants on mouse behavior

The mechanisms by which odors induce instinctive behaviors are largely unknown. Odor detection in the mouse nose is mediated by >1, 000 different odorant receptors (ORs) and trace amine-associated receptors (TAARs). Odor perceptions are encoded combinatorially by ORs and can be altered by slight changes in the combination of activated receptors. However, the stereotyped nature of instinctive odor responses suggests the involvement of specific receptors and genetically programmed neural circuits relatively immune to extraneous odor stimuli and receptor inputs. Here, we report that, contrary to expectation, innate odor-induced behaviors can be context-dependent. First, different ligands for a given TAAR can vary in behavioral effect. Second, when combined, some attractive and aversive odorants neutralize one another's behavioral effects. Both a TAAR ligand and a common odorant block aversion to a predator odor, indicating that this ability is not unique to TAARs and can extend to an aversive response of potential importance to survival. In vitro testing of single receptors with binary odorant mixtures indicates that behavioral blocking can occur without receptor antagonism in the nose. Moreover, genetic ablation of a single receptor prevents its cognate ligand from blocking predator odor aversion, indicating that the blocking requires sensory input from the receptor. Together, these findings indicate that innate odor-induced behaviors can depend on context, that signals from a single receptor can block innate odor aversion, and that instinctive behavioral responses to odors can be modulated by interactions in the brain among signals derived from different receptors.

Molecular Evolution and Functional Divergence of Trace Amine-Associated Receptors

Trace amine-associated receptors (TAARs) are a member of the G-protein-coupled receptor superfamily and are known to be expressed in olfactory sensory neurons. A limited number of molecular evolutionary studies have been done for TAARs so far. To elucidate how lineage-specific evolution contributed to their functional divergence, we examined 30 metazoan genomes. In total, 493 TAAR gene candidates (including 84 pseudogenes) were identified from 26 vertebrate genomes. TAARs were not identified from non-vertebrate genomes. An ancestral-type TAAR-like gene appeared to have emerged in lamprey. We found four therian-specific TAAR subfamilies (one eutherian-specific and three metatherian-specific) in addition to previously known nine subfamilies. Many species-specific TAAR gene duplications and losses contributed to a large variation of TAAR gene numbers among mammals, ranging from 0 in dolphin to 26 in flying fox. TAARs are classified into two groups based on binding preferences for primary or tertiary amines as well as their sequence similarities. Primary amine-detecting TAARs (TAAR1-4) have emerged earlier, generally have single-copy orthologs (very few duplication or loss), and have evolved under strong functional constraints. In contrast, tertiary amine-detecting TAARs (TAAR5-9) have emerged more recently and the majority of them experienced higher rates of gene duplications. Protein members that belong to the tertiary amine-detecting TAAR group also showed the patterns of positive selection especially in the area surrounding the ligand-binding pocket, which could have affected ligand-binding activities and specificities. Expansions of the tertiary amine-detecting TAAR gene family may have played important roles in terrestrial adaptations of therian mammals. Molecular evolution of the TAAR gene family appears to be governed by a complex, species-specific, interplay between environmental and evolutionary factors.

DoOR 2.0--Comprehensive Mapping of Drosophila melanogaster Odorant Responses

Odors elicit complex patterns of activated olfactory sensory neurons. Knowing the complete olfactome, i.e. the responses in all sensory neurons for all relevant odorants, is desirable to understand olfactory coding. The DoOR project combines all available Drosophila odorant response data into a single consensus response matrix. Since its first release many studies were published: receptors were deorphanized and several response profiles were expanded. In this study, we add unpublished data to the odor-response profiles for four odorant receptors (Or10a, Or42b, Or47b, Or56a). We deorphanize Or69a, showing a broad response spectrum with the best ligands including 3-hydroxyhexanoate, alpha-terpineol, 3-octanol and linalool. We include all of these datasets into DoOR, provide a comprehensive update of both code and data, and new tools for data analyses and visualizations. The DoOR project has a web interface for quick queries (http://neuro.uni.kn/DoOR), and a downloadable, open source toolbox written in R, including all processed and original datasets. DoOR now gives reliable odorant-responses for nearly all Drosophila olfactory responding units, listing 693 odorants, for a total of 7381 data points.

Contrasting Patterns of Evolutionary Diversification in the Olfactory Repertoires of Reptile and Bird Genomes

Olfactory receptors (ORs) are membrane proteins that mediate the detection of odorants in the environment, and are the largest vertebrate gene family. Comparative studies of mammalian genomes indicate that OR repertoires vary widely, even between closely related lineages, as a consequence of frequent OR gains and losses. Several studies also suggest that mammalian OR repertoires are influenced by life history traits. Sauropsida is a diverse group of vertebrates group that is the sister group to mammals, and includes birds, testudines, squamates, and crocodilians, and represents a natural system to explore predictions derived from mammalian studies. In this study, we analyzed olfactory receptor (OR) repertoire variation among several representative species and found that the number of intact OR genes in sauropsid genomes analyzed ranged over an order of magnitude, from 108 in the green anole to over 1,000 in turtles. Our results suggest that different sauropsid lineages have highly divergent OR repertoire composition that derive from lineage-specific combinations of gene expansions, losses, and retentions of ancestral OR genes. These differences also suggest that varying degrees of adaption related to life history have shaped the unique OR repertoires observed across sauropsid lineages.

Tetrapod V1R-like ora genes in an early-diverging ray-finned fish species: the canonical six ora gene repertoire of teleost fish resulted from gene loss in a larger ancestral repertoire

BACKGROUND

Chemical senses serve a multitude of essential functions across the animal kingdom. Vertebrates employ four GPCR families to detect odors, among them the v1r/ora gene family. The V1R family is known to evolve rapidly in the lobe-finned lineage giving rise to tetrapods, but the homologous ORA family consists of just six highly conserved genes in teleost fish, with direct orthologs in the lobe-finned fish coelacanth. Thus, the teleost repertoire of six canonical ora genes was assumed to be the ancestral feature before the divergence of ray-finned and lobe-finned fish. So far, this hypothesis has not been tested with earlier diverging ray-finned fish.

RESULTS

We have newly identified the complete ora gene repertoires of five teleost species, and of spotted gar, a basal ray-finned fish, using thorough data mining and extensive phylogenetic analysis. The genomes of eight further teleost species were re-analyzed for their ORA repertoires. We report that direct orthologs of the six canonical ora genes (ora1-6) were present in all newly analyzed species, with faithfully preserved exon/intron structure and mostly preserved genomic arrangement in symmetric pairs for ora1-4. In four teleost species including medaka and cave fish we observe species-specific gene duplication events. Thus, the ora gene repertoire in teleost fish is not quite as strictly conserved as previously assumed. In fact, the examination of non-synonymous vs. synonymous substitution rates (dN/dS) shows pronounced negative selection in five of the six ora genes, but also rare occurrence of positive selection in ora3 and ora6. Surprisingly, spotted gar possesses beyond the six canonical genes three additional genes, ora7-8b, orthologous to coelacanth genes v1r07-10. No orthologs for these genes were found in teleosts and cartilaginous fish.

CONCLUSIONS

Early diverging ray-finned fish such as the spotted gar possess several v1r-like genes previously assumed to be restricted to the lobe-finned lineage, but now found to be already present in the most recent common ancestor of lobe- and ray-finned fish. Thus, the presence of just six canonical ora genes in many teleost species is not the ancestral feature of the ray-finned lineage, but caused by loss of two ancestral genes in teleosts.

The Mouse Solitary Odorant Receptor Gene Promoters as Models for the Study of Odorant Receptor Gene Choice

BACKGROUND

In vertebrates, several anatomical regions located within the nasal cavity mediate olfaction. Among these, the main olfactory epithelium detects most conventional odorants. Olfactory sensory neurons, provided with cilia exposed to the air, detect volatile chemicals via an extremely large family of seven-transmembrane chemoreceptors named odorant receptors. Their genes are expressed in a monogenic and monoallelic fashion: a single allele of a single odorant receptor gene is transcribed in a given mature neuron, through a still uncharacterized molecular mechanism known as odorant receptor gene choice.

AIM

Odorant receptor genes are typically arranged in genomic clusters, but a few are isolated (we call them solitary) from the others within a region broader than 1 Mb upstream and downstream with respect to their transcript's coordinates. The study of clustered genes is problematic, because of redundancy and ambiguities in their regulatory elements: we propose to use the solitary genes as simplified models to understand odorant receptor gene choice.

PROCEDURES

Here we define number and identity of the solitary genes in the mouse genome (C57BL/6J), and assess the conservation of the solitary status in some mammalian orthologs. Furthermore, we locate their putative promoters, predict their homeodomain binding sites (commonly present in the promoters of odorant receptor genes) and compare candidate promoter sequences with those of wild-caught mice. We also provide expression data from histological sections.

RESULTS

In the mouse genome there are eight intact solitary genes: Olfr19 (M12), Olfr49, Olfr266, Olfr267, Olfr370, Olfr371, Olfr466, Olfr1402; five are conserved as solitary in rat. These genes are all expressed in the main olfactory epithelium of three-day-old mice. The C57BL/6J candidate promoter of Olfr370 has considerably varied compared to its wild-type counterpart. Within the putative promoter for Olfr266 a homeodomain binding site is predicted. As a whole, our findings favor Olfr266 as a model gene to investigate odorant receptor gene choice.

Evaluation of Giant African Pouched Rats for Detection of Pulmonary Tuberculosis in Patients from a High-Endemic Setting

Background

This study established evidence about the diagnostic performance of trained giant African pouched rats for detecting Mycobacterium tuberculosis in sputum of well-characterised patients with presumptive tuberculosis (TB) in a high-burden setting.

Methods

The TB detection rats were evaluated using sputum samples of patients with presumptive TB enrolled in two prospective cohort studies in Bagamoyo, Tanzania. The patients were characterised by sputum smear microscopy and culture, including subsequent antigen or molecular confirmation of Mycobacterium tuberculosis, and by clinical data at enrolment and for at least 5-months of follow-up to determine the reference standard. Seven trained giant African pouched rats were used for the detection of TB in the sputum samples after shipment to the APOPO project in Morogoro, Tanzania.

Results

Of 469 eligible patients, 109 (23.2%) were culture-positive for Mycobacterium tuberculosis and 128 (27.3%) were non-TB controls with sustained recovery after 5 months without anti-TB treatment. The HIV prevalence was 46%. The area under the receiver operating characteristic curve of the seven rats for the detection of culture-positive pulmonary tuberculosis was 0.72 (95% CI 0.66–0.78). An optimal threshold could be defined at ≥2 indications by rats in either sample with a corresponding sensitivity of 56.9% (95% CI 47.0–66.3), specificity of 80.5% (95% CI 72.5–86.9), positive and negative predictive value of 71.3% (95% CI 60.6–80.5) and 68.7% (95% CI 60.6–76.0), and an accuracy for TB diagnosis of 69.6%. The diagnostic performance was negatively influenced by low burden of bacilli, and independent of the HIV status.

Conclusion

Giant African pouched rats have potential for detection of tuberculosis in sputum samples. However, the diagnostic performance characteristics of TB detection rats do not currently meet the requirements for high-priority, rapid sputum-based TB diagnostics as defined by the World Health Organization.

Gene-pseudogene evolution: a probabilistic approach

Over the last decade, methods have been developed for the reconstruction of gene trees that take into account the species tree. Many of these methods have been based on the probabilistic duplication-loss model, which describes how a gene-tree evolves over a species-tree with respect to duplication and losses, as well as extension of this model, e.g., the DLRS (Duplication, Loss, Rate and Sequence evolution) model that also includes sequence evolution under relaxed molecular clock. A disjoint, almost as recent, and very important line of research has been focused on non protein-coding, but yet, functional DNA. For instance, DNA sequences being pseudogenes in the sense that they are not translated, may still be transcribed and the thereby produced RNA may be functional.

We extend the DLRS model by including pseudogenization events and devise an MCMC framework for analyzing extended gene families consisting of genes and pseudogenes with respect to this model, i.e., reconstructing gene-trees and identifying pseudogenization events in the reconstructed gene-trees. By applying the MCMC framework to biologically realistic synthetic data, we show that gene-trees as well as pseudogenization points can be inferred well. We also apply our MCMC framework to extended gene families belonging to the Olfactory Receptor and Zinc Finger superfamilies. The analysis indicate that both these super families contains very old pseudogenes, perhaps so old that it is reasonable to suspect that some are functional. In our analysis, the sub families of the Olfactory Receptors contains only lineage specific pseudogenes, while the sub families of the Zinc Fingers contains pseudogene lineages common to several species.

Exome sequencing of a colorectal cancer family reveals shared mutation pattern and predisposition circuitry along tumor pathways

The molecular basis of cancer and cancer multiple phenotypes are not yet fully understood. Next Generation Sequencing promises new insight into the role of genetic interactions in shaping the complexity of cancer. Aiming to outline the differences in mutation patterns between familial colorectal cancer cases and controls we analyzed whole exomes of cancer tissues and control samples from an extended colorectal cancer pedigree, providing one of the first data sets of exome sequencing of cancer in an African population against a background of large effective size typically with excess of variants. Tumors showed hMSH2 loss of function SNV consistent with Lynch syndrome. Sets of genes harboring insertions-deletions in tumor tissues revealed, however, significant GO enrichment, a feature that was not seen in control samples, suggesting that ordered insertions-deletions are central to tumorigenesis in this type of cancer. Network analysis identified multiple hub genes of centrality. ELAVL1/HuR showed remarkable centrality, interacting specially with genes harboring non-synonymous SNVs thus reinforcing the proposition of targeted mutagenesis in cancer pathways. A likely explanation to such mutation pattern is DNA/RNA editing, suggested here by nucleotide transition-to-transversion ratio that significantly departed from expected values (p-value 5e-6). NFKB1 also showed significant centrality along with ELAVL1, raising the suspicion of viral etiology given the known interaction between oncogenic viruses and these proteins.

Baton pass hypothesis: successive incorporation of unconserved endogenous retroviral genes for placentation during mammalian evolution

It is well accepted that numerous RNAs derived from endogenous retroviruses (ERVs) are expressed in mammalian reproductive structures, particularly in the uterus, trophoblast, and placenta. Syncytin 1 and syncytin 2 in humans and syncytin A and syncytin B in mice are membrane proteins originating from Env genes of ERVs. These ERVs are involved in the fusion of trophoblast cells, resulting in multinucleated syncytiotrophoblast formation. Evidence accumulated indicates that syncytin-like fusogenic proteins are expressed in the placenta of rabbits, dogs/cats, ruminant ungulates, tenrecs, and opossums. The syncytin genes so far characterized are known to be endogenized to the host genome only within the past 12-80 million years, more recently than the appearance of mammalian placentas, estimated to be 160-180 million years ago. We speculate that ERVs including syncytin-like gene variants integrated into mammalian genomes in a locus-specific manner have replaced the genes previously responsible for cell fusion. We therefore propose the 'baton pass' hypothesis, in which multiple successive ERV variants 'take over' cell-fusion roles, resulting in increased trophoblast cell fusion, morphological variations in placental structures, and enhanced reproductive success in placental mammals.

Large-scale transcriptional profiling of chemosensory neurons identifies receptor-ligand pairs in vivo

In mammals, olfactory perception is based on the combinatorial activation of G protein-coupled receptors. Identifying the full repertoire of receptors activated by a given odorant in vivo, a quest that has been hampered for over 20 years by technical difficulties, would represent an important step in deciphering the rules governing chemoperception. We found that odorants induced a fast and reversible concentration-dependent decrease in the transcription of genes corresponding to activated receptors in intact mice. On the basis of this finding, we developed a large-scale transcriptomic approach to uncover receptor-ligand pairs in vivo. We identified the mouse and rat odorant receptor signatures corresponding to specific odorants. Finally, we found that this approach, which can be used for species for which no genomic sequence is available, is also applicable to non-vertebrate species such as Drosophila.

Genomic Variants Revealed by Invariably Missing Genotypes in Nelore Cattle

High density genotyping panels have been used in a wide range of applications. From population genetics to genome-wide association studies, this technology still offers the lowest cost and the most consistent solution for generating SNP data. However, in spite of the application, part of the generated data is always discarded from final datasets based on quality control criteria used to remove unreliable markers. Some discarded data consists of markers that failed to generate genotypes, labeled as missing genotypes. A subset of missing genotypes that occur in the whole population under study may be caused by technical issues but can also be explained by the presence of genomic variations that are in the vicinity of the assayed SNP and that prevent genotyping probes from annealing. The latter case may contain relevant information because these missing genotypes might be used to identify population-specific genomic variants. In order to assess which case is more prevalent, we used Illumina HD Bovine chip genotypes from 1,709 Nelore (Bos indicus) samples. We found 3,200 missing genotypes among the whole population. NGS re-sequencing data from 8 sires were used to verify the presence of genomic variations within their flanking regions in 81.56% of these missing genotypes. Furthermore, we discovered 3,300 novel SNPs/Indels, 31% of which are located in genes that may affect traits of importance for the genetic improvement of cattle production.

Repertoire and evolution of TNF superfamily in Crassostrea gigas: implications for expansion and diversification of this superfamily in Mollusca

Tumor necrosis factor superfamily (TNFSF) members represent a group of cytokines participating in diverse immunological, pathological and developmental pathways. However, compared with deuterostomia and cnidaia, the composition and evolution of TNF homologous in protostomia are still not well understood. In the present study, a total of 81 TNF superfamily (TNFSF) genes from 15 mollusk species, including 23 TNFSF genes from Crassostrea gigas, were surveyed by genome-wide bioinformatics analysis. The phylogenetic analysis showed that 14 out of 23 C. gigas TNFSF genes in five clades exhibited orthologous relationships with Pinctada fucata TNFSF genes. Moreover, there were 15 C. gigas TNFSF genes located in oyster-specific clusters, which were contributed by small-scaled tandem and/or segmental duplication events in oyster. By comparing the sequences of duplicated TNFSF pairs, exon loss and variant in exon/intron length were revealed as the major modes of divergence in gene structure. Most of the duplicated C. gigas TNFSF pairs were evolved under purifying selection with consistent tissue expression patterns, implying functional constraint shaped diversification. This study demonstrated the expansion and early divergence of TNF superfamily in C. gigas, which provides potential insight into revealing the evolution and function of this superfamily in mollusk.

Kiwi genome provides insights into evolution of a nocturnal lifestyle

BACKGROUND

Kiwi, comprising five species from the genus Apteryx, are endangered, ground-dwelling bird species endemic to New Zealand. They are the smallest and only nocturnal representatives of the ratites. The timing of kiwi adaptation to a nocturnal niche and the genomic innovations, which shaped sensory systems and morphology to allow this adaptation, are not yet fully understood.

RESULTS

We sequenced and assembled the brown kiwi genome to 150-fold coverage and annotated the genome using kiwi transcript data and non-redundant protein information from multiple bird species. We identified evolutionary sequence changes that underlie adaptation to nocturnality and estimated the onset time of these adaptations. Several opsin genes involved in color vision are inactivated in the kiwi. We date this inactivation to the Oligocene epoch, likely after the arrival of the ancestor of modern kiwi in New Zealand. Genome comparisons between kiwi and representatives of ratites, Galloanserae, and Neoaves, including nocturnal and song birds, show diversification of kiwi's odorant receptors repertoire, which may reflect an increased reliance on olfaction rather than sight during foraging. Further, there is an enrichment of genes influencing mitochondrial function and energy expenditure among genes that are rapidly evolving specifically on the kiwi branch, which may also be linked to its nocturnal lifestyle.

CONCLUSIONS

The genomic changes in kiwi vision and olfaction are consistent with changes that are hypothesized to occur during adaptation to nocturnal lifestyle in mammals. The kiwi genome provides a valuable genomic resource for future genome-wide comparative analyses to other extinct and extant diurnal ratites.

Quantifying the cribriform plate: influences of allometry, function, and phylogeny in Carnivora

The small, perforated bony cup of the anterior cranial fossa called the cribriform plate (CP) is perhaps the best-preserved remnant of olfactory anatomy in fossil mammal skulls. The CP and its myriad foramina record the passage of peripheral olfactory nerves from nasal cavity to olfactory bulb. Previous work has suggested that CP surface area reflects aspects of olfactory capacity (as inferred from habitat and observed behavior) in mammals. To further explore the utility of CP as a proxy for olfactory function, we designed novel, nondestructive digital methods to quantify CP morphology from dry skulls. Using CT scans and 3-D imaging software, we quantified CP features from 42 species of Carnivora, a group that represents a wide spectrum of ecologies and sensory demands. Two metrics, CP surface area (CPSA) and cumulative CP foramina area (FXSA), scaled to skull length with negative allometry, and differed between aquatic and terrestrial species, with the former having reduced areas. Number of foramina (NF) was not correlated with skull length but tended to be greater in caniforms than feliforms. Both CPSA and FXSA are well correlated with ethmoturbinal surface area, a known osteological correlate of olfactory function. This suggests that CPSA and FXSA are useful proxies for olfactory ability, especially when studying fossils or skulls in which turbinals are not preserved. Total area of CP foramina (FXSA), an exacting measure of olfactory nerve endocasts, is tightly correlated with CPSA. Because of this, it may be desirable to use CPSA alone as a proxy given that it is easier to measure than FXSA.

Organization and distribution of glomeruli in the bowhead whale olfactory bulb

Although modern baleen whales (Mysticeti) retain a functional olfactory system that includes olfactory bulbs, cranial nerve I and olfactory receptor genes, their olfactory capabilities have been reduced to a great degree. This reduction likely occurred as a selective response to their fully aquatic lifestyle. The glomeruli that occur in the olfactory bulb can be divided into two non-overlapping domains, a dorsal domain and a ventral domain. Recent molecular studies revealed that all modern whales have lost olfactory receptor genes and marker genes that are specific to the dorsal domain. Here we show that olfactory bulbs of bowhead whales (Balaena mysticetus) lack glomeruli on the dorsal side, consistent with the molecular data. In addition, we estimate that there are more than 4,000 glomeruli elsewhere in the bowhead whale olfactory bulb, which is surprising given that bowhead whales possess only 80 intact olfactory receptor genes. Olfactory sensory neurons that express the same olfactory receptors in rodents generally project to two specific glomeruli in an olfactory bulb, implying an approximate 1:2 ratio of the number of olfactory receptors to the number of glomeruli. Here we show that this ratio does not apply to bowhead whales, reiterating the conceptual limits of using rodents as model organisms for understanding the initial coding of odor information among mammals.

Identification of novel post-transcriptional features in olfactory receptor family mRNAs

Olfactory receptor (Olfr) genes comprise the largest gene family in mice. Despite their importance in olfaction, how most Olfr mRNAs are regulated remains unexplored. Using RNA-seq analysis coupled with analysis of pre-existing databases, we found that Olfr mRNAs have several atypical features suggesting that post-transcriptional regulation impacts their expression. First, Olfr mRNAs, as a group, have dramatically higher average AU-content and lower predicted secondary structure than do control mRNAs. Second, Olfr mRNAs have a higher density of AU-rich elements (AREs) in their 3'UTR and upstream open reading frames (uORFs) in their 5 UTR than do control mRNAs. Third, Olfr mRNAs have shorter 3' UTR regions and with fewer predicted miRNA-binding sites. All of these novel properties correlated with higher Olfr expression. We also identified striking differences in the post-transcriptional features of the mRNAs from the two major classes of Olfr genes, a finding consistent with their independent evolutionary origin. Together, our results suggest that the Olfr gene family has encountered unusual selective forces in neural cells that have driven them to acquire unique post-transcriptional regulatory features. In support of this possibility, we found that while Olfr mRNAs are degraded by a deadenylation-dependent mechanism, they are largely protected from this decay in neural lineage cells.

Ion transporter NKCC1, modulator of neurogenesis in murine olfactory neurons

Olfaction is one of the most crucial senses for vertebrates regarding foraging and social behavior. Therefore, it is of particular interest to investigate the sense of smell, its function on a molecular level, the signaling proteins involved in the process and the mechanism of required ion transport. In recent years, the precise role of the ion transporter NKCC1 in olfactory sensory neuron (OSN) chloride accumulation has been a controversial subject. NKCC1 is expressed in OSNs and is involved in chloride accumulation of dissociated neurons, but it had not been shown to play a role in mouse odorant sensation. Here, we present electro-olfactogram recordings (EOG) demonstrating that NKCC1-deficient mice exhibit significant defects in perception of a complex odorant mixture (Henkel100) in both air-phase and submerged approaches. Using next generation sequencing (NGS) and RT-PCR experiments of NKCC1-deficient and wild type mouse transcriptomes, we confirmed the absence of a highly expressed ion transporter that could compensate for NKCC1. Additional histological investigations demonstrated a reduced number of cells in the olfactory epithelium (OE), resulting in a thinner neuronal layer. Therefore, we conclude that NKCC1 is an important transporter involved in chloride ion accumulation in the olfactory epithelium, but it is also involved in OSN neurogenesis.

Aquatic adaptation and the evolution of smell and taste in whales

INTRODUCTION

While olfaction is one of the most important senses in most terrestrial mammals, it is absent in modern toothed whales (Odontoceti, Cetacea). Furthermore, behavioral evidence suggests that gustation is very limited. In contrast, their aquatic sistergroup, baleen whales (Mysticeti) retain small but functional olfactory organs, and nothing is known about their gustation. It is difficult to investigate mysticete chemosensory abilities because experiments in a controlled setting are impossible.

RESULTS

Here, we use the functional regionalization of the olfactory bulb (OB) to identify the loss of specific olfactory functions in mysticetes. We provide the whole-genome sequence of a mysticete and show that mysticetes lack the dorsal domain of the OB, an area known to induce innate avoidance behavior against odors of predators and spoiled foods. Genomic and fossil data suggest that mysticetes lost the dorsal domain of the OB before the Odontoceti-Mysticeti split. Furthermore, we found that all modern cetaceans are revealed to have lost the functional taste receptors.

CONCLUSION

These results strongly indicate that profound changes in the chemosensory capabilities had occurred in the cetacean lineage during the period when ancestral whales migrated from land to water.

Deep sequencing of the murine olfactory receptor neuron transcriptome

The ability of animals to sense and differentiate among thousands of odorants relies on a large set of olfactory receptors (OR) and a multitude of accessory proteins within the olfactory epithelium (OE). ORs and related signaling mechanisms have been the subject of intensive studies over the past years, but our knowledge regarding olfactory processing remains limited. The recent development of next generation sequencing (NGS) techniques encouraged us to assess the transcriptome of the murine OE. We analyzed RNA from OEs of female and male adult mice and from fluorescence-activated cell sorting (FACS)-sorted olfactory receptor neurons (ORNs) obtained from transgenic OMP-GFP mice. The Illumina RNA-Seq protocol was utilized to generate up to 86 million reads per transcriptome. In OE samples, nearly all OR and trace amine-associated receptor (TAAR) genes involved in the perception of volatile amines were detectably expressed. Other genes known to participate in olfactory signaling pathways were among the 200 genes with the highest expression levels in the OE. To identify OE-specific genes, we compared olfactory neuron expression profiles with RNA-Seq transcriptome data from different murine tissues. By analyzing different transcript classes, we detected the expression of non-olfactory GPCRs in ORNs and established an expression ranking for GPCRs detected in the OE. We also identified other previously undescribed membrane proteins as potential new players in olfaction. The quantitative and comprehensive transcriptome data provide a virtually complete catalogue of genes expressed in the OE and present a useful tool to uncover candidate genes involved in, for example, olfactory signaling, OR trafficking and recycling, and proliferation.

Olfactory receptor patterning in a higher primate

The mammalian olfactory system detects a plethora of environmental chemicals that are perceived as odors or stimulate instinctive behaviors. Studies using odorant receptor (OR) genes have provided insight into the molecular and organizational strategies underlying olfaction in mice. One important unanswered question, however, is whether these strategies are conserved in primates. To explore this question, we examined the macaque, a higher primate phylogenetically close to humans. Here we report that the organization of sensory inputs in the macaque nose resembles that in mouse in some respects, but not others. As in mouse, neurons with different ORs are interspersed in the macaque nose, and there are spatial zones that differ in their complement of ORs and extend axons to different domains in the olfactory bulb of the brain. However, whereas the mouse has multiple discrete band-like zones, the macaque appears to have only two broad zones. It is unclear whether the organization of OR inputs in a rodent/primate common ancestor degenerated in primates or, alternatively became more sophisticated in rodents. The mouse nose has an additional small family of chemosensory receptors, called trace amine-associated receptors (TAARs), which may detect social cues. Here we find that TAARs are also expressed in the macaque nose, suggesting that TAARs may also play a role in human olfactory perception. We further find that one human TAAR responds to rotten fish, suggesting a possible role as a sentinel to discourage ingestion of food harboring pathogenic microorganisms.

Three crocodilian genomes reveal ancestral patterns of evolution among archosaurs

To provide context for the diversification of archosaurs--the group that includes crocodilians, dinosaurs, and birds--we generated draft genomes of three crocodilians: Alligator mississippiensis (the American alligator), Crocodylus porosus (the saltwater crocodile), and Gavialis gangeticus (the Indian gharial). We observed an exceptionally slow rate of genome evolution within crocodilians at all levels, including nucleotide substitutions, indels, transposable element content and movement, gene family evolution, and chromosomal synteny. When placed within the context of related taxa including birds and turtles, this suggests that the common ancestor of all of these taxa also exhibited slow genome evolution and that the comparatively rapid evolution is derived in birds. The data also provided the opportunity to analyze heterozygosity in crocodilians, which indicates a likely reduction in population size for all three taxa through the Pleistocene. Finally, these data combined with newly published bird genomes allowed us to reconstruct the partial genome of the common ancestor of archosaurs, thereby providing a tool to investigate the genetic starting material of crocodilians, birds, and dinosaurs.

Increased rate of hair keratin gene loss in the cetacean lineage

BACKGROUND

Hair represents an evolutionary innovation that appeared early on mammalian evolutionary history, and presumably contributed significantly to the rapid radiation of the group. An interesting event in hair evolution has been its secondary loss in some mammalian groups, such as cetaceans, whose hairless phenotype appears to be an adaptive response to better meet the environmental conditions. To determine whether different repertoire of keratin genes among mammals can potentially explain the phenotypic hair features of different lineages, we characterized the type I and II clusters of alpha keratins from eight mammalian species, including the hairless dolphin and minke whale representing the order Cetacea.

RESULTS

We combined the available genomic information with phylogenetic analysis to conduct a comprehensive analysis of the evolutionary patterns of keratin gene clusters. We found that both type I and II gene clusters are fairly conserved among the terrestrial mammals included in this study, with lineage specific gene duplication and gene loss. Nevertheless, there is also evidence for an increased rate of pseudogenization in the cetacean lineage when compared to their terrestrial relatives, especially among the hair type keratins.

CONCLUSIONS

Here we present a comprehensive characterization of alpha-keratin genes among mammals and elucidate the mechanisms involved in the evolution of this gene family. We identified lineage-specific gene duplications and gene loss among the Laurasiatherian and Euarchontoglires species included in the study. Interestingly, cetaceans present an increased loss of hair-type keratin genes when compared to other terrestrial mammals. As suggested by the 'less-is-more' hypothesis, we do not rule out the possibility that the gene loss of hair-type keratin genes in these species might be associated to the hairless phenotype and could have been adaptive in response to new selective pressures imposed by the colonization of a new habitat. Our study provides support for the idea that pseudogenes are not simply 'genomic fossils' but instead have adaptive roles during the evolutionary process.

Extreme expansion of the olfactory receptor gene repertoire in African elephants and evolutionary dynamics of orthologous gene groups in 13 placental mammals

Olfactory receptors (ORs) detect odors in the environment, and OR genes constitute the largest multigene family in mammals. Numbers of OR genes vary greatly among species—reflecting the respective species' lifestyles—and this variation is caused by frequent gene gains and losses during evolution. However, whether the extent of gene gains/losses varies among individual gene lineages and what might generate such variation is unknown. To answer these questions, we used a newly developed phylogeny-based method to classify >10,000 intact OR genes from 13 placental mammal species into 781 orthologous gene groups (OGGs); we then compared the OGGs. Interestingly, African elephants had a surprisingly large repertoire (∼2000) of functional OR genes encoded in enlarged gene clusters. Additionally, OR gene lineages that experienced more gene duplication had weaker purifying selection, and Class II OR genes have evolved more dynamically than those in Class I. Some OGGs were highly expanded in a lineage-specific manner, while only three OGGs showed complete one-to-one orthology among the 13 species without any gene gains/losses. These three OGGs also exhibited highly conserved amino acid sequences; therefore, ORs in these OGGs may have physiologically important functions common to every placental mammal. This study provides a basis for inferring OR functions from evolutionary trajectory.

Inventory of the cichlid olfactory receptor gene repertoires: identification of olfactory genes with more than one coding exon

BACKGROUND

To help understand the molecular mechanisms underlying the remarkable phenotypic diversity displayed by cichlids, the genome sequences of O. niloticus, P. nyererei, H. burtoni, N. brichardi and M. zebra were recently determined. Here, we present the contents of the olfactory receptor (OR) repertoires in the genomes of these five fishes.

RESULTS

We performed an exhaustive TBLASTN search of the five cichlid genomes to identify their OR repertoires as completely as possible. We used as bait a set of ORs described in the literature. The cichlid repertoires thereby extracted contained large numbers of complete genes (O. niloticus 158; H. burtoni 90; M. zebra 102; N. brichardi 69; P. nyererei 88), a small numbers of pseudogenes and many "edge genes" corresponding to incomplete genes located at the ends of contigs. A phylogenetic tree was constructed and showed these repertoires include a large number of families and subfamilies. It also allowed the identification of a large number of OR analogues between cichlids with very high amino-acid identity (≥ 99%). Nearly 9% of the full-length cichlid OR genes are composed of several coding exons. This is very unusual for vertebrate OR genes. Nevertheless, the evidence is strong, and includes the donor and acceptor splice junction sequences; also, the positions of these genes in the phylogenetic tree indicate that they constitute subfamilies well apart from non-OR G protein-coupled receptor families.

CONCLUSIONS

Cichlid OR repertoires are made up of a larger number of genes and fewer pseudogenes than those in other teleosts except zebrafish. These ORs share all identified properties common to all fish ORs; however, the large number of families and subfamilies, each containing few ORs implies that they have evolved more rapidly. This high level of OR diversity is consistent with the substantial phenotypic diversity that characterizes cichlids.

Not so pseudo anymore: pseudogenes as therapeutic targets

Pseudogenes are junk DNA gene remnants generated by inactivating mutations or the loss of regulatory sequences, often following gene duplication or retrotransposition events. These pseudogenes have previously been considered to be molecular fossils derived from once-coding genes. In many cases, pseudogenes confer no observable selective advantage to the host organism and may be on a path towards removal from the genome. However, pseudogenes can also serve as raw material for the exaptation of novel functions, particularly in relation to the regulation of gene expression. Many pseudogenes are resurrected as noncoding RNA genes, which function in RNA-based gene regulatory circuits. As such, functional pseudogenes might simply be considered as 'genes'. Here, we discuss the role of these pseudogene-derived RNAs as regulators of gene expression in the context of human disease. In particular, we consider the manipulation of pseudogene transcripts through the use of antisense oligonucleotides, siRNAs, aptamers or classical gene therapy approaches as novel pharmacological strategies.

Rapid birth-and-death evolution of imprinted snoRNAs in the Prader-Willi syndrome locus: implications for neural development in Euarchontoglires

Imprinted small nucleolar RNAs (snoRNAs) are only found in eutherian genomes and closely related to brain functions. A complex human neurological disease, Prader-Willi syndrome (PWS), is primarily attributed to the deletion of imprinted snoRNAs in chromosome 15q11-q13. Here we investigated the snoRNA repertoires in the PWS locus of 12 mammalian genomes and their evolution processes. A total of 613 imprinted snoRNAs were identified in the PWS homologous loci and the gene number was highly variable across lineages, with a peak in Euarchontoglires. Lineage-specific gene gain and loss events account for most extant genes of the HBII-52 (SNORD115) and the HBII-85 (SNORD116) gene family, and remarkable high gene-birth rates were observed in the primates and the rodents. Meanwhile, rapid sequence substitution occurred only in imprinted snoRNA genes, rather than their flanking sequences or the protein-coding genes located in the same imprinted locus. Strong selective constraints on the functional elements of these imprinted snoRNAs further suggest that they are subjected to birth-and-death evolution. Our data suggest that the regulatory role of HBII-52 on 5-HT_2CR pre-mRNA might originate in the Euarchontoglires through adaptive process. We propose that the rapid evolution of PWS-related imprinted snoRNAs has contributed to the neural development of Euarchontoglires.

Stereological quantification of olfactory receptor neurons in mice

The total number of olfactory receptor neurons (ORNs) in the mouse main olfactory epithelium (MOE) was estimated using stereological sampling. Noses and skulls of male and female 8-week-old C57BL/6J mice were de-calcified, embedded in paraffin, cut into 10-μm-thick sections serially at 100-μm intervals, and processed for immunohistochemistry for the olfactory marker protein (OMP), a specific marker for ORNs. The number of OMP (+) receptor neurons was measured using an optical fractionator with the Stereo-Investigator software. The mean values of the total number of OMP (+) receptor neurons in the unilateral MOE were 5,140,000±380,000 in males and 5,210,000±380,000 in females, with no significant differences between the sexes. We concluded that the total number of ORNs in the unilateral MOE is approximately 5×10(6) in mice.

The membrane proteome of sensory cilia to the depth of olfactory receptors

In the nasal cavity, the nonmotile cilium of olfactory sensory neurons (OSNs) constitutes the chemosensory interface between the ambient environment and the brain. The unique sensory organelle facilitates odor detection for which it includes all necessary components of initial and downstream olfactory signal transduction. In addition to its function in olfaction, a more universal role in modulating different signaling pathways is implicated, for example, in neurogenesis, apoptosis, and neural regeneration. To further extend our knowledge about this multifunctional signaling organelle, it is of high importance to establish a most detailed proteome map of the ciliary membrane compartment down to the level of transmembrane receptors. We detached cilia from mouse olfactory epithelia via Ca(2+)/K(+) shock followed by the enrichment of ciliary membrane proteins at alkaline pH, and we identified a total of 4,403 proteins by gel-based and gel-free methods in conjunction with high resolution LC/MS. This study is the first to report the detection of 62 native olfactory receptor proteins and to provide evidence for their heterogeneous expression at the protein level. Quantitative data evaluation revealed four ciliary membrane-associated candidate proteins (the annexins ANXA1, ANXA2, ANXA5, and S100A5) with a suggested function in the regulation of olfactory signal transduction, and their presence in ciliary structures was confirmed by immunohistochemistry. Moreover, we corroborated the ciliary localization of the potassium-dependent Na(+)/Ca(2+) exchanger (NCKX) 4 and the plasma membrane Ca(2+)-ATPase 1 (PMCA1) involved in olfactory signal termination, and we detected for the first time NCKX2 in olfactory cilia. Through comparison with transcriptome data specific for mature, ciliated OSNs, we finally delineated the membrane ciliome of OSNs. The membrane proteome of olfactory cilia established here is the most complete today, thus allowing us to pave new avenues for the study of diverse molecular functions and signaling pathways in and out of olfactory cilia and thus to advance our understanding of the biology of sensory organelles in general.

Strong links between genomic and anatomical diversity in both mammalian olfactory chemosensory systems

Mammalian olfaction comprises two chemosensory systems: the odorant-detecting main olfactory system (MOS) and the pheromone-detecting vomeronasal system (VNS). Mammals are diverse in their anatomical and genomic emphases on olfactory chemosensation, including the loss or reduction of these systems in some orders. Despite qualitative evidence linking the genomic evolution of the olfactory systems to specific functions and phenotypes, little work has quantitatively tested whether the genomic aspects of the mammalian olfactory chemosensory systems are correlated to anatomical diversity. We show that the genomic and anatomical variation in these systems is tightly linked in both the VNS and the MOS, though the signature of selection is different in each system. Specifically, the MOS appears to vary based on absolute organ and gene family size while the VNS appears to vary according to the relative proportion of functional genes and relative anatomical size and complexity. Furthermore, there is little evidence that these two systems are evolving in a linked fashion. The relationships between genomic and anatomical diversity strongly support a role for natural selection in shaping both the anatomical and genomic evolution of the olfactory chemosensory systems in mammals.

Genetic signatures for enhanced olfaction in the African mole-rats

The Olfactory Receptor (OR) superfamily, the largest in the vertebrate genome, is responsible for vertebrate olfaction and is traditionally subdivided into 17 OR families. Recent studies characterising whole-OR subgenomes revealed a 'birth and death' model of evolution for a range of species, however little is known about fine-scale evolutionary dynamics within single-OR families. This study reports the first assessment of fine-scale OR evolution and variation in African mole-rats (Bathyergidae), a family of subterranean rodents endemic to sub-Saharan Africa. Because of the selective pressures of life underground, enhanced olfaction is proposed to be fundamental to the evolutionary success of the Bathyergidae, resulting in a highly diversified OR gene-repertoire. Using a PCR-sequencing approach, we analysed variation in the OR7 family across 14 extant bathyergid species, which revealed enhanced levels of functional polymorphisms concentrated across the receptors' ligand-binding region. We propose that mole-rats are able to recognise a broad range of odorants and that this diversity is reflected throughout their OR7 gene repertoire. Using both classic tests and tree-based methods to test for signals of selection, we investigate evolutionary forces across the mole-rat OR7 gene tree. Four well-supported clades emerged in the OR phylogeny, with varying signals of selection; from neutrality to positive and purifying selection. Bathyergid life-history traits and environmental niche-specialisation are explored as possible drivers of adaptive OR evolution, emerging as non-exclusive contributors to the positive selection observed at OR7 genes. Our results reveal unexpected complexity of evolutionary mechanisms acting within a single OR family, providing insightful perspectives into OR evolutionary dynamics.

Deleted copy number variation of Hanwoo and Holstein using next generation sequencing at the population level

BACKGROUND

Copy number variation (CNV), a source of genetic diversity in mammals, has been shown to underlie biological functions related to production traits. Notwithstanding, there have been few studies conducted on CNVs using next generation sequencing at the population level.

RESULTS

Illumina NGS data was obtained for ten Holsteins, a dairy cattle, and 22 Hanwoo, a beef cattle. The sequence data for each of the 32 animals varied from 13.58-fold to almost 20-fold coverage. We detected a total of 6,811 deleted CNVs across the analyzed individuals (average length = 2732.2 bp) corresponding to 0.74% of the cattle genome (18.6 Mbp of variable sequence). By examining the overlap between CNV deletion regions and genes, we selected 30 genes with the highest deletion scores. These genes were found to be related to the nervous system, more specifically with nervous transmission, neuron motion, and neurogenesis. We regarded these genes as having been effected by the domestication process. Further analysis of the CNV genotyping information revealed 94 putative selected CNVs and 954 breed-specific CNVs.

CONCLUSIONS

This study provides useful information for assessing the impact of CNVs on cattle traits using NGS at the population level.

Comparative genomics provide insights into evolution of trichoderma nutrition style

Saprotrophy on plant biomass is a recently developed nutrition strategy for Trichoderma. However, the physiology and evolution of this new nutrition strategy is still elusive. We report the deep sequencing and analysis of the genome of Trichoderma longibrachiatum, an efficient cellulase producer. The 31.7-Mb genome, smallest among the sequenced Trichoderma species, encodes fewer nutrition-related genes than saprotrophic T. reesei (Tr), including glycoside hydrolases and nonribosomal peptide synthetase-polyketide synthase. Homology and phylogenetic analyses suggest that a large number of nutrition-related genes, including GH18 chitinases, β-1,3/1,6-glucanases, cellulolytic enzymes, and hemicellulolytic enzymes, were lost in the common ancestor of T. longibrachiatum (Tl) and Tr. dN/dS (ω) calculation indicates that all the nutrition-related genes analyzed are under purifying selection. Cellulolytic enzymes, the key enzymes for saprotrophy on plant biomass, are under stronger purifying selection pressure in Tl and Tr than in mycoparasitic species, suggesting that development of the nutrition strategy of saprotrophy on plant biomass has increased the selection pressure. In addition, aspartic proteases, serine proteases, and metalloproteases are subject to stronger purifying selection pressure in Tl and Tr, suggesting that these enzymes may also play important roles in the nutrition. This study provides insights into the physiology and evolution of the nutrition strategy of Trichoderma.

A cluster of olfactory receptor genes linked to frugivory in bats

Diversity of the mammalian olfactory receptor (OR) repertoire has been globally reshaped by niche specialization. However, little is known about the variability of the OR repertoire at a shallower evolutionary timeframe. The vast bat radiation exhibits an extraordinary variety of trophic and sensory specializations. Unlike other mammals, bats possess a unique and diverse OR gene repertoire. We elucidated whether the evolution of the OR gene repertoire can be linked to ecological niche specializations, such as sensory modalities and diet. The OR gene repertoires of 27 bat species spanning the chiropteran radiation were amplified and sequenced. For each species, intact and nonfunctional genes were assessed, and the OR gene abundances in each gene family were analyzed and compared. We identified a unique OR pattern linked to the frugivorous diet of New World fruit-eating bats and a similar convergent pattern in the Old World fruit-eating bats. Our results show a strong association between niche specialization and OR repertoire diversity even at a shallow evolutionary timeframe.

New genes as drivers of phenotypic evolution

During the course of evolution, genomes acquire novel genetic elements as sources of functional and phenotypic diversity, including new genes that originated in recent evolution. In the past few years, substantial progress has been made in understanding the evolution and phenotypic effects of new genes. In particular, an emerging picture is that new genes, despite being present in the genomes of only a subset of species, can rapidly evolve indispensable roles in fundamental biological processes, including development, reproduction, brain function and behaviour. The molecular underpinnings of how new genes can develop these roles are starting to be characterized. These recent discoveries yield fresh insights into our broad understanding of biological diversity at refined resolution.

Age distribution patterns of human gene families: divergent for Gene Ontology categories and concordant between different subcellular localizations

The age distribution of gene duplication events within the human genome exhibits two waves of duplications along with an ancient component. However, because of functional constraint differences, genes in different functional categories might show dissimilar retention patterns after duplication. It is known that genes in some functional categories are highly duplicated in the early stage of vertebrate evolution. However, the correlations of the age distribution pattern of gene duplication between the different functional categories are still unknown. To investigate this issue, we developed a robust pipeline to date the gene duplication events in the human genome. We successfully estimated about three-quarters of the duplication events within the human genome, along with the age distribution pattern in each Gene Ontology (GO) slim category. We found that some GO slim categories show different distribution patterns when compared to the whole genome. Further hierarchical clustering of the GO slim functional categories enabled grouping into two main clusters. We found that human genes located in the duplicated copy number variant regions, whose duplicate genes have not been fixed in the human population, were mainly enriched in the groups with a high proportion of recently duplicated genes. Moreover, we used a phylogenetic tree-based method to date the age of duplications in three signaling-related gene superfamilies: transcription factors, protein kinases and G-protein coupled receptors. These superfamilies were expressed in different subcellular localizations. They showed a similar age distribution as the signaling-related GO slim categories. We also compared the differences between the age distributions of gene duplications in multiple subcellular localizations. We found that the distribution patterns of the major subcellular localizations were similar to that of the whole genome. This study revealed the whole picture of the evolution patterns of gene functional categories in the human genome.