Distribution of events of positive selection and population differentiation in a metabolic pathway: the case of asparagine N-glycosylation

BMC Evolutionary Biology

Jaume Bertranpetit²  Ludovica Montanucci²  Martin Sikora¹  Pierre Luisi²  Hafid Laayouni²  Giovanni Marco Dall’Olio² 

[1]Department of Genetics, Stanford University School of Medicine, Stanford, USA

[2]IBE, Institut de Biologia Evolutiva (UPF-CSIC), Parc de Recerca Biomèdica de Barcelona (PRBB), Dr. Aiguader, 88, 08003, Barcelona, Catalonia, Spain

关键词: Adaptation to environment;    Calnexin/calreticulin cycle;    Pathway analysis;    Glycosylation;    Asparagine N-Glycosylation;    Population differentiation;    Positive selection;    Homo sapiens;   

Others  :  1140966
DOI  :  10.1186/1471-2148-12-98

 received in 2011-12-31, accepted in 2012-06-25,  发布年份 2012

【 摘 要 】

Background

Asparagine N-Glycosylation is one of the most important forms of protein post-translational modification in eukaryotes. This metabolic pathway can be subdivided into two parts: an upstream sub-pathway required for achieving proper folding for most of the proteins synthesized in the secretory pathway, and a downstream sub-pathway required to give variability to trans-membrane proteins, and involved in adaptation to the environment and innate immunity. Here we analyze the nucleotide variability of the genes of this pathway in human populations, identifying which genes show greater population differentiation and which genes show signatures of recent positive selection. We also compare how these signals are distributed between the upstream and the downstream parts of the pathway, with the aim of exploring how forces of population differentiation and positive selection vary among genes involved in the same metabolic pathway but subject to different functional constraints.

Results

Our results show that genes in the downstream part of the pathway are more likely to show a signature of population differentiation, while events of positive selection are equally distributed among the two parts of the pathway. Moreover, events of positive selection are frequent on genes that are known to be at bifurcation points, and that are identified as being in key position by a network-level analysis such as MGAT3 and GCS1.

Conclusions

These findings indicate that the upstream part of the Asparagine N-Glycosylation pathway has lower diversity among populations, while the downstream part is freer to tolerate diversity among populations. Moreover, the distribution of signatures of population differentiation and positive selection can change between parts of a pathway, especially between parts that are exposed to different functional constraints. Our results support the hypothesis that genes involved in constitutive processes can be expected to show lower population differentiation, while genes involved in traits related to the environment should show higher variability. Taken together, this work broadens our knowledge on how events of population differentiation and of positive selection are distributed among different parts of a metabolic pathway.

【 授权许可】

   
2012 Dall'Olio et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

【 参考文献 】

• [1]Wright S: The genetical structure of populations. Ann Eugen 1951, 15:323-353.
• [2]Weir BS, Hill WG: Estimating F-statistics. Annu Rev Genet 2002, 36:721-750.
• [3]Voight BF, Kudaravalli S, Wen X, Pritchard JK: A map of recent positive selection in the human genome. PLoS Biol 2006, 4(3):e72.
• [4]Grossman I, Avidan N, Singer C, Paperna T, Lancet D, Beckmann JS, Miller A: Genomic profiling of interpopulation diversity guides prioritization of candidate-genes for autoimmunity. Genes Immun 2004, 5:493-504.
• [5]Ryan AW, Mapp J, Moyna S, Mattiangeli V, Kelleher D, Bradley DG, McManus R: Levels of interpopulation differentiation among different functional classes of immunologically important genes. Genes Immun 2006, 7:179-183.
• [6]Caramelo JJ, Parodi AJ: Getting in and out from calnexin/calreticulin cycles. J Biol Chem 2008, 283:10221-10225.
• [7]Parodi AJ: Protein glucosylation and its role in protein folding. Annu Rev Biochem 2000, 69:69-93.
• [8]Lederkremer GZ: Glycoprotein folding, quality control and ER-associated degradation. Curr Opin Struct Biol 2009, 19:515-523.
• [9]Shental-Bechor D, Levy Y: Folding of glycoproteins: toward understanding the biophysics of the glycosylation code. Curr Opin Struct Biol 2009, 19:524-533.
• [10]Roth J, Yam GH, Fan J, et al.: Protein quality control: the who’s who, the where’s and therapeutic escapes. Histochem Cell Biol 2008, 129:163-177.
• [11]Stanley P, Schachter H, Taniguchi N: N-Glycans. Cold Spring Harbor Laboratory Press, In Essentials of Glycobiology; 2009.
• [12]Oriol R, Martinez-Duncker I, Chantret I, Mollicone R, Codogno P: Common origin and evolution of glycosyltransferases using Dol-P-monosaccharides as donor substrate. Mol Biol Evol 2002, 19:1451-1463.
• [13]Haeuptle MA, Hennet T: Congenital disorders of glycosylation: an update on defects affecting the biosynthesis of dolichol-linked oligosaccharides. Hum Mutat 2009, 30:1628-1641.
• [14]Jaeken J: Congenital disorders of glycosylation. Ann N Y Acad Sci 2010, 1214:190-198.
• [15]Murrell MP, Yarema KJ, Levchenko A: The systems biology of glycosylation. Chembiochem: a European journal of chemical biology 2004, 5:1334-1347.
• [16]Hossler P, Goh LT, Lee MM, Hu WS: GlycoVis: Visualizing Glycan Distribution in the Protein N-Glycosylation Pathway in Mammalian Cells. Biotechnol Bioeng 2006, 95:946-960.
• [17]Hossler P, Mulukutla BC, Hu WS: Systems analysis of N-glycan processing in mammalian cells. PLoS One 2007, 2(1):e713.
• [18]Kim PM, Korbel JO, Gerstein MB: Positive selection at the protein network periphery: evaluation in terms of structural constraints and cellular context. Proc Natl Acad Sci USA 2007, 104:20274-20279.
• [19]Dennis JW, Lau KS, Demetriou M, Nabi IR: Adaptive regulation at the cell surface by N-glycosylation. Traffic 2009, 10:1569-1578.
• [20]Dennis JW, Nabi IR, Demetriou M: Metabolism, cell surface organization, and disease. Cell 2009, 139:1229-1241.
• [21]Varki A, Freeze H: Glycans in Acquired Human Diseases. Cold Spring Harbor Laboratory Press, In Essentials of Glycobiology; 2009.
• [22]Harduin-Lepers A, Mollicone R, Delannoy P, Oriol R: The animal sialyltransferases and sialyltransferase-related genes: a phylogenetic approach. Glycobiology 2005, 15:805-817.
• [23]Varki A: Nothing in glycobiology makes sense, except in the light of evolution. Cell 2006, 126:841-845.
• [24]Gagneux P, Varki A: Evolutionary considerations in relating oligosaccharide diversity to biological function. Glycobiology 1999, 9:747-755.
• [25]Montanucci L, Laayouni H, Dall’Olio GM, Bertranpetit J: Molecular evolution and network-level analysis of the N-glycosylation metabolic pathway across primates. Mol Biol Evol 2011, 28:813-823.
• [26]Scardoni G, Petterlini M, Laudanna C: Analyzing biological network parameters with CentiScaPe. Bioinformatics (Oxford, England) 2009, 25:2857-2859.
• [27]Peng G, Luo L, Siu H, et al.: Gene and pathway-based second-wave analysis of genome-wide association studies. European journal of human genetics: EJHG 2009, 18:111-117.
• [28]Song Y, Aglipay JA, Bernstein JD, Goswami S, Stanley P: The bisecting GlcNAc on N-glycans inhibits growth factor signaling and retards mammary tumor progression. Cancer Res 2010, 70:3361-3371.
• [29]Makino T, Suzuki Y, Gojobori T: Differential evolutionary rates of duplicated genes in protein interaction network. Gene 2006, 385:57-63.
• [30]Vitkup D, Kharchenko P, Wagner A: Influence of metabolic network structure and function on enzyme evolution. Genome Biol 2006, 7:R39. BioMed Central Full Text
• [31]Yang YH, Zhang FM, Ge S: Evolutionary rate patterns of the Gibberellin pathway genes. BMC Evol Biol 2009, 9:206. BioMed Central Full Text
• [32]Alvarez-Ponce D, Aguadé M, Rozas J: Network-level molecular evolutionary analysis of the insulin/TOR signal transduction pathway across 12 Drosophila genomes. Genome Res 2009, 19:234-242.
• [33]Hahn MW, Conant GC, Wagner A: Molecular evolution in large genetic networks: does connectivity equal constraint? J Mol Evol 2004, 58:203-211.
• [34]Luisi P, Alvarez-Ponce D, Dall'olio GM, Sikora M, Bertranpetit J, Laayouni H: Network-level and population genetics analysis of the insulin/TOR signal transduction pathway across human populations. Mol Biol Evol 2011. Epub ahead of print
• [35]Teshima KM, Coop G, Przeworski M: How reliable are empirical genomic scans for selective sweeps? Genome Res 2006, 16:702-12.
• [36]Pickrell JK, Coop G, Novembre J, et al.: Signals of recent positive selection in a worldwide sample of human populations. Genome Res 2009, 19:826-37.
• [37]Li JZ, Absher DM, Tang H, et al.: Worldwide human relationships inferred from genome-wide patterns of variation. Science 2008, 319:1100-4.
• [38]Cann H, de Toma C, Cazes L, et al.: A human genome diversity cell line panel. Science 2002, 296:261-2.
• [39]Rosenberg N: Standardized subsets of the HGDP-CEPH Human Genome Diversity Cell Line Panel, accounting for atypical and duplicated samples and pairs of close relatives. Ann Hum Genet 2006, 70:841-7.
• [40]Scheet P, Stephens M: A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet 2006, 78:629-644.
• [41]Dall’Olio GM, Jassal B, Montanucci L, Gagneux P, Bertranpetit J, Laayouni H: The annotation of the Asparagine N-linked Glycosylation pathway in the Reactome Databas. Glycobiology 2011, 21:1395-1400.
• [42]Stajich JE, Block D, Boulez K, et al.: The Bioperl toolkit: Perl modules for the life sciences. Genome Res 2002, 12:1611-8.
• [43]Kelley JL, Madeoy J, Calhoun JC, Swanson W, Akey JM: Genomic signatures of positive selection in humans and the limits of outlier approaches. Genome Res 2006, 16:980-989.
• [44]iHS software download. http://hgdp.uchicago.edu/Software webcite
• [45]Cline MS, Smoot M, Cerami E, et al.: Integration of biological networks and gene expression data using Cytoscape. Nat Protoc 2007, 2:2366-82.
• [46]Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T: Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics 2011, 27:431-432.
• [47]Cytoscape session file used in this analysis. https://bitbucket.org/dalloliogm/pathways/src/tip/cytoscape webcite
• [48]Script used to plot report on centralities. https://bitbucket.org/dalloliogm/cytoscape-compare-node-centralities webcite
• [49]Repository for all code used in this analysis. https://bitbucket.org/dalloliogm/pathways webcite