Landscape genomics analysis reveals the genetic basis underlying cashmere goats and dairy goats adaptation to frigid environments

Abstract

Understanding the genetic mechanism of cold adaptation in cashmere goats and dairy goats is very important to improve their production performance. The purpose of this study was to comprehensively analyze the genetic basis of goat adaptation to cold environments, clarify the impact of environmental factors on genome diversity, and lay the foundation for breeding goat breeds to adapt to climate change. A total of 240 dairy goats were subjected to genome resequencing, and the whole genome sequencing data of 57 individuals from 6 published breeds were incorporated. By integrating multiple approaches such as phylogenetic analysis, population structure analysis, gene flow and population history exploration, selection signal analysis, and genome-environment association analysis, an in-depth investigation was carried out. Phylogenetic analysis unraveled the genetic relationships and differentiation patterns among dairy goats and other goat breeds. Through signal analysis (θπ, FST, XP-CLR), we identified numerous candidate genes associated with cold adaptation in dairy goats (STRIP1, ALX3, HTR4, NTRK2, MRPL11, PELI3, DPP3, BBS1) and cashmere goats (MED12L, MARC2, MARC1, DSG3, C6H4orf22, CHD7, MYPN, KIAA0825, MITF). Genome-environment association (GEA) analysis confirmed the link between these genes and environmental factors. Moreover, a detailed analysis of the critical genes C6H4orf22 and STRIP1 demonstrated their significant roles in the geographical variations of cold adaptation and allele frequency differences among different breeds. This study contributes to understanding the genetic basis of cold adaptation, providing crucial theoretical support for precision breeding programs aimed at improving production performance in cold regions by leveraging adaptive alleles, thereby ensuring sustainable animal husbandry.

Keywords

Genome-environment association / Environmental adaptation / Frigid environments / Whole genome resequencing / Select signal analysis

References

[1] AbbasiT, AbbasiSA. Reducing the global environmental impact of livestock production: the minilivestock option. J Clean Prod Prev Smog Crises, 2016, 112: 1754-1766.

[2] AlexanderDH, LangeK. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinformatics, 2011, 12246.

[3] AmillsM, CapoteJ, Tosser-KloppG. Goat domestication and breeding: a jigsaw of historical, biological and molecular data with missing pieces. Anim Genet, 2017, 48: 631-644.

[4] Amiri GhanatsamanZ, Ayatolahi MehrgardiA, Asadollahpour NanaeiH, EsmailizadehA. Comparative genomic analysis uncovers candidate genes related with milk production and adaptive traits in goat breeds. Sci Rep, 2023, 138722.

[5] AriantiR, VinnaiBÁ, GyőryF, GubaA, CsőszÉ, KristófE, FésüsL. Availability of abundant thiamine determines efficiency of thermogenic activation in human neck area derived adipocytes. J Nutr Biochem, 2023, 119. 109385

[6] AzzuV, JastrochM, DivakaruniAS, BrandMD. The regulation and turnover of mitochondrial uncoupling proteins. Biochim Biophys Acta, 2010, 1797: 785-791.

[7] BenjellounB, LeempoelK, BoyerF, StuckiS, StreeterI, Orozco-terWengelP, AlbertoFJ, ServinB, BiscariniF, AlbertiA, EngelenS, StellaA, ColliL, CoissacE, BrufordMW, Ajmone-MarsanP, NegriniR, ClarkeL, FlicekP, ChikhiA, JoostS, TaberletP, PompanonF. Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints. Mol Ecol, 2024, 33. e17257

[8] BolgerAM, LohseM, UsadelB. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014, 30: 2114-2120.

[9] BuD, LuoH, HuoP, WangZ, ZhangS, HeZ, WuY, ZhaoL, LiuJ, GuoJ, FangS, CaoW, YiL, ZhaoY, KongL. KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res, 2021, 49: W317-W325.

[10] CayeK, JumentierB, LepeuleJ, FrançoisO. LFMM 2: Fast and Accurate Inference of Gene-Environment Associations in Genome-Wide Studies. Mol Biol Evol, 2019, 36: 852-860.

[11] ChangL, ZhengY, LiS, NiuX, HuangS, LongQ, RanX, WangJ. Identification of genomic characteristics and selective signals in Guizhou black goat. BMC Genomics, 2024, 25164.

[12] ChenH, PattersonN, ReichD. Population differentiation as a test for selective sweeps. Genome Res, 2010, 20: 393-402.

[13] CoumouD, RahmstorfS. A decade of weather extremes. Nature Clim Change, 2012, 2: 491-496.

[14] Danecek, P., Auton, A., Abecasis, G., Albers, C.A., Banks, E., DePristo, M.A., Handsaker, R.E., Lunter, G., Marth, G.T., Sherry, S.T., McVean, G., Durbin, R., 1000 Genomes Project Analysis Group. The variant call format and VCFtools. Bioinformatics, 2011, 27: 2156-2158.

[15] DanmaigoroA, MuhammadMA, AbubakarK, MagiriRB, BakareAG, IjiPA. Morphological and physiological features in small ruminants: an adaptation strategy for survival under changing climatic conditions. Int J Biometeorol, 2024, 68: 1497-1505.

[16] DansoF, IddrisuL, LunguSE, ZhouG, JuX. Effects of Heat Stress on Goat Production and Mitigating Strategies: A Review. Animals, 2024, 141793.

[17] Li C, Cao H, Ren Y, Jia J, Yang G, Jin J, Shi X (2024) Eicosapentaenoic acid-mediated activation of PGAM2 regulates skeletal muscle growth and development via the PI3K/AKT pathway. Int J Biol Macromol 268:131547. https://doi.org/10.1016/j.ijbiomac.2024.131547

[18] ExcoffierL, DupanloupI, Huerta-SánchezE, SousaVC, FollM. Robust demographic inference from genomic and SNP data. PLoS Genet, 2013, 9. e1003905

[19] FilomenaMC, YamamotoDL, CaremaniM, KadarlaVK, MastrototaroG, SerioS, VydyanathA, MutarelliM, GarofaloA, PerticiI, KnöllR, NigroV, LutherPK, LieberRL, BeckMR, LinariM, BangM-L. Myopalladin promotes muscle growth through modulation of the serum response factor pathway. J Cachexia Sarcopenia Muscle, 2020, 11: 169-194.

[20] FrichotE, SchovilleSD, BouchardG, FrançoisO. Testing for Associations between Loci and Environmental Gradients Using Latent Factor Mixed Models. Mol Biol Evol, 2013, 30: 1687-1699.

[21] Gaitskell-PhillipsG, Martín-CanoFE, Ortiz-RodríguezJM, Silva-RodríguezA, da Silva-ÁlvarezE, GilMC, Ortega-FerrusolaC, PeñaFJ. The seminal plasma proteins Peptidyl arginine deaminase 2, rRNA adenine N (6)-methyltransferase and KIAA0825 are linked to better motility post thaw in stallions. Theriogenology, 2022, 177: 94-102.

[22] GoddeCM, Mason-D’CrozD, MayberryDE, ThorntonPK, HerreroM. Impacts of climate change on the livestock food supply chain; a review of the evidence. Glob Food Sec, 2021, 28. 100488

[23] GuoH, XiaL, WangW, XuW, ShenX, WuX, HeT, JiangX, XuY, ZhaoP, TanD, ZhangX, ZhangY. Hypoxia induces alterations in tRNA modifications involved in translational control. BMC Biol, 2023, 2139.

[24] GuoJ, ZhongJ, LiL, ZhongT, WangL, SongT, ZhangH. Comparative genome analyses reveal the unique genetic composition and selection signals underlying the phenotypic characteristics of three Chinese domestic goat breeds. Genet Sel Evol, 2019, 5170.

[25] Habeeb, A.A., Osman, S.F., Teama, F.E.I., Gad, A.E.. The detrimental impact of high environmental temperature on physiological response, growth, milk production, and reproductive efficiency of ruminants. Trop Anim Health Prod, 2023, 55388.

[26] HarelS, HigginsCA, CeriseJE, DaiZ, ChenJC, ClynesR, ChristianoAM. Pharmacologic inhibition of JAK-STAT signaling promotes hair growth. Sci Adv, 2015, 1. e1500973

[27] HeJ, ZhaoB, HuangX, FuX, LiuG, TianY, WuC, MaoJ, LiuJ, GunS, TianK. Gene network analysis reveals candidate genes related with the hair follicle development in sheep. BMC Genomics, 2022, 23428.

[28] HoffmannAA, WeeksAR, SgròCM. Opportunities and challenges in assessing climate change vulnerability through genomics. Cell, 2021, 184: 1420-1425.

[29] HuH, SunS-C. Ubiquitin signaling in immune responses. Cell Res, 2016, 26: 457-483.

[30] HuX-J, YangJ, XieX-L, LvF-H, CaoY-H, LiW-R, LiuM-J, WangY-T, LiJ-Q, LiuY-G, RenY-L, ShenZ-Q, WangF, HehuaEe, HanJ-L, LiM-H. The Genome Landscape of Tibetan Sheep Reveals Adaptive Introgression from Argali and the History of Early Human Settlements on the Qinghai-Tibetan Plateau. Mol Biol Evol, 2019, 36: 283-303.

[31] HwangJ, PallasDC. STRIPAK complexes: structure, biological function, and involvement in human diseases. Int J Biochem Cell Biol, 2014, 47: 118-148.

[32] IwenKA, OelkrugR, BrabantG. Effects of thyroid hormones on thermogenesis and energy partitioning. J Mol Endocrinol, 2018, 60: R157-R170.

[33] JonesAK, BajramiB, CampbellMK, ErzurumluogluAM, GuoQ, ChenH, ZhangX, ZevelevaS, KvaskoffD, BrunnerA-D, MullerS, GatheyV, DaveRM, TannerJW, RixenS, StruweMA, PhoenixK, KlumphKJ, RobinsonH, VeyelD, MullerA, NoyvertB, BartholdyBA, Steixner-KumarAA, StutzkiJ, DrichelD, OmlandS, SheehanR, HillJ, BretschneiderT, GottschlingD, ScheidigAJ, ClementB, GieraM, DingZ, BroadwaterJ, WarrenCR. mARC1 in MASLD: Modulation of lipid accumulation in human hepatocytes and adipocytes. Hepatology Communications, 2024, 8. e0365

[34] KlingensporM. Cold-induced recruitment of brown adipose tissue thermogenesis. Exp Physiol, 2003, 88: 141-148.

[35] KorneliussenTS, MoltkeI, AlbrechtsenA, NielsenR. Calculation of Tajima’s D and other neutrality test statistics from low depth next-generation sequencing data. BMC Bioinformatics, 2013, 14289.

[36] KorytinaGF, AkhmadishinaLZ, KochetovaOV, AznabaevaYG, ZagidullinShZ, VictorovaTV. Polymorphic variants of glutamate receptor (GRIK5, GRIN2B) and serotonin receptor (HTR2A) genes are associated with chronic. Mol Biol, 2017, 51: 533-542.

[37] KumarS, StecherG, LiM, KnyazC, TamuraK. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol, 2018, 35: 1547-1549.

[38] LakeIR, HooperL, AbdelhamidA, BenthamG, BoxallABA, DraperA, Fairweather-TaitS, HulmeM, HunterPR, NicholsG, WaldronKW. Climate Change and Food Security: Health Impacts in Developed Countries. Environ Health Perspect, 2012, 120: 1520-1526.

[39] LetunicI, BorkP. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res, 2021, 49: W293-W296.

[40] LiH, DurbinR. Inference of human population history from individual whole-genome sequences. Nature, 2011, 475: 493-496.

[41] LiH, DurbinR. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics, 2009, 25: 1754-1760.

[42] Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., Marth, G., Abecasis, G., Durbin, R., 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics, 2009, 25: 2078-2079.

[43] LimaARC, SilveiraRMF, CastroMSM, De VecchiLB, da R FernandesMHM, de ResendeKT. Relationship between thermal environment, thermoregulatory responses and energy metabolism in goats: A comprehensive review. J Therm Biol, 2022, 109103324.

[44] LiZ, ChenL, ChenX, ZhongY, YangY, XiaW, et al.. Biological adaptations in the Arctic cervid, the reindeer (Rangifer tarandus). Science, 2019, 364eaav6312.

[45] LiuY, ShepherdEG, NelinLD. MAPK phosphatases — regulating the immune response. Nat Rev Immunol, 2007, 7: 202-212.

[46] LobellDB, SchlenkerW, Costa-RobertsJ. Climate Trends and Global Crop Production Since 1980. Science, 2011, 333: 616-620.

[47] Ma Y, Zhu S, Lv T, Gu X, Feng H, Zhen J, Xin W, Wan Q (2020) SQSTM1/p62 Controls mtDNA Expression and Participates in Mitochondrial Energetic Adaption via MRPL12. iScience 23. https://doi.org/10.1016/j.isci.2020.101428

[48] MacHughDE, BradleyDG. Livestock genetic origins: Goats buck the trend. Proc Natl Acad Sci USA, 2001, 985382.

[49] ManzettiS, ZhangJ, Van Der SpoelD. Thiamin Function, Metabolism, Uptake, and Transport. Biochemistry, 2014, 53: 821-835.

[50] MasekM, EtardC, HofmannC, HülsmeierAJ, ZangJ, TakamiyaM, GesemannM, NeuhaussSCF, HornemannT, SträhleU, Bachmann-GagescuR. Loss of the Bardet-Biedl protein Bbs1 alters photoreceptor outer segment protein and lipid composition. Nat Commun, 2022, 131282.

[51] McKennaA, HannaM, BanksE, SivachenkoA, CibulskisK, KernytskyA, GarimellaK, AltshulerD, GabrielS, DalyM, DePristoMA. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res, 2010, 20: 1297-1303.

[52] MendozaMN, RaudseppT, AlshanbariF, GutiérrezG, Ponce de LeónFA. Chromosomal Localization of Candidate Genes for Fiber Growth and Color in Alpaca (Vicugna pacos). Front Genet, 2019, 10583.

[53] NaderiS, RezaeiH-R, PompanonF, BlumMGB, NegriniR, NaghashH-R, BalkızÖ, MashkourM, GaggiottiOE, Ajmone-MarsanP, KenceA, VigneJ-D, TaberletP. The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals. Proc Natl Acad Sci, 2008, 105: 17659-17664.

[54] NiuY, LiY, ZhaoY, HeX, ZhaoQ, PuY, MaY, JiangL. Whole-genome sequencing identifies functional genes for environmental adaptation in Chinese sheep. J Genet Genomics, 2024, 51: 1278-1285.

[55] PattersonN, PriceAL, ReichD. Population structure and eigenanalysis. PLoS Genet, 2006, 2. e190

[56] PeriasamyM, HerreraJL, ReisFCG. Skeletal Muscle Thermogenesis and Its Role in Whole Body Energy Metabolism. Diabetes Metab J, 2017, 41327.

[57] PickrellJK, PritchardJK. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet, 2012, 8. e1002967

[58] PravenecM, SabaLM, ZídekV, LandaV, MlejnekP, ŠilhavýJ, ŠimákováM, StrnadH, TrnovskáJ, ŠkopV, HüttlM, MarkováI, OliyarnykO, MalínskáH, KazdováL, SmithH, TabakoffB. Systems genetic analysis of brown adipose tissue function. Physiol Genomics, 2018, 50: 52-66.

[59] QiuQ, ZhangG, MaT, QianW, WangJ, YeZ, CaoC, HuQ, KimJ, LarkinDM, AuvilL, CapitanuB, MaJ, LewinHA, QianX, LangY, ZhouR, WangL, WangK, XiaJ, LiaoS, PanS, LuX, HouH, WangY, ZangX, YinY, MaH, ZhangJ, WangZ, ZhangY, ZhangD, YonezawaT, HasegawaM, ZhongY, LiuW, ZhangY, HuangZ, ZhangS, LongR, YangH, WangJ, LenstraJA, CooperDN, WuY, WangJ, ShiP, WangJ, LiuJ. The yak genome and adaptation to life at high altitude. Nat Genet, 2012, 44: 946-949.

[60] QuinlanAR, HallIM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics, 2010, 26: 841-842.

[61] RellstabC, GugerliF, EckertAJ, HancockAM, HoldereggerR. A practical guide to environmental association analysis in landscape genomics. Mol Ecol, 2015, 24: 4348-4370.

[62] RopiquetA, HassaninA. Hybrid origin of the Pliocene ancestor of wild goats. Mol Phylogenet Evol, 2006, 41: 395-404.

[63] RossJ, ChauM, MillerB, MukherjeeS, ZhangC, KongJ, KaissiE, NewsamA, MahboubiD, BerryJ, Tucker-BurdenC, BratD. TMIC-14. Hypoxia Inducible Olig2 Mediates Glioma Stem Cell Migration. Neuro-Oncology, 2017, 19vi246.

[64] RoulinA. Melanin-based colour polymorphism responding to climate change. Glob Change Biol, 2014, 20: 3344-3350.

[65] RowlandLA, BalNC, PeriasamyM. The role of skeletal-muscle-based thermogenic mechanisms in vertebrate endothermy. Biol Rev, 2015, 90: 1279-1297.

[66] SilvaJE. Thermogenic Mechanisms and Their Hormonal Regulation. Physiol Rev, 2006, 86: 435-464.

[67] SinhaRA, YenPM. Metabolic Messengers: Thyroid Hormones. Nat Metab, 2024, 6: 639-650.

[68] SlominskiA, TobinDJ, ShibaharaS, WortsmanJ. Melanin Pigmentation in Mammalian Skin and Its Hormonal Regulation. Physiol Rev, 2004.

[69] SmithRE. Thermoregulatory and Adaptive Behavior of Brown Adipose Tissue. Science, 1964, 146: 1686-1689.

[70] Teixeira IAMA, Härter CJ, Vargas JAC, Souza AP, Fernandes MHMR (2024) Review: Update of nutritional requirements of goats for growth and pregnancy in hot environments. Animal 18 Suppl 2:101219. https://doi.org/10.1016/j.animal.2024.101219

[71] ThorntonPK, EricksenPJ, HerreroM, ChallinorAJ. Climate variability and vulnerability to climate change: a review. Glob Change Biol, 2014, 20: 3313-3328.

[72] ThorntonPK, van de SteegJ, NotenbaertA, HerreroM. The impacts of climate change on livestock and livestock systems in developing countries: A review of what we know and what we need to know. Agric Syst, 2009, 101: 113-127.

[73] VijhRK, SharmaU, KapoorP, RahejaM, AroraR, AhlawatS, DurejaV. Design and validation of high-density SNP array of goats and population stratification of Indian goat breeds. Gene, 2023, 885. 147691

[74] VillarroyaF, CereijoR, VillarroyaJ, GiraltM. Brown adipose tissue as a secretory organ. Nat Rev Endocrinol, 2017, 13: 26-35.

[75] WangD, TianY, WangQ, ZhangY, YeB, ZuoZ, HeJ, PanZ, SunD, ZouJ, XieS. Cold stress-induced autophagy and apoptosis disorders are mainly mediated by AMPK/PPAR/PI3K/AKT/mTOR pathways. Aquaculture, 2024, 583. 740574

[76] WangK, LiM, HakonarsonH. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res, 2010, 38. e164

[77] YangR, HanZ, ZhouW, LiX, ZhangX, ZhuL, WangJ, LiX, ZhangC-L, HanY, LiL, LiuS. Population structure and selective signature of Kirghiz sheep by Illumina Ovine SNP50 BeadChip. PeerJ, 2024, 12. e17980

[78] Ye P, Duan W, Leng YQ, Wang YK, Tan X, Wang WZ (2022) DPP3: From biomarker to therapeutic target of cardiovascular diseases. Front Cardiovasc Med 9. https://doi.org/10.3389/fcvm.2022.974035

[79] ZhangC, Asadollahpour NanaeiH, Jafarpour NegariN, Amiri RoudbarM, Amiri GhanatsamanZ, NiyazbekovaZ, YangX. Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats. BMC Genomics, 2024, 25477.

[80] ZhangC, DongS-S, XuJ-Y, HeW-M, YangT-L. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics, 2019, 35: 1786-1788.

[81] ZhaoH, ZhuS, GuoT, HanM, ChenB, QiaoG, WuY, YuanC, LiuJ, LuZ, SunW, WangT, LiF, ZhangY, HouF, YueY, YangB. Whole-genome re-sequencing association study on yearling wool traits in Chinese fine-wool sheep. J Anim Sci, 2021, 99skab210.

[82] Zhao J, Mu Y, Gong P, Liu B, Zhang F, Zhu L, Shi C, Lv X, Luo J (2024) Whole-genome resequencing of native and imported dairy goat identifies genes associated with productivity and immunity. Front Vet Sci 11. https://doi.org/10.3389/fvets.2024.1409282

[83] ZhengZ, WangX, LiM, LiY, YangZ, WangX, et al.. The origin of domestication genes in goats. Sci Adv, 2020, 6eaaz5216.

[84] ZhouW, CaiH, LiJ, XuH, WangX, MenH, ZhengY, CaiL. Potential roles of mediator Complex Subunit 13 in Cardiac Diseases. Int J Biol Sci, 2021, 17: 328-338.

Funding

the Hohhot City Science and Technology Plan Project (Major Science and Technology Special Project) in China(2023150103000025)

the National Key Research and Development Program of China(2021YFD1600704)

the Tianchi Talent Special Appointed Expert Project

RIGHTS & PERMISSIONS

The Author(s)