Gene Expression Analysis Provides New Insights into the Mechanism of Intramuscular Fat Formation in Japanese Black Cattle

Circular RNAs and host genes act synergistically in regulating cellular processes and functions in skeletal myogenesis

Circular RNAs (circRNAs) are post-transcriptional regulators generated from backsplicing of pre-mRNAs of host genes. A major circRNA regulatory mechanism involves microRNA (miRNA) sequestering, relieving miRNA-blocked mRNAs for translation and functions. To investigate possible circRNA-host gene relationship, skeletal myogenesis is chosen as a study model for its developmental importance and for readily available muscle tissues from farm animals for studies at different myogenic stages. This review aims to provide an integrated interpretations on methodologies, regulatory mechanisms and possible host gene-circRNA synergistic functional relationships in skeletal myogenesis, focusing on myoblast differentiation and proliferation, core drivers of muscle formation in myogenesis, while other myogenic processes that play supportive roles in the structure, maintenance and function of muscle tissues are also briefly discussed. On literature review,thirty-two circRNAs derived from thirty-one host genes involved in various myogenic stages are identified; twenty-two (68.6 %) of these circRNAs regulate myogenesis by sequestering miRNAs to engage PI3K/AKT and other signaling pathways while four (12.5 %) are translated into proteins for functions. In circRNA-host gene relationship,ten (32.3 %) host genes are shown to regulate myogenesis,nine (29.0 %) are specific to skeletal muscle functions,and twelve (38.8 %) are linked to skeletal muscle disorders.Our analysis of skeletal myogenesis suggests that circRNAs and host genes act synergistically to regulate cellular functions. Such circRNA-host gene functional synergism may also be found in other major cellular processes. CircRNAs may have evolved later than miRNAs to counteract the suppressive effects of miRNAs and to augment host gene functions to further fine-tune gene regulation.

Genome-Wide Association Study of Conformation Traits in Brazilian Holstein Cattle

The linear conformation of animals exerts an influence on health, reproduction, production, and welfare, in addition to longevity, which directly affects the profitability of milk-producing farms. The objectives of this study were (1) to perform genome-wide association studies (GWASs) of conformation traits, namely the Rump, Feet and Legs, Mammary System, Dairy Strength, and Final Classification traits, and (2) to identify genes and related pathways involved in physiological processes associated with conformation traits in Brazilian Holstein cattle. Phenotypic and genotypic data from 2339 Holstein animals distributed across the states of Rio Grande do Sul, Paraná, São Paulo, and Minas Gerais were used. The genotypic data were obtained with a 100 K SNP marker panel. The single-step genome-wide association study (ssGWAS) method was employed in the analyses. Genes close to a significant SNP were identified in an interval of 100 kb up- and downstream using the Ensembl database available in the BioMart tool. The DAVID database was used to identify the main metabolic pathways and the STRING program was employed to create the gene regulatory network. In total, 36 significant SNPs were found on 15 chromosomes; 27 of these SNPs were linked to genes that may influence the traits studied. Fourteen genes most closely related to the studied traits were identified, as well as four genes that showed interactions in important metabolic pathways such as myogenesis, adipogenesis, and angiogenesis. Among the total genes, four were associated with myogenesis (TMOD2, TMOD3, CCND2, and CTBP2), three with angiogenesis (FGF23, FGF1, and SCG3), and four with adipogenesis and body size and development (C5H12orf4, CCND2, EMILIN1, and FGF6). These results contribute to a better understanding of the biological mechanisms underlying phenotypic variability in conformation traits in Brazilian Holstein cattle.

A review of the role of epigenetic studies for intramuscular fat deposition in beef cattle

Intramuscular fat (IMF) deposition profoundly influences meat quality and economic value in beef cattle production. Meanwhile, contemporary developments in epigenetics have opened new outlooks for understanding the molecular basics of IMF regulation, and it has become a key area of research for world scholars. Therefore, the aim of this paper was to provide insight and synthesis into the intricate relationship between epigenetic mechanisms and IMF deposition in beef cattle. The methodology involves a thorough analysis of existing literature, including pertinent books, academic journals, and online resources, to provide a comprehensive overview of the role of epigenetic studies in IMF deposition in beef cattle. This review summarizes the contemporary studies in epigenetic mechanisms in IMF regulation, high-resolution epigenomic mapping, single-cell epigenomics, multi-omics integration, epigenome editing approaches, longitudinal studies in cattle growth, environmental epigenetics, machine learning in epigenetics, ethical and regulatory considerations, and translation to industry practices from perspectives of IMF deposition in beef cattle. Moreover, this paper highlights DNA methylation, histone modifications, acetylation, phosphorylation, ubiquitylation, non-coding RNAs, DNA hydroxymethylation, epigenetic readers, writers, and erasers, chromatin immunoprecipitation followed by sequencing, whole genome bisulfite sequencing, epigenome-wide association studies, and their profound impact on the expression of crucial genes governing adipogenesis and lipid metabolism. Nutrition and stress also have significant influences on epigenetic modifications and IMF deposition. The key findings underscore the pivotal role of epigenetic studies in understanding and enhancing IMF deposition in beef cattle, with implications for precision livestock farming and ethical livestock management. In conclusion, this review highlights the crucial significance of epigenetic pathways and environmental factors in affecting IMF deposition in beef cattle, providing insightful information for improving the economics and meat quality of cattle production.

A single-cell atlas of bovine skeletal muscle reveals mechanisms regulating intramuscular adipogenesis and fibrogenesis

BACKGROUND

Intramuscular fat (IMF) and intramuscular connective tissue (IMC) are often seen in human myopathies and are central to beef quality. The mechanisms regulating their accumulation remain poorly understood. Here, we explored the possibility of using beef cattle as a novel model for mechanistic studies of intramuscular adipogenesis and fibrogenesis.

METHODS

Skeletal muscle single-cell RNAseq was performed on three cattle breeds, including Wagyu (high IMF), Brahman (abundant IMC but scarce IMF), and Wagyu/Brahman cross. Sophisticated bioinformatics analyses, including clustering analysis, gene set enrichment analyses, gene regulatory network construction, RNA velocity, pseudotime analysis, and cell-cell communication analysis, were performed to elucidate heterogeneities and differentiation processes of individual cell types and differences between cattle breeds. Experiments were conducted to validate the function and specificity of identified key regulatory and marker genes. Integrated analysis with multiple published human and non-human primate datasets was performed to identify common mechanisms.

RESULTS

A total of 32 708 cells and 21 clusters were identified, including fibro/adipogenic progenitor (FAP) and other resident and infiltrating cell types. We identified an endomysial adipogenic FAP subpopulation enriched for COL4A1 and CFD (log2FC = 3.19 and 1.92, respectively; P < 0.0001) and a perimysial fibrogenic FAP subpopulation enriched for COL1A1 and POSTN (log2FC = 1.83 and 0.87, respectively; P < 0.0001), both of which were likely derived from an unspecified subpopulation. Further analysis revealed more progressed adipogenic programming of Wagyu FAPs and more advanced fibrogenic programming of Brahman FAPs. Mechanistically, NAB2 drives CFD expression, which in turn promotes adipogenesis. CFD expression in FAPs of young cattle before the onset of intramuscular adipogenesis was predictive of IMF contents in adulthood (R2  = 0.885, P < 0.01). Similar adipogenic and fibrogenic FAPs were identified in humans and monkeys. In aged humans with metabolic syndrome and progressed Duchenne muscular dystrophy (DMD) patients, increased CFD expression was observed (P < 0.05 and P < 0.0001, respectively), which was positively correlated with adipogenic marker expression, including ADIPOQ (R2  = 0.303, P < 0.01; and R2  = 0.348, P < 0.01, respectively). The specificity of Postn/POSTN as a fibrogenic FAP marker was validated using a lineage-tracing mouse line. POSTN expression was elevated in Brahman FAPs (P < 0.0001) and DMD patients (P < 0.01) but not in aged humans. Strong interactions between vascular cells and FAPs were also identified.

CONCLUSIONS

Our study demonstrates the feasibility of beef cattle as a model for studying IMF and IMC. We illustrate the FAP programming during intramuscular adipogenesis and fibrogenesis and reveal the reliability of CFD as a predictor and biomarker of IMF accumulation in cattle and humans.

Accurate haplotype construction and detection of selection signatures enabled by high quality pig genome sequences

High-quality whole-genome resequencing in large-scale pig populations with pedigree structure and multiple breeds would enable accurate construction of haplotype and robust selection-signature detection. Here, we sequence 740 pigs, combine with 149 of our previously published resequencing data, retrieve 207 resequencing datasets, and form a panel of worldwide distributed wild boars, aboriginal and highly selected pigs with pedigree structures, amounting to 1096 genomes from 43 breeds. Combining with their haplotype-informative reads and pedigree structure, we accurately construct a panel of 1874 haploid genomes with 41,964,356 genetic variants. We further demonstrate its valuable applications in GWAS by identifying five novel loci for intramuscular fat content, and in genomic selection by increasing the accuracy of estimated breeding value by 36.7%. In evolutionary selection, we detect MUC13 gene under a long-term balancing selection, as well as NPR3 gene under positive selection for pig stature. Our study provides abundant genomic variations for robust selection-signature detection and accurate haplotypes for deciphering complex traits in pigs.

Genome-Wide Association Analysis of Muscle pH in Texel Sheep × Altay Sheep F2 Resource Population

pH was one of the important meat quality traits, which was an important factor affecting the storage/shelf life and quality of meat in meat production. In order to find a way to extend the storage/shelf life, the pH values (pH_45min, pH_24h, pH_48h and pH_72h) of the longissimus dorsi muscles in F₂ individuals of 462 Texel sheep × Altay sheep were determined, genotyping was performed using Illumina Ovine SNP 600 K BeadChip and whole genome resequencing technology, a genome-wide association analysis (GWAS) was used to screen the candidate genes and molecular markers for pH values related to the quality traits of mutton, and the effects of population stratification were detected by Q-Q plots. The results showed that the pH population stratification analysis did not find significant systemic bias, and there was no obvious population stratification effect. The results of the association analysis showed that 28 SNPs significantly associated with pH reached the level of genomic significance. The candidate gene associated with pH_45min was identified as the CCDC92 gene by gene annotation and a search of the literature. Candidate genes related to pH_24h were KDM4C, TGFB2 and GOT2 genes. The candidate genes related to pH_48h were MMP12 and MMP13 genes. The candidate genes related to pH_72h were HILPDA and FAT1 genes. Further bioinformatics analyses showed 24 gene ontology terms and five signaling pathways that were significantly enriched (p ≤ 0.05). Many terms and pathways were related to cellular components, processes of protein modification, the activity of protein dimerization and hydrolase activity. These identified SNPs and genes could provide useful information about meat and the storage/shelf life of meat, thereby extending the storage/shelf life and quality of meat.

Diversity in Cell Morphology, Composition, and Function among Adipose Depots in River Buffaloes

Fat deposition is a significant economic trait in livestock animals. Adipose tissues (ATs) developed in subcutaneous and visceral depots are considered waste whereas those within muscle are highly valued. In river buffaloes, lipogenesis is highly active in subcutaneous (especially in the sternum subcutaneous) and visceral depots but not in muscle tissue. Revealing the features and functions of ATs in different depots is significant for the regulation of their development. Here, we characterize the cell size, composition, and function of six AT depots in river buffaloes. Our data support that the subcutaneous AT depots have a larger cell size than visceral AT depots, and the subcutaneous AT depots, especially the sternum subcutaneous AT, are mainly associated with the extracellular matrix whereas the visceral AT depots are mainly associated with immunity. We found that sternum subcutaneous AT is significantly different from ATs in other depots, due to the high unsaturated fatty acid content and the significant association with metabolic protection. The perirenal AT is more active in FA oxidation for energy supply. In addition, the expression of HOX paralogs supports the variable origins of ATs in different depots, indicating that the development of ATs in different depots is mediated by their progenitor cells. The present study enhances our understanding of the cellular and molecular features, metabolism, and origin of AT depots in buffaloes, which is significant for the regulation of fat deposition and provides new insights into the features of AT depots in multiple discrete locations.

Candidate Genes and Gene Networks Change with Age in Japanese Black Cattle by Blood Transcriptome Analysis

Age is an important physiological factor that affects the metabolism and immune function of beef cattle. While there have been many studies using the blood transcriptome to study the effects of age on gene expression, few have been reported on beef cattle. To this end, we used the blood transcriptomes of Japanese black cattle at different ages as the study subjects and screened 1055, 345, and 1058 differential expressed genes (DEGs) in the calf vs. adult, adult vs. old, and calf vs. old comparison groups, respectively. The weighted co-expression network consisted of 1731 genes. Finally, blue, brown, and yellow age-specific modules were obtained, in which genes were enriched in signaling pathways related to growth and development and immune metabolic dysfunction, respectively. Protein-protein interaction (PPI) analysis showed gene interactions in each specific module, and 20 of the highest connectivity genes were chosen as potential hub genes. Finally, we identified 495, 244, and 1007 genes by exon-wide selection signature (EWSS) analysis of different comparison groups. Combining the results of hub genes, we found that VWF, PARVB, PRKCA, and TGFB1I1 could be used as candidate genes for growth and development stages of beef cattle. CORO2B and SDK1 could be used as candidate marker genes associated with aging. In conclusion, by comparing the blood transcriptome of calves, adult cattle, and old cattle, the candidate genes related to immunity and metabolism affected by age were identified, and the gene co-expression network of different age stages was constructed. It provides a data basis for exploring the growth, development, and aging of beef cattle.

Integrated Analysis of Long Non-Coding RNA and mRNA to Reveal Putative Candidate Genes Associated with Backfat Quality in Beijing Black Pig

Pigs' backfat quality has an important impact on the quality of pork and pork products and has a strong relationship with nutrition and sensory characteristics. This study aimed to identify the related candidate genes of backfat quality and to preliminary clarify the molecular regulatory mechanism underlying pig backfat quality phenotypes. Expression assessments of long non-coding RNA (lncRNA) and mRNA profiling in backfat from high-quality (firm and white) and low-quality (soft and yellow) Beijing Black pigs were performed by RNA sequencing. Significantly different expressions were observed in 610 protein-coding genes and 290 lncRNAs between the two groups. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway annotation showed that some candidate differentially expressed genes that participate in lipid-related pathways and pigmentation terms may play a role in backfat quality in pigs. The cis-target and trans-target genes were predicted to explore the regulatory function of lncRNAs, and integrative analyses of different expression lncRNAs targets and different expression genes were performed. The results showed the regulatory networks of lncRNA-mRNA related to backfat quality, and our study obtained strong candidate genes for backfat quality: ELOVL5, SCD, DGAT2, SLC24A5, and TYRP1, which were involved in fat metabolism, adipogenesis regulation, and pigmentation. To our knowledge, this study is the first to demonstrate the molecular genetic mechanisms of backfat quality in pigs, and these findings improve the current understanding of backfat quality mechanisms and provide a foundation for further studies.

Assessing Genetic Diversity and Searching for Selection Signatures by Comparison between the Indigenous Livni and Duroc Breeds in Local Livestock of the Central Region of Russia

Indigenous pig breeds are mainly associated with the adaptive capacity that is necessary to respond adequately to climate change, food security, and livelihood needs, and natural resources conservation. Livni pigs are an indigenous fat-type breed farmed in a single farm in the Orel region and located in the Central European part of the Russian Federation. To determine the genomic regions and genes that are affected by artificial selection, we conducted the comparative study of two pig breeds with different breeding histories and breeding objectives, i.e., the native fat-type Livni and meat-type Duroc breeds using the Porcine GGP HD BeadChip, which contains ~80,000 SNPs. To check the Livni pigs for possible admixture, the Landrace and the Large White breeds were included into the study of genetic diversity as these breeds participated in the formation of the Livni pigs. We observed the highest level of genetic diversity in Livni pigs compared to commercial breeds (UHE = 0.409 vs. 0.319–0.359, p < 0.001; AR = 1.995 vs. 1.894–1.964, p < 0.001). A slight excess of heterozygotes was found in all of the breeds. We identified 291 candidate genes, which were localized within the regions under putative selection, including 22 and 228 genes, which were specific for Livni and Duroc breeds, respectively, and 41 genes common for both breeds. A detailed analysis of the molecular functions identified the genes, which were related to the formation of meat and fat traits, and adaptation to environmental stress, including extreme temperatures, which were different between breeds. Our research results are useful for conservation and sustainable breeding of Livni breed, which shows a high level of genetic diversity. This makes Livni one of the valuable national pig genetic resources.

Polo-Like Kinase 2 Plays an Essential Role in Cytoprotection against MG132-Induced Proteasome Inhibition via Phosphorylation of Serine 19 in HSPB5

Protein homeostasis, including protein folding, refolding, and degradation, is thought to decline with aging. HSPB5 (also known as αB-crystallin) prevents target protein aggregation as a molecular chaperone and exhibits a cytoprotective function against various cell stresses. To elucidate the effect of HSPB5 on endoplasmic reticulum (ER) stress, we searched for novel binding proteins of HSPB5 using the proximity-dependent biotin labeling method. Proteins presumed to interact with HSPB5 in cells treated with the proteasome inhibitor MG132 were identified by a reversible biotin-binding capacity method combining tamavidin2-REV magnetic beads and mass spectrometry. We discovered a new binding protein for HSPB5, polo-like kinase 2 (PLK2), which is an apoptosis-related enzyme. The expression of PLK2 was upregulated by MG132 treatment, and it was co-localized with HSPB5 near the ER in L6 muscle cells. Inhibition of PLK2 decreased ER stress-induced phosphorylation of serine 19 in HSPB5 and increased apoptosis by activation of caspase 3 under ER stress. Overexpression of HSPB5 (WT) suppressed the ER stress-induced caspase 3 activity, but this was not observed with phospho-deficient HSPB5 (3A) mutants. These results clarify the role of HSPB5 phosphorylation during ER stress and suggest that the PLK2/HSPB5 pathway plays an essential role in cytoprotection against proteasome inhibition-induced ER stress.

Application of Mass Spectrometry for Determining the Geographic Production Area of Wagyu Beef

Japanese Black cattle (Japanese Wagyu) beef is attracting attention for its aroma and marbling, and its handling is increasing worldwide. Here, we focused on the origin discrimination of Wagyu beef and analyzed the nutritional components of Japanese Wagyu (produced in multiple prefectures of Japan), Hybrid Wagyu (a cross between Angus and Wagyu cattle born in Australia and transported to Japan), and Australian Wagyu beef using mass spectrometry (MS). Triple-quadrupole liquid chromatography-MS was used to clarify the molecular species of lipids in Wagyu beef. Fourteen classes of lipids were separated, and 128 different triacylglycerides (TGs) were detected. A simple comparative analysis of these TGs using high-performance liquid chromatography revealed significantly higher levels of triolein (C18:1/C18:1/C18:1; abbreviated OOO) and C18:1/C18:1/C16:1 (OOPo) in Japanese Wagyu. Wagyu elements beef were comprehensively analyzed using inductively coupled plasma (ICP)-MS and ICP-optical emission spectrometry. We found significant differences in the rubidium, cesium, and lithium levels of Japanese and Australian Wagyu beef. On comparing metabolites using gas chromatography-MS, we identified significant differences in the levels of amino acids and other components of the Japanese and Australian Wagyu beef. These results suggest the possibility of determining the origin of Wagyu cattle breeds using MS and genetic discrimination.

Genome-Wide Selection Signatures and Human-Mediated Introgression Events in Bos taurus indicus-influenced Composite Beef Cattle

Genetic introgression from interbreeding hybridization of European Bos taurus taurus (EBT) and Indian Bos taurus indicus (IBI) cattle breeds have been widely used to combine the climatic resilience of the IBI cattle and the higher productivity of EBT when forming new composite beef cattle (CB) populations. The subsequent breeding strategies have shifted their initial genomic compositions. To uncover population structure, signatures of selection, and potential introgression events in CB populations, high-density genotypes [containing 492,954 single nucleotide polymorphisms (SNPs) after the quality control] of 486 individuals from 15 cattle breeds, including EBT, IBI, and CB populations, along with two Bos grunniens genotypes as outgroup were used in this study. Then, in-depth population genetics analyses were performed for three CB breeds of Beefmaster, Brangus, and Santa Gertrudis. Neighbor-joining, principal components, and admixture analyses confirmed the historical introgression of EBT and IBI haplotypes into CB breeds. The f_dM statistics revealed that only 12.9% of CB populations' genetic components are of IBI origin. The results of signatures of selection analysis indicated different patterns of selection signals in the three CB breeds with primary pressure on pathways involved in protein processing and stress response in Beefmaster, cell proliferation regulation and immune response in Brangus, and amino acids and glucose metabolisms in Santa Gertrudis. An average of >90% of genomic regions underlying selection signatures were of EBT origin in the studied CB populations. Investigating the CB breeds' genome allows the estimation of EBT and IBI ancestral proportions and the locations within the genome where either taurine or indicine origin alleles are under selective pressure. Such findings highlight various opportunities to control the selection process more efficiently and explore complementarity at the genomic level in CB populations.

Production of Hydroxy Fatty Acids, Precursors of γ-Hexalactone, Contributes to the Characteristic Sweet Aroma of Beef

Aroma is an essential factor for meat quality. The meat of Japanese Black cattle exhibits fine marbling and a rich and sweet aroma with a characteristic lactone composition. The mechanism of lactone formation associated with beef aroma has not been elucidated. In this study, we examined the precursors of γ-hexalactone, an indicator of the sweet aroma of beef and identified the mechanism underlying γ-hexalactone production. A low-temperature vacuum system was used to prepare beef tallow from Japanese Black cattle and Holstein cattle. The odor components were identified using headspace–gas chromatography. The analysis revealed that γ-hexalactone, γ-dodecalactone, δ-tetradecalactone, and δ-hexadecalactone were present as sweet aroma components of beef tallow prepared from marbling and muscle. Since we previously reported that γ-hexalactone formation correlates with linoleic acid content in beef, we analyzed ten oxidized fatty acids derived from linoleic acid by liquid chromatography–triple quadrupole mass spectrometry and detected two hydroxy-octadecadienoic acids (9S-HODE and 13S-HODE) in beef tallow. Significant differences in arachidonic acid 15-lipoxygenase and cyclooxygenase protein expression levels among subcutaneous fat, intramuscular fat, and muscle tissue were observed. Our results suggest that the combination of linoleic acid and the expression of lipid oxidase derived from beef muscle and intramuscular fat produce hydroxy fatty acids that result in a sweet aroma.

Transcriptome Analysis of Bovine Rumen Tissue in Three Developmental Stages

Rumen development is a crucial physiological challenge for ruminants. However, the molecular mechanism regulating rumen development has not been clearly elucidated. In this study, we investigated genes involved in rumen development in 13 rumen tissues from three developmental stages (birth, youth, and adult) using RNA sequencing. We identified that 6,048 genes were differentially expressed among three developmental stages. Using weighted correlation network analysis, we found that 12 modules were significantly associated with developmental stages. Functional annotation and protein–protein interaction (PPI) network analysis revealed that CCNB1, CCNB2, IGF1, IGF2, HMGCL, BDH1, ACAT1, HMGCS2, and CREBBP involved in rumen development. Integrated transcriptome with GWAS information of carcass weight (CW), stomach weight (SW), marbling score (MS), backfat thickness (BFT), ribeye area (REA), and lean meat weight (LMW), we found that upregulated DEGs (fold change 0∼1) in birth–youth comparison were significantly enriched with GWAS signals of MS, downregulated DEGs (fold change >3) were significantly enriched with GWAS signals of SW, and fold change 0∼1 up/downregulated DEGs in birth–adult comparison were significantly enriched with GWAS signals of CW, LMW, REA, and BFT. Furthermore, we found that GWAS signals for CW, LMW, and REA were enriched in turquoise module, and GWAS signals for CW was enriched in lightgreen module. Our study provides novel insights into the molecular mechanism underlying rumen development in cattle and highlights an integrative analysis for illustrating the genetic architecture of beef complex traits.

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.