Category Archives: Journal of Animal Breeding and Genetics

Investigation of phenotypic, genetic and genomic background of Milk spectra in Sarda dairy sheep

Posted on by

Abstract

Aim of this study was to analyse the genetic background of milk Fourier transform infrared (FTIR) spectra in dairy sheep. Individual milk FTIR spectra, with 1060 wavenumbers each, were available for 793 adult Sarda breed ewes genotyped at 45,813 SNP. The absorbance values of each wavenumber was analysed using a linear mixed model that included dim class, parity and lambing month as fixed effects and flock-test date and animal as random effects. The model was applied to estimate variance components and heritability and to perform a genome-wide association study for each wavenumber. Average h2 of wavenumbers absorbance was 0.13 ± 0.08, with the largest values observed in the regions associated with the characteristic bonds of carbonylic and methylenic groups of milk fat (h2 = 0.57 at 1724–1728 cm−1; and h2 = 0.34 at 2811–2834 cm−1, respectively). The absorbance values of wavenumbers were moderately correlated with the estimated heritabilities. After the Bonferroni correction, a total of nine markers were found to be significantly associated with 32 different wavenumbers. Of particular interest was the SNP s63269.1, mapped on chromosome 2, that was found to be associated with 27 wavenumbers. Genes previously found to be related to traits of interest (e.g. disease resistance, milk yield and quality, cheese firmness) are located close to the significant SNP. As expected, the heritability estimated for the absorbance of each wavenumbers seems to be associated with the related milk components.

Genome‐wide association studies for diarrhoea outcomes identified genomic regions affecting resistance to a severe enteropathy in suckling rabbits

Posted on by

Abstract

Selection and breeding strategies to improve resistance to enteropathies are essential to reaching the sustainability of the rabbit production systems. However, disease heterogeneity (having only as major visible symptom diarrhoea) and low disease heritability are two barriers for the implementation of these strategies. Diarrhoea condition can affect rabbits at different life stages, starting from the suckling period, with large negative economic impacts. In this study, from a commercial population of suckling rabbits (derived from 133 litters) that experienced an outbreak of enteropathy, we first selected a few animals that died with severe symptoms of diarrhoea and characterized their microbiota, using 16S rRNA gene sequencing data. Clostridium genus was consistently present in all affected specimens. In addition, with the aim to identify genetic markers in the rabbit genome that could be used as selection tools, we performed genome-wide association studies for symptoms of diarrhoea in the same commercial rabbit population. These studies were also complemented with F ST analyses between the same groups of rabbits. A total of 332 suckling rabbits (151 with severe symptoms of diarrhoea, 42 with mild symptoms and 129 without any symptoms till the weaning period), derived from 45 different litters (a subset of the 133 litters) were genotyped with the Affymetrix Axiom OrcunSNP Array. In both genomic approaches, rabbits within litters were paired to constitute two groups (susceptible and resistant, including the mildly affected in one or the other group) and run case and control genome-wide association analyses. Genomic heritability estimated in the designed experimental structure integrated in a commercial breeding scheme was 0.19–0.21 (s.e. 0.09–0.10). A total of eight genomic regions on rabbit chromosome 2 (OCU2), OCU3, OCU7, OCU12, OCU13, OCU16 and in an unassembled scaffold had significant single nucleotide polymorphisms (SNPs) and/or markers that trespassed the FST percentile distribution. Among these regions, three main peaks of SNPs were identified on OCU12, OCU13 and OCU16. The QTL region on OCU13 encompasses several genes that encode members of a family of immunoglobulin Fc receptors (FCER1G, FCRLA, FCRLB and FCGR2A) involved in the immune innate system, which might be important candidate genes for this pathogenic condition. The results obtained in this study demonstrated that resistance to an enteropathy occurring in suckling rabbits is in part genetically determined and can be dissected at the genomic level, providing DNA markers that could be used in breeding programmes to increase resistance to enteropathies in meat rabbits.

Effects of reference population size and structure on genomic prediction of maternal traits in two pig lines using whole‐genome sequence‐, high‐density‐ and combined annotation‐dependent depletion genotypes

Posted on by

Abstract

The aim of this study was to investigate the reference population size required to obtain substantial prediction accuracy within- and across-lines and the effect of using a multi-line reference population for genomic predictions of maternal traits in pigs. The data consisted of two nucleus pig populations, one pure-bred Landrace (L) and one Synthetic (S) Yorkshire/Large White line. All animals were genotyped with up to 30 K animals in each line, and all had records on maternal traits. Prediction accuracy was tested with three different marker data sets: High-density SNP (HD), whole genome sequence (WGS), and markers derived from WGS based on pig combined annotation dependent depletion-score (pCADD). Also, two different genomic prediction methods (GBLUP and Bayes GC) were compared for four maternal traits; total number piglets born (TNB), total number of stillborn piglets (STB), Shoulder Lesion Score and Body Condition Score. The main results from this study showed that a reference population of 3 K–6 K animals for within-line prediction generally was sufficient to achieve high prediction accuracy. However, when the number of animals in the reference population was increased to 30 K, the prediction accuracy significantly increased for the traits TNB and STB. For multi-line prediction accuracy, the accuracy was most dependent on the number of within-line animals in the reference data. The S-line provided a generally higher prediction accuracy compared to the L-line. Using pCADD scores to reduce the number of markers from WGS data in combination with the GBLUP method generally reduced prediction accuracies relative to GBLUP using HD genotypes. The BayesGC method benefited from a large reference population and was less dependent on the different genotype marker datasets to achieve a high prediction accuracy.

Prevalence of sex‐chromosome aneuploidy estimated using SNP genotype intensity information in a large population of juvenile dairy and beef cattle

Posted on by

Abstract

Aneuploidy is a genetic condition characterized by the loss or gain of one or more chromosomes. Aneuploidy affecting the sex chromosomes can lead to infertility in otherwise externally phenotypically normal cattle. Early identification of cattle with sex chromosomal aneuploidy is important to minimize the costs associated with rearing infertile cattle and futile breeding attempts. As most livestock breeding programs routinely genotype their breeding populations using single nucleotide polymorphism (SNP) arrays, this study aimed to assess the feasibility of integrating an aneuploidy screening tool into the existing pipelines that handle dense SNP genotype data. A further objective was to estimate the prevalence of sex chromosome aneuploidy in a population of 146,431 juvenile cattle using available genotype intensity data. Three genotype intensity statistics were used: the LogR Ratio (LRR), R-value (the sum of X and Y SNP probe intensities), and B-allele frequency (BAF) measurements. Within the female-verified population of 124,958 individuals, the estimated prevalence rate was 0.0048% for XO, 0.0350% for XXX, and 0.0004% for XXY. The prevalence of XXY in the male-verified population was 0.0870% (i.e., 18 out of 20,670 males). Cytogenetic testing was used to verify 2 of the XXX females who were still alive. The proposed approach can be readily integrated into existing genomic pipelines, serving as an efficient, large-scale screening tool for aneuploidy. Its implementation could enable the early identification of infertile animals with sex-chromosome aneuploidy.

Estimation of heritability with genomic information by method R

Posted on by

Abstract

Estimating heritabilities with large genomic models by established methods such as restricted maximum likelihood (REML) or Bayesian via Gibbs sampling is computationally expensive. Alternatively, heritability can be estimated indirectly by method R and by maximum predictivity, referred to as MaxPred here, at a much lower computing cost. By method R, the heritability used for predictions with whole and partial data is considered the best estimate when the predictions based on partial data are unbiased relative to those with the complete data. By MaxPred, the heritability estimate is the one that maximizes predictivity. This study compared heritability estimation with genomic information using average information REML (AI–REML), method R and MaxPred. A simulated population was generated with ten generations of 5000 animals each and an effective population size of 80. Each animal had one record for a trait with a heritability of 0.3, a phenotypic variance of 10.0 and was genotyped at 50 k SNP. In method R, the heritability estimate is found when the expectation of a regression coefficient is equal to one. The regression is the EBV of selection candidates calculated with the whole dataset regressed on the EBV of candidates calculated from a partial dataset. In this study, we used the GBLUP framework and therefore, GEBV was calculated. The partial dataset was created by removing the last generation of phenotypes. Predictivity was defined as the correlation between the adjusted phenotypes of the selection candidates and their GEBV calculated from the partial data. We estimated the heritability for populations that included between three and 10 generations. In every scenario, predictivity increased as more data was used and was the highest at the simulated heritability. However, the predictivity for all data subsets and all heritabilities compared did not differ more than 0.01, suggesting MaxPred is not the best indication for heritability estimation. For the whole dataset, the heritability was estimated as 0.30 ± 0.01, 0.26 ± 0.01 and 0.30 ± 0.04 for AI–REML without genomics, AI–REML with genomics and method R with genomics, respectively. Heritability estimation with genomics by method R reduced timing by 83%, implying a reduction in computing time from 9.5 to 1.6 h, on average, compared to AI–REML with genomics. Method R has the potential to estimate heritabilities with large genomic information at a low cost when many generations of animals are present; however, the standard error can be high when only a few iterations are used.

Genomic analyses of nitrogen utilization efficiency, its indicator trait blood urea nitrogen and the relationship to classical growth performance and feed efficiency traits in a Landrace × Piétrain crossbred population

Posted on by

Abstract

Improving the nutrient efficiency in pork production is required to reduce the resource competition between human food and animal feed regarding diet components edible for humans and to minimize emissions relevant to climate or the environment. Thereby, protein utilization efficiency and its equivalent nitrogen utilization efficiency (NUE) play a major role. Breeding for more nitrogen (N) efficient pigs bears a promising strategy to improve such traits, however, directly phenotyping NUE based on N balance data is neither cost-efficient nor straightforward and not applicable for routine evaluations. Blood urea nitrogen (BUN) levels in the pig are suitable to predict the NUE and, therefore, might be an indicator trait for NUE because BUN is a relatively easy-to-measure trait. This study investigated the suitability of NUE as a selection trait in future breeding programs. The relationships to classical growth performance and feed efficiency traits were analysed as well as the relationship to BUN to infer the role of BUN as an indicator trait to improve NUE via breeding. The analyzes were based on a Landrace F1 cross population consisting of 502 individuals who descended from 20 Piétrain sires. All animals were genotyped for 48,525 SNPs. They were phenotyped in two different fattening phases, i.e., FP1 and FP2, during the experiment. Uni- and bivariate analyses were run to estimate variance components and to determine the genetic correlation between different traits or between the same trait measured at different time points. Moderate heritabilities were estimated for all traits, whereby the heritability for NUE was h 2 = 0.293 in FP1 and h 2 = 0.163 in FP2 and BUN had the by far highest heritability (h 2 = 0.415 in FP1 and h 2 = 0.460 in FP2). The significant genetic correlation between NUE and BUN showed the potential of BUN to be considered an indicator trait for NUE. This was particularly pronounced when NUE was measured in FP1 (genetic correlations rg=−0.631$$ {r}_g=-0.631 $$ and rg=−0.688$$ {r}_g=-0.688 $$ between NUE and BUN measured in FP1 and FP2, respectively). The genetic correlations of NUE and BUN with important production traits suggest selecting pigs with high growth rates and low BUN levels to breed more efficient pigs in future breeding programs.

Genetic diversity, population structure and origin of the native goats in Central Laos

Posted on by

Abstract

Maintaining genetic diversity and variation in livestock populations is critical for natural and artificial selection promoting genetic improvement while avoiding problems due to inbreeding. In Laos, there are concerns that there has been a decline in genetic diversity and a rise in inbreeding among native goats in their village-based smallholder system. In this study, we investigated the genetic diversity of Lao native goats in Phin, Songkhone and Sepon districts in Central Laos for the first time using Illumina's Goat SNP50 BeadChip. We also explored the genetic relationships between Lao goats with 163 global goat populations from 36 countries. Our results revealled a close genetic relationship between Lao native goats and Chinese, Mongolian and Pakistani goats, sharing ancestries with Guangfen, Jining Grey and Luoping Yellow breeds (China) and Teddi goats (Pakistan). The observed (Ho) and expected (He) heterozygosity were 0.292 and 0.303 (Laos), 0.288 and 0.288 (Sepon), 0.299 and 0.308 (Phin) and 0.289 and 0.305 (Songkhone), respectively. There was low to moderate genetic differentiation (F ST: 0.011–0.043) and negligible inbreeding coefficients (F IS: −0.001 to 0.052) between goat districts. The runs of homozygosity (ROH) had an average length of 5.92–6.85 Mb, with short ROH segments (1–5 Mb length) being the most prevalent (66.34%). Longer ROH segments (20–40 and >40 Mb length categories) were less common, comprising only 4.81% and 1.01%, respectively. Lao goats exhibit moderate genetic diversity, low-inbreeding levels and adequate effective population size. Some genetic distinctions between Lao goats may be explained by geographic and cultural features.

Genetic parameters for different types of medullated fibre in Alpacas

Posted on by

Abstract

The quality of alpaca textile fibre has great potential, especially if objectionable fibres (coarse and medullated fibres) that cause itching are reduced, considering that objectionable fibres can be identified by diameter and medullation types. The objective of this study was to estimate genetic parameters for medullar types and their respective diameters to evaluate the possibility of incorporating them as selection criteria in alpaca breeding programmes. The research used 3149 alpaca fibre samples collected from 2020 to 2022, from a population of 1626 Huacaya type alpacas. The heritability and correlations of the percentages of non-medullated (NM), fragmented medulle (FM), uncontinuous medullated (UM), continuous medullated (CM), and strongly medullated (SM) fibres were analysed, also the fibre diameter (FD) for each of the medullation types. The heritability estimated for medullation types were 0.25 ± 0.01, 0.18 ± 0.01, 0.10 ± 0.01, 0.20 ± 0.01 and 0.11 ± 0.01 for NM, FM, UM, CM and SM, respectively. The genetic correlations for medullation categories ranged from 0.15 ± 0.03 to 0.66 ± 0.02 (in absolute values). The heritabilility estimated for fibre diameter (FD) of each of the medullation types were 0.29 ± 0.03, 0.27 ± 0.02, 0.35 ± 0.02, 0.30 ± 0.02, 0.25 ± 0.02 and 0.10 ± 0.02 for FD, FD_NM, FD_FM, FD_UM, FD_CM and FD_SM, respectively. The genetic correlations for fibre diameter of the medullation types ranged from 0.04 ± 0.04 to 0.97 ± 0.01. FD, NM and FM are the main traits to be used as selection criteria under a genetic index, since they would reduce fibre diameter, and also increase NM and FM, and, in addition reducing indirectly CM, SM, and SM_FD. Therefore, the quality of alpaca fibre could be improved.

Characterization of runs of homozygosity islands in American mink using whole‐genome sequencing data

Posted on by

Abstract

The genome-wide analysis of runs of homozygosity (ROH) islands can be an effective strategy for identifying shared variants within a population and uncovering important genomic regions related to complex traits. The current study performed ROH analysis to characterize the genome-wide patterns of homozygosity, identify ROH islands and annotated genes within these candidate regions using whole-genome sequencing data from 100 American mink (Neogale vison). After sequence processing, variants were called using GATK and Samtools pipelines. Subsequent to quality control, 8,373,854 bi-allelic variants identified by both pipelines remained for further analysis. A total of 34,652 ROH segments were identified in all individuals, among which shorter segments (0.3–1 Mb) were abundant throughout the genome, approximately accounting for 84.39% of all ROH. Within these segments, we identified 63 ROH islands housing 156 annotated genes. The genes located in ROH islands were associated with fur quality (EDNRA, FGF2, FOXA2 and SLC24A4), body size/weight (MYLK4, PRIM2, FABP2, EYS and PHF3), immune capacity (IL2, IL21, PTP4A1, SEMA4C, JAK2, CCNA2 and TNIP3) and reproduction (ADAD1, KHDRBS2, INSL6, PGRMC2 and HSPA4L). Furthermore, Gene Ontology and KEGG pathway enrichment analyses revealed 56 and 9 significant terms (FDR-corrected p-value < 0.05), respectively, among which cGMP-PKG signalling pathway, regulation of actin cytoskeleton, and calcium signalling pathway were highlighted due to their functional roles in growth and fur characteristics. This is the first study to present ROH islands in American mink. The candidate genes from ROH islands and functional enrichment analysis suggest possible signatures of selection in response to the mink breeding targets, such as increased body length, reproductive performance and fur quality. These findings contribute to our understanding of genetic characteristics, and provide complementary information to assist with implementation of breeding strategies for genetic improvement in American mink.

Identification of genomic regions and pathways associated with traits related to rumen acidosis in feedlot Nellore cattle

Posted on by

Abstract

There may be an increased risk of metabolic disorders, such as rumen acidosis, in cattle fed high-concentrate diets, particularly those from Bos taurus indicus genotypes, which have shown to be more sensitive to ruminal acidification. Therefore, this study aimed to estimate (co)variance components and identify genomic regions and pathways associated with ruminal acidosis in feedlot Nellore cattle fed high-concentrate diets. It was utilized a dataset containing a total of 642 Nellore bulls that were genotyped from seven feedlot nutrition studies. The GGP Indicus 35k panel was used with the single step genome-wide association study methodology in which the effects of the markers were obtained from the genomic values estimated by the GBLUP model. A bivariate model to estimate genetic correlations between the economically important traits and indicator traits for acidosis was used. The traits evaluated in this study that were nutritionally related to rumen acidosis included average daily gain (ADG), final body weight, time spent eating (TSE), time spent ruminating, rumenitis score (RUM), rumen absorptive surface area (ASA), rumen keratinized layer thickness (KER) and hot carcass weight (HCW). The identified candidate genes were mainly involved in the negative or non-regulation of the apoptotic process, salivary secretion, and transmembrane transport. The genetic correlation between HCW and ASA was low positive (0.27 ± 0.23), and between ADG and ASA was high moderate (0.58 ± 0.59). A positive genetic correlation between RUM and all performance traits was observed, and TSE correlated negatively with HCW (−0.33 ± 0.21), ASA (−0.75 ± 0.48), and KER (−0.40 ± 0.27). The genetic association between economically important traits and indicator traits for acidosis suggested that Nellore cattle may be more sensitive to acidosis in feedlot systems.