A uniquely human consequence of domain-specific functional adaptation in a sialic acid-binding receptor

Most mammalian cell surfaces display two major sialic acids (Sias), N-acetylneuraminic acid (Neu5Ac) and N-glycolylneuraminic acid (Neu5Gc). Humans lack Neu5Gc due to a mutation in CMP-Neu5Ac hydroxylase, which occurred after evolutionary divergence from great apes. We describe an apparent consequence of human Neu5Gc loss: domain-specific functional adaptation of Siglec-9, a member of the family of sialic acid-binding receptors of innate immune cells designated the CD33-related Siglecs (CD33rSiglecs). Binding studies on recombinant human Siglec-9 show recognition of both Neu5Ac and Neu5Gc. In striking contrast, chimpanzee and gorilla Siglec-9 strongly prefer binding Neu5Gc. Simultaneous probing of multiple endogenous CD33rSiglecs on circulating blood cells of human, chimp, or gorilla suggests that the binding differences observed for Siglec-9 are representative of multiple CD33rSiglecs. We conclude that Neu5Ac-binding ability of at least some human CD33rSiglecs is a derived state selected for following loss of Neu5Gc in the hominid lineage. These data also indicate that endogenous Sias (rather than surface Sias of bacterial pathogens) are the functional ligands of CD33rSiglecs and suggest that the endogenous Sia landscape is the major factor directing evolution of CD33rSiglec binding specificity. Exon-1-encoded Sia-recognizing domains of human and ape Siglec-9 share only approximately 93-95% amino acid identity. In contrast, the immediately adjacent intron and exon 2 have the approximately 98-100% identity typically observed among these species. Together, our findings suggest ongoing adaptive evolution specific to the Sia-binding domain, possibly of an episodic nature. Such domain-specific divergences should also be considered in upcoming comparisons of human and chimpanzee genomes.     

Terrestrial vertebrates have two keratin gene clusters; striking differences in teleost fish

Keratins I and II form the largest subgroups of mammalian intermediate filament (IF) proteins and account as obligatory heteropolymers for the keratin filaments of epithelia. All human type I genes except for the K18 gene are clustered on chromosome 17q21, while all type II genes form a cluster on chromosome 12q13, that ends with the type I gene K18. Highly related keratin gene clusters are found in rat and mouse. Since fish seem to lack a keratin II cluster we screened the recently established draft genomes of a bird (chicken) and an amphibian (Xenopus). The results show that keratin I and II gene clusters are a feature of all terrestrial vertebrates. Because hair with its multiple hair keratins and inner root sheath keratins is a mammalian acquisition, the keratin gene clusters of chicken and Xenopus tropicalis have only about half the number of genes found in mammals. Within the type I clusters all genes have the same orientation. In type II clusters there is a rare gene of opposite orientation. Finally we show that the genes for keratins 8 and 18, which are the first expression pair in embryology, are not only adjacent in mammals, but also in Xenopus and three different fish. Thus neighboring K8 and K18 genes seem a feature shared by all vertebrates. In contrast to the two well defined keratin gene clusters of terrestrial vertebrates, three teleost fish show an excess of type I over type II genes, the lack of a keratin type II gene cluster and a striking dispersal of type I genes, that are probably the result of the teleost-specific whole genome duplication followed by a massive gene loss. This raises the question whether keratin gene clusters extend beyond the ancestral bony vertebrate to cartilage fish and lamprey. We also analyzed the complement of non-keratin IF genes of the chicken. Surprisingly, an additional nuclear lamin gene, previously overlooked by cDNA cloning, is documented on chromosome 10. The two splice variants closely resemble the lamin LIII a + b of amphibia and fish. This lamin gene is lost on the mammalian lineage.     

Conserved synteny of mammalian imprinted genes in chicken, frog, and fish genomes

Conservation of synteny of mammalian imprinted genes between chicken and human suggested that highly conserved gene clusters were selected long before these genes were recruited for genomic imprinting in mammals. Here we have applied in silico mapping of orthologous genes in pipid frog, zebrafish, spotted green and Japanese pufferfish to show considerable conservation of synteny in lower vertebrates. More than 400 million years ago in a common ancestor of teleost fish and tetrapods, 'preimprinted' chromosome regions homologous to human 6q25, 7q21, 7q32, 11p15, and 15q11-->q12 already contained most present-day mammalian imprinted genes. Interestingly, some imprinted gene orthologues which are isolated from imprinted clusters in mouse and human could be linked to preimprinted regions in lower vertebrates, indicating that separation occurred during mammalian evolution. On the contrary, newly arisen genes by segmental duplication in the mammalian lineage, i.e. SNRPN and FRAT3, were transposed or translocated to imprinted clusters and recruited for parent-specific activity. By analysis of currently available sequences of non-mammalian vertebrates, the imprinted gene clusters homologous to human chromosomes 14q32 and 19q12 are only poorly conserved in chicken, frog, and fish and, therefore, may not have evolved from ancestral preimprinted gene arrays. Evidently, evolution of imprinted gene clusters is an ongoing and dynamic process in mammals. In general, imprinted gene orthologues do not show a higher degree of synteny conservation in vertebrates than non-imprinted genes interspersed with or adjacent to an imprinted cluster.     

Identification and Comparative Analysis of the Protocadherin Cluster in a Reptile,the Green Anole Lizard

Background

The vertebrate protocadherins are a subfamily of cell adhesion molecules that are predominantly expressed in the nervous system and are believed to play an important role in establishing the complex neural network during animal development. Genes encoding these molecules are organized into a cluster in the genome. Comparative analysis of the protocadherin subcluster organization and gene arrangements in different vertebrates has provided interesting insights into the history of vertebrate genome evolution. Among tetrapods, protocadherin clusters have been fully characterized only in mammals. In this study, we report the identification and comparative analysis of the protocadherin cluster in a reptile, the green anole lizard (Anolis carolinensis).

Methodology/Principal Findings

We show that the anole protocadherin cluster spans over a megabase and encodes a total of 71 genes. The number of genes in the anole protocadherin cluster is significantly higher than that in the coelacanth (49 genes) and mammalian (54–59 genes) clusters. The anole protocadherin genes are organized into four subclusters: the δ, α, β and γ. This subcluster organization is identical to that of the coelacanth protocadherin cluster, but differs from the mammalian clusters which lack the δ subcluster. The gene number expansion in the anole protocadherin cluster is largely due to the extensive gene duplication in the γb subgroup. Similar to coelacanth and elephant shark protocadherin genes, the anole protocadherin genes have experienced a low frequency of gene conversion.

Conclusions/Significance

Our results suggest that similar to the protocadherin clusters in other vertebrates, the evolution of anole protocadherin cluster is driven mainly by lineage-specific gene duplications and degeneration. Our analysis also shows that loss of the protocadherin δ subcluster in the mammalian lineage occurred after the divergence of mammals and reptiles. We present a model for the evolutionary history of the protocadherin cluster in tetrapods.     

Group B streptococcal capsular sialic acids interact with siglecs (immunoglobulin-like lectins) on human leukocytes

Group B Streptococcus (GBS) is classified into nine serotypes that vary in capsular polysaccharide (CPS) architecture but share in common the presence of a terminal sialic acid (Sia) residue. This position and linkage of GBS Sia closely resembles that of cell surface glycans found abundantly on human cells. CD33-related Siglecs (CD33rSiglecs) are a family of Sia-binding lectins expressed on host leukocytes that engage host Sia-capped glycans and send signals that dampen inflammatory gene activation. We hypothesized that GBS evolved to display CPS Sia as a form of molecular mimicry limiting the activation of an effective innate immune response. In this study, we applied a panel of immunologic and cell-based assays to demonstrate that GBS of several serotypes interacts in a Sia- and serotype-specific manner with certain human CD33rSiglecs, including hSiglec-9 and hSiglec-5 expressed on neutrophils and monocytes. Modification of GBS CPS Sia by O acetylation has recently been recognized, and we further show that the degree of O acetylation can markedly affect the interaction between GBS and hSiglec-5, -7, and -9. Thus, production of Sia-capped bacterial polysaccharide capsules that mimic human cell surface glycans in order to engage CD33rSiglecs may be an example of a previously unrecognized bacterial mechanism of leukocyte manipulation.     

Analysis of the genomic structure of the porcine CD1 gene cluster

CD1 is an MHC class I-like protein that presents lipid antigens to T cell receptors. We determined 470,187 bp of the genomic sequence encompassing the region encoding porcine CD1 genes. We identified 16 genes in this region and newly identified CD1A2, CD1B, CD1C, CD1D, and CD1E. Porcine CD1 genes were located in clusters between KIRREL and olfactory receptor (OR) genes, as observed in humans, although they were divided into two regions by a region encoding OR genes. Comparison of the genomic sequences of CD1 gene loci in pigs with other mammals showed that separation of the CD1 gene cluster by ORs was observed only in pigs. CD1A duplication in the porcine genome was estimated to have occurred after the divergence of the human and porcine. This analysis of the genomic sequence of the porcine CD1 family will contribute to our understanding of the evolution of mammalian CD1 genes.     

Chicken CD69 and CD94/NKG2-like genes in a chromosomal region syntenic to mammalian natural killer gene complex

In mammals, natural killer (NK) cell C-type lectin receptors were encoded in a gene cluster called natural killer gene complex (NKC). The NKC is not reported in chicken yet. Instead, NK receptor genes were found in the major histocompatibility complex. In this study, two novel chicken C-type lectin-like receptor genes were identified in a region on chromosome 1 that is syntenic to mammalian NKC region. The chromosomal locations were validated with fluorescent in situ hybridization. Based on 3D structure modeling, sequence homology, chromosomal location, and phlylogenetic analysis, one receptor is the orthologue of mammalian cluster of differentiation 69 (CD69), and the other is highly homologous to CD94 and NKG2. Like CD94/NKG2 gene found in teleostean fishes, chicken CD94/NKG2 has the features of both human CD94 and NKG2A. Unlike mammalian NKC, these two chicken C-type lectin receptors are not closely linked but separated by 42 million base pairs according to the chicken draft genome sequence. The arrangement of several other genes that are located outside the mammalian NKC is conserved among chicken, human, and mouse. The chicken NK C-type lectin-like receptors in the NKC syntenic region indicate that this chromosomal region existed before the divergence between mammals and aves. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequences have been submitted to the GenBank nucleotide sequence database under the accession number chicken CD69 (DQ156495), CD94/NKG2 (DQ156496), and CD94/NKG2 variant (DQ241793).     

Cloning and characterization of human Siglec-11. A recently evolved signaling molecule that can interact with SHP-1 and SHP-2 and is expressed by tissue macrophages, including brain microglia

Siglecs are sialic acid-recognizing animal lectins of the immunoglobulin superfamily. We have cloned and characterized a novel human molecule, Siglec-11, that belongs to the subgroup of CD33/Siglec-3-related Siglecs. As with others in this subgroup, the cytosolic domain of Siglec-11 is phosphorylated at tyrosine residue(s) upon pervanadate treatment of cells and then recruits the protein-tyrosine phosphatases SHP-1 and SHP-2. However, Siglec-11 has several novel features relative to the other CD33/Siglec-3-related Siglecs. First, it binds specifically to alpha2-8-linked sialic acids. Second, unlike other CD33/Siglec-3-related Siglecs, Siglec-11 was not found on peripheral blood leukocytes. Instead, we observed its expression on macrophages in various tissues, such as liver Kupffer cells. Third, it was also expressed on brain microglia, thus becoming the second Siglec to be found in the nervous system. Fourth, whereas the Siglec-11 gene is on human chromosome 19, it lies outside the previously described CD33/Siglec-3-related Siglec cluster on this chromosome. Fifth, analyses of genome data bases indicate that Siglec-11 has no mouse ortholog and that it is likely to be the last canonical human Siglec to be reported. Finally, although Siglec-11 shows marked sequence similarity to human Siglec-10 in its extracellular domain, the cytosolic tail appears only distantly related. Analysis of genomic regions surrounding the Siglec-11 gene suggests that it is actually a chimeric molecule that arose from relatively recent gene duplication and recombination events, involving the extracellular domain of a closely related ancestral Siglec gene (which subsequently became a pseudogene) and a transmembrane and cytosolic tail derived from another ancestral Siglec.     

Identification of duplicated fourth alpha2-adrenergic receptor subtype by cloning and mapping of five receptor genes in zebrafish

The alpha(2)-adrenergic receptors (alpha(2)-ARs) belong to the large family of rhodopsinlike G-protein-coupled receptors that share a common structure of seven transmembrane (TM) alpha-helices. The aims of this study were (1) to determine the number of alpha(2)-AR genes in a teleost fish, the zebrafish (Danio rerio), (2) to study the gene duplication events that generated the alpha(2)-AR subtypes, and (3) to study changes in receptor structure that have occurred since the divergence of the mammalian and fish lineages. Here, we report the cloning and chromosomal mapping of fish orthologs for all three mammalian alpha(2)-ARs. In addition, we identified a fourth alpha(2)-AR subtype with two duplicates in zebrafish. Chromosomal mapping showed that the zebrafish alpha(2)-AR genes are located within conserved chromosomal segments, consistent with the origin of the four alpha(2)-AR subtypes by two rounds of chromosome or block duplication before the divergence of the ray fin fish and tetrapod lineages. Thus, the fourth subtype has apparently been present in the common ancestor of vertebrates but has been deleted or is yet to be identified in mammals. The overall percentage identity between the fish and mammalian orthologs is 53% to 67%, and in the TM regions 80% to 87%. These values are clearly lower than what is observed between mammalian orthologs. Still, all of the residues thought to be important for alpha(2)-adrenergic ligand binding are conserved across species and subtypes, and even the most divergent regions of the fish receptors show clear "molecular fingerprints" typical for orthologs of a given subtype.     

Cloning and characterization of a novel mouse Siglec, mSiglec-F: differential evolution of the mouse and human (CD33) Siglec-3-related gene clusters

A novel mouse Siglec (mSiglec-F) belonging to the subfamily of Siglec-3-related Siglecs has been cloned and characterized. Unlike most human Siglec-3 (hSiglec-3)-related Siglecs with promiscuous linkage specificity, mSiglec-F shows a strong preference for alpha2-3-linked sialic acids. It is predominantly expressed in immature cells of the myelomonocytic lineage and in a subset of CD11b (Mac-1)-positive cells in some tissues. As with previously cloned Siglec-3-related mSiglecs, the lack of strong sequence similarity to a singular hSiglec made identification of the human ortholog difficult. We therefore conducted a comprehensive comparison of Siglecs between the human and mouse genomes. The mouse genome contains eight Siglec genes, whereas the human genome contains 11 Siglec genes and a Siglec-like gene. Although a one-to-one orthologous correspondence between human and mouse Siglecs 1, 2, and 4 is confirmed, the Siglec-3-related Siglecs showed marked differences between human and mouse. We found only four Siglec genes and two pseudogenes in the mouse chromosome 7 region syntenic to the Siglec-3-related gene cluster on human chromosome 19, which, in contrast, contains seven Siglec genes, a Siglec-like gene, and thirteen pseudogenes. Although analysis of gene maps and exon structures allows tentative assignments of mouse-human Siglec ortholog pairs, the possibility of unequal genetic recombination makes the assignments inconclusive. We therefore support a temporary lettered nomenclature for additional mouse Siglecs. Current information suggests that mSiglec-F is likely a hSiglec-5 ortholog. The previously reported mSiglec-3/CD33 and mSiglec-E/MIS are likely orthologs of hSiglec-3 and hSiglec-9, respectively. The other Siglec-3-like gene in the cluster (mSiglec-G) is probably a hSiglec-10 ortholog. Another mouse gene (mSiglec-H), without an apparent human ortholog, lies outside of the cluster. Thus, although some duplications of Siglec-3-related genes predated separation of the primate and rodent lineages (about 80-100 million years ago), this gene cluster underwent extensive duplications in the primate lineage thereafter.     

Duplicated KOX zinc finger gene clusters flank the centromere of human chromosome 10: evidence for a pericentric inversion during primate evolution.

Two related zinc finger (ZNF) gene clusters from the pericentromeric region of human chromosome 10, defined by cDNAs of the KOX series, have been cloned in yeast artificial chromosomes (YACs). The two clusters evolved by duplication of an ancestral gene cluster before the divergence of the human and great ape lineages. Included in cluster A are the ZNF gene sequences ZNF11A, ZNF33A, and ZNF37A, while cluster B comprises the related sequences ZNF11B, ZNF33B and ZNF37B. Genes from both clusters are expressed: cDNAs KOX2, KOX31 and KOX21 derive from ZNF11B, ZNF33A and ZNF37A, respectively. Further YACs have been isolated which link ZNF11A and ZNF33A to another gene, ZNF25, defined by cDNA clone KOX19. Therefore ZNF25 also forms part of cluster A, but has no counterpart in cluster B. Surprisingly, the KOX ZNF gene clusters are located on opposite sides of the centromere: cluster A maps to 10p11.2, while cluster B is in 10q11.2. This suggests the occurrence during primate evolution of a previously undescribed pericentric inversion subsequent to the cluster duplication. The evolution of this subset of KOX ZNF genes has therefore involved three types of genetic event: local gene duplication, gene cluster duplication, and chromosome rearrangement.     

Evolution of immunoglobulin kappa chain variable region genes in vertebrates

The major source of immunoglobulin diversity is variation in DNA sequence among multiple copies of variable region (V) genes of the heavy- and light-chain multigene families. In order to clarify the evolutionary pattern of the multigene family of immunoglobulin light kappa chain V region (V kappa) genes, phylogenetic analyses of V kappa genes from humans and other vertebrate species were conducted. The results obtained indicate that the V kappa genes so far sequenced can be grouped into three major monophyletic clusters, the cartilaginous fish, bony fish and amphibian, and mammalian clusters, and that the cartilaginous fish cluster first separated from the rest of the V kappa genes and then the remaining two clusters diverged. The mammalian V kappa genes can further be divided into 10 V kappa groups, 7 of which are present in the human genome. Human and mouse V kappa genes from different V kappa groups are intermingled rather than clustered on the chromosome, and there are a large number of pseudogenes scattered on the chromosome. This indicates that the chromosomal locations of V kappa genes have been shuffled many times by gene duplication, deletion, and transposition in the evolutionary process and that many genes have become nonfunctional during this process. This mode of evolution is consistent with the model of birth-and-death evolution rather than with the model of concerted evolution. An analysis of duplicate V kappa functional genes and pseudogenes in the human genome has indicated that pseudogenes evolve faster than functional genes but that the rate of nonsynonymous nucleotide substitution in the complementarity-determining regions of V kappa genes has been enhanced by positive Darwinian selection.      

Organization of<Emphasis Type="Italic"> Iroquois</Emphasis> genes in fish

In mammals, a total of six iroquois (Irx) genes exist, which are organized into two clusters. Here we report on the organization of all iroquois genes present in fish, using zebrafish (Danio rerio) and pufferfish (Fugu rubripes and Tetraodon nigroviridis) as examples. A total of 10 Irx genes were found in pufferfish, and 11 in zebrafish; all but one of these genes are organized into clusters (four clusters plus one isolated gene locus). The extra fish clusters result from chromosome duplication in the fish lineage, after its divergence from tetrapod vertebrates. Two of the four fish clusters are highly conserved to the ones in mammals, with regard to similarity of genes and cluster architecture. Irx genes within the other two clusters have diverged in sequence and cluster organization, suggesting functional divergence. These results will allow us to use the zebrafish system for functional and comparative studies of iroquois genes in vertebrate development.Electronic Supplementary Material Supplementary material is available in the online version of this article at Edited by D. Tautz     

Evolution of glucagon genes

Statistical analyses of DNA sequences of the preproglucagon genes from bovine, human, hamster, and anglerfish suggest that a gene duplication creating two anglerfish genes (AF I and II) occurred about 160 Myr ago, long after the separation of fish and mammals. The analyses further suggest that the internal duplication producing the glucagon and glucagon-like peptide II (GLP-II) regions occurred about 1.2 billion years ago, which would indicate that the GLP-II region was present in the ancestral anglerfish sequence but was silenced or deleted before the gene duplication separating AF I and II. The glucagon-like peptide I (GLP-I) was derived from a duplication of the ancestral glucagon region about 800 Myr ago. The rate of synonymous substitution in these genes is approximately 4.3 x 10(-9) substitutions per year per synonymous site. The rate of nonsynonymous substitution in the signal peptide region is about 1.1 x 10(-9) substitutions per year per nonsynonymous site, a high rate comparable to that in the C-peptide region of preproinsulin. The rate of nonsynonymous substitution in the glicentin-related pancreatic polypeptide (GRPP) region is 0.63 x 10(-9) for the comparisons between mammalian species and 1.8 x 10(-9) for the comparisons between fish and mammals; the moderate rate in mammals suggests a physiological role for GRPP. The glucagon region is extremely conservative; no nonsynonymous substitution is observed in the mammalian genes, and a nonsynonymous rate of 0.18 x 10(-9) was obtained from the comparisons between fish and mammals. In the GLP-I region, the rate of nonsynonymous substitution was estimated to be 0.08 x 10(-9) for the comparisons between mammalian species and 0.30 x 10(- 9) for the comparisons between fish and mammals. In the GLP-II region, the rate was estimated to be 0.25 x 10(-9) for the comparisons between mammalian species. Thus, GLP-I and II are also very conservative, which suggests an important physiological role for these peptides.      

Differential evolution of voltage-gated sodium channels in tetrapods and teleost fishes

The voltage-gated sodium channel (SCN) alpha subunits are large proteins with central roles in the generation of action potentials. They consist of approximately 2,000 amino acids encoded by 24-27 exons. Previous evolutionary studies have been unable to reconcile the proposed gene duplication schemes with the species distribution and molecular phylogeny of the genes. We have carefully annotated the complete SCN gene sequences, correcting numerous database errors, for a broad range of vertebrate species and analyzed their phylogenetic relationships. We have also compared the chromosomal positions of the SCN genes relative to adjacent gene families. Our studies show that the ancestor of the vertebrates probably had a single sodium channel gene with two characteristic AT-AC introns, the second of which is unique to vertebrate SCN genes. This ancestral gene, located close to a HOX gene cluster, was quadrupled along with HOX in the two rounds of basal vertebrate tetraploidizations to generate the ancestors of the four channels SCN1A, SCN4A, SCN5A, and SCN8A. The third tetraploidization in the teleost fish ancestor doubled this set of genes and all eight are still present in at least three of four investigated teleost fish genomes. In tetrapods, the gene family expanded by local duplications before the radiation of amniotes, generating the cluster SCN5A, SCN10A, and SCN11A on one chromosome and the cluster SCN1A, SCN2A, SCN3A, and SCN9A on a different chromosome. In eutherian mammals, a tenth gene, SCN7A, arose in a local duplication in the SCN1A gene cluster. The SCN7A gene has undergone rapid evolution and has lost the ability to cause action potentials-instead, it functions as a sodium sensor. The three genes in the SCN5A cluster were translocated from the HOX-bearing chromosome in a mammalian ancestor along with several adjacent genes. This evolutionary scenario is supported by the adjacent TGF-β receptor superfamily (comprised of five distinct families) and the cysteine-serine-rich nuclear protein gene family as well as the HOX clusters. The independent expansions of the SCN repertoires in tetrapods and teleosts suggest that the functional diversification may differ between the two lineages.     

Differential evolution of the 13 Atlantic salmon Hox clusters

Hox cluster organization represents a valuable marker to study the effects of recent genome duplication in salmonid fish (25-100 Mya). Using polymerase chain reaction amplification of cDNAs, BAC library screening, and genome walking, we reconstructed 13 Hox clusters in the Atlantic salmon containing 118 Hox genes including 8 pseudogenes. Hox paralogs resulting from the genome duplication preceding the radiation of ray-finned fish have been much better preserved in salmon than in other model teleosts. The last genome duplication in the salmon lineage has been followed by the loss of 1 of the 4 HoxA clusters. Four rounds of genome duplication after the vertebrate ancestor salmon Hox clusters display the main organizational features of vertebrate Hox clusters, with Hox genes exclusively that are densely packed in the same orientation. Recently, duplicated Hox clusters have engaged a process of divergence, with several cases of pseudogenization or asymmetrical evolution of Hox gene duplicates, and a marked erosion of identity in noncoding sequences. Strikingly, the level of divergence attained strongly depends on the Hox cluster pairs rather than on the Hox genes within each cluster. It is particularly high between both HoxBb clusters and both HoxDa clusters, whereas both HoxBa clusters remained virtually identical. Positive selection on the Hox protein-coding sequences could not be detected.     

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.     

基因倍增和脊椎动物起源

有机体基因复制导致基因复杂性增加及其和脊椎动物起源的关系已经成为进化生物学研究的热点。20世纪70年代由Ohno提出后经Holland等修正的原始脊索动物经两轮基因组复制产生脊椎动物的假设目前已被广泛接受。脊椎动物起源和进化过程中发生过两轮基因组复制的主要证据有三点：（1）据估计脊椎动物基因组内编码基因数目大约相当于果蝇、海鞘等无脊椎动物的4倍;原口动物如果蝇和后口动物如头索动物文昌鱼的基因组大都只有单拷贝的基因,而脊椎动物的基因组则通常有4个同属于一个家族的基因。（2）无脊椎动物如节肢动物、海胆和头索动物文昌鱼都只有一个Hox基因簇,而脊椎动物除鱼类外,有7个具有Hox基因簇,其余都具有4个Hox基因簇。（3）基因作图证明,不但在鱼类和哺乳动物染色体广大片段上基因顺序相似,而且有证据显示哺乳动物基因组不同染色体之间存在相似性。据认为第一次基因倍增发生在脊椎动物与头索动物分开之后,第二次基因倍增发生在有颌类脊椎动物和无颌类脊椎动物分开以后。但是,基因是逐个发生倍增,还是通过基因组内某些DNA片段抑或整个基因组的加倍而实现的,目前还颇有争议。