A genome wide study of genetic adaptation to high altitude in feral Andean Horses of the páramo
© Hendrickson; licensee BioMed Central Ltd. 2013
Received: 19 July 2013
Accepted: 3 December 2013
Published: 17 December 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.
Keywordshypoxia Equus Hypoxia Inducing Factor Adaptation Genomes wide association
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.
Genome wide allelic differentiation
Most significant allelic divergences between feral Andean horses and breed out-groups
7.90 × 10-14
7.90 × 10-14
3.01 × 10-8
1.14 × 10-13
1.48 × 10-11
2.02 × 10-13
2.19 × 10-8
1.36 × 10-9
2.08 × 10-12
2.42 × 10-10
3.48 × 10-12
1.33 × 10-8
3.59 × 10-11
4.73 × 10-11
5.67 × 10-11
TENM2 and other nervous system genes
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 .
Neurological system process genes represented in a Gene ontology biological processes analysis (GO:0050877, p = 5.1 × 10 -5 )
Entrez gene ID
C1q and tumor necrosis factor related protein 5
solute carrier family 12 (sodium/potassium/chloride transporters), member 2
endothelial PAS domain protein 1
Kv channel interacting protein 1
nuclear receptor subfamily 4, group A, member 3
glutamate receptor, ionotropic, delta 2
ubiquitin protein ligase E3 component n-recognin 3 (putative)
MAS-related GPR, member X2
olfactory receptor, family 13, subfamily C, member 5
olfactory receptor, family 13, subfamily F, member 1
olfactory receptor, family 13, subfamily C, member 4
olfactory receptor, family 13, subfamily C, member 3
olfactory receptor, family 13, subfamily C, member 8
olfactory receptor, family 13, subfamily C, member 2
olfactory receptor, family 51, subfamily Q, member 1
Cytochrome P450 genes
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
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 .
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=.
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.
- 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.Google Scholar
- 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.Google Scholar
- Simpson GG: Splendid Isolation: The Curious History of South American Mammals. 1980, New Haven, CT: New Haven Yale University PressGoogle Scholar
- McFadden B: Fossil Horses. 1992, Cambridge: Cambridge University PressGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Zhou F, You G: Molecular insights into the structure-function relationship of organic anion transporters OATs. Pharm Res. 2007, 24: 28-36.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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 issuePubMed CentralPubMedView ArticleGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.