Database of cattle candidate genes and genetic markers for milk production and mastitis

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

SNP Analysis of TLR4 Promoter and Its Transcriptional Factor Binding Profile in Relevance to Bovine Subclinical Mastitis

Bovine mastitis is a complex infectious disease that develops in the mammary gland, predominantly caused by a bacterial infection of mammary tissue. Genetic variability of mastitis is well established and depends upon different quantitative trait loci (QTL) related to mastitis resistance or susceptibility. The susceptibility is often attributed to single-nucleotide polymorphisms (SNPs) in the variable cow breed genomes. Several global investigative attempts have resulted in studies mapping mastitis to the variations in the relevant genes. Reports have been attributed to dramatic genetic expression changes in Toll-Like Receptor 4 (TLR4) genes in mastitis-positive cows. However, the mechanism behind this variable genetic expression of TLR4 genes has been studied poorly. The present study aims to investigate SCM through various screening tests like somatic cell count (SCC), electric conductivity (EC), pH, and California mastitis test (CMT) in milk samples. This study also aims to investigate possible mechanisms behind this variable expression of TLR4 by comparative SNP evaluation and transcriptional factor profile mining. So that the important genetic mutations and effects thereof can be exploited in selecting specific breeds with higher mastitis resistance and milk yield. Seventy Holstein Frisian (HF) crossbred dairy cows were selected in the present study. The animals were screened based on various diagnostic tests (SCC, pH, EC, and CMT). Blood samples (5 mL) were collected for extraction of DNA followed by amplification of PPR1 and PPR2 of the promoter region and 5'UTR of the bovine TLR4 gene using specific primers. Sanger's enzymatic DNA sequencing technique sequenced the amplified PCR products. Further, the identification of SNPs was done through various bioinformatic tools used in this study. The findings of the present study revealed that CMT, EC, pH, and SCC could be used for the early detection of subclinical mastitis. In the present study, a significant increase in the EC, pH, and SCC in milk samples of animals affected with SCM was found in comparison to the healthy animals. The present study also revealed 16 SNPs falling in TLR4 promoter and 5' untranslated region (5'UTR) sequences in mastitis-positive genotypes compared to reference genomes. The study also investigates the potential transcriptional factor program deployed in response to variable mastitis development resistance. In the present study, the allelic and genotype frequencies of all SNP variants in the three regions viz., PPR1, PPR2, and 5'UTR, were the same indicating the absence of heterozygous condition at the respective loci. The present study has wide applicability for researchers developing mastitis-resistant breeding programs and the data generated may aid in the selection of better genetic breeds. The transcription factor binding profiles can serve as concrete leads about the studies on bovine mastitis at the molecular level and may also aid global research groups working on transcription factor (TF)-based molecular pathology of mastitis.

Whole-genome SNP markers reveal runs of homozygosity in indigenous cattle breeds of Pakistan

The runs of homozygosity (ROH) were identified in 14 Pakistani cattle breeds (n = 105) by genotyping with the Illumina 50 K SNP BeadChip. These breeds were categorized into Dairy, Dual, and Draft breeds based on their utility and production performance. We identified a total of 10,936 ROHs which mainly consisted of a high number of shorter segments (1-4 Mb). Dairy group exhibited the highest level of inbreeding (FROH: 0.078 ± 0.028) while the lowest (FROH: 0.002 ± 0.008) was observed in Dual group. In 48 genomic regions identified with a high frequency of ROH, 207 genes were detected in the three breed groups. A substantially higher number of ROH islands detected in dairy breeds indicated the impact of the positive selection pressure over the years. Important candidate genes and QTL were detected in the ROH islands associated with economic traits like milk production, reproduction, meat, carcass, and health traits in dairy cattle.

Genomic Characterization and Initial Insight into Mastitis-Associated SNP Profiles of Local Latvian Bos taurus Breeds

Latvia has two local Bos taurus breeds-Latvian Brown (LBG) and Latvian Blue (LZG)-characterized by a good adaptation to the local climate, longevity, and high fat and protein contents in milk. Since these are desired traits in the dairy industry, this study investigated the genetic background of the LBG and LZG breeds and identified the genetic factors associated with mastitis. Blood and semen samples were acquired, and whole genome sequencing was then performed to acquire a genomic sequence with at least 35× or 10× coverage. The heterozygosity, nucleotide diversity, and LD analysis indicated that LBG and LZG cows have similar levels of genetic diversity compared to those of other breeds. An analysis of the population structure revealed that each breed clustered together, but the overall differentiation between the breeds was small. The highest genetic variance was observed in the LZG breed compared with the LBG breed. Our results show that SNP rs721295390 is associated with mastitis in the LBG breed, and SNPs rs383806754, chr29:43998719CG>C, and rs462030680 are associated with mastitis in the LZG breed. This study shows that local Latvian LBG and LZG breeds have a pronounced genetic differentiation, with each one suggesting its own mastitis-associated SNP profile.

Comparative genomic analysis uncovers candidate genes related with milk production and adaptive traits in goat breeds

During the process of animal domestication, both natural and artificial selection cause variation in allele frequencies among populations. Identifying genomic areas of selection in domestic animals may aid in the detection of genomic areas linked to ecological and economic traits. We studied genomic variation in 140 worldwide goat individuals, including 75 Asian, 30 African and 35 European goats. We further carried out comparative population genomics to detect genomic regions under selection for adaptability to harsh conditions in local Asian ecotypes and also milk production traits in European commercial breeds. In addition, we estimated the genetic distances among 140 goat individuals. The results showed that among all studied goat groups, local breeds from West and South Asia emerged as an independent group. Our search for selection signatures in local goats from West and South Asia revealed candidate genes related to adaptation to hot climate (HSPB6, HSF4, VPS13A and NBEA genes) and immune response (IL7, IL5, IL23A and LRFN5) traits. Furthermore, selection signatures in European commercial goats involved several milk production related genes, such as VPS13C, NCAM2, TMPRSS15, CSN3 and ABCG2. The identified candidate genes could be the fundamental genetic resource for enhancement of goat production and environmental-adaptive traits, and as such they should be used in goat breeding programs to select more efficient breeds.

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.

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.

Non-Synonymous Variants in Fat QTL Genes among High- and Low-Milk-Yielding Indigenous Breeds

The effect of breed on milk components-fat, protein, lactose, and water-has been observed to be significant. As fat is one of the major price-determining factors for milk, exploring the variations in fat QTLs across breeds would shed light on the variable fat content in their milk. Here, on whole-genome sequencing, 25 differentially expressed hub or bottleneck fat QTLs were explored for variations across indigenous breeds. Out of these, 20 genes were identified as having nonsynonymous substitutions. A fixed SNP pattern in high-milk-yielding breeds in comparison to low-milk-yielding breeds was identified in the genes GHR, TLR4, LPIN1, CACNA1C, ZBTB16, ITGA1, ANK1, and NTG5E and, vice versa, in the genes MFGE8, FGF2, TLR4, LPIN1, NUP98, PTK2, ZTB16, DDIT3, and NT5E. The identified SNPs were ratified by pyrosequencing to prove that key differences exist in fat QTLs between the high- and low-milk-yielding breeds.

Candidate genes for mastitis resistance in dairy cattle: a data integration approach

BACKGROUND

Inflammation of the mammary tissue (mastitis) is one of the most detrimental health conditions in dairy ruminants and is considered the most economically important infectious disease of the dairy sector. Improving mastitis resistance is becoming an important goal in dairy ruminant breeding programmes. However, mastitis resistance is a complex trait and identification of mastitis-associated alleles in livestock is difficult. Currently, the only applicable approach to identify candidate loci for complex traits in large farm animals is to combine different information that supports the functionality of the identified genomic regions with respect to a complex trait.

METHODS

To identify the most promising candidate loci for mastitis resistance we integrated heterogeneous data from multiple sources and compiled the information into a comprehensive database of mastitis-associated candidate loci. Mastitis-associated candidate genes reported in association, expression, and mouse model studies were collected by searching the relevant literature and databases. The collected data were integrated into a single database, screened for overlaps, and used for gene set enrichment analysis.

RESULTS

The database contains candidate genes from association and expression studies and relevant transgenic mouse models. The 2448 collected candidate loci are evenly distributed across bovine chromosomes. Data integration and analysis revealed overlaps between different studies and/or with mastitis-associated QTL, revealing promising candidate genes for mastitis resistance.

CONCLUSION

Mastitis resistance is a complex trait influenced by numerous alleles. Based on the number of independent studies, we were able to prioritise candidate genes and propose a list of the 22 most promising. To our knowledge this is the most comprehensive database of mastitis associated candidate genes and could be helpful in selecting genes for functional validation studies.

Comprehensive RNA-seq profiling to evaluate the rabbit mammary gland transcriptome after mastitis

Mastitis is a relatively common disease in rabbit does. The aim of this study was to investigate a relationship between the severity of clinical signs and pathological observations and to analyze differentially expressed genes (DEGs) in the mammary gland with mastitis versus healthy mammary gland. The result showed that rectal temperatures of the rabbits with both mild mastitis and severe mastitis were higher than that of control. Cell counting results showed that the somatic cell count (SCC) only in milk of the rabbit with severe mastitis was significantly higher than that in the control group. However, the number of heterophils in the histological sections of mammary glands with mild mastitis was significantly higher than that of control. A total of 1,096 DEGs between the control and mastitis mammary glands was identified by RNA-sequencing (RNA-seq). Gene ontology (GO) showed that most of up-regulated genes were enriched in terms such as response to stimulus, signal transduction, and cell communication. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these genes were mostly enriched in the pathways such as Rap1 signaling pathway, proteoglycans in cancer, and PI3K-Akt signaling pathway. However, the downregulated genes were mainly enriched in metabolic processes and significantly involved in metabolic pathways. The data provides useful information to further dissect the molecular genetic mechanisms underlying rabbit mastitis, which is a prerequisite for designing effective intervention strategies.

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

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

How Can We Advance Integrative Biology Research in Animal Science in 21st Century? Experience at University of Ljubljana from 2002 to 2022

In this perspective analysis, we strive to answer the following question: how can we advance integrative biology research in the 21st century with lessons from animal science? At the University of Ljubljana, Biotechnical Faculty, Department of Animal Science, we share here our three lessons learned in the two decades from 2002 to 2022 that we believe could inform integrative biology, systems science, and animal science scholarship in other countries and geographies. Cultivating multiomics knowledge through a conceptual lens of integrative biology is crucial for life sciences research that can stand the test of diverse biological, clinical, and ecological contexts. Moreover, in an era of the current COVID-19 pandemic, animal nutrition and animal science, and the study of their interactions with human health (and vice versa) through integrative biology approaches hold enormous prospects and significance for systems medicine and ecosystem health.

Immune mechanisms, resistance genes, and their roles in the prevention of mastitis in dairy cows

Mastitis is one of the most important diseases of the mammary gland. The increased incidence of this disease in cows is due to the breeding of dairy cattle for higher yields, which is accompanied by an increased susceptibility to mastitis. Therefore, the difficulty involved with preventing this disease has increased. An integral part of current research is the elimination of mastitis in order to reduce the consumption of antibiotic drugs, thereby reducing the resistance of microorganisms and decreasing companies' economic losses due to mastitis (i.e. decreased milk yield, increased drug costs, and reduced milk supply). Susceptibility to mastitis is based on dairy cows' immunity, health, nutrition, and welfare. Thus, it is important to understand the immune processes in the body in order to increase the resistance of animals. Recently, various studies have focused on the selection of mastitis resistance genes. An important point is also the prevention of mastitis. This publication aims to describe the physiology of the mammary gland along with its immune mechanisms and to approximate their connection with potential mastitis resistance genes. This work describes various options for mastitis elimination and focuses on genetic selection and a closer specification of resistance genes to mastitis. Among the most promising resistance genes for mastitis, we consider CD14, CXCR1, lactoferrin, and lactoglobulin.

Genetic parameters and genome-wide association for milk production traits and somatic cell score in different lactation stages of Shanghai Holstein population

The aim of this study was to investigate the genetic parameters and genetic architectures of six milk production traits in the Shanghai Holstein population. The data used to estimate the genetic parameters consisted of 1,968,589 test-day records for 305,031 primiparous cows. Among the cows with phenotypes, 3,016 cows were genotyped with Illumina Bovine SNP50K BeadChip, GeneSeek Bovine 50K BeadChip, GeneSeek Bovine LD BeadChip v4, GeneSeek Bovine 150K BeadChip, or low-depth whole-genome sequencing. A genome-wide association study was performed to identify quantitative trait loci and genes associated with milk production traits in the Shanghai Holstein population using genotypes imputed to whole-genome sequences and both fixed and random model circulating probability unification and a mixed linear model with rMVP software. Estimated heritabilities (h2) varied from 0.04 to 0.14 for somatic cell score (SCS), 0.07 to 0.22 for fat percentage (FP), 0.09 to 0.27 for milk yield (MY), 0.06 to 0.23 for fat yield (FY), 0.09 to 0.26 for protein yield (PY), and 0.07 to 0.35 for protein percentage (PP), respectively. Within lactation, genetic correlations for SCS, FP, MY, FY, PY, and PP at different stages of lactation estimated in random regression model were ranged from -0.02 to 0.99, 0.18 to 0.99, 0.04 to 0.99, 0.04 to 0.99, 0.01 to 0.99, and 0.33 to 0.99, respectively. The genetic correlations were highest between adjacent DIM but decreased as DIM got further apart. Candidate genes included those related to production traits (DGAT1, MGST1, PTK2, and SCRIB), disease-related (LY6K, COL22A1, TECPR2, and PLCB1), heat stress–related (ITGA9, NDST4, TECPR2, and HSF1), and reproduction-related (7SK and DOCK2) genes. This study has shown that there are differences in the genetic mechanisms of milk production traits at different stages of lactation. Therefore, it is necessary to conduct research on milk production traits at different stages of lactation as different traits. Our results can also provide a theoretical basis for subsequent molecular breeding, especially for the novel genetic loci.

Low-density SNP panel for efficient imputation and genomic selection of milk production and technological traits in dairy sheep

The present study aimed to ascertain how different strategies for leveraging genomic information enhance the accuracy of estimated breeding values for milk and cheese-making traits and to evaluate the implementation of a low-density (LowD) SNP chip designed explicitly for that aim. Thus, milk samples from a total of 2,020 dairy ewes from 2 breeds (1,039 Spanish Assaf and 981 Churra) were collected and analyzed to determine 3 milk production and composition traits and 2 traits related to milk coagulation properties and cheese yield. The 2 studied populations were genotyped with a customized 50K Affymetrix SNP chip (Affymetrix Inc.) containing 55,627 SNP markers. The prediction accuracies were obtained using different multitrait methodologies, such as the BLUP model based on pedigree information, the genomic BLUP (GBLUP), and the BLUP at the SNP level (SNP-BLUP), which are based on genotypic data, and the single-step GBLUP (ssGBLUP), which combines both sources of information. All of these methods were analyzed by cross-validation, comparing predictions of the whole population with the test population sets. Additionally, we describe the design of a LowD SNP chip (3K) and its prediction accuracies through the different methods mentioned previously. Furthermore, the results obtained using the LowD SNP chip were compared with those based on the 50K SNP chip data sets. Finally, we conclude that implementing genomic selection through the ssGBLUP model in the current breeding programs would increase the accuracy of the estimated breeding values compared with the BLUP methodology in the Assaf (from 0.19 to 0.39) and Churra (from 0.27 to 0.44) dairy sheep populations. The LowD SNP chip is cost-effective and has proven to be an accurate tool for estimating genomic breeding values for milk and cheese-making traits, microsatellite imputation, and parentage verification. The results presented here suggest that the routine use of this LowD SNP chip could potentially increase the genetic gains of the breeding selection programs of the 2 Spanish dairy sheep breeds considered here.

Inputs for optimizing selection platform for milk production traits of dairy Sahiwal cattle

The premises for the potential success of molecular breeding is the ability to identify major genes associated with important dairy related traits. The present study was taken up with the objectives to identify single nucleotide polymorphism (SNP) of bovine MASP2 and SIRT1 genes and its effect on estimated breeding values (EBVs) and to estimate genetic parameters for lactation milk yield (LMY), 305-day milk yield (305dMY), 305-day fat yield (305dFY), 305-day solid not fat yield (305dSNFY) and lactation length (LL) in Sahiwal dairy cattle to devise a promising improvement strategy. Genetic parameters and breeding values of milk production traits were estimated from 935 Sahiwal cattle population (1979–2019) reared at National Dairy Research Institute at Karnal, India. A total of 7 SNPs, where one SNP (g.499C>T) in exon 2 and four SNPs (g.576G>A, g.609T>C, g.684G>T and g.845A>G) in exon 3 region of MASP2 gene and 2 SNPs (g.-306T>C and g.-274G>C) in the promoter region of SIRT1 gene were identified in Sahiwal cattle population. Five of these identified SNPs were chosen for further genotyping by PCR-RFLP and association analysis. Association analysis was performed using estimated breeding values (n = 150) to test the effect of SNPs on LMY, 305dMY, 305dFY, 305dSNFY and LL. Association analysis revealed that, three SNP markers (g.499C>T, g.609T>C and g.-306T>C) were significantly associated with all milk yield traits. The estimates for heritability using repeatability model for LMY, 305dMY, 305dFY, 305dSNFY and LL were low, however the corresponding estimates from first parity were 0.20±0.08, 0.17±0.08, 0.13±0.09, 0.13±0.09 and 0.24, respectively. The repeatability estimates were moderate to high indicating consistency of performance over the parities and hence reliability of first lactation traits. Genetic correlations among the traits of first parity were high (0.55 to 0.99). From the results we could conclude that optimum strategy to improve the Sahiwal cattle further would be selecting the animals based on their first lactation 305dMY. Option top include the significant SNP in selection criteria can be explored. Taken together, a 2-stage selection approach, select Sahiwal animals early for the SNP and then on the basis of first lactation 305dMY will help to save resources.

Deciphering the Molecular Mechanism Underlying African Animal Trypanosomiasis by Means of the 1000 Bull Genomes Project Genomic Dataset

Simple Summary

Climate change is increasing the risk of spreading vector-borne diseases such as African Animal Trypanosomiasis (AAT), which is causing major economic losses, especially in sub-Saharan African countries. Mainly considering this disease, we have investigated transcriptomic and genomic data from two cattle breeds, namely Boran and N‘Dama, where the former is known for its susceptibility and the latter one for its tolerance to the AAT. Despite the rich literature on this disease, there is still a need to investigate underlying genetic mechanisms to decipher the complex interplay of regulatory SNPs (rSNPs), their corresponding gene expression profiles and the downstream effectors associated with the AAT disease. The findings of this study complement our previous results, which mainly involve the upstream events, including transcription factors (TFs) and their co-operations as well as master regulators. Moreover, our investigation of significant rSNPs and effectors found in the liver, spleen and lymph node tissues of both cattle breeds could enhance the understanding of distinct mechanisms leading to either resistance or susceptibility of cattle breeds.

Abstract

African Animal Trypanosomiasis (AAT) is a neglected tropical disease and spreads by the vector tsetse fly, which carries the infectious Trypanosoma sp. in their saliva. Particularly, this parasitic disease affects the health of livestock, thereby imposing economic constraints on farmers, costing billions of dollars every year, especially in sub-Saharan African countries. Mainly considering the AAT disease as a multistage progression process, we previously performed upstream analysis to identify transcription factors (TFs), their co-operations, over-represented pathways and master regulators. However, downstream analysis, including effectors, corresponding gene expression profiles and their association with the regulatory SNPs (rSNPs), has not yet been established. Therefore, in this study, we aim to investigate the complex interplay of rSNPs, corresponding gene expression and downstream effectors with regard to the AAT disease progression based on two cattle breeds: trypanosusceptible Boran and trypanotolerant N’Dama. Our findings provide mechanistic insights into the effectors involved in the regulation of several signal transduction pathways, thereby differentiating the molecular mechanism with regard to the immune responses of the cattle breeds. The effectors and their associated genes (especially MAPKAPK5, CSK, DOK2, RAC1 and DNMT1) could be promising drug candidates as they orchestrate various downstream regulatory cascades in both cattle breeds.

Maternal Variants in the MFGE8 Gene are Associated with Perceived Breast Milk Supply

Background: The World Health Organization recommends exclusive breastfeeding for ≥6 months, but many mothers are unable to meet this goal. A major reason why mothers undergo early, unplanned breastfeeding cessation is perceived inadequate of milk supply (PIMS). We hypothesized that defining genetic polymorphisms associated with PIMS could aid early identification of at-risk mothers, providing an opportunity for targeted lactation support. Materials and Methods: This prospective observational cohort study followed 221 breastfeeding mothers for 12 months, collecting medical, demographic, and breastfeeding characteristics. Eighteen mammary secretory genes were assessed for single-nucleotide polymorphisms in 88 women (45 with PIMS and 43 with perceived adequate milk supply [PAMS]), matched by age/race/parity. Hierarchical regressions were used to assess the ability of genotype to aid PIMS prediction. Results: Mothers with PIMS exclusively breastfed for a shorter period (7 ± 12 weeks; p = 0.001) and reported lower milk production (17.6 ± 13.3 oz/day; p = 0.001), and their infants displayed reduced weight-for-length Z-score gains (0.74 ± 1.4; p = 0.038) relative to mothers with PAMS (22 ± 19 weeks; 27.03 ± 12.2 oz/day; 1.4 ± 1.5). Maternal genotype for the rs2271714 variant within milk fat globule EGF and factor V/VIII domain containing gene (MFGE8) was associated with PIMS status (p = 0.009, adjusted p = 0.09, likelihood ratio = 9.33) and duration of exclusive breastfeeding (p = 0.009, adjusted p = 0.09, χ² = 9.39). Addition of MFGE8 genotype to a model employing maternal characteristics (age, parity, previous breast-feeding duration, body mass index, education, and depression status) significantly increased predictive accuracy for PIMS status (p = 0.001; χ² = 13.5; area under the curve = 0.813 versus 0.725). Conclusions: Genotyping one lactogenic gene aided identification of mothers at risk for PIMS. If validated in a larger cohort, such an approach could be used to identify mothers who may benefit from increased lactation support.

Challenging Sustainable and Innovative Technologies in Cheese Production: A Review

It is well known that cheese yield and quality are affected by animal genetics, milk quality (chemical, physical, and microbiological), production technology, and the type of rennet and dairy cultures used in production. Major differences in the same type of cheese (i.e., hard cheese) are caused by the rennet and dairy cultures, which affect the ripening process. This review aims to explore current technological advancements in animal genetics, methods for the isolation and production of rennet and dairy cultures, along with possible applications of microencapsulation in rennet and dairy culture production, as well as the challenge posed to current dairy technologies by the preservation of biodiversity. Based on the reviewed scientific literature, it can be concluded that innovative approaches and the described techniques can significantly improve cheese production.

Single nucleotide polymorphisms and metabolic biochemical profile of productive markers characterize three European breeds of dairy cattle

The objective of this study was to investigate polymorphisms of DGAT1, FABP, OLR1 and ATP1A1 genes using PCR-DNA sequencing, and to associate these genetic structures to changes in metabolic biochemical markers and milk composition indicators in a total of 90 dairy cows of the Holstein, Simmental, and Brown Swiss breeds (30 cows each). PCR was carried out for amplification of 411-bp of DGAT1, 525-bp of FABP, 582-bp of OLR1, and 300-bp of ATP1A1 genes. Three breeds’ nucleotide sequence variations in the form of single nucleotide polymorphisms (SNPs) were detailed by DNA sequencing analysis. Chisquare analysis showed that the distribution of all discovered SNPs varied significantly (P < 0.001). Biochemical indices in cow’s serum revealed no significant difference in serum total protein, albumin, and total cholesterol among the three breeds. However, triglyceride showed a significant increase in Simmental compared to either Holsteins or Brown Swiss, while the highest mean value of triiodothyronine (T3) and tetraiodothyronine (T4) was detected in Holstein dairy cows The milk composition indicators analysis revealed that milk protein, sugar, and density were significantly higher in Holsteins than both Simmental and Brown Swiss. Meanwhile, milk fat and total solids revealed a significantly higher increase in Simmental than both brown Swiss and Holstein. As a result, the metabolic biochemical markers profile along with the identified SNPs could be used as a candidate and a reference guide for effective characterization of the Holstein, Simmental, and Brown Swiss breeds, leading to the creation of a marker-assisted selection system for production traits in dairy cattle breeds.

Genomewide Association Analyses of Lactation Persistency and Milk Production Traits in Holstein Cattle Based on Imputed Whole-Genome Sequence Data

Lactation persistency and milk production are among the most economically important traits in the dairy industry. In this study, we explored the association of over 6.1 million imputed whole-genome sequence variants with lactation persistency (LP), milk yield (MILK), fat yield (FAT), fat percentage (FAT%), protein yield (PROT), and protein percentage (PROT%) in North American Holstein cattle. We identified 49, 3991, 2607, 4459, 805, and 5519 SNPs significantly associated with LP, MILK, FAT, FAT%, PROT, and PROT%, respectively. Various known associations were confirmed while several novel candidate genes were also revealed, including ARHGAP35, NPAS1, TMEM160, ZC3H4, SAE1, ZMIZ1, PPIF, LDB2, ABI3, SERPINB6, and SERPINB9 for LP; NIM1K, ZNF131, GABRG1, GABRA2, DCHS1, and SPIDR for MILK; NR6A1, OLFML2A, EXT2, POLD1, GOT1, and ETV6 for FAT; DPP6, LRRC26, and the KCN gene family for FAT%; CDC14A, RTCA, HSTN, and ODAM for PROT; and HERC3, HERC5, LALBA, CCL28, and NEURL1 for PROT%. Most of these genes are involved in relevant gene ontology (GO) terms such as fatty acid homeostasis, transporter regulator activity, response to progesterone and estradiol, response to steroid hormones, and lactation. The significant genomic regions found contribute to a better understanding of the molecular mechanisms related to LP and milk production in North American Holstein cattle.

Single nucleotide polymorphism in STAT5A could not endorse variation in milk production traits in Indian bovine population

The Signal Transducer and Activator of Transcription 5A (STAT5A) gene involved in activating the transcription of milk protein genes was predicted to be influencing milk production traits. The present study was undertaken to investigate the suitability of the polymorphism of STAT5A as a marker for milk traits in Ongole, crossbred cattle and Murrah buffaloes from Southern India. Blood samples (n = 502) for DNA isolation and milk samples (n = 222) from different genetic groups were collected from various farms. The gene variants upon polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) on the exon 7 region of STAT5A were subjected to GLM analysis to evaluate their association with milk production traits. The frequencies of C and T alleles at the STAT5A/AvaI locus were 0.98 and 0.02 (Jersey crossbred), 0.94 and 0.06 [Holstein-Friesian (HF) crossbred], 0.97 and 0.03 (Ongole). T allele was not observed in Murrah buffaloes. The least squares mean lactation milk yield of CC and CT genotypes of STAT5A were 2,096.90 ± 48.63 and 2,294.41 ± 215.85 kg in Jersey crossbred, 2,312.92 ± 91.01 and 2,392.82 ± 207.66 kg in HF crossbred and 528.40 ± 22.10 and 396.37 ± 76.17 kg in Ongole cattle, respectively. The milk fat content of the CC genotype was higher (P > 0.05) in Jersey crossbred cattle. The CT genotypes of Ongole and HF crossbred cattle recorded a higher fat per cent than the CC genotypes. Significant associations were not observed in support of STAT5A as a marker for milk production traits in either Ongole or crossbred cattle of indicine admixture and no reason could be found to consider this locus as universal markers for milk production traits in indicine cattle and buffaloes. Considering the monomorphic nature of the gene in buffaloes and their higher milk fat content as compared to bovine milk, much remains to be explored regarding the underlying differences across the bovine and the bubaline species.

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.

Top Trends in Multiomics Research: Evaluation of 52 Published Studies and New Ways of Thinking Terminology and Visual Displays

Multiomics study designs have significantly increased understanding of complex biological systems. The multiomics literature is rapidly expanding and so is their heterogeneity. However, the intricacy and fragmentation of omics data are impeding further research. To examine current trends in multiomics field, we reviewed 52 articles from PubMed and Web of Science, which used an integrated omics approach, published between March 2006 and January 2021. From studies, data regarding investigated loci, species, omics type, and phenotype were extracted, curated, and streamlined according to standardized terminology, and summarized in a previously developed graphical summary. Evaluated studies included 21 omics types or applications of omics technology such as genomics, transcriptomics, metabolomics, epigenomics, environmental omics, and pharmacogenomics, species of various phyla including human, mouse, Arabidopsis thaliana, Saccharomyces cerevisiae, and various phenotypes, including cancer and COVID-19. In the analyzed studies, diverse methods, protocols, results, and terminology were used and accordingly, assessment of the studies was challenging. Adoption of standardized multiomics data presentation in the future will further buttress standardization of terminology and reporting of results in systems science. This shall catalyze, we suggest, innovation in both science communication and laboratory medicine by making available scientific knowledge that is easier to grasp, share, and harness toward medical breakthroughs.

Estimation of the associations between GH and DGAT1 genes and growth traits by using decision tree in Awassi sheep

GH and DGAT1 are candidate genes associated with growth traits in sheep breeds. This study aims to determine the association between growth traits and GH and DGAT1 gene polymorphism using three statistical methods in Awassi sheep. The polymorphism was detected by the PCR-RFLP method. Three genotypes (AA, AB and BB) were observed for GH gene locus with allele and genotype frequency 0.70(A) and 0.30(B); 0.60(AA), 0.20(AB) and 0.20(BB). Three genotypes (CC, CT and TT) were found for DGAT1 gene locus with allele and genotype frequency 0.58(C) and 0.42(T); 0.47(CC), 0.23(CT) and 0.30(TT). The genes were in agreement with Hardy-Weinberg equilibrium (p > 0.05). CHAID, CRT and GLM were used to identify the association between growth traits and gene polymorphism. The results showed an association between GH locus and body weight and tail length. Also, found an association between DGAT1 locus and tail length. The three methods showed similar results in determining the association between genes and growth traits. Thus, the CHAID and CRT methods can evaluate the association between genes and growth traits.

Epigenetic states of genes controlling immune responsiveness in bovine chronic mastitis

Mastitis is a common disease in dairy cows, causing substantial economic losses. The leading cause of elevated milk somatic cell count (SCC), which is the best indicator for detecting mastitis, is the invasion of pathogens. A major pathogen responsible for bovine mastitis is Staphylococcus aureus, a member of the coagulase-positive staphylococci. Some strains of coagulase-negative staphylococci can also be a cause of clinical or subclinical mastitis. Our study used bisulfite sequencing PCR (BSP) to detect the methylation status of nine candidate genes (CCL2, HCK, F11R, CD8A, PDIA3, LGMN, HSPA1A, IL18 and NFKBIA). We investigated the mechanisms associated with overexpression of these genes, in the mammary gland secretory tissue of cows diagnosed with mastitis and infected with coagulase-positive or coagulase-negative staphylococci. The results showed no changes at the DNA methylation level between the mastitis (CoPS and CoNS) and control groups (H), except for in the HCK region, where the observed differences between the CoPS and H groups were statistically significant. The low methylation level of the CpG sequence seems not to correspond to the previously observed increased activity of these genes, suggesting that mechanisms other than DNA methylation may control mRNA expression at the analyzed loci.

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 studies (GWAS) and post-GWAS analyses for technological traits in Assaf and Churra dairy breeds

This study aimed to perform a GWAS to identify genomic regions associated with milk and cheese-making traits in Assaf and Churra dairy sheep breeds; second, it aimed to identify possible positional and functional candidate genes and their interactions through post-GWAS studies. For 2,020 dairy ewes from 2 breeds (1,039 Spanish Assaf and 981 Churra), milk samples were collected and analyzed to determine 6 milk production and composition traits and 6 traits related to milk coagulation properties and cheese yield. The genetic profiles of the ewes were obtained using a genotyping chip array that included 50,934 SNP markers. For both milk and cheese-making traits, separate single-breed GWAS were performed using GCTA software. The set of positional candidate genes identified via GWAS was subjected to guilt-by-association-based prioritization analysis with ToppGene software. Totals of 84 and 139 chromosome-wise significant associations for the 6 milk traits and the 6 cheese-making traits were identified in this study. No significant SNPs were found in common between the 2 studied breeds, possibly due to their genetic heterogeneity of the phenotypes under study. Additionally, 63 and 176 positional candidate genes were located in the genomic intervals defined as confidence regions in relation to the significant SNPs identified for the analyzed traits for Assaf and Churra breeds. After the functional prioritization analysis, 71 genes were identified as promising positional and functional candidate genes and proposed as targets of future research to identify putative causative variants in relation to the traits under examination. In addition, this multitrait study allowed us to identify variants that have a pleiotropic effect on both milk production and cheese-related traits. The incorporation of variants among the proposed functional and positional candidate genes into genomic selection strategies represent an interesting approach for achieving rapid genetic gains, specifically for those traits difficult to measure, such as cheese-making traits.

Evaluation of candidate genotype of immune gene MBL1 associated with udder health and performance traits in dairy cattle and buffalo of India

An investigation was conducted to identify polymorphism in mannose-binding lectin 1 (MBL1) gene and its effect on udder health and performance traits in dairy cattle and buffalo of India. Candidate single-nucleotide polymorphism (SNP) c.2534G > A of MBL1 gene was genotyped by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). All the possible genotypes for SNP c.2534G > A (GG, AG, and AA) were observed in the studied population. However, Sahiwal cows revealed dimorphic pattern (AG and GG). The effect of targeted SNP on incidence of mastitis was evaluated and found to be significant. Animals with GG genotype were less susceptible to clinical mastitis and had comparatively lower somatic cell score (SCS) in Hardhenu cattle (P < 0.01) and Murrah buffalo (P < 0.05). Animals having GG genotype also exhibited significantly (P < 0.05) lower age at first calving (AFC). AG genotyped Murrah buffalo animals revealed significantly higher second lactation milk yield (P < 0.01). GG genotype with SCS and AFC could therefore be exploited as a promising candidate marker for the genetic improvement of udder health and AFC in dairy animals.

Sequence-based genome-wide association study of individual milk mid-infrared wavenumbers in mixed-breed dairy cattle

BACKGROUND

Fourier-transform mid-infrared (FT-MIR) spectroscopy provides a high-throughput and inexpensive method for predicting milk composition and other novel traits from milk samples. While there have been many genome-wide association studies (GWAS) conducted on FT-MIR predicted traits, there have been few GWAS for individual FT-MIR wavenumbers. Using imputed whole-genome sequence for 38,085 mixed-breed New Zealand dairy cattle, we conducted GWAS on 895 individual FT-MIR wavenumber phenotypes, and assessed the value of these direct phenotypes for identifying candidate causal genes and variants, and improving our understanding of the physico-chemical properties of milk.

RESULTS

Separate GWAS conducted for each of 895 individual FT-MIR wavenumber phenotypes, identified 450 1-Mbp genomic regions with significant FT-MIR wavenumber QTL, compared to 246 1-Mbp genomic regions with QTL identified for FT-MIR predicted milk composition traits. Use of mammary RNA-seq data and gene annotation information identified 38 co-localized and co-segregating expression QTL (eQTL), and 31 protein-sequence mutations for FT-MIR wavenumber phenotypes, the latter including a null mutation in the ABO gene that has a potential role in changing milk oligosaccharide profiles. For the candidate causative genes implicated in these analyses, we examined the strength of association between relevant loci and each wavenumber across the mid-infrared spectrum. This revealed shared association patterns for groups of genomically-distant loci, highlighting clusters of loci linked through their biological roles in lactation and their presumed impacts on the chemical composition of milk.

CONCLUSIONS

This study demonstrates the utility of FT-MIR wavenumber phenotypes for improving our understanding of milk composition, presenting a larger number of QTL and putative causative genes and variants than found from FT-MIR predicted composition traits. Examining patterns of significance across the mid-infrared spectrum for loci of interest further highlighted commonalities of association, which likely reflects the physico-chemical properties of milk constituents.

Genome Wide Prediction, Mapping and Development of Genomic Resources of Mastitis Associated Genes in Water Buffalo

Water buffalo (Bubalus bubalis) are an important animal resource that contributes milk, meat, leather, dairy products, and power for plowing and transport. However, mastitis, a bacterial disease affecting milk production and reproduction efficiency, is most prevalent in populations having intensive selection for higher milk yield, especially where the inbreeding level is also high. Climate change and poor hygiene management practices further complicate the issue. The management of this disease faces major challenges, like antibiotic resistance, maximum residue level, horizontal gene transfer, and limited success in resistance breeding. Bovine mastitis genome wide association studies have had limited success due to breed differences, sample sizes, and minor allele frequency, lowering the power to detect the diseases associated with SNPs. In this work, we focused on the application of targeted gene panels (TGPs) in screening for candidate gene association analysis, and how this approach overcomes the limitation of genome wide association studies. This work will facilitate the targeted sequencing of buffalo genomic regions with high depth coverage required to mine the extremely rare variants potentially associated with buffalo mastitis. Although the whole genome assembly of water buffalo is available, neither mastitis genes are predicted nor TGP in the form of web-genomic resources are available for future variant mining and association studies. Out of the 129 mastitis associated genes of cattle, 101 were completely mapped on the buffalo genome to make TGP. This further helped in identifying rare variants in water buffalo. Eighty-five genes were validated in the buffalo gene expression atlas, with the RNA-Seq data of 50 tissues. The functions of 97 genes were predicted, revealing 225 pathways. The mastitis proteins were used for protein-protein interaction network analysis to obtain additional cross-talking proteins. A total of 1,306 SNPs and 152 indels were identified from 101 genes. Water Buffalo-MSTdb was developed with 3-tier architecture to retrieve mastitis associated genes having genomic coordinates with chromosomal details for TGP sequencing for mining of minor alleles for further association studies. Lastly, a web-genomic resource was made available to mine variants of targeted gene panels in buffalo for mastitis resistance breeding in an endeavor to ensure improved productivity and the reproductive efficiency of water buffalo.

Genomic Loci Affecting Milk Production in German Black Pied Cattle (DSN)

German Black Pied cattle (DSN) is an endangered population of about 2,550 dual-purpose cattle in Germany. Having a milk yield of about 2,500 kg less than the predominant dairy breed Holstein, the preservation of DSN is supported by the German government and the EU. The identification of the genomic loci affecting milk production in DSN can provide a basis for selection decisions for genetic improvement of DSN in order to increase market chances through the improvement of milk yield. A genome-wide association analysis of 30 milk traits was conducted in different lactation periods and numbers. Association using multiple linear regression models in R was performed on 1,490 DSN cattle genotyped with BovineSNP50 SNP-chip. 41 significant and 20 suggestive SNPs affecting milk production traits in DSN were identified, as well as 15 additional SNPs for protein content which are less reliable due to high inflation. The most significant effects on milk yield in DSN were detected on chromosomes 1, 6, and 20. The region on chromosome 6 was located nearby the casein gene cluster and the corresponding haplotype overlapped the CSN3 gene (casein kappa). Associations for fat and protein yield and content were also detected. High correlation between traits of the same lactation period or number led to some SNPs being significant for multiple investigated traits. Half of all identified SNPs have been reported in other studies, previously. 15 SNPs were associated with the same traits in other breeds. The other associated SNPs have been reported previously for traits such as exterior, health, meat and carcass, production, and reproduction traits. No association could be detected between DGAT1 and other known milk genes with milk production traits despite the close relationship between DSN and Holstein. The results of this study confirmed that many SNPs identified in other breeds as associated with milk traits also affect milk traits in dual-purpose DSN cattle and can be used for further genetic analysis to identify genes and causal variants that affect milk production in DSN cattle.

Evaluation of non-synonym mutation in DGAT1 K232A as a marker for milk production traits in Ongole cattle and Murrah buffalo from Southern India

Various candidate genes have been reported to affect milk yield and composition in dairy cattle. A non-synonymous mutation in the DGAT1 gene, i.e., K232A was reported to have a strong association with milk yield and milk composition of Bos taurus. A study has been undertaken on 502 unrelated individuals belonging to indigenous Ongole cattle, crossbred cattle, and Murrah buffaloes from the Indian sub-continent with the objective to determine the polymorphism of the K232A locus and their association with milk yield and composition. Typing DGAT1 K232A allelic variation by PCR-RFLP using CfrI restriction enzyme revealed three genotypes in crossbred cattle. Genotype KK was more prevalent (0.60) in Jersey crossbred, whereas in Holstein Friesian crossbred it was KA genotype (0.48). In Ongole cattle and Murrah buffaloes, the locus did not exhibit polymorphism. The least-square mean of milk yields pooled over lactations across the DGAT1 variants was significantly (P < 0.05) higher among the homozygous (AA) genotypes, both in Jersey crossbred and HF crossbred cattle after adjusting for the effects of farm, parity, and season. The fat, SNF, and protein content values of AA genotypes were less than the KK genotypes in both the genetic groups (P > 0.05). The fixation of the DGAT1^K allele at the locus in Bos indicus cattle and Bubalus bubalis in the present study did not support its use as a reliable universal marker for milk production and composition traits.

The bioactivity of colostrum and milk exosomes of high, average, and low immune responder cows on human intestinal epithelial cells

Bovine milk contains bioactive components that are nutritionally and immunologically important to calves and humans. Dairy cows classified as high (H) immune responders using the patented high immune response technology have higher concentrations of immunoglobulin and specific antibodies in sera and milk compared with average (A) and low (L) responders. MicroRNA post-transcriptionally regulate expression of milk bioactive components and are enriched in extracellular vesicles known as exosomes, which protect them from degradation. The bioactivity of colostrum and milk exosomes at the human intestinal epithelial barrier remains to be explored, particularly in the context of the high immune response technology. Therefore, the purpose of this study was to evaluate the functional role of bovine milk exosomes compared with colostrum exosomes from H, A, and L responders at the intestinal interface using human colorectal adenocarcinoma epithelial (Caco-2) cells. Exosomes were isolated by successive ultracentrifugation and confirmed by western blot analysis for the presence of common exosomal proteins (CD9, CD63, and heat shock protein 70). Fluorescent labeling of exosomes using PKH67 dye confirmed their uptake by Caco-2 cells, demonstrating their potential bioavailability. The MTT assays showed that colostrum and milk exosomes maintain Caco-2 metabolic activity and are not cytotoxic to these cells. Specifically, metabolic activity after co-incubation with colostrum and milk exosomes from H responder cows was significantly greater than after co-incubation with exosomes from L responders. Caspase 3 activity, an indicator of apoptosis, was significantly lower after co-incubation of Caco-2 cells with milk exosomes compared with colostrum exosomes, suggesting that unlike colostrum exosomes, particularly those from L responders, milk exosomes do not activate the caspase 3 pathway in Caco-2 cells. This study helps us better understand the functional importance of colostrum and milk exosomes from dairy cows and emphasizes differences in functionality among exosomes from H, A, and L immune responders.

Significance of molecular chaperones and micro RNAs in acquisition of thermo-tolerance in dairy cattle

Ambient temperature is considered as the major abiotic factor which regulates body physiological mechanisms of all living creatures across the globe. Variation in ambient temperature which emulates thermoneutral zone culminates in heat stress. Heat stress has been emerged as major ultimatum to livestock's growth, development, production and reproduction across the world. Livestock's responds to the heat stress via different mechanisms such as behavioral, physiological, biochemical, endocrine and molecular mechanisms. Amongst the aforementioned mechanisms, molecular mechanism plays crucial role to achieve thermo-tolerance via expression of highly conserved family of proteins known as heat shock proteins (HSPs) across livestock species. HSPs serve as molecular chaperones to ameliorate the menace of heat stress in domestic species. In addition, microRNAs are small non-coding RNA which down regulates post-transcriptional gene expression by targeting various HSPs to regulate the thermoregulatory responses in livestock species. Despite of thermal adaptation mechanisms, heat stress breaches animal body homeostasis thereby depresses their production and productivity. Therefore, veterinary researches have been targeting to explore different repertoire of HSPs and microRNAs expression to counteract the rigors of heat stress thereby confer thermo-tolerance in livestock species. The present review highlights the significance of molecular chaperones and microRNAs in the acquisition of thermo-tolerance in dairy cattle.

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.

Association of GH, STAT5A, MYF5 gene polymorphisms with milk somatic cell count, EC and pH levels of Holstein dairy cattle

This study was conducted to ivnestigate the associations of GH-AluI, STAT5A-AvaI and MYF5-TaqI gene polymorphisms with milk somatic cell count (SCC), electrical conductivity (EC) and pH levels in Holstein dairy cows. For this purpose, 167 blood and 1670 milk samples of 167 Holstein cows in their 2nd lactation were used. There were significant relationships between GH-AluI genotypes and milk EC (p < 0.001) and between STAT5A-AvaI genotypes and milk EC (p = 0.007), but there were not any significant relationships between MYF5 gene polymorphism and the investigated traits (p > 0.05). The greatest EC values were observed in GH-AluI-LV and STAT5A-AvaI-TT-genotyped individuals. Just because of association of EC with mastitis, it was concluded that present GH-AluI and STAT5A-AvaI polymorphisms could be used in further studies to be conducted to improve mastitis resistance and milk quality traits of Holstein dairy cows.

Association Analysis between SPP1, POFUT1 and PRLR Gene Variation and Milk Yield, Composition and Coagulation Traits in Sarda Sheep

Many studies focus on the identification of genomic regions that undergo selective processes, where evidence of selection is revealed and positional candidate genes are identified. The aim of the research was to evaluate the association between positional candidate genes, namely secreted phosphoprotein 1 (SPP1, sheep chromosome Ovis aries OAR6, 36.651-36.658 Mb), protein O-fucosyltransferase 1 (POFUT1, OAR13, 61.006-61.027 Mb) and prolactin receptor (PRLR, OAR16, 38.969-39.028 Mb) with milk yield, composition and coagulation traits. Eight single nucleotide polymorphisms (SNPs) mapping to the three genes were genotyped in 380 Sarda dairy sheep. Statistical analysis revealed an association between SNP rs161844011 at SPP1 (chromosome position Oar_v3 OAR6:36651870, gene region exon 7) and somatic cell score, while POFUT1 SNP rs424501869 (OAR13:61007495, intron 1) was associated with curd firmness both 45 and 60 min after rennet addition (p = 0.015 and p = 0.007, respectively). SNP rs400874750 at PRLR gene (OAR16:39004070, intron 2) had a significant association with lactose content (p = 0.020), somatic cell score (p = 0.038), rennet coagulation time (p = 0.018) and curd firming time (p = 0.047). The outcome of this research confirmed predictions based on genomic studies, producing new information regarding the SPP1, POFUT1 and PRLR genes, which may be useful for future breeding schemes.

Meta-analysis for milk fat and protein percentage using imputed sequence variant genotypes in 94,321 cattle from eight cattle breeds

Background

Sequence-based genome-wide association studies (GWAS) provide high statistical power to identify candidate causal mutations when a large number of individuals with both sequence variant genotypes and phenotypes is available. A meta-analysis combines summary statistics from multiple GWAS and increases the power to detect trait-associated variants without requiring access to data at the individual level of the GWAS mapping cohorts. Because linkage disequilibrium between adjacent markers is conserved only over short distances across breeds, a multi-breed meta-analysis can improve mapping precision.

Results

To maximise the power to identify quantitative trait loci (QTL), we combined the results of nine within-population GWAS that used imputed sequence variant genotypes of 94,321 cattle from eight breeds, to perform a large-scale meta-analysis for fat and protein percentage in cattle. The meta-analysis detected (p ≤ 10⁻⁸) 138 QTL for fat percentage and 176 QTL for protein percentage. This was more than the number of QTL detected in all within-population GWAS together (124 QTL for fat percentage and 104 QTL for protein percentage). Among all the lead variants, 100 QTL for fat percentage and 114 QTL for protein percentage had the same direction of effect in all within-population GWAS. This indicates either persistence of the linkage phase between the causal variant and the lead variant across breeds or that some of the lead variants might indeed be causal or tightly linked with causal variants. The percentage of intergenic variants was substantially lower for significant variants than for non-significant variants, and significant variants had mostly moderate to high minor allele frequencies. Significant variants were also clustered in genes that are known to be relevant for fat and protein percentages in milk.

Conclusions

Our study identified a large number of QTL associated with fat and protein percentage in dairy cattle. We demonstrated that large-scale multi-breed meta-analysis reveals more QTL at the nucleotide resolution than within-population GWAS. Significant variants were more often located in genic regions than non-significant variants and a large part of them was located in potentially regulatory regions.

Genome-Wide Identification of Discriminative Genetic Variations in Beef and Dairy Cattle via an Information-Theoretic Approach

Analyzing the associations between genotypic changes and phenotypic traits on a genome-wide scale can contribute to understanding the functional roles of distinct genetic variations during breed development. We performed a whole-genome analysis of Angus and Jersey cattle breeds using conditional mutual information, which is an information-theoretic method estimating the conditional independency among multiple factor variables. The proposed conditional mutual information-based approach allows breed-discriminative genetic variations to be explicitly identified from tens of millions of SNP (single nucleotide polymorphism) positions on a genome-wide scale while minimizing the usage of prior knowledge. Using this data-driven approach, we identified biologically relevant functional genes, including breed-specific variants for cattle traits such as beef and dairy production. The identified lipid-related genes were shown to be significantly associated with lipid and triglyceride metabolism, fat cell differentiation, and muscle development. In addition, we confirmed that milk-related genes are involved in mammary gland development, lactation, and mastitis-associated processes. Our results provide the distinct properties of Angus and Jersey cattle at a genome-wide level. Moreover, this study offers important insights into discovering unrevealed genetic variants for breed-specific traits and the identification of genetic signatures of diverse cattle breeds with respect to target breed-specific properties.

Distinguishing pleiotropy from linked QTL between milk production traits and mastitis resistance in Nordic Holstein cattle

BACKGROUND

Production and health traits are central in cattle breeding. Advances in next-generation sequencing technologies and genotype imputation have increased the resolution of gene mapping based on genome-wide association studies (GWAS). Thus, numerous candidate genes that affect milk yield, milk composition, and mastitis resistance in dairy cattle are reported in the literature. Effect-bearing variants often affect multiple traits. Because the detection of overlapping quantitative trait loci (QTL) regions from single-trait GWAS is too inaccurate and subjective, multi-trait analysis is a better approach to detect pleiotropic effects of variants in candidate genes. However, large sample sizes are required to achieve sufficient power. Multi-trait meta-analysis is one approach to deal with this problem. Thus, we performed two multi-trait meta-analyses, one for three milk production traits (milk yield, protein yield and fat yield), and one for milk yield and mastitis resistance.

RESULTS

For highly correlated traits, the power to detect pleiotropy was increased by multi-trait meta-analysis compared with the subjective assessment of overlapping of single-trait QTL confidence intervals. Pleiotropic effects of lead single nucleotide polymorphisms (SNPs) that were detected from the multi-trait meta-analysis were confirmed by bivariate association analysis. The previously reported pleiotropic effects of variants within the DGAT1 and MGST1 genes on three milk production traits, and pleiotropic effects of variants in GHR on milk yield and fat yield were confirmed. Furthermore, our results suggested that variants in KCTD16, KCNK18 and ENSBTAG00000023629 had pleiotropic effects on milk production traits. For milk yield and mastitis resistance, we identified possible pleiotropic effects of variants in two genes, GC and DGAT1.

CONCLUSIONS

Multi-trait meta-analysis improves our ability to detect pleiotropic interactions between milk production traits and identifies variants with pleiotropic effects on milk production traits and mastitis resistance. In particular, this should contribute to better understand the biological mechanisms that underlie the unfavorable genetic correlation between milk yield and mastitis.

The effect of immunogenetic factors on productive longevity and reproductive qualities of cows

Increasing the productive longevity of cows with high reproductive ability is an urgent task in dairy cattle breeding both in Russia and foreign countries. Research was conducted in Krasnoyarsk region. The object was retired cows (n = 484). The blood groups of cows, indicators of productive longevity, and reproductive ability of cows were studied. The number of antigen carriers and non-carriers (frequency 30.0% and more) was analyzed. Out of 21 antigens with a frequency of more than 30.0%, the most frequently encountered antigens were E (80.6%), H` (77.5%), L (58.5%), W (57.6 %), A2 (57.0%), Z (54.1%), C (52.5%) and X2 (50.4%). The cows carrying O1 antigen had the highest productive longevity, – the difference from the cows not carrying this antigen by the retirement and lactation age was 0.8 years (p <0.001), and the period of economic use was 9.5 months (p <0.001), life-long milk yield was 3852 kg (p <0.001). Cows carrying antigens G3, O1, and L` had the best reproductive ability. In carriers of antigens O1, G3, and E`3, the service period was lower by 22.9; 20.3 and 18.9 days (p <0.05), and the fertility index was higher by 1.8; 1.6 and 1.7 units (p <0.05). In cows carrying L` antigen, the service period was shorter by 16.9 days (p <0.05).