Comparison of rumen contents' characteristics in Nguni and Bonsmara cows raised under two different grazing systems.

Objective: This study aimed to evaluate rumen fermentation parameters influenced by both grazing system and breed. Materials and Methods: A 2 × 2 factorial design was employed, involving 40 cows with matched age, parity, and physiological status. The cows were evenly divided between Bonsmara and Nguni breeds, as well as communal and commercial grazing systems. Rumen fluid samples were collected and analyzed for parameters including ammonia-nitrogen (NH3-N), pH, temperature, and volatile fatty acids (VFAs). Results: Nguni cows exhibited significantly higher ruminal NH3-N levels (p < 0.05) compared to Bonsmara, ranging from 69.05 to 96.78 mg/l. Commercial grazing demonstrated significantly higher NH3-N concentrations (p < 0.05) than communal grazing. Ruminal pH, temperature, total VFAs, and specific VFAs (Iso-butyrate, valeric, and iso-valeric) did not show significant differences (p > 0.05). However, total VFAs were slightly lower in communal grazing (78.87 mmol/l) than in commercial grazing (89.80 mmol/l). Acetate, propionate, butyrate, and the acetate to propionate ratio did not display significant differences (p > 0.05) between breeds but varied between grazing systems. Communal systems had higher acetate and acetate to propionate ratio (p < 0.05), while commercial systems showed higher propionate and butyrate levels (p < 0.05). Conclusion: Grazing conditions significantly influenced rumen fermentation parameters, irrespective of breed. Further research is necessary to explore the relationship between forage conditions, diversity, and rumen fermentation within different grazing systems.

Genetic differentiation and structuration of the Gobra zebu cattle breeds reared in Senegal.

The Gobra zebu genetic breeding program has resulted in the genetic improvement of a new population. This population gained genetic characteristics that set them apart from the other cattle populations reared in Senegal. The cause of these differences might be the reproductive isolation and selection to which this population of the "Centre de Recherches Zootechniques" of Dahra has been subjected since the 1950s. This study aimed to assess the genetic differentiation and structuration of this population compared to the main cattle breeds used in Senegal. A total of 180 individuals, selected from the Gobra selection nucleus and bovine populations from four main breeds in Senegal, were included in this study. We used a panel of 21 microsatellite markers among those recommended by the Food Agriculture Organization, to conduct the molecular genotyping of our sampled populations. The basic genetic parameters of differentiation and structuration were calculated using various bioinformatics software. The results of this study, particularly the degrees of genetic differentiation (Fst), the coefficient of genetic homogeneity (Gst), and the gene flow (Nm), show a significant genetic differentiation of the Gobra from the station compared to the other populations studied. Structuring and phylogeny analyses reveal a micro-structuring within the Gobra population as a novelty. This micro-structuring clearly identifies the Gobra individuals from Dahra's station among the other Gobra populations studied. The main causes of these observations would be reproductive isolation and the selection pressure exerted on this population for several decades.

Genome-wide screening for selection signatures in native and cosmopolitan cattle breeds reared in Türkiye.

Via long-term natural and artificial selection pressure, homozygosity may extend across the genome, leaving genomic patterns called selection signatures. This study is the first attempt to assess genome-wide selection signatures in six native Turkish and two cosmopolitan cattle breeds by 211.119 bi-allelic SNPs recovered using the double digest restriction associated DNA sequencing method. The integrated haplotype score (iHS) statistic was utilised to reveal selection signatures within populations, whereas the cross-population extended haplotype homozygosity (XP-EHH) and fixation index (FST ) approaches were preferred to reveal differently fixed genomic regions between native Turkish and cosmopolitan cattle breeds. Selection signatures in 142 genomic regions containing 305 genes were detected within eight cattle breeds by iHS statistics. The XP-EHH and FST approaches revealed that 197 and 114 SNPs were under selection pressure, respectively, which overlapped with 144 and 190 genes, respectively. A total of 18 genes were detected by at least two approaches. Six genes related to disease resistance (TTP2), meat yield (DIAPH3 and METTL21C), meat quality (ZNF24 and ZNF397) and first calving interval (ZSCAN30) turned out to be differently fixed between native Turkish and cosmopolitan cattle breeds, as they were identified by both XP-EHH and FST approaches. In addition, the iHS approach revealed that eight genes associated with visual modality (LSGN), olfaction (MOXD2, OR4C1F and OR4C1F), and immune response (TRBV3-1 and CLDN10) were under selection pressure in both native and cosmopolitan cattle breeds. Owing to their being significantly related to survival traits, these regions may have played a key role in cattle genome evolution. Future studies utilising denser genetic data are required to obtain deeper knowledge on effects of natural and artificial selection in Anatolian cattle breeds. © 2023 Stichting International Foundation for Animal Genetics.

Genome-wide association study for growth traits in Blanco Orejinegro and Romosinuano cattle.

Growth traits are economically important characteristics for the genetic improvement of local cattle breeds. Genome-wide association studies (GWAS) provide valuable information to enhance the understanding on the genetics of complex traits. The aim of this study was to perform a GWAS to identify genomic regions and genes associated to birth weight, weaning weight adjusted for 240 days, 16 months, and 24 months weight in Romosinuano (ROMO) and Blanco Orejinegro (BON) cattle. A single-step genomic-BLUP was implemented using 596 BON and 569 ROMO individuals that were genotyped with an Illumina BovineSNP50 BeadChip. There were 25 regions of interest identified on different chromosomes, with few of them simultaneously associated with two or more growth traits and some were common to both breeds. The gene mapping allowed to find 173 annotations on these regions, from which 49 represent potential candidate genes with known growth-related functions in cattle and other species. Among the regions that were associated with several growth traits, that at 24 - 27 MB of BTA14, has important candidate genes such as LYPLA1, XKR4, TMEM68 and PLAG1. Another region of interest at 0.40-0.77 Mb of BTA23 was identified in both breeds, containing KHDRBS2 as a potential candidate gene influencing body weight. Future studies targeting these regions could provide more knowledge to uncover the genetic architecture underlying growth traits in BON and ROMO cattle. The genomic regions and genes identified in this study could be used to improve the prediction of genetic merit for growth traits in these creole cattle breeds.

Genome-wide analysis emancipates genomic diversity and signature of selection in Altay white-headed cattle of Xinjiang, China.

Altay white-headed cattle have not received enough attention for several reasons. Due to irrational breeding and selection practices, the number of pure Altay white-headed cattle has decreased significantly and the breed is now on the eve of extinction. The genomic characterization will be a crucial step towards understanding the genetic basis of productivity and adaptability to survival under native Chinese agropastoral systems; nevertheless, no attempt has been made in Altay white-headed cattle. In the current study, we compared the genomes of 20 Altay white-headed cattle to the genomes of 144 individuals in representative breeds. Population genetic diversity revealed that the nucleotide diversity of Altay white-headed cattle was less than that of indicine breeds and comparable to that of Chinese taurus cattle. Using population structure analysis, we also found that Altay white-headed cattle carried the ancestry of the European and East Asian cattle lineage. In addition, we used three different methods (F ST, θπ ratio and XP-EHH) to investigate the adaptability and white-headed phenotype of Altay white-headed cattle and compared it with Bohai black cattle. We found EPB41L5, SCG5 and KIT genes on the list of the top one percent genes, these genes might have an association with environmental adaptability and the white-headed phenotype for this breed. Our research reveals the distinctive genomic features of Altay white-headed cattle at the genome-wide level.

Genetic diversity and population structure of a Peruvian cattle herd using SNP data.

New-generation sequencing technologies, among them SNP chips for massive genotyping, are useful for the effective management of genetic resources. To date, molecular studies in Peruvian cattle are still scarce. For the first time, the genetic diversity and population structure of a reproductive nucleus cattle herd of four commercial breeds from a Peruvian institution were determined. This nucleus comprises Brahman (N = 9), Braunvieh (N = 9), Gyr (N = 5), and Simmental (N = 15) breeds. Additionally, samples from a locally adapted creole cattle, the Arequipa Fighting Bull (AFB, N = 9), were incorporated. Female individuals were genotyped with the GGPBovine100K and males with the BovineHD. Quality control, and the proportion of polymorphic SNPs, minor allele frequency, expected heterozygosity, observed heterozygosity, and inbreeding coefficient were estimated for the five breeds. Admixture, principal component analysis (PCA), and discriminant analysis of principal components (DAPC) were performed. Also, a dendrogram was constructed using the Neighbor-Joining clustering algorithm. The genetic diversity indices in all breeds showed a high proportion of polymorphic SNPs, varying from 51.42% in Gyr to 97.58% in AFB. Also, AFB showed the highest expected heterozygosity estimate (0.41 ± 0.01), while Brahman the lowest (0.33 ± 0.01). Besides, Braunvieh possessed the highest observed heterozygosity (0.43 ± 0.01), while Brahman the lowest (0.37 ± 0.02), indicating that Brahman was less diverse. According to the molecular variance analysis, 75.71% of the variance occurs within individuals, whereas 24.29% occurs among populations. The pairwise genetic differentiation estimates (FST) between breeds showed values that ranged from 0.08 (Braunvieh vs. AFB) to 0.37 (Brahman vs. Braunvieh). Similarly, pairwise Reynold's distance ranged from 0.09 (Braunvieh vs. AFB) to 0.46 (Brahman vs. Braunvieh). The dendrogram, similar to the PCA, identified two groups, showing a clear separation between Bos indicus (Brahman and Gyr) and B. taurus breeds (Braunvieh, Simmental, and AFB). Simmental and Braunvieh grouped closely with the AFB cattle. Similar results were obtained for the population structure analysis with K = 2. The results from this study would contribute to the appropriate management, avoiding loss of genetic variability in these breeds and for future improvements in this nucleus. Additional work is needed to speed up the breeding process in the Peruvian cattle system.

Analysis of runs of homozygosity of cattle living in different climate zones.

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

Domestication and evolution of cattle originated different modern breeds in different places worldwide. The interaction between natural and artificial selection and the adaptation to environment shaped the genome, and the three different types of cattle here considered (taurus, indicus, and composite) may harbor different selection signatures. To study the difference among types and region of origin (tropics, temperate, and continental zones), Regions of Homozygosity (ROH) were used. ROH are continuous homozygous chromosomal segments identical by descendent, which characteristics can give information about inbreeding occurrence and natural and artificial selection. Moreover, it had been investigated which genes were mapped in these regions, and if interesting differences pertaining to environmental adaptation or fitness in general, could be observed.

Genomic rotational crossbreeding with advanced optimum contribution selection methods applied to simulated German Angler and German Holstein dairy cattle populations.

Many local dairy cattle breeds are facing genetic extinction due to a large proportion of foreign genes, which have been introgressed in the past. In addition, the performance gap to popular high-yielding breeds is increasing, resulting in a risk of numeric extinction. In the present simulation study, a genomic rotational crossbreeding scheme with the high-yielding German Holstein breed and the numerically small German Angler breed was analysed with the aim to utilize heterosis effects in the crossbred animals. Simultaneously inbreeding was controlled, and the amount of Holstein introgression observed in the Angler breed was reduced. Different scenarios of implementing OCS methods for Angler individuals were evaluated, which differed in their restrictions regarding kinship, native kinship, as well as the amount of genetic contributions from German Holstein. The results showed that rotational crossbreeding can result in superior crossbred offspring compared to the purebred parental lines, whereby OCS methods can simultaneously restrict the increase in inbreeding and keep the Holstein contributions at their current level. However, reducing the amount of migrant contributions while restricting the increase in the native kinship in Angler turned out to be a costly restriction. The reason was that Angler with low genetic contributions from Holsteins tended to have similar Angler ancestors. Consequently, reducing Holstein contributions would considerably increase the native kinship in Angler if it were not constrained. The constraint on the native kinship made a constraint on the conventional kinship superfluous and caused it to increase at a much lower rate than envisaged. This led to both, a high genetic diversity and a low genetic gain. The high genetic diversity in Angler also resulted in lower and oscillating heterosis effects in the crossbred animals. Thus, the reduction of migrant contribution did not increase heterosis effects in the crossbred offspring, and did not result in superior crossbred offspring in general.

Comparative analysis of buffalo, local and continental cattle carcasses with the European Union classification system in Greece.

The carcasses produced in Greece from Greek Buffalo, indigenous and local crossbred cattle were evaluated using the European Union (EU) bovine classification system and compared with carcasses of continental breeds reared in the country. The main carcass production (67.8%) came from animals born in Greece with a relatively low percentage from local breeds (17%). The lightest carcass weight was found in Greek Brachyceros [145.2 ± 2.3 kg (males); 146.0 ± 3.1 kg (females)], while among locals the Greek Blonde rendered the heaviest carcass weights [303.2 ± 2.9 kg (males); 265.3 ± 3.5 kg (females)]. Muscle conformation scoring was the lowest in Buffalo [4.1 ± 0.1 (males); 4.9 ± 0.0 (females)] and the highest for males of Charolais and Blonde d'Aquitaine (11.2 ± 0.0 and 10.7 ± 0.0, respectively). The meat-purposed breeds produced carcasses of higher quality than dual-purposed, local breeds and dairy Holstein-Friesian. With the discriminant analysis, the correct classification of carcasses varied from 53.4% (males) to 47.8% (females) with high assignment percentage of Greek Buffalo males (81.6%), Simmental males (80.7%) and Holstein-Friesian females (98.5%).

Comprehensive assessment of genetic diversity, structure, and relationship in four Japanese cattle breeds by Illumina 50 K SNP array analysis.

There are four unique cattle breeds in Japan: Japanese Black, Japanese Brown, Japanese Polled, and Japanese Shorthorn. The objective of this study was to comprehensively assess the genetic diversity, structure, relationship, and the degree of influence from foreign breeds (Angus, Simmental, Hanwoo, Shorthorn, Ayrshire, Brown Swiss, and Devon) in the Japanese cattle breeds using Illumina 50 K SNP array. In principal component analysis, each Japanese breed was separately clustered except for Japanese Shorthorn and Shorthorn. Japanese cattle breeds also showed different genetic components from each other at K ≥ 5 in population structure analysis. Japanese Shorthorn, on the other hand, had a very similar structure to Shorthorn at K ≤ 9, and Japanese Polled had a partially similar component with Angus at K = 3-7. Such close relationships were also observed in the phylogenetic tree. These findings imply that Japanese cattle breeds share genetic components with European cattle breeds to some extent while they have been almost differentiated. In population structure analysis, Japanese Black cattle shared little genetic component (3.5%) with European breeds. This is the first study to determine the extent to which European breeds impact Japanese breeds.

Genetic Diversity and Selective Signature in Dabieshan Cattle Revealed by Whole-Genome Resequencing.

Dabieshan cattle are a typical breed of southern Chinese cattle that have the characteristics of muscularity, excellent meat quality and tolerance to temperature and humidity. Based on 148 whole-genome data, our analysis disclosed the ancestry components of Dabieshan cattle with Chinese indicine (0.857) and East Asian taurine (0.139). The Dabieshan genome demonstrated a higher genomic diversity compared with the other eight populations, supported by the observed nucleotide diversity, linkage disequilibrium decay and runs of homozygosity. The candidate genes were detected by a selective sweep, which might relate to the fertility (GPX5, GPX6), feed efficiency (SLC2A5), immune response (IGLL1, BOLA-DQA2, BOLA-DQB), heat resistance (DnaJC1, DnaJC13, HSPA4), fat deposition (MLLT10) and the coat color (ASIP). We also identified the "East Asian taurine-like" segments in Dabieshan cattle, which might contribute to meat quality traits. The results revealed by the unique and valuable genomic data can build a foundation for the genetic improvement and conservation of genetic resources for indigenous cattle breeds.

Genetic diversity of DGAT1 gene linked to milk production in cattle populations of Ethiopia.

BACKGROUND: Diacylglycerol acyl-CoA acyltransferase 1 (DGAT1) has become a promising candidate gene for milk production traits because of its important role as a key enzyme in catalyzing the final step of triglyceride synthesis. Thus use of bovine DGAT1 gene as milk production markers in cattle is well established. However, there is no report on polymorphism of the DGAT1 gene in Ethiopian cattle breeds. The present study is the first comprehensive report on diversity, evolution, neutrality evaluation and genetic differentiation of DGAT1 gene in Ethiopian cattle population. The aim of this study was to characterize the genetic variability of exon 8 region of DGAT1 gene in Ethiopian cattle breeds. RESULTS: Analysis of the level of genetic variability at the population and sequence levels with genetic distance in the breeds considered revealed that studied breeds had 11, 0.615 and 0.010 haplotypes, haplotype diversity and nucleotide diversity respectively. Boran-Holstein showed low minor allele frequency and heterozygosity, while Horro showed low nucleotide and haplotype diversities. The studied cattle DGAT1 genes were under purifying selection. The neutrality test statistics in most populations were negative and statistically non-significant (p > 0.10) and consistent with a populations in genetic equilibrium or in expansion. Analysis for heterozygosity, polymorphic information content and inbreeding coefficient revealed sufficient genetic variation in DGAT1 gene. The pairwise FST values indicated significant differentiation among all the breeds (FST = 0.13; p ≤ 0.05), besides the rooting from the evolutionary or domestication history of the cattle inferred from the phylogenetic tree based on the neighbourhood joining method. There was four separated cluster among the studied cattle breeds, and they shared a common node from the constructed tree. CONCLUSION: The cattle populations studied were polymorphic for DGAT1 locus. The DGAT1 gene locus is extremely crucial and may provide baseline information for in-depth understanding, exploitation of milk gene variation and could be used as a marker in selection programmes to enhance the production potential and to accelerate the rate of genetic gain in Ethiopian cattle populations exposed to different agro ecology condition.

Seasonal survey, risk factor's analysis and genotyping of Theileria annulata infecting cattle in Punjab province, Pakistan.

Theileriosis is one of the most frequently reported tick borne diseases in tropical and subtropical regions and leads to annual economic losses, such as the reduced dairy products and increased casualties. Tropical theileriosis is caused by Theileria annulata and the present study was designed to improve our knowledge of Theileria annulata infection in Pakistani cattle. In order to assess the prevalence of Theileria annulata on cattle from Multan district in the Punjab province (Pakistan) according to seasons and other risk factors, a total of 1020 blood samples (340 samples each from cross, Holstein Frisian and Sahiwal breed) were collected between 2020 and 2022. Based on Tams-1 partial sequence amplification, the overall T. annulata prevalence was estimated at 11.3% (115/1020). The highest prevalence was observed in autumn season (14.1%), followed by winter (12.9%), summer (11.4%) and spring (6.7%) season. Sahiwal cattle were most susceptible to T. annulata infection (13.2%) followed by Crossbred (11.8%) and Holstein Frisian (8.8%). Epidemiological factor analysis revealed that female cattle, cattle rose with other dairy animals at farm, tick infested cattle, and cattle raised with dogs infested with ticks were associated with the prevalence of T. annulata. White blood cells, lymphocyte (%), Monocyte (%) hemoglobin, mean cell hemoglobin, mean corpuscular hemoglobin concentration, and platelet count were significantly affected blood parameters in T. annulata positive cattle of all three breeds. Representative partial Tams-1 sequences of four Pakistani T. annulata isolates revealed a single genotype genetically close to those infecting cattle from neighboring countries like Iran, Turkey and Egypt. Longitudinal survey and phylogenetic positioning of T. annulata is recommended for epidemiological correlation, diagnosis and treatment of theileriosis in such an agricultural region of Pakistan.

Selection and Drift: A Comparison between Historic and Recent Dutch Friesian Cattle and Recent Holstein Friesian Using WGS Data.

Over the last century, genetic diversity in many cattle breeds has been affected by the replacement of traditional local breeds with just a few milk-producing breeds. In the Netherlands, the local Dutch Friesian breed (DF) has gradually been replaced by the Holstein Friesian breed (HF). The objective of this study is to investigate genomewide genetic diversity between a group of historically and recently used DF bulls and a group of recently used HF bulls. Genetic material of 12 historic (hDF), 12 recent DF bulls (rDF), and 12 recent HF bulls (rHF) in the Netherlands was sequenced. Based on the genomic information, different parameters-e.g., allele frequencies, inbreeding coefficient, and runs of homozygosity (ROH)-were calculated. Our findings showed that a large amount of diversity is shared between the three groups, but each of them has a unique genetic identity (12% of the single nucleotide polymorphisms were group-specific). The rDF is slightly more diverged from rHF than hDF. The inbreeding coefficient based on runs of homozygosity (Froh) was higher for rDF (0.24) than for hDF (0.17) or rHF (0.13). Our results also displayed the presence of several genomic regions that differentiated between the groups. In addition, thirteen, forty-five, and six ROH islands were identified in hDF, rDF, and rHF, respectively. The genetic diversity of the DF breed reduced over time, but this did not lead to higher inbreeding levels-especially, inbreeding due to recent ancestors was not increased.

Naturally Occurring Glycosidases in Milk from Native Cattle Breeds: Activity and Consequences on Free and Protein Bound-Glycans.

Little is known about the extent of variation and activity of naturally occurring milk glycosidases and their potential to degrade milk glycans. A multi-omics approach was used to investigate the relationship between glycosidases and important bioactive compounds such as free oligosaccharides and O-linked glycans in bovine milk. Using 4-methylumbelliferone (4-MU) assays activities of eight indigenous glycosidases were determined, and by mass spectrometry and 1H NMR spectroscopy various substrates and metabolite products were quantified in a subset of milk samples from eight native North European cattle breeds. The results showed a clear variation in glycosidase activities among the native breeds. Interestingly, negative correlations between some glycosidases including ß-galactosidase, N-acetyl-ß-d-glucosaminidase, certain oligosaccharide isomers as well as O-linked glycans of κ-casein were revealed. Further, a positive correlation was found for free fucose content and α-fucosidase activity (r = 0.37, p-value < 0.001) indicating cleavage of fucosylated glycans in milk at room temperature. The results obtained suggest that milk glycosidases might partially degrade valuable glycans, which would result in lower recovery of glycans and thus represent a loss for the dairy ingredients industry if these activities are pronounced.

Access to ultra-long IgG CDRH3 bovine antibody sequences using short read sequencing technology.

The advances in high-throughput DNA sequencing and recombinant antibody technologies has presented new methods for characterizing antibody repertoires and significantly increased our understanding on the functional role of antibodies in immunity and their use in diagnostics, vaccine antigen design and as biological therapeutics. A subset of Bos taurus antibodies possesses unique ultra-long third complementary-determining region of the heavy chain (CDRH3) and are of special interest because they are thought to have unique functional abilities of broadly neutralizing properties - a functional role that has not been fully explored in vaccine development. Next generation sequencing technologies that are widely used to profile immunoglobulin (Ig) repertoires are based on short-read methods such as the Illumina technology. Although this technology has worked well in sequencing Ig V-D-J regions of most jawed vertebrates, it has faced serious technical challenges with sequencing regions in bovine Ig bearing ultra-long CDRH3 sequences, which are longer than 120 bp. To overcome this limitation, we have developed a sequencing strategy based on nested PCR products that allows sequence assembly of full-length bovine Ig heavy-chain (IgH) V-D-J regions. We have used this strategy to sequence IgH V-D-J regions of two Bos indicus breeds, Ankole and Boran. We confirm the presence of ultra-long CDRH3 sequences in IgG transcripts in both African cattle breeds, and provide preliminary evidence for differences and preferences in germline VH, DH and JH allele gene usage as well as differences in the length of the VH region in the two bovine breeds. Our method provides tools that should allow more robust analyses of ultra-long CDRH3 sequences aiding antibody and epitope discovery in different cattle breeds and their role in mediating immunity.

Genomic Prediction in Local Breeds: The Rendena Cattle as a Case Study.

The maintenance of local cattle breeds is key to selecting for efficient food production, landscape protection, and conservation of biodiversity and local cultural heritage. Rendena is an indigenous cattle breed from the alpine North-East of Italy, selected for dual purpose, but with lesser emphasis given to beef traits. In this situation, increasing accuracy for beef traits could prevent detrimental effects due to the antagonism with milk production. Our study assessed the impact of genomic information on estimated breeding values (EBVs) in Rendena performance-tested bulls. Traits considered were average daily gain, in vivo EUROP score, and in vivo estimate of dressing percentage. The final dataset contained 1691 individuals with phenotypes and 8372 animals in pedigree, 1743 of which were genotyped. Using the cross-validation method, three models were compared: (i) Pedigree-BLUP (PBLUP); (ii) single-step GBLUP (ssGBLUP), and (iii) weighted single-step GBLUP (WssGBLUP). Models including genomic information presented higher accuracy, especially WssGBLUP. However, the model with the best overall properties was the ssGBLUP, showing higher accuracy than PBLUP and optimal values of bias and dispersion parameters. Our study demonstrated that integrating phenotypes for beef traits with genomic data can be helpful to estimate EBVs, even in a small local breed.

Harmful mutation load in the mitochondrial genomes of cattle breeds.

OBJECTIVE: Domestication of wild animals results in a reduction in the effective population size, and this could affect the deleterious mutation load of domesticated breeds. Furthermore, artificial selection will also contribute to the accumulation of deleterious mutations due to the increased rate of inbreeding among these animals. The process of domestication, founder population size, and artificial selection differ between cattle breeds, which could lead to a variation in their deleterious mutation loads. We investigated this using mitochondrial genome data from 364 animals belonging to 18 cattle breeds of the world. RESULTS: Our analysis revealed more than a fivefold difference in the deleterious mutation load among cattle breeds. We also observed a negative correlation between the breed age and the proportion of deleterious amino acid-changing polymorphisms. This suggests a proportionally higher deleterious SNPs in young breeds compared to older breeds. Our results highlight the magnitude of difference in the deleterious mutations present in the mitochondrial genomes of various breeds. The results of this study could be useful in predicting the rate of incidence of genetic diseases in different breeds.

Genetic assessment of BoLA-DRB3 polymorphisms by comparing Bangladesh, Ethiopian, and Korean cattle.

Attributable to their major function in pathogen recognition, the use of bovine leukocyte antigens (BoLA) as disease markers in immunological traits in cattle is well established. However, limited report exists on polymorphism of the BoLA gene in zebu cattle breeds by high resolution typing methods. Thus, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 100 animals (Boran, n = 13; Sheko, n = 20; Fogera, n = 16; Horro, n = 19), Hanwoo cattle (n = 18) and Bangladesh Red Chittagong zebu (n = 14). Out of the 59 detected alleles, 43 were already deposited under the Immuno Polymorphism Database for major histocompatibility complex (IPD-MHC) while 16 were unique to this study. Assessment of the level of genetic variability at the population and sequence levels with genetic distance in the breeds considered in this study showed that Zebu breeds had a gene diversity score greater than 0.752, nucleotide diversity score greater than 0.152, and mean number of pairwise differences higher than 14, being very comparable to those investigated for other cattle breeds. Regarding neutrality tests analyzed, we investigated that all the breeds except Hanwoo had an excess number of alleles and could be expected from a recent population expansion or genetic hitchhiking. Howbeit, the observed heterozygosity was not significantly (p < 0.05) higher than the expected heterozygosity. The Hardy Weinberg equilibrium (HWE) analysis revealed non-significant excess of heterozygote animals, indicative of plausible over-dominant selection. The pairwise FST values suggested a low genetic variation among all the breeds (FST = 0.056; p < 0.05), besides the rooting from the evolutionary or domestication history of the cattle. No detached clade was observed in the evolutionary divergence study of the BoLA-DRB3 gene, inferred from the phylogenetic tree based on the maximum likelihood model. The investigation herein indicated the clear differences in BoLA-DRB3 gene variability between African and Asian cattle breeds.

Genome-Wide Analysis Reveals Selection Signatures Involved in Meat Traits and Local Adaptation in Semi-Feral Maremmana Cattle.

The Maremmana cattle is an ancient Podolian-derived Italian breed raised in semi-wild conditions with distinctive morphological and adaptive traits. The aim of this study was to detect potential selection signatures in Maremmana using medium-density single nucleotide polymorphism array. Putative selection signatures were investigated combining three statistical approaches designed to quantify the excess of haplotype homozygosity either within (integrated haplotype score, iHS) or among pairs of populations (Rsb and XP-EHH), and contrasting the Maremmana with a single reference population composed of a pool of seven Podolian-derived Italian breeds. Overall, the three haplotype-based analyses revealed selection signatures distributed over 19 genomic regions. Of these, six relevant candidate regions were identified by at least two approaches. We found genomic signatures of selective sweeps spanning genes related to mitochondrial function, muscle development, growth, and meat traits (SCIN, THSD7A, ETV1, UCHL1, and MYOD1), which reflects the different breeding schemes between Maremmana (semi-wild conditions) and the other Podolian-derived Italian breeds (semi-extensive). We also identified several genes linked to Maremmana adaptation to the environment of the western-central part of Italy, known to be hyperendemic for malaria and other tick-borne diseases. These include several chemokine (C-C motif) ligand genes crucially involved in both innate and adaptive immune responses to intracellular parasite infections and other genes playing key roles in pulmonary disease (HEATR9, MMP28, and ASIC2) or strongly associated with malaria resistance/susceptibility (AP2B1). Our results provide a glimpse into diverse selection signatures in Maremmana cattle and can be used to enhance our understanding of the genomic basis of environmental adaptation in cattle.