The acquisition of novel N-glycosylation sites in conserved proteins during human evolution
BMC Bioinformatics volume 16, Article number: 29 (2015)
N-linked protein glycosylation plays an important role in various biological processes, including protein folding and trafficking, and cell adhesion and signaling. The acquisition of a novel N-glycosylation site may have significant effect on protein structure and function, and therefore, on the phenotype.
We analyzed the human glycoproteome data set (2,534 N-glycosylation sites in 1,027 proteins) and identified 112 novel N-glycosylation sites in 91 proteins that arose in the human lineage since the last common ancestor of Euarchonta (primates and treeshrews). Three of them, Asn-196 in adipocyte plasma membrane-associated protein (APMAP), Asn-91 in cluster of differentiation 166 (CD166/ALCAM), and Asn-76 in thyroglobulin, are human-specific. Molecular evolutionary analysis suggested that these sites were under positive selection during human evolution. Notably, the Asn-76 of thyroglobulin might be involved in the increased production of thyroid hormones in humans, especially thyroxine (T4), because the removal of the glycan moiety from this site was reported to result in a significant decrease in T4 production.
We propose that the novel N-glycosylation sites described in this study may be useful candidates for functional analyses to identify innovative genetic modifications for beneficial phenotypes acquired in the human lineage.
N-linked glycosylation of the Asn residue in the consensus motif Asn-X-Ser/Thr, where X is any amino acid except Pro, is one of the most well-studied protein posttranslational modifications (PTMs) . N-glycosylation, which mainly occurs in secreted or cell membrane proteins, plays important roles in protein folding, quality control, and trafficking , as well as cell adhesion and signalling [3,4]. The emergence of a new N-glycosylation site may alter protein function either positively or negatively. For example, missense mutations in factor VIII created novel N-glycosylation sites that cause severe hemophilia A . Similarly, a missense mutation in the interferon γ receptor 2 induces novel N-glycosylation, which results in a Mendelian susceptibility to mycobacterial disease [6,7]. The abolishment of N-glycosylation sites often causes disrupted protein folding, trafficking, or activity; thus, proper N-glycosylation is crucial for normal protein function [8,9]. A proteome-wide analysis of nonsynonymous single-nucleotide variations in the N-glycosylation motifs of human genes showed that more than 1,000 human proteins had either lost or gained N-glycosylation sites due to missense substitutions, some of which may be implicated in diseases .
The gain of new N-glycosylation sites during evolution may affect the structure and molecular function of proteins; when these novel modifications confer beneficial traits, they will be fixed during evolution. Previously, we identified a large variety of genetic changes that could have been involved in the acquisition of human traits, including gene inactivation [11,12], exon evolution [13,14], and gains of phosphorylation or ubiquitylation [15,16]. Therefore, it would be of great interest to collect information on novel N-glycosylation sites that arose during human evolution, as the sites might have been involved in the development of some human phenotypes.
In order to study associations between the acquisition of an N-glycosylation site and its phenotypic outcome, a large amount of N-glycosylation site data and mammalian orthologous protein sequence data are required. Recent developments and advances in various high-throughput proteomics techniques for N-glycoproteome identification using immunoaffinity subtraction, hydrazide chemistry, and mass spectrometry have made it possible to access massive amounts of N-glycosylation site data from human proteomes [17-19]. These data are available at the UniProt database (http://www.uniprot.org), which is a universal protein sequence database, as well as some specialized PTM databases such as PHOSIDA (http://www.phosida.com/) .
Since human genome sequences were completed [21,22], a large amount of nucleotide and protein sequence data have become available not only from humans but also from many other organisms. Comparative sequence data, including alignments of mammalian orthologous protein sequences, are available at the University of California Santa Cruz (UCSC) Genome Browser Database (http://genome.ucsc.edu) .
In this study, a bioinformatics method was devised to identify novel N-glycosylated Asn residues that are located in the consensus motif Asn-X-Ser/Thr and arose during human evolution after the Euarchonta lineage diverged from the Glires lineage. Both a comprehensive literature survey and extensive data mining were conducted to examine the possible functional implications of novel N-glycosylation sites, especially in cases of human-specific gains.
Identification of novel N-glycosylation sites acquired during human evolution and determination of the timing of acquisition
We developed a bioinformatics procedure to identify the acquisition by proteins of the N-glycosylation motif Asn-X-Ser/Thr, where the Asn residue was experimentally verified to be N-glycosylated, during human evolution (Figure 1). A novel N-glycosylation site can arise by the emergence of not only an Asn residue but also a Ser or Thr residue to form the consensus motif . The overall procedure devised in this study is similar to that used to identify novel ubiquitylation sites in a previous study . Initially, there were 2,534 experimentally verified human N-glycosylation sites from 1,027 proteins in the UniProt database, and 57,289 orthologous protein sequence alignments from 62 mammalian species, including species from Euarchonta, Glires, Laurasiatheria, Afrotheria, Xenarthra, Marsupialia, and Monotremata, extracted from the UCSC “multiz100way” data  (see Additional file 1 for the list of mammalian species). These data were analyzed to collect N-glycosylation sites in human proteins that newly appeared during the evolution from the common ancestor, Euarchonta (primates and treeshrews); as the result, 112 novel N-glycosylation sites from 91 proteins were identified. A summary of the results are presented in Additional file 2, and detailed alignments are provided in Additional file 3. Of the 91 proteins, one protein (CFH) had acquired four N-glycosylation sites (Nos. 28 to 31 in Additional file 2 and Additional file 3; two proteins (PTPRC and PTPRJ) had acquired three sites each (Nos. 85 to 87, and 88 to 90, respectively); 14 proteins had acquired two sites each; and the remaining 74 proteins had acquired one site each. Figure 2 shows the number of the N-glycosylation sites that are shared by each of the Euarchonta clades along the human lineage: humans, three; ancestor of humans and chimpanzees, two; African great apes, 10; great apes, four; apes, 12; catarrhines, 16; simians, 45; primates, 15; and euarchonts, 5.
Most of the novel N-glycosylation sites were generated by the emergence of an Asn residue in an existing X-X-Ser/Thr motif. However, in some cases, the emergence of a Ser or a Thr residue in an Asn-X-X sequence created a novel N-glycosylation site: for example, a change from Asn-Glu-Ile to Asn-Glu-Thr generated a novel N-glycosylation site at Asn-911 in complement factor H (CFH) in apes (No. 30 in Additional file 3).
Of the 112 novel N-glycosylation sites, three sites in three proteins were human specific (Table 1 and Figure 3); therefore, these Asn residues subject to the N-glycosylation evolved and were fixed in human proteins after the divergence of humans and chimpanzees. The residues are Asn-196 in adipocyte plasma membrane-associated protein (APMAP), Asn-91 in cluster of differentiation 166 (CD166), and Asn-76 in thyroglobulin.
Human-specific N-glycosylation site Asn-196 in APMAP
The human APMAP (also known as C20orf3) has two N-glycosylation sites: Asn-160 and Asn-196, the latter of which is human specific: almost all the other mammals examined have a Lys residue at this position (Figure 3A). Full-length protein and coding sequences of APMAP orthologs were determined from human, chimpanzee, gorilla, orangutan, gibbon, and rhesus macaque genomes (see Methods section for details). Multiple alignment of these proteins showed that human APMAP protein has two human-specific amino acid changes, Val-100 and Asn-196, where all the other five primates have Ile and Lys, respectively (Additional file 4). The Asn-196 is a human-specific N-glycosylation site.
To test if the APMAP protein has been under positive selection during human evolution, the ratios of nonsynonymous to synonymous rates (dN/dS, ω) across different branches and sites of the selected primate phylogeny were estimated [26-28]. First, we used “branch models”, M0 (one ω ratio for all branches), free ratio (one ω ratio for each branch), and two ratio (ω1 for the human branch and ω0 for other branches) models (Table 2 and Additional file 4). The likelihood ratio test (LRT) comparing M0 (one ratio) and free ratio model was not significant. However, the LRT comparing M0 and two ratio model was highly significant (P = 0.006943), suggesting the human APMAP has evolved at different rate compared to other primates. The estimated dN and dS rates for the human branch using two ratio model were 0.0024 and 0.0000 (see Additional file 4 for details), respectively, indicating possible accelerated nonsynonymous substitution during human evolution. Next, we used “branch-site models”, model A (ω ratio is left to vary) and null model A (ω ratio is fixed to 1), to infer positively selected sites in human APMAP. The two aforementioned human-specific amino acid positions, Val-100 and Asn-196, were detected to be under positive selection with overall probability of 0.828 and 0.953, respectively, using the Bayes empirical Bayes (BEB) test . However, the LRT comparing model A and null model A was not significant. Although it is statistically insignificant, the acquisition of Asn-196 and its subsequent N-glycosylation might have a significant effect on the structure and function of APMAP in humans.
APMAP is an adipocyte plasma membrane-associated protein, which is induced during adipocyte differentiation . It is ubiquitously expressed in human embryonic and adult tissues, with the highest levels in liver, placenta, and kidney . APMAP may exhibit calcium-dependent hydrolase activity and is regulated by the peroxisome proliferator activated receptor γ protein that is a master regulator of adipocyte differentiation . Expression of APMAP was reported to be strongly correlated with hepatic-specific metastasis in patients with metastatic colorectal cancer . A recent study demonstrated that APMAP is a negative regulator of amyloid-beta (Aβ) production through its interaction with amyloid precursor protein and γ-secretase . Although APMAP seems to be involved in various biological processes in humans, the molecular function directly associated with the human-specific N-glycosylation site in APMAP is yet to be determined.
Human-specific N-glycosylation site Asn-91 in CD166
The human CD166, which is also known as activated leukocyte cell adhesion molecule (ALCAM), has 10 N-glycosylation sites. The residue Asn-91 is found in humans but not in other euarchonts; thus, this residue evolved after the divergence of humans and chimpanzees (Figure 3B). Most other mammals have a Ser residue at this position. However, some Glires and Laurasiatheria species, especially whales and ruminants, independently acquired a consensus sequence for N-glycosylation at this position.
Multiple alignment of full-length CD166 proteins from human, chimpanzee, gorilla, orangutan, gibbon, and rhesus macaque genomes revealed that the Asn-91 is only residue that differs between humans and chimpanzees (Additional file 5). Although the Asn-91 is human-specific, all the LRTs (M0 versus free ratio model, M0 versus two ratio model, and model A versus null model A) were insignificant (Table 2 and Additional file 5), implying that there has been no statistically noticeable positive selection on the human CD166. Interestingly, the human-specific N-glycosylation site Asn-91 was still inferred to be under positive selection with overall probability of 0.613, suggesting the acquisition of this site and its N-glycosylation might have an effect on the function of human CD166 proteins.
CD166 binds to the T-cell differentiation antigen CD6 and may play a role in the binding of T and B cells to activated leukocytes, as well as in interactions between cells of the nervous system . CD166 is composed of five extracellular immunoglobulin (Ig)-like domains: two Ig-like V-type domains and three Ig-like C2-type domains. The human-specific N-glycosylation site Asn-91 is located in the first Ig-like V-type domain, which mediates the CD166–CD6 interaction [35-37]. Most functional studies on CD166 have focused on its cancer-related functions such as invasion, migration, and adhesion [37,38]. However, recent studies show that CD166 is also involved in axon growth in neuronal cells such as retinal ganglion cells and dorsal root ganglion cells [39,40]. Therefore, it is possible that the gain of the Asn-91 N-glycosylation site in CD166 might be involved in the evolution of novel phenotypes in nervous system development, as well as in immune response and cell adhesion processes, which must be validated experimentally.
Human-specific N-glycosylation site Asn-76 in thyroglobulin
Thyroglobulin is the precursor of the thyroid hormones T4 and triiodothyronine (T3), both of which regulate metabolism in humans [41-43]. The human thyroglobulin has 17 N-glycosylation sites; Asn-76, which becomes Asn-57 in the mature form of thyroglobulin, is a human-specific N-glycosylation site: most other mammals have Asp or Glu at this position (Figure 3C). Interestingly, multiple alignment of thyroglobulin orthologs from humans, chimpanzees, gorillas, orangutans, gibbons, and rhesus macaques revealed additional 20 amino acid positions with human-specific substitution (Additional file 6).
To test if the thyroglobulin has been under positive selection during human evolution, ω ratios across different branches and sites of the selected primate phylogeny were estimated (Table 2 and Additional file 6). The LRT comparing M0 and free ratio model was not significant. However, the LRT comparing M0 and two ratio model was highly significant (P = 0.006511), implying the human thyroglobulin has evolved at different rate compared to other primates. The estimated ω ratio for the human branch (ω1) using two ratio model was 0.78473, while ω0 for other branches was 0.33454 (see Additional file 6 for details), suggesting a slightly accelerated nonsynonymous substitution during human evolution.
Inference of positively selected sites in human thyroglobulin using model A showed that the 21 aforementioned amino acid positions might have been under positive selection although the LRT comparing model A and null model A was not significant. In spite of statistical insignificance, the acquisition of novel N-glycosylation site Asn-76, together with other 20 human-specific amino acid substitutions, might have a significant effect on the thyroid hormone metabolism in humans.
The thyroglobulin protein precursor itself has no biological function but serves as a chemical platform for thyroid hormone production. When the two N-glycosylation sites in its N-terminal region, including the human-specific Asn-76, were deglycosylated by peptide-N4-(N-acetyl-β-glucosaminyl)-asparagine amidase, T4 production decreased by half compared to that seen with the normal protein . Therefore, proper N-glycosylation modifications, including those at the human-specific site, in thyroglobulin are crucial for normal T4 production and the control of metabolism. It has been suggested that humans and chimpanzees differ with respect to their thyroid hormone metabolism . Interestingly, compared to chimpanzees, humans have a higher T4 plasma concentration, which may be implicated in the origins of human intelligence .
Novel N-glycosylation sites shared by other animals
Of the 112 novel N-glycosylation sites in human proteins, 109 sites were shared by other animals. For example, 12 N-glycosylation sites were shared by all apes, indicating that this site appeared in the common ancestor of apes. The Asn-480 of pappalysin-1 (also known as pregnancy-associated plasma protein-A), which is N-glycosylated in humans , is shared by all apes examined (Figure 4A); in contrast, all the other mammals, even Marsupialia and Monotremata species, have an Asp residue at this position. Pappalysin-1 has metalloproteinase activity and specifically cleaves insulin-like growth factor-binding proteins [48,49]; it is present at high concentration in maternal blood during pregnancy and is essential for normal fetal development . The serum pappalysin-1 concentration frequently increases in patients with severe sepsis and appears to be associated with sepsis-related myocardial dysfunction . However, there is no comprehensive study of whether the gain of Asn-480 in apes is associated with these phenotypes.
There are 12 N-glycosylation sites that might have arisen in the common ancestor of apes. One is Asn-113 in the thyrotropin receptor, or thyroid-stimulating hormone receptor (TSHR) (Figure 4B). The THSR responds to thyroid-stimulating hormone (also known as thyrotropin) and stimulates the production of T4 and T3 in the thyroid gland . Human THSR has six N-glycosylation sites: Asn-113 is specific to apes, and the other five are conserved in other mammals. However, a mutated TSHR in which N-glycosylation at Asn-113 had been disrupted had the same expression level and function as the wildtype TSHR; thus, Asn-113 N-glycosylation may not be important for TSHR function . Therefore, gain of N-glycosylation at Asn-113 may be neutral or have a function yet to be determined.
The UDP-glucuronosyltransferase 1–9 (UD19) Asn-344 is one of 45 N-glycosylation sites that are shared by simians (apes and monkeys; Figure 4C). Nonsimian mammals have a Lys residue at this position. UD19, which is also known as UDP glucuronosyltransferase 1 family, polypeptide A9, is involved in the conjugation and elimination of toxic xenobiotic and endogenous compounds [54,55]. Unglycosylation of UD19 resulted in the inhibition of proper protein folding and the impairment of glucuronidation activity; thus, N-glycosylation plays a role in folding the human UD19 protein . UD19 is one of nine functional isoforms produced by the alternative utilization of the first nine exons in the UGT1A gene locus . Because Asn-344 is located in the common exon 5, not only UD19 but also eight other isoforms of UDP glucuronosyltransferase 1 enzyme have this novel N-glycosylation site .
Previously, it has been suggested that the gain of novel protein PTM sites such as ubiquitylation sites may be associated with the acquisition of novel phenotypes during human evolution by modulating the activity or network of proteins . It is also highly probable that gains of novel N-glycosylation sites may result in functional modification of proteins and phenotypic changes in an organism. In this study, 1,027 human glycoproteins containing experimentally verified N-glycosylation sites and their orthologous mammalian proteins were systematically compared. As a result, 112 novel N-glycosylation sites were identified in 91 proteins that newly appeared during human evolution after the Euarchonta lineage diverged from the Glires lineage. It must be noted that most of these novel N-glycosylation sites were obtained by high-throughput mass spectrometry. The presence of these modifications must be further scrutinized by conventional molecular biology techniques.
Not all the novel N-glycosylation sites described in this study may have resulted in functional innovation. Some of them might have appeared as a result of random genetic drift and be functionally neutral. Nevertheless, some of them could have conferred selective advantage during human evolution and be fixed in the human genome. One such example identified in this study is the novel N-glycosylation site in UD19, which is involved in the elimination of potentially toxic xenobiotics and endogenous compounds. UD19 acquired the novel N-glycosylation site Asn-344 during the evolution of the common ancestor of apes and monkeys (see Figure 4C). When the N-glycosylation at Asn-344 is abolished, folding is inhibited in UD19, and its glucuronidase activity is reduced . Therefore, glycosylation at Asn-344 is required for proper folding and activity of UD19. It is possible that ancestral simian primates required better defense mechanisms against toxic compounds introduced into their systems by environmental or dietary shifts. The acquisition of a new N-glycosylation site in UD19 might have conferred improved xenobiotics metabolism to apes and monkeys, although there is no direct evidence for this hypothesis.
The three human-specific N-glycosylation sites are particularly interesting (see Table 1 and Figure 3). The residue Asn-196 in APMAP is the first of the three human-specific N-glycosylation sites, which was inferred to be positively selected with an extremely high probability in humans (see Table 2 and Additional file 4). The human APMAP has been reported to be involved in a variety of biological processes including adipocyte differentiation, hepatic-specific metastasis in cancer, and inhibition of Aβ production [32-34]. The fact that APMAP is implicated in adipocyte differentiation is particularly interesting because humans and great apes exhibit large differences in adipose tissue and fatty acid storage, and these differences may be associated with the development of subcutaneous fat and even in brain development [58,59]. Therefore, the molecular functional study of human-specific sequence changes in proteins such as APMAP, which are associated with adipose tissue and lipid metabolism, may reveal the molecular mechanisms for the evolution of these traits.
The human CD166 protein has two Ig-like V-type domains and three Ig-like C2-type domains (Figure 5A) and functions as a cell adhesion molecule. The human-specific N-glycosylation site Asn-91, which was inferred to be positively selected (see Table 2 and Additional file 5), is located within the first Ig-like V-type domain, which is responsible for protein–protein interactions [35-37]. The addition of a bulky glycan moiety to this domain may change its structural profile and thus affect cell–cell adhesion activity or ligand specificity. The most interesting function of CD166 is its involvement in axon growth in neuronal cells [39,40]. Of the 583 amino acid residues in CD166, only the residue at position 91 differs between humans and chimpanzees; therefore, the emergence of Asn-91 and its N-glycosylation might be associated with evolution of human-specific phenotypes, probably in the nervous system, which must be determined experimentally.
The human thyroglobulin, which serves as a precursor molecule for the thyroid hormones T4 and T3, has a human-specific N-glycosylation site Asn-76. The human thyroglobulin contains 11 thyroglobulin type-1 domains, which are involved in the control of proteolytic degradation . The human-specific Asn-76 is located within the first thyroglobulin type-1 domain (Figure 5B). The Asn-76 was inferred to be positively selected during human evolution, along with the other 20 positions (see Table 2 and Additional file 6). Removal of the glycan group from the Asn-76 reduced thyroid hormone production, especially T4 production ; thus, the gain of Asn-76 and its N-glycosylation, together with the other 20 putatively adaptive amino acid changes, may be implicated in the increased T4 concentration present in humans as compared to chimpanzees. It is possible that the additional glycan moiety may confer increased resistance to the proteolytic degradation of thyroglobulin proteins and thus lead to increased thyroid hormone production. It is very interesting that the T4 concentration in humans is higher than that in chimpanzees , as elevated T4 production may have caused the modification of human physiology in response to selection pressures in a specific environment: specifically, it has been proposed that an altered thyroid hormone metabolism might have been beneficial for early humans in the savannah environment, as they practiced persistence hunting and thus had large energy requirements .
Losses of N-glycosylation sites during human evolution are also very interesting. Some human diseases have been reported to be caused by the loss of N-glycosylation sites . To find cases where ancestrally conserved N-glycosylation sites were lost during human evolution, a large amount of N-glycosylation data collected from animals distantly related to humans is required. The N-glycoproteome data obtained from mouse tissues and plasma using high-throughput mass spectrometry would be an ideal dataset for this analysis . With a simple modification, the procedure described in this study could be used to analyze these data for the identification of N-glycosylation sites that were lost during human evolution and their possible phenotypic implications.
We have devised and applied a bioinformatics method to identify the acquisition of N-glycosylation sites during human evolution. We propose that the acquisition of novel N-glycosylation sites may play a role in the development of lineage-specific phenotypes during evolution. Thus, the cases identified in this study may provide a useful resource for molecular functional analyses in search of human traits acquired during evolution.
Human N-glycosylation site data
The N-glycosylation sites in human proteins were obtained from the UniProt database (as of 13 November, 2013). The feature table of the UniProt records was scanned to collect entries with experimentally identified N-glycosylation sites. Specifically, the lines starting with “FT” followed by the “CARBOHYD” tag were examined for whether they contained a term “N-linked (GlcNAc…)”, which would indicate that the protein was N-glycosylated. Sites without experimental evidence, labeled as “potential”, “by similarity”, “partial”, or “probable”, were excluded. As a result, 2,534 N-glycosylation sites from 1,027 human proteins were obtained.
Mammalian orthologous proteins
Mammalian orthologs of the human glycosylated proteins were obtained from the UCSC Genome Browser Database (http://genome.ucsc.edu). The “CDS FASTA alignment from multiple alignments” data, derived from the “multiz100way” alignment data prepared from 100 vertebrate genomes , were downloaded using the Table Browser tool of the UCSC Genome Browser. Protein sequences from 62 mammalian species were extracted from these alignment datasets. The selected mammalian species include humans, 12 other Euarchonta species (chimpanzees, gorillas, orangutans, gibbons, rhesus macaques, crab-eating macaques, baboons, green monkeys, marmosets, squirrel monkeys, bushbabies, and treeshrews), 13 Glires species (lesser Egyptian jerboas, prairie voles, Chinese hamsters, golden hamsters, mice, rats, naked mole-rats, guinea pigs, chinchillas, brush-tailed rats, rabbits, and pikas), 25 Laurasiatheria species (pigs, alpacas, Bactrian camels, dolphins, killer whales, Tibetan antelopes, cows, sheep, goats, horses, white rhinoceroses, cats, dogs, ferrets, pandas, Pacific walruses, Weddell seals, black flying-foxes, megabats, David’s myotis bats, microbats, big brown bats, hedgehogs, shrews, and star-nosed moles), six Afrotheria species (elephants, cape elephant shrews, manatees, cape golden moles, tenrecs, and aardvarks), one Xenarthra species (armadillos), three Marsupialia species (opossums, Tasmanian devils, and wallabies), and one Monotremata species (platypuses). Additional file 1 contains detailed information on species and genome assemblies.
Computational screening for candidate novel N-glycosylation sites in human proteins
The total number of experimentally identified N-glycosylation sites collected from human proteins was 2,534. To identify mammalian proteins that were orthologous to each of the human N-glycosylated proteins, the “multiz100way” alignment data containing 57,289 alignment sets were analyzed (see Figure 1 for the overall procedure). There were 1,027 orthologous protein datasets comprising 2,534 human N-glycosylation sites. From each dataset, sequences of 62 mammalian species were extracted and realigned using MUSCLE (http://www.drive5.com/muscle) . Then, each modification site in the alignment was analyzed, and cases where more than 30% of non-Euarchonta species had an N-glycosylation motif, which might represent ancestrally conserved sites, were discarded; cases where only a small number of sequences were aligned were also discarded. A total of 184 sites in 130 protein alignments were retained after this computational screening step and subjected to in-depth semimanual inspection.
Manual inspection to select novel N-glycosylation sites in human proteins
As the final step, extensive manual inspection and curation on the 184 candidate sites was carried out to identify highly plausible cases of gains of N-glycosylation sites in the human lineage. Datasets showing the following conditions were filtered out: cases where the human sequence of UniProt database was different from that of UCSC in three or more amino acid sequence residues because of a possible paralogous relationship; cases where the N-glycosylated site was different from the consensus Asn-X-Ser/Thr; cases where the human N-glycosylation occurred only in a rare variant or mutant allele; or cases where the chimpanzee protein sequence was not included. In each dataset, sequences containing many gaps in alignment were removed from the dataset to retain only high quality sequences.
As the final result, 112 novel N-glycosylation sites in 91 human proteins were identified. Then, multiple alignments were constructed to determine when the N-glycosylation motifs first appeared. The possible functional consequences of the novel N-glycosylation site were then assessed by comprehensive literature survey and sequence analysis.
Molecular evolutionary analysis
Full-length protein and coding sequences of APMAP, CD166, and thyroglobulin were collected from humans, chimpanzees, gorillas, orangutans, gibbons, and rhesus macaques. Human cDNA RefSeq sequences were obtained from the National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov/refseq): accession numbers are [NCBI:NM_020531] (APMAP), [NCBI:NM_001627] (CD166), and [NCBI:NM_003235] (thyroglobulin). For chimpanzee, gorilla, orangutan, gibbon, and rhesus macaque orthologs, genome assemblies (panTro4, gorGor3, ponAbe2, nomLeu3, and rheMac3, respectively) were searched using each of human cDNA sequences at the UCSC Genome Browser Database. Exons, which were predicted from genomic segments, were assembled into a virtual cDNA and then conceptually translated to get a protein sequence. Some exons, which were missing in the current genome assembly, were obtained by assembling whole genome shotgun reads by searching the NCBI Sequence Read Archive (SRA) with SRA-BLAST server (http://www.ncbi.nlm.nih.gov/sra).
The ratio of nonsynonymous to synonymous substitution rates (dN/dS, ω) was estimated by a likelihood method implemented in the codeml program of the PAML package (version 4.8a) . To detect possible accelerated evolution in human proteins, we employed “branch models” that allow the ω ratio to vary among branches in phylogeny ; M0 (one ω ratio for all lineages), free ratio (one ω ratio for each branch), and two ratio (ω1 for the human branch and ω0 for the other branches). To infer positively selected sites in human proteins, we used “branch-site models” that allow the ω ratio to vary among both sites and lineages [26,28]; model A (ω is left to vary) and null model A (ω is fixed to 1). To compare the fit of nested models, the likelihood ratio test (LRT) was performed . P values were obtained using the “chi2” program in the PAML package. Protein and coding sequences, tree files, control files, and major result files for APMAP, CD166, and thyroglobulin are provided in Additional files 4, 5, and 6, respectively.
Adipocyte plasma membrane-associated protein
Cluster of differentiation 166
Cluster of differentiation 6
Thyroid-stimulating hormone receptor
Schwarz F, Aebi M. Mechanisms and principles of N-linked protein glycosylation. Curr Opin Struct Biol. 2011;21(5):576–82.
Helenius A, Aebi M. Intracellular functions of N-linked glycans. Science. 2001;291(5512):2364–9.
Dennis JW, Nabi IR, Demetriou M. Metabolism, cell surface organization, and disease. Cell. 2009;139(7):1229–41.
Scott H, Panin VM. The role of protein N-glycosylation in neural transmission. Glycobiology. 2014;24(5):407–17.
Aly AM, Higuchi M, Kasper CK, Kazazian Jr HH, Antonarakis SE, Hoyer LW. Hemophilia A due to mutations that create new N-glycosylation sites. Proc Natl Acad Sci U S A. 1992;89(11):4933–7.
Vogt G, Chapgier A, Yang K, Chuzhanova N, Feinberg J, Fieschi C, et al. Gains of glycosylation comprise an unexpectedly large group of pathogenic mutations. Nat Genet. 2005;37(7):692–700.
Vogt G, Vogt B, Chuzhanova N, Julenius K, Cooper DN, Casanova JL. Gain-of-glycosylation mutations. Curr Opin Genet Dev. 2007;17(3):245–51.
Winterpacht A, Hilbert K, Stelzer C, Schweikardt T, Decker H, Segerer H, et al. A novel mutation in FGFR-3 disrupts a putative N-glycosylation site and results in hypochondroplasia. Physiol Genomics. 2000;2(1):9–12.
Wujek P, Kida E, Walus M, Wisniewski KE, Golabek AA. N-glycosylation is crucial for folding, trafficking, and stability of human tripeptidyl-peptidase I. J Biol Chem. 2004;279(13):12827–39.
Mazumder R, Morampudi KS, Motwani M, Vasudevan S, Goldman R. Proteome-wide analysis of single-nucleotide variations in the N-glycosylation sequon of human genes. PLoS One. 2012;7(5):e36212.
Hahn Y, Lee B. Identification of nine human-specific frameshift mutations by comparative analysis of the human and the chimpanzee genome sequences. Bioinformatics. 2005;21 Suppl 1:i186–194.
Kim DS, Wang Y, Oh HJ, Lee K, Hahn Y. Frequent loss and alteration of the MOXD2 gene in catarrhines and whales: a possible connection with the evolution of olfaction. PLoS One. 2014;9(8):e104085.
Kim DS, Hahn Y. Identification of human-specific transcript variants induced by DNA insertions in the human genome. Bioinformatics. 2011;27(1):14–21.
Kim DS, Hahn Y. Human-specific protein isoforms produced by novel splice sites in the human genome after the human-chimpanzee divergence. BMC Bioinformatics. 2012;13:299.
Kim DS, Hahn Y. Identification of novel phosphorylation modification sites in human proteins that originated after the human-chimpanzee divergence. Bioinformatics. 2011;27(18):2494–501.
Kim DS, Hahn Y. Gains of ubiquitylation sites in highly conserved proteins in the human lineage. BMC Bioinformatics. 2012;13:306.
Chen R, Jiang X, Sun D, Han G, Wang F, Ye M, et al. Glycoproteomics analysis of human liver tissue by combination of multiple enzyme digestion and hydrazide chemistry. J Proteome Res. 2009;8(2):651–61.
Liu T, Qian WJ, Gritsenko MA, Camp 2nd DG, Monroe ME, Moore RJ, et al. Human plasma N-glycoproteome analysis by immunoaffinity subtraction, hydrazide chemistry, and mass spectrometry. J Proteome Res. 2005;4(6):2070–80.
Wollscheid B, Bausch-Fluck D, Henderson C, O'Brien R, Bibel M, Schiess R, et al. Mass-spectrometric identification and relative quantification of N-linked cell surface glycoproteins. Nat Biotechnol. 2009;27(4):378–86.
Gnad F, Gunawardena J, Mann M. PHOSIDA 2011: the posttranslational modification database. Nucleic Acids Res. 2011;39(Database issue):D253–260.
International Human Genome Sequencing Consortium. Finishing the euchromatic sequence of the human genome. Nature. 2004;431(7011):931–45.
Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, et al. The sequence of the human genome. Science. 2001;291(5507):1304–51.
Mangan ME, Williams JM, Kuhn RM, Lathe 3rd WC. The UCSC Genome Browser: What every molecular biologist should know. Curr Protoc Mol Biol. 2014;107:19 19 11–36.
Williams R, Ma X, Schott RK, Mohammad N, Ho CY, Li CF, et al. Encoding asymmetry of the N-glycosylation motif facilitates glycoprotein evolution. PLoS One. 2014;9(1):e86088.
Blanchette M, Kent WJ, Riemer C, Elnitski L, Smit AF, Roskin KM, et al. Aligning multiple genomic sequences with the threaded blockset aligner. Genome Res. 2004;14(4):708–15.
Dasmeh P, Serohijos AW, Kepp KP, Shakhnovich EI. Positively selected sites in cetacean myoglobins contribute to protein stability. PLoS Comput Biol. 2013;9(3):e1002929.
Yang Z. Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998;15(5):568–73.
Yang Z, dos Reis M. Statistical properties of the branch-site test of positive selection. Mol Biol Evol. 2011;28(3):1217–28.
Yang Z, Wong WS, Nielsen R. Bayes empirical Bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005;22(4):1107–18.
Albrektsen T, Richter HE, Clausen JT, Fleckner J. Identification of a novel integral plasma membrane protein induced during adipocyte differentiation. Biochem J. 2001;359(Pt 2):393–402.
Ilhan A, Gartner W, Nabokikh A, Daneva T, Majdic O, Cohen G, et al. Localization and characterization of the novel protein encoded by C20orf3. Biochem J. 2008;414(3):485–95.
Bogner-Strauss JG, Prokesch A, Sanchez-Cabo F, Rieder D, Hackl H, Duszka K, et al. Reconstruction of gene association network reveals a transmembrane protein required for adipogenesis and targeted by PPARγ. Cell Mol Life Sci. 2010;67(23):4049–64.
Mekenkamp LJ, Haan JC, Koopman M, Vink-Borger ME, Israeli D, Teerenstra S, et al. Chromosome 20p11 gains are associated with liver-specific metastasis in patients with colorectal cancer. Gut. 2013;62(1):94–101.
Mosser S, Alattia JR, Dimitrov M, Matz A, Pascual J, Schneider BL, et al. The adipocyte differentiation protein APMAP is an endogenous suppressor of Aβ production in the brain. Hum Mol Genet. 2015;24(2):371–82.
Bowen MA, Patel DD, Li X, Modrell B, Malacko AR, Wang WC, et al. Cloning, mapping, and characterization of activated leukocyte-cell adhesion molecule (ALCAM), a CD6 ligand. J Exp Med. 1995;181(6):2213–20.
Bowen MA, Bajorath J, D'Egidio M, Whitney GS, Palmer D, Kobarg J, et al. Characterization of mouse ALCAM (CD166): the CD6-binding domain is conserved in different homologs and mediates cross-species binding. Eur J Immunol. 1997;27(6):1469–78.
Weidle UH, Eggle D, Klostermann S, Swart GW. ALCAM/CD166: cancer-related issues. Cancer Genomics Proteomics. 2010;7(5):231–43.
Jannie KM, Stipp CS, Weiner JA. ALCAM regulates motility, invasiveness, and adherens junction formation in uveal melanoma cells. PLoS One. 2012;7(6):e39330.
Thelen K, Jaehrling S, Spatz JP, Pollerberg GE. Depending on its nano-spacing, ALCAM promotes cell attachment and axon growth. PLoS One. 2012;7(12):e40493.
Thelen K, Maier B, Faber M, Albrecht C, Fischer P, Pollerberg GE. Translation of the cell adhesion molecule ALCAM in axonal growth cones - regulation and functional importance. J Cell Sci. 2012;125(Pt 4):1003–14.
Malthiery Y, Lissitzky S. Primary structure of human thyroglobulin deduced from the sequence of its 8448-base complementary DNA. Eur J Biochem. 1987;165(3):491–8.
van de Graaf SA, Ris-Stalpers C, Pauws E, Mendive FM, Targovnik HM, de Vijlder JJ. Up to date with human thyroglobulin. J Endocrinol. 2001;170(2):307–21.
Yen PM. Physiological and molecular basis of thyroid hormone action. Physiol Rev. 2001;81(3):1097–142.
Mallet B, Lejeune PJ, Baudry N, Niccoli P, Carayon P, Franc JL. N-glycans modulate in vivo and in vitro thyroid hormone synthesis. Study at the N-terminal domain of thyroglobulin. J Biol Chem. 1995;270(50):29881–8.
Gagneux P, Amess B, Diaz S, Moore S, Patel T, Dillmann W, et al. Proteomic comparison of human and great ape blood plasma reveals conserved glycosylation and differences in thyroid hormone metabolism. Am J Phys Anthropol. 2001;115(2):99–109.
Previc FH. Thyroid hormone production in chimpanzees and humans: implications for the origins of human intelligence. Am J Phys Anthropol. 2002;118(4):402–3.
Overgaard MT, Sorensen ES, Stachowiak D, Boldt HB, Kristensen L, Sottrup-Jensen L, et al. Complex of pregnancy-associated plasma protein-A and the proform of eosinophil major basic protein. Disulfide structure and carbohydrate attachment. J Biol Chem. 2003;278(4):2106–17.
Laursen LS, Overgaard MT, Soe R, Boldt HB, Sottrup-Jensen L, Giudice LC, et al. Pregnancy-associated plasma protein-A (PAPP-A) cleaves insulin-like growth factor binding protein (IGFBP)-5 independent of IGF: implications for the mechanism of IGFBP-4 proteolysis by PAPP-A. FEBS Lett. 2001;504(1–2):36–40.
Lawrence JB, Oxvig C, Overgaard MT, Sottrup-Jensen L, Gleich GJ, Hays LG, et al. The insulin-like growth factor (IGF)-dependent IGF binding protein-4 protease secreted by human fibroblasts is pregnancy-associated plasma protein-A. Proc Natl Acad Sci U S A. 1999;96(6):3149–53.
Kalousova M, Muravska A, Zima T. Pregnancy-associated plasma protein A (PAPP-A) and preeclampsia. Adv Clin Chem. 2014;63:169–209.
Zhang Z, Dai H, Yu Y, Yang J, Chen J, Wu L. Elevated pregnancy-associated plasma protein A predicts myocardial dysfunction and death in severe sepsis. Ann Clin Biochem. 2014;51(Pt 1):22–9.
Farid NR, Szkudlinski MW. Minireview: structural and functional evolution of the thyrotropin receptor. Endocrinology. 2004;145(9):4048–57.
Nagayama Y, Nishihara E, Namba H, Yamashita S, Niwa M. Identification of the sites of asparagine-linked glycosylation on the human thyrotropin receptor and studies on their role in receptor function and expression. J Pharmacol Exp Ther. 2000;295(1):404–9.
Burchell B, Coughtrie MW. UDP-glucuronosyltransferases. Pharmacol Ther. 1989;43(2):261–89.
Rowland A, Miners JO, Mackenzie PI. The UDP-glucuronosyltransferases: their role in drug metabolism and detoxification. Int J Biochem Cell Biol. 2013;45(6):1121–32.
Nakajima M, Koga T, Sakai H, Yamanaka H, Fujiwara R, Yokoi T. N-glycosylation plays a role in protein folding of human UGT1A9. Biochem Pharmacol. 2010;79(8):1165–72.
Gong QH, Cho JW, Huang T, Potter C, Gholami N, Basu NK, et al. Thirteen UDP-glucuronosyltransferase genes are encoded at the human UGT1 gene complex locus. Pharmacogenetics. 2001;11(4):357–68.
Cunnane SC, Crawford MA. Survival of the fattest: fat babies were the key to evolution of the large human brain. Comp Biochem Physiol A Mol Integr Physiol. 2003;136(1):17–26.
Varki A, Altheide TK. Comparing the human and chimpanzee genomes: searching for needles in a haystack. Genome Res. 2005;15(12):1746–58.
Mihelic M, Turk D. Two decades of thyroglobulin type-1 domain research. Biol Chem. 2007;388(11):1123–30.
Li S, Iakoucheva LM, Mooney SD, Radivojac P. Loss of post-translational modification sites in disease. Pac Symp Biocomput. 2010;15:337–47.
Zielinska DF, Gnad F, Wisniewski JR, Mann M. Precision mapping of an in vivo N-glycoproteome reveals rigid topological and sequence constraints. Cell. 2010;141(5):897–907.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.
This work was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2012R1A1B3001513), Republic of Korea.
The authors declare that they have no competing interests.
YH conceived of the study. DSK and YH conducted the analysis. DSK and YH prepared the manuscript. Both authors read and approved the final manuscript.
List of mammalian species and genome assemblies.
List of novel N-glycosylation sites.
Sequence alignments of novel N-glycosylation sites.
Molecular evolutionary analysis of APMAP.
Molecular evolutionary analysis of CD166.
Molecular evolutionary analysis of thyroglobulin.
About this article
Cite this article
Kim, D.S., Hahn, Y. The acquisition of novel N-glycosylation sites in conserved proteins during human evolution. BMC Bioinformatics 16, 29 (2015). https://doi.org/10.1186/s12859-015-0468-5