Improving selection decisions with mating information by accounting for Mendelian sampling variances looking two generations ahead.

BACKGROUND: Breeding programs are judged by the genetic level of animals that are used to disseminate genetic progress. These animals are typically the best ones of the population. To maximise the genetic level of very good animals in the next generation, parents that are more likely to produce top performing offspring need to be selected. The ability of individuals to produce high-performing progeny differs because of differences in their breeding values and gametic variances. Differences in gametic variances among individuals are caused by differences in heterozygosity and linkage. The use of the gametic Mendelian sampling variance has been proposed before, for use in the usefulness criterion or Index5, and in this work, we extend existing approaches by not only considering the gametic Mendelian sampling variance of individuals, but also of their potential offspring. Thus, the criteria developed in this study plan one additional generation ahead. For simplicity, we assumed that the true quantitative trait loci (QTL) effects, genetic map and the haplotypes of all animals are known. RESULTS: In this study, we propose a new selection criterion, ExpBVSelGrOff, which describes the genetic level of selected grand-offspring that are produced by selected offspring of a particular mating. We compare our criterion with other published criteria in a stochastic simulation of an ongoing breeding program for 21 generations for proof of concept. ExpBVSelGrOff performed better than all other tested criteria, like the usefulness criterion or Index5 which have been proposed in the literature, without compromising short-term gains. After only five generations, when selection is strong (1%), selection based on ExpBVSelGrOff achieved 5.8% more commercial genetic gain and retained 25% more genetic variance without compromising inbreeding rate compared to selection based only on breeding values. CONCLUSIONS: Our proposed selection criterion offers a new tool to accelerate genetic progress for contemporary genomic breeding programs. It retains more genetic variance than previously published criteria that plan less far ahead. Considering future gametic Mendelian sampling variances in the selection process also seems promising for maintaining more genetic variance.

Comparison of genomic prediction methods for early growth traits of Inner Mongolia cashmere goats based on multi trait models.

Inner Mongolia cashmere goat is an excellent livestock breed formed through long-term natural selection and artificial breeding, and is currently a world-class dual-purpose breed producing cashmere and meat. Multi trait animal model is considered to significantly improve the accuracy of genetic evaluation in livestock and poultry, enabling indirect selection between traits. In this study, the pedigree, genotype, environment, and phenotypic records of early growth traits of Inner Mongolia cashmere goats were used to build multi trait animal model., Then three methods including ABLUP, GBLUP, and ssGBLUP wereused to estimate the genetic parameters and genomic breeding values of early growth traits (birth weight, weaning weight, average daily weight gain before weaning, and yearling weight). The accuracy and reliability of genomic estimated breeding value are further evaluated using the five fold cross validation method. The results showed that the heritability of birth weight estimated by three methods was 0.13-0.15, the heritability of weaning weight was 0.13-0.20, heritability of daily weight gain before weaning was 0.11-0.14, and the heritability of yearling weight was 0.09-0.14, all of which belonged to moderate to low heritability. There is a strong positive genetic correlation between weaning weight and daily weight gain before weaning, daily weight gain before weaning and yearling weight, with correlation coefficients of 0.77-0.79 and 0.56-0.67, respectively. The same pattern was found in phenotype correlation among traits. The accuracy of the estimated breeding values by ABLUP, GBLUP, and ssGBLUP methods for birth weight is 0.5047, 0.6694, and 0.7156, respectively; the weaning weight is 0.6207, 0.6456, and 0.7254, respectively; the daily weight gain before weaning was 0.6110, 0.6855, and 0.7357 respectively; and the yearling weight was 0.6209, 0.7155, and 0.7756, respectively. In summary, the early growth traits of Inner Mongolia cashmere goats belong to moderate to low heritability, and the speed of genetic improvement is relatively slow. The genetic improvement of other growth traits can be achieved through the selection of weaning weight. The ssGBLUP method has the highest accuracy and reliability in estimating genomic breeding value of early growth traits in Inner Mongolia cashmere goats, and is significantly higher than that from ABLUP method, indicating that it is the best method for genomic breeding of early growth weight in Inner Mongolia cashmere goats.

Deciphering genetic diversity in conserved cattle bulls to achieve sustainable development goals.

The primary objective of Sustainable Development Goal target 2.5 established by the United Nations is to ensure the preservation of genetic diversity in domesticated animals. The ICAR-National Bureau of Animal Genetic Resources in India has been actively engaged in the conservation of cattle and buffalo bull semen for long-term storage. This present study aimed to assess the genetic diversity present in the conserved cattle bull semen, which would aid in determining the most suitable strategy for future conservation management. A total of 192 bull semen belonging to 19 cattle breeds were selected to evaluate genetic diversity using 17 pairs of FAO recommended microsatellite primers. Total 267 alleles were detected across all the samples which indicates substantial amount of allelic variation is being maintained in conserved bulls. Further, all cattle bulls semen conserved showed higher observed heterozygosity than expected heterozygosity which indicates excess genetic diversity in all the populations. The FST, F IT and FIS value across the loci and population is 0.146 ± 0.009, 0.054 ± 0.038, and - 0.105 ± 0.035, respectively, which suggests lack of inbreeding in conserved cattle bull semen. This study has established genetic diversity in conserved cattle semen samples to achieve sustainable development goals. In addition, it provides compelling evidence that the current approach for conserving cattle bull semen is heading in the correct direction.

Development and verification of a 10K liquid chip for Hainan black goat based on genotyping by pinpoint sequencing of liquid captured targets.

BACKGROUND: China has thousands years of goat breeding and abundant goat genetic resources. Additionally, the Hainan black goat is one of the high-quality local goat breeds in China. In order to conserve the germplasm resources of the Hainan black goat, facilitate its genetic improvement and further protect the genetic diversity of goats, it is urgent to develop a single nucleotide polymorphism (SNP) chip for Hainan black goat. RESULTS: In this study, we aimed to design a 10K liquid chip for Hainan black goat based on genotyping by pinpoint sequencing of liquid captured targets (cGPS). A total of 45,588 candidate SNP sites were obtained, 10,677 of which representative SNP sites were selected to design probes, which finally covered 9,993 intervals and formed a 10K cGPS liquid chip for Hainan black goat. To verify the 10K cGPS liquid chip, some southern Chinese goat breeds and a sheep breed with similar phenotype to the Hainan black goat were selected. A total of 104 samples were used to verify the clustering ability of the 10K cGPS liquid chip for Hainan black goat. The results showed that the detection rate of sites was 97.34% -99.93%. 84.5% of SNP sites were polymorphic. The heterozygosity rate was 3.08%-36.80%. The depth of more than 99.4% sites was above 10X. The repetition rate was 99.66%-99.82%. The average consistency between cGPS liquid chip results and resequencing results was 85.58%. In addition, the phylogenetic tree clustering analysis verified that the SNP sites on the chip had better clustering ability. CONCLUSION: These results indicate that we have successfully realized the development and verification of the 10K cGPS liquid chip for Hainan black goat, which provides a useful tool for the genome analysis of Hainan black goat. Moreover, the 10K cGPS liquid chip is conducive to the research and protection of Hainan black goat germplasm resources and lays a solid foundation for its subsequent breeding work.

Lifetime milk production of Holstein cattle in the humid tropics compared to Holstein-Gyr and Holstein-Brahman crosses.

Crossbred cattle are commonly used for milk production in the tropics, combining the potential benefits of pure breeds with the heterosis effects of the offspring. However, no comprehensive assessment of lifetime productivity for crossbred versus purebred cattle in low-altitude tropical environments has been carried out. The present study compares the lifetime productivity of purebred Holstein (HO, n = 17,269), Gyr (GY4, n = 435), and Brahman (BR4, n = 622) with crossbreds Gyr × Holstein (GY × HO, n = 5521) and Brahman×Holstein (BR × HO, n = 5429) cows from dairy farms located in low and medium altitude tropical regions in Costa Rica. The production traits of interest were age at first calving (AFC), days open (DO), milk production per lactation (TMP), lactation length (LLEN), age at culling (ACUL), and number of lactations (NLAC). Estimates of heterosis were also calculated. The AFC for GY × HO crosses (33-34 months) was not significantly different (p > .05) from HO (33.8 months). For BR × HO crosses, a significant (p < .05) decrease in AFC (BR3HO1 35.6 months, BR2HO2 34.5 months, and BR1H03 33.3 months) was observed as the fraction of HO breed increased. Estimates of heterosis for AFC were favourable for both crosses, of a magnitude close to 3%. The DO for F1 crosses (GY2HO2 94 days; BR2HO2 96 days) was significantly (p < .05) lower than HO (123 days). Estimates of heterosis for DO were also favourable and above 15% for both crosses. The TMP and LLEN were higher for HO (TMP = 5003 kg; LLEN = 324 days) compared with GY × HO (TMP = 4428 to 4773 kg; LLEN = 298 to 312 days) and BR × HO (TMP = 3950 to 4761 kg; LLEN = 273 to 313 days) crosses. Heterosis for TMP was favourable but low for both crosses, with a magnitude below 3.0%. The NLAC for HO (4.6 lactations) was significantly (p < .05) lower than F1 (GY2HO2, 5.8 lactations; BR2HO2, 5.4 lactations). Heterosis for NLAC was above 6.0% for both crosses. Overall, estimates of lifetime income over feed costs per cow on average were USD 2637 (30.3%) and USD 734 (8.4%) higher in F1 GY × HO and BR × HO, respectively, compared to HO. In conclusion, crossbred animals, specifically those with Gyr and Brahman genetics, extend the productive lifespan, increasing economic returns.

Equivalence of variance components between standard and recursive genetic models using LDL' transformations.

BACKGROUND: Recursive models are a category of structural equation models that propose a causal relationship between traits. These models are more parameterized than multiple trait models, and they require imposing restrictions on the parameter space to ensure statistical identification. Nevertheless, in certain situations, the likelihood of recursive models and multiple trait models are equivalent. Consequently, the estimates of variance components derived from the multiple trait mixed model can be converted into estimates under several recursive models through LDL' or block-LDL' transformations. RESULTS: The procedure was employed on a dataset comprising five traits (birth weight-BW, weight at 90 days-W90, weight at 210 days-W210, cold carcass weight-CCW and conformation-CON) from the Pirenaica beef cattle breed. These phenotypic records were unequally distributed among 149,029 individuals and had a high percentage of missing data. The pedigree used consisted of 343,753 individuals. A Bayesian approach involving a multiple-trait mixed model was applied using a Gibbs sampler. The variance components obtained at each iteration of the Gibbs sampler were subsequently used to estimate the variance components within three distinct recursive models. CONCLUSIONS: The LDL' or block-LDL' transformations applied to the variance component estimates achieved from a multiple trait mixed model enabled inference across multiple sets of recursive models, with the sole prerequisite of being likelihood equivalent. Furthermore, the aforementioned transformations simplify the handling of missing data when conducting inference within the realm of recursive models.

Estimating genomic relationships of metafounders across and within breeds using maximum likelihood, pseudo-expectation-maximization maximum likelihood and increase of relationships.

BACKGROUND: The theory of "metafounders" proposes a unified framework for relationships across base populations within breeds (e.g. unknown parent groups), and base populations across breeds (crosses) together with a sensible compatibility with genomic relationships. Considering metafounders might be advantageous in pedigree best linear unbiased prediction (BLUP) or single-step genomic BLUP. Existing methods to estimate relationships across metafounders Γ are not well adapted to highly unbalanced data, genotyped individuals far from base populations, or many unknown parent groups (within breed per year of birth). METHODS: We derive likelihood methods to estimate Γ . For a single metafounder, summary statistics of pedigree and genomic relationships allow deriving a cubic equation with the real root being the maximum likelihood (ML) estimate of Γ . This equation is tested with Lacaune sheep data. For several metafounders, we split the first derivative of the complete likelihood in a term related to Γ , and a second term related to Mendelian sampling variances. Approximating the first derivative by its first term results in a pseudo-EM algorithm that iteratively updates the estimate of Γ by the corresponding block of the H-matrix. The method extends to complex situations with groups defined by year of birth, modelling the increase of Γ using estimates of the rate of increase of inbreeding ( Δ F ), resulting in an expanded Γ and in a pseudo-EM+ Δ F algorithm. We compare these methods with the generalized least squares (GLS) method using simulated data: complex crosses of two breeds in equal or unsymmetrical proportions; and in two breeds, with 10 groups per year of birth within breed. We simulate genotyping in all generations or in the last ones. RESULTS: For a single metafounder, the ML estimates of the Lacaune data corresponded to the maximum. For simulated data, when genotypes were spread across all generations, both GLS and pseudo-EM(+ Δ F ) methods were accurate. With genotypes only available in the most recent generations, the GLS method was biased, whereas the pseudo-EM(+ Δ F ) approach yielded more accurate and unbiased estimates. CONCLUSIONS: We derived ML, pseudo-EM and pseudo-EM+ Δ F methods to estimate Γ in many realistic settings. Estimates are accurate in real and simulated data and have a low computational cost.

Skull asymmetry in various sheep breeds: Directional asymmetry and fluctuating asymmetry.

Sheep (Ovis aries) play an important role in the economy of Turkey and the Balkan Peninsula due to their use in farming. As a domesticated species, sheep's morphometric and morphological diversity is likely determined by selective breeding practices rather than geographic distribution. This study aimed to analyse four different sheep breed skulls and reveal skull asymmetry using geometric morphometric methods. For this purpose, 2D images of 52 sheep skulls from different breeds were analysed from the dorsal view of the skull, using 28 landmarks. In the comparison of sheep skulls from the dorsal view, the first principal components for directional asymmetry (DA) and fluctuating asymmetry (FA) were 32.98% and 39.62% of the total variation, respectively. Sharri and Ivesi (Awassi) sheep breeds had the broadest distribution of skull shapes among the breeds, while Lara e Polisit was the most conservative breed. DA was used as a measure of biomechanical constraints, and FA was used as an indicator of environmental stress. Consistent with DA, both differences in centroid size and shape between breeds were statistically significant. No differences between males and females related to asymmetry were revealed. Ivesi sheep revealed the highest fluctuating asymmetry. Geometric morphometric methods proved to be a useful tool for distinguishing differences in the shape of the skull of different sheep breeds and also can be useful for taxonomic purposes.

Morphological variation of tail bone among two chicken breeds and their F1 progeny.

Fancy breeds of Japanese indigenous chicken display extensive morphological diversity, particularly in tail feathers. Although marked differences in tail and bone traits have been reported between Tosa-jidori (wild type) and Minohikichabo (rich type) breeds, little is known about the pattern of genetic inheritance in cross experiments. Therefore, this study aimed to investigate the strain and sex effects, and inheritance patterns, in the morphometric variation of pygostyle bones among Tosa-jidori, Minohikichabo, and their F1 hybrids. Five morphological traits, angle of the apex of the pygostyle, pygostyle length, margo cranialis length, tail feather number, and body weight, were evaluated at the adult stage. A significant strain difference was detected in all traits, whereas significant sex differences were observed in only three traits, but not in the angle of the apex of the pygostyle and tail feather number. In F1 hybrids, the angle of the apex of the pygostyle was significantly different to that of Tosa-jidori but not that of Minohikichabo, whereas the pygostyle length and tail number of F1 hybrids were significantly different from those of Minohikichabo but not those of Tosa-jidori. A significant heterosis effect was found in the margo cranialis length and body weight. All five traits showed nonadditive inheritance patterns but varied in each trait between partial dominance (angle of the apex of pygostyle), full dominance (pygostyle length and tail feather number), and over-dominance (margo cranialis length and body weight). Interestingly, different patterns of genetic inheritance in the F1 hybrid were observed at different locations, even within the same pygostyle bone. Using the Japanese indigenous chicken model, these results provide a substantial step toward understanding the genetic architecture of morphology in chickens.

Behavioural synchronisation between different groups of dogs and wolves and their owners/handlers: Exploring the effect of breed and human interaction.

Dogs have previously been shown to synchronise their behaviour with their owner and the aim of this study was to test the effect of immediate interactions, breed, and the effects of domestication. The behavioural synchronisation test was conducted in outdoor enclosures and consisted of 30 s where the owner/handler was walking and 30 s of standing still. Three studies were conducted to explore the effect of immediate interaction (study A), the effect of breed group (study B), and the effect of domestication (study C). In study A, a group of twenty companion dogs of various breeds were tested after three different human interaction treatments: Ignore, Pet, and Play. The results showed that dogs adjusted their movement pattern to align with their owner's actions regardless of treatment. Furthermore, exploration, eye contact, and movement were all influenced by the owners moving pattern, and exploration also decreased after the Play treatment. In study B, the synchronisation test was performed after the Ignore treatment on three groups: 24 dogs of ancient dog breeds, 17 solitary hunting dogs, and 20 companion dogs (data from study A). Irrespective of the group, all dogs synchronised their moving behaviour with their owner. In addition, human walking positively influenced eye contact behaviour while simultaneously decreasing exploration behaviour. In study C, a group of six socialised pack-living wolves and six similarly socialised pack-living dogs were tested after the Ignore treatment. Interestingly, these animals did not alter their moving behaviour in response to their handler. In conclusion, dogs living together with humans synchronise with their owner's moving behaviour, while wolves and dogs living in packs do not. Hence, the degree of interspecies behavioural synchronisation may be influenced by the extent to which the dogs are immersed in everyday life with humans.

Phenotypic characterization of Guraghe and Jimma cattle breeds in Ethiopia: Implications for breed differentiation and in-situ conservation.

The aim of this study was to phenotypically characterize Guraghe and Jimma cattle breeds at the farm level and quantify the relationships between them. Eight morphometric measurements and sixteen morphological traits were recorded for a total of 313 (221 females and 92 males) randomly selected adult cattle from three purposively selected districts. Univariate and multivariate analysis procedures of Statistical Analysis Software (SAS 9.0) were used to analyze the data. Clear morphological and morphometric variations were not observed between the two cattle breeds. All of the studied cattle populations possessed straight-edged ears, a sloppy rump profile, and straight face and back profiles. Moreover, the majority of the studied cattle possess widely spaced curved horns, red-colored uniform body color patterns, and erected small humps located at the cervicothoracic position. In addition to their phenotypic similarities, multivariate analysis also failed to reveal significant differences between the two breeds. These results suggest the inseparable nature of the two cattle breeds. However, such similarities in phenotypic traits between the two cattle breeds do not necessarily indicate genetic similarities. Therefore, further genetic characterization is recommended to quantify the degree of genetic relationship between the breeds. In the meantime, it is recommended to design breed-specific in situ conservation as well as genetic improvement programs that consider cattle breeds as one. Furthermore, an inclusive and uniform breed name that can represent the two cattle populations is obtained from the country's steering committee for indigenous animal genetic resources.

Copy number variant scan in more than four thousand Holstein cows bred in Lombardy, Italy.

Copy Number Variants (CNV) are modifications affecting the genome sequence of DNA, for instance, they can be duplications or deletions of a considerable number of base pairs (i.e., greater than 1000 bp and up to millions of bp). Their impact on the variation of the phenotypic traits has been widely demonstrated. In addition, CNVs are a class of markers useful to identify the genetic biodiversity among populations related to adaptation to the environment. The aim of this study was to detect CNVs in more than four thousand Holstein cows, using information derived by a genotyping done with the GGP (GeneSeek Genomic Profiler) bovine 100K SNP chip. To detect CNV the SVS 8.9 software was used, then CNV regions (CNVRs) were detected. A total of 123,814 CNVs (4,150 non redundant) were called and aggregated into 1,397 CNVRs. The PCA results obtained using the CNVs information, showed that there is some variability among animals. For many genes annotated within the CNVRs, the role in immune response is well known, as well as their association with important and economic traits object of selection in Holstein, such as milk production and quality, udder conformation and body morphology. Comparison with reference revealed unique CNVRs of the Holstein breed, and others in common with Jersey and Brown. The information regarding CNVs represents a valuable resource to understand how this class of markers may improve the accuracy in prediction of genomic value, nowadays solely based on SNPs markers.

Genetic parameters for reproductive traits and somatic cell score in U.S. dairy goats.

There are no studies regarding the estimation of genetic parameters and genetic trends for reproductive traits and somatic cells in goats. Their knowledge allows optimization of selection schemes. The objective of this study was to estimate genetic parameters and genetic and phenotypic trends for age at first kidding (AFK), kidding interval (KIN) and somatic cell score (SCS). Analyses were conducted within and across seven US goat breeds, namely, Nubian (NU), Alpine (AL), LaMancha (LM), Toggenburg (TO), Saanen (SA), Nigerian Dwarf (ND) and Oberhasli (OB), and a set of all of these breeds (AB). The restricted maximum likelihood methodology and trivariate animal models were used. Genetic and phenotypic trends were estimated using regression models. The average and standard deviation of AFK, KIN and SCS for AB were 573.6 ± 178.5 days, 418.8 ± 125.5 days and 4.67 ± 2.23 Log2, respectively. The heritabilities (h2) and standard errors of AFK, KIN and SCS for AB were 0.28 ± 0.02, 0.04 ± 0.02 and 0.22 ± 0.01, respectively. The h2 ranged from 0.15 (SA) to 0.37 (NU) for AFK, from 0.04 (AB) to 0.10 (AL) for KIN, and from 0.11 (TO) to 0.26 (LM and ND) for SCS. Genetic correlations between AFK and KIN and between AFK and SCS for AB were positive and weak (0.07 and 0.12, respectively) but significant (P < 0.01). Genetic correlations between SCS and KIN were significant (P < 0.01) for all the breeds and ranged from -0.15 (NU) to 0.44 (AL). Genetic correlations between AFK and SCS in the NU and AL breeds were similar (approximately 0.21). A positive genetic trend was found for KIN in the SA breed, which caused an increase in the number of days between consecutive kiddings. The genetic trend of SCS for the NU, AL and ND breeds was negative and decreased annually, which is beneficial for producers. These first results show the intensity and direction of some favorable/unfavorable relationships between AFK or KIN and SCS Log2 in some U.S. goat genetic groups.

Conservation and upgradation of indigenous Cholistani cattle breed of Pakistan: a pathway to sustainable livestock.

Pakistan is endowed with many established indigenous zebu Bos indicus type (humped) cattle breeds including Sahiwal, Red Sindhi, Bhagnari and Cholistani. Amongst these indigenous cattle breeds, Sahiwal and Red Sindhi have extensively been navigated and hence these two are acclaimed as internationally recognized breeds. However, research work on Cholistani cattle breed actually initiated in 2010 and has attained a steady pace. This breed was a new entrant in Livestock Census of Pakistan since 2006. Cholistani is a hardy, tick-resistant, adaptable cattle breed being reared under pastoral nomadism of the Cholistan desert, Pakistan. The present narrative review is the first of its kind intended to sum-up all the research work conducted about this indigenous cattle breed, and to put forth research gaps for this formerly neglected cattle breed. The review discusses the research work conducted on Cholistani cattle breed under five major research subjects/domains i.e. production attributes, theriogenology-related attributes, hematochemical attributes, disease, epidemiologic and therapeutic attributes, and genetic attributes. Future horizon for research avenues has also been given. It is the dire need of time that specific breed-oriented conservation and propagation programs may be initiated in the country so that sustained livestock and enhance socioeconomic profiling of rural communities may be attained.

Genome-wide scans for selection signatures in indigenous chickens reveal candidate genes associated with local adaptation.

Population growth and climate change pose challenges to the sustainability of poultry farming. The emphasis on high-yield traits in commercialized breeds has led to a decline in their adaptability. Chicken varieties adapted to the local environment, possessing traits that facilitate adaptation to climate change, such as disease resistance and tolerance to extreme weather conditions, can improve hybridization outcomes. In this study, we conducted an analysis of the population structure and genetic diversity of 110 chickens representing indigenous breeds from southern China and two different commercial breeds. Further, we performed comparative population genomics, utilizing nucleotide diversity and fixation statistics, to characterize genomic features of natural selection and to identify unique genetic traits and potential selection markers developed by indigenous breeds after adapting to the local environment. Results based on genetic diversity and population structure analyses showed that indigenous varieties exhibited high levels of genetic diversity. Commercial breeds that have been indigenously bred demonstrated higher levels of genetic diversity than those that have not, and breeds with different selection practices displayed significant differences in genetic structure. Additionally, we further searched for potential genomic regions in native chicken ecotypes, uncovering several candidate genes related to ecological adaptations affecting local breeds, such as IKBKB, S1PR1, TSHR, IL1RAPL1 and AMY2A, which are involved in disease resistance, heat tolerance, immune regulation and behavioral traits. This work provides important insights into the genomic characterization of ecotypes of native chicken in southern China. The identification of candidate genes associated with traits such as disease resistance, heat tolerance, immunomodulation, and behavioral traits is a significant outcome. These candidate genes may contribute to the understanding of the molecular basis of the adaptation of the southern native chicken to the local environment. It is recommended that these genes be integrated into chicken breeding programs to enhance sustainable agriculture and promote effective conservation and utilization strategies.

Evaluation of non-additive genetic effects on carcass and meat quality traits in Korean Hanwoo cattle using genomic models.

The traditional genetic evaluation methods generally consider additive genetic effects only and often ignore non-additive (dominance and epistasis) effects that may have contributed to genetic variation of complex traits of livestock species. The available dense single nucleotide polymorphisms (SNPs) panels offer to investigate the potential benefits of including non-additive genetic effects in the genomic evaluation models. Data from 16 971 genotyped (Illumina Bovine 50 K SNP chip) Korean Hanwoo cattle were used to estimate genetic variance components and prediction accuracy of genomic breeding values (GEBVs) for four carcass and meat quality traits: carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT) and marbling score (MS). Five different genetic models were evaluated through including additive, dominance and epistatic interactions (additive by additive, A × A; additive by dominance, A × D and dominance by dominance, D × D) successively in the models. The estimates of additive genetic variances and narrow sense heritabilities (ha2) were found similar across the evaluated models and traits except when additive interaction (A × A) was included. The dominance variance estimates relative to phenotypic variance ranged from 1.7-3.4% for CWT and MS traits, whereas, they were close to zero for EMA and BFT traits. The magnitude of A × A epistatic heritability (haa2) ranged between 14.8 and 27.7% in all traits. However, heritability estimates for A × D and D × D epistatic interactions (had2 and hdd2) were quite low compared to haa2 and were contributed only 0.0-9.7% of the total phenotypic variation. In general, broad sense heritability (hG2) estimates were almost twice (ranging between 0.54 and 0.68) the ha2 for all of the investigated traits. The inclusion of dominance effects did not improve the prediction accuracy of GEBV but improved 2.0-3.0% when epistatic effects were included in the model. More importantly, rank correlation revealed that partitioning of variance components considering dominance and epistatic effects in the model would enable to re-rank of top animals with better prediction of GEBV. The present result suggests that dominance and epistatic effects could be included in the genomic evaluation model for better estimates of variance components and more accurate prediction of GEBV for carcass and meat quality traits in Korean Hanwoo cattle.

Genetic characterization of the Latvian local goat breed and genetic traits associated with somatic cell count.

The Latvian local goat (LVK) breed represents the only native domestic goat breed in Latvia, but its limited population places it within the endangered category. However, the LVK breed has not yet undergone a comprehensive genetic characterization. Therefore, we completed whole genome sequencing to reveal the genetic foundation of the LVK breed while identifying genetic traits linked to the somatic cell count (SCC) levels. The study included 40 genomes of LVK goats sequenced to acquire at least 35x or 10x coverage. A Principal component analysis, a genetic distance tree, and an admixture analysis showed LVK's similarity to some European breeds, such as Finnish Landrace, Alpine, and Saanen, which aligns with the breed's history. An analysis of genome-wide heterozygosity, nucleotide diversity, and LD analysis indicated that the LVK population exhibits substantial levels of genetic diversity. LVK genome was dominated by short runs of homozygosity (ROHs, ≤ 500 kb) with a median length of 25 kb. With FROH 2.49%, average inbreeding levels were low; however, FROH ranged broadly from 0.13 to 12.2%. With the exception of one pure-blood breeding buck exhibiting FROH of 9.3% and FSNP of 8.5%, animals with at least 66% LVK ancestry showed moderate or no inbreeding. Overall, this study demonstrated that the LVK goats can be differentiated from imported breeds, although the population has a complex genetic structure. We were able to identify potential genetic traits associated with SCC levels, although the kinship of the animals and the heterogenic substructure of the population might have largely influenced the association analysis. We identified 26 genetic variants associated with SCC levels, which included the potentially relevant SNP rs662053371 in the OSBPL8 gene, indicating a potential signal linked to lipid metabolism in goats. To conclude, these findings present valuable insight into the genetic structure of the LVK breed for the conservation of local genetic resources.

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.

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats.

BACKGROUND: Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate. RESULTS: To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds. CONCLUSIONS: The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

Accuracy of genomic prediction using multiple Atlantic salmon populations.

BACKGROUND: The accuracy of genomic prediction is partly determined by the size of the reference population. In Atlantic salmon breeding programs, four parallel populations often exist, thus offering the opportunity to increase the size of the reference set by combining these populations. By allowing a reduction in the number of records per population, multi-population prediction can potentially reduce cost and welfare issues related to the recording of traits, particularly for diseases. In this study, we evaluated the accuracy of multi- and across-population prediction of breeding values for resistance to amoebic gill disease (AGD) using all single nucleotide polymorphisms (SNPs) on a 55K chip or a selected subset of SNPs based on the signs of allele substitution effect estimates across populations, using both linear and nonlinear genomic prediction (GP) models in Atlantic salmon populations. In addition, we investigated genetic distance, genetic correlation estimated based on genomic relationships, and persistency of linkage disequilibrium (LD) phase across these populations. RESULTS: The genetic distance between populations ranged from 0.03 to 0.07, while the genetic correlation ranged from 0.19 to 0.99. Nonetheless, compared to within-population prediction, there was limited or no impact of combining populations for multi-population prediction across the various models used or when using the selected subset of SNPs. The estimates of across-population prediction accuracy were low and to some extent proportional to the genetic correlation estimates. The persistency of LD phase between adjacent markers across populations using all SNP data ranged from 0.51 to 0.65, indicating that LD is poorly conserved across the studied populations. CONCLUSIONS: Our results show that a high genetic correlation and a high genetic relationship between populations do not guarantee a higher prediction accuracy from multi-population genomic prediction in Atlantic salmon.