Single nucleotide polymorphisms, haplotypes and combined genotypes of LAP3 gene in bovine and their association with milk production traits

Gene expression analysis of Aquaporin genes in ruminants during growth phase in response to heat stress

Aquaporins (AQPs) are trans-membrane protein involved in water transport and different cellular functions such as cell adhesion, signalling and proliferation. These membrane proteins are essential for key physiological functions such as organ development, osmoregulation, tissue regeneration and metabolism. The regulation of AQP5 gene expression in ruminants during growth phase has not been analysed in-vivo. Therefore, the gene expression pattern was analysed in Jamunapari goats during 3 months to 12 month of age and adult age group in response to heat stress. The genotyping of the AQP5 gene was carried out by High-Resolution Melting (HRM) in four different goat breeds, which indicated four distinct genotypes in the population. The nucleotide diversity for the AQP5 gene ranged from 0.315 and 0.524 across the breeds. Additionally, a close evolutionary relationship between AQP5 and the HSP70 gene was observed, indicating a shared pathway for heat stress regulation. The m-RNA expression level of AQP5 at 3, 9, 12 month and adult age group exhibited 47.24, 1140, 43.17 and 12.55-fold higher expression than control. The m-RNA expression level of the AQP5 gene was up-regulated and significantly higher (P < 0.05) at 9-month age as compared to the other age groups. Heat stress phenotypes were classified based on respiration rate and heart rate, and the m-RNA expression of AQP5 was higher in heat stress-susceptible (HSS) individuals than heat stress-tolerant (HST) individuals at 3, 9, and 12 months of age. The AQP5 plays a significant role in thermoregulation during growth phases in response to heat stress in goats, however, it is required to understand the role of aquaporins at cellular level as well as to establish the association with production performance in ruminant system in-vivo.

Study on Single Nucleotide Polymorphism of LAP3 Gene and Its Correlation with Dairy Quality Traits of Gannan Yak

This study explored the polymorphism of the leucine aminopeptidase (LAP3) gene and its relationship with milk quality characteristics in Gannan yak. A cohort of 162 Gannan yak was genotyped utilizing the Illumina Yak cGPS 7K BeadChip, and the identified single nucleotide polymorphisms (SNPs) were evaluated for their association with milk protein, casein, lactose, and fat concentrations. The results showed that four SNPs (g.4494G > A, g.5919A > G, g.8033G > C, and g.15,615A > G) in the LAP3 gene exhibited polymorphism with information content values of 0.267, 0.267, 0.293, and 0.114, respectively. All four SNPs were in Hardy-Weinberg equilibrium (p > 0.05). The g.4494G > A and g.5919A > G SNPs were significantly associated with protein content (p < 0.05), with homozygous genotypes showing significantly higher protein content than heterozygous genotypes (p < 0.05). The g.8033G > C SNP was significantly associated with casein content, protein content, non-fat solids, and acidity (p < 0.05), with the CC genotype having significantly higher casein, protein, and non-fat solids content than the GG and GC genotypes (p < 0.05). The g.15,615A > G SNP was significantly associated with average fat globule diameter (p < 0.05). In general, the mutations within the LAP3 gene demonstrated a positive impact on milk quality traits in Gannan yak, with mutated genotypes correlating with enhanced milk quality. These results indicate that the LAP3 gene could be a significant or candidate gene affecting milk quality traits in Gannan yak and offer potential genetic markers for molecular breeding programs in this species.

Genome-wide association analysis identify candidate genes for feed efficiency and growth traits in Wenchang chickens

BACKGROUND

Wenchang chickens are one of the most popular local chicken breeds in the Chinese chicken industry. However, the low feed efficiency is the main shortcoming of this breed. Therefore, there is a need to find a more precise breeding method to improve the feed efficiency of Wenchang chickens. In this study, we explored important candidate genes and variants for feed efficiency and growth traits through genome-wide association study (GWAS) analysis.

RESULTS

Estimates of genomic heritability for growth and feed efficiency traits, including residual feed intake (RFI) of 0.05, average daily food intake (ADFI) of 0.21, average daily weight gain (ADG) of 0.24, body weight (BW) at 87, 95, 104, 113 days of age (BW87, BW95, BW104 and BW113) ranged from 0.30 to 0.44. Important candidate genes related to feed efficiency and growth traits were identified, such as PLCE1, LAP3, MED28, QDPR, LDB2 and SEL1L3 genes.

CONCLUSION

The results identified important candidate genes for feed efficiency and growth traits in Wenchang chickens and provide a theoretical basis for the development of new molecular breeding technology.

Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics

BACKGROUND

Rendena is a dual-purpose cattle breed, which is primarily found in the Italian Alps and the eastern areas of the Po valley, and recognized for its longevity, fertility, disease resistance and adaptability to steep Alpine pastures. It is categorized as 'vulnerable to extinction' with only 6057 registered animals in 2022, yet no comprehensive analyses of its molecular diversity have been performed to date. The aim of this study was to analyse the origin, genetic diversity, and genomic signatures of selection in Rendena cattle using data from samples collected in 2000 and 2018, and shed light on the breed's evolution and conservation needs.

RESULTS

Genetic analysis revealed that the Rendena breed shares genetic components with various Alpine and Po valley breeds, with a marked genetic proximity to the Original Braunvieh breed, reflecting historical restocking efforts across the region. The breed shows signatures of selection related to both milk and meat production, environmental adaptation and immune response, the latter being possibly the result of multiple rinderpest epidemics that swept across the Alps in the eighteenth century. An analysis of the Rendena cattle population spanning 18 years showed an increase in the mean level of inbreeding over time, which is confirmed by the mean number of runs of homozygosity per individual, which was larger in the 2018 sample.

CONCLUSIONS

The Rendena breed, while sharing a common origin with Brown Swiss, has developed distinct traits that enable it to thrive in the Alpine environment and make it highly valued by local farmers. Preserving these adaptive features is essential, not only for maintaining genetic diversity and enhancing the ability of this traditional animal husbandry to adapt to changing environments, but also for guaranteeing the resilience and sustainability of both this livestock system and the livelihoods within the Rendena valley.

Genetic variation in bovine LAP3 and SIRT1 genes associated with fertility traits in dairy cattle

BACKGROUND

The genetic progress of fertility and reproduction traits in dairy cattle has been constrained by the low heritability of these traits. Identifying candidate genes and variants associated with fertility and reproduction could enhance the accuracy of genetic selection and expedite breeding process of dairy cattle with low-heritability traits. While the bovine LAP3 and SIRT1 genes exhibit well-documented associations with milk production traits in dairy cattle, their effect on cow fertility have not yet been explored. Eleven single nucleotide polymorphisms (SNPs), comprising five in the promoter (rs717156555: C > G, rs720373055: T > C, rs516876447: A > G, rs461857269: C > T and rs720349928: G > A), two in 5'UTR (rs722359733: C > T and rs462932574: T > G), two in intron 12 (rs110932626: A > G and rs43702363: C > T), and one in 3'UTR of exon 13 (rs41255599: C > T) in LAP3 and one in SIRT1 (rs718329990:T > C) genes, have previously been reported to be associated with various traits of milk production and clinical mastitis in Sahiwal and Karan Fries dairy cattle. In this study, the analysis primarily aimed to assess the impact of SNPs within LAP3 and SIRT1 genes on fertility traits in Sahiwal and Karan Fries cattle. Association studies were conducted using mixed linear models, involving 125 Sahiwal and 138 Karan Fries animals in each breed. The analysis utilized a designated PCR-RFLP panel.

RESULTS

In the promoter region of the LAP3 gene, all variants demonstrated significant (P < 0.05) associations with AFC, except for rs722359733: C > T. However, specific variants with the LAP3 gene's promoter region, namely rs722359733: C > T, rs110932626: A > G, rs43702363: C > T, and rs41255599: C > T, showed significant associations with CI and DO in Sahiwal and Karan Fries cows, respectively. The SNP rs718329990: T > C in the promoter region of SIRT1 gene exhibited a significant association with CI and DO in Sahiwal cattle. Haplotype-based association analysis revealed significant associations between haplotype combinations and AFC, CI and DO in the studied dairy cattle population. Animals with H2H3 and H2H4 haplotype combination exhibited higher AFC, CI and DO than other combinations.

CONCLUSIONS

These results affirm the involvement of the LAP3 and SIRT1 genes in female fertility traits, indicating that polymorphisms within these genes are linked to the studied traits. Overall, the significant SNPs and haplotypes identified in this study could have the potential to enhance herd profitability and ensure long-term sustainability on dairy farms by enabling the selection of animals with early age first calving and enhance reproductive performance in the dairy cattle breeding program.

Identify Candidate Genes Associated with the Weight and Egg Quality Traits in Wenshui Green Shell-Laying Chickens by the Copy Number Variation-Based Genome-Wide Association Study

Copy number variation (CNV), as an essential source of genetic variation, can have an impact on gene expression, genetic diversity, disease susceptibility, and species evolution in animals. To better understand the weight and egg quality traits of chickens, this paper aimed to detect CNVs in Wenshui green shell-laying chickens and conduct a copy number variation regions (CNVRs)-based genome-wide association study (GWAS) to identify variants and candidate genes associated with their weight and egg quality traits to support related breeding efforts. In our paper, we identified 11,035 CNVRs in Wenshui green shell-laying chickens, which collectively spanned a length of 13.1 Mb, representing approximately 1.4% of its autosomal genome. Out of these CNVRs, there were 10,446 loss types, 491 gain types, and 98 mixed types. Notably, two CNVRs showed significant correlations with egg quality, while four CNVRs exhibited significant associations with body weight. These significant CNVRs are located on chromosome 4. Further analysis identified potential candidate genes that influence weight and egg quality traits, including FAM184B, MED28, LAP3, ATOH8, ST3GAL5, LDB2, and SORCS2. In this paper, the CNV map of the Wenshui green shell-laying chicken genome was constructed for the first time through population genotyping. Additionally, CNVRs can be employed as molecular markers to genetically improve chickens' weight and egg quality traits.

The implications of various gene variant combinations on breeding values for Awassi ewe milk production

Background and Aim

Milk production is an important factor to consider in selecting Awassi sheep. This trait is influenced by various genes that can be managed to boost production. The breeding values (BVs) for milk production in Jordan's Awassi sheep flocks have been established. This study determined how combined gene variants of BLG, PRL, CSN3, CSN1S1, and CSN2 affect the estimated BVs (EBVs) for milk production in Awassi sheep.

Materials and Methods

Milk and blood samples were collected from 391 Awassi sheep, which was then subjected to molecular analysis through sequencing in order to identify potential alleles and genotypes that could be linked to the EBVs of milk.

Results

The predicted BVs for milk were significantly influenced by PRL and CSN3 gene variants. Through performing the act of epistasis, the interactions of BLG with CSN3 and -CSN1S1 greatly impacted EBVs for milk production. Likewise, the three-way interaction among PRL, CSN3, and CSN1S1, as well as the combined effect of CSN3 with CSN1S1 and CNS1 significantly improved BVs for milk production. When the breeding selection program incorporates the polymorphisms of these genes, gains in milk production can be obtained.

Conclusion

Alleles within the examined genomic areas are crucial for evaluating BVs and maximizing genetic gain in milk production.

Comparative expression analysis of water buffalo (Bubalus bubalis) to identify genes associated with economically important traits

The milk, meat, skins, and draft power of domestic water buffalo (Bubalus bubalis) provide substantial contributions to the global agricultural economy. The world's water buffalo population is primarily found in Asia, and the buffalo supports more people per capita than any other livestock species. For evaluating the workflow, output rate, and completeness of transcriptome assemblies within and between reference-free (RF) de novo transcriptome and reference-based (RB) datasets, abundant bioinformatics studies have been carried out to date. However, comprehensive documentation of the degree of consistency and variability of the data produced by comparing gene expression levels using these two separate techniques is lacking. In the present study, we assessed the variations in the number of differentially expressed genes (DEGs) attained with RF and RB approaches. In light of this, we conducted a study to identify, annotate, and analyze the genes associated with four economically important traits of buffalo, viz., milk volume, age at first calving, post-partum cyclicity, and feed conversion efficiency. A total of 14,201 and 279 DEGs were identified in RF and RB assemblies. Gene ontology (GO) terms associated with the identified genes were allocated to traits under study. Identified genes improve the knowledge of the underlying mechanism of trait expression in water buffalo which may support improved breeding plans for higher productivity. The empirical findings of this study using RNA-seq data-based assembly may improve the understanding of genetic diversity in relation to buffalo productivity and provide important contributions to answer biological issues regarding the transcriptome of non-model organisms.

Incorporating genome-wide and transcriptome-wide association studies to identify genetic elements of longissimus dorsi muscle in Huaxi cattle

Locating the genetic variation of important livestock and poultry economic traits is essential for genetic improvement in breeding programs. Identifying the candidate genes for the productive ability of Huaxi cattle was one crucial element for practical breeding. Based on the genotype and phenotype data of 1,478 individuals and the RNA-seq data of 120 individuals contained in 1,478 individuals, we implemented genome-wide association studies (GWAS), transcriptome-wide association studies (TWAS), and Fisher's combined test (FCT) to identify the candidate genes for the carcass trait, the weight of longissimus dorsi muscle (LDM). The results indicated that GWAS, TWAS, and FCT identified seven candidate genes for LDM altogether: PENK was located by GWAS and FCT, PPAT was located by TWAS and FCT, and XKR4, MTMR3, FGFRL1, DHRS4, and LAP3 were only located by one of the methods. After functional analysis of these candidate genes and referring to the reported studies, we found that they were mainly functional in the progress of the development of the body and the growth of muscle cells. Combining advanced breeding techniques such as gene editing with our study will significantly accelerate the genetic improvement for the future breeding of Huaxi cattle.

Genetic variation in promoter region of the bovine LAP3 gene associated with estimated breeding values of milk production traits and clinical mastitis in dairy cattle

The purpose of this study was to identify genetic variants in the promoter and 5'UTR regions of bovine leucine amino peptidase three (LAP3) gene and analysed their associations with estimated breeding values (EBVs) of milk production traits and clinical mastitis in Sahiwal and Karan Fries cattle. Eleven SNPs were identified within the region under study of the LAP3 gene, including seven promoter variants (rs717156555: C>G, rs720373055: T>C, rs715189731: A>G, rs516876447: A>G, rs461857269: C>T, rs136548163: C>T, and rs720349928: G>A) and four 5'UTR variants (rs717884982: C>T, rs722359733: C>T, rs481631804: C>T and rs462932574: T>G). Out of them, 10 SNPs variants were found in both Sahiwal and Karan Fries cattle, with one SNP variant (rs481631804: C>T) being unique to Karan Fries cattle. Seven of these identified SNPs were chosen for association analyses. Individual SNP based association analysis revealed that two SNPs (rs720373055: T>C and rs720349928: G>A) were significantly associated with EBVs of lactation milk yield (LMY), 305-day milk yield (305dMY), and one significant association of SNP rs722359733: C>T with lactation length (LL) was observed. Haplotype based association analysis indicated that diplotypes are significantly associated with EBVs of LMY, 305dMY, and LL, individuals with H1H3 (CTACGCT/GCGTACG) being linked to higher lactation performance than other diplotypes. Further logistic regression analysis revealed that, animals with diplotype H1H3 was less susceptible to the incidence of clinical mastitis than other cows, as the odds ratio for the non-incidence of clinical mastitis was found to be low. Altogether, variations in the LAP3 gene promoter could be used as a genetic marker, most notably diplotype H1H3, may greatly benefit the simultaneous improvement of mastitis resistance and milk yield traits in dairy cattle. Moreover, bioinformatics analysis predicted that the SNPs rs720373055: T>C, rs715189731:A>G and rs720349928: G>A is located in the core promoter region and in TFBs, play key role in regulation of studied phenotypes.

Assessing genomic diversity and signatures of selection in Pinan cattle using whole-genome sequencing data

BACKGROUND

Crossbreeding is an important way to improve production beef cattle performance. Pinan cattle is a new hybrid cattle obtained from crossing Piedmontese bulls with Nanyang cows. After more than 30 years of cross-breeding, Pinan cattle show a variety of excellent characteristics, including fast growth, early onset of puberty, and good meat quality. In this study, we analyzed the genetic diversity, population structure, and genomic region under the selection of Pinan cattle based on whole-genome sequencing data of 30 Pinan cattle and 169 published cattle genomic data worldwide.  RESULTS: Estimating ancestry composition analysis showed that the composition proportions for our Pinan cattle were mainly Piedmontese and a small amount of Nanyang cattle. The analyses of nucleotide diversity and linkage disequilibrium decay indicated that the genomic diversity of Pinan cattle was higher than that of European cattle and lower than that of Chinese indigenous cattle. De-correlated composite of multiple selection signals, which combines four different statistics including θπ, CLR, F_ST, and XP-EHH, was computed to detect the signatures of selection in the Pinan cattle genome. A total of 83 genes were identified, affecting many economically important traits. Functional annotation revealed that these selected genes were related to immune (BOLA-DQA2, BOLA-DQB, LSM14A, SEC13, and NAALADL2), growth traits (CYP4A11, RPL26, and MYH10), embryo development (REV3L, NT5E, CDX2, KDM6B, and ADAMTS9), hornless traits (C1H21orf62), and climate adaptation (ANTXR2).

CONCLUSION

In this paper, we elucidated the genomic characteristics, ancestry composition, and selective signals related to important economic traits in Pinan cattle. These results will provide the basis for further genetic improvement of Pinan cattle and reference for other hybrid cattle related studies.

Weighted Single-Step Genome-Wide Association Study Uncovers Known and Novel Candidate Genomic Regions for Milk Production Traits and Somatic Cell Score in Valle del Belice Dairy Sheep

Simple Summary

Milk production is the most economically crucial dairy sheep trait and constitutes the major genetic enhancement purpose via selective breeding. Also, mastitis is one of the most frequently encountered diseases, having a significant impact on animal welfare, milk yield, and quality. The aim of this study was to identify genomic region(s) associated with the milk production traits and somatic cell score (SCS) in Valle del Belice sheep using single-step genome-wide association (ssGWA) and genotyping data from medium density SNP panels. We identified several genomic regions (OAR1, OAR2, OAR3, OAR4, OAR6, OAR9, and OAR25) and candidate genes implicated in milk production traits and SCS. Our findings offer new insights into the genetic basis of milk production traits and SCS in dairy sheep.

Abstract

The objective of this study was to uncover genomic regions explaining a substantial proportion of the genetic variance in milk production traits and somatic cell score in a Valle del Belice dairy sheep. Weighted single-step genome-wide association studies (WssGWAS) were conducted for milk yield (MY), fat yield (FY), fat percentage (FAT%), protein yield (PY), protein percentage (PROT%), and somatic cell score (SCS). In addition, our aim was also to identify candidate genes within genomic regions that explained the highest proportions of genetic variance. Overall, the full pedigree consists of 5534 animals, of which 1813 ewes had milk data (15,008 records), and 481 ewes were genotyped with a 50 K single nucleotide polymorphism (SNP) array. The effects of markers and the genomic estimated breeding values (GEBV) of the animals were obtained by five iterations of WssGBLUP. We considered the top 10 genomic regions in terms of their explained genomic variants as candidate window regions for each trait. The results showed that top ranked genomic windows (1 Mb windows) explained 3.49, 4.04, 5.37, 4.09, 3.80, and 5.24% of the genetic variances for MY, FY, FAT%, PY, PROT%, and total SCS, respectively. Among the candidate genes found, some known associations were confirmed, while several novel candidate genes were also revealed, including PPARGC1A, LYPLA1, LEP, and MYH9 for MY; CACNA1C, PTPN1, ROBO2, CHRM3, and ERCC6 for FY and FAT%; PCSK5 and ANGPT1 for PY and PROT%; and IL26, IFNG, PEX26, NEGR1, LAP3, and MED28 for SCS. These findings increase our understanding of the genetic architecture of six examined traits and provide guidance for subsequent genetic improvement through genome selection.

Single and combined effects of CSN1S1 and CSN2-casein genes on Awassi sheep milk quantity and quality

Background and Aim: Milk produced from Awassi sheep is of high nutritive value; its production is relatively low in Awassi sheep, so the genetic improvement programs targeted milk production and its components are of high importance, especially when using genes that have an important signal to milk traits. This study was aimed at assessing the influence of alpha S1 (CSN1S1) and beta-casein (CSN2) genes genotypes interaction on Awassi ewes milk productivity. Materials and Methods: A total number of 391 milk yield and its composition records (taken through five consecutive years, 2007-2011) of 167 ewes were utilized for this study. DNA samples were extracted from the ewe's blood samples, then the polymerase chain reaction products of alpha S1 (CSN1S1) and beta-casein (CSN2) genes were sequenced. The obtained sequences were analyzed; thereafter, the detected variants were tested for their possible association with milk traits. Results: The CSN1S1 and CSN2 variants allelic frequencies were 0.85 and 0.15, and 0.95 and 0.05, respectively. Lactose and solid not fat (SNF) % were associated with TC CSN1S1 genotypes. No association was found among CSN1S1 polymorphic genotypes with milk production, lactose, and SNF % were associated with TC CSN1S1 genotypes. Ewes of CSN2 AC genotype showed higher milk production traits, while no association was found between milk composition traits and CNS2 genotypes. Nevertheless, CSN1S1∗CSN2 interaction showed the highest SNF, fat percentages, and milk production. Conclusion: The substantial interaction effects between CSN1S1×CSN2 genes were significantly affected the amount of milk, fat, and SNF% produced. The detected variants should be included in the breeding programs of Awassi sheep that are designed for improving their milk quantity and quality.

Genome-wide association study on milk production and somatic cell score for Thai dairy cattle using weighted single-step approach with random regression test-day model

Genome-wide association studies are a powerful tool to identify genomic regions and variants associated with phenotypes. However, only limited mutual confirmation from different studies is available. The objectives of this study were to identify genomic regions as well as genes and pathways associated with the first-lactation milk, fat, protein, and total solid yields; fat, protein, and total solid percentage; and somatic cell score (SCS) in a Thai dairy cattle population. Effects of SNPs were estimated by a weighted single-step GWAS, which back-solved the genomic breeding values predicted using single-step genomic BLUP (ssGBLUP) fitting a single-trait random regression test-day model. Genomic regions that explained at least 0.5% of the total genetic variance were selected for further analyses of candidate genes. Despite the small number of genotyped animals, genomic predictions led to an improvement in the accuracy over the traditional BLUP. Genomic predictions using weighted ssGBLUP were slightly better than the ssGBLUP. The genomic regions associated with milk production traits contained 210 candidate genes on 19 chromosomes [Bos taurus autosome (BTA) 1 to 7, 9, 11 to 16, 20 to 21, 26 to 27 and 29], whereas 21 candidate genes on 3 chromosomes (BTA 11, 16, and 21) were associated with SCS. Many genomic regions explained a small fraction of the genetic variance, indicating polygenic inheritance of the studied traits. Several candidate genes coincided with previous reports for milk production traits in Holstein cattle, especially a large region of genes on BTA14. We identified 141 and 5 novel genes related to milk production and SCS, respectively. These novel genes were also found to be functionally related to heat tolerance (e.g., SLC45A2, IRAG1, and LOC101902172), longevity (e.g., SYT10 and LOC101903327), and fertility (e.g., PAG1). These findings may be attributed to indirect selection in our population. Identified biological networks including intracellular cell transportation and protein catabolism implicate milk production, whereas the immunological pathways such as lymphocyte activation are closely related to SCS. Further studies are required to validate our findings before exploiting them in genomic selection.

Genome-Wide Patterns of Homozygosity Reveal the Conservation Status in Five Italian Goat Populations

Simple Summary

In the local populations, the increase in inbreeding is a relevant problem for the reduction in production, reproduction, and adaptive traits. The application of genomic technologies has facilitated the assessment of inbreeding in these populations. The current study aims to investigate the patterns of homozygosity in five Italian local goat populations. The results showed the different selection histories and breeding schemes of these goat populations. The analysis also indicated the importance of this information to avoid future loss of diversity and to produce information for designing optimal breeding and conservation programs.

Abstract

The application of genomic technologies has facilitated the assessment of genomic inbreeding based on single nucleotide polymorphisms (SNPs). In this study, we computed several runs of homozygosity (ROH) parameters to investigate the patterns of homozygosity using Illumina Goat SNP50 in five Italian local populations: Argentata dell’Etna (N = 48), Derivata di Siria (N = 32), Girgentana (N = 59), Maltese (N = 16) and Messinese (N = 22). The ROH results showed well-defined differences among the populations. A total of 3687 ROH segments >2 Mb were detected in the whole sample. The Argentata dell’Etna and Messinese were the populations with the lowest mean number of ROH and inbreeding coefficient values, which reflect admixture and gene flow. In the Girgentana, we identified an ROH pattern related with recent inbreeding that can endanger the viability of the breed due to reduced population size. The genomes of Derivata di Siria and Maltese breeds showed the presence of long ROH (>16 Mb) that could seriously impact the overall biological fitness of these breeds. Moreover, the results confirmed that ROH parameters are in agreement with the known demography of these populations and highlighted the different selection histories and breeding schemes of these goat populations. In the analysis of ROH islands, we detected harbored genes involved with important traits, such as for milk yield, reproduction, and immune response, and are consistent with the phenotypic traits of the studied goat populations. Finally, the results of this study can be used for implementing conservation programs for these local populations in order to avoid further loss of genetic diversity and to preserve the production and fitness traits. In view of this, the availability of genomic data is a fundamental resource.

Polymorphism in DNA repair gene BRCA1 associated with clinical mastitis and production traits in indigenous dairy cattle

The present study was conducted on Karan Fries cows to identify SNPs in breast cancer 1 (BRCA1) gene and their association with inflammation of mammary gland and production traits. Chi square analysis showed significant association of both SNPs i.e. G43737229T (P < 0.01) and G43761121A (P < 0.05) with incidence of clinical mastitis. Logistic modeling showed that GG genotype cows of G43737229T SNP was less susceptible to mastitis in comparison to TT genotype cows as the odds ratio of TT against GG for the non-incidence of mastitis was less i.e. 0.23. Whereas, cows with genotype AA of G43761121A SNP were less prone to clinical mastitis in comparison to GG genotype cows. G-T haplotype and GGTT combined genotype exhibited maximum mastitis incidences. Least squares analysis revealed significant association of G43737229T SNP with total milk yield (P < 0.01), 305 days milk yield (P < 0.01) and lactation length (P < 0.05). The TT and GT genotyped Karan Fries cows were associated with higher milk yields during first and second lactation, as well as with first lactation length. These observations on SNPs, genotypes and combined genotypes of BRCA1 gene and their differential association with the incidence of mastitis and production traits can be utilized as an aid to selection for simultaneous improvement of both antagonistic traits.

Identification and characterization of trait-specific SNPs using ddRAD sequencing in water buffalo

Single Nucleotide Polymorphism (SNP) is one of the important molecular markers widely used in animal breeding program for improvement of any desirable genetic traits. Considering this, the present study was carried out to identify, annotate and analyze the SNPs related to four important traits of buffalo viz. milk volume, age at first calving, post-partum cyclicity and feed conversion efficiency. We identified 246,495, 168,202, 74,136 and 194,747 genome-wide SNPs related to mentioned traits, respectively using ddRAD sequencing technique based on 85 samples of Murrah Buffaloes. Distribution of these SNPs were highest (61.69%) and lowest (1.78%) in intron and exon regions, respectively. Under coding regions, the SNPs for the four traits were further classified as synonymous (4697) and non-synonymous (3827). Moreover, Gene Ontology (GO) terms of identified genes assigned to various traits. These characterized SNPs will enhance the knowledge of cellular mechanism for enhancing productivity of water buffalo through molecular breeding.

Polymorphisms of LAP3 gene and their association with the growth traits in the razor clam Sinonovacula constricta

Leucine aminopeptidase 3 (LAP3) is an important proteolytic enzyme that catalyzes the hydrolysis of leucine residues from the amino termini of protein or peptide substrates and plays a critical role in protein metabolism and growth. In the present study, we investigated the full-length cDNA sequence of the LAP3 gene in Sinonovacula constricta (ScLAP3) using expressed sequence tags and rapid amplification of cDNA ends. The full-length ScLAP3 cDNA was 2885 bp, with a 1560 bp open reading frame encoding 519 amino acids. Sequence analysis revealed that ScLAP3 shared 70.9% identity with LAP3 from the blood clam Tegillarca granosa and 62.0-68.0% with other species. ScLAP3 was expressed in all six tested tissues, with significantly higher expression levels in the foot compared with mantle, adductor muscle, liver, gills, and siphon tissues in adults (P < 0.01). In the eight developmental stages, ScLAP3 expression gradually increased, with significantly higher levels in D-shaped larvae compared with other developmental stages (P < 0.01), suggesting that it may be involved in the formation of certain organs during early development. Association analysis identified three shared single nucleotide polymorphisms (SNPs), c.1073A > G, c.1139C > T and c.1154A > G in exons of ScLAP3 gene from 177 individuals of two groups, one selective strain and one wild population, which had significant effects on growth traits of S. constricta. The results provided candidate genetic markers to assist selective breeding of razor clams toward improved growth.

Genome-wide association and pathway analysis of carcass and meat quality traits in Piemontese young bulls

A key concern in beef production is how to improve carcass and meat quality traits. Identifying the genomic regions and biological pathways that contribute to explaining variability in these traits is of great importance for selection purposes. In this study, genome wide-association (GWAS) and pathway-based analyses of carcass traits (age at slaughter (AS), carcass weight (CW), carcass daily gain (CDG), conformation score and rib-eye muscle area) and meat quality traits (pH, Warner-Bratzler shear force, purge loss, cooking loss and colour parameters (lightness, redness, yellowness, chroma, hue)) were conducted using genotype data from the 'GeneSeek Genomic Profiler Bovine LD' array in a cohort of 1166 double-muscled Piemontese beef cattle. The genome wide-association analysis was based on the GRAMMAR-GC approach and identified 37 significant single nucleotide polymorphisms (SNPs), which were associated with 12 traits (P<5 × 10-5). In particular, 14 SNPs associated with CW, CDG and AS were located at 38.57 to 38.94 Mb on Bos taurus autosome 6 and mapped within four genes, that is, Leucine Aminopeptidase 3, Family with Sequence Similarity 184 Member B, Non-SMC Condensin I Complex Subunit G and Ligand-Dependent Nuclear Receptor Corepressor-Like. Strong pairwise linkage disequilibrium was found in this region. For meat quality traits, most associations were 1 SNP per trait, except for a signal on BTA25 (at ~11.96 Mb), which was significant for four of the five meat colour parameters assessed. Gene-set enrichment analyses yielded significant results for six traits (right-sided hypergeometric test, false discovery rate <0.05). In particular, several pathways related to transmembrane transport (i.e., oxygen, calcium, ion and cation) were overrepresented for meat colour parameters. The results obtained provide useful information for genomic selection for beef production and quality in the Piemontese breed.

Effect and Interaction of β-Lactoglobulin, Kappa Casein, and Prolactin Genes on Milk Production and Composition of Awassi Sheep

Simple Summary

Beta-lactoglobulin (β-LG), prolactin (PRL), and Kappa casein (CSN3) all contribute to the determination of milk production and composition, but have not been assessed in local Awassi sheep. Therefore, our aim was to analyze the contribution of these genes in milk production and composition traits in commercial Awassi ewe population by genotyping and sequencing these genes. Our results showed the prevalence of the different variations (alleles) of the tested genes in the Awassi population, and no association among β-LG and CSN3 polymorphic genotypes and milk production, or PRL and fat%. Also, all 3 genes help determine the milk production potential of Awassi ewes and help assess milk components, and thus can be used in breeding programs to select for milk potential.

Abstract

A participatory animal-breeding program was applied to 9 commercial Awassi sheep flocks in Jordan. This study aimed to assess the influence of Beta-lactoglobulin (β-LG), Prolactin (PRL), and Kappa casein (CSN3) genes, genotypes and their interaction on milk production and composition traits of 167 genotyped Awassi ewes via Polymerase Chain Reaction (PCR) followed by sequencing. Allele frequencies for the two variants were 0.42 and 0.58 for β-LG, 0.82 and 0.18 for PRL, and 0.92 and 0.08 for CSN3. No association was found among β-LG and CSN3 polymorphic genotypes with milk production traits. However, ewes with PRL AA genotype showed higher milk production, β-LG AB was associated with lowest fat%, high solid not fat (SNF)%, protein%, and lactose%. β-LG BB was associated with highest milk density. PRL, β-LG, and CSN3 polymorphic genotypes were differentially associated with milk production and component traits. Furthermore, β-LG × PRL interaction showed the highest milk production and fat%; β-LG × PRL recorded the highest SNF%, protein%, lactose%, and milk density, while the PRL × CSN3 had the highest fat% and SNF%. The enhancing effects of these gene interactions can be incorporated in Awassi breeding programs to improve milk production and composition.

Genome-Wide Association Study for Milk Protein Composition Traits in a Chinese Holstein Population Using a Single-Step Approach

Genome-wide association studies (GWASs) have been widely used to determine the genetic architecture of quantitative traits in dairy cattle. In this study, with the aim of identifying candidate genes that affect milk protein composition traits, we conducted a GWAS for nine such traits (α_s1-casein, α_s2-casein, β-casein, κ-casein, α-lactalbumin, β-lactoglobulin, casein index, protein percentage, and protein yield) in 614 Chinese Holstein cows using a single-step strategy. We used the Illumina BovineSNP50 Bead chip and imputed genotypes from high-density single-nucleotide polymorphisms (SNPs) ranging from 50 to 777 K, and subsequent to genotype imputation and quality control, we screened a total of 586,304 informative high-quality SNPs. Phenotypic observations for six major milk proteins (α_s1-casein, α_s2-casein, β-casein, κ-casein, α-lactalbumin, and β-lactoglobulin) were evaluated as weight proportions of the total protein fraction (wt/wt%) using a commercial enzyme-linked immunosorbent assay kit. Informative windows comprising five adjacent SNPs explaining no < 0.5% of the genomic variance per window were selected for gene annotation and gene network and pathway analyses. Gene network analysis performed using the STRING Genomics 10.0 database revealed a co-expression network comprising 46 interactions among 62 of the most plausible candidate genes. A total of 178 genomic windows and 194 SNPs on 24 bovine autosomes were significantly associated with milk protein composition or protein percentage. Regions affecting milk protein composition traits were mainly observed on chromosomes BTA 1, 6, 11, 13, 14, and 18. Of these, several windows were close to or within the CSN1S1, CSN1S2, CSN2, CSN3, LAP3, DGAT1, RPL8, and HSF1 genes, which have well-known effects on milk protein composition traits of dairy cattle. Taken together with previously reported quantitative trait loci and the biological functions of the identified genes, we propose 19 novel candidate genes affecting milk protein composition traits: ARL6, SST, EHHADH, PCDHB4, PCDHB6, PCDHB7, PCDHB16, SLC36A2, GALNT14, FPGS, LARP4B, IDI1, COG4, FUK, WDR62, CLIP3, SLC25A21, IL5RA, and ACADSB. Our findings provide important insights into milk protein synthesis and indicate potential targets for improving milk quality.

Scans for signatures of selection in Russian cattle breed genomes reveal new candidate genes for environmental adaptation and acclimation

Domestication and selective breeding has resulted in over 1000 extant cattle breeds. Many of these breeds do not excel in important traits but are adapted to local environments. These adaptations are a valuable source of genetic material for efforts to improve commercial breeds. As a step toward this goal we identified candidate regions to be under selection in genomes of nine Russian native cattle breeds adapted to survive in harsh climates. After comparing our data to other breeds of European and Asian origins we found known and novel candidate genes that could potentially be related to domestication, economically important traits and environmental adaptations in cattle. The Russian cattle breed genomes contained regions under putative selection with genes that may be related to adaptations to harsh environments (e.g., AQP5, RAD50, and RETREG1). We found genomic signatures of selective sweeps near key genes related to economically important traits, such as the milk production (e.g., DGAT1, ABCG2), growth (e.g., XKR4), and reproduction (e.g., CSF2). Our data point to candidate genes which should be included in future studies attempting to identify genes to improve the extant breeds and facilitate generation of commercial breeds that fit better into the environments of Russia and other countries with similar climates.

Human-Mediated Introgression of Haplotypes in a Modern Dairy Cattle Breed

Domestic animals can serve as model systems of adaptive introgression and their genomic signatures. In part, their usefulness as model systems is due to their well-known histories. Different breeding strategies such as introgression and artificial selection have generated numerous desirable phenotypes and superior performance in domestic animals. The modern Danish Red Dairy Cattle is studied as an example of an introgressed population. It originates from crossing the traditional Danish Red Dairy Cattle with the Holstein and Brown Swiss breeds, both known for high milk production. This crossing happened, among other things due to changes in the production system, to raise milk production and overall performance. The genomes of modern Danish Red Dairy Cattle are heavily influenced by regions introgressed from the Holstein and Brown Swiss breeds and under subsequent selection in the admixed population. The introgressed proportion of the genome was found to be highly variable across the genome. Haplotypes introgressed from Holstein and Brown Swiss contained or overlapped known genes affecting milk production, as well as protein and fat content (CD14, ZNF215, BCL2L12, and THRSP for Holstein origin and ITPR2, BCAT1, LAP3, and MED28 for Brown Swiss origin). Genomic regions with high introgression signals also contained genes and enriched QTL associated with calving traits, body confirmation, feed efficiency, carcass, and fertility traits. These introgressed signals with relative identity-by-descent scores larger than the median showing Holstein or Brown Swiss introgression are mostly significantly correlated with the corresponding test statistics from signatures of selection analyses in modern Danish Red Dairy Cattle. Meanwhile, the putative significant introgressed signals have a significant dependency with the putative significant signals from signatures of selection analyses. Artificial selection has played an important role in the genomic footprints of introgression in the genome of modern Danish Red Dairy Cattle. Our study on a modern cattle breed contributes to an understanding of genomic consequences of selective introgression by demonstrating the extent to which adaptive effects contribute to shape the specific genomic consequences of introgression.

Genomic selection signatures in sheep from the Western Pyrenees

BACKGROUND

The current large spectrum of sheep phenotypic diversity results from the combined product of sheep selection for different production traits such as wool, milk and meat, and its natural adaptation to new environments. In this study, we scanned the genome of 25 Sasi Ardi and 75 Latxa sheep from the Western Pyrenees for three types of regions under selection: (1) regions underlying local adaptation of Sasi Ardi semi-feral sheep, (2) regions related to a long traditional dairy selection pressure in Latxa sheep, and (3) regions experiencing the specific effect of the modern genetic improvement program established for the Latxa breed during the last three decades.

RESULTS

Thirty-two selected candidate regions including 147 annotated genes were detected by using three statistical parameters: pooled heterozygosity H, Tajima's D, and Wright's fixation index F_st. For Sasi Ardi sheep, chromosomes Ovis aries (OAR)4, 6, and 22 showed the strongest signals and harbored several candidate genes related to energy metabolism and morphology (BBS9, ELOVL3 and LDB1), immunity (NFKB2), and reproduction (H2AFZ). The major genomic difference between Sasi Ardi and Latxa sheep was on OAR6, which is known to affect milk production, with highly selected regions around the ABCG2, SPP1, LAP3, NCAPG, LCORL, and MEPE genes in Latxa sheep. The effect of the modern genetic improvement program on Latxa sheep was also evident on OAR15, on which several olfactory genes are located. We also detected several genes involved in reproduction such as ESR1 and ZNF366 that were affected by this selection program.

CONCLUSIONS

Natural and artificial selection have shaped the genome of both Sasi Ardi and Latxa sheep. Our results suggest that Sasi Ardi traits related to energy metabolism, morphological, reproductive, and immunological features have been under positive selection to adapt this semi-feral sheep to its particular environment. The highly selected Latxa sheep for dairy production showed clear signatures of selection in genomic regions related to milk production. Furthermore, our data indicate that the selection criteria applied in the modern genetic improvement program affect immunity and reproduction traits.

Effect of single nucleotide polymorphisms on intramuscular fat content in Hungarian Simmental cattle

Objective

To estimate effect of single nucleotide polymorphisms on the intramuscular fat content (IMF) of Hungarian Simmental bulls.

Methods

Genotypes were determined on high-density Illumina Bovine DNA Chip. After slaughtering of animals, chemical percentage of intramuscular fat was determined from longissimus dorsi muscle. A multi-locus mixed-model was applied for statistical analyses.

Results

Analyses revealed four loci (rs43284251, rs109210955, rs41630030, and rs41642251) to be highly associated (−log₁₀P>12) with IMF located on chromosome 1, 6, 13, and 17, respectively. The frequency of their minor alleles was 0.426, 0.221, 0.162, and 0.106.

Conclusion

The loci above can be useful in selection programs and gives the possibility to assist selection by molecular tools.

Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium-density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (F_ROH  = 0.099), whereas Comisana showed the lowest (F_ROH  = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.

Effect of DGAT1, FASN and PRL genes on milk production and milk composition traits in Simmental and crossbred Holstein cattle

The aim of this study was to evaluate the effects of polymorphisms in the acylCoA diacylglycerol acyltransferase (DGAT1), fatty acid synthase (FASN) and prolactin (PRL) genes on milk production traits and the potential influence of the PRL gene on SCC of Simmental and crossbred Holstein cattle as well as comparatively among them. A total of eighty-three cows were genotyped by PCR-RFLP method. Crossbred Holsteins with DGAT1 genotype KA hadsignificantly higher fat content than cows with AA genotype and for FASN diplotype TW/AR had significantly higher fat content compared to diplotype AR/AR. In the observed group of animals, there was no association found between PRL SNP G8398A and milk production traits and SCC (SCS). For all three genes, cross bred Holstein hadsignificantly higher milk yield, fat and protein content in certain genotypes/diplotypes. The results of combined genotype analysis of the three SNPs showed that GG/TWAR/KK genotype combination had the lowest SCS.

Searching for new loci and candidate genes for economically important traits through gene-based association analysis of Simmental cattle

Single-marker genome-wide association study (GWAS) is a convenient strategy of genetic analysis that has been successful in detecting the association of a number of single-nucleotide polymorphisms (SNPs) with quantitative traits. However, analysis of individual SNPs can only account for a small proportion of genetic variation and offers only limited knowledge of complex traits. This inadequacy may be overcome by employing a gene-based GWAS analytic approach, which can be considered complementary to the single-SNP association analysis. Here we performed an initial single-SNP GWAS for bone weight (BW) and meat pH value with a total of 770,000 SNPs in 1141 Simmental cattle. Additionally, 21836 cattle genes collected from the Ensembl Genes 83 database were analyzed to find supplementary evidence to support the importance of gene-based association study. Results of the single SNP-based association study showed that there were 11 SNPs significantly associated with bone weight (BW) and two SNPs associated with meat pH value. Interestingly, all of these SNPs were located in genes detected by the gene-based association study.

Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight

Background

Body weight (BW) is an important trait for meat production in sheep. Although over the past few years, numerous quantitative trait loci (QTL) have been detected for production traits in cattle, few QTL studies have been reported for sheep, with even fewer on meat production traits. Our objective was to perform a genome-wide association study (GWAS) with the medium-density Illumina Ovine SNP50 BeadChip to identify genomic regions and corresponding haplotypes associated with BW in Australian Merino sheep.

Methods

A total of 1781 Australian Merino sheep were genotyped using the medium-density Illumina Ovine SNP50 BeadChip. Among the 53 862 single nucleotide polymorphisms (SNPs) on this array, 48 640 were used to perform a GWAS using a linear mixed model approach. Genotypes were phased with hsphase; to estimate SNP haplotype effects, linkage disequilibrium blocks were identified in the detected QTL region.

Results

Thirty-nine SNPs were associated with BW at a Bonferroni-corrected genome-wide significance threshold of 1 %. One region on sheep (Ovis aries) chromosome 6 (OAR6) between 36.15 and 38.56 Mb, included 13 significant SNPs that were associated with BW; the most significant SNP was OAR6_41936490.1 (P = 2.37 × 10⁻¹⁶) at 37.69 Mb with an allele substitution effect of 2.12 kg, which corresponds to 0.248 phenotypic standard deviations for BW. The region that surrounds this association signal on OAR6 contains three genes: leucine aminopeptidase 3 (LAP3), which is involved in the processing of the oxytocin precursor; NCAPG non-SMC condensin I complex, subunit G (NCAPG), which is associated with foetal growth and carcass size in cattle; and ligand dependent nuclear receptor corepressor-like (LCORL), which is associated with height in humans and cattle.

Conclusions

The GWAS analysis detected 39 SNPs associated with BW in sheep and a major QTL region was identified on OAR6. In several other mammalian species, regions that are syntenic with this region have been found to be associated with body size traits, which may reflect that the underlying biological mechanisms share a common ancestry. These findings should facilitate the discovery of causative variants for BW and contribute to marker-assisted selection.

Electronic supplementary material

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

An interpretive review of selective sweep studies in Bos taurus cattle populations: identification of unique and shared selection signals across breeds

This review compiles the results of 21 genomic studies of European Bos taurus breeds and thus provides a general picture of the selection signatures in taurine cattle identified by genome-wide selection-mapping scans. By performing a comprehensive summary of the results reported in the literature, we compiled a list of 1049 selection sweeps described across 37 cattle breeds (17 beef breeds, 14 dairy breeds, and 6 dual-purpose breeds), and four different beef-vs.-dairy comparisons, which we subsequently grouped into core selective sweep (CSS) regions, defined as consecutive signals within 1 Mb of each other. We defined a total of 409 CSSs across the 29 bovine autosomes, 232 (57%) of which were associated with a single-breed (Single-breed CSSs), 134 CSSs (33%) were associated with a limited number of breeds (Two-to-Four-breed CSSs) and 39 CSSs (9%) were associated with five or more breeds (Multi-breed CSSs). For each CSS, we performed a candidate gene survey that identified 291 genes within the CSS intervals (from the total list of 5183 BioMart-extracted genes) linked to dairy and meat production, stature, and coat color traits. A complementary functional enrichment analysis of the CSS positional candidates highlighted other genes related to pathways underlying behavior, immune response, and reproductive traits. The Single-breed CSSs revealed an over-representation of genes related to dairy and beef production, this was further supported by over-representation of production-related pathway terms in these regions based on a functional enrichment analysis. Overall, this review provides a comparative map of the selection sweeps reported in European cattle breeds and presents for the first time a characterization of the selection sweeps that are found in individual breeds. Based on their uniqueness, these breed-specific signals could be considered as “divergence signals,” which may be useful in characterizing and protecting livestock genetic diversity.

Designing of an artificial neural network model to evaluate the association of three combined Y-specific microsatellite loci on the actual and predicted postthaw motility in crossbred bull semen

The freezing of bull semen significantly hamper the motility of sperm which reduces the conception rate in dairy cattle. The prediction of postthaw motility (PTM) before freezing will be useful to take the decision on discarding or freezing of the germplasm. The artificial neural network (ANN) methodology found to be useful in prediction and classification problems related to animal science, and hence, the present study was undertaken to compare the efficiency of ANN in prediction of PTM on the basis of the number of ejaculates, volume, and concentration of sperms. The combined effect of Y-specific microsatellite alleles on the actual and predicted PTM was also studied. The results revealed that the prediction accuracy of PTM based on the semen quality parameters was comparatively lower because of higher variability in the data set. The ANN gave better prediction accuracy (34.88%) than the multiple regression analysis models (32.04%). The root mean square error was lower for ANN (8.4353) than that in the multiple regression analysis (8.6168). The haplotype or combined effect of microsatellite alleles on actual and predicted PTM was found to be highly significant (P < 0.01). On the basis of results, it was concluded that the ANN methodology can be used for prediction of PTM in crossbred bulls.

QTLs associated with dry matter intake, metabolic mid-test weight, growth and feed efficiency have little overlap across 4 beef cattle studies

Background

The identification of genetic markers associated with complex traits that are expensive to record such as feed intake or feed efficiency would allow these traits to be included in selection programs. To identify large-effect QTL, we performed a series of genome-wide association studies and functional analyses using 50 K and 770 K SNP genotypes scored in 5,133 animals from 4 independent beef cattle populations (Cycle VII, Angus, Hereford and Simmental × Angus) with phenotypes for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake.

Results

A total of 5, 6, 11 and 10 significant QTL (defined as 1-Mb genome windows with Bonferroni-corrected P-value <0.05) were identified for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake, respectively. The identified QTL were population-specific and had little overlap across the 4 populations. The pleiotropic or closely linked QTL on BTA 7 at 23 Mb identified in the Angus population harbours a promising candidate gene ACSL6 (acyl-CoA synthetase long-chain family member 6), and was the largest effect QTL associated with dry matter intake and mid-test body weight explaining 10.39% and 14.25% of the additive genetic variance, respectively. Pleiotropic or closely linked QTL associated with average daily gain and mid-test body weight were detected on BTA 6 at 38 Mb and BTA 7 at 93 Mb confirming previous reports. No QTL for residual feed intake explained more than 2.5% of the additive genetic variance in any population. Marker-based estimates of heritability ranged from 0.21 to 0.49 for residual feed intake across the 4 populations.

Conclusions

This GWAS study, which is the largest performed for feed efficiency and its component traits in beef cattle to date, identified several large-effect QTL that cumulatively explained a significant percentage of additive genetic variance within each population. Differences in the QTL identified among the different populations may be due to differences in power to detect QTL, environmental variation, or differences in the genetic architecture of trait variation among breeds. These results enhance our understanding of the biology of growth, feed intake and utilisation in beef cattle.

Association analyses of FGFR2 gene polymorphisms with femoral neck bone mineral density in Chinese Han population

Femoral neck (FN) bone mineral density (BMD) is the most important risk phenotype for osteoporosis and has been used as a reference standard for describing osteoporosis. Identification of genetic variations associated with FN BMD may provide potential targets for therapeutic studies. Given the important biological role of FGFR2 gene involved in bone, we tested the associations between FGFR2 polymorphisms and FN BMD in 1,300 Chinese Han subjects. Of the 28 total SNPs, 2 SNPs, namely rs11200014 and rs1078806, were significantly associated with FN BMD under dominant model (P = 0.0014 and 0.0012, respectively) after conservative Bonferroni correction. The two SNPs were in complete linkage disequilibrium. In addition, haplotype-based association tests identified two haplotypes significantly associated with FN BMD, including one haplotype in block 4 where the two SNPs located. However, different from previous studies in white older men, we did not detect any significant association in sex-stratified analyses. In summary, our findings suggest that the FGFR2 gene may play an important role in variation in FN BMD in Chinese Han population, independent of gender effects. Further studies performed in multiple and large samples are needed to elucidate the underlying molecular mechanism and pathophysiology of osteoporosis.

Large-effect pleiotropic or closely linked QTL segregate within and across ten US cattle breeds

Background

The availability of high-density SNP assays including the BovineSNP50 (50 K) enables the identification of novel quantitative trait loci (QTL) and improvement of the resolution of the locations of previously mapped QTL. We performed a series of genome-wide association studies (GWAS) using 50 K genotypes scored in 18,274 animals from 10 US beef cattle breeds with observations for twelve body weights, calving ease and carcass traits.

Results

A total of 159 large-effects QTL (defined as 1-Mb genome windows explaining more than 1% of additive genetic variance) were identified. In general, more QTL were identified in analyses with bigger sample sizes. Four large-effect pleiotropic or closely linked QTLs located on BTA6 at 37–42 Mb (primarily at 38 Mb), on BTA7 at 93 Mb, on BTA14 at 23–26 Mb (primarily at 25 Mb) and on BTA20 at 4 Mb were identified in more than one breed. Several breed-specific large-effect pleiotropic or closely linked QTL were also identified. Some identified QTL regions harbor genes known to have large effects on a variety of traits in cattle such as PLAG1 and MSTN and others harbor promising candidate genes including NCAPG, ARRDC3, ERGIC1, SH3PXD2B, HMGA2, MSRB3, LEMD3, TIGAR, SEPT7, and KIRREL3. Gene ontology analysis revealed that genes involved in ossification and in adipose tissue development were over-represented in the identified pleiotropic QTL. Also, the MAPK signaling pathway was identified as a common pathway affected by the genes located near the pleiotropic QTL.

Conclusions

This largest GWAS ever performed in beef cattle, led us to discover several novel across-breed and breed-specific large-effect pleiotropic QTL that cumulatively account for a significant percentage of additive genetic variance (e.g. more than a third of additive genetic variance of birth and mature weights; and calving ease direct in Hereford). These results will improve our understanding of the biology of growth and body composition in cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-442) contains supplementary material, which is available to authorized users.

Signatures of selection in five Italian cattle breeds detected by a 54K SNP panel

In this study we used a medium density panel of SNP markers to perform population genetic analysis in five Italian cattle breeds. The BovineSNP50 BeadChip was used to genotype a total of 2,935 bulls of Piedmontese, Marchigiana, Italian Holstein, Italian Brown and Italian Pezzata Rossa breeds. To determine a genome-wide pattern of positive selection we mapped the F st values against genome location. The highest F st peaks were obtained on BTA6 and BTA13 where some candidate genes are located. We identified selection signatures peculiar of each breed which suggest selection for genes involved in milk or meat traits. The genetic structure was investigated by using a multidimensional scaling of the genetic distance matrix and a Bayesian approach implemented in the STRUCTURE software. The genotyping data showed a clear partitioning of the cattle genetic diversity into distinct breeds if a number of clusters equal to the number of populations were given. Assuming a lower number of clusters beef breeds group together. Both methods showed all five breeds separated in well defined clusters and the Bayesian approach assigned individuals to the breed of origin. The work is of interest not only because it enriches the knowledge on the process of evolution but also because the results generated could have implications for selective breeding programs.

A genome-wide scan for signatures of differential artificial selection in ten cattle breeds

Background

Since the times of domestication, cattle have been continually shaped by the influence of humans. Relatively recent history, including breed formation and the still enduring enormous improvement of economically important traits, is expected to have left distinctive footprints of selection within the genome. The purpose of this study was to map genome-wide selection signatures in ten cattle breeds and thus improve the understanding of the genome response to strong artificial selection and support the identification of the underlying genetic variants of favoured phenotypes. We analysed 47,651 single nucleotide polymorphisms (SNP) using Cross Population Extended Haplotype Homozygosity (XP-EHH).

Results

We set the significance thresholds using the maximum XP-EHH values of two essentially artificially unselected breeds and found up to 229 selection signatures per breed. Through a confirmation process we verified selection for three distinct phenotypes typical for one breed (polledness in Galloway, double muscling in Blanc-Bleu Belge and red coat colour in Red Holstein cattle). Moreover, we detected six genes strongly associated with known QTL for beef or dairy traits (TG, ABCG2, DGAT1, GH1, GHR and the Casein Cluster) within selection signatures of at least one breed. A literature search for genes lying in outstanding signatures revealed further promising candidate genes. However, in concordance with previous genome-wide studies, we also detected a substantial number of signatures without any yet known gene content.

Conclusions

These results show the power of XP-EHH analyses in cattle to discover promising candidate genes and raise the hope of identifying phenotypically important variants in the near future. The finding of plausible functional candidates in some short signatures supports this hope. For instance, MAP2K6 is the only annotated gene of two signatures detected in Galloway and Gelbvieh cattle and is already known to be associated with carcass weight, back fat thickness and marbling score in Korean beef cattle. Based on the confirmation process and literature search we deduce that XP-EHH is able to uncover numerous artificial selection targets in subpopulations of domesticated animals.

Genome-wide association study identifies Loci and candidate genes for body composition and meat quality traits in Beijing-You chickens

Body composition and meat quality traits are important economic traits of chickens. The development of high-throughput genotyping platforms and relevant statistical methods have enabled genome-wide association studies in chickens. In order to identify molecular markers and candidate genes associated with body composition and meat quality traits, genome-wide association studies were conducted using the Illumina 60 K SNP Beadchip to genotype 724 Beijing-You chickens. For each bird, a total of 16 traits were measured, including carcass weight (CW), eviscerated weight (EW), dressing percentage, breast muscle weight (BrW) and percentage (BrP), thigh muscle weight and percentage, abdominal fat weight and percentage, dry matter and intramuscular fat contents of breast and thigh muscle, ultimate pH, and shear force of the pectoralis major muscle at 100 d of age. The SNPs that were significantly associated with the phenotypic traits were identified using both simple (GLM) and compressed mixed linear (MLM) models. For nine of ten body composition traits studied, SNPs showing genome wide significance (P<2.59E-6) have been identified. A consistent region on chicken (Gallus gallus) chromosome 4 (GGA4), including seven significant SNPs and four candidate genes (LCORL, LAP3, LDB2, TAPT1), were found to be associated with CW and EW. Another 0.65 Mb region on GGA3 for BrW and BrP was identified. After measuring the mRNA content in beast muscle for five genes located in this region, the changes in GJA1 expression were found to be consistent with that of breast muscle weight across development. It is highly possible that GJA1 is a functional gene for breast muscle development in chickens. For meat quality traits, several SNPs reaching suggestive association were identified and possible candidate genes with their functions were discussed.

Identification of a short region on chromosome 6 affecting direct calving ease in Piedmontese cattle breed

Calving in cattle is affected by calf morphology and by dam characteristics. It is described by two different traits: maternal calving ease, which is the ability to generate dams with good physiological predisposition to calving, and direct calving ease, which is the ability to generate calves that are easily born. The aim of this study was to identify regions of cattle genome harboring genes possibly affecting direct calving ease in the Piedmontese cattle breed. A population of 323 bulls scored for direct calving ease (EBV) was analyzed by a medium-density SNP marker panel (54,001 SNPs) to perform a genome-wide scan. The strongest signal was detected on chromosome 6 between 37.8 and 38.7 Mb where 13 SNPs associated to direct calving ease were found. Three genes are located in this region: LAP3, encoding for a leucine aminopeptidase involved in the oxytocin hydrolysis; NCAPG, encoding for a non-SMC condensin I complex, which has been associated in cattle with fetal growth and carcass size; and LCORL, which has been associated to height in humans and cattle. To further confirm the results of the genome-wide scan we genotyped additional SNPs within these genes and analyzed their association with direct calving ease. The results of this additional analysis fully confirmed the findings of the GWAS and particularly indicated LAP3 as the most probable gene involved. Linkage Disequilibrium (LD) analysis showed high correlation between SNPs located within LAP3 and LCORL indicating a possible selection signature due either to increased fitness or breeders' selection for the trait.

Functional characterization of genetic polymorphisms identified in the promoter region of the bovine PEPS gene

Peptidase S (PEPS) is a metallopeptidase that cleaves N-terminal residues from proteins and peptides. PEPS is used as a cell maintenance enzyme with critical roles in peptide turnover. The promoter region located upstream of the initiation site plays an important role in regulating gene expression. Polymorphism in the promoter region can alter gene expression and lead to biological changes. In the current study, polymorphisms in the promoter region of the PEPS gene were investigated. Polymerase chain reaction (PCR)-restriction fragment length polymorphism and DNA sequencing methods were used to screen sequence variations in the promoter region of DNA samples from 743 Chinese Holstein cattle. Two polymorphisms (g. -534 T>C and g. -2545 G>A) were identified and eight haplotypes were classified by haplotype analysis. The two genetic polymorphisms and haplotypes were associated with fat percentage and somatic cell score in Chinese Holstein cattle. The results of real-time PCR showed that cow kidneys exhibit the highest PEPS expression level. Moreover, bioinformatics analysis predicted that the single-nucleotide polymorphism g. -534 T>C is located in the core promoter region and in the transcription factor binding sites. The promoter activities of the polymorphism of -543 T>C were measured by luciferase assay in the human kidney epithelial cell line 293T. Transcriptional activity is significantly lower in cell lines transfected with the reporter construct containing 2.5 kb upstream fragments with -543 C than in those with wild-type -543 T. The results indicated that genetic variation at locus -543 influences PEPS promoter activity. The genetic variation in the promoter region of PEPS gene may regulate PEPS gene transcription and might have consequences at a regulatory level.

BRCA1: a new candidate gene for bovine mastitis and its association analysis between single nucleotide polymorphisms and milk somatic cell score

Bovine mastitis is a very complex and common disease of dairy cattle and a major source of economic losses to the dairy industry worldwide. In this study, the bovine breast cancer 1, early onset gene (BRCA1) was taken as a candidate gene for mastitis resistance. The main object of this study was to investigate whether the BRCA1 gene was associated with mastitis in cattle. Through DNA sequencing, Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP) and Created Restriction Site PCR (CRS-PCR) methods, three SNPs (G22231T, T25025A, and C28300A) were detected and twenty-four combinations of these SNPs were observed. The single SNP and their genetic effects on somatic cell score (SCS) were evaluated and a significant association with SCS was found in C28300A. The mean of genotype EE was significantly lower than those of genotypes EF and FF. The results of combined genotypes analysis of three SNPs showed that BBDDFF genotype with the highest SCS were easily for the mastitis susceptibility, whereas AACCEE genotype with the lowest SCS were favorable for the mastitis resistance. The information provided in the present study will be very useful for improving mastitis resistance in dairy cattle by marker-assisted selection.