Genetic factors influencing milk and fat yields in tropically adapted dairy cattle: insights from quantitative trait loci analysis and gene associations

OBJECTIVE

The objective of this study was to identify genes associated with 305-day milk yield (MY) and fat yield (FY) that also influence the adaptability of the Thai multibreed dairy cattle population to tropical conditions.

METHODS

A total of 75,776 imputed and actual single nucleotide polymorphisms (SNPs) from 2,661 animals were used to identify genomic regions associated with MY and FY using the single-step genomic best linear unbiased predictions. Fixed effects included herd-yearseason, breed regression, heterosis regression and calving age regression effects. Random effects were animal additive genetic and residual. Individual SNPs with a p-value smaller than 0.05 were selected for gene mapping, function analysis, and quantitative trait loci (QTL) annotation analysis.

RESULTS

A substantial number of QTLs associated with MY (9,334) and FY (8,977) were identified by integrating SNP genotypes and QTL annotations. Notably, we discovered 17 annotated QTLs within the health and exterior QTL classes, corresponding to nine unique genes. Among these genes, Rho GTPase activating protein 15 (ARHGAP15) and catenin alpha 2 (CTNNA2) have previously been linked to physiological traits associated with tropical adaptation in various cattle breeds. Interestingly, these two genes also showed signs of positive selection, indicating their potential role in conferring tolerance to trypanosomiasis, a prevalent tropical disease.

CONCLUSION

Our findings provide valuable insights into the genetic basis of MY and FY in the Thai multibreed dairy cattle population, shedding light on the underlying mechanisms of tropical adaptation. The identified genes represent promising targets for future breeding strategies aimed at improving milk and fat production while ensuring resilience to tropical challenges. This study significantly contributes to our understanding of the genetic factors influencing milk production and adaptability in dairy cattle, facilitating the development of sustainable genetic selection strategies and breeding programs in tropical environments.

Genetic parameters and genome-wide association studies for mozzarella and milk production traits, lactation length, and lactation persistency in Murrah buffaloes

Genetic and genomic analyses of longitudinal traits related to milk production efficiency are paramount for optimizing water buffaloes breeding schemes. Therefore, this study aimed to (1) compare single-trait random regression models under a single-step genomic BLUP setting based on alternative covariance functions (i.e., Wood, Wilmink, and Ali and Schaeffer) to describe milk (MY), fat (FY), protein (PY), and mozzarella (MZY) yields, fat-to-protein ratio (FPR), somatic cell score (SCS), lactation length (LL), and lactation persistency (LP) in Murrah dairy buffaloes (Bubalus bubalis); (2) combine the best functions for each trait under a multiple-trait framework; (3) estimate time-dependent SNP effects for all the studied longitudinal traits; and (4) identify the most likely candidate genes associated with the traits. A total of 323,140 test-day records from the first lactation of 4,588 Murrah buffaloes were made available for the study. The model included the average curve of the population nested within herd-year-season of calving, systematic effects of number of milkings per day, and age at first calving as linear and quadratic covariates, and additive genetic, permanent environment, and residual as random effects. The Wood model had the best goodness of fit based on the deviance information criterion and posterior model probabilities for all traits. Moderate heritabilities were estimated over time for most traits (0.30 ± 0.02 for MY; 0.26 ± 0.03 for FY; 0.45 ± 0.04 for PY; 0.28 ± 0.05 for MZY; 0.13 ± 0.02 for FPR; and 0.15 ± 0.03 for SCS). The heritability estimates for LP ranged from 0.38 ± 0.02 to 0.65 ± 0.03 depending on the trait definition used. Similarly, heritabilities estimated for LL ranged from 0.10 ± 0.01 to 0.14 ± 0.03. The genetic correlation estimates across days in milk (DIM) for all traits ranged from -0.06 (186-215 DIM for MY-SCS) to 0.78 (66-95 DIM for PY-MZY). The SNP effects calculated for the random regression model coefficients were used to estimate the SNP effects throughout the lactation curve (from 5 to 305 d). Numerous relevant genomic regions and candidate genes were identified for all traits, confirming their polygenic nature. The candidate genes identified contribute to a better understanding of the genetic background of milk-related traits in Murrah buffaloes and reinforce the value of incorporating genomic information in their breeding programs.

Characterization of runs of Homozygosity revealed genomic inbreeding and patterns of selection in indigenous sahiwal cattle

Runs of homozygosity (ROH) are contiguous genomic regions, homozygous across all sites which arise in an individual due to the parents transmitting identical haplotypes to their offspring. The genetic improvement program of Sahiwal cattle after decades of selection needs re-assessment of breeding strategy and population phenomena. Hence, the present study was carried out to optimize input parameters in PLINK for ROH estimates, to explore ROH islands and assessment of pedigree and genome-based inbreeding in Sahiwal cattle. The sliding window approach with parameters standardized to define ROH for the specific population under study was used for the identification of runs. The optimum maximum gap, density, window-snp and window-threshold were 250 Kb, 120 Kb/SNP, 10, 0.05 respectively and ROH patterns were also characterized. ROH islands were defined as the short homozygous genomic regions shared by a large proportion of individuals in a population, containing significantly higher occurrences of ROH than the population specific threshold level. These were identified using the -homozyg-group function of the PLINK v1.9 program. Our results indicated that the Islands of ROH harbor a few candidate genes, ACAD11, RFX4, BANP, UBA5 that are associated with major economic traits. The average FPED (Pedigree based inbreeding coefficient), FROH (Genomic inbreeding coefficient), FHOM (Inbreeding estimated as the ratio of observed and expected homozygous genotypes), FGRM (Inbreeding estimated on genomic relationship method) and FGRM0.5 (Inbreeding estimated from the diagonal of a GRM with allele frequencies near to 0.5) were 0.009, 0.091, 0.035, -0.104 and -0.009, respectively. Our study revealed the optimum parameter setting in PLINK viz. maximal gaps between two SNPs, minimal density of SNPs in a segment (in kb/SNP) and scanning window size to identify ROH segments, which will enable ROH estimation more efficient and comparable across various SNP genotyping-based studies. The result further emphasized the significant role of genomics in unraveling population diversity, selection signatures and inbreeding in the ongoing Sahiwal breed improvement programs.

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.

Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.

Signatures of selection in riverine buffalo populations revealed by genome-wide SNP data

The detection of selection signatures assists in understanding domestication, evolution, and the identification of genomic regions related to adaptation and production traits in buffaloes. The emergence of high-throughput technologies like Next Generation Sequencing and SNP genotyping had expanded our ability to detect these signatures of selection. In this study, we sought to identify signatures of selection in five buffalo populations (Brazilian Murrah, Bulgarian Murrah, Indian Murrah, Nili-Ravi, and Kundi) using Axiom Buffalo 90 K Genotyping Array data. Using seven different methodologies (Tajima's D, CLR, ROH, iHS, FST, FLK and hapFLK), we identified selection signatures in 374 genomic regions, spanning a total of 381 genes and 350 quantitative trait loci (QTLs). Among these, several candidate genes were associated with QTLs for milk production, reproduction, growth and carcass traits. The genes and QTLs reported in this study provide insight into selection signals shaping the genome of buffalo breeds. Our findings can aid in further genomic association studies, genomic prediction, and the implementation of breeding programmes in Indian buffaloes.

Methylation Levels in the Promoter Region of FHIT and PIAS1 Genes Associated with Mastitis Resistance in Xinjiang Brown Cattle

Mastitis causes serious economic losses in the dairy industry, but there are no effective treatments or preventive measures. In this study, the ZRANB3, PIAS1, ACTR3, LPCAT2, MGAT5, and SLC37A2 genes in Xinjiang brown cattle, which are associated with mastitis resistance, were identified using a GWAS. Pyrosequencing analysis showed that the promoter methylation levels of the FHIT and PIAS1 genes in the mastitis group were higher and lower, respectively, than those in the healthy group (65.97 ± 19.82% and 58.00 ± 23.52%). However, the methylation level of the PIAS1 gene promoter region in the mastitis group was lower than that in the healthy group (11.48 ± 4.12% and 12.17 ± 4.25%). Meanwhile, the methylation levels of CpG3, CpG5, CpG8, and CpG15 in the promoter region of the FHIT and PIAS1 genes in the mastitis group were significantly higher than those in the healthy group (p < 0.01), respectively. RT-qPCR showed that the expression levels of the FHIT and PIAS1 genes were significantly higher in the healthy group than those in the mastitis group (p < 0.01). Correlation analysis showed that the promoter methylation level of the FHIT gene was negatively correlated with its expression. Hence, increased methylation in the promoter of the FHIT gene reduces the mastitis resistance in Xinjiang brown cattle. Finally, this study provides a reference for the molecular-marker-assisted selection of mastitis resistance in dairy cattle.

The association of gene polymorphisms with milk production and mastitis resistance phenotypic traits in dairy cattle

The aim of this study was to evaluate the association between gene polymorphisms (SNPs) and mastitis indicators and their relationship with milk production profitability in dairy herd.A functional analysis was also performed of five genes containing the studied SNPs and those located close by. DNA was isolated from the hair bulb of 320 dairy cows kept in three herds and SNP-microarray analysis was performed. The data on 299 cows was subjected to final statistical analysis using AI-REML method with one-trait repeatability test-day animal model and pedigree information using the DMU4 package. Five from 35 SNPs significantly associated with mastitis indicators or production traits and located within a gene or no more than 500,000 nucleotides from the gene were selected for the functional and economic analysis. A questionnaire was also developed to collect associated economic data of 219 cows from three herds, such as the value of milk production and direct costs incurred over three years; this allowed the gross margin, direct profitability index and direct costs incurred to produce one liter of milk to be determined, among others. None of the five studied SNPs were related to protein content. The rs110785912(T/A), found near CXCR4, and rs136813430(T/C), located in the TLR4 gene exon, were associated with lnSCC, while rs110455063(C/G), located near IGFI, was associated with milk yield, fat and total solid contents. rs109421300(T/C), associated with fat/protein content ratio, as well as fat and total solid content, is located in the DGAT1 gene intron. rs41587003(A/C), located in the DLG2 gene intron, was associated with lactose content. The economic analysis revealed differences between the variants of the three tested SNPs. The T/C variant of the rs136813430(T/C) SNP was characterized by the highest gross margin, the highest direct profitability index and the lowest costs incurred to produce 1 liter of milk. The T/A variant of rs110785912(T/A) was related to low lnSCC and was characterized by the highest direct profitability index. In turn, the C/C variant of the rs41587003(T/C) was related to the lowest level of lactose and the highest costs of milk production. It appears that rs136813430(T/C) may be the most promising of the tested SNPs for increasing the profitability of milk production. To our knowledge, it is the first effort to assess directly a correlation between the DNA polymorphism and economic output of a dairy enterprise.

Differential gene expression in peripheral leukocytes of pre-weaned Holstein heifer calves with respiratory disease

Bovine respiratory disease (BRD) is a leading cause of calf morbidity and mortality, and prevalence remains high despite current management practices. Differential gene expression (DGE) provides detailed insight into individual immune responses and can illuminate enriched pathways and biomarkers that contribute to disease susceptibility and outcomes. The aims of this study were to investigate differences in peripheral leukocyte gene expression in Holstein preweaned heifer calves 1) with and without BRD, and 2) across weeks of age. Calves were enrolled for this short-term longitudinal study on two commercial dairies in Washington State. Calves were assessed every two weeks throughout the pre-weaning period using clinical respiratory scoring (CRS) and thoracic ultrasonography (TUS), and blood samples were collected. Calves were selected that were either healthy (n = 10) or had BRD diagnosed by CRS (n = 7), TUS (n = 6), or both (n = 6) in weeks 5 or 7 of life). Three consecutive time point samples were analyzed for each BRD calf consisting of PRE, ONSET, and POST samples. Nineteen genes of interest were selected based on previous gene expression studies in cattle: ALOX15, BPI, CATHL6, CXCL8, DHX58, GZMB, HPGD, IFNG, IL17D, IL1R2, ISG15, LCN2, LIF, MX1, OAS2, PGLYRP1, S100A8, SELP, and TNF. Comparisons were made between age and disease time point matched BRD and healthy calves as well as between calf weeks of age. No DGE was observed between diseased and healthy calves; however, DGE was observed between calf weeks of age regardless of disease state. Developmental differences in leukocyte gene expression, phenotype, and functionality make pre-weaned calves immunologically distinct from mature cattle, and early life shifts in calf leukocyte populations likely contribute to the age-related gene expression differences we observed. Age overshadows disease impacts to influence gene expression in young calves, and immune development progresses upon a common trajectory regardless of disease during the preweaning period.

Identification of Key Genes and Potential Pathways Associated with Mastitis Induced by E. coli

Mastitis is one of the most important diseases of dairy cattle. It is an infectious disease leading to an inflammatory reaction in the cow's mammary gland. Escherichia coli is one of the common bacteria which induce mastitis in cows. The aim of this study was to identify key genes and potential pathways associated with mastitis induced by E. coli in dairy cattle using bioinformatics analysis. The gene expression profile of ten samples including five adjacent tissues from a quarter infected with Escherichia coli and five tissues from a healthy quarter of dairy cattle was assessed using GEO2R. Gene ontology and pathway analysis were performed using bioinformatics tools. A total of 156 differentially expressed genes were detected which 95 genes were upregulated and 61 genes were downregulated in adjacent tissue of quarter infected compared with healthy tissue. Cellular oxidant detoxification and oxidation-reduction process were the most significant biological process terms in gene ontology analysis. The most important pathways of DEGs were the biosynthesis of amino acids, p53 signaling pathway, and Metabolic pathways. Three important modules were identified and their path enrichment analysis was performed. There are 10 core genes, among which SOD2, COL1A2, COL3A1, POSTN, ALDH18A1, and CBS may be the main genes associated with mastitis, which can be considered as candidate genes in the prevention and carly diagnosis program of mastitis.

Proportion of Concentrate in the Diet of Early Lactation Dairy Cows Has Contrasting Effects on Circulating Leukocyte Global Transcriptomic Profiles, Health and Fertility According to Parity

The functionality of circulating leukocytes in dairy cows is suppressed after calving, with negative energy balance as a risk factor. Leukocyte transcriptomic profiles were compared separately in 44 multiparous (MP) and 18 primiparous (PP) Holstein-Friesian cows receiving diets differing in concentrate proportion to test whether immune dysfunction could be mitigated by appropriate nutrition. After calving, cows were offered either (1) low concentrate (LC); (2) medium concentrate (MC) or (3) high concentrate (HC) diets with proportions of concentrate to grass silage of 30%:70%, 50%:50% and 70%:30%, respectively. Cow phenotype data collected included circulating metabolites, milk yield and health and fertility records. RNA sequencing of circulating leukocytes at 14 days in milk was performed. The HC diet improved energy balance in both age groups. There were more differentially expressed genes in PP than MP cows (460 vs. 173, HC vs. LC comparison) with few overlaps. The MP cows on the LC diet showed upregulation of the complement and coagulation cascade and innate immune defence mechanisms against pathogens and had a trend of more cases of mastitis and poorer fertility. In contrast, the PP cows on the HC diet showed greater immune responses based on both gene expression and phenotypic data and longer interval of calving to conception. The leukocytes of MP and PP cows therefore responded differentially to the diets between age, nutrient supply and immunity affecting their health and subsequent fertility.

Utilization of 2D and 3D cell cultures for the modelling of intramammary infection in sheep

Ovine mastitis is defined as the inflammation of the sheep udder, most commonly caused in response to intramammary infections. Based on the occurrence of clinical signs, mastitis is characterized as either clinical or subclinical (SCM). The impact of ovine SCM on the overall sustainability of dairy sheep farms has been substantially documented underpinning the significance of efficient diagnosis. Although SCM can be detected in cows, the performance and the validity of the methods used do not transfer in dairy sheep. This fact challenges the development of evidence-based ovine udder health management protocols and renders the detection and control of ovine mastitis rather problematic. Currently, cell culture-based models are being successfully used in biomedical studies and have also been effectively used in the case of bovine mastitis. The objective of the present study was to culture ovine primary mammary cells for the development of 2D and 3D cell culture-based models for the study of ovine mammary gland and to focus on the first stages of the intramammary infection by common mastitis-inducing pathogens. Cells were infected by E. coli and S. aureus mimicking the first stages of natural intramammary infections. The secreted proteins were subjected to mass-spectrometry resulting in the identification of a total of 79 distinct proteins. Among those, several had already been identified in healthy or mastitic milk, while others had not been previously detected for in the ovine mammary secretome. Our results suggest that the development of cell-based models for studying specific stages of intramammary infection has the potential to be beneficial for the udder health management in dairy sheep.

Methylome and transcriptome data integration reveals potential roles of DNA methylation and candidate biomarkers of cow Streptococcus uberis subclinical mastitis

Background

Mastitis caused by different pathogens including Streptococcus uberis (S. uberis) is responsible for huge economic losses to the dairy industry. In order to investigate the potential genetic and epigenetic regulatory mechanisms of subclinical mastitis due to S. uberis, the DNA methylome (whole genome DNA methylation sequencing) and transcriptome (RNA sequencing) of milk somatic cells from cows with naturally occurring S. uberis subclinical mastitis and healthy control cows (n = 3/group) were studied.

Results

Globally, the DNA methylation levels of CpG sites were low in the promoters and first exons but high in inner exons and introns. The DNA methylation levels at the promoter, first exon and first intron regions were negatively correlated with the expression level of genes at a whole-genome-wide scale. In general, DNA methylation level was lower in S. uberis-positive group (SUG) than in the control group (CTG). A total of 174,342 differentially methylated cytosines (DMCs) (FDR < 0.05) were identified between SUG and CTG, including 132,237, 7412 and 34,693 DMCs in the context of CpG, CHG and CHH (H = A or T or C), respectively. Besides, 101,612 methylation haplotype blocks (MHBs) were identified, including 451 MHBs that were significantly different (dMHB) between the two groups. A total of 2130 differentially expressed (DE) genes (1378 with up-regulated and 752 with down-regulated expression) were found in SUG. Integration of methylome and transcriptome data with MethGET program revealed 1623 genes with significant changes in their methylation levels and/or gene expression changes (MetGDE genes, MethGET P-value < 0.001). Functional enrichment of genes harboring ≥ 15 DMCs, DE genes and MetGDE genes suggest significant involvement of DNA methylation changes in the regulation of the host immune response to S. uberis infection, especially cytokine activities. Furthermore, discriminant correlation analysis with DIABLO method identified 26 candidate biomarkers, including 6 DE genes, 15 CpG-DMCs and 5 dMHBs that discriminated between SUG and CTG.

Conclusion

The integration of methylome and transcriptome of milk somatic cells suggests the possible involvement of DNA methylation changes in the regulation of the host immune response to subclinical mastitis due to S. uberis. The presented genetic and epigenetic biomarkers could contribute to the design of management strategies of subclinical mastitis and breeding for mastitis resistance.

LncRNA GAS5 Suppresses Colorectal Cancer Progress by Target miR-21/LIFR Axis

GAS5 is abnormally high in colorectal cancer tissues, which is a specific expression of lncRNA in colorectal cancer (CRC). Nevertheless, its biological function in CRC has not been elucidated. The abnormal high expression of GAS5 in CRC is the specific expression of lncRNA in CRC. The purpose of our study is to explore the effect of GAS5 on CRC and its mechanism. The expression of GAS5 in 53 paired normal and colorectal cancer tissues and colorectal cancer cell lines was detected by real-time PCR. The biological effects of GAS5, miR-21, and LIFR were measured by functional assays, including wound healing, transwell assays, and in vivo assays. We ensured the carcinogenesis role of GAS5 in CRC in the xenograft nude model. The dual-luciferase reporter assay system and chromatin immunoprecipitation method were used for target evaluation and Western blot for verification. GAS5 was significantly decreased in tumor tissues and CRC cells, and the low expression of CAS5 in CRC promoted tumor metastasis and decreased the survival of patients. GAS5 knockdown increases the cell viability, inhibits apoptosis, and promotes migration. Xenografted tumors in nude mice studies showed that GAS5 knockdown promoted tumor growth and caused worse lesions in colorectal. Furthermore, GAS5 increases the expression level of target gene LIFR to promote the apoptosis of CRC cells by binding to miR-21. Our study revealed that a novel pathway about lncRNA GAS5 inhibited the proliferation and metastasis of CRC cells by targeting miR-21/LIFR which provides a new strategy to treat CRC.

Genome-wide elucidation of CNV regions and their association with production and reproduction traits in composite Vrindavani cattle

The present study was aimed to analyze the genome-wide copy number variations (CNVs) in Vrindavani composite cattle and concatenate them into CNV regions (CNVRs), and finally test the association of CNVRs with different production and reproduction traits. Genotypic data, generated on BovineSNP50 Beadchip (v3) array for 96 Vrindavani animals, was used to elucidate the CNVs at the genome level. Intensity data covering over 53,218 SNP genotypes on bovine genome was used. Algorithm based on Hidden Markov Model was employed in PennCNV program to detect, normalize and filter CNVs across the genome. 252 putative CNVs, detected via PennCNV program, in different individuals were concatenated into 71 CNV regions (CNVRs) using CNVRuler program. Association of CNVRs with important (re)production traits in Vrindavani animals was assessed using linear regression. Five CNVRs were found to be significantly associated with ten important (re)production traits. The genes harbored in these regions provided useful insights into the association of CNVRs with genes and ultimately the variation at phenotype level. Important genes that overlapped with CNVRs included WASHC4, HS6ST3, MBNL2, TOLLIP, PIDD1 and TSPAN4. Furthermore, the CNVRs were found to overlap with important QTLs available in AnimalQTL database which affect milk yield and composition along with reproduction and immune function traits. The copy number states of three enes were validated using digital droplet PCR technique. The results from the present study significantly enhance the understanding about CNVs in Vrindavani cattle and should help establish its CNV map. The study will also enable further investigation on association of these variants with important traits of economic interest including disease incidence.

Identifying biomarkers for evaluating wound extent and age in the contused muscle of rats using microarray analysis: a pilot study

Wound age estimation is still one of the most important and significant challenges in forensic practice. The extent of wound damage greatly affects the accuracy and reliability of wound age estimation, so it is important to find effective biomarkers to help diagnose wound degree and wound age. In the present study, the gene expression profiles of both mild and severe injuries in 33 rats were assayed at 0, 1, 3, 24, 48, and 168 hours using the Affymetrix microarray system to provide biomarkers for the evaluation of wound age and the extent of the wound. After obtaining thousands of differentially expressed genes, a principal component analysis, the least absolute shrinkage and selection operator, and a time-series analysis were used to select the most predictive prognostic genes. Finally, 15 genes were screened for evaluating the extent of wound damage, and the top 60 genes were also screened for wound age estimation in mild and severe injury. Selected indicators showed good diagnostic performance for identifying the extent of the wound and wound age in a Fisher discriminant analysis. A function analysis showed that the candidate genes were mainly related to cell proliferation and the inflammatory response, primarily IL-17 and the Hematopoietic cell lineage signalling pathway. The results revealed that these genes play an essential role in wound-healing and yield helpful and valuable potential biomarkers for further targeted studies.

Detection of candidate genes affecting milk production traits in sheep using whole-genome sequencing analysis

BACKGROUND

Artificial and natural selection for important economic traits and genetic adaptation of the populations to specific environments have led to the changes on the sheep genome. Recent advances in genome sequencing methods have made it possible to use comparative genomics tools to identify genes under selection for traits of economic interest in domestic animals.

OBJECTIVES

In this study, we compared the genomes of Assaf and Awassi sheep breeds with those of the Cambridge, Romanov and British du cher sheep breeds to explore positive selection signatures for milk traits using nucleotide diversity (Pi) and FST statistical methods.

METHODS

Genome sequences from fourteen sheep with a mean sequence depth of 9.32X per sample were analysed, and a total of 23 million single nucleotide polymorphisms (SNPs) were called and applied for this study. Genomic clustering of breeds was identified using ADMIXTURE software. The FST and Pi values for each SNP were computed between population A (Assaf and Awassi) and population B (Cambridge, British du cher, and Romanov).

RESULTS

The results of the PCA grouped two classes for these five dairy sheep breeds. The selection signatures analysis displayed 735 and 515 genes from FST and nucleotide diversity (Pi) statistical methods, respectively. Among all these, 12 genes were shared between the two approaches. The most conspicuous genes were related to milk traits, including ST3GAL1 (the synthesis of oligosacáridos), CSN1S1 (milk protein), CSN2 (milk protein), OSBPL8 (fatty acid traits), SLC35A3 (milk fat and protein percentage), VPS13B (total milk production, fat yield, and protein yield), DPY19L1 (peak yield), CCDC152 (lactation persistency and somatic cell count), NT5DC1 (lactation persistency), P4HTM (test day protein), CYTH4 (FAT Production) and METRNL (somatic cell), U1 (milk traits), U6 (milk traits) and 5S_RRNA (milk traits).

CONCLUSIONS

The findings provide new insight into the genetic basis of sheep milk properties and can play a role in designing sheep breeding programs incorporating genomic information.

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.

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.

Effects of Maternal Nutrition on Female Offspring Weight Gain and Sexual Development

Maternal nutrition during pregnancy influences postnatal life of animals; nevertheless, few studies have investigated its effects on the productive performance and reproductive development of heifers. This study evaluated the performance, reproductive development, and correlation between reproduction × fat thickness and performance × ribeye area (REA) traits of heifers. We also performed an exploratory genomic association during the rearing period in heifers submitted to fetal programming. The study comprised 55 Nellore heifers born to dams exposed to one of the following nutritional planes: control, without protein-energy supplementation; PELT, protein-energy last trimester, protein-energy supplementation offered in the final third of pregnancy; and PEWG, protein-energy whole gestation, protein-energy supplementation upon pregnancy confirmation. Protein-energy supplementation occurred at the level of 0.3% live weight. After weaning, heifers were submitted to periodic evaluations of weight and body composition by ultrasonography. From 12 to 18 months, we evaluated the reproductive tract of heifers to monitor its development for sexual precocity and ovarian follicle population. The treatments had no effect (p > 0.05) on average daily gain; however, the weight of the animals showed a significant difference over time (p = 0.017). No differences were found between treatments for REA, backfat, and rump fat thickness, nor for puberty age, antral follicular count, and other traits related to reproductive tract development (p > 0.05). The correlation analysis between performance traits and REA showed high correlations (r > 0.37) between REA at weaning and year versus weight from weaning until yearling; however, no correlation was found for reproductive development traits versus fat thickness (p > 0.05). The exploratory genomic association study showed one single-nucleotide polymorphism (SNP) for each treatment on an intergenic region for control and PEWG, and the one for PELT on an intronic region of RAPGEF1 gene. Maternal nutrition affected only the weight of the animals throughout the rearing period.

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 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.

Integrative Systems Biology Analysis Elucidates Mastitis Disease Underlying Functional Modules in Dairy Cattle

Background: Mastitis is the most prevalent disease in dairy cattle and one of the most significant bovine pathologies affecting milk production, animal health, and reproduction. In addition, mastitis is the most common, expensive, and contagious infection in the dairy industry. Methods: A meta-analysis of microarray and RNA-seq data was conducted to identify candidate genes and functional modules associated with mastitis disease. The results were then applied to systems biology analysis via weighted gene coexpression network analysis (WGCNA), Gene Ontology, enrichment analysis for the Kyoto Encyclopedia of Genes and Genomes (KEGG), and modeling using machine-learning algorithms. Results: Microarray and RNA-seq datasets were generated for 2,089 and 2,794 meta-genes, respectively. Between microarray and RNA-seq datasets, a total of 360 meta-genes were found that were significantly enriched as "peroxisome," "NOD-like receptor signaling pathway," "IL-17 signaling pathway," and "TNF signaling pathway" KEGG pathways. The turquoise module (n = 214 genes) and the brown module (n = 57 genes) were identified as critical functional modules associated with mastitis through WGCNA. PRDX5, RAB5C, ACTN4, SLC25A16, MAPK6, CD53, NCKAP1L, ARHGEF2, COL9A1, and PTPRC genes were detected as hub genes in identified functional modules. Finally, using attribute weighting and machine-learning methods, hub genes that are sufficiently informative in Escherichia coli mastitis were used to optimize predictive models. The constructed model proposed the optimal approach for the meta-genes and validated several high-ranked genes as biomarkers for E. coli mastitis using the decision tree (DT) method. Conclusion: The candidate genes and pathways proposed in this study may shed new light on the underlying molecular mechanisms of mastitis disease and suggest new approaches for diagnosing and treating E. coli mastitis in dairy cattle.

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.

Global transcriptomic profiles of circulating leucocytes in early lactation cows with clinical or subclinical mastitis

Bovine mastitis, an inflammatory disease of the mammary gland, is classified as subclinical or clinical. Circulating neutrophils are recruited to the udder to combat infection. We compared the transcriptomic profiles in circulating leukocytes between healthy cows and those with naturally occurring subclinical or clinical mastitis. Holstein Friesian dairy cows from six farms in EU countries were recruited. Based on milk somatic cell count and clinical records, cows were classified as healthy (n = 147), subclinically (n = 45) or clinically mastitic (n = 22). Circulating leukocyte RNA was sequenced with Illumina NextSeq single end reads (30 M). Differentially expressed genes (DEGs) between the groups were identified using CLC Genomics Workbench V21, followed by GO enrichment analysis. Both subclinical and clinical mastitis caused significant changes in the leukocyte transcriptome, with more intensive changes attributed to clinical mastitis. We detected 769 DEGs between clinical and healthy groups, 258 DEGs between subclinical and healthy groups and 193 DEGs between clinical and subclinical groups. Most DEGs were associated with cell killing and immune processes. Many upregulated DEGs in clinical mastitis encoded antimicrobial peptides (AZU1, BCL3, CAMP, CATHL1, CATHL2, CATHL4,CATHL5, CATHL6, CCL1, CXCL2, CXCL13, DEFB1, DEFB10, DEFB4A, DEFB7, LCN2, PGLYRP1, PRTN3, PTX3, S100A8, S100A9, S100A12, SLC11A1, TF and LTF) which were not upregulated in subclinical mastitis. The use of transcriptomic profiles has identified a much greater up-regulation of genes encoding antimicrobial peptides in circulating leukocytes of cows with naturally occurring clinical compared with subclinical mastitis. These could play a key role in combatting disease organisms.

Effect of diet and nonesterified fatty acid levels on global transcriptomic profiles in circulating peripheral blood mononuclear cells in early lactation dairy cows

After calving, lipid mobilization caused by increased nutrient demands for lactation leads to elevated circulating concentrations of nonesterified fatty acids (NEFA). Excessive NEFA levels have previously been identified as a major risk factor for postpartum immunosuppression. The aim of this study was to investigate changes in global transcriptomic gene expression of peripheral blood mononuclear cells (PBMC) in dairy cows offered different early lactation diets (high concentrate, n = 7; medium, n = 8; or low, n = 9) and with differing circulating levels of NEFA. Cows were classified as having NEFA concentrations of either <500 µM (low, n = 6), 500 to 750 µM (medium, n = 8) or >750 µM (high, n = 10) at 14 d in milk. Plasma urea concentrations were greater for cows on the high concentrate diet but β-hydroxybutyrate and glucose concentrations did not differ significantly between either dietary treatments or NEFA groups. Cows with high NEFA weighed more at drying off and suffered greater body condition score loss after calving. The PBMC were isolated at 14 d in milk, and RNA was extracted for RNA sequencing. Differential gene expression was analyzed with DESeq2 with q-value for false discovery rate control followed by Gene Ontology Enrichment. Although there were no differentially expressed genes associated with lactation diet, 304 differentially expressed genes were identified between cows with high and low circulating NEFA, with 118 upregulated and 186 downregulated. Gene Ontology enrichment analysis demonstrated that biological adhesion and immune system process were foremost among various PBMC functions which were altered relating to body defenses and immunity. High NEFA concentrations were associated with inhibited cellular adhesion function by downregulating 20 out of 26 genes (by up to 17-fold) related to this process. Medium NEFA concentrations altered a similar set of functions as high NEFA, but with smaller enrichment scores. Localization and immune system process were most significant, with biological adhesion ranking only eleventh. Our results demonstrated that increased circulating NEFA concentrations, but not diet, were associated with immune system processes in PBMC in early lactation cows. Leukocyte cell-to-cell adhesion was inhibited when the NEFA concentration exceeded 750 µM, which would reduce the efficiency of diapedesis and so contribute to decreased body defense mechanisms and predispose animals to infection.

Structural equation modeling for unraveling the multivariate genomic architecture of milk proteins in dairy cattle

The aims of this study were to investigate potential functional relationships among milk protein fractions in dairy cattle and to carry out a structural equation model (SEM) GWAS to provide a decomposition of total SNP effects into direct effects and effects mediated by traits that are upstream in a phenotypic network. To achieve these aims, we first fitted a mixed Bayesian multitrait genomic model to infer the genomic correlations among 6 milk nitrogen fractions [4 caseins (CN), namely κ-, β-, α_S1-, and α_S2-CN, and 2 whey proteins, namely β-lactoglobulin (β-LG) and α-lactalbumin (α-LA)], in a population of 989 Italian Brown Swiss cows. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2 (Illumina Inc.). A Bayesian network approach using the max-min hill-climbing (MMHC) algorithm was implemented to model the dependencies or independence among traits. Strong and negative genomic correlations were found between β-CN and α_S1-CN (-0.706) and between β-CN and κ-CN (-0.735). The application of the MMHC algorithm revealed that κ-CN and β-CN seemed to directly or indirectly influence all other milk protein fractions. By integrating multitrait model GWAS and SEM-GWAS, we identified a total of 127 significant SNP for κ-CN, 89 SNP for β-CN, 30 SNP for α_S1-CN, and 14 SNP for α_S2-CN (mostly shared among CN and located on Bos taurus autosome 6) and 15 SNP for β-LG (mostly located on Bos taurus autosome 11), whereas no SNP passed the significance threshold for α-LA. For the significant SNP, we assessed and quantified the contribution of direct and indirect paths to total marker effect. Pathway analyses confirmed that common regulatory mechanisms (e.g., energy metabolism and hormonal and neural signals) are involved in the control of milk protein synthesis and metabolism. The information acquired might be leveraged for setting up optimal management and selection strategies aimed at improving milk quality and technological characteristics in dairy cattle.

Genomic studies of milk-related traits in water buffalo (Bubalus bubalis) based on single-step genomic best linear unbiased prediction and random regression models

Genomic selection has been widely implemented in many livestock breeding programs, but it remains incipient in buffalo. Therefore, this study aimed to (1) estimate variance components incorporating genomic information in Murrah buffalo; (2) evaluate the performance of genomic prediction for milk-related traits using single- and multitrait random regression models (RRM) and the single-step genomic best linear unbiased prediction approach; and (3) estimate longitudinal SNP effects and candidate genes potentially associated with time-dependent variation in milk, fat, and protein yields, as well as somatic cell score (SCS) in multiple parities. The data used to estimate the genetic parameters consisted of a total of 323,140 test-day records. The average daily heritability estimates were moderate (0.35 ± 0.02 for milk yield, 0.22 ± 0.03 for fat yield, 0.42 ± 0.03 for protein yield, and 0.16 ± 0.03 for SCS). The highest heritability estimates, considering all traits studied, were observed between 20 and 280 d in milk (DIM). The genetic correlation estimates at different DIM among the evaluated traits ranged from -0.10 (156 to 185 DIM for SCS) to 0.61 (36 to 65 DIM for fat yield). In general, direct selection for any of the traits evaluated is expected to result in indirect genetic gains for milk yield, fat yield, and protein yield but also increase SCS at certain lactation stages, which is undesirable. The predicted RRM coefficients were used to derive the genomic estimated breeding values (GEBV) for each time point (from 5 to 305 DIM). In general, the tuning parameters evaluated when constructing the hybrid genomic relationship matrices had a small effect on the GEBV accuracy and a greater effect on the bias estimates. The SNP solutions were back-solved from the GEBV predicted from the Legendre random regression coefficients, which were then used to estimate the longitudinal SNP effects (from 5 to 305 DIM). The daily SNP effect for 3 different lactation stages were performed considering 3 different lactation stages for each trait and parity: from 5 to 70, from 71 to 150, and from 151 to 305 DIM. Important genomic regions related to the analyzed traits and parities that explain more than 0.50% of the total additive genetic variance were selected for further analyses of candidate genes. In general, similar potential candidate genes were found between traits, but our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the traits across parities. These results contribute to a better understanding of the genetic architecture of milk production traits in dairy buffalo and reinforce the relevance of incorporating genomic information to genetically evaluate longitudinal traits in dairy buffalo. Furthermore, the candidate genes identified can be used as target genes in future functional genomics studies.

Genetic Parameters and Genome-Wide Association Studies of Eight Longevity Traits Representing Either Full or Partial Lifespan in Chinese Holsteins

Due to the complexity of longevity trait in dairy cattle, two groups of trait definitions are widely used to measure longevity, either covering the full lifespan or representing only a part of it to achieve an early selection. Usually, only one group of longevity definition is used in breeding program for one population, and genetic studies on the comparisons of two groups of trait definitions are scarce. Based on the data of eight traits well representing the both groups of trait definitions, the current study investigated genetic parameters and genetic architectures of longevity in Holsteins. Heritabilities and correlations of eight longevity traits were estimated using single-trait and multi-trait animal models, with the data from 103,479 cows. Among the cows with phenotypes, 2,630 cows were genotyped with the 150K-SNP panel. A single-trait fixed and random Circuitous Probability Unification model was performed to detect candidate genes for eight longevity traits. Generally, all eight longevity traits had low heritabilities, ranging from 0.038 for total productive life and herd life to 0.090 for days from the first calving to the end of first lactation or culling. High genetic correlations were observed among the traits within the same definition group: from 0.946 to 0.997 for three traits reflecting full lifespan and from 0.666 to 0.997 for five traits reflecting partial productive life. Genetic correlations between two groups of traits ranged from 0.648 to 0.963, and increased gradually with the extension of lactations number regarding the partial productive life traits. A total of 55 SNPs located on 25 chromosomes were found genome-wide significantly associated with longevity, in which 12 SNPs were associated with more than one trait, even across traits of different definition groups. This is the first study to investigate the genetic architecture of longevity representing both full and the partial lifespan simultaneously, which will assist the selection of an appropriate trait definition for genetic improvement of longevity. Because of high genetic correlations with the full lifespan traits and higher heritability, the partial productive life trait measured as the days from the first calving to the end of the third lactation or culling could be a good alternative for early selection on longevity. The candidate genes identified by this study, such as RPRM, GRIA3, GTF2H5, CA5A, CACNA2D1, FGF10, and DNAJA3, could be used to pinpoint causative mutations for longevity and further benefit the genomic improvement of longevity in dairy cattle.

Weighted gene co-expression network analysis identifies modules and functionally enriched pathways in the lactation process

The exponential growth in knowledge has resulted in a better understanding of the lactation process in a wide variety of animals. However, the underlying genetic mechanisms are not yet clearly known. In order to identify the mechanisms involved in the lactation process, various mehods, including meta-analysis, weighted gene co-express network analysis (WGCNA), hub genes identification, gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment at before peak (BP), peak (P), and after peak (AP) stages of the lactation processes have been employed. A total of 104, 85, and 26 differentially expressed genes were identified based on PB vs. P, BP vs. AP, and P vs. AP comparisons, respectively. GO and KEGG pathway enrichment analysis revealed that DEGs were significantly enriched in the "ubiquitin-dependent ERAD" and the "chaperone cofactor-dependent protein refolding" in BP vs. P and P vs. P, respectively. WGCNA identified five significant functional modules related to the lactation process. Moreover, GJA1, AP2A2, and NPAS3 were defined as hub genes in the identified modules, highlighting the importance of their regulatory impacts on the lactation process. The findings of this study provide new insights into the complex regulatory networks of the lactation process at three distinct stages, while suggesting several candidate genes that may be useful for future animal breeding programs. Furthermore, this study supports the notion that in combination with a meta-analysis, the WGCNA represents an opportunity to achieve a higher resolution analysis that can better predict the most important functional genes that might provide a more robust bio-signature for phenotypic traits, thus providing more suitable biomarker candidates for future studies.

Signatures of Selection in Composite Vrindavani Cattle of India

Vrindavani is an Indian composite cattle breed developed by crossbreeding taurine dairy breeds with native indicine cattle. The constituent breeds were selected for higher milk production and adaptation to the tropical climate. However, the selection response for production and adaptation traits in the Vrindavani genome is not explored. In this study, we provide the first overview of the selection signatures in the Vrindavani genome. A total of 96 Vrindavani cattle were genotyped using the BovineSNP50 BeadChip and the SNP genotype data of its constituent breeds were collected from a public database. Within-breed selection signatures in Vrindavani were investigated using the integrated haplotype score (iHS). The Vrindavani breed was also compared to each of its parental breeds to discover between-population signatures of selection using two approaches, cross-population extended haplotype homozygosity (XP-EHH) and fixation index (F_ST). We identified 11 common regions detected by more than one method harboring genes such as LRP1B, TNNI3K, APOB, CACNA2D1, FAM110B, and SPATA17 associated with production and adaptation. Overall, our results suggested stronger selective pressure on regions responsible for adaptation compared to milk yield.

Inbreeding coefficients and runs of homozygosity islands in Brazilian water buffalo

Characterization of autozygosity is relevant to monitor genetic diversity and manage inbreeding levels in breeding programs. Identification of autozygosity hotspots can unravel genomic regions targeted by selection for economically important traits and can help identify candidate genes for selection. In this study, we estimated the inbreeding levels of a Brazilian population of Murrah buffalo undergoing selection for milk production traits, particularly milk yield. We also studied the distribution of runs of homozygosity (ROH) islands and identified putative genes and quantitative trait loci (QTL) under selection. We genotyped 422 Murrah buffalo for 51,611 SNP; 350 of these had ROH longer than 10 Mb, indicating the occurrence of inbreeding in the last 5 generations. The mean length of the ROH per animal was 4.28 ± 1.85 Mb. Inbreeding coefficients were calculated from the genomic relationship matrix, the pedigree, and the ROH, with estimates varying between 0.242 and 0.035. Inbreeding estimates from the pedigree had a low correlation with the genomic estimates, and estimates from the genomic relationship matrix were much higher than those from the pedigree or the ROH. Signatures of selection were identified in 6 genomic regions, located on chromosomes 1, 2, 3, 5, 16, and 18, encompassing a total of 190 genes and 174 QTL. Many of the genes (e.g., APRT and ACSF3) and QTL identified are related to milk production traits, such as milk yield, milk fat yield and percentage, and milk protein yield and percentage. Other genes are associated with reproduction and immune response traits as well as morphological aspects of the buffalo species. Inbreeding levels in this population are still low but are increasing due to selection and should be managed to avoid future losses due to inbreeding depression. The proximity of genes linked to milk production traits with genes associated with reproduction and immune system traits suggests the need to include these latter genes in the breeding program to avoid negatively affecting them due to selection for production traits.

A Fixed Cohort Field Study of Gene Expression in Circulating Leukocytes From Dairy Cows With and Without Mastitis

Specifically designed gene expression studies can be used to prioritize candidate genes and identify novel biomarkers affecting resilience against mastitis and other diseases in dairy cattle. The primary goal of this study was to assess whether specific peripheral leukocyte genes expressed differentially in a previous study of dairy cattle with postpartum disease, also would be expressed differentially in peripheral leukocytes from a diverse set of different dairy cattle with moderate to severe clinical mastitis. Four genes were selected for this study due to their differential expression in a previous transcriptomic analysis of circulating leukocytes from dairy cows with and without evidence of early postpartum disease. An additional 15 genes were included based on their cellular, immunologic, and inflammatory functions associated with resistance and tolerance to mastitis. This fixed cohort study was conducted on a conventional dairy in Washington state. Cows >50 days in milk (DIM) with mastitis (n = 12) were enrolled along with healthy cows (n = 8) selected to match the DIM and lactation numbers of mastitic cows. Blood was collected for a complete blood count (CBC), serum biochemistry, leukocyte isolation, and RNA extraction on the day of enrollment and twice more at 6 to 8-days intervals. Latent class analysis was performed to discriminate healthy vs. mastitic cows and to describe disease resolution. RNA samples were processed by the Primate Diagnostic Services Laboratory (University of Washington, Seattle, WA). Gene expression analysis was performed using the Nanostring System (Nanostring Technologies, Seattle, Washington, USA). Of the four genes (C5AR1, CATHL6, LCN2, and PGLYRP1) with evidence of upregulation in cows with mastitis, three of those genes (CATHL6, LCN2, and PGLYRP1) were investigated due to their previously identified association with postpartum disease. These genes are responsible for immunomodulatory molecules that selectively enhance or alter host innate immune defense mechanisms and modulate pathogen-induced inflammatory responses. Although further research is warranted to explain their functional mechanisms and bioactivity in cattle, our findings suggest that these conserved elements of innate immunity have the potential to bridge disease states and target tissues in diverse dairy populations.

Identification of markers associated with estimated breeding value and horn colour in Hungarian Grey cattle

OBJECTIVE

This study was conducted to estimate effect of single nucleotide polymorphisms (SNP) on the estimated breeding value of Hungarian Grey (HG) bulls and to find markers associated with horn colour.

METHODS

Genotypes 136 HG animals were determined on Geneseek high-density Bovine SNP 150K BeadChip. A multi-locus mixed-model was applied for statistical analyses.

RESULTS

Six SNPs were identified to be associated (-log10P>10) with green and white horn. These loci are located on chromosome 1, 3, 9, 18, and 25. Seven loci (on chromosome 1, 3, 6, 9, 10, 28) showed considerable association (-log10P>10) with the estimated breeding value.

CONCLUSION

Analysis provides markers for further research of horn colour and supplies markers to achieve more effective selection work regarding estimated breeding value of HG.

Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle

Highlighting genomic profiles for geographically distinct subpopulations of the same breed may provide insights into adaptation mechanisms to different environments, reveal genomic regions divergently selected, and offer initial guidance to joint genomic analysis. Here, we characterized similarities and differences between the genomic patterns of Angus subpopulations, born and raised in Canada (N = 382) and Brazil (N = 566). Furthermore, we systematically scanned for selection signatures based on the detection of autozygosity islands common between the two subpopulations, and signals of divergent selection, via FST and varLD tests. The principal component analysis revealed a sub-structure with a close connection between the two subpopulations. The averages of genomic relationships, inbreeding coefficients, and linkage disequilibrium at varying genomic distances were rather similar across them, suggesting non-accentuated differences in overall genomic diversity. Autozygosity islands revealed selection signatures common to both subpopulations at chromosomes 13 (63.77-65.25 Mb) and 14 (22.81-23.57 Mb), which are notably known regions affecting growth traits. Nevertheless, further autozygosity islands along with FST and varLD tests unravel particular sites with accentuated population subdivision at BTAs 7 and 18 overlapping with known QTL and candidate genes of reproductive performance, thermoregulation, and resistance to infectious diseases. Our findings indicate overall genomic similarity between Angus subpopulations, with noticeable signals of divergent selection in genomic regions associated with the adaptation in different environments.

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.

Genome-wide scan of selection signatures in Dehong humped cattle for heat tolerance and disease resistance

Dehong humped cattle (DHH) is an indigenous zebu breed from southwestern China that possesses characteristics of heat tolerance and strong disease resistance and adapts well to the local tropical and subtropical climatic conditions. However, information on selection signatures of DHH is scarce. Herein, we compared the genomes of DHH and each of Diqing and Zhaotong cattle breeds using the population differentiation index (F_ST ), cross-population extended haplotype homozygosity (XP-EHH) and cross-population composite likelihood ratio (XP-CLR) methods to explore the genomic signatures of heat tolerance and disease resistance in DHH. Several pathways and genes carried selection signatures, including thermal sweating (calcium signaling pathway), heat shock (HSF1) and oxidative stress response (PLCB1, PLCB4), coat color (RAB31), feed intake (ATP8A1, SHC3) and reproduction (TP63, MAP3K13, PTPN4, PPP3CC, ADAMTSL1, SS18L1, OSBPL2, TOX, RREB1, GRK2). These identified pathways and genes may contribute to heat tolerance in DHH. Simultaneously, we also identified LIPH, TP63 and CBFA2T3 genes under positive selection that were associated with immunity.

Transcriptomic analysis on the promoter regions discover gene networks involving mastitis in cattle

Mastitis is one of the costliest diseases in dairy farms caused by infection of different microorganisms such as Escherichia coli, Streptococcus uberis and Staphylococcus aureus. Promoters are significantly involved in regulating gene expression and shedding light on the mechanisms of transcriptional regulation in physiological and immunological processes of the infections. Exploiting regulatory elements such as transcription factor binding sites (TFBSs modules) on the promoter region could reveal co-regulated genes, which allow screating regulatory models and executing a cross-sectional analysis on several databases. In this study, the promoter regions of 11 genes associated with contagious mastitis including CCL4, CXCL8, STAT3, IKBKB, MAPK14, NFKBIA, NFKB1, TNF, IL18, IL6, and HCK were investigated to predict the activating regulatory modules on promoters and to discover the key related transcription factors. By exploring the promoter regions, 228 genes were discovered comprising the same transcription factors modules. Out of 228 genes, 36 were validated using five microarray datasets. The promoter research of these genes revealed that as many as 7 down-regulated and 12 up-regulated genes are predictable in the network. The genes whose functions were associated with the initial gene list (11 genes), were identified by DAVID queries with TFBSs models implying that the approach provides a clear image of the underlying regulatory mechanism of gene expression profile and offers a novel approach in designing gene networks in cattle.

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.

Eph/ephrin signalling serves a bidirectional role in lipopolysaccharide‑induced intestinal injury

A growing body of evidence has demonstrated that Eph/ephrin signalling may serve a central role in intestinal diseases. However, whether erythropoietin-producing hepatocellular (Eph)/ephrin signalling is associated with the development of post-infectious irritable bowel syndrome (PI-IBS) is still unknown. In the present study, the role of Eph/Ephrin signalling in lipopolysaccharide (LPS)-induced intestinal injury was evaluated in vivo and in vitro. LPS treatment significantly increased the levels of proinflammatory mediators [monocyte chemoattractant protein-1, tumour necrosis factor α, interleukin (IL)-1β, IL-6, intercellular adhesion molecule 1 and vascular cell adhesion molecule-1], activated the EphA2-Ephrin A1, protein kinase B (Akt)-nuclear factor (NF)-κB, Src-NF-κB and Wnt/β-catenin signalling pathways, and inhibited EphB1-Ephrin B3 signalling in colon tissues, and primary cultured enteric neuronal and glial cells. Notably, EphA2 monoclonal antibody (mAb) treatment or Ephrin B3 overexpression could partially alleviate the LPS-induced upregulation of proinflammatory mediators, and Akt-NF-κB, Src-NF-κB and Wnt/β-catenin signalling pathways. In addition, EphA2 mAb treatment could partially inhibit LPS-induced inactivation of EphB-Ephrin B3 signalling, while Ephrin B3 overexpression could abrogate LPS-induced activation of EphA2-Ephrin A1 signalling. EphB1/Ephrin B3 signalling may antagonise the EphA2/Ephrin A1-dependent pathway following LPS treatment. The results associated with the EphA2 signaling pathway, indicated that Eph/ephrin signalling may serve a bidirectional role in LPS-induced intestinal injury. Eph/ephrin signalling may be a novel therapeutic target for LPS-induced intestinal injury and potentially PI-IBS.

Polymorphisms in the selectin gene cluster are associated with fertility and survival time in a population of Holstein Friesian cows

Selectins are adhesion molecules, which mediate attachment between leucocytes and endothelium. They aid extravasation of leucocytes from blood into inflamed tissue during the mammary gland’s response to infection. Selectins are also involved in attachment of the conceptus to the endometrium and subsequent placental development. Poor fertility and udder health are major causes for culling dairy cows. The three identified bovine selectin genes SELP, SELL and SELE are located in a gene cluster. SELP is the most polymorphic of these genes. Several SNP in SELP and SELE are associated with human vascular disease, while SELP SNP rs6127 has been associated with recurrent pregnancy loss in women. This study describes the results of a gene association study for SNP in SELP (n = 5), SELL (n = 2) and SELE (n = 1) with fertility, milk production and longevity traits in a population of 337 Holstein Friesian dairy cows. Blood samples for PCR-RFLP were collected at 6 months of age and animals were monitored until either culling or 2,340 days from birth. Three SNP in SELP_Ex4-6 formed a haplotype block containing a Glu/Ala substitution at rs42312260. This region was associated with poor fertility and reduced survival times. SELP_Ex8 (rs378218397) coded for a Val475Met variant locus in the linking region between consensus repeats 4 and 5, which may influence glycosylation. The synonymous SNP rs110045112 in SELE_Ex14 deviated from Hardy Weinberg equilibrium. For both this SNP and rs378218397 there were too few AA homozygotes present in the population and AG heterozygotes had significantly worse fertility than GG homozygotes. Small changes in milk production associated with some SNP could not account for the reduced fertility and only SELP_Ex6 showed any association with somatic cell count. These results suggest that polymorphisms in SELP and SELE are associated with the likelihood of successful pregnancy, potentially through compromised implantation and placental development.

Determination of Genetic Structure and Signatures of Selection in Three Strains of Tanzania Shorthorn Zebu, Boran and Friesian Cattle by Genome-Wide SNP Analyses

Background

More than 90 percent of cattle in Tanzania belong to the indigenous Tanzania Short Horn Zebu (TSZ) population which has been classified into 12 strains based on historical evidence, morphological characteristics, and geographic distribution. However, specific genetic information of each TSZ population has been lacking and has caused difficulties in designing programs such as selection, crossbreeding, breed improvement or conservation. This study was designed to evaluate the genetic structure, assess genetic relationships, and to identify signatures of selection among cattle of Tanzania with the main goal of understanding genetic relationship, variation and uniqueness among them.

Methodology/Principal findings

The Illumina Bos indicus SNP 80K BeadChip was used to genotype genome wide SNPs in 168 DNA samples obtained from three strains of TSZ cattle namely Maasai, Tarime and Sukuma as well as two comparative breeds; Boran and Friesian. Population structure and signatures of selection were examined using principal component analysis (PCA), admixture analysis, pairwise distances (F_ST), integrated haplotype score (iHS), identical by state (IBS) and runs of homozygosity (ROH). There was a low level of inbreeding (F~0.01) in the TSZ population compared to the Boran and Friesian breeds. The analyses of F_ST, IBS and admixture identified no considerable differentiation between TSZ trains. Importantly, common ancestry in Boran and TSZ were revealed based on admixture and IBD, implying gene flow between two populations. In addition, Friesian ancestry was found in Boran. A few common significant iHS were detected, which may reflect influence of recent selection in each breed or strain.

Conclusions

Population admixture and selection signatures could be applied to develop conservation plan of TSZ cattle as well as future breeding programs in East African cattle.

Bovine P-selectin mediates leukocyte adhesion and is highly polymorphic in dairy breeds

Bovine P-selectin (SELP) mediates leukocyte rolling and primes leukocyte adhesion to endothelium, both essential for leukocyte recruitment to an infection site. We investigated SELP-mediated adhesion between bovine peripheral blood leukocytes (PBLs) and cultured bovine aortic endothelial cells pre-activated with lipopolysaccharide (LPS). We examined gene polymorphism for bovine selectins SELP, l-selectin (SELL) and E-selectin (SELE) and compared their SNP frequency between five dairy breeds (Holstein, Friesian, Jersey, Ayrshire and Brown Swiss). LPS treatment caused a rapid (10min) and slower (4h) enhancement of PBL adhesion (P<0.01). Antibody blocking of SELP inhibited LPS induced cell adhesion. SELP was highly polymorphic, with 9 of the 13 SNPs in its exons, whereas only three synonymous SNPs in SELL and one in SELE. The resulting amino acid changes for the three missense SELP SNP were located in the lectin domain and in two consensus repeat (CR) regions, CR2 and CR5. The Val475Met variant locus in the CR4 and CR5 linking region was very close to a predicted N-acetyl-d-glucosamine glycosylation site, which is likely to influence SELP function. The AA genotype was under-represented, only being found in 1% of 373 heifers genotyped from the 5 breeds (P=0.056), suggesting that AA homozygous animals carrying the Val475Met substitution for SELP may have compromised development. Our study thus confirmed that SELP mediates the attachment of PBL to endothelium and provides novel evidence that its high polymorphism is likely to affect biological function. This may potentially influence leukocyte migration and fertility, both key to successful performance in dairy cows.