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.

Electronic supplementary material

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

Comparative genomics search for losses of long-established genes on the human lineage

Taking advantage of the complete genome sequences of several mammals, we developed a novel method to detect losses of well-established genes in the human genome through syntenic mapping of gene structures between the human, mouse, and dog genomes. Unlike most previous genomic methods for pseudogene identification, this analysis is able to differentiate losses of well-established genes from pseudogenes formed shortly after segmental duplication or generated via retrotransposition. Therefore, it enables us to find genes that were inactivated long after their birth, which were likely to have evolved nonredundant biological functions before being inactivated. The method was used to look for gene losses along the human lineage during the approximately 75 million years (My) since the common ancestor of primates and rodents (the euarchontoglire crown group). We identified 26 losses of well-established genes in the human genome that were all lost at least 50 My after their birth. Many of them were previously characterized pseudogenes in the human genome, such as GULO and UOX. Our methodology is highly effective at identifying losses of single-copy genes of ancient origin, allowing us to find a few well-known pseudogenes in the human genome missed by previous high-throughput genome-wide studies. In addition to confirming previously known gene losses, we identified 16 previously uncharacterized human pseudogenes that are definitive losses of long-established genes. Among them is ACYL3, an ancient enzyme present in archaea, bacteria, and eukaryotes, but lost approximately 6 to 8 Mya in the ancestor of humans and chimps. Although losses of well-established genes do not equate to adaptive gene losses, they are a useful proxy to use when searching for such genetic changes. This is especially true for adaptive losses that occurred more than 250,000 years ago, since any genetic evidence of the selective sweep indicative of such an event has been erased.     

Gene loss through pseudogenization contributes to the ecological diversification of a generalist Roseobacter lineage

Ecologically relevant genes generally show patchy distributions among related bacterial genomes. This is commonly attributed to lateral gene transfer, whereas the opposite mechanism—gene loss—has rarely been explored. Pseudogenization is a major mechanism underlying gene loss, and pseudogenes are best characterized by comparing closely related genomes because of their short life spans. To explore the role of pseudogenization in microbial ecological diversification, we apply rigorous methods to characterize pseudogenes in the 279 newly sequenced Ruegeria isolates of the globally abundant Roseobacter group collected from two typical coastal habitats in Hong Kong, the coral Platygyra acuta and the macroalga Sargassum hemiphyllum. Pseudogenes contribute to ~16% of the accessory genomes of these strains. Ancestral state reconstruction reveals that many pseudogenization events are correlated with ancestral niche shifts. Specifically, genes related to resource scavenging and energy acquisition were often pseudogenized when roseobacters inhabiting carbon-limited and energy-poor coral skeleton switched to other resource-richer niches. For roseobacters inhabiting the macroalgal niches, genes for nitrogen regulation and carbohydrate utilization were important but became dispensable upon shift to coral skeleton where nitrate is abundant but carbohydrates are less available. Whereas low-energy-demanding secondary transporters are more favorable in coral skeleton, ATP-driven primary transporters are preferentially kept in the energy-replete macroalgal niches. Moreover, a large proportion of these families mediate organismal interactions, suggesting their rapid losses by pseudogenization as a potential response to host and niche shift. These findings illustrate an important role of pseudogenization in shaping genome content and driving ecological diversification of marine roseobacters.Subject terms: Ecology, Evolution, Microbiology     

Similar numbers but different repertoires of olfactory receptor genes in humans and chimpanzees

Animals recognize their external world through the detection of tens of thousands of chemical odorants. Olfactory receptor (OR) genes encode proteins for detecting odorant molecules and form the largest multigene family in mammals. It is known that humans have fewer OR genes and a higher fraction of OR pseudogenes than mice or dogs. To investigate whether these features are human specific or common to all higher primates, we identified nearly complete sets of OR genes from the chimpanzee and macaque genomes and compared them with the human OR genes. In contrast to previous studies, here we show that the number of OR genes ( approximately 810) and the fraction of pseudogenes (51%) in chimpanzees are very similar to those in humans, though macaques have considerably fewer OR genes. The pseudogenization rates and the numbers of genes affected by positive selection are also similar between humans and chimpanzees. Moreover, the most recent common ancestor between humans and chimpanzees had a larger number of functional OR genes (>500) and a lower fraction of pseudogenes (41%) than its descendents, suggesting that the OR gene repertoires are in a phase of deterioration in both lineages. Interestingly, despite the close evolutionary relationship between the 2 species, approximately 25% of their functional gene repertoires are species specific due to massive gene losses. These findings suggest that the tempo of evolution of OR genes is similar between humans and chimpanzees, but the OR gene repertoires are quite different between them. This difference might be responsible for the species-specific ability of odor perception.     

Divergence,demography and gene loss along the human lineage

Genomic DNA sequences are an irreplaceable source for reconstructing the vanished past of living organisms. Based on updated sequence data, this paper summarizes our studies on species divergence time, ancient population size and functional loss of genes in the primate lineage leading to modern humans (Homo sapiens sapiens). The inter- and intraspecific comparisons of DNA sequences suggest that the human lineage experienced a rather severe bottleneck in the Middle Pleistocene, throughout which period the subdivided African population played a predominant role in shaping the genetic architecture of modern humans. Also, published and newly identified human-specific pseudogenes (HSPs) are enumerated in order to infer their significance for human evolution. Of the 121 candidate genes obtained, authentic HSPs turn out to comprise only 25 olfactory receptor genes, four T cell receptor genes and nine other genes. The fixation of HSPs has been too rare over the past 6–7 Myr to account for species differences between humans and chimpanzees.     

Evolution of Human-Specific Alleles Protecting Cognitive Function of Grandmothers

The myelomonocytic receptor CD33 (Siglec-3) inhibits innate immune reactivity by extracellular V-set domain recognition of sialic acid (Sia)-containing “self-associated molecular patterns” (SAMPs). We earlier showed that V-set domain-deficient CD33-variant allele, protective against late-onset Alzheimer’s Disease (LOAD), is derived and specific to the hominin lineage. We now report multiple hominin-specific CD33 V-set domain mutations. Due to hominin-specific, fixed loss-of-function mutation in the CMAH gene, humans lack N-glycolylneuraminic acid (Neu5Gc), the preferred Sia-ligand of ancestral CD33. Mutational analysis and molecular dynamics (MD)-simulations indicate that fixed change in amino acid 21 of hominin V-set domain and conformational changes related to His45 corrected for Neu5Gc-loss by switching to N-acetylneuraminic acid (Neu5Ac)-recognition. We show that human-specific pathogens Neisseria gonorrhoeae and Group B Streptococcus selectively bind human CD33 (huCD33) as part of immune-evasive molecular mimicry of host SAMPs and that this binding is significantly impacted by amino acid 21 modification. In addition to LOAD-protective CD33 alleles, humans harbor derived, population-universal, cognition-protective variants at several other loci. Interestingly, 11 of 13 SNPs in these human genes (including CD33) are not shared by genomes of archaic hominins: Neanderthals and Denisovans. We present a plausible evolutionary scenario to compile, correlate, and comprehend existing knowledge about huCD33-evolution and suggest that grandmothering emerged in humans.     

Loss of gene function and evolution of human phenotypes

Humans have acquired many distinct evolutionary traits after the human-chimpanzee divergence. These phenotypes have resulted from genetic changes that occurred in the human genome and were retained by natural selection. Comparative primate genome analyses reveal that loss-of-function mutations are common in the human genome. Some of these gene inactivation events were revealed to be associated with the emergence of advantageous phenotypes and were therefore positively selected and fixed in modern humans (the “less-ismore” hypothesis). Representative cases of human gene inactivation and their functional implications are presented in this review. Functional studies of additional inactive genes will provide insight into the molecular mechanisms underlying acquisition of various human-specific traits. [BMB Reports 2015; 48(7): 373-379]     

The use of racial, ethnic, and ancestral categories in human genetics research

The global dispersal of anatomically modern humans over the past 100,000 years has produced patterns of phenotypic variation that have exerted—and continue to exert—powerful influences on the lives of individuals and the experiences of groups. The recency of our common ancestry and continued gene flow among populations have resulted in less genetic differentiation among geographically distributed human populations than is observed in many other mammalian species. Nevertheless, differences in appearance have contributed to the development of ideas about “race” and “ethnicity” that often include the belief that significant inherited differences distinguish humans. The use of racial, ethnic, and ancestral categories in genetics research can imply that group differences arise directly through differing allele frequencies, with little influence from socially mediated mechanisms. At the same time, careful investigations of the biological, environmental, social, and psychological attributes associated with these categories will be an essential component of cross-disciplinary research into the origins, prevention, and treatment of common diseases, including those diseases that differ in prevalence among groups.     

Lineage-specific loss of function of bitter taste receptor genes in humans and nonhuman primates

Since the process of becoming dead genes or pseudogenes (pseudogenization) is irreversible and can occur rather rapidly under certain environmental circumstances, it is one plausible determinant for characterizing species specificity. To test this evolutionary hypothesis, we analyzed the tempo and mode of duplication and pseudogenization of bitter taste receptor (T2R) genes in humans as well as in 12 nonhuman primates. The results show that primates have accumulated more pseudogenes than mice after their separation from the common ancestor and that lineage-specific pseudogenization becomes more conspicuous in humans than in nonhuman primates. Although positive selection has operated on some amino acids in extracellular domains, functional constraints against T2R genes are more relaxed in primates than in mice and this trend has culminated in the rapid deterioration of the bitter-tasting capability in humans. Since T2R molecules play an important role in avoiding generally bitter toxic and harmful substances, substantial modification of the T2R gene repertoire is likely to reflect different responses to changes in the environment and to result from species-specific food preference during primate evolution.     

MicroRNA expression and regulation in human, chimpanzee, and macaque brains

Among other factors, changes in gene expression on the human evolutionary lineage have been suggested to play an important role in the establishment of human-specific phenotypes. However, the molecular mechanisms underlying these expression changes are largely unknown. Here, we have explored the role of microRNA (miRNA) in the regulation of gene expression divergence among adult humans, chimpanzees, and rhesus macaques, in two brain regions: prefrontal cortex and cerebellum. Using a combination of high-throughput sequencing, miRNA microarrays, and Q-PCR, we have shown that up to 11% of the 325 expressed miRNA diverged significantly between humans and chimpanzees and up to 31% between humans and macaques. Measuring mRNA and protein expression in human and chimpanzee brains, we found a significant inverse relationship between the miRNA and the target genes expression divergence, explaining 2%-4% of mRNA and 4%-6% of protein expression differences. Notably, miRNA showing human-specific expression localize in neurons and target genes that are involved in neural functions. Enrichment in neural functions, as well as miRNA-driven regulation on the human evolutionary lineage, was further confirmed by experimental validation of predicted miRNA targets in two neuroblastoma cell lines. Finally, we identified a signature of positive selection in the upstream region of one of the five miRNA with human-specific expression, miR-34c-5p. This suggests that miR-34c-5p expression change took place after the split of the human and the Neanderthal lineages and had adaptive significance. Taken together these results indicate that changes in miRNA expression might have contributed to evolution of human cognitive functions.     

Natural selection in the evolution of SARS-CoV-2 in bats created a generalist virus and highly capable human pathogen

Virus host shifts are generally associated with novel adaptations to exploit the cells of the new host species optimally. Surprisingly, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has apparently required little to no significant adaptation to humans since the start of the Coronavirus Disease 2019 (COVID-19) pandemic and to October 2020. Here we assess the types of natural selection taking place in Sarbecoviruses in horseshoe bats versus the early SARS-CoV-2 evolution in humans. While there is moderate evidence of diversifying positive selection in SARS-CoV-2 in humans, it is limited to the early phase of the pandemic, and purifying selection is much weaker in SARS-CoV-2 than in related bat Sarbecoviruses. In contrast, our analysis detects evidence for significant positive episodic diversifying selection acting at the base of the bat virus lineage SARS-CoV-2 emerged from, accompanied by an adaptive depletion in CpG composition presumed to be linked to the action of antiviral mechanisms in these ancestral bat hosts. The closest bat virus to SARS-CoV-2, RmYN02 (sharing an ancestor about 1976), is a recombinant with a structure that includes differential CpG content in Spike; clear evidence of coinfection and evolution in bats without involvement of other species. While an undiscovered “facilitating” intermediate species cannot be discounted, collectively, our results support the progenitor of SARS-CoV-2 being capable of efficient human–human transmission as a consequence of its adaptive evolutionary history in bats, not humans, which created a relatively generalist virus.

A study of the natural origins of SARS-CoV-2 reveals very little adaptive evolution occurring since it emerged in humans, but strong evolutionary signals in the bat virus lineage from which SARS-CoV-2 arose. Evolution in bats involved lineage-specific depletion of CpG nucleotides (linked to host anti-viral molecules), and clear evidence of recombination across these lineages, supporting bat host species’ influence on the ancestral viruses.     

Human-specific amino acid changes found in 103 protein-coding genes

We humans have many characteristics that are different from those of the great apes. These human-specific characters must have arisen through mutations accumulated in the genome of our direct ancestor after the divergence of the last common ancestor with chimpanzee. Gene trees of human and great apes are necessary for extracting these human-specific genetic changes. We conducted a systematic analysis of 103 protein-coding genes for human, chimpanzee, gorilla, and orangutan. Nucleotide sequences for 18 genes were newly determined for this study, and those for the remaining genes were retrieved from the DDBJ/EMBL/GenBank database. The total number of amino acid changes in the human lineage was 147 for 26,199 codons (0.56%). The total number of amino acid changes in the human genome was, thus, estimated to be about 80,000. We applied the acceleration index test and Fisher's synonymous/nonsynonymous exact test for each gene tree to detect any human-specific enhancement of amino acid changes compared with ape branches. Six and two genes were shown to have significantly higher nonsynonymous changes at the human lineage from the acceleration index and exact tests, respectively. We also compared the distribution of the differences of the nonsynonymous substitutions on the human lineage and those on the great ape lineage. Two genes were more conserved in the ape lineage, whereas one gene was more conserved in the human lineage. These results suggest that a small proportion of protein-coding genes started to evolve differently in the human lineage after it diverged from the ape lineage.     

Harvesting of novel polyhydroxyalkanaote (PHA) synthase encoding genes from a soil metagenome library using phenotypic screening

Salmonella enterica serovar Typhi and Typhimurium are closely related serovars. However, S. Typhi, a human-specific pathogen, has 5% of genes as pseudogenes, far more than S. Typhimurium, which only has 1%. One of these pseudogenes corresponds to sopD2, which in S. Typhimurium encodes an effector protein involved in Salmonella-containing vacuole biogenesis in human epithelial cell lines, which is needed for full virulence of the pathogen. We investigated whether S. Typhi trans-complemented with the functional sopD2 gene from S. Typhimurium (sopD2(STM) ) would reduce the invasion of human epithelial cell lines. Our results showed that the presence of sopD2(STM) in S. Typhi significantly modified the bacterial ability to alter cellular permeability and decrease the CFUs recovered after cell invasion of human epithelial cell line. These results add to mounting evidence that pseudogenes contribute to S. Typhi adaptation to humans.     

High-resolution copy-number variation map reflects human olfactory receptor diversity and evolution

Olfactory receptors (ORs), which are involved in odorant recognition, form the largest mammalian protein superfamily. The genomic content of OR genes is considerably reduced in humans, as reflected by the relatively small repertoire size and the high fraction ( approximately 55%) of human pseudogenes. Since several recent low-resolution surveys suggested that OR genomic loci are frequently affected by copy-number variants (CNVs), we hypothesized that CNVs may play an important role in the evolution of the human olfactory repertoire. We used high-resolution oligonucleotide tiling microarrays to detect CNVs across 851 OR gene and pseudogene loci. Examining genomic DNA from 25 individuals with ancestry from three populations, we identified 93 OR gene loci and 151 pseudogene loci affected by CNVs, generating a mosaic of OR dosages across persons. Our data suggest that approximately 50% of the CNVs involve more than one OR, with the largest CNV spanning 11 loci. In contrast to earlier reports, we observe that CNVs are more frequent among OR pseudogenes than among intact genes, presumably due to both selective constraints and CNV formation biases. Furthermore, our results show an enrichment of CNVs among ORs with a close human paralog or lacking a one-to-one ortholog in chimpanzee. Interestingly, among the latter we observed an enrichment in CNV losses over gains, a finding potentially related to the known diminution of the human OR repertoire. Quantitative PCR experiments performed for 122 sampled ORs agreed well with the microarray results and uncovered 23 additional CNVs. Importantly, these experiments allowed us to uncover nine common deletion alleles that affect 15 OR genes and five pseudogenes. Comparison to the chimpanzee reference genome revealed that all of the deletion alleles are human derived, therefore indicating a profound effect of human-specific deletions on the individual OR gene content. Furthermore, these deletion alleles may be used in future genetic association studies of olfactory inter-individual differences.     

GLADX: an automated approach to analyze the lineage-specific loss and pseudogenization of genes

A well-established ancestral gene can usually be found, in one or multiple copies, in different descendant species. Sometimes during the course of evolution, all the representatives of a well-established ancestral gene disappear in specific lineages; such gene losses may occur in the genome by deletion of a DNA fragment or by pseudogenization. The loss of an entire gene family in a given lineage may reflect an important phenomenon, and could be due either to adaptation, or to a relaxation of selection that leads to neutral evolution. Therefore, the lineage-specific gene loss analyses are important to improve the understanding of the evolutionary history of genes and genomes. In order to perform this kind of study from the increasing number of complete genome sequences available, we developed a unique new software module called GLADX in the DAGOBAH framework, based on a comparative genomic approach. The software is able to automatically detect, for all the species of a phylum, the presence/absence of a representative of a well-established ancestral gene, and by systematic steps of re-annotation, confirm losses, detect and analyze pseudogenes and find novel genes. The approach is based on the use of highly reliable gene phylogenies, of protein predictions and on the analysis of genomic mutations. All the evidence associated to evolutionary approach provides accurate information for building an overall view of the evolution of a given gene in a selected phylum. The reliability of GLADX has been successfully tested on a benchmark analysis of 14 reported cases. It is the first tool that is able to fully automatically study the lineage-specific losses and pseudogenizations. GLADX is available at http://ioda.univ-provence.fr/IodaSite/gladx/.     

Accelerated evolution of the pituitary adenylate cyclase-activating polypeptide precursor gene during human origin

Pituitary adenylate cyclase-activating polypeptide (PACAP) is a neuropeptide abundantly expressed in the central nervous system and involved in regulating neurogenesis and neuronal signal transduction. The amino acid sequence of PACAP is extremely conserved across vertebrate species, indicating a strong functional constraint during the course of evolution. However, through comparative sequence analysis, we demonstrated that the PACAP precursor gene underwent an accelerated evolution in the human lineage since the divergence from chimpanzees, and the amino acid substitution rate in humans is at least seven times faster than that in other mammal species resulting from strong Darwinian positive selection. Eleven human-specific amino acid changes were identified in the PACAP precursors, which are conserved from murine to African apes. Protein structural analysis suggested that a putative novel neuropeptide might have originated during human evolution and functioned in the human brain. Our data suggested that the PACAP precursor gene underwent adaptive changes during human origin and may have contributed to the formation of human cognition.     

Irish Cepaea nemoralis Land Snails Have a Cryptic Franco-Iberian Origin That Is Most Easily Explained by the Movements of Mesolithic Humans

The origins of flora and fauna that are only found in Ireland and Iberia, but which are absent from intervening countries, is one of the enduring questions of biogeography. As Southern French, Iberian and Irish populations of the land snail Cepaea nemoralis sometimes have a similar shell character, we used mitochondrial phylogenies to begin to understand if there is a shared “Lusitanian” history. Although much of Europe contains snails with A and D lineages, by far the majority of Irish individuals have a lineage, C, that in mainland Europe was only found in a restricted region of the Eastern Pyrenees. A past extinction of lineage C in the rest of Europe cannot be ruled out, but as there is a more than 8000 year continuous record of Cepaea fossils in Ireland, the species has long been a food source in the Pyrenees, and the Garonne river that flanks the Pyrenees is an ancient human route to the Atlantic, then we suggest that the unusual distribution of the C lineage is most easily explained by the movements of Mesolithic humans. If other Irish species have a similarly cryptic Lusitanian element, then this raises the possibility of a more widespread and significant pattern.