Genome-wide association study in Chinese Holstein cows reveal two candidate genes for somatic cell score as an indicator for mastitis susceptibility

Genome-Wide Association Study for Milk Somatic Cell Score in Holstein Friesian Cows in Slovenia

Mastitis is a serious challenge for the dairy industry, leading to economic losses and affecting milk quality. The aim of this study is to identify genetic factors associated with mastitis resistance by conducting a genome-wide association study (GWAS) for the somatic cell score (SCS). Phenotypic records of 350 Holstein Friesian cows were obtained from the Slovenian Cattle Recording Scheme Database and consisted of around 1500 lactation data from 2012 to 2023 collected on a single farm in Slovenia. Corresponding genotypic data were also retrieved from the same database and genotyped using the Illumina BovineSNP50 BeadChip (Illumina, Inc., San Diego, CA, USA). For the association study, three SCS parameters were considered, including lactation mean somatic cell score (LM_SCS), maximum SCS value (SCSMAX), and top three mean value of SCS (TOP3). After performing a GWAS using FarmCPU and BLINK models, five significant SNPs associated with the TOP3 trait were found on BTA 14, 15, 22, and 29. The identified SNP markers were closely linked to six known candidate genes (DNASE1L3, SLC36A4, ARMC1, PDE7A, MMP13, CD44). These results indicate potential genetic markers associated with SCS in the Slovenian Holstein Friesian population.

African Local Pig Genetic Resources in the Context of Climate Change Adaptation

Africa is home to a wide diversity of locally adapted pig breeds whose genetic architecture offers important insights into livestock adaptation to climate change. However, the majority of these inherent traits have not been fully highlighted. This review presents an overview of the current state of African pig genetic resources, providing highlights on their population and production statistics, production system, population diversity indices, and genomic evidence underlying their evolutionary potential. The study results reveal an incomplete characterization of local pig genotypes across the continent. The characterized population, however, demonstrates moderate to high levels of genetic diversity, enough to support breeding and conservation programs. Owing to low genetic differentiation and limited evidence of distinct population structures, it appears that most local pig populations are strains within larger breeds. Genomic evidence has shown a higher number of selection signatures associated with various economically important traits, thus making them potential candidates for climate change adaptation. The reportedly early evidence of hybridization with wild suid groups further suggests untapped insights into disease resistance and resilience traits that need to be illuminated using higher-density markers. Nevertheless, gene introgression from commercial breeds is prevalent across Africa; thus, efforts to realize and utilize these traits must increase before they are permanently depleted.

Genome-Wide Association Study Reveals Quantitative Trait Loci and Candidate Genes Associated with High Interferon-gamma Production in Holstein Cattle Naturally Infected with Mycobacterium Bovis

Mycobacterium bovis (Mb) is the causative agent of bovine tuberculosis (bTb). Genetic selection aiming to identify less susceptible animals has been proposed as a complementary measure in ongoing programs toward controlling Mb infection. However, individual animal phenotypes for bTb based on interferon-gamma (IFNɣ) and its use in bovine selective breeding programs have not been explored. In the current study, IFNɣ production was measured using a specific IFNɣ ELISA kit in bovine purified protein derivative (bPPD)-stimulated blood samples collected from Holstein cattle. DNA isolated from the peripheral blood samples collected from the animals included in the study was genotyped with the EuroG Medium Density bead Chip, and the genotypes were imputed to whole-genome sequences. A genome-wide association analysis (GWAS) revealed that the IFNɣ in response to bPPD was associated with a specific genetic profile (heritability = 0.23) and allowed the identification of 163 SNPs, 72 quantitative trait loci (QTLs), 197 candidate genes, and 8 microRNAs (miRNAs) associated with this phenotype. No negative correlations between this phenotype and other phenotypes and traits included in the Spanish breeding program were observed. Taken together, our results define a heritable and distinct immunogenetic profile associated with strong production of IFNɣ in response to Mb.

Genome-wide association study for stayability at different calvings in Nellore beef cattle

BACKGROUNDING

Stayability, which may be defined as the probability of a cow remaining in the herd until a reference age or at a specific number of calvings, is usually measured late in the animal's life. Thus, if used as selection criteria, it will increase the generation interval and consequently might decrease the annual genetic gain. Measuring stayability at an earlier age could be a reasonable strategy to avoid this problem. In this sense, a better understanding of the genetic architecture of this trait at different ages and/or at different calvings is important. This study was conducted to identify possible regions with major effects on stayability measured considering different numbers of calvings in Nellore cattle as well as pathways that can be involved in its expression throughout the female's productive life.

RESULTS

The top 10 most important SNP windows explained, on average, 17.60% of the genetic additive variance for stayability, varying between 13.70% (at the eighth calving) and 21% (at the fifth calving). These SNP windows were located on 17 chromosomes (1, 2, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 18, 19, 20, 27, and 28), and they harbored a total of 176 annotated genes. The functional analyses of these genes, in general, indicate that the expression of stayability from the second to the sixth calving is mainly affected by genetic factors related to reproductive performance, and nervous and immune systems. At the seventh and eighth calvings, genes and pathways related to animal health, such as density bone and cancer, might be more relevant.

CONCLUSION

Our results indicate that part of the target genomic regions in selecting for stayability at earlier ages (from the 2th to the 6th calving) would be different than selecting for this trait at later ages (7th and 8th calvings). While the expression of stayability at earlier ages appeared to be more influenced by genetic factors linked to reproductive performance together with an overall health/immunity, at later ages genetic factors related to an overall animal health gain relevance. These results support that selecting for stayability at earlier ages (perhaps at the second calving) could be applied, having practical implications in breeding programs since it could drastically reduce the generation interval, accelerating the genetic progress.

Genomic signatures of selection, local adaptation and production type characterisation of East Adriatic sheep breeds

BACKGROUND

The importance of sheep breeding in the Mediterranean part of the eastern Adriatic has a long tradition since its arrival during the Neolithic migrations. Sheep production system is extensive and generally carried out in traditional systems without intensive systematic breeding programmes for high uniform trait production (carcass, wool and milk yield). Therefore, eight indigenous Croatian sheep breeds from eastern Adriatic treated here as metapopulation (EAS), are generally considered as multipurpose breeds (milk, meat and wool), not specialised for a particular type of production, but known for their robustness and resistance to certain environmental conditions. Our objective was to identify genomic regions and genes that exhibit patterns of positive selection signatures, decipher their biological and productive functionality, and provide a "genomic" characterization of EAS adaptation and determine its production type.

RESULTS

We identified positive selection signatures in EAS using several methods based on reduced local variation, linkage disequilibrium and site frequency spectrum (eROHi, iHS, nSL and CLR). Our analyses identified numerous genomic regions and genes (e.g., desmosomal cadherin and desmoglein gene families) associated with environmental adaptation and economically important traits. Most candidate genes were related to meat/production and health/immune response traits, while some of the candidate genes discovered were important for domestication and evolutionary processes (e.g., HOXa gene family and FSIP2). These results were also confirmed by GO and QTL enrichment analysis.

CONCLUSIONS

Our results contribute to a better understanding of the unique adaptive genetic architecture of EAS and define its productive type, ultimately providing a new opportunity for future breeding programmes. At the same time, the numerous genes identified will improve our understanding of ruminant (sheep) robustness and resistance in the harsh and specific Mediterranean environment.

Karrow N, Mao Y. Identification and Classification of Long Non-Coding RNAs in the Mammary Gland of the Holstein Cow

The mammary glands, responsible for milk secretion, are regulated at a local level by various hormones, growth factors, non-coding RNAs, and other elements. Recent research has discovered the presence of lncRNAs in these glands, with suggestions that they may be essential for the maintenance and function of mammary glands. Besides directly controlling the gene and protein expression, lncRNAs are believed to play a significant part in numerous physiological and pathological processes. This study focused on examining the mammary gland tissues of Chinese Holstein cows, to identify and categorize long non-coding RNAs (lncRNAs). The research intended to distinguish lncRNAs in the mammary tissues of Holstein cows and contrast them between lactation and non-lactation periods. In this study, mammary gland tissues were sampled from three Holstein cows in early lactation (n = 3, 30 days postpartum) and non-lactation (n = 3, 315 days postpartum) on a large dairy farm in Jiangsu province. Mammary tissue samples were collected during early lactation and again during non-lactation. In total, we detected 1905 lncRNAs, with 57.3% being 500 bp and 612 intronic lncRNAs. The exon count for lncRNAs varied from 2 to 10. It was observed that 96 lncRNA expressions markedly differed between the two stages, with 83 genes being upregulated and 53 downregulated. Enrichment analysis results revealed that Gene Ontology (GO) analysis was primarily abundant in cellular processes. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that target genes were predominantly abundant in metabolic pathways, fatty acid biosynthesis, the immune system, and glycosphingolipid biosynthesis. This study analyzed the expression profile and characteristics of lncRNAs in the mammary gland tissues of Holstein cows during both lactation and non-lactation stages, forming a foundation for further investigation into the functional roles of lncRNAs in Holstein cows throughout lactation.

Kompetitive Allele Specific PCR Genotyping of 89 SNPs in Romanian Spotted and Romanian Brown Cattle Breeds and Their Association with Clinical Mastitis

Mastitis is the most common production disease in the dairy sector worldwide, its incidence being associated with both cows' exposure to bacteria and the cows' genetic make-up for resistance to pathogens. The objective of our study was to analyse 89 missense SNPs belonging to six genes (CXCR2, CXCL8, TLR4, BRCA1, LTF, BOLA-DRB3), which were found to be associated with genetic resistance or susceptibility to mastitis. A total of 298 cattle (250 Romanian Spotted and 48 Romanian Brown) were genotyped by Kompetitive Allele Specific PCR (KASP) and a chi-squared test was used for genetic association studies with clinical mastitis. A total of 35 SNPs (39.3%) among the selected 89 SNPs were successfully genotyped, of which 31 markers were monomorphic. The polymorphic markers were found in two genes: TLR4 (rs460053411) and BOLA-DRB3 (rs42309897, rs208816121, rs110124025). The polymorphic SNPs with MAF > 5% and call rates > 95% were used for the association study. The results showed that rs110124025 in the BOLA-DRB3 gene was significantly associated with mastitis prevalence (p ≤ 0.05) in both investigated breeds. Current results show that the SNP rs110124025 in the BOLA-DRB3 gene can be used as a candidate genetic marker in selection for mastitis resistance in Romanian dairy cattle.

Genome-wide association study for selected cheese-making properties in Dual-Purpose Belgian Blue cows

This study aimed to estimate genetic parameters and identify genomic region(s) associated with selected cheese-making properties (CMP) in Dual-Purpose Belgian Blue (DPBB) cows. Edited data were 46,301 test-day records of milk yield, fat percentage, protein percentage, casein percentage, milk calcium content (CC), coagulation time (CT), curd firmness after 30 min from rennet addition (a30), and milk titratable acidity (MTA) collected from 2014 to 2020 on 4,077 first-parity (26,027 test-day records), and 3,258 second-parity DPBB cows (20,274 test-day records) distributed in 124 herds in the Walloon Region of Belgium. Data of 28,266 SNP, located on 29 Bos taurus autosomes (BTA) of 1,699 animals were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 25 consecutive SNPs (with an average size of ∼2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Heritability estimates for the included CMP ranged from 0.19 (CC) to 0.50 (MTA), and 0.24 (CC) to 0.41 (MTA) in the first and second parity, respectively. The genetic correlation estimated between CT and a30 varied from -0.61 to -0.41 and from -0.55 to -0.38 in the first and second lactations, respectively. Negative genetic correlations were found between CT and milk yield and composition, while those estimated between curd firmness and milk composition were positive. Genome-wide association analyses results identified 4 genomic regions (BTA1, BTA3, BTA7, and BTA11) associated with the considered CMP. The identified genomic regions showed contrasting results between parities and among the different stages of each parity. It suggests that different sets of candidate genes underlie the phenotypic expression of the considered CMP between parities and lactation stages of each parity. The findings of this study can be used for future implementation and use of genomic evaluation to improve the cheese-making traits in DPBB cows.

Genetic polymorphisms of TRAPPC9 and CD4 genes and their association with milk production and mastitis resistance phenotypic traits in Chinese Holstein

The present study was designed to evaluate the association of polymorphisms in bovine trafficking protein particle complex subunit 9 (TRAPPC9) and cluster of differentiation 4 (CD4) genes with milk production and mastitis resistance phenotypic traits in a different cattle population. Three single nucleotide polymorphisms (SNPs) (SNP1 Position: Chr14:2484891, SNP2 (rs110017379), SNP3 Position: Chr14:2525852) in bovine TRAPPC9 and one SNP (Position: Chr5:104010752) in CD4 were screened through Chinese Cow's SNPs Chip-I (CCSC-I) and genotyped in a population of 312 Chinese Holsteins (156: Mastitis, 156: Healthy). The results were analyzed using the general linear model in SAS 9.4. Our analysis revealed that milk protein percentage, somatic cell count (SCC), somatic cell score (SCS), serum cytokines interleukin 6 (IL-6) and interferon-gamma (IFN-γ) were significantly (P < 0.05) associated with at least one or more identified SNPs of TRAPPC9 and CD4 genes. Furthermore, the expression status of SNPs in CD4 and TRAPPC9 genes were verified through RT-qPCR. The expression analysis showed that genotypes GG in SNP3 of TRAPPC9 and TT genotype in SNP4 of CD4 showed higher expression level compared to other genotypes. The GG genotype in SNP2 and TT genotype in SNP3 of TRAPPC9 were associated with higher bovine milk SCC and lower IL6. Altogether, our findings suggested that the SNPs of TRAPPC9 and CD4 genes could be useful genetic markers in selection for milk protein improvement and mastitis resistance phenotypic traits in dairy cattle. The CCSC-I used in current study is proposed to be validate in different and large population of dairy cattle not only in China but also in other countries. Moreover, our analyses recommended that besides SCC and SCS, the association of genetic markers could also be considered with the serum cytokines (IL-6, IFN-γ) while selecting genetically mastitis resistance dairy cattle.

Genome-Wide Association Study for Udder Conformation Traits in Chinese Holstein Cattle

Udder conformation traits are one of the most economic traits in dairy cows, greatly affecting animal health, milk production, and producer profitability in the dairy industry. Genetic analysis of udder structure and scores have been developed in Holstein cattle. In our research, we conducted a genome-wide association study for five udder traits, including anterior udder attachment (AUA), central suspensory ligament (CSL), posterior udder attachment height (PUAH), posterior udder attachment width (PUAW), and udder depth (UD), in which the fixed and random model circulating probability unification (FarmCPU) model was applied for the association analysis. The heritability and the standard errors of these five udder traits ranged from 0.04 ± 0.00 to 0.49 ± 0.03. Phenotype data were measured from 1000 Holstein cows, and the GeneSeek Genomic Profiler (GGP) Bovine 100 K SNP chip was used to analyze genotypic data in Holstein cattle. For GWAS analysis, 984 individual cows and 84,407 single-nucleotide polymorphisms (SNPs) remained after quality control; a total of 18 SNPs were found at the GW significant threshold (p < 5.90 × 10−7). Many candidate genes were identified within 200kb upstream or downstream of the significant SNPs, which include MGST1, MGST2, MTUS1, PRKN, STXBP6, GRID2, E2F8, CDH11, FOXP1, SLF1, TMEM117, SBF2, GC, ADGRB3, and GCLC. Pathway analysis revealed that 58 Gene Ontology (GO) terms and 18 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched with adjusted p values, and these GO terms and the KEGG pathway analysis were associated with biological information, metabolism, hormonal growth, and development processes. These results could give valuable biological information for the genetic architecture of udder conformation traits in dairy Holstein cattle.

Emerging Roles of Noncoding RNAs in Bovine Mastitis Diseases

Non-coding RNAs (ncRNAs) are an abundant class of RNA with varying nucleotide lengths. They have been shown to have great potential in eutherians/human disease diagnosis and treatments and are now gaining more importance for the improvement of diseases in livestock. To date, thousands of ncRNAs have been discovered in the bovine genome and the continuous advancement in deep sequencing technologies and various bioinformatics tools has enabled the elucidation of their roles in bovine health. Among farm animals' diseases, mastitis, a common inflammatory disease in cattle, has caused devastating economic losses to dairy farmers over the last few decades. Here, we summarize the biology of bovine mastitis and comprehensively discuss the roles of ncRNAs in different types of mastitis infection. Based on our findings and relevant literature, we highlighted various evidence of ncRNA roles in mastitis. Different approaches (in vivo versus in vitro) for exploring ncRNA roles in mastitis are emphasized. More particularly, the potential applications of emerging genome editing technologies, as well as integrated omics platforms for ncRNA studies and implications for mastitis are presented.

Identification of TRAPPC9 and BAIAP2 Gene Polymorphisms and Their Association With Fat Deposition-Related Traits in Hu Sheep

Fat deposition is an important economic trait that is closely related to feed efficiency and carcass performance in livestock. In this study, the fat deposition-related traits of 1,293 Hu sheep were measured and descriptive statistical analysis was conducted. The results showed that the coefficient of variation of all fat deposition-related traits was higher than 24%. In addition, single nucleotide polymorphisms and the expression characteristics of TRAPPC9 (encoding trafficking protein particle complex subunit 9) and BAIAP2 (encoding brain-specific Angiogenesis inhibitor 1-associated protein 2) genes in Hu sheep were detected using PCR amplification, Sanger sequencing, KASPar genotyping, and quantitative real-time reverse transcription PCR (qRT-PCR). The associations between SNPs and fat deposition-related traits were also analyzed. Two intronic mutations, TRAPPC9 g.57654 A > G and BAIAP2 g.46061 C > T, were identified in Hu sheep. The result of association analysis showed that TRAPPC9 g.57654 A > G and BAIAP2 g.46061 C > T were both significantly associated with the weight of tail fat, tail fat relative weight (body weight), and tail fat relative weight (carcass) (P < 0.05). Comprehensive effects analysis showed that there were significant differences between the combined genotypes and tail fat and perirenal fat deposition. Moreover, qRT-PCR analysis showed that TRAPPC9 and BAIAP2 are widely expressed, and their expression levels were significantly higher in the small-tail group compared with those in the big-tail group (P < 0.01). These results provided important candidate molecular markers that could be used in strategies to reduce tail fat deposition in Hu sheep.

Genome-wide associations for heat stress response suggest potential candidate genes underlying milk fatty acid composition in dairy cattle

Contents of milk fatty acids (FA) display remarkable alterations along climatic gradients. Detecting candidate genes underlying such alterations might be beneficial for the exploration of climate sensitivity in dairy cattle. Consequently, we aimed on the definition of FA heat stress indicators, considering FA breeding values in response to temperature-humidity index (THI) alterations. Indicators were used in GWAS, in ongoing gene annotations and for the estimation of chromosome-wide variance components. The phenotypic data set consisted of 39,600 test-day milk FA records from 5,757 first-lactation Holstein dairy cows kept in 16 large-scale German cooperator herds. The FA traits were C18:0, polyunsaturated fatty acids (PUFA), saturated fatty acids (SFA), and unsaturated fatty acids (UFA). After genotype quality control, 40,523 SNP markers from 3,266 cows and 930 sires were considered. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI, which were allocated to 10 different THI classes. The same FA from 3 stages of lactation were considered as different, but genetically correlated traits. Consequently, a 3-trait reaction norm model was used to estimate genetic parameters and breeding values for FA along THI classes, considering either pedigree (A) or genomic (G) relationship matrices. De-regressed proofs and genomic estimated breeding values at the intermediate THI class 5 and at the extreme THI class 10 were used as pseudophenotypes in ongoing genomic analyses for thermoneutral (TNC) and heat stress conditions (HSC), respectively. The differences in de-regressed proofs and in genomic estimated breeding values from both THI classes were pseudophenotypes for heat stress response (HSR). Genetic correlations between the same FA under TNC and HSC were smallest in the first lactation stage and ranged from 0.20 for PUFA to 0.87 for SFA when modeling with the A matrix, and from 0.35 for UFA to 0.86 for SFA when modeling with the G matrix. In the first lactation stage, larger additive genetic variances under HSC compared with TNC indicate climate sensitivity for C18:0, PUFA, and UFA. Climate sensitivity was also reflected by pronounced chromosome-wide genetic variances for HSR of PUFA and UFA in the first stage of lactation. For all FA under TNC, HSC, and HSR, quite large genetic variance proportions were explained by BTA14. In GWAS, 30 SNP (within or close to 38 potential candidate genes) overlapped for HSR of the different FA. One unique potential candidate gene (AMFR) was detected for HSR of PUFA, 15 for HSR of SFA (ADGRB1, DENND3, DUSP16, EFR3A, EMP1, ENSBTAG00000003838, EPS8, MGP, PIK3C2G, STYK1, TMEM71, GSG1, SMARCE1, CCDC57, and FASN) and 3 for HSR of UFA (ENSBTAG00000048091, PAEP, and EPPK1). The identified unique genes play key roles in milk FA synthesis and are associated with disease resistance in dairy cattle. The results suggest consideration of FA in combination with climatic responses when inferring genetic mechanisms of heat stress in dairy cows.

Multi-Omics Integration and Network Analysis Reveal Potential Hub Genes and Genetic Mechanisms Regulating Bovine Mastitis

Mastitis, inflammation of the mammary gland, is the most prevalent disease in dairy cattle that has a potential impact on profitability and animal welfare. Specifically designed multi-omics studies can be used to prioritize candidate genes and identify biomarkers and the molecular mechanisms underlying mastitis in dairy cattle. Hence, the present study aimed to explore the genetic basis of bovine mastitis by integrating microarray and RNA-Seq data containing healthy and mastitic samples in comparative transcriptome analysis with the results of published genome-wide association studies (GWAS) using a literature mining approach. The integration of different information sources resulted in the identification of 33 common and relevant genes associated with bovine mastitis. Among these, seven genes—CXCR1, HCK, IL1RN, MMP9, S100A9, GRO1, and SOCS3—were identified as the hub genes (highly connected genes) for mastitis susceptibility and resistance, and were subjected to protein-protein interaction (PPI) network and gene regulatory network construction. Gene ontology annotation and enrichment analysis revealed 23, 7, and 4 GO terms related to mastitis in the biological process, molecular function, and cellular component categories, respectively. Moreover, the main metabolic-signalling pathways responsible for the regulation of immune or inflammatory responses were significantly enriched in cytokine–cytokine-receptor interaction, the IL-17 signaling pathway, viral protein interaction with cytokines and cytokine receptors, and the chemokine signaling pathway. Consequently, the identification of these genes, pathways, and their respective functions could contribute to a better understanding of the genetics and mechanisms regulating mastitis and can be considered a starting point for future studies on bovine mastitis.

Genome-Wide Association Study Candidate Genes on Mammary System-Related Teat-Shape Conformation Traits in Chinese Holstein Cattle

In the dairy industry, mammary system traits are economically important for dairy animals, and it is important to explain their fundamental genetic architecture in Holstein cattle. Good and stable mammary system-related teat traits are essential for producer profitability in animal fitness and in the safety of dairy production. In this study, we conducted a genome-wide association study on three traits—anterior teat position (ATP), posterior teat position (PTP), and front teat length (FTL)—in which the FarmCPU method was used for association analyses. Phenotypic data were collected from 1000 Chinese Holstein cattle, and the GeneSeek Genomic Profiler Bovine 100K single-nucleotide polymorphisms (SNP) chip was used for cattle genotyping data. After the quality control process, 984 individual cattle and 84,406 SNPs remained for GWAS work analysis. Nine SNPs were detected significantly associated with mammary-system-related teat traits after a Bonferroni correction (p < 5.92 × 10⁻⁷), and genes within a region of 200 kb upstream or downstream of these SNPs were performed bioinformatics analysis. A total of 36 gene ontology (GO) terms and 3 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly enriched (p < 0.05), and these terms and pathways are mainly related to metabolic processes, immune response, and cellular and amino acid catabolic processes. Eleven genes including MMS22L, E2F8, CSRP3, CDH11, PEX26, HAL, TAMM41, HIVEP3, SBF2, MYO16 and STXBP6 were selected as candidate genes that might play roles in the teat traits of cows. These results identify SNPs and candidate genes that give helpful biological information for the genetic architecture of these teat traits, thus contributing to the dairy production, health, and genetic selection of Chinese Holstein cattle.

Genetic Diversity and Population Structure for Resistance and Susceptibility to Mastitis in Braunvieh Cattle

Mastitis is a disease that causes significant economic losses, since resistance to mastitis is a difficult trait to be improved due to its multifactorial occurrence. Therefore, our objective was to characterize a Mexican Braunvieh cattle population for genetic resistance and susceptibility to mastitis. We used 66 SNP markers for 45 candidate genes in 150 animals. The average heterozygosity was 0.445 ± 0.076, a value higher than those reported for some European breeds. The inbreeding coefficient was slightly negative for resistance to subclinical (−0.058 ± 0.055) and clinical (−0.034 ± 0.076) mastitis, possibly due to low selection for the immunological candidate genes that influence these traits. The genotypic profiles for the candidate loci per K-means group were obtained, as well as the group distribution through the graphics of the principal component analysis. The genotypic profiles showed high genetic diversity among groups. Resistance to clinical mastitis had the lowest presence of the heterozygous genotypes. Although the percentage of highly inbred animals (>50%) is up to 13.3%, there are highly heterozygous groups in terms of the studied traits, a favorable indicator of the presence of genetic diversity. The results of this study constitute evidence of the genetic potential of the Mexican Braunvieh population to improve mastitis-related traits.

Understanding the genomic architecture of clinical mastitis in Bos indicus

This study elucidated potential genetic variants and QTLs associated with clinical mastitis incidence traits in Bos indicus breed, Sahiwal. Estimated breeding values for the traits (calculated using Bayesian inference) were used as pseudo-phenotypes for association with genome-wide SNPs and further QTL regions underlying the traits were identified. In all, 25 SNPs were found to be associated with the traits at the genome-wide suggestive threshold (p ≤ 5 × 10^-4) and these SNPs were used to define QTL boundaries based on the linkage disequilibrium structure. A total of 16 QTLs were associated with the trait EBVs including seven each for clinical mastitis incidence (CMI) in first and second lactations and two for CMI in third lactation. Nine out of sixteen QTLs overlapped with the already reported QTLs for mastitis traits, whereas seven were adjudged as novel ones. Important candidates for clinical mastitis in the identified QTL regions included DNAJB9, ELMO1, ARHGAP26, NR3C1, CACNB2, RAB4A, GRB2, NUP85, SUMO2, RBPJ, and RAB33B genes. These findings shed light on the genetic architecture of the disease in Bos indicus, and present potential regions for fine mapping and downstream analysis in future.

Genome-Wide Association Studies for Milk Somatic Cell Score in Romanian Dairy Cattle

Mastitis is one of the most frequently encountered diseases in dairy cattle, negatively affecting animal welfare and milk production. For this reason, contributions to understanding its genomic architecture are of great interest. Genome-wide association studies (GWAS) have identified multiple loci associated with somatic cell score (SCS) and mastitis in cattle. However, most of the studies have been conducted in different parts of the world on various breeds, and none of the investigations have studied the genetic architecture of mastitis in Romanian dairy cattle breeds up to this point in time. In this study, we report the first GWAS for SCS in dairy cattle breeds from Romania. For GWAS, we used an Axiom Bovine v3 SNP-chip (>63,000 Single Nucleotide Polymorphism -SNPs) and 33,330 records from 690 cows belonging to Romanian Spotted (RS) and Romanian Brown (RB) cattle. The results found one SNP significantly associated with SCS in the RS breed and 40 suggestive SNPs with -log₁₀ (p) from 4 to 4.9 for RS and from 4 to 5.4 in RB. From these, 14 markers were located near 12 known genes (AKAP8, CLHC1, MEGF10, SATB2, GATA6, SPATA6, COL12A1, EPS8, LUZP2, RAMAC, IL12A and ANKRD55) in RB cattle, 3 markers were close to ZDHHC19, DAPK1 and MMP7 genes, while one SNP overlapped the HERC3 gene in RS cattle. Four genes (HERC3, LUZP2, AKAP8 and MEGF10) associated with SCS in this study were previously reported in different studies. The most significant SNP (rs110749552) associated with SCS was located within the HERC3 gene. In both breeds, the SNPs and position of association signals were distinct among the three parities, denoting that mastitis is controlled by different genes that are dependent according to parity. The current results contribute to an expansion in the body of knowledge regarding the proportion of genetic variability explained by SNPs for SCS in dairy cattle.

Mining for the association of bovine mastitis linked genes to pathological signatures and Pathways

Background: Bovine mastitis is a common infectious disease with a serious threat to the dairy industry and public health. Mastitis is a polygenetic trait under the control of many genes. In the current study, our research attempted to address the role of mastitis-associated genes in various signalings including parasitic, viral, cancer and fungal diseases by using online bioinformatics software. Methods: We selected mastitis-associated genes from already published data and using online bioinformatics tools including DAVID and String classify the pathological role of relevant genes. A Venn diagram was used to show the status of overlapping genes among different biological function processes. Result: This study revealed that the genes gathered in published resources of mastitis were significantly correlated with Influenza A, Chagas disease, Leishmaniasis, Toxoplasmosis, Tuberculosis, Cancer signaling, Hepatitis B, Type I &II diabetes mellitus and Prion diseases biological pathways. Based on our findings, we concluded that mastitis-linked genes could be used as markers for many other diseases. Moreover, the Bioinformatics tools applied in the current study might be helpful in screening the genes involved in one disease and their association with other diseases as well.

Genome-Wide Association Study Identifies Candidate Genes Associated with Feet and Leg Conformation Traits in Chinese Holstein Cattle

Simple Summary

Feet and leg problems are among the major reasons for dairy cows leaving the herd, as well as having direct association with production and reproduction efficiency, health (e.g., claw disorders and lameness) and welfare. Hence, understanding the genetic architecture underlying feet and conformation traits in dairy cattle offers new opportunities toward the genetic improvement and long-term selection. Through a genome-wide association study on Chinese Holstein cattle, we identified several candidate genes associated with feet and leg conformation traits. These results could provide useful information about the molecular breeding basis of feet and leg traits, thus improving the longevity and productivity of dairy cattle.

Abstract

Feet and leg conformation traits are considered one of the most important economical traits in dairy cattle and have a great impact on the profitability of milk production. Therefore, identifying the single nucleotide polymorphisms (SNPs), genes and pathways analysis associated with these traits might contribute to the genomic selection and long-term plan selection for dairy cattle. We conducted genome-wide association studies (GWASs) using the fixed and random model circulating probability unification (FarmCPU) method to identify SNPs associated with bone quality, heel depth, rear leg side view and rear leg rear view of Chinese Holstein cows. Phenotypic measurements were collected from 1000 individuals of Chinese Holstein cattle and the GeneSeek Genomic Profiler Bovine 100 K SNP chip was utilized for individual genotyping. After quality control, 984 individual cows and 84,906 SNPs remained for GWAS work; as a result, we identified 20 significant SNPs after Bonferroni correction. Several candidate genes were identified within distances of 200 kb upstream or downstream to the significant SNPs, including ADIPOR2, INPP4A, DNMT3A, ALDH1A2, PCDH7, XKR4 and CADPS. Further bioinformatics analyses showed 34 gene ontology terms and two signaling pathways were significantly enriched (p ≤ 0.05). Many terms and pathways are related to biological quality, metabolism and development processes; these identified SNPs and genes could provide useful information about the genetic architecture of feet and leg traits, thus improving the longevity and productivity of Chinese Holstein dairy cattle.

Genome-Wide Association Studies of Somatic Cell Count in the Assaf Breed

Simple Summary

Mastitis causes economic loss due to discarded milk and reduced milk production and quality, increased medical care costs and somatic cell count (SCC) penalties. The use of genetic markers associated with the variability of this trait through marker-assisted selection (MAS) could help traditional methods. Our objectives were to identify new single nucleotide polymorphisms (SNPs) and genes associated with mastitis resistance in Assaf sheep by using the Illumina Ovine Infinium® HD SNP BeadChip (680K). Firstly, corrected phenotype estimates for somatic cell score (SCS) were calculated using 6173 records from 1894 multiparous Assaf ewes, and were used to select 192 extreme animals (low SCS group: n = 96; and high SCS group: n = 96) for the genome-wide association study (GWAS). Four SNPs (rs419096188, rs415580501, rs410336647, and rs424642424), three of them totally linked, were found to be significant at the chromosome level (FDR 10%) in two different regions of OAR19 close to genes related to the immune system response. Validation studies of two SNPs (rs419096188 and rs424642424) by Kompetitive Allele-Specific PCR (KASP) genotyping in the total population (n = 1894) confirmed previous GWAS association results for the SCS trait. Finally, the SNP rs419096188 was also associated with lactose content trait.

Abstract

A genome-wide association study (GWAS) was performed to identify new single nucleotide polymorphisms (SNPs) and genes associated with mastitis resistance in Assaf sheep by using the Illumina Ovine Infinium^® HD SNP BeadChip (680K). In total, 6173 records from 1894 multiparous Assaf ewes with at least three test day records and aged between 2 and 7 years old were used to estimate a corrected phenotype for somatic cell score (SCS). Then, 192 ewes were selected from the top (n = 96) and bottom (n = 96) tails of the corrected SCS phenotype distribution to be used in a GWAS. Although no significant SNPs were found at the genome level, four SNPs (rs419096188, rs415580501, rs410336647, and rs424642424) were significant at the chromosome level (FDR 10%) in two different regions of OAR19. The SNP rs419096188 was located in intron 1 of the NUP210 and close to the HDAC11 genes (61 kb apart), while the other three SNPs were totally linked and located 171 kb apart from the ARPP21 gene. These three genes were related to the immune system response. These results were validated in two SNPs (rs419096188 and rs424642424) in the total population (n = 1894) by Kompetitive Allele-Specific PCR (KASP) genotyping. Furthermore, rs419096188 was also associated with lactose content.

Genome-Wide Identification of Candidate Genes for Milk Production Traits in Korean Holstein Cattle

Simple Summary

Milk production traits that are economically important in the dairy industry have been considered the main selection criteria for breeding. The present genome-wide association study was performed to identify chromosomal loci and candidate genes with potential effects on milk production phenotypes in a Korean Holstein population. A total of eight significant quantitative trait locus regions were identified for milk yield (Bos taurus autosome (BTA) 7 and 14), adjusted 305-d fat yield (BTA 3, 5, and 14), adjusted 305-d protein yield (BTA 8), and somatic cell score (BTA 8 and 23) of milk production traits. Furthermore, we discovered three main candidate genes (diacylglycerol O-acyltransferase 1 (DGAT1), phosphodiesterase 4B (PDE4B), and anoctamin 2 (ANO2)) through bioinformatics analysis. These genes may help to understand better the underlying genetic and molecular mechanisms for milk production phenotypes in the Korean Holstein population.

Abstract

We performed a genome-wide association study and fine mapping using two methods (single marker regression: frequentist approach and Bayesian C (BayesC): fitting selected single nucleotide polymorphisms (SNPs) in a Bayesian framework) through three high-density SNP chip platforms to analyze milk production phenotypes in Korean Holstein cattle (n = 2780). We identified four significant SNPs for each phenotype in the single marker regression model: AX-311625843 and AX-115099068 on Bos taurus autosome (BTA) 14 for milk yield (MY) and adjusted 305-d fat yield (FY), respectively, AX-428357234 on BTA 18 for adjusted 305-d protein yield (PY), and AX-185120896 on BTA 5 for somatic cell score (SCS). Using the BayesC model, we discovered significant 1-Mb window regions that harbored over 0.5% of the additive genetic variance effects for four milk production phenotypes. The concordant significant SNPs and 1-Mb window regions were characterized into quantitative trait loci (QTL). Among the QTL regions, we focused on a well-known gene (diacylglycerol O-acyltransferase 1 (DGAT1)) and newly identified genes (phosphodiesterase 4B (PDE4B), and anoctamin 2 (ANO2)) for MY and FY, and observed that DGAT1 is involved in glycerolipid metabolism, fat digestion and absorption, metabolic pathways, and retinol metabolism, and PDE4B is involved in cAMP signaling. Our findings suggest that the candidate genes in QTL are strongly related to physiological mechanisms related to the fat production and consequent total MY in Korean Holstein cattle.

Bovine mastitis prevention and control in the post-antibiotic era

Mastitis is the most important disease in the dairy industry. Antibiotics are considered to be the first choice in the treatment of the disease. However, the problem of antibiotic residue and antimicrobial resistance, in addition to the impact of antibiotic abuse on public health, leads to many restrictions on uncontrolled antibiotic therapy in the dairy sector worldwide. Researchers have investigated novel therapeutic approaches to replace the use of antibiotics in mastitis control. These efforts, supported by the revolutionary development of nanotechnology, stem cell assays, molecular biological tools, and genomics, enabled the development of new approaches for mastitis-treatment and control. The present review discusses recent concepts to control mastitis such as breeding of mastitis-resistant dairy cows, the development of novel diagnostic and therapeutic tools, the application of communication technology as an educational and epidemiological tool, application of modern mastitis vaccines, cow drying protocols, teat disinfection, housing, and nutrition. These include the application of nanotechnology, stem cell technology, photodynamic and laser therapy or the use of traditional herbal medical plants, nutraceuticals, antibacterial peptides, bacteriocins, antibodies therapy, bacteriophages, phage lysins, and probiotics as alternatives to antibiotics.

SNPs with intermediate minor allele frequencies facilitate accurate breed assignment of Indian Tharparkar cattle

Tharparkar cattle breed is widely known for its superior milch quality and hardiness attributes. This study aimed to develop an ultra-low density breed-specific single nucleotide polymorphism (SNP) genotype panel to accurately quantify Tharparkar populations in biological samples. In this study, we selected and genotyped 72 Tharparkar animals randomly from Cattle & Buffalo Farm of IVRI, India. This Bovine SNP50 BeadChip genotypic datum was merged with the online data from six indigenous cattle breeds and five taurine breeds. Here, we used a combination of pre-selection statistics and the MAF-LD method developed in our laboratory to analyze the genotypic data obtained from 317 individuals of 12 distinct breeds to identify breed-informative SNPs for the selection of Tharparkar cattle. This methodology identified 63 unique Tharparkar-specific SNPs near intermediate gene frequencies. We report several informative SNPs in genes/QTL regions affecting phenotypes or production traits that might differentiate the Tharparkar breed.

Identification of SNPs Associated with Somatic Cell Score in Candidate Genes in Italian Holstein Friesian Bulls

Mastitis is an infectious disease affecting the mammary gland, leading to inflammatory reactions and to heavy economic losses due to milk production decrease. One possible way to tackle the antimicrobial resistance issue stemming from antimicrobial therapy is to select animals with a genetic resistance to this disease. Therefore, aim of this study was to analyze the genetic variability of the SNPs found in candidate genes related to mastitis resistance in Holstein Friesian bulls. Target regions were amplified, sequenced by Next-Generation Sequencing technology on the Illumina® MiSeq, and then analyzed to find correlation with mastitis related phenotypes in 95 Italian Holstein bulls chosen with the aid of a selective genotyping approach. On a total of 557 detected mutations, 61 showed different genotype distribution in the tails of the deregressed EBVs for SCS and 15 were identified as significantly associated with the phenotype using two different approaches. The significant SNPs were identified in intergenic or intronic regions of six genes, known to be key components in the immune system (namely CXCR1, DCK, NOD2, MBL2, MBL1 and M-SAA3.2). These SNPs could be considered as candidates for a future genetic selection for mastitis resistance, although further studies are required to assess their presence in other dairy cattle breeds and their possible negative correlation with other traits.

Discovery of Genomic Characteristics and Selection Signatures in Southern Chinese Local Cattle

Chinese local cattle with a high level of genetic diversity mainly originate from two subspecies; the cattle in northern China are primarily Bos Taurus, and the cattle in southern China are primarily Bos indicus. Cattle from southern China are characterized by a specific phenotype and adapted to the local environment. This study explored the genetic diversity, degree of admixture, and selection signature in eight local cattle breeds in southern China. The lowest level of heterozygosity was found in Hainan and Nandan cattle from Hainan and Guangxi province, respectively, whereas the highest level of heterozygosity was detected in Zhaotong cattle from Yunnan province. A neighbor-joining phylogenetic tree analysis clearly separated Lufeng cattle from other breeds, whereas Leiqiong and Hainan cattle have some crossover. Based on linkage disequilibrium-filtered single nucleotide polymorphisms (SNPs), the admixture analysis revealed two clusters corresponding to the taurine and indicine cattle lineages, and the local cattle breeds from southern China showed a certain degree of admixture. When K = 4 and 9, we found a slight separation among Leiqiong, Lufeng, and Hainan cattle. Meanwhile, we performed a selection signature analysis in Hainan, Leiqiong, and Lufeng cattle distributed in the extreme south of China, using the integrated haplotype score (iHS), Rsb statistic, and BayeScan software. Using the iHS approach, we identified 251, 270, and 256 candidate regions in Lufeng, Leiqiong, and Hainan cattle, respectively. Moreover, we identified 184, 174, and 146 candidate regions in pairwise comparisons of Leiqiong vs. Lufeng, Leiqiong vs. Hainan, and Hainan vs. Lufeng cattle using the Rsb approach. In addition, we identified 76 loci with a total of 48 genes under selection, based on the F_ST approach. Several candidate genes under selection were found to be related to meat quality, immunity, and adaptation to the local environment in southern China. Our results provide significant information about the genetic differences among the cattle breeds from southern China and the possible cause of difference in breed-specific characteristics. Selection signature analysis identified a few candidate SNPs and genes related to certain important traits of these cattle. In general, our results provide valuable insights into the genetic basis of specific traits under selection in certain local cattle breeds.

Genome-Wide DNA Methylation Analysis of Mammary Gland Tissues From Chinese Holstein Cows With Staphylococcus aureus Induced Mastitis

Staphylococcus aureus intramammary infection is one of the most common causes of chronic mastitis in dairy cows, whose development may be associated with epigenetic changes in the expression of important host defense genes. This study aimed to construct a genome-wide DNA methylation profile of the mammary gland of Chinese Holstein cows (n = 3) following experimentally induced S. aureus mastitis, and to explore the potential gene regulatory mechanisms affected by DNA methylation during S. aureus mastitis. DNA was extracted from S. aureus-positive (n = 3) and S. aureus-negative (n = 3) mammary gland quarters and subjected to methylation-dependent restriction-site associated DNA sequencing (Methyl-RAD Seq). Results showed that CmCGG/CmCWGG DNA methylation sites were unevenly distributed and concentrated on chromosomes 5, 11, and 19, and within intergenic regions and intron regions of genes. Compared with healthy control quarters, 9,181 significantly differentially methylated (DM) CmCGG sites and 1,790 DM CmCWGG sites were found in the S. aureus-positive quarters (P < 0.05, |log2FC| > 1). Furthermore, 363 CmCGG differently methylated genes (DMGs) and 301 CmCWGG DMGs (adjusted P < 0.05, |log2FC| > 1) were identified. Gene ontology and KEGG enrichment analysis indicated that CmCGG DMGs are involved in immune response pathways, while the CmCWGG DMGs were mainly enriched in gene ontology terms related to metabolism. The mRNAs of 526 differentially methylated CmCGG genes and 124 differentially methylated CmCWGG genes were also significantly differentially expressed (RNA-Seq data) in the same samples, herein denoted differentially methylated and expressed genes (DMEGs) (P < 0.05). Functional enrichment analysis of DMEGs revealed roles related to biological processes, especially the regulation of immune response to diseases. CmCGG DMEGs like IL6R, TNF, BTK, IL1R2, and TNFSF8 enriched in several immune-related GO terms and pathways indicated their important roles in host immune response and their potential as candidate genes for S. aureus mastitis. These results suggest potential regulatory roles for DNA methylation in bovine mammary gland processes during S. aureus mastitis and serves as a reference for future epigenetic regulation and mechanistic studies.

Identification of novel molecular markers of mastitis caused by Staphylococcus aureus using gene expression profiling in two consecutive generations of Chinese Holstein dairy cattle

Background

Mastitis in dairy cows caused by Staphylococcus aureus is a major problem hindering economic growth in dairy farms worldwide. It is difficult to prevent or eliminate due to its asymptomatic nature and long persistence of infection. Although transcriptomic responses of bovine mammary gland cells to pathogens that cause mastitis have been studied, the common responses of peripheral blood leukocytes to S. aureus infection across two consecutive generations of dairy cattle have not been investigated.

Methods

In the current study, RNA-Seq was used to profile the transcriptomes of peripheral blood leukocytes sampled from S. aureus-infected mothers and their S. aureus-infected daughters, and also healthy non-infected mothers and their healthy daughters. Differential gene expression was evaluated as follows: 1) S. aureus-infected cows versus healthy non-infected cows (S vs. H, which include all the mothers and daughters), 2) S. aureus-infected mothers versus healthy non-infected mothers (SM vs. HM), and 3) S. aureus-infected daughters versus healthy non-infected daughters (SMD vs. HMD).

Results

Analysis of all identified expressed genes in the four groups (SM, SMD, HM, and HMD) showed that EPOR, IL9, IFNL3, CCL26, IL26 were exclusively expressed in both the HM and HMD groups, and that they were significantly (P <  0.05) enriched for the cytokine-cytokine receptor interaction pathway. A total of 17, 13 and 10 differentially expressed genes (DEGs) (FDR P_adj. < 0.1 and |FC| > 1.2) were detected in the three comparisons, respectively. DEGs with P <  0.05 and |FC| > 2 were used for functional enrichment analyses. For the S vs. H comparison, DEGs detected included CCL20, IL13 and MMP3, which are associated with the IL-17 signaling pathway. In the SM vs. HM and SMD vs. HMD comparisons, five (BLA-DQB, C1R, C2, FCGR1A, and KRT10) and six (BLA-DQB, C3AR1, CFI, FCAR, FCGR3A, and LOC10498484) genes, respectively, were involved in the S. aureus infection pathway.

Conclusions

Our study provides insights into the transcriptomic responses of bovine peripheral blood leukocytes across two generations of cattle naturally infected with S. aureus. The genes highlighted in this study could serve as expression biomarkers for mastitis and may also contain sequence variation that can be used for genetic improvement of dairy cattle for resilience to mastitis.

Overview of Research Development on the Role of NF-κB Signaling in Mastitis

Mastitis is the inflammation of the mammary gland. Escherichia coli and Staphylococcus aureus are the most common bacteria responsible for mastitis. When mammary epithelial cells are infected by microorganisms, this activates an inflammatory response. The bacterial infection is recognized by innate pattern recognition receptors (PRRs) in the mammary epithelial cells, with the help of Toll-like receptors (TLRs). Upon activation by lipopolysaccharides, a virulent agent of bacteria, the TLRs further trigger nuclear factor-κB (NF-κB) signaling to accelerate its pathogenesis. The NF-κB has an essential role in many biological processes, such as cell survival, immune response, inflammation and development. Therefore, the NF-κB signaling triggered by the TLRs then regulates the transcriptional expression of specific inflammatory mediators to initiate inflammation of the mammary epithelial cells. Thus, any aberrant regulation of NF-κB signaling may lead to many inflammatory diseases, including mastitis. Hence, the inhibiting of NF-κB signaling has potential therapeutic applications in mastitis control strategies. In this review, we highlighted the regulation and function of NF-κB signaling in mastitis. Furthermore, the role of NF-κB signaling for therapeutic purposes in mastitis control has been explored in the current review.

Metabolite Genome-Wide Association Study (mGWAS) and Gene-Metabolite Interaction Network Analysis Reveal Potential Biomarkers for Feed Efficiency in Pigs

Metabolites represent the ultimate response of biological systems, so metabolomics is considered the link between genotypes and phenotypes. Feed efficiency is one of the most important phenotypes in sustainable pig production and is the main breeding goal trait. We utilized metabolic and genomic datasets from a total of 108 pigs from our own previously published studies that involved 59 Duroc and 49 Landrace pigs with data on feed efficiency (residual feed intake (RFI)), genotype (PorcineSNP80 BeadChip) data, and metabolomic data (45 final metabolite datasets derived from LC-MS system). Utilizing these datasets, our main aim was to identify genetic variants (single-nucleotide polymorphisms (SNPs)) that affect 45 different metabolite concentrations in plasma collected at the start and end of the performance testing of pigs categorized as high or low in their feed efficiency (based on RFI values). Genome-wide significant genetic variants could be then used as potential genetic or biomarkers in breeding programs for feed efficiency. The other objective was to reveal the biochemical mechanisms underlying genetic variation for pigs' feed efficiency. In order to achieve these objectives, we firstly conducted a metabolite genome-wide association study (mGWAS) based on mixed linear models and found 152 genome-wide significant SNPs (p-value < 1.06 × 10-6) in association with 17 metabolites that included 90 significant SNPs annotated to 52 genes. On chromosome one alone, 51 significant SNPs associated with isovalerylcarnitine and propionylcarnitine were found to be in strong linkage disequilibrium (LD). SNPs in strong LD annotated to FBXL4, and CCNC consisted of two haplotype blocks where three SNPs (ALGA0004000, ALGA0004041, and ALGA0004042) were in the intron regions of FBXL4 and CCNC. The interaction network revealed that CCNC and FBXL4 were linked by the hub gene N6AMT1 that was associated with isovalerylcarnitine and propionylcarnitine. Moreover, three metabolites (i.e., isovalerylcarnitine, propionylcarnitine, and pyruvic acid) were clustered in one group based on the low-high RFI pigs. This study performed a comprehensive metabolite-based genome-wide association study (GWAS) analysis for pigs with differences in feed efficiency and provided significant metabolites for which there is significant genetic variation as well as biological interaction networks. The identified metabolite genetic variants, genes, and networks in high versus low feed efficient pigs could be considered as potential genetic or biomarkers for feed efficiency.

Genome Wide Assessment of Genetic Variation and Population Distinctiveness of the Pig Family in South Africa

Genetic diversity is of great importance and a prerequisite for genetic improvement and conservation programs in pigs and other livestock populations. The present study provides a genome wide analysis of the genetic variability and population structure of pig populations from different production systems in South Africa relative to global populations. A total of 234 pigs sampled in South Africa and consisting of village (n = 91), commercial (n = 60), indigenous (n = 40), Asian (n = 5) and wild (n = 38) populations were genotyped using Porcine SNP60K BeadChip. In addition, 389 genotypes representing village and commercial pigs from America, Europe, and Asia were accessed from a previous study and used to compare population clustering and relationships of South African pigs with global populations. Moderate heterozygosity levels, ranging from 0.204 for Warthogs to 0.371 for village pigs sampled from Capricorn municipality in Eastern Cape province of South Africa were observed. Principal Component Analysis of the South African pigs resulted in four distinct clusters of (i) Duroc; (ii) Vietnamese; (iii) Bush pig and Warthog and (iv) a cluster with the rest of the commercial (SA Large White and Landrace), village, Wild Boar and indigenous breeds of Koelbroek and Windsnyer. The clustering demonstrated alignment with genetic similarities, geographic location and production systems. The PCA with the global populations also resulted in four clusters that where populated with (i) all the village populations, wild boars, SA indigenous and the large white and landraces; (ii) Durocs (iii) Chinese and Vietnamese pigs and (iv) Warthog and Bush pig. K = 10 (The number of population units) was the most probable ADMIXTURE based clustering, which grouped animals according to their populations with the exception of the village pigs that showed presence of admixture. AMOVA reported 19.92%-98.62% of the genetic variation to be within populations. Sub structuring was observed between South African commercial populations as well as between Indigenous and commercial breeds. Population pairwise F ST analysis showed genetic differentiation (P ≤ 0.05) between the village, commercial and wild populations. A per marker per population pairwise F ST analysis revealed SNPs associated with QTLs for traits such as meat quality, cytoskeletal and muscle development, glucose metabolism processes and growth factors between both domestic populations as well as between wild and domestic breeds. Overall, the study provided a baseline understanding of porcine diversity and an important foundation for porcine genomics of South African populations.

Short communication: Time-dependent genetic parameters and single-step genome-wide association analyses for predicted milk fatty acid composition in Ayrshire and Jersey dairy cattle

Milk fat content and fatty acid (FA) composition have great economic value to the dairy industry as they are directly associated with taste and chemical-physical characteristics of milk and dairy products. In addition, consumers' choices are not only based on the nutritional aspects of food, but also on products known to promote better health. Milk FA composition is also related to the metabolic status and physiological stages of cows and thus can also be used as indicator for other novel traits of interest (e.g., metabolic diseases and methane yield). Genetic selection is a promising alternative to manipulate milk FA composition. In this study, we aimed to (1) estimate time-dependent genetic parameters for 5 milk FA groups (i.e., short-chain, medium-chain, long-chain, saturated, and unsaturated) predicted based on milk mid-infrared spectroscopy, for Canadian Ayrshire and Jersey breeds, and (2) conduct a time-dependent, single-step genome-wide association study to identify genomic regions, candidate genes, and metabolic pathways associated with milk FA. We analyzed 31,709 test-day records of 9,648 Ayrshire cows from 268 herds, and 34,341 records of 11,479 Jersey cows from 883 herds. The genomic database contained a total of 2,330 Ayrshire and 1,019 Jersey animals. The average daily heritability ranged from 0.18 (long-chain FA) to 0.34 (medium-chain FA) in Ayrshire, and from 0.25 (long-chain and unsaturated FA) to 0.52 (medium-chain and saturated FA) in Jersey. Important genomic regions were identified in Bos taurus autosomes BTA3, BTA5, BTA12, BTA13, BTA14, BTA16, BTA18, BTA20, and BTA21. The proportion of the variance explained by 20 adjacent SNP ranged from 0.71% (saturated FA) to 1.11% (long-chain FA) in Ayrshire, and from 0.70% (unsaturated FA) to 3.09% (medium-chain FA) in Jersey cattle. Important candidate genes and pathways were also identified, such as the PTK2 and TRAPPC9 genes, associated with milk fat percentage, and HMGCS, FGF10, and C6 genes, associated with fertility traits and immune response. Our findings on the genetic parameters and candidate genes contribute to a better understanding of the genetic architecture of milk FA composition in Ayrshire and Jersey dairy cattle.

Copy Number Variations and Expression Levels of Guanylate-Binding Protein 6 Gene Associated with Growth Traits of Chinese Cattle

Association studies have indicated profound effects of copy number variations (CNVs) on various phenotypes in different species. In this study, we identified the CNV distributions and expression levels of guanylate-binding protein 6 (GBP6) associated with the growth traits of Chinese cattle. The results showed that the phenotypic values of body size and weight of Xianan (XN) cattle were higher than those of Nanyang (NY) cattle. The medium CNV types were mostly identified in the XN and NY breeds, but their CNV distributions were significantly different (adjusted p < 0.05). The association analysis revealed that the body weight, cannon circumference and chest circumference of XN cattle had significantly different values in different CNV types (p < 0.05), with CNV gain types (Log22-ΔΔCt > 0.5) displaying superior phenotypic values. We also found that transcription levels varied in different tissues (p < 0.001) and the CNV gain types showed the highest relative gene expression levels in the muscle tissue, consistent with the highest phenotypic values of body weight and cannon circumference among the three CNV types. Consequently, our results suggested that CNV gain types of GBP6 could be used as the candidate markers in the cattle-breeding program for growth traits.

Significant genetic effects of JAK2 and DGAT1 mutations on milk fat content and mastitis resistance in Holsteins

Improving the production traits and resistance against mastitis in dairy cattle is a challenge for animal scientists across the globe. The present study was designed to investigate the genetic effects of single nucleotide polymorphisms (SNPs) in Janus kinase 2 (JAK2) and diacylglycerol acyltransferase (DGAT1) genes with production and mastitis-related traits. Four SNPs in JAK2 and one in DGAT1 were analyzed through Chinese Cow's SNPs Chip-I (CCSC-I) and genotyped in a population of 312 Chinese Holsteins. Our findings demonstrated that milk fat percentage, somatic cell count (SCC), somatic cell score (SCS), serum cytokines interleukin 6 (IL-6) and interferon gamma (IFN-γ) showed significant associations (P < 0.05) with at least one or more identified SNPs. Consequently, the analysis based on haplotypes amongst the SNPs in JAK2 revealed noteworthy (P < 0.05) association with SCC and IL-6. Collectively, our results verified the pleiotropic ability of detected SNPs in bovine JAK2 and DGAT1 for milk fat percentage as well as mastitis-related traits. The significant SNPs in both the genes could serve as powerful genetic markers to minimize mastitis risk. In addition, besides SCC and SCS, the IFN-γ and IL-6 could also be used as indicators of improved genetic resistance against mastitis.

Single-step genome-wide association for longitudinal traits of Canadian Ayrshire, Holstein, and Jersey dairy cattle

Estimating single nucleotide polymorphism (SNP) effects over time is essential to identify and validate candidate genes (or quantitative trait loci) associated with time-dependent variation of economically important traits and to better understand the underlying mechanisms of lactation biology. Therefore, in this study, we aimed to estimate time-dependent effects of SNP and identifying candidate genes associated with milk (MY), fat (FY), and protein (PY) yields, and somatic cell score (SCS) in the first 3 lactations of Canadian Ayrshire, Holstein, and Jersey breeds, as well as suggest their potential pattern of phenotypic effect over time. Random regression coefficients for the additive direct genetic effect were estimated for each animal using single-step genomic BLUP, based on 2 random regression models: one considering MY, FY, and PY in the first 3 lactations and the other considering SCS in the first 3 lactations. Thereafter, SNP solutions were obtained for random regression coefficients, which were used to estimate the SNP effects over time (from 5 to 305 d in lactation). The top 1% of SNP that showed a high magnitude of SNP effect in at least 1 d in lactation were selected as relevant SNP for further analyses of candidate genes, and clustered according to the trajectory of their SNP effects over time. The majority of SNP selected for MY, FY, and PY increased the magnitude of their effects over time, for all breeds. In contrast, for SCS, most selected SNP decreased the magnitude of their effects over time, especially for the Holstein and Jersey breeds. In general, we identified a different set of candidate genes for each breed, and similar genes were found across different lactations for the same trait in the same breed. For some of the candidate genes, the suggested pattern of phenotypic effect changed among lactations. Among the lactations, candidate genes (and their suggested phenotypic effect over time) identified for the second and third lactations were more similar to each other than for the first lactation. Well-known candidate genes with major effects on milk production traits presented different suggested patterns of phenotypic effect across breeds, traits, and lactations in which they were identified. The candidate genes identified in this study can be used as target genes in studies of gene expression.

On the genomic regions associated with milk lactose in Fleckvieh cattle

Lactose is a sugar uniquely found in mammals' milk and it is the major milk solid in bovines. Lactose yield (LY, kg/d) is responsible for milk volume, whereas lactose percentage (LP) is thought to be more related to epithelial integrity and thus to udder health. There is a paucity of studies that have investigated lactose at the genomic level in dairy cows. This paper aimed to improve our knowledge on LP and LY, providing new insights into the significant genomic regions affecting these traits. A genome-wide association study for LP and LY was carried out in Fleckvieh cattle by using bulls' deregressed estimated breeding values of first lactation as pseudo-phenotypes. Heritabilities of first-lactation test-day LP and LY estimated using linear animal models were 0.38 and 0.25, respectively. A total of 2,854 bulls genotyped with a 54K SNP chip were available for the genome-wide association study; a linear mixed model approach was adopted for the analysis. The significant SNP of LP were scattered across the whole genome, with signals on chromosomes 1, 2, 3, 7, 12, 16, 18, 19, 20, 28, and 29; the top 4 significant SNP explained 4.90% of the LP genetic variance. The signals were mostly in regions or genes with involvement in molecular intra- or extracellular transport; for example, CDH5, RASGEF1C, ABCA6, and SLC35F3. A significant region within chromosome 20 was previously shown to affect mastitis or somatic cell score in cattle. As regards LY, the significant SNP were concentrated in fewer regions (chromosomes 6 and 14), related to mastitis/somatic cell score, immune response, and transport mechanisms. The 5 most significant SNP for LY explained 8.45% of genetic variance and more than one-quarter of this value has to be attributed to the variant within ADGRB1. Significant peaks in target regions remained even after adjustment for the 2 most significant variants previously detected on BTA6 and BTA14. The present study is a prelude for deeper investigations into the biological role of lactose for milk secretion and volume determination, stressing the connection with genes regulating intra- or extracellular trafficking and immune and inflammatory responses in dairy cows. Also, these results improve the knowledge on the relationship between lactose and udder health; they support the idea that LP and its derived traits are potential candidates as indicators of udder health in breeding programs aimed to enhance cows' resistance to mastitis.

Identification of genetic markers associated with milk production traits in Chinese Holstein cattle based on post genome-wide association studies

With the rapid development of dairy industry, the breeding process of dairy cows has been accelerated. In previous genome-wide association studies (GWAS), a large number of genetic markers have been reported which may contribute to the selection of Holstein populations with superior milk-producing traits, but they remain to be further verified before practical application. In this study, 90 single nucleotide polymorphisms (SNPs) were selected, which were reported to be significantly associated with five milk production traits, including 305-day milk yield (305MY), 305-day milk fat percent (305FC), 305-day milk protein percent (305PC), 305-day milk fat yield (305FY) and 305-day milk protein yield (305PY). Effective 305-day data and fresh DNA samples were obtained from 295 healthy cows with gestational age of 1-4. Matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) was used to perform precise genotyping of these loci, followed by site association and haplotype analysis. Results showed that 36 out of 90 loci were supported to be used as genetic markers. In particular, several novel and effective haplotypes were also presented. Overall, our results verified tens of useful markers and provided a basis for further development of breeding strategies.

Three Novel Players: PTK2B, SYK, and TNFRSF21 Were Identified to Be Involved in the Regulation of Bovine Mastitis Susceptibility via GWAS and Post-transcriptional Analysis

Bovine mastitis is a common inflammatory disease caused by multiple factors in early lactation or dry period. Genome wide association studies (GWAS) can provide a convenient and effective strategy for understanding the biological basis of mastitis and better prevention. 2b-RADseq is a high-throughput sequencing technique that offers a powerful method for genome-wide genetic marker development and genotyping. In this study, single nucleotide polymorphisms (SNPs) of the immune-regulated gene correlative with mastitis were screened and identified by two stage association analysis via GWAS-2b-RADseq in Chinese Holstein cows. We have screened 10,058 high quality SNPs from 7,957,920 tags and calculated their allele frequencies. Twenty-seven significant SNPs were co-labeled in two GWAS analysis models [Bayesian (P < 0.001) and Logistic regression (P < 0.01)], and only three SNPs (rs75762330, C > T, PIC = 0.2999; rs88640083, A > G, PIC = 0.1676; rs20438858, G > A, PIC = 0.3366) were annotated to immune-regulated genes (PTK2B, SYK, and TNFRSF21). Identified three SNPs are located in non-coding regions with low or moderate genetic polymorphisms. However, independent sample population validation (Case-control study) data showed that three important SNPs (rs75762330, P < 0.025, OR > 1; rs88640083, P < 0.005, OR > 1; rs20438858, P < 0.001, OR < 1) were significantly associated with clinical mastitis trait. Importantly, PTK2B and SYK expression was down-regulated in both peripheral blood leukocytes (PBLs) of clinical mastitis cows and in vitro LPS (E. coli)-stimulated bovine mammary epithelial cells, while TNFRSF21 was up-regulated. Under the same conditions, expression of Toll-like receptor 4 (TLR4), AKT1, and pro-inflammatory factors (IL-1β and IL-8) were also up-regulated. Interestingly, network analysis indicated that PTK2B and SYK are co-expressed in innate immune signaling pathway of Chinese Holstein. Taken together, these results provided strong evidence for the study of SNPs in bovine mastitis, and revealed the role of SYK, PTK2B, and TNFRSF21 in bovine mastitis susceptibility/tolerance.

Identification of QTLs for low somatic cell count in Murrah buffaloes

Mastitis, the most frequent and costly disease in buffalo, is the major cause of morbidity. The somatic cell count, an indirect indicator of susceptibility/resistance to mastitis, is a low heritable trait and thus a perfect candidate for marker assisted selection. Half sib families (12) were created and the somatic cell count was recorded at 3 stages of lactation during the first lactation of the 2,422 daughters belonging to 12 sires. Partial genome scan was carried out using interval mapping with different algorithms. The QTLs obtained for each half sib family were further subjected to meta analysis to identify chromosomal regions associated with somatic cell count on 8 chromosomes of buffalo. Four metaQTL regions were identified on chromosomes BBU1q, BBU8, and BBU10; 3 metaQTL regions on BBU2q, BBU9 and BBU15; 2 metaQTL regions on BBU6 and 1 on BBU7 of buffalo. Comparative genomics was used for finding out genes underlying the metaQTL regions; 1,065 genes were underlying the metaQTL regions in buffaloes assuming buffalo–cattle–human synteny. Genes (78) mapped to immune response. These genes are supposedly important candidate genes for further analysis. Gene ontology and network analysis was carried out on these genes. The genes identified belonged to immune response and defense mechanism. The QTL markers identified in the present analysis can be used in the breeding programs of buffalo to select the bulls, which are less susceptible to mastitis.

Genome-wide association analysis for β-hydroxybutyrate concentration in Milk in Holstein dairy cattle

BACKGROUND

Ketosis in dairy cattle has been shown to cause a high morbidity in the farm and substantial financial losses to dairy farmers. Ketosis symptoms, however, are difficult to identify, therefore, the amount of ketone bodies (mainly β-hydroxybutyric acid, BHB) is used as an indicator of subclinical ketosis in cows. It has also been shown that milk BHB concentrations have a strong correlation with ketosis in dairy cattle. Mid-infrared spectroscopy (MIR) has recently became a fast, cheap and high-throughput method for analyzing milk components. The aim of this study was to perform a genome-wide association study (GWAS) on the MIR-predicted milk BHB to identify genomic regions, genes and pathways potentially affecting subclinical ketosis in North American Holstein dairy cattle.

RESULTS

Several significant regions were identified associated with MIR-predicted milk BHB concentrations (indicator of subclinical ketosis) in the first lactation (SCK1) and second and later lactations (SCK2) in Holstein dairy cows. The strongest association was located on BTA6 for SCK1 and BTA14 on SCK2. Several SNPs on BTA6 were identified in regions and variants reported previously to be associated with susceptibility to ketosis and clinical mastitis in Jersey and Holstein dairy cattle, respectively. One highly significant SNP on BTA14 was found within the DGAT1 gene with known functions on fat metabolism and inflammatory response in dairy cattle. A region on BTA6 and three SNPs on BTA20 were found to overlap between SCK1 and SCK2. However, a novel region on BTA20 (55-63 Mb) for SCK2 was also identified, which was not reported in previous association studies. Enrichment analysis of the list of candidate genes within the identified regions for MIR-predicted milk BHB concentrations yielded molecular functions and biological processes that may be involved in the inflammatory response and lipid metabolism in dairy cattle.

CONCLUSIONS

The results of this study confirmed several SNPs and genes identified in previous studies as associated with ketosis susceptibility and immune response, and also found a novel region that can be used for further analysis to identify causal variations and key regulatory genes that affect clinical/ subclinical ketosis.

Identifying Loci Under Positive Selection in Yellow Korean Cattle (Hanwoo)

Jeju Black cattle is one of the aboriginal Korean cattle breeds that has been isolated in Jeju island for a long time, while Yellow Hanwoo cattle has been extensively selected for beef production traits for the last several decades. Aiming to investigate broader patterns of selection, we genotyped 352 Yellow Hanwoo and 169 Jeju Black cattle using a customized 150 K bovine chip. Our composite selection signals’ analysis to identify selection signatures (cross-population extended haplotype homozygosity [XP-EHH], ΔSAF, and F_ST) identified recent and strong signature of selection near many loci with mutations affecting the traits under strong selection as outlier in Yellow Hanwoo, including SCP2 (P = 8.41 × 10⁻¹⁰) that may be involved in the meat quality. We found nine candidate regions with significant clusters of selection signals, and further bioinformatics analyses of the genes located within these regions revealed mainly genes involved in G-protein coupled receptor signaling pathway (GO:0007186) or olfactory transduction (bta04740), which may be due to adaptation to natural environments in Jeju island. Based on the stronger correlation of Ne₁₀/Ne₁₀₀ ratio between Yellow Hanwoo (0.61) and Jeju Black (0.66) cattle, our results suggest that the difference of chromosomal regions of selection signature between the 2 cattle breeds was due to a consequence of selection processes to adapt to environmental differences between Jeju island and the main inland, Korean peninsula.

A Meta-Analysis Including Pre-selected Sequence Variants Associated With Seven Traits in Three French Dairy Cattle Populations

A within-breed genome-wide association study (GWAS) is useful when identifying the QTL that segregates in a breed. However, an across-breed meta-analysis can be used to increase the power of identification and precise localization of QTL that segregate in multiple breeds. Precise localization will allow including QTL information from other breeds in genomic prediction due to the persistence of the linkage phase between the causal variant and the marker. This study aimed to identify and confirm QTL detected in within-breed GWAS through a meta-analysis in three French dairy cattle breeds. A set of sequence variants selected based on their functional annotations were imputed into 50 k genotypes for 46,732 Holstein, 20,096 Montbeliarde, and 11,944 Normande cows to identify QTL for milk production, the success rate at insemination of cows (fertility) and stature. We conducted within-breed GWAS followed by across-breed meta-analysis using a weighted Z-scores model on the GWAS summary data (i.e., P-values, effect direction, and sample size). After Bonferroni correction, the GWAS result identified 21,956 significantly associated SNP (P _FWER < 0.05), while meta-analysis result identified 9,604 significant SNP (P _FWER < 0.05) associated with the phenotypes. The meta-analysis identified 36 QTL for milk yield, 48 QTL for fat yield and percentage, 29 QTL for protein yield and percentage, 13 QTL for fertility, and 16 QTL for stature. Some of these QTL were not significant in the within-breed GWAS. Some previously identified causal variants were confirmed, e.g., BTA14:1802265 (fat percentage, P = 1.5 × 10^-760; protein percentage, P = 7.61 × 10^-348) both mapping the DGAT1-K232A mutation and BTA14:25006125 (P = 8.58 × 10^-140) mapping PLAG1 gene was confirmed for stature in Montbeliarde. New QTL lead SNP shared between breeds included the intronic variant rs109205829 (NFIB gene), and the intergenic variant rs41592357 (1.38 Mb upstream of the CNTN6 gene and 0.65 Mb downstream of the CNTN4 gene). Rs110425867 (ZFAT gene) was the top variant associated with fertility, and new QTL lead SNP included rs109483390 (0.1 Mb upstream of the TNFAIP3 gene and 0.07 Mb downstream of PERP gene), and rs42412333 (0.45 Mb downstream of the RPL10L gene). An across-breed meta-analysis had greater power to detect QTL as opposed to a within breed GWAS. The QTL detected here can be incorporated in routine genomic predictions.

Prioritizing candidate genes post-GWAS using multiple sources of data for mastitis resistance in dairy cattle

Background

Improving resistance to mastitis, one of the costliest diseases in dairy production, has become an important objective in dairy cattle breeding. However, mastitis resistance is influenced by many genes involved in multiple processes, including the response to infection, inflammation, and post-infection healing. Low genetic heritability, environmental variations, and farm management differences further complicate the identification of links between genetic variants and mastitis resistance. Consequently, studies of the genetics of variation in mastitis resistance in dairy cattle lack agreement about the responsible genes.

Results

We associated 15,552,968 imputed whole-genome sequencing markers for 5147 Nordic Holstein cattle with mastitis resistance in a genome-wide association study (GWAS). Next, we augmented P-values for markers in genes in the associated regions using Gene Ontology terms, Kyoto Encyclopedia of Genes and Genomes pathway analysis, and mammalian phenotype database. To confirm results of gene-based analyses, we used gene expression data from E. coli-challenged cow udders. We identified 22 independent quantitative trait loci (QTL) that collectively explained 14% of the variance in breeding values for resistance to clinical mastitis (CM). Using association test statistics with multiple pieces of independent information on gene function and differential expression during bacterial infection, we suggested putative causal genes with biological relevance for 12 QTL affecting resistance to CM in dairy cattle.

Conclusion

Combining information on the nearest positional genes, gene-based analyses, and differential gene expression data from RNA-seq, we identified putative causal genes (candidate genes with biological evidence) in QTL for mastitis resistance in Nordic Holstein cattle. The same strategy can be applied for other traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5050-x) contains supplementary material, which is available to authorized users.

Genome-wide microRNA profiling of bovine milk-derived exosomes infected with Staphylococcus aureus

Bovine milk is rich in exosomes, which contain abundant miRNAs and play important roles in the regulation of neonatal growth and development of adaptive immunity. Here, we analyzed miRNA expression profiles of bovine milk exosomes from three healthy and three mastitic cows, and then six miRNA libraries were constructed. Interestingly, we detected no scRNAs and few snRNAs in milk exosomes; this result indicated a potential preference for RNA packaging in milk exosomes. A total of 492 known and 980 novel exosomal miRNAs were detected, and the 10 most expressed miRNAs in the six samples accounted for 80-90% of total miRNA-associated reads. Expression analyses identified 18 miRNAs with significantly different expression between healthy and infected animals; the predicted target genes of differentially expressed miRNAs were significantly enriched in immune system process, response to stimulus, growth, etc. Moreover, target genes were significantly enriched in several Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways including inflammatory, immune, and cancer pathways. Our survey provided comprehensive information about milk exosomes and exosomal miRNAs involved in mastitis. Moreover, the differentially expressed miRNAs, especially miR-223 and miR-142-5p, could be considered as potential candidates for mastitis.

Assessment of runs of homozygosity islands and estimates of genomic inbreeding in Gyr (Bos indicus) dairy cattle

Background

Runs of homozygosity (ROH) are continuous homozygous segments of the DNA sequence. They have been applied to quantify individual autozygosity and used as a potential inbreeding measure in livestock species. The aim of the present study was (i) to investigate genome-wide autozygosity to identify and characterize ROH patterns in Gyr dairy cattle genome; (ii) identify ROH islands for gene content and enrichment in segments shared by more than 50% of the samples, and (iii) compare estimates of molecular inbreeding calculated from ROH (F_ROH), genomic relationship matrix approach (F_GRM) and based on the observed versus expected number of homozygous genotypes (F_HOM), and from pedigree-based coefficient (F_PED).

Results

ROH were identified in all animals, with an average number of 55.12 ± 10.37 segments and a mean length of 3.17 Mb. Short segments (ROH_{1–2 Mb}) were abundant through the genomes, which accounted for 60% of all segments identified, even though the proportion of the genome covered by them was relatively small. The findings obtained in this study suggest that on average 7.01% (175.28 Mb) of the genome of this population is autozygous. Overlapping ROH were evident across the genomes and 14 regions were identified with ROH frequencies exceeding 50% of the whole population. Genes associated with lactation (TRAPPC9), milk yield and composition (IRS2 and ANG), and heat adaptation (HSF1, HSPB1, and HSPE1), were identified. Inbreeding coefficients were estimated through the application of F_ROH, F_GRM, F_HOM, and F_PED approaches. F_PED estimates ranged from 0.00 to 0.327 and F_ROH from 0.001 to 0.201. Low to moderate correlations were observed between F_PED-F_ROH and F_GRM-F_ROH, with values ranging from −0.11 to 0.51. Low to high correlations were observed between F_ROH-F_HOM and moderate between F_PED-F_HOM and F_GRM-F_HOM. Correlations between F_ROH from different lengths and F_PED gradually increased with ROH length.

Conclusions

Genes inside ROH islands suggest a strong selection for dairy traits and enrichment for Gyr cattle environmental adaptation. Furthermore, low F_PED-F_ROH correlations for small segments indicate that F_PED estimates are not the most suitable method to capture ancient inbreeding. The existence of a moderate correlation between larger ROH indicates that F_ROH can be used as an alternative to inbreeding estimates in the absence of pedigree records.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4365-3) contains supplementary material, which is available to authorized users.