Searching new signals for production traits through gene‐based association analysis in three Italian cattle breeds

Genome‐wide association studies (GWAS) have been widely applied to disentangle the genetic basis of complex traits. In cattle breeds, classical GWAS approaches with medium‐density marker panels are far from conclusive, especially for complex traits. This is due to the intrinsic limitations of GWAS and the assumptions that are made to step from the association signals to the functional variations. Here, we applied a gene‐based strategy to prioritize genotype–phenotype associations found for milk production and quality traits with classical approaches in three Italian dairy cattle breeds with different sample sizes (Italian Brown n = 745; Italian Holstein n = 2058; Italian Simmental n = 477). Although classical regression on single markers revealed only a single genome‐wide significant genotype–phenotype association, for Italian Holstein, the gene‐based approach identified specific genes in each breed that are associated with milk physiology and mammary gland development. As no standard method has yet been established to step from variation to functional units (i.e., genes), the strategy proposed here may contribute to revealing new genes that play significant roles in complex traits, such as those investigated here, amplifying low association signals using a gene‐centric approach.     

Systematic analyses for candidate genes of milk production traits in water buffalo (Bubalus Bubalis)

Water buffalo (Bubalus bubalis) is of great economic importance as a provider of milk and meat in many countries. However, the milk yield of buffalo is much lower than that of Holstein cows. Selection of candidate genes related to milk production traits can be applied to improve buffalo milk performance. A systematic review of studies of these candidate genes will be greatly beneficial for researchers to timely and efficiently understand the research development of molecular markers for buffalo milk production traits. Here, we identified and classified the candidate genes associated with buffalo milk production traits. A total of 517 candidate genes have been identified as being associated with milk performance in different buffalo breeds. Nineteen candidate genes containing 47 mutation sites have been identified using the candidate gene approach. In addition, 499 candidate genes have been identified in six genome‐wide association studies (GWASes) including two studies performed with the bovine SNP chip and four studies with the buffalo SNP chip. Genes CTNND2 (catenin delta 2), APOB (apolipoprotein B), FHIT (fragile histidine triad) and ESRRG (estrogen related receptor gamma) were identified in at least two GWASes. These four genes, especially APOB, deserve further study to explore regulatory roles in buffalo milk production. With growth in the number of buffalo genomic studies, more candidate genes associated with buffalo milk production traits will be identified. Therefore, future studies, such as those investigating gene location and functional analyses, are necessary to facilitate the exploitation of genetic potential and the improvement of buffalo milk performance.     

Genome-wide association study of milk production traits in a crossbred dairy sheep population using three statistical models

Milk production is one of the most important characteristics of dairy sheep, and the identification of genes affecting milk production traits is critical to understanding the genetics and improve milk production in future generations. Three statistical techniques, namely GWAS, ridge-regression BLUP and BayesC , were used to identify SNPs in significant association with three milk production traits (milk yield, fat yield and protein yield) in a crossbred dairy sheep population. The results suggested that chromosomes 1, 3, 4, 5, 7 and 11 were likely to harbor genes important to milk production because these chromosomes had the greatest top-100-SNP variance contributions on the three milk production traits. The GWAS analysis identified between 74 and 288 genome-wide significant SNP (P < 0.05) whereas the BayesCπ model revealed between six and 63 SNPs, each with >95% posterior probability of inclusion as having a non-zero association effect on at least one of the three milk production traits. Positional candidate genes for milk production in sheep were searched, based on the sheep genomic assembly OAR version 3.1, such as those which map position coincided with or was located within 0.1 Mbp of a genome-wide suggestive or significant SNP. These identified SNPs and candidate genes supported some previous findings and also added new information about genetic markers for genetic improvement of lactation in dairy sheep, but keeping in mind that the majority of these positional candidate genes are not necessarily true causative loci for these traits and future validations are thus necessary.     

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.     

Genome‐wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305‐day milk yield in Holstein dairy cows

Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.     

Genetic diversity and genomic signatures of selection among cattle breeds from Siberia,eastern and northern Europe

Domestication in the near eastern region had a major impact on the gene pool of humpless taurine cattle (Bos taurus). As a result of subsequent natural and artificial selection, hundreds of different breeds have evolved, displaying a broad range of phenotypic traits. Here, 10 Eurasian B. taurus breeds from different biogeographic and production conditions, which exhibit different demographic histories and have been under artificial selection at various intensities, were investigated using the Illumina BovineSNP50 panel to understand their genetic diversity and population structure. In addition, we scanned genomes from eight breeds for signatures of diversifying selection. Our population structure analysis indicated six distinct breed groups, the most divergent being the Yakutian cattle from Siberia. Selection signals were shared (experimental P‐value < 0.01) with more than four breeds on chromosomes 6, 7, 13, 16 and 22. The strongest selection signals in the Yakutian cattle were found on chromosomes 7 and 21, where a miRNA gene and genes related to immune system processes are respectively located. In general, genomic regions indicating selection overlapped with known QTL associated with milk production (e.g. on chromosome 19), reproduction (e.g. on chromosome 24) and meat quality (e.g. on chromosome 7). The selection map created in this study shows that native cattle breeds and their genetic resources represent unique material for future breeding.     

Genetic diversity analysis of French goat populations reveals selective sweeps involved in their differentiation

After domestication 11 000 years ago in Asia Minor, the goat followed human migration to Europe and Asia. It was then introduced in Africa and is now raised all over the world. In this study, we exploited a dataset composed of 54 000 SNPs (Illumina goat DNA chip) to analyze the genetic diversity of 223 individuals belonging to eight French breeds (Alpine, Angora, Corse, Fossés, Poitevine, Provençale, Pyrénées and Saanen). Analyses carried out included individual‐based approaches (principal component analysis and population structure) and population‐based approaches (phylogenetic tree constructions). The results of the genetic diversity analyses revealed that French breeds are clearly differentiated, in particular, the Angora breed that originates from south west Asia. The Provençale breed shows a very original genetic pattern that could be the result of ancient admixture. Then, selection signatures were detected by identifying regions of outlying genetic differentiation between populations. Five genomic regions were detected under selection on chromosomes 5, 6, 11, 13 and 20, revealing mainly soft selective sweeps and a few hard selective sweeps and highlighting candidate genes that had been selected for during the evolutionary history of these breeds. Among them, two coat coloration genes (ADAMTS20 and ASIP) and one gene related to milk composition (CSN1S1) were involved.     

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.     

Genome-wide associations for investigating time-dependent genetic effects for milk production traits in dairy cattle

Phenotypic variation in milk production traits has been described over the course of a lactation as well as between different parities. The objective of this study was to investigate whether variation in production is affected by different loci across lactations. A genome-wide association study (GWAS) using a 50-k SNP chip was conducted in 152 divergent German Holstein Friesian cows to test for association with milk production traits over different lactations. The first four lactations were analysed regarding milk yield, fat, protein, lactose, milk urea nitrogen yield and content as well as somatic cell score. Two approaches were used: (i) Wilmink curve parameters were used to assess the genetic effects over the course of a lactation and (ii) test-day yield deviations (YD) were used as a normative approach for a GWAS. The significant effects were largest for markers affecting curve parameters for which there was a statistical power <0.8 of detection even in this small design. While significant markers for YDs were detected in this study, the power to detect effects of a similar magnitude was only 0.11, suggesting that many loci may have been missed with this approach in the present design. Furthermore, all significant effects were specific for a single lactation, leading to the conclusion that the variance explained by a certain locus changes from lactation to lactation. We confirm the common evidence that most production traits vary in the degree of persistency after the peak as a result of genetic influence.     

Detection of selection signatures of population‐specific genomic regions selected during domestication process in Jinhua pigs

Chinese pigs have been undergoing both natural and artificial selection for thousands of years. Jinhua pigs are of great importance, as they can be a valuable model for exploring the genetic mechanisms linked to meat quality and other traits such as disease resistance, reproduction and production. The purpose of this study was to identify distinctive footprints of selection between Jinhua pigs and other breeds utilizing genome‐wide SNP data. Genotyping by genome reducing and sequencing was implemented in order to perform cross‐population extended haplotype homozygosity to reveal strong signatures of selection for those economically important traits. This work was performed at a 2% genome level, which comprised 152 006 SNPs genotyped in a total of 517 individuals. Population‐specific footprints of selective sweeps were searched for in the genome of Jinhua pigs using six native breeds and three European breeds as reference groups. Several candidate genes associated with meat quality, health and reproduction, such as GH1, CRHR2, TRAF4 and CCK, were found to be overlapping with the significantly positive outliers. Additionally, the results revealed that some genomic regions associated with meat quality, immune response and reproduction in Jinhua pigs have evolved directionally under domestication and subsequent selections. The identified genes and biological pathways in Jinhua pigs showed different selection patterns in comparison with the Chinese and European breeds.     

Genetic diversity analyses reveal first insights into breed‐specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH) with pedigree‐based inbreeding coefficients (F_PED). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.     

Candidate gene polymorphisms for behavioural adaptations during urbanization in blackbirds

Successful urban colonization by formerly rural species represents an ideal situation in which to study adaptation to novel environments. We address this issue using candidate genes for behavioural traits that are expected to play a role in such colonization events. We identified and genotyped 16 polymorphisms in candidate genes for circadian rhythms, harm avoidance and migratory and exploratory behaviour in 12 paired urban and rural populations of the blackbird Turdus merula across the Western Palaearctic. An exonic microsatellite in the SERT gene, a candidate gene for harm avoidance behaviour, exhibited a highly significant association with habitat type in an analysis conducted across all populations. Genetic divergence at this locus was consistent in 10 of the 12 population pairs; this contrasts with previously reported stochastic genetic divergence between these populations at random markers. Our results indicate that behavioural traits related to harm avoidance and associated with the SERT polymorphism experience selection pressures during most blackbird urbanization events. These events thus appear to be influenced by homogeneous adaptive processes in addition to previously reported demographic founder events.     

Detecting selection signatures in three Iranian sheep breeds

The objective of genome mapping is to achieve valuable insight into the connection between gene variants (genotype) and observed traits (phenotype). Part of that objective is to understand the selective forces that have operated on a population. Finding links between genotype–phenotype changes makes it possible to identify selective sweeps by patterns of genetic variation and linkage disequilibrium. Based on Illumina 50KSNP chip data, two approaches, XP‐EHH (cross‐population extend haplotype homozygosity) and F_ST (fixation index), were carried out in this research to identify selective sweeps in the genome of three Iranian local sheep breeds: Baluchi (n = 86), Lori‐Bakhtiari (n = 45) and Zel (n = 45). Using both methods, 93 candidate genomic regions were identified as harboring putative selective sweeps. Bioinformatics analysis of the genomic regions showed that signatures of selection related to multiple candidate genes, such as HOXB9, HOXB13, ACAN, NPR2, TRIL, AOX1, CSF2, GHR, TNS2, SPAG8, HINT2, ALS2, AAAS, RARG, SYCP2, CAV1, PPP1R3D, PLA2G7, TTLL7 and C20orf10, that play a role in skeletal system and tail, sugar and energy metabolisms, growth, reproduction, immune and nervous system traits. Our findings indicated diverse genomic selection during the domestication of Iranian sheep breeds.     

A genome‐wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle

A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (P < 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.     

Gelsolin expression in sheep milk somatic cells during lactation

The identification of genes involved in phenotypes related to milk quality is important for both economic and health aspects in livestock production. The aim of this study was to assess the level of gelsolin gene expression in two breeds of dairy sheep – Sarda and Gentile – with pronounced differences in quantitative and qualitative milk traits. Gelsolin, a type of actin-modulating proteins is involved in the processes of actin remodeling during cell growth and apoptosis; therefore a role of this protein in mammary changes during lactation was here hypothesized. Individual milk samples were collected three times during lactation from 26 ewes of the two breeds. The differential gene expression of gelsolin in the two breeds and the three lactation times was estimated by quantitative PCR on RNA extracted from milk somatic cells. Correlations of gelsolin gene expression with milk yield and quality and days of lactation were also estimated. The results showed that gelsolin gene expression was significantly higher in the Sarda compared to the Gentile at each lactation stage, in agreement with the longer lactation duration and the higher daily milk yield of the first breed. Significant correlations of gelsolin gene expression were found with milk fat content in Sarda breed (−0.46, P<0.05). Gelsolin expression analysis confirmed the link between gelsolin gene function and milk fat content of sheep.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.     

Increasing the accuracy of genomic prediction in pure‐bred Limousin beef cattle by including cross‐bred Limousin data and accounting for an F94L variant in MSTN

Explicitly fitting effects for major genes or QTL that account for a large percentage of variation in a whole genomic prediction model may increase prediction accuracy. This study compared approaches to account for a major effect of an F94L variant in the MSTN gene within the genomic prediction using bovine whole‐genomic SNP markers. Among the beef cattle breeds, Limousin have been known to have an F94L variant that is not present in Angus. The reference population in this study consisted of 3060 beef cattle including pure‐bred Limousin (PL), cross‐bred Limousin with Angus (LF) and pure‐bred Angus, genotyped using a BovineSNP50 BeadChip and directly for the MSTN‐F94L variant. We compared prediction accuracies in PL animals using the three datasets from only the PL population, admixed PL and LF (AL) or multibreed analysis using all of the PL, LF and Angus (MB) population according to four‐fold cross‐validation after K‐means clustering. The MSTN‐F94L variant was the most strongly associated with five traits (birth weight, calving ease direct, milk, weaning weight and yield grade) among the 13 measured traits in PL and AL populations. Fitting the MSTN‐F94L variant as a random effect, the genomic prediction accuracies for birth weight increased by 2.7% in PL, by 2.2% in AL and by 3.2% in MB. Prediction accuracies for five traits increased in the MB analysis. Fitting MSTN‐F94L as a fixed effect in PL, AL and MB analyses resulted in increased prediction accuracy in PL for eight traits. Prediction accuracies can be improved by including a causal variant in genomic evaluation compared with simply using whole‐genome SNP markers. Fitting the causal variant as a fixed effect along with markers fitted as random effects resulted in greater prediction accuracies for most traits. Causal variants should be genotyped along with SNP markers.     

Taking advantage from phenotype variability in a local animal genetic resource: identification of genomic regions associated with the hairless phenotype in Casertana pigs

Casertana is an endangered autochthonous pig breed (raised in south‐central Italy) that is considered to be the descendant of the influential Neapolitan pig population that was used to improve British breeds in the 19th century. Casertana pigs are characterized by a typical, almost complete, hairless phenotype, even though a few Casertana pigs are normal haired. In this work, using Illumina PorcineSNP60 BeadChip data, we carried out a genome‐wide association study and an F_ST analysis with this breed by comparing animals showing the classical hairless phenotype (n = 81) versus pigs classified as haired (n = 15). Combining the results obtained with the two approaches, we identified two significant regions: one on porcine chromosome (SSC) 7 and one on SSC15. The SSC7 region contains the forkhead box N3 (FOXN3) gene, the most plausible candidate gene of this region, considering that mutations in another gene of the same family (forkhead box N1; Foxn1 or FOXN1) are responsible for the nude locus in rodents and alopecia in humans. Another potential candidate gene, rho guanine nucleotide exchange factor 10 (ARHGEF10), is located in the SSC15 region. FOXN3 and ARHGEF10 have been detected as differentially expressed in androgenetic and senescent alopecia respectively. This study on an autochthonous pig breed contributes to shed some light on novel genes potentially involved in hair development and growth and demonstrates that local animal breeds can be valuable genetic resources for disclosing genetic factors affecting unique traits, taking advantage of phenotype variability segregating in small populations.