Assessing the genetic background and genomic relatedness of red cattle populations originating from Northern Europe

Unraveling inbreeding patterns and selection signals in Alpine Grey cattle

Inbreeding plays a crucial role in livestock breeding, influencing genetic diversity and phenotypic traits. Genomic data have helped address limitations posed by incomplete pedigrees, providing deeper insights into breed genetic diversity. This study assesses inbreeding levels via pedigree and genomic approaches and analyzes old and recent inbreeding using runs of homozygosity (ROH), and selection signals in Alpine Grey cattle. Pedigree data from 165 575 individuals, analyzed with INBUPGF90 software, computed inbreeding coefficients. Genomic-based coefficients derived from PLINK v1.9. or DetectRUNS R package analyses of 1 180 individuals' genotypes. Common single nucleotide polymorphisms within ROH pinpointed genomic regions, aggregating into "ROH islands" indicative of selection pressure. Overlaps with USCS Genome Browser unveiled gene presence. Moderate correlations (0.20-0.54) existed between pedigree and genomic coefficients, with most genomic estimators having higher (>0.8) correlation values. Inbreeding averaged 0.04 in < 8 Mb ROH segments, and 0.03 in > 16 Mb segments; > 90% of ROHs were < 8 Mb, indicating ancient inbreeding prevalence. Recent inbreeding proved less detrimental than in cosmopolitan breeds. Two major ROH islands on chromosomes 6 and 7 harbored genes linked to immune response, disease resistance (PYURF, HERC3), and fertility (EIF4EBP3, SRA1). This study underscores the need for detailed inbreeding analyses to understand genetic characteristics and historical changes in local breeds like Alpine Grey cattle. Genomic insights, especially from ROH, facilitated overcoming pedigree limitations, illuminating breed genetic diversity. Our findings reveal ancient inbreeding's enduring genetic impact and ROH islands potential for selective sweeps, elucidating traits in Alpine Grey cattle.

Genetic Variation and Composition of Two Commercial Estonian Dairy Cattle Breeds Assessed by SNP Data

The aims of this study were to assess the genomic relatedness of Estonian and selected European dairy cattle breeds and to examine the within-breed diversity of two Estonian dairy breeds using genome-wide SNP data. This study was based on a genotyped heifer population of the Estonian Red (ER) and Estonian Holstein (EH) breeds, including about 10% of all female cattle born in 2017-2020 (sample sizes n = 215 and n = 2265, respectively). The within-breed variation study focused on the level of inbreeding using the ROH-based inbreeding coefficient. The genomic relatedness analyses were carried out among two Estonian and nine European breeds from the WIDDE database. Admixture analysis revealed the heterogeneity of ER cattle with a mixed pattern showing several ancestral populations containing a relatively low proportion (1.5-37.0%) of each of the reference populations used. There was a higher FROH in EH (FROH = 0.115) than in ER (FROH = 0.044). Compared to ER, the long ROHs of EH indicated more closely related parents. The paternal origin of the genetic material used in breeding had a low effect on the inbreeding level. However, among EH, the highest genomic inbreeding was estimated in daughters of USA-born sires.

Exploring the Effects of Robertsonian Translocation 1/29 (Rob (1;29)) on Genetic Diversity in Minor Breeds of Spanish Berrenda Cattle via Genome-Wide Analysis

Most of the previous studies on the genetic variability in Spanish "Berrenda" breeds have been carried out using DNA microsatellites. The present work aimed to estimate the genetic diversity, population structure, and potential genetic differences among individuals of both Berrenda breeds and groups based on the presence of the Robertsonian chromosomal translocation, rob (1;29). A total of 373 samples from animals belonging to the two breeds, including 169 cases diagnosed as rob (1;29)-positive, were genotyped using an SNP50K chip. The genetic diversity at the breed level did not show significant differences, but it was significantly lower in those subpopulations containing the rob (1;29). Runs of homozygosity identified a region of homozygosity on chromosome 6, where the KIT (KIT proto-oncogene, receptor tyrosine kinase) gene, which determines the typical spotted coat pattern in both breeds, is located. The four subpopulations considered showed minor genetic differences. The regions of the genome that most determined the differences between the breeds were observed on chromosomes 4, 6, 18, and 22. The presence of this Robertsonian translocation did not result in sub-structuring within each of the breeds considered. To improve the reproductive performance of Berrenda breeds, it would be necessary to implement strategies considering the involvement of potential breeding stock carrying rob (1;29).

Genetic diversity and positive signatures of selection in indigenous cattle breeds of Iran

Animal domestication, climate changes over time, and artificial selection have played significant roles in shaping the genome structure of various animal species, including cattle. These processes have led to the emergence of several indigenous cattle breeds with distinct genetic characteristics. This study focused on unraveling the genetic diversity and identifying candidate genomic regions in eight indigenous cattle breeds of Iran. The data consisted of ∼777 962 single nucleotide polymorphisms (SNPs) of 89 animals from Iranian indigenous cattle scattered throughout the country. We employed various methods, including integrated haplotype score, FST, and cross-population composite likelihood ratio, to conduct a genome scan for detecting selection signals within and between cattle populations. Average observed heterozygosity across the populations was 0.36, with a range of 0.32-0.40. In addition, negative and low rates of inbreeding (FIS) in the populations were observed. The genome-wide analysis revealed several genomic regions that harbored candidate genes associated with production traits (e.g., MFSD1, TYW5, ADRB2, BLK, and CRTC3), adaptation to local environmental constraints (CACNA2D1, CXCL3, and GRO1), and coat color (DYM). Finally, the study of the reported quantitative trait loci (QTL) regions in the cattle genome demonstrated that the identified regions were associated with QTL related to important traits such as milk composition, body weight, daily gain, feed conversion, and residual feed intake. Overall, this study contributes to a better understanding of the genetic diversity and potential candidate genes underlying important traits in Iranian indigenous cattle breeds, which can inform future breeding and conservation efforts.

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.

Estimation of inbreeding, between-breed genomic relatedness and definition of sub-populations in red-pied cattle breeds

Currently, enhancing the collaboration between related breeds is of main importance to increase the competitivity and the sustainability of local breeds. One type of collaboration is the development of an across-breed reference population that will allow a better management of local breeds. For this purpose, the genomic relatedness between the local target breed and possible breeds to be included in the reference population should be estimated. In Europe, there are several local red-pied cattle breeds that would benefit from this kind of collaboration. However, how different red-pied cattle breeds from the Benelux are related to each other and can collaborate is still unclear. The objectives of this study were therefore: (1) to estimate the level of inbreeding of the East Belgian Red and White (EBRW), the Red-Pied of the Ösling (RPO) and Dutch red-pied cattle breeds; (2) to determine the genomic relatedness of several red-pied cattle breeds, with a special focus on two endangered breeds: the EBRW and the RPO, and (3) based on the second objective, to detect animals from other breeds that were genomically close enough to be considered as advantageous in the creation of an across-breed reference population of EBRW or RPO. The estimated inbreeding levels based on runs of homozygosity were relatively low for almost all the studied breeds and especially for the EBRW and RPO. This would imply that inbreeding is currently not an issue in these two endangered breeds and that their sustainability is not threatened by their level of inbreeding. The results from the principal component analysis, the phylogenetic tree and the clustering all highlighted that the EBRW and RPO breeds were included in the genomic continuum of the studied red-pied cattle breeds and can be therefore considered as genomically close to Dutch red-pied cattle breeds, highlighting the possibility of a collaboration between these breeds. Especially, EBRW animals were closely related to Deep Red and Improved Red animals while, to a lesser extent, the RPO animals were closely related to the Meuse-Rhine-Yssel breed. Based on these results, we could use distance measures, based either on the principal component analysis or clustering, to detect animals from Dutch breeds that were genomically closest to the EBRW or RPO breeds. This will finally allow the building of an across-breed reference population for EBRW or RPO for further genomic evaluations, considering these genomically closest animals from other breeds.

Genomic diversity and relationship analyses of endangered German Black Pied cattle (DSN) to 68 other taurine breeds based on whole-genome sequencing

German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (F_ST), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (F_Hom) and genomic inbreeding (F_RoH) through runs of homozygosity (RoHs). Region-wide F_ST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The F_Hom and F_RoH values in DSN were among the lowest. Regions with high F_ST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. F_ST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN.

Gene Expression of Aquaporins (AQPs) in Cumulus Oocytes Complex and Embryo of Cattle

Aquaporins (AQPs) are proteins with various functions related to proper cell function and early development in mammals. The aim of this study was to evaluate the presence of AQPs and determine their mRNA levels in the cumulus oocyte complex (COC) of four bovine breeds and in blastocysts of five bovine crosses. Grade I, II and III COCs were collected by ovum pick up from non-lactating heifers of the Brahaman, Holstein, Gir and Romosinuano breeds. Embryos were produced in vitro up to the blastocyst stage of the bovine ♀Gir × ♂Holstein, ♀Holstein × ♂Gir, ♀Brahman × ♂Holstein, ♀Holstein × ♂Brahman, and ♀Romosinuano × ♂Holstein crosses. mRNA expression of AQP1-AQP12b was estimated in COC and embryos by real-time-PCR. The presence of the twelve AQPs in the COCs and bovine embryos was established. Additionally, significant differences were determined in the expression of AQP6 and AQP12b in COCs, as well as in transcripts levels of AQP4, AQP8 and AQP9 from bovine embryos. Gene expression of AQPs in COCs and bovine embryos is consistent with the previously described biological functions. This is the first report of AQPs in COC of Gir, Brahman, Holstein and Romosinuano and embryos of five crossbreeds between Bos indicus and B. taurus.

Opportunities of Genomics for the Use of Semen Cryo-Conserved in Gene Banks

Shortly after the introduction of cryo-conserved semen in the main farm animal species, gene banks were founded. Safeguarding farm animal genetic diversity for future use was and is the main objective. A sampling of sires was based on their pedigree and phenotypic information. Nowadays, DNA information from cryo-conserved sires and from animals in the living populations has become available. The combination of their DNA information can be used to realize three opportunities: 1) to make the gene bank a more complete archive of genetic diversity, 2) to determine the history of the genetic diversity from the living populations, and 3) to improve the performance and genetic diversity of living populations. These three opportunities for the use of gene bank sires in the genomic era are outlined in this study, and relevant recent literature is summarized to illustrate the great value of a gene bank as an archive of genetic diversity.

Genetic Differentiation among Livestock Breeds—Values for Fst

Simple Summary

The degree of relationship among livestock breeds can be quantified by the F_st statistic, which measures the extent of genetic differentiation between them. An F_st value of 0.1 has often been taken as indicating that two breeds are indeed genetically distinct, but this concept has not been evaluated critically. Here, F_st values have been collated for the six major livestock species: cattle, sheep, goats, pigs, horses, and chickens. These values are remarkably variable both within and between species, demonstrating that F_st > 0.1 is not a reliable criterion for breed distinctiveness. However, the large body of F_st data accumulated in the last 20–30 years represents an untapped database that could contribute to the development of interdisciplinary research involving livestock breeds.

Abstract

(1) Background: The F_st statistic is widely used to characterize between-breed relationships. F_st = 0.1 has frequently been taken as indicating genetic distinctiveness between breeds. This study investigates whether this is justified. (2) Methods: A database was created of 35,080 breed pairs and their corresponding F_st values, deduced from microsatellite and SNP studies covering cattle, sheep, goats, pigs, horses, and chickens. Overall, 6560 (19%) of breed pairs were between breeds located in the same country, 7395 (21%) between breeds of different countries within the same region, 20,563 (59%) between breeds located far apart, and 562 (1%) between a breed and the supposed wild ancestor of the species. (3) Results: General values for between-breed F_st were as follows, cattle: microsatellite 0.06–0.12, SNP 0.08–0.15; sheep: microsatellite 0.06–0.10, SNP 0.06–0.17; horses: microsatellite 0.04–0.11, SNP 0.08–0.12; goats: microsatellite 0.04–0.14, SNP 0.08–0.16; pigs: microsatellite 0.06–0.27, SNP 0.15–0.22; chickens: microsatellite 0.05–0.28, SNP 0.08–0.26. (4) Conclusions: (1) Large amounts of F_st data are available for a substantial proportion of the world’s livestock breeds, (2) the value for between-breed F_st of 0.1 is not appropriate owing to its considerable variability, and (3) accumulated F_st data may have value for interdisciplinary research.

Improving the Accuracy of Multi-Breed Prediction in Admixed Populations by Accounting for the Breed Origin of Haplotype Segments

Numerically small breeds have often been upgraded with mainstream breeds. This historic introgression predisposes the breeds for joint genomic evaluations with mainstream breeds. The linkage disequilibrium structure differs between breeds. The marker effects of a haplotype segment may, therefore, depend on the breed from which the haplotype segment originates. An appropriate method for genomic evaluation would account for this dependency. This study proposes a method for the computation of genomic breeding values for small admixed breeds that incorporate phenotypic and genomic information from large introgressed breeds by considering the breed origin of alleles (BOA) in the evaluation. The proposed BOA model classifies haplotype segments according to their origins and assumes different but correlated SNP effects for the different origins. The BOA model was compared in a simulation study to conventional within-breed genomic best linear unbiased prediction (GBLUP) and conventional multi-breed GBLUP models. The BOA model outperformed within-breed GBLUP as well as multi-breed GBLUP in most cases.

Genomics Confirm an Alarming Status of the Genetic Diversity of Belgian Red and Belgian White Red Cattle

Simple Summary

Genetic diversity of livestock is vitally important as it enables the adaptation of future populations to changing environments. Therefore, preserving a sufficiently large genetic diversity is key. However, for many local cattle populations, little is known about their genetic diversity such as inbreeding level, effective size etc. We studied the genetic diversity of two local Belgian red cattle populations (Belgian Red and Belgian White Red cattle) using state-of-the-art genomic techniques. These tools assessed diversity at the population and individual level, and allowed the positioning of these two breeds in an international context of 52 other (European) cattle breeds. Accordingly, we contribute to the general knowledge of European red cattle, and more specifically we help the breeders, breed organization and the government to manage the genetic diversity of both breeds.

Abstract

Genetic diversity is increasingly important for researchers and society. Small and local populations deserve more attention especially, as they may harbor important characteristics. Moreover, small populations are at greater risk and their genetic management is often more challenging. Likewise, European red cattle populations are threatened, as they are outcompeted by more specialized cattle breeds. In this study, we investigate the genetic diversity of two local Belgian red cattle breeds: Belgian Red and Belgian White Red cattle. A total of 270 animals were genotyped via medium density SNP arrays. Genetic diversity was assessed using runs of homozygosity screening, effective population size estimation and F_st analyses. Genomic inbreeding coefficients based on runs of homozygosity were estimated at 7.0% for Belgian Red and 6.1% for Belgian White Red cattle, and both populations had a low effective population size (68 and 86, respectively). PCA, F_st and admixture analyses revealed the relationship to 52 other international breeds, where they were closest related to some Belgian, French, Scandinavian and Dutch breeds. Moreover, F_st analyses revealed for Belgian Red cattle a signature of selection on BTA6, adjacent to the KIT gene. This study gains important knowledge on the genetic diversity of these two small local red cattle breeds, and will aid in their (genetic) management.

Identification of Reference Genes for Expression Studies in the Whole-Blood from Three Cattle Breeds under Two States of Livestock Weather Safety

Simple Summary

Reductions in the fertility, body weight, and growth rate of cattle across the world are associated with the global warming phenomenon. Developing optimal management strategies is an important aspect of breeding programs for different breeds. Blood tissue undergoes dramatic physiological and metabolic changes during heat stress conditions, which involves the expression and regulation of a great number of genes across species. Real-time quantitative PCR (qPCR) is a method for the rapid and reliable quantification of mRNA transcription. Reference genes are used to normalize mRNA levels between different samples. Thus, the selection of high-quality reference genes is necessary for the interpretation of data generated by real-time PCR.

Abstract

Real-time PCR is widely used to study the relative abundance of mRNA due to its specificity, sensitivity, and repeatability quantification. However, relative quantification requires a reference gene, which should be stable in its expression, showing lower variation by experimental conditions or tissues. The aim of this study was to evaluate the stability of the expression of five commonly used reference genes (actb, ywhaz, b2m, sdha, and 18s rRNA) at different physiological stages (alert and emergency) in three different cattle breeds. In this study, five genes (actb, ywhaz, b2m, sdha, and 18s rRNA) were selected as candidate reference genes for expression studies in the whole blood from three cattle breeds (Romosinuano, Gyr, and Brahman) under heat stress conditions. The transcription stability of the candidate reference genes was evaluated using geNorm and NormFinder. The results showed that actb, 18SrRNA, and b2m expression were the most stable reference genes for whole blood of Gyr and Brahman breeds under two states of livestock weather safety (alert and emergency). Meanwhile, actb, b2m, and ywhaz were the most stable reference genes for the Romosinuano breed.

Genome-wide SNP analysis clearly distinguished the Belarusian Red cattle from other European cattle breeds

Local breeds can serve as an important source of genetic variability in domestic animal species. This study aimed to assess the genetic diversity and population structure of Belarusian Red cattle and their differentiation from other European cattle populations based on genome-wide SNP genotypes. Twenty pedigree-recorded non-closely related cows of Belarusian Red cattle were genotyped using the Illumina BovineHD BeadChip. Genotypes of 22 other European cattle breeds were included in the study for comparison. A total of 28 562 SNPs passed through the quality control checks and were selected for analysis. The Belarusian Red cattle displayed a moderate level of genetic variability (_U H_E  = 0.341, H_O  = 0.368), and the highest heterozygote excess (_U F_IS  = -0.066), among the studied breeds; this reflects the contribution of multiple breeds to their formation. The principal component analysis, F_ST -based Neighbor-Net tree and Admixture clustering, clearly distinguished the Belarusian Red cattle from the other European cattle breeds. Moreover, the presence of ancestral genomic components of Danish Red and Brown Swiss breeds were clearly visible, which agrees with the breed's history and its recent development. Our study highlights the importance of maintaining the specific genomic components, which makes a significant contribution to the global genetic diversity in the modern population of Belarusian Red cattle, allowing us to consider them a valuable national genetic resource. Our research results will be useful for the development of conservation programs for this local cattle breed.