- Research article
- Open Access
A genome wide study of genetic adaptation to high altitude in feral Andean Horses of the páramo
BMC Evolutionary Biology volume 13, Article number: 273 (2013)
Life at high altitude results in physiological and metabolic challenges that put strong evolutionary pressure on performance due to oxidative stress, UV radiation and other factors dependent on the natural history of the species. To look for genes involved in altitude adaptation in a large herbivore, this study explored genome differentiation between a feral population of Andean horses introduced by the Spanish in the 1500s to the high Andes and their Iberian breed relatives.
Using allelic genetic models and Fst analyses of ~50 K single nucleotide polymorphisms (SNPs) across the horse genome, 131 candidate genes for altitude adaptation were revealed (Bonferoni of p ≤ 2 × 10–7). Significant signals included the EPAS1 in the hypoxia-induction-pathway (HIF) that was previously discovered in human studies (p = 9.27 × 10-8); validating the approach and emphasizing the importance of this gene to hypoxia adaptation. Strong signals in the cytochrome P450 3A gene family (p = 1.5 ×10-8) indicate that other factors, such as highly endemic vegetation in altitude environments are also important in adaptation. Signals in tenuerin 2 (TENM2, p = 7.9 × 10-14) along with several other genes in the nervous system (gene categories representation p = 5.1 × 10-5) indicate the nervous system is important in altitude adaptation.
In this study of a large introduced herbivore, it becomes apparent that some gene pathways, such as the HIF pathway are universally important for high altitude adaptation in mammals, but several others may be selected upon based on the natural history of a species and the unique ecology of the altitude environment.
In addition to hypoxia and high levels of ultraviolet radiation, many aspects of high altitude environments put strong evolutionary pressures on resident species. Extreme temperature and humidity fluctuation, highly endemic vegetation and fauna, and other biological factors relative to a species’ natural history, such as length or timing of the breeding season or level of population isolation can create strong adaptive pressure. Understanding patterns and timing of genetic adaptations through organisms of different life histories is important to elucidate the commonalities in unique adaptation pathways that occur in extreme environments.
The feral horses of the high Andean pàramo were originally brought with the Spanish conquistadors in the 1500s. Although most herds have been extirpated, small groups now called the párameros or cerreros persist in Ecuador’s eastern range between the Cotopaxi volcano and Quilindaña. These horses came predominantly from Andalusia, and had a mixture of Jennet, Andalusia, and Berber ancestry. They quickly adapted to the Andes, interbreeding naturally in small bands in isolation , with the consequence that they present an ideal natural experiment in adaptation to high altitude conditions.
In addition to hypoxia, the highly inclement pàramo presents other physiologic and metabolic challenges. The pàramo is cold and humid with yearly rainfalls between 500 to over 3000 mm. Extreme temperature fluctuations, with diurnal ranges from below freezing to up to 30C and periods of fog, heavy rain, hail or even snow dominate a good portion of the day during the rainy season. The dry season is very short. Regular subzero temperatures, high UV radiation, and low pH soils have selected for an highly endemic alpine vegetation primarily consisting of tussock grasses, ground rosettes, dwarf shrubs, cushion plants and giant rosettes such as Espeletia and Puya . Although Equus fossils have been found throughout the Andes, all New World Equids went extinct in the Americas at the end of the Pleistocene  and in the high Andes, the remaining large ungulates (i.e. Camelids) shifted their center of distribution to the more temperate puna grasslands to the south . The success of the introduced horses in the northern pàramo provides a unique model for looking at extreme hypoxic, thermoregulatory and metabolic adaptive pressure over a relatively short time span of ~200 generations.
Several recent genome wide studies have focused on human adaptation to hypoxia. Comparisons of high-altitude Tibetan populations with low-altitude Han Chinese populations uncovered strong signatures of selection in the hypoxia inducing factor (HIF) response pathway in the region of endothelial PAS domain 1 (EPAS1) [5–7]. Several other genes in this pathway were also associated including peroxi-some proliferator-activated receptor-a (PPARA) (Tibetans) and egl 9 homolog 1 (EGLN1) (Tibetans and Andeans [5–9]), and the protein kinase, AMP-activated a 1 catalytic subunit (PRKAA1) (Andeans [10, 11]). Selection in the HIF pathway is detected in most human studies; however, several other pathways such as methylation, pathogen resistance, or DNA damage repair do not clearly overlap between study groups . This could reflect differences in evolutionary paths, time scale, or other natural history variables that have not been elucidated yet.
As a large-bodied herbivore, the introduced Andean horse provides a unique comparative perspective for further understanding whether the HIF pathway is a universal target of natural selection for hypoxia, and how time, altitudinal ecosystem, and species life history influence genetic adaptation to extreme environments. This study explores which genomic regions appear to be under selection in the Andean horse, and how these genes relate to previous studies and to the natural history of the horse and the páramo.
Results and discussion
The feral population of horses is made up of several small bands located on lands adjacent to the Cotopaxi National Park, including Hacienda Yanahurco, a private wildlife reserve of ~7300 acres, and other unpopulated land to the south. The populations size fluctuates from ~100 to over 200 individuals due to extreme weather fluctuations and removal of animals by wranglers.
Horses are small in stature and many individuals exhibit the heritable “gaited” trait found in breeds such as the Paso fino or Mangalarga Marchador. A principal component analysis indicated that the feral horses form a cluster close to these two breeds in PC1 (Figure 1A), and they are a distinct cluster as defined by PC3 (Figure 1B). Several breeds form tight clusters (e.g. Spanish Barb Breed Association, North American Peruvian horses, Galicinos), indicating either under strong artificial selection and/or a small number of founders. The general arrangement of the clusters reflected relationships illustrated in a larger prior analysis of breeds .
Genome wide allelic differentiation
Using Allelic Genetic models for association, 129 SNPs in the feral Andean horses show genetic divergence from the mixed breed sample below the Bonferoni cut-off (p < 2 × 10-7) as shown in the Manhattan plot in Figure 2 (p-values given in Additional file 1: Table S2). 131 genes are located within the regions identified by these SNPs that may have been under strong selection once horses were introduced into the Andes. Table 1 lists the most significant ten associations, and each region of interest is discussed in turn below.
TENM2 and other nervous system genes
The most significant divergence in allele frequencies in the Andean horse is found in the TENM2 gene, a member of the teneurin family of type II tramsmembrane glycoproteins (p = 7.9 × 10–14, Figure 3). This signal is also supported by Fst analysis (Figure 3B). The genes for teneurins are highly conserved across invertebrate and vertebrate species. The region on chromosome 14 in Equus shares 93.7% sequence identity with human TENM2, and 81.7% with chicken TENM2 (formerly ODZ2).
Teneurin proteins are proposed to regulate gene expression in the nervous system during development . TEMN2 is involved in development of neuronal circuits in the visual system  and has been shown to be expressed in the developing limbs, somites, and craniofacial mesenchyme .
Further support for the nervous system as a focus of adaptation in the feral Andean horse population is that "neurological system process" was the most significantly represented category in a Gene Ontology Biological Processes analysis (GO:0050877, p = 5.1 × 10-5) (Additional file 2: Table S3). Eighteen of 87 genes categorized in GO fell into this category (Table 2), including EPAS1. Hypoxia is implicated in several human central nervous system pathologies such as stroke or neurodegenerative disease and a connection between the HIF pathway and neuronal response to hypoxia has been recognized . 12 of the genes were in sensory perception, particularly olfaction, which may indicate that better sensory abilities are a benefit in the wild where individuals must locate mates, maintain herd structure, avoid predators and find appropriate food sources. In addition developmental and functional genes within the nervous system were significant.
Cytochrome P450 genes
Significant frequency divergences occur on chromosome 13 in a region containing several cytochrome P450 (CYP) 3A genes (Figure 4). Three SNPs show significant allele frequency divergence meeting the Bonferoni correction of p <2 × 10-7 and three additional SNPS meet a lower cut-off of 10-6 (Bonferoni of p < 0.05, BIEC2-207561, BIEC2-207603, BIEC2-207605). The SNP BIEC2-207449 is located within CYP3A93; however the entire horse CYP3A gene cluster consists of seven genes and one pseudogene in strong linkage disequilibrium covering a region of ~115 Kb. Fst analyses also supported this region as a candidate of selection (Figure 4B).
Cytochrome P450 enzymes are a superfamily of membrane bound heme containing monooxygenases whose function is to catalyze the oxidation of organic substances. Best known for their role in drug metabolism, CYPs also are involved in hormone synthesis and breakdown, cholesterol synthesis, vitamin D metabolism, and metabolism of toxic compounds. The horse CYP3A gene cluster is more extensive than the human CYP cluster. It has been suggested that the horse CYP gene cluster is more complex because horses must digest a range of plant toxins in different habitats . In the páramo ~ 60% of plant species are endemic and include largely bunchgrasses, shrubs, mosses and lichens which are highly evolved, likely in response to high UV, low pH soil and temperature and humidity fluctuation. This vegetation would differ greatly from that originally found in the habitats of the founding horses. It should be noted that members of this gene family were found in genome-wide hypoxia adaptation studies of Tibetan human populations  and frogs  and experimental data from rabbits indicate hypoxemia up-regulates the expression of CYP3A6 .
Replication of associations between EPAS1 and life at high altitude
EPAS1 has been found to be associated with altitude in several studies of Tibetans [5–7, 21]. EPAS1 encodes a transcription factor involved in the induction of genes regulated by oxygen, which is induced as oxygen levels fall. In the horse, BIEC2–310909 (rs69041973) was highly significant (p = 9.27 × 10-8) as shown in Figure 5. This SNP is an intronic SNP with no known function. It has been suggested that in Tibetans, the association between EPAS1 expression and lower Hb levels is an adaptive response that lowers hypoxia-induced erythropoeisis in Tibetans such that the negative effects of this hypoxia response do not occur [6, 22]. A more in-depth follow-up study will be needed to confirm the result and elucidate the exact mechanism for the difference, but the commonality of this gene association in horses with the studies in humans is intriguing.
Other genes with high allele frequency differences within the Andean herd include LINGO2, which was recently associated with BMI  and Parkinson disease [24, 25], and a region containing ubiquitin protein ligase E3 component n-recognin 3 (UBR3), Sp5 transcription factor, glutamate decarboxylase 1 (GAD1) and myosin IIIB (MYO3B). Solute carrier family 22 (organic anion/urate transporter), member 11 (SLC22A11) (p = 5.7 × 10–8) is part of the SLC22A gene family that mediates the absorption, distribution, and excretion of a diverse array of environmental toxins .
Most Fst signals reflected the significant association signals. The highest region of allele frequency divergence as measured by Fst was observed on chromosome 11 (Figure 6), and contained the following genes: TOM1L1, COX11, STXBP4, HLF, MMD, TMEM100, and PCTP. The Fst data was corroborated by a highly significant association found between COX11/STXBP4 (p = 1.44 × 10-9). These two genes have been explored for their association to human breast cancer in several studies [27–30]. They are in strong LD in both humans and horses. COX11 is a nuclear encoded mitochondrial protein and the terminal component of the mitochondrial respiratory chain. COX11 has been found to be up-regulated in chronic hypoxia of cyanotic patients undergoing repair of heart defects, suggesting a role in dealing with lack of oxygen, possibly by acting as an heme biosynthetic enzyme which transports copper to heme A [31, 32].
Several large allelic divergences between the feral Andean horses and a comparison group of primarily Iberian horse breeds have been observed in this study. Many of the genes identified in this study appear to be biologically relevant, however one cannot rule out type 1 error due to genetic drift or the composition of the comparison group. A recent study of artificial selection in 33 horse breeds found strong selection for aesthetics and performance traits resulted in high homozygosity within breeds and large divergence between breeds . The Andean population was likely founded by a mixed sample of Iberian breeds; however, by combining the descendants of those breeds in the comparison group, many of which have undergone strong subsequent artificial selection since their introduction to the Americas, one effectively makes a comparison to a more heterogeneous out-group. This may increase the potential for false positives. The replication of the significance of the EPAS1 loci suggests that many genes found in this study are the result of strong selection in the wild population. The replication of this signal in a large herbivore in a relatively short evolutionary time span further implicates HIF- response pathway as essential to high altitude. Additional signals in the nervous system and the cytochrome-P450 gene family, which may be in response to the local endemic plants, suggest that high altitude adaptation is complicated, and that life history of the species in question and local ecology of the environment are both important factors in evolutionary adaptation.
Population samples from the feral horse herd (n = 97) were collected on the southeastern side of the Cotopaxi National Park in Ecuador (0°41´28.97¨ S 78°16´44.69¨ W) in November 2006 and March 2007 during the annual “round-up” which was started in the late 1990s. DNA in this study is from pulled hair because it is non-evasive and can be collected rapidly; however, blood and small skin clippings were collected when possible and preserved in RNALater (Ambion, Inc). No samples were taken from foals born during the current season to avoid excess stress. The round-up employs Chagras from adjacent villages who ride horses originally captured from the feral herd in previous years, therefore samples were also obtained from riding and pack horses.
For a comparison group, registered breed samples (N = 55) were collected in the United States and South America primarily from Iberian breeds (Andalusian, Lipizzaner) or breeds that have recently been defined and founded in the Americas from original Spanish-mixed stock; including Colonial Spanish Horse, Sorraia Mustang, American Paint Horse, American Heritage Horse, Horse of the Americas, Galiceno, North American Peruvian Horse, Paso Fino, Quarter Horse, US Mangalarga Marchador and Spanish Barb (details given in Additional file 3: Table S1). More distant breed samples included Arabian (Bask and Polish), Shetland Pony, Miniature Horse, Thoroughbred, and Sandohit Oldenburg. The Andean populations were likely founded from horses from more than one Iberian breed. Further, a recent study of artificial selection showed that many breeds included in this study have been under strong selection for aesthetic and performance traits resulting in high level of homogeneity within breeds . Artificial selection would have continued after the Iberian breeds were established in the Americas. Therefore breeds were combined to lesson artifacts from artificial selection and to provide a more heterogeneous out-group similar to the potential founding population.
DNA isolation and genotyping
DNA was isolated with the Qiagen DNeasy Blood and Tissue Kit (Valencia, CA) following the manufacturer’s protocol with the exception of over-night digestion at 56C for all samples. The Equine SNP50 BeadChip (Illumina, Inc., San Diego, CA) was used to genotype 50,023 SNPs across the horse genome in 151 horses from the Andean population and related breeds. Average spacing between functional SNPs on the chip is 43.1kB with lower coverage on chromosome X (49.44kB) . The average genotyping completeness in this study was 99.7% (92.0–99.9%).
Initial quality control (QC) was performed in Genome Studio version 2010.3 (Illumina, Inc., San Diego, CA). Of the 54,602 SNPs typed, 2200 were excluded due to <60% Gentrain Scores (n = 917), <95% Call frequency (n = 625), or >3 clusters or low intensity scores (n = 657). No SNPs had >2 Mendel Errors based on the 3 trios included in the genotyping. 108 SNPs were excluded with missingness >5% between case and controls. Of the remaining 52,294, 50,023 had a minor allele frequency of >0.01 and were included in subsequent analyses. Of these, 584 (0.01%) were monomorphic within Andeans.
Estimations of relatedness and population structure
Because the samples in this study were collected during a round–up of an unmarked herd over a large montane area, it was necessary to estimate relatedness between samples, and then remove all first order relatives from population analyses. Kinship coefficients were based on 33,483 SNPs pruned to remove those in high linkage disequilibrium (independent pairwise pruning, window size = 50, r2 = 0.5) in PLINK version 1.07 . 18 feral individuals had a kinship coefficient > 0.5, of which 8 individuals were removed from subsequent analyses. Kinship coefficients also identified 9 offspring of an Arabian stallion introduced to the riding herd adjacent to the hacienda in 2006, which were not included in downstream analyses. Kinship coefficient distributions from the samples collected indicate the population is not extremely inbred (Additional file 4: Figure S1).
The population structure of the feral Andean horses and several related breeds were analyzed by principle component analyses in Eigensoft 3.0  with LD-pruned dataset (32, 672 SNPs). These data confirmed the half Arabian offspring as well 6 horses originally from the wild herd that were now owned by local Chagras.
Genome–wide statistical analyses
Selection was analyzed on two levels. First a SNP-by-SNP allelic model was performed in PLINK (version 1.07)  between feral Andean horses (n = 78) and the comparison group to look for significant differences in allele frequency between the horses that have been living at altitude (3500-4500 m) for the past 500 years and founding breeds. P–values were based on the Fisher’s Exact test. Allelic models for all SNPs in the analyses are shown in Additional file 5. A Bonferroni correction of p < 2 × 10–7 used to define genome-wide significance.
A second analysis using sliding window Fst was applied to look for signals of selection  within each chromosome with window sizes of 100, 50 and 25 SNPs, which were on average 4.7 ± 0.2, 2.3 ±0.2 and 1.2 ±0.2 Mb in size respectively.
Genes within statistically significant regions were identified in the Equus assembly equCab2.0 (2007) within the UCSC browser [38, 39]. Genes were included if they were within ±10Kb from the significant SNP. In cases where large blocks of LD surrounded the SNP, or several adjacent SNPs were significant, these areas were included in for gene identification. Human and other available reference assemblies were aligned to the horse assembly to assist in identification of candidate genes. A GO functional analysis was performed on the genes found with the GWAS method with Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 . P-values were based on the Fisher’s Exact test.
Availability of supporting data
Genotyping data can be found at https://mynotebook.labarchives.com/share/Horse%2520SNP%2520chip%2520/My45fDI0OTgzLzMvVHJlZU5vZGUvMjQ0NjE5MjkxN3w5Ljk=.
Moscoso LS: Horses and Tack. Chagras: Ecuadors Andean Coyboys. Edited by: Oxford P, Bish R, Quito EC. 2004, Imprenta Mariscal, Quito Ecuador: The Nature Conservancy, 61-90.
Luteyn JL: Paramos: Why Study Them?. Páramo an Andean Ecosyst Under Hum Influ. Edited by: Balslev H, Luteyn JL. 1992, San Diego, CA, USA: Academic Press, 1-14.
Simpson GG: Splendid Isolation: The Curious History of South American Mammals. 1980, New Haven, CT: New Haven Yale University Press
McFadden B: Fossil Horses. 1992, Cambridge: Cambridge University Press
Beall CM, Cavalleri GL, Deng L, Elston RC, Gao Y, Knight J, Li C, Li JC, Liang Y, McCormack M, Montgomery HE, Pan H, Robbins PA, Shianna KV, Tam SC, Tsering N, Veeramah KR, Wang W, Wangdui P, Weale ME, Xu Y, Xu Z, Yang L, Zaman MJ, Zeng C, Zhang L, Zhang X, Zhaxi P, Zheng YT: Natural selection on EPAS1 (HIF2alpha) associated with low hemoglobin concentration in Tibetan highlanders. Proc Natl Acad Sci U S A. 2010, 107: 11459-11464. 10.1073/pnas.1002443107.
Peng Y, Yang Z, Zhang H, Cui C, Qi X, Luo X, Tao X, Wu T, Chen H, Shi H, Su B, Basang , Ouzhuluobu , Danzengduojie , Ciwangsangbu : Genetic variations in Tibetan populations and high-altitude adaptation at the Himalayas. Mol Biol Evol. 2011, 28: 1075-1081. 10.1093/molbev/msq290.
Yi X, Liang Y, Huerta-Sanchez E, Jin X, Cuo ZXP, Pool JE, Xu X, Jiang H, Vinckenbosch N, Korneliussen TS, Zheng HH, Liu T, He W, Li K, Luo R, Nie X, Wu H, Zhao M, Cao H, Zou J, Shan Y, Li SS, Yang Q, Yi X, Liang Y, Huerta-Sanchez E, Jin X, Cuo ZXP, Pool JE, Xu X, et al: Sequencing of 50 human exomes reveals adaptation to high altitude. Science. 2010, 329 (Zhao M): 75-78.
Simonson TS, Yang Y, Huff CD, Yun H, Qin G, Witherspoon DJ, Bai Z, Lorenzo FR, Xing J, Jorde LB, Prchal JT, Ge R: Genetic evidence for high-altitude adaptation in Tibet. Science. 2010, 329: 72-75. 10.1126/science.1189406.
Wang B, Zhang Y-B, Zhang F, Lin H, Wang X, Wan N, Ye Z, Weng H, Zhang L, Li X, Yan J, Wang P, Wu T, Cheng L, Wang J, Wang D-M, Ma X, Yu J: On the origin of Tibetans and their genetic basis in adapting high-altitude environments. PLoS One. 2011, 6: e17002-10.1371/journal.pone.0017002.
Bigham A, Bauchet M, Pinto D, Mao X, Akey JM, Mei R, Scherer SW, Julian CG, Wilson MJ, López Herráez D, Brutsaert T, Parra EJ, Moore LG, Shriver MD: Identifying signatures of natural selection in Tibetan and Andean populations using dense genome scan data. PLoS Genet. 2010, 6: 1-14.
Bigham AW, Mao X, Mei R, Brutsaert T, Wilson MJ, Julian CG, Parra EJ, Akey JM, Moore LG, Shriver MD: Identifying positive selection candidate loci for high-altitude adaptation in Andean populations. Hum Genomics. 2009, 4: 79-90.
Alkorta-Aranburu G, Beall CM, Witonsky DB, Gebremedhin A, Pritchard JK, Di Rienzo A: The genetic architecture of adaptations to high altitude in Ethiopia. PLoS Genet. 2012, 8: e1003110-10.1371/journal.pgen.1003110.
Petersen JL, Mickelson JR, Cothran EG, Andersson LS, Axelsson J, Bailey E, Bannasch D, Binns MM, Borges AS, Brama P, da Câmara Machado A, Distl O, Felicetti M, Fox-Clipsham L, Graves KT, Guérin G, Haase B, Hasegawa T, Hemmann K, Hill EW, Leeb T, Lindgren G, Lohi H, Lopes MS, McGivney BA, Mikko S, Orr N, Penedo MCT, Piercy RJ, Raekallio M, et al: Genetic Diversity in the Modern Horse Illustrated from Genome-Wide SNP Data. PLoS One. 2013, 8: e54997-10.1371/journal.pone.0054997.
Tucker RP, Chiquet-Ehrismann R: Teneurins: a conserved family of transmembrane proteins involved in intercellular signaling during development. Dev Biol. 2006, 290: 237-245. 10.1016/j.ydbio.2005.11.038.
Rubin BP, Tucker RP, Brown-Luedi M, Martin D, Chiquet-Ehrismann R: Teneurin 2 is expressed by the neurons of the thalamofugal visual system in situ and promotes homophilic cell-cell adhesion in vitro. Development. 2002, 129: 4697-4705.
Tucker RP, Chiquet-Ehrismann R, Chevron MP, Martin D, Hall RJ, Rubin BP: Teneurin-2 is expressed in tissues that regulate limb and somite pattern formation and is induced in vitro and in situ by FGF8. Dev Dyn. 2001, 220: 27-39. 10.1002/1097-0177(2000)9999:9999<::AID-DVDY1084>3.0.CO;2-B.
Corcoran A, O’Connor JJ: Hypoxia inducible factor signalling mechanisms in the central nervous system. Acta Physiol (Oxf). 2013, 208: 298-310. 10.1111/apha.12117.
Schmitz A, Demmel S, Peters LM, Leeb T, Mevissen M, Haase BLA: Comparative human-horse sequence analysis of the CYP3A subfamily gene cluster. Anim Genet. 2010, 41 (Suppl 2): 72-79.
Yang W, Qi Y, Bi K, Fu J: Toward understanding the genetic basis of adaptation to high-elevation life in poikilothermic species: a comparative transcriptomic analysis of two ranid frogs, Rana chensinensis and R. kukunoris. BMC Genomics. 2012, 13: 588-10.1186/1471-2164-13-588.
Fradette C, Batonga J, Teng S, Piquette-Miller M, du Souich P: Animal models of acute moderate hypoxia are associated with a down-regulation of CYP1A1, 1A2, 2B4, 2C5, and 2C16 and up-regulation of CYP3A6 and P-glycoprotein in liver. Drug Metab Dispos. 2007, 35: 765-771. 10.1124/dmd.106.013508.
Xu S, Li S, Yang Y, Tan J, Lou H, Jin W, Yang L, Pan X, Wang J, Shen Y, Wu B, Wang H, Jin L: A genome-wide search for signals of high-altitude adaptation in Tibetans. Mol Biol Evol. 2011, 28: 1003-1011. 10.1093/molbev/msq277.
Storz JF, Scott GR, Cheviron ZA: Phenotypic plasticity and genetic adaptation to high-altitude hypoxia in vertebrates. J Exp Biol. 2010, 213 (Pt 24): 4125-4136.
Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, Jackson AU, Lango Allen H, Lindgren CM, Luan J, Mägi R, Randall JC, Vedantam S, Winkler TW, Qi L, Workalemahu T, Heid IM, Steinthorsdottir V, Stringham HM, Weedon MN, Wheeler E, Wood AR, Ferreira T, Weyant RJ, Segrè AV, Estrada K, Liang L, Nemesh J, Park J-H, Gustafsson S, Kilpeläinen TO, et al: Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet. 2010, 42: 937-948. 10.1038/ng.686.
Vilariño-Güell C, Wider C, Ross OA, Jasinska-Myga B, Kachergus J, Cobb SA, Soto-Ortolaza AI, Behrouz B, Heckman MG, Diehl NN, Testa CM, Wszolek ZK, Uitti RJ, Jankovic J, Louis ED, Clark LN, Rajput A, Farrer MJ: LINGO1 and LINGO2 variants are associated with essential tremor and Parkinson disease. Neurogenetics. 2010, 11: 401-408. 10.1007/s10048-010-0241-x.
Wu Y-W, Prakash KM, Rong T-Y, Li H-H, Xiao Q, Tan LC, Au W-L, Ding J, Chen S, Tan E-K: Lingo2 variants associated with essential tremor and Parkinson’s disease. Hum Genet. 2011, 129: 611-615. 10.1007/s00439-011-0955-3.
Zhou F, You G: Molecular insights into the structure-function relationship of organic anion transporters OATs. Pharm Res. 2007, 24: 28-36.
Tang L, Xu J, Wei F, Wang L, Nie W-W, Chen L-B, Guan X-X: Association of STXBP4/COX11 rs6504950 (G>A) polymorphism with breast cancer risk: evidence from 17,960 cases and 22,713 controls. Arch Med Res. 2012, 43: 383-388. 10.1016/j.arcmed.2012.07.008.
Campa D, Kaaks R, Le Marchand L, Haiman CA, Travis RC, Berg CD, Buring JE, Chanock SJ, Diver WR, Dostal L, Fournier A, Hankinson SE, Henderson BE, Hoover RN, Isaacs C, Johansson M, Kolonel LN, Kraft P, Lee I-M, McCarty CA, Overvad K, Panico S, Peeters PHM, Riboli E, Sanchez MJ, Schumacher FR, Skeie G, Stram DO, Thun MJ, Trichopoulos D, et al: Interactions between genetic variants and breast cancer risk factors in the breast and prostate cancer cohort consortium. J Natl Cancer Inst. 2011, 103: 1252-1263. 10.1093/jnci/djr265.
Antoniou AC, Beesley J, McGuffog L, Sinilnikova OM, Healey S, Neuhausen SL, Ding YC, Rebbeck TR, Weitzel JN, Lynch HT, Isaacs C, Ganz PA, Tomlinson G, Olopade OI, Couch FJ, Wang X, Lindor NM, Pankratz VS, Radice P, Manoukian S, Peissel B, Zaffaroni D, Barile M, Viel A, Allavena A, Dall’Olio V, Peterlongo P, Szabo CI, Zikan M, Claes K, et al: Common breast cancer susceptibility alleles and the risk of breast cancer for BRCA1 and BRCA2 mutation carriers: implications for risk prediction. Cancer Res. 2010, 70: 9742-9754. 10.1158/0008-5472.CAN-10-1907.
Bhatti P, Doody MM, Rajaraman P, Alexander BH, Yeager M, Hutchinson A, Burdette L, Thomas G, Hunter DJ, Simon SL, Weinstock RM, Rosenstein M, Stovall M, Preston DL, Linet MS, Hoover RN, Chanock SJ, Sigurdson AJ: Novel breast cancer risk alleles and interaction with ionizing radiation among U.S. radiologic technologists. Radiat Res. 2010, 173: 214-224. 10.1667/RR1985.1.
Hiser L, Di Valentin M, Hamer AG, Hosler JP: Cox11p is required for stable formation of the Cu(B) and magnesium centers of cytochrome c oxidase. J Biol Chem. 2000, 275: 619-623. 10.1074/jbc.275.1.619.
Ghorbel MT, Cherif M, Jenkins E, Mokhtari A, Kenny D, Angelini GD, Caputo M: Transcriptomic analysis of patients with tetralogy of Fallot reveals the effect of chronic hypoxia on myocardial gene expression. J Thorac Cardiovasc Surg. 2010, 140: 37–345-e26.
Petersen JL, Mickelson JR, Rendahl AK, Valberg SJ, Andersson LS, Axelsson J, Bailey E, Bannasch D, Binns MM, Borges AS, Brama P, da Câmara MA, Capomaccio S, Cappelli K, Cothran EG, Distl O, Fox-Clipsham L, Graves KT, Guérin G, Haase B, Hasegawa T, Hemmann K, Hill EW, Leeb T, Lindgren G, Lohi H, Lopes MS, McGivney BA, Mikko S, Orr N, et al: Genome-wide analysis reveals selection for important traits in domestic horse breeds. PLoS Genet. 2013, 9: e1003211-10.1371/journal.pgen.1003211.
McCue ME, Bannasch DL, Petersen JL, Gurr J, Bailey E, Binns MM, Distl O, Guérin G, Hasegawa T, Hill EW, Leeb T, Lindgren G, Penedo MCT, Røed KH, Ryder OA, Swinburne JE, Tozaki T, Valberg SJ, Vaudin M, Lindblad-Toh K, Wade CM, Mickelson JR: A high density SNP array for the domestic horse and extant Perissodactyla: utility for association mapping, genetic diversity, and phylogeny studies. PLoS Genet. 2012, 8: e1002451-10.1371/journal.pgen.1002451.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira M A R, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC): 559-575.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38: 904-909. 10.1038/ng1847.
Weir BS, Cardon LR, Anderson AD, Nielsen DM, Hill WG: Measures of human population structure show heterogeneity among genomic regions. Genome Res. 2005, 15: 1468-1476. 10.1101/gr.4398405.
Wade CM, Giulotto E, Sigurdsson S, Zoli M, Gnerre S, Imsland F, Lear TL, Adelson DL, Bailey E, Bellone RR, Blöcker H, Distl O, Edgar RC, Garber M, Leeb T, Mauceli E, MacLeod JN, Penedo MCT, Raison JM, Sharpe T, Vogel J, Andersson L, Antczak DF, Biagi T, Binns MM, Chowdhary BP, Coleman SJ, Della Valle G, Fryc S, Guérin G, et al: Genome sequence, comparative analysis, and population genetics of the domestic horse. Science. 2009, 326: 865-867. 10.1126/science.1178158.
Dreszer TR, Karolchik D, Zweig AS, Hinrichs AS, Raney BJ, Kuhn RM, Meyer LR, Wong M, Sloan CA, Rosenbloom KR, Roe G, Rhead B, Pohl A, Malladi VS, Li CH, Learned K, Kirkup V, Kirkup V, Hsu F, Harte RA, Guruvadoo L, Goldman M, Giardine BM, Fujita PA, Diekhans M, Cline MS, Clawson H, Barber GP, Haussler D, James Kent W: The UCSC Genome Browser database: extensions and updates 2011. Nucleic Acids Res. 2012, 40: D918-D923. 10.1093/nar/gkr1055. Database issue
Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.
This study was funded by the Morris Animal Foundation and in part by the Laboratory of Genomic Diversity at the National Cancer Institute, and Shepherd University and the West Virginia Space Grant Consortium.
Fernando Cobo of Hacienda Yanahurco provided access to the lands and field support through both his own work and that of his Chagras. Eric Hendrickson-Lambert provided much needed field assistance. Juan Carlos and Lilly Matheus (F.E.I.C.E.) assisted with logistics and permits in Ecuador. Juan Carlos died unexpectedly during this study—a great loss to science. Laboratory assistance was provided by Joanna Martinez and Taylor Donley through the Werner Student Intern Program and Cancer Research Training Awards (lab of SLH), and Michael Malasky of the National Cancer Institute. Russ Hanson and Janice Martensen arranged U.S. permits. Gila Kahila Bar-Gal, Sagi Polani, Marilyn Menotti-Raymond and Victor David gave valuable discussions. Stephen J. O’Brien provided laboratory and scholarly support.
The following people provided breed samples: Haley L. Kostinas of Horizon Equine Vet Service LLC, James Kehler, Diane Gates, Marye Ann Thompson, Ed Ogg, Michaele Srock, Janice Brennan, Barbara Windom, Lucia Roda, Silke Schneider, Amy Heuter-Imming, Amy Jo Maurer, Simrat Khalsa, Karen Moore, James R. Dildine, John O. Banion, Richard M. Blaney, Anne Beggs, Mike And Sheri Olson, and Pat Miller. Two anonymous reviewers provided valuable comments.
The author declares that she has no competing interests.
Electronic supplementary material
Additional file 1: Table S2: SNPs found to be associated with the Andean horse population association using the allelic model at p<0.05 (Bonferoni p<1 ×10-6). Additional analyses were limited to p<0.01 (Bonferoni p<2 ×10-7). * = SNP located in gene. (XLS 269 KB)
Additional file 3: Table S1: Horse samples included in the study. No.=The number of individuals within that group. (DOCX 64 KB)
Kinship coefficients distribution for the Andean horses collected during the study based on 33,483 LD-pruned SNPs (independent pairwise pruning, window size = 50, r
Additional file 4: 2 = 0.5, PLINK version 1.07).(PNG 24 KB)
Additional file 5: Table S4: All results from association tests between the Andean horse population and the outgroup using the allelic model. Fisher’s p-value. A1=Allele 1, A2= Allele 2. (CSV 2 MB)
Authors’ original submitted files for images
About this article
Cite this article
Hendrickson, S.L. A genome wide study of genetic adaptation to high altitude in feral Andean Horses of the páramo. BMC Evol Biol 13, 273 (2013). https://doi.org/10.1186/1471-2148-13-273
- Hypoxia Inducing Factor
- Genomes wide association