Genome-wide association study uncovers four QTL predisposing to supernumerary teats in cattle

Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints

The question of how local adaptation takes place remains a fundamental question in evolutionary biology. The variation of allele frequencies in genes under selection over environmental gradients remains mainly theoretical and its empirical assessment would help understanding how adaptation happens over environmental clines. To bring new insights to this issue we set up a broad framework which aimed to compare the adaptive trajectories over environmental clines in two domesticated mammal species co-distributed in diversified landscapes. We sequenced the genomes of 160 sheep and 161 goats extensively managed along environmental gradients, including temperature, rainfall, seasonality and altitude, to identify genes and biological processes shaping local adaptation. Allele frequencies at putatively adaptive loci were rarely found to vary gradually along environmental gradients, but rather displayed a discontinuous shift at the extremities of environmental clines. Of the 430 candidate adaptive genes identified, only 6 were orthologous between sheep and goats and those responded differently to environmental pressures, suggesting different putative mechanisms involved in local adaptation in these two closely related species. Interestingly, the genomes of the 2 species were impacted differently by the environment, genes related to signatures of selection were most related to altitude, slope and rainfall seasonality for sheep, and summer temperature and spring rainfall for goats. The diversity of candidate adaptive pathways may result from a high number of biological functions involved in the adaptations to multiple eco-climatic gradients, and a differential role of climatic drivers on the two species, despite their co-distribution along the same environmental gradients. This study describes empirical examples of clinal variation in putatively adaptive alleles with different patterns in allele frequency distributions over continuous environmental gradients, thus showing the diversity of genetic responses in adaptive landscapes and opening new horizons for understanding genomics of adaptation in mammalian species and beyond.

Transcriptome analysis to identify candidate genes related to mammary gland development of Bactrian camel (Camelus bactrianus)

Introduction

The demand for camel milk, which has unique therapeutic properties, is increasing. The mammary gland is the organ in mammals responsible for the production and quality of milk. However, few studies have investigated the genes or pathways related to mammary gland growth and development in Bactrian camels. This study aimed to compare the morphological changes in mammary gland tissue and transcriptome expression profiles between young and adult female Bactrian camels and to explore the potential candidate genes and signaling pathways related to mammary gland development.

Methods

Three 2  years-old female camels and three 5  years-old adult female camels were maintained in the same environment. The parenchyma of the mammary gland tissue was sampled from the camels using percutaneous needle biopsy. Morphological changes were observed using hematoxylin-eosin staining. High-throughput RNA sequencing was performed using the Illumina HiSeq platform to analyze changes in the transcriptome between young and adult camels. Functional enrichment, pathway enrichment, and protein-protein interaction networks were also analyzed. Gene expression was verified using quantitative real-time polymerase chain reaction (qRT-PCR).

Results

Histomorphological analysis showed that the mammary ducts and mammary epithelial cells in adult female camels were greatly developed and differentiated from those in young camels. Transcriptome analysis showed that 2,851 differentially expressed genes were obtained in the adult camel group compared to the young camel group, of which 1,420 were upregulated, 1,431 were downregulated, and 2,419 encoded proteins. Functional enrichment analysis revealed that the upregulated genes were significantly enriched for 24 pathways, including the Hedgehog signaling pathway which is closely related to mammary gland development. The downregulated genes were significantly enriched for seven pathways, among these the Wnt signaling pathway was significantly related to mammary gland development. The protein-protein interaction network sorted the nodes according to the degree of gene interaction and identified nine candidate genes: PRKAB2, PRKAG3, PLCB4, BTRC, GLI1, WIF1, DKK2, FZD3, and WNT4. The expression of fifteen genes randomly detected by qRT-PCR showed results consistent with those of the transcriptome analysis.

Discussion

Preliminary findings indicate that the Hedgehog, Wnt, oxytocin, insulin, and steroid biosynthesis signaling pathways have important effects on mammary gland development in dairy camels. Given the importance of these pathways and the interconnections of the involved genes, the genes in these pathways should be considered potential candidate genes. This study provides a theoretical basis for elucidating the molecular mechanisms associated with mammary gland development and milk production in Bactrian camels.

Design and characterization of a high-resolution multiple-SNP capture array by target sequencing for sheep

The efficiency of molecular breeding largely depends on inexpensive genotyping arrays. In this study, we aimed to develop an ovine high-resolution multiple-single-nucleotide polymorphism (SNP) capture array, based on genotyping by target sequencing (GBTS) system with capture-in-solution (liquid chip) technology. All the markers were from 40K captured regions, including genes located within selective sweep regions, breed-specific regions, quantitative trait loci (QTL), and the potential functional SNPs on the sheep genome. The results showed that a total of 210K high-quality SNPs were identified in the 40K regions, indicating a high average capture ratio (99.7%) for the target genomic regions. Using genotyped data (n = 317) from liquid chip technology, we further performed genome-wide association studies (GWAS) to detect the genetic loci affecting sheep hair types and teat number. A single significant association signal for hair types was identified on 6.7-7.1 Mb of chromosome 25. The IRF2BP2 gene (chr25: 7,067,974-7,071,785), which is located within this genomic region, has been previously known to be involved in hair/wool traits in sheep. The results further showed a new candidate region around 26.4 Mb of chromosome 13, between the ARHGAP21 and KIAA1217 genes, that was significantly related to teat number in sheep. The haplotype patterns of this region also showed differences in animals with 2, 3, or 4 teats. Advances in using the high-accuracy and low-cost liquid chip are expected to accelerate sheep genomic and breeding studies in the coming years.

Variation and frequency of supernumerary teats, litter size, histological features and the fibroblast growth factor 2 (FGF-2) gene expression pattern in goats

Historically, female domestic goats carrying multiple kids are mostly unable to express sufficient nursing ability due to a limited number of functional teats. Therefore, the functional teat is an important component in prolific goat breeding, and plays a key role in the future health of their kids. With this motivation, we wanted to investigate the phenotypic features, litter size, histology of adult female mammary glands, and the gene expression profile of the fibroblast growth factor 2 (FGF-2) gene in goats. To illustrate this, the initial dataset of the current study consists of an electronic questionnaire that includes 697 individuals (548 does and 149 bucks) of five endemic and three exotic goats from 2015 to 2020 in different geographic areas of Iran, from 59 Markhoz (MARG), 50 Azari (AZAR), 73 Busheri (BUSH), 69 Sarbisheh (SARB), 165 Mahabadi (MOHA) indigenous goats and also exotic breeds, including 183 Saanen (SANN), 39 Alpine (ALPN), and 59 Boer (BORE) goats. The results of this study confirmed that MOHA goats (4.16%), BORE (4.43%) and SANN goat breeds (5.75%) have larger litter sizes. Interestingly, the evidence gathering when SNTs occurred showed that both the BUSH and BORE goat breeds had the highest frequency of SNTs. Moreover, under the same physiological and lactation conditions, there was no statistically significant difference in histological features between the three compared does class consist of the two teats, SNTs, and four functional teats. In addition, the results of the gene expression profile significantly highlight the FGF-2 gene pattern in two teat groups compared to other SNT groups (P < 0.01). In summary, this scenario can be used to generate further research and facts on responsible candidate genes, the variations in teat numbers in goats, examining both the incidence of SNT and litter size.

Genetic parameters and weighted single-step genome-wide association study for supernumerary teats in Holstein cattle

Supernumerary teats (SNT) are a common epidermal abnormality of udders in mammals. The SNT negatively affect machine milking ability, udder health, and animal welfare and sometimes act as reservoirs for undesirable bacteria, resulting in economic losses on calves and lactating cows due to the cost of SNT removal surgery, early culling, and low milk yield. This study aimed to analyze the incidence and genetic parameter of SNT and detect SNT-related genes in Chinese Holstein cattle. In this study, the incidence of SNT was recorded in 4,670 Chinese Holstein cattle (born between 2008 and 2017) from 2 farms, including 734 genotyped cows with 114,485 SNPs. The SNT had a total frequency of 9.8% and estimated heritability of 0.22 (SE = 0.07), which were obtained using a threshold model in the studied Chinese Holstein population. Furthermore, we calculated approximate genetic correlations between SNT and the following indicator traits: 12 milk production, 28 body conformation, 5 fertility and reproduction, 5 health, and 9 longevity. Generally, the estimated correlations, such as 305-d milk yield for third parity (-0.55; SE = 0.02) and age at first calving in heifer (0.19; SE = 0.03), were low to moderate. A single-step GWAS was implemented, and 10 genes associated with SNT located in BTA4 were identified. The region (112.70-112.90 Mb) on BTA4 showed the highest genetic variance for SNT. The quantitative trait loci on BTA4 was mapped into the RARRES2 gene, which was previously shown to affect adipogenesis and hormone secretion. The WIF1 gene, which was located in BTA5, was also considered as a candidate gene for SNT. Overall, these findings provide useful information for breeders who are interested in reducing SNT.

Exploiting within-breed variability in the autochthonous Reggiana breed identified several candidate genes affecting pigmentation-related traits, stature and udder defects in cattle

Autochthonous cattle breeds constitute important reservoirs of genetic diversity. Reggiana is an Italian local cattle breed reared in the north of Italy for the production of a mono‐breed Parmigiano–Reggiano cheese. Reggiana cattle usually have a classical solid red coat colour and pale muzzle. As part of the strategies designed for the sustainable conservation of this genetic resource, we investigated at the genome‐wise level the within‐breed detected variability of three pigmentation‐related traits (intensity of red coat colour, based on three classes – light/diluted, normal and dark; spotted patterns/piebaldism that sometime emerge in the breed; muzzle colour – pink/pale, grey and black), stature, presence/absence and number of supernumerary teats and teat length. A total of 1776 Reggiana cattle (about two‐thirds of the extant breed population) were genotyped with the GeneSeek GGP Bovine 150k SNP array and single‐marker and haplotype‐based GWASs were carried out. The results indicated that two main groups of genetic factors affect the intensity of red coat colour: darkening genes (including EDN3 and a few other genes) and diluting genes (including PMEL and a few other genes). Muzzle colour was mainly determined by MC1R gene markers. Piebaldism was mainly associated with KIT gene markers. Stature was associated with BTA6 markers upstream of the NCAPG–LCORL genes. Teat defects were associated with TBX3/TBX5, MCC and LGR5 genes. Overall, the identified genomic regions not only can be directly used in selection plans in the Reggiana breed, but also contribute to clarifying the genetic mechanisms involved in determining exterior traits in cattle.

Genome-wide associations and detection of potential candidate genes for direct genetic and maternal genetic effects influencing dairy cattle body weight at different ages

Background

Body weight (BW) at different ages are of increasing importance in dairy cattle breeding schemes, because of their strong correlation with energy efficiency traits, and their impact on cow health, longevity and farm economy. In total, 15,921 dairy cattle from 56 large-scale test-herds with BW records were genotyped for 45,613 single nucleotide polymorphisms (SNPs). This dataset was used for genome-wide association studies (GWAS), in order to localize potential candidate genes for direct and maternal genetic effects on BW recorded at birth (BW0), at 2 to 3 months of age (BW23), and at 13 to 14 months of age (BW1314).

Results

The first 20 principal components (PC) of the genomic relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G) grouped the genotyped cattle into three clusters. In the statistical models used for GWAS, correction for population structure was done by including polygenic effects with various genetic similarity matrices, such as the pedigree-based relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A), the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G-matrix, the reduced \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G-matrix LOCO (i.e. exclusion of the chromosome on which the candidate SNP is located), and LOCO plus chromosome-wide PC. Inflation factors for direct genetic effects using \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A and LOCO were larger than 1.17. For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G and LOCO plus chromosome-wide PC, inflation factors were very close to 1.0. According to Bonferroni correction, ten, two and seven significant SNPs were detected for the direct genetic effect on BW0, BW23, and BW1314, respectively. Seventy-six candidate genes contributed to direct genetic effects on BW with four involved in growth and developmental processes: FGF6, FGF23, TNNT3, and OMD. For maternal genetic effects on BW0, only three significant SNPs (according to Bonferroni correction), and four potential candidate genes, were identified. The most significant SNP on chromosome 19 explained only 0.14% of the maternal de-regressed proof variance for BW0.

Conclusions

For correction of population structure in GWAS, we suggest a statistical model that considers LOCO plus chromosome-wide PC. Regarding direct genetic effects, several SNPs had a significant effect on BW at different ages, and only two SNPs on chromosome 5 had a significant effect on all three BW traits. Thus, different potential candidate genes regulate BW at different ages. Maternal genetic effects followed an infinitesimal model.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0444-4) contains supplementary material, which is available to authorized users.

A genome-wide association study reveals candidate genes for the supernumerary nipple phenotype in sheep (Ovis aries)

Genome-wide association studies (GWASs) have been widely applied in livestock to identify genes associated with traits of economic interest. Here, we conducted the first GWAS of the supernumerary nipple phenotype in Wadi sheep, a native Chinese sheep breed, based on Ovine Infinium HD SNP BeadChip genotypes in a total of 144 ewes (75 cases with four teats, including two normal and two supernumerary teats, and 69 control cases with two teats). We detected 63 significant SNPs at the chromosome-wise threshold. Additionally, one candidate region (chr1: 170.723-170.734 Mb) was identified by haplotype-based association tests, with one SNP (rs413490006) surrounding functional genes BBX and CD47 on chromosome 1 being commonly identified as significant by the two mentioned analyses. Moreover, Gene Ontology enrichment for the significant SNPs identified by the GWAS analysis was functionally clustered into the categories of receptor activity and synaptic membrane. In addition, pathway mapping revealed four promising pathways (Wnt, oxytocin, MAPK and axon guidance) involved in the development of the supernumerary nipple phenotype. Our results provide novel and important insights into the genetic mechanisms underlying the phenotype of supernumerary nipples in mammals, including humans. These findings may be useful for future breeding and genetics in sheep and other livestock.

Heritability and genome-wide association mapping for supernumerary teats in French Alpine and Saanen dairy goats

This paper reports a quantitative genetics and genomic analysis of undesired presence of supernumerary teats (SNT) in goats. Supernumerary teats are a problem in goat breeding as they can considerably impede machine milking efficiency, leading to increased milking time and injury. This phenotype has routinely been recorded for the past 15 yr in French Alpine and Saanen goats. Around 4% of the females had been assigned the SNT phenotype and consequently could not be included in the breeding program as elite animals. The heritability of this binary trait, estimated by applying linear logistic polygenic models to 32,908 Alpine and 23,217 Saanen females, was 0.40 and 0.44, respectively. A genome-wide association study was implemented using a daughter design composed of 810 Saanen goats sired by 9 artificial insemination bucks and 1,185 Alpine goats sired by 11 bucks, genotyped with the goatSNP50 chip (Illumina Inc., San Diego, CA). This association study was based on logistic polygenic models, one with separately taken single nucleotide polymorphisms and the other with haplotypes as fixed effects. The 2 breeds were analyzed together and separately. No region was found to be significant at the genome level, but 17 regions on 10 chromosomes were significant at the chromosome level. These signals were always only slightly above the chromosome significance threshold and only a few of them overlapped across analyses. No evidence of segregation of a major gene in our Saanen and Alpine populations was observed, suggesting that SNT presence is inherited in a polygenic fashion. This conclusion regarding SNT determinism agrees with recent association analyses in cattle, and one locus was even found in an orthologous region. The possibility of applying markers-based selection on the SNT trait is therefore unlikely, but, as this trait is heritable and routinely recorded, it could be managed by attributing a dedicated estimated breeding value.

A multi-trait meta-analysis with imputed sequence variants reveals twelve QTL for mammary gland morphology in Fleckvieh cattle

Background

The availability of whole-genome sequence data from key ancestors in bovine populations provides an exhaustive catalogue of polymorphic sites that segregate within and across cattle breeds. Sequence variants identified from the sequenced genome of key ancestors can be imputed into animals that have been genotyped using medium- and high-density genotyping arrays. Association analysis with imputed sequences, particularly when applied to multiple traits simultaneously, is a very powerful approach to detect candidate causal variants that underlie complex phenotypes.

Results

We used whole-genome sequence data from 157 key ancestors of the German Fleckvieh cattle population to impute 20,561,798 sequence variants into 10,363 animals that had (partly imputed) genotypes based on 634,109 single nucleotide polymorphisms (SNPs). Rare variants were more frequent among the sequence-derived than the array-derived genotypes. Association studies with imputed sequence variants were performed using seven correlated udder conformation traits as response variables. The calculation of an approximate multi-trait test statistic enabled us to detect 12 quantitative trait loci (QTL) (P < 2.97 × 10⁻⁹) that affect different morphological features of the mammary gland. Among the tested variants, the most significant associations were found for imputed sequence variants at 11 QTL, whereas the top association signal was observed for an array-derived variant at a QTL on bovine chromosome 14. Seven QTL were associated with multiple phenotypes. Most QTL were located in non-coding regions of the genome but in close proximity of candidate genes that could be involved in mammary gland morphology (SP5, GC, NPFFR2, CRIM1, RXFP2, TBX5, RBM19 and ADAM12).

Conclusions

Using imputed sequence variants in association analyses allows the detection of QTL at maximum resolution. Multi-trait approaches can reveal QTL that are not detected in single-trait association studies. Most QTL for udder conformation traits were located in non-coding regions of the genome, which suggests that mutations in regulatory sequences are the major determinants of variation in mammary gland morphology in cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0190-4) contains supplementary material, which is available to authorized users.

Association analysis for udder health based on SNP-panel and sequence data in Danish Holsteins

BACKGROUND

The sensitivity of genome-wide association studies for the detection of quantitative trait loci (QTL) depends on the density of markers examined and the statistical models used. This study compares the performance of three marker densities to refine six previously detected QTL regions for mastitis traits: 54 k markers of a medium-density SNP (single nucleotide polymorphism) chip (MD), imputed 777 k markers of a high-density SNP chip (HD), and imputed whole-genome sequencing data (SEQ). Each dataset contained data for 4496 Danish Holstein cattle. Comparisons were performed using a linear mixed model (LM) and a Bayesian variable selection model (BVS).

RESULTS

After quality control, 587, 7825, and 78 856 SNPs in the six targeted regions remained for MD, HD, and SEQ data, respectively. In general, the association patterns between SNPs and traits were similar for the three marker densities when tested using the same statistical model. With the LM model, 120 (MD), 967 (HD), and 7209 (SEQ) SNPs were significantly associated with mastitis, whereas with the BVS model, 43 (MD), 131 (HD), and 1052 (SEQ) significant SNPs (Bayes factor > 3.2) were observed. A total of 26 (MD), 75 (HD), and 465 (SEQ) significant SNPs were identified by both models. In addition, one, 16, and 33 QTL peaks for MD, HD, and SEQ data were detected according to the QTL intensity profile of SNP bins by post-analysis of the BVS model.

CONCLUSIONS

The power to detect significant associations increased with increasing marker density. The BVS model resulted in clearer boundaries between linked QTL than the LM model. Using SEQ data, the six targeted regions were refined to 33 candidate QTL regions for udder health. The comparison between these candidate QTL regions and known genes suggested that NPFFR2, SLC4A4, DCK, LIFR, and EDN3 may be considered as candidate genes for mastitis susceptibility.

Genome-Wide Specific Selection in Three Domestic Sheep Breeds

BACKGROUND

Commercial sheep raised for mutton grow faster than traditional Chinese sheep breeds. Here, we aimed to evaluate genetic selection among three different types of sheep breed: two well-known commercial mutton breeds and one indigenous Chinese breed.

RESULTS

We first combined locus-specific branch lengths and di statistical methods to detect candidate regions targeted by selection in the three different populations. The results showed that the genetic distances reached at least medium divergence for each pairwise combination. We found these two methods were highly correlated, and identified many growth-related candidate genes undergoing artificial selection. For production traits, APOBR and FTO are associated with body mass index. For meat traits, ALDOA, STK32B and FAM190A are related to marbling. For reproduction traits, CCNB2 and SLC8A3 affect oocyte development. We also found two well-known genes, GHR (which affects meat production and quality) and EDAR (associated with hair thickness) were associated with German mutton merino sheep. Furthermore, four genes (POL, RPL7, MSL1 and SHISA9) were associated with pre-weaning gain in our previous genome-wide association study.

CONCLUSIONS

Our results indicated that combine locus-specific branch lengths and di statistical approaches can reduce the searching ranges for specific selection. And we got many credible candidate genes which not only confirm the results of previous reports, but also provide a suite of novel candidate genes in defined breeds to guide hybridization breeding.

A genome-wide association study in large white and landrace pig populations for number piglets born alive

The number of piglets born alive (NBA) per litter is one of the most important traits in pig breeding due to its influence on production efficiency. It is difficult to improve NBA because the heritability of the trait is low and it is governed by a high number of loci with low to moderate effects. To clarify the biological and genetic background of NBA, genome-wide association studies (GWAS) were performed using 4,012 Large White and Landrace pigs from herdbook and commercial breeding companies in Germany (3), Austria (1) and Switzerland (1). The animals were genotyped with the Illumina PorcineSNP60 BeadChip. Because of population stratifications within and between breeds, clusters were formed using the genetic distances between the populations. Five clusters for each breed were formed and analysed by GWAS approaches. In total, 17 different significant markers affecting NBA were found in regions with known effects on female reproduction. No overlapping significant chromosome areas or QTL between Large White and Landrace breed were detected.