scholarly journals Genome-wide scans for signatures of selection in Mangalarga Marchador horses using high-throughput SNP genotyping

Author(s):  
Wellington Bizarria dos Santos ◽  
Gustavo Schettini ◽  
Amanda Marchi Maiorano ◽  
Fernando Oliveira Bussiman ◽  
Júlio Carvalho Balieiro ◽  
...  

Abstract Background: The Mangalarga Marchador horse (MM) is one of the breeds shaped over generations by local adaptations and specific preferences of Brazilian breeders to morphology, functionality, and locomotion. The animals genetically have “batida” or “picada” natural gait trait, which is a trademark of the breed. The movement biomechanics of this breed promote stability during the execution, comfort, and softness of the ride. The detection of signatures of selection in genomic regions provides insights into the evolutionary process, enabling discoveries regarding complex phenotypic traits. In this research, we focused on the identification of genomic regions affected by different selection pressures, mainly highlighting recent selection, as well as understanding the candidate genes and functional pathways associated with the signatures of selection in the MM genome. A broader discussion of genes and traits of importance in this breed, especially traits related to the type and quality of the gait, temperament, conformation, and locomotor system, was also provided.Results: Three different methods were used to search for signals of selection: Tajima’s D (TD), the integrated haplotype score (iHS), and runs of homozygosity (ROH). The samples were composed of males (n=62) and females (n=130) that were initially chosen considering well-defined phenotypes for gait: picada (n=86) and batida (n=106). All horses were genotyped using a 670k Axiom ® Equine Genotyping Array​ (Axiom MNEC670). In total, 169 pruned candidate genes harboring important biological processes were found, highlighting the following: anterior/posterior pattern for the set of genes (GLI3, HOXC9, HOXC6, HOXC5, HOXC4, HOXC13, HOXC11, and HOXC10); limb morphogenesis, skeletal system, proximal/distal pattern formation, JUN kinase activity (CCL19 and MAP3K6); and muscle stretch response (MAPK14). Other candidate genes were associated with energy metabolism, bronchodilator response, NADH regeneration, reproduction, keratinization, and the immunological system.Conclusions: Collectively, our findings revealed for the first time some pieces of evidence for selection signals acting on athletic performance, gait type, and energy to muscle activity in the MM breed, with the particular involvement of the HOX family of genes. Gene network analysis showed no relationship between the signals observed in this population and the DMRT3 gene, a relevant gene associated with gait trait, however candidate genes that were located close (~28Mb) to DMRT3 were identified.

BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Wellington B. Santos ◽  
Gustavo P. Schettini ◽  
Amanda M. Maiorano ◽  
Fernando O. Bussiman ◽  
Júlio C. C. Balieiro ◽  
...  

Abstract Background The detection of signatures of selection in genomic regions provides insights into the evolutionary process, enabling discoveries regarding complex phenotypic traits. In this research, we focused on identifying genomic regions affected by different selection pressures, mainly highlighting the recent positive selection, as well as understanding the candidate genes and functional pathways associated with the signatures of selection in the Mangalarga Marchador genome. Besides, we seek to direct the discussion about genes and traits of importance in this breed, especially traits related to the type and quality of gait, temperament, conformation, and locomotor system. Results Three different methods were used to search for signals of selection: Tajima’s D (TD), the integrated haplotype score (iHS), and runs of homozygosity (ROH). The samples were composed of males (n = 62) and females (n = 130) that were initially chosen considering well-defined phenotypes for gait: picada (n = 86) and batida (n = 106). All horses were genotyped using a 670 k Axiom® Equine Genotyping Array​ (Axiom MNEC670). In total, 27, 104 (chosen), and 38 candidate genes were observed within the signatures of selection identified in TD, iHS, and ROH analyses, respectively. The genes are acting in essential biological processes. The enrichment analysis highlighted the following functions: anterior/posterior pattern for the set of genes (GLI3, HOXC9, HOXC6, HOXC5, HOXC4, HOXC13, HOXC11, and HOXC10); limb morphogenesis, skeletal system, proximal/distal pattern formation, JUN kinase activity (CCL19 and MAP3K6); and muscle stretch response (MAPK14). Other candidate genes were associated with energy metabolism, bronchodilator response, NADH regeneration, reproduction, keratinization, and the immunological system. Conclusions Our findings revealed evidence of signatures of selection in the MM breed that encompass genes acting on athletic performance, limb development, and energy to muscle activity, with the particular involvement of the HOX family genes. The genome of MM is marked by recent positive selection. However, Tajima’s D and iHS results point also to the presence of balancing selection in specific regions of the genome.


2021 ◽  
Vol 12 ◽  
Author(s):  
Pedro H. F. Freitas ◽  
Yachun Wang ◽  
Ping Yan ◽  
Hinayah R. Oliveira ◽  
Flavio S. Schenkel ◽  
...  

Understanding the biological mechanisms of climatic adaptation is of paramount importance for the optimization of breeding programs and conservation of genetic resources. The aim of this study was to investigate genetic diversity and unravel genomic regions potentially under selection for heat and/or cold tolerance in thirty-two worldwide cattle breeds, with a focus on Chinese local cattle breeds adapted to divergent climatic conditions, Datong yak (Bos grunniens; YAK), and Bali (Bos javanicus) based on dense SNP data. In general, moderate genetic diversity levels were observed in most cattle populations. The proportion of polymorphic SNP ranged from 0.197 (YAK) to 0.992 (Mongolian cattle). Observed and expected heterozygosity ranged from 0.023 (YAK) to 0.366 (Sanhe cattle; SH), and from 0.021 (YAK) to 0.358 (SH), respectively. The overall average inbreeding (±SD) was: 0.118 ± 0.028, 0.228 ± 0.059, 0.194 ± 0.041, and 0.021 ± 0.004 based on the observed versus expected number of homozygous genotypes, excess of homozygosity, correlation between uniting gametes, and runs of homozygosity (ROH), respectively. Signatures of selection based on multiple scenarios and methods (FST, HapFLK, and ROH) revealed important genomic regions and candidate genes. The candidate genes identified are related to various biological processes and pathways such as heat-shock proteins, oxygen transport, anatomical traits, mitochondrial DNA maintenance, metabolic activity, feed intake, carcass conformation, fertility, and reproduction. This highlights the large number of biological processes involved in thermal tolerance and thus, the polygenic nature of climatic resilience. A comprehensive description of genetic diversity measures in Chinese cattle and YAK was carried out and compared to 24 worldwide cattle breeds to avoid potential biases. Numerous genomic regions under positive selection were detected using three signature of selection methods and candidate genes potentially under positive selection were identified. Enriched function analyses pinpointed important biological pathways, molecular function and cellular components, which contribute to a better understanding of the biological mechanisms underlying thermal tolerance in cattle. Based on the large number of genomic regions identified, thermal tolerance has a complex polygenic inheritance nature, which was expected considering the various mechanisms involved in thermal stress response.


Biology ◽  
2021 ◽  
Vol 10 (9) ◽  
pp. 876
Author(s):  
Alexandra S. Abdelmanova ◽  
Arsen V. Dotsev ◽  
Michael N. Romanov ◽  
Olga I. Stanishevskaya ◽  
Elena A. Gladyr ◽  
...  

Comparison of genomic footprints in chicken breeds with different selection history is a powerful tool in elucidating genomic regions that have been targeted by recent and more ancient selection. In the present work, we aimed at examining and comparing the trajectories of artificial selection in the genomes of the native egg-type Russian White (RW) and meat-type White Cornish (WC) breeds. Combining three different statistics (top 0.1% SNP by FST value at pairwise breed comparison, hapFLK analysis, and identification of ROH island shared by more than 50% of individuals), we detected 45 genomic regions under putative selection including 11 selective sweep regions, which were detected by at least two different methods. Four of such regions were breed-specific for each of RW breed (on GGA1, GGA5, GGA8, and GGA9) and WC breed (on GGA1, GGA5, GGA8, and GGA28), while three remaining regions on GGA2 (two sweeps) and GGA3 were common for both breeds. Most of identified genomic regions overlapped with known QTLs and/or candidate genes including those for body temperatures, egg productivity, and feed intake in RW chickens and those for growth, meat and carcass traits, and feed efficiency in WC chickens. These findings were concordant with the breed origin and history of their artificial selection. We determined a set of 188 prioritized candidate genes retrieved from the 11 overlapped regions of putative selection and reviewed their functions relative to phenotypic traits of interest in the two breeds. One of the RW-specific sweep regions harbored the known domestication gene, TSHR. Gene ontology and functional annotation analysis provided additional insight into a functional coherence of genes in the sweep regions. We also showed a greater candidate gene richness on microchromosomes relative to macrochromosomes in these genomic areas. Our results on the selection history of RW and WC chickens and their key candidate genes under selection serve as a profound information for further conservation of their genomic diversity and efficient breeding.


2021 ◽  
Vol 12 ◽  
Author(s):  
Satish Kumar Illa ◽  
Sabyasachi Mukherjee ◽  
Sapna Nath ◽  
Anupama Mukherjee

BackgroundIn the evolutionary time scale, selection shapes the genetic variation and alters the architecture of genome in the organisms. Selection leaves detectable signatures at the genomic coordinates that provide clues about the protein-coding regions. Sahiwal is a valuable indicine cattle adapted to tropical environments with desirable milk attributes. Insights into the genomic regions under putative selection may reveal the molecular mechanisms affecting the quantitative and other important traits. To understand this, the present investigation was undertaken to explore signatures of selection in the genome of Sahiwal cattle using a medium-density genotyping INDUS chip.ResultDe-correlated composite of multiple selection signals (DCMS), which combines five different univariate statistics, was computed in the dataset to detect the signatures of selection in the Sahiwal genome. Gene annotations, Quantitative Trait Loci (QTL) enrichment, and functional analyses were carried out for the identification of significant genomic regions. A total of 117 genes were identified, which affect a number of important economic traits. The QTL enrichment analysis highlighted 14 significant [False Discovery Rate (FDR)-corrected p-value ≤ 0.05] regions on chromosomes BTA 1, 3, 6, 11, 20, and 21. The top three enriched QTLs were found on BTA 6, 20, and 23, which are associated with exterior, health, milk production, and reproduction traits. The present study on selection signatures revealed some key genes related with coat color (PDGFRA, KIT, and KDR), facial pigmentation (LEF), milk fat percent (MAP3K1, HADH, CYP2U1, and SGMS2), sperm membrane integrity (OSTC), lactation persistency (MRPS30, NNT, CCL28, HMGCS1, NIM1K, ZNF131, and CCDC152), milk yield (GHR and ZNF469), reproduction (NKX2-1 and DENND1A), and bovine tuberculosis susceptibility (RNF144B and PAPSS1). Further analysis of candidate gene prioritization identified four hub genes, viz., KIT, KDR, MAP3K1, and LEF, which play a role in coat color, facial pigmentation, and milk fat percentage in cattle. Gene enrichment analysis revealed significant Gene ontology (GO) terms related to breed-specific coat color and milk fat percent.ConclusionThe key candidate genes and putative genomic regions associated with economic traits were identified in Sahiwal using single nucleotide polymorphism data and the DCMS method. It revealed selection for milk production, coat color, and adaptability to tropical climate. The knowledge about signatures of selection and candidate genes affecting phenotypes have provided a background information that can be further utilized to understand the underlying mechanism involved in these traits in Sahiwal cattle.


Animals ◽  
2019 ◽  
Vol 10 (1) ◽  
pp. 33 ◽  
Author(s):  
Gebremedhin Gebreselassie ◽  
Haile Berihulay ◽  
Lin Jiang ◽  
Yuehui Ma

Sheep (Ovis aries) is one of the most economically, culturally, and socially important domestic animals. They are reared primarily for meat, milk, wool, and fur production. Sheep were reared using natural selection for a long period of time to offer these traits. In fact, this production system has been slowing the productivity and production potential of the sheep. To improve production efficiency and productivity of this animal through genetic improvement technologies, understanding the genetic background of traits such as body growth, weight, carcass quality, fat percent, fertility, milk yield, wool quality, horn type, and coat color is essential. With the development and utilization of animal genotyping technologies and gene identification methods, many functional genes and genetic variants associated with economically important phenotypic traits have been identified and annotated. This is useful and presented an opportunity to increase the pace of animal genetic gain. Quantitative trait loci and genome wide association study have been playing an important role in identifying candidate genes and animal characterization. This review provides comprehensive information on the identified genomic regions and candidate genes associated with production and reproduction traits, and gene function in sheep.


Genes ◽  
2021 ◽  
Vol 12 (3) ◽  
pp. 456
Author(s):  
Hewa Bahithige Pavithra Chathurangi Ariyarathne ◽  
Martin Correa-Luna ◽  
Hugh Thomas Blair ◽  
Dorian John Garrick ◽  
Nicolas Lopez-Villalobos

The objective of this study was to identify genomic regions associated with milk fat percentage (FP), crude protein percentage (CPP), urea concentration (MU) and efficiency of crude protein utilization (ECPU: ratio between crude protein yield in milk and dietary crude protein intake) using grazing, mixed-breed, dairy cows in New Zealand. Phenotypes from 634 Holstein Friesian, Jersey or crossbred cows were obtained from two herds at Massey University. A subset of 490 of these cows was genotyped using Bovine Illumina 50K SNP-chips. Two genome-wise association approaches were used, a single-locus model fitted to data from 490 cows and a single-step Bayes C model fitted to data from all 634 cows. The single-locus analysis was performed with the Efficient Mixed-Model Association eXpedited model as implemented in the SVS package. Single nucleotide polymorphisms (SNPs) with genome-wide association p-values ≤ 1.11 × 10−6 were considered as putative quantitative trait loci (QTL). The Bayes C analysis was performed with the JWAS package and 1-Mb genomic windows containing SNPs that explained > 0.37% of the genetic variance were considered as putative QTL. Candidate genes within 100 kb from the identified SNPs in single-locus GWAS or the 1-Mb windows were identified using gene ontology, as implemented in the Ensembl Genome Browser. The genes detected in association with FP (MGST1, DGAT1, CEBPD, SLC52A2, GPAT4, and ACOX3) and CPP (DGAT1, CSN1S1, GOSR2, HERC6, and IGF1R) were identified as candidates. Gene ontology revealed six novel candidate genes (GMDS, E2F7, SIAH1, SLC24A4, LGMN, and ASS1) significantly associated with MU whose functions were in protein catabolism, urea cycle, ion transportation and N excretion. One novel candidate gene was identified in association with ECPU (MAP3K1) that is involved in post-transcriptional modification of proteins. The findings should be validated using a larger population of New Zealand grazing dairy cows.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Pablo Cáceres ◽  
Agustín Barría ◽  
Kris A. Christensen ◽  
Liane N. Bassini ◽  
Katharina Correa ◽  
...  

AbstractSea lice (Caligus rogercresseyi) is an ectoparasite which causes major production losses in the salmon aquaculture industry worldwide. Atlantic salmon (Salmo salar) and rainbow trout (Oncorhynchus mykiss) are two of the most susceptible salmonid species to sea lice infestation. The objectives of this study were to: (1) identify genomic regions associated with resistance to Caligus rogercresseyi in Atlantic salmon and rainbow trout by performing single-step Genome-Wide Association studies (ssGWAS), and (2) identify candidate genes related to trait variation based on exploring orthologous genes within the associated regions across species. A total of 2626 Atlantic salmon and 2643 rainbow trout were challenged and genotyped with 50 K and 57 K SNP panels, respectively. We ran two independent ssGWAS for sea lice resistance on each species and identified 7 and 13 regions explaining more than 1% of the genetic variance for the trait, with the most important regions explaining 3% and 2.7% for Atlantic salmon and rainbow trout, respectively. We identified genes associated with immune response, cytoskeleton function, and cell migration when focusing on important genomic regions for each species. Moreover, we found 15 common orthogroups which were present in more than one associated genomic region, within- or between-species; however, only one orthogroup showed a clear potential biological relevance in the response against sea lice. For instance, dual-specificity protein phosphatase 10-like (dusp10) and dual-specificity protein phosphatase 8 (dusp8) were found in genomic regions associated with lice density in Atlantic salmon and rainbow trout, respectively. Dusp10 and dusp8 are modulators of the MAPK pathway and might be involved in the differences of the inflammation response between lice resistant and susceptible fish from both species. Our results provide further knowledge on candidate genes related to sea lice resistance and may help establish better control for sea lice in fish populations.


2021 ◽  
Vol 11 (1) ◽  
pp. 59
Author(s):  
Kirsten Voorhies ◽  
Joanne E. Sordillo ◽  
Michael McGeachie ◽  
Elizabeth Ampleford ◽  
Alberta L. Wang ◽  
...  

An unaddressed and important issue is the role age plays in modulating response to short acting β2-agonists in individuals with asthma. The objective of this study was to identify whether age modifies genetic associations of single nucleotide polymorphisms (SNPs) with bronchodilator response (BDR) to β2-agonists. Using three cohorts with a total of 892 subjects, we ran a genome wide interaction study (GWIS) for each cohort to examine SNP by age interactions with BDR. A fixed effect meta-analysis was used to combine the results. In order to determine if previously identified BDR SNPs had an age interaction, we also examined 16 polymorphisms in candidate genes from two published genome wide association studies (GWAS) of BDR. There were no significant SNP by age interactions on BDR using the genome wide significance level of 5 × 10−8. Using a suggestive significance level of 5 × 10−6, three interactions, including one for a SNP within PRAG1 (rs4840337), were significant and replicated at the significance level of 0.05. Considering candidate genes from two previous GWAS of BDR, three SNPs (rs10476900 (near ADRB2) [p-value = 0.009], rs10827492 (CREM) [p-value = 0.02], and rs72646209 (NCOA3) [p-value = 0.02]) had a marginally significant interaction with age on BDR (p < 0.05). Our results suggest age may be an important modifier of genetic associations for BDR in asthma.


2021 ◽  
Vol 22 (7) ◽  
pp. 3477
Author(s):  
Julia Zaborowska ◽  
Bartosz Łabiszak ◽  
Annika Perry ◽  
Stephen Cavers ◽  
Witold Wachowiak

Mountain plants, challenged by vegetation time contractions and dynamic changes in environmental conditions, developed adaptations that help them to balance their growth, reproduction, survival, and regeneration. However, knowledge regarding the genetic basis of species adaptation to higher altitudes remain scarce for most plant species. Here, we attempted to identify such corresponding genomic regions of high evolutionary importance in two closely related European pines, Pinus mugo and P. uncinata, contrasting them with a reference lowland relative—P. sylvestris. We genotyped 438 samples at thousands of single nucleotide polymorphism (SNP) markers, tested their genetic differentiation and population structure followed by outlier detection and gene ontology annotations. Markers clearly differentiated the species and uncovered patterns of population structure in two of them. In P. uncinata three Pyrenean sites were grouped together, while two outlying populations constituted a separate cluster. In P. sylvestris, Spanish population appeared distinct from the remaining four European sites. Between mountain pines and the reference species, 35 candidate genes for altitude-dependent selection were identified, including such encoding proteins responsible for photosynthesis, photorespiration and cell redox homeostasis, regulation of transcription, and mRNA processing. In comparison between two mountain pines, 75 outlier SNPs were found in proteins involved mainly in the gene expression and metabolism.


Sign in / Sign up

Export Citation Format

Share Document