Reaction-norm animal model analysis of average daily gain heat tolerance in purebred Duroc pigs.

The present study aimed to genetically improve growth performance under high-heat environments by specifically designing a reaction-norm animal model (RNAM) for purebred Duroc pigs in Japan. A total of 54,750 records of average daily gain (ADG) measured for pigs reared at four farms in different prefectures were analyzed. Estimated maximum daily temperatures at the respective farm locations were used to calculate the average cumulative thermal load (TL). The TL values served as an indicator of high-heat environments for pigs. The plausible cumulative period length and threshold temperature for calculating TL were determined to be 8 weeks until just before shipping and 25°C, respectively. Variance components were estimated via RNAM analysis using TL as a linear covariate. The estimated additive genetic variances under both responsive and non-responsive to TL were found to be significant. Moreover, the estimated heritability of ADG ranged from 0.38 to 0.73 for TL values of 0-8. These results suggest that the RNAM developed holds the potential for improving the genetic ability of growth under high-heat environments in pigs.

Plateau-linear reaction norm model analysis of number born alive in purebred Landrace pigs using meteorological data in Japan.

We performed a plateau-linear reaction norm model (RNM) analysis of number born alive (NBA) in purebred Landrace pigs, where breeding value changes according to maximum temperature at mating day, using public meteorological observation data in Japan. We analysed 52,668 NBA records obtained from 10,320 Landrace sows. Pedigree data contained 99,201 animals. Off-farm daily temperature data at the nearest weather station from each of the farms were downloaded from the Japan Meteorological Agency website. A plateau-linear RNM analysis based on daily maximum temperature on mating day (threshold temperature of 16.6°C) was performed. The percentage of the records with daily maximum temperatures at mating days of ≤16.6, ≥25.0, ≥30.0 and ≥35.0°C were 34.3%, 33.6%, 14.0% and 0.8%, respectively. The value of Akaike's information criterion for the plateau-linear RNM was lower than that for a simple repeatability model (RM). With the plateau-linear RNM, estimated value of heritability ranged from 0.14 to 0.15, while that from the RM analysis was 0.15. Additive genetic correlation between intercept and slope terms was estimated to be -0.52 from the plateau-linear RNM analysis. Estimated additive genetic correlations were >0.9 between NBA at different temperatures ranging from 16.6 to 37.6°C. For the 10,320 sows, average values of prediction reliability of the intercept and slope terms for breeding values in the plateau-linear RNM were 0.47 and 0.16, respectively. Increasing weight for slope term in linear selection index could bring positive genetic gain in the slope part, but prediction accuracy would decrease. Our results imply that genetically improving heat tolerance in sows reared in Japan focusing on NBA using RNM is possible, while RNM is more complex to implement and interpret. Therefore, further study should be encouraged to make genetic improvement for heat tolerance in sows more efficient.

Estimation of off-farm temperatures and temperature-humidity index from meteorological observation records.

The objective of this study was to devise an optimal method for estimating air temperatures outside pig farms to be able to evaluate the genetic performance of pigs. Using daily temperature data from Japan Meteorological Agency meteorological stations, we investigated the optimal number of observation weather stations (number of records), and methods of estimating outside temperature when temperature records are missing. We also considered the possibility of using relative humidity data. Our results showed that it is possible to use records from the three nearest weather stations to estimate off-farm ambient temperatures. We also concluded that estimates of outside temperatures when records are missing can be made by using data from at least one weather station that holds a full set of data. The correlation coefficients between the true THI (temperature-humidity index) and the estimated THI and the average daily temperature were almost the same, indicating that the daily average temperature can be used instead of estimated THI.

Application of linear and machine learning models to genomic prediction of fatty acid composition in Japanese Black cattle.

We collected 3180 records of oleic acid (C18:1) and monounsaturated fatty acid (MUFA) measured using gas chromatography (GC) and 6960 records of C18:1 and MUFA measured using near-infrared spectroscopy (NIRS) in intermuscular fat samples of Japanese Black cattle. We compared genomic prediction performance for four linear models (genomic best linear unbiased prediction [GBLUP], kinship-adjusted multiple loci [KAML], BayesC, and BayesLASSO) and five machine learning models (Gaussian kernel [GK], deep kernel [DK], random forest [RF], extreme gradient boost [XGB], and convolutional neural network [CNN]). For GC-based C18:1 and MUFA, KAML showed the highest accuracies, followed by BayesC, XGB, DK, GK, and BayesLASSO, with more than 6% gain of accuracy by KAML over GBLUP. Meanwhile, DK had the highest prediction accuracy for NIRS-based C18:1 and MUFA, but the difference in accuracies between DK and KAML was slight. For all traits, accuracies of RF and CNN were lower than those of GBLUP. The KAML extends GBLUP methods, of which marker effects are weighted, and involves only additive genetic effects; whereas machine learning methods capture non-additive genetic effects. Thus, KAML is the most suitable method for breeding of fatty acid composition in Japanese Black cattle.

Population structure of Vietnamese pigs using mitochondrial DNA.

The D-loop region on mitochondrial DNA (mtDNA) is frequently used for analyses of maternal lineages within domestic animal species. There are many native pig breeds in Vietnam, but their origins remain unclear. This study investigated maternal lineages using the D-loop region on mtDNA of 260 samples collected from native pigs in 20 provinces across Vietnam. The D-loop region of all samples was amplified and sequenced. We obtained 713 bp sequences of the D-loop region for each sample excluding the repeat region, and variants on this region were used to construct a phylogenetic tree. We detected 50 haplotypes from Vietnamese native pigs, with 27 novel haplotypes. Phylogenetic tree analysis showed two haplotype groups: one for the MTSEA group, frequently found in domestic pigs in the mountainous areas of Cambodia and Laos; and the D2 group, found in pigs originating from Chinese pigs. No European haplotype was found. Haplotypes in northeast Vietnam comprised only haplotypes of the D2 group, whereas in areas from the northwest mountains to the south, we found haplotypes belonging to both the D2 and MTSEA groups. This study suggested that both origins contributed to maternal lineages of current populations of Vietnamese native pigs.

Comparing pedigree and genomic inbreeding coefficients, and inbreeding depression of reproductive traits in Japanese Black cattle.

BACKGROUND: Pedigree-based inbreeding coefficients have been generally included in statistical models for genetic evaluation of Japanese Black cattle. The use of genomic data is expected to provide precise assessment of inbreeding level and depression. Recently, many measures have been used for genome-based inbreeding coefficients; however, with no consensus on which is the most appropriate. Therefore, we compared the pedigree- ([Formula: see text]) and multiple genome-based inbreeding coefficients, which were calculated from the genomic relationship matrix with observed allele frequencies ([Formula: see text]), correlation between uniting gametes ([Formula: see text]), the observed vs expected number of homozygous genotypes ([Formula: see text]), runs of homozygosity (ROH) segments ([Formula: see text]) and heterozygosity by descent segments ([Formula: see text]). We quantified inbreeding depression from estimating regression coefficients of inbreeding coefficients on three reproductive traits: age at first calving (AFC), calving difficulty (CD) and gestation length (GL) in Japanese Black cattle. RESULTS: The highest correlations with [Formula: see text] were for [Formula: see text] (0.86) and [Formula: see text] (0.85) whereas [Formula: see text] and [Formula: see text] provided weak correlations with [Formula: see text], with range 0.33-0.55. Except for [Formula: see text] and [Formula: see text], there were strong correlations among genome-based inbreeding coefficients ([Formula: see text] 0.94). The estimates of regression coefficients of inbreeding depression for [Formula: see text] was 2.1 for AFC, 0.63 for CD and -1.21 for GL, respectively, but [Formula: see text] had no significant effects on all traits. Genome-based inbreeding coefficients provided larger effects on all reproductive traits than [Formula: see text]. In particular, for CD, all estimated regression coefficients for genome-based inbreeding coefficients were significant, and for GL, that for [Formula: see text] had a significant.. Although there were no significant effects when using overall genome-level inbreeding coefficients for AFC and GL, [Formula: see text] provided significant effects at chromosomal level in four chromosomes for AFC, three chromosomes for CD, and two chromosomes for GL. In addition, similar results were obtained for [Formula: see text]. CONCLUSIONS: Genome-based inbreeding coefficients can capture more phenotypic variation than [Formula: see text]. In particular, [Formula: see text] and [Formula: see text] can be considered good estimators for quantifying inbreeding level and identifying inbreeding depression at the chromosome level. These findings might improve the quantification of inbreeding and breeding programs using genome-based inbreeding coefficients.

Evaluating the performance of genomic prediction accounting for effects of single nucleotide polymorphism markers in reproductive traits of Japanese Black cattle.

We examined the prediction accuracies of genomic best linear unbiased prediction (GBLUP), various weighted GBLUP according to the degrees of marker effects (WGBLUP) and machine learning (ML) methods, and compared them with traditional BLUP for age at first calving (AFC), calving difficulty (CD), and gestation length in Japanese Black cattle. For WGBLUP, firstly, BayesC and FarmCPU were used to estimate marker effects. Then, we constructed three weighted genomic relationship matrices from information of estimated marker effects in the first step: absolute value of the estimated marker-effect WGBLUP, estimated marker-variance WGBLUP, and genomic-feature WGBLUP. For ML, we applied Gaussian kernel, random forest, extreme gradient boost, and support vector regression. We collected a total of 2583 animals having both phenotypic records and genotypes with 30,105 markers and 16,406 pedigree records. For AFC, prediction accuracies of WGBLUP methods using FarmCPU exceeded BLUP by 25.7%-39.5%. For CD, estimated marker-variance WGBLUP using BayesC achieved the highest prediction accuracy. Among ML methods, extreme gradient boost, support vector regression, and Gaussian kernel increased prediction accuracies by 28.4%, 19.0%, and 36.4% for AFC, CD, and gestation length compared with BLUP, respectively. Thus, prediction performance could be improved using suitable WGBLUP and ML methods according to target reproductive traits for the population used.

Genome-wide detection of changes in allelic frequency in Landrace pigs selected for resistance to mycoplasma pneumonia of swine.

Closed-pig line breeding could change the genetic structure at a genome-wide scale because of the selection in a pig breeding population. We investigated the changes in population structure among generations at a genome-wide scale and the selected loci across the genome by comparing the observed and expected allele frequency changes in mycoplasma pneumonia of swine (MPS)-selected pigs. Eight hundred and seventy-four Landrace pigs, selected for MPS resistance without reducing average daily gain over five generations, had 37,299 single nucleotide polymorphisms (SNPs) and were used for genomic analyses. Regarding population structure, individuals in the first generation were the most widely distributed and then converged into a specific group, as they were selected over five generations. For allele frequency changes, 96 and 14 SNPs had higher allele frequency changes than the 99.9% and 99.99% thresholds of the expected changes, respectively. These SNPs were evenly spread across the genome, and a few of these selected regions overlapped with previously detected quantitative trait loci for MPS and immune-related traits. Our results indicated that the considerable changes in allele frequency were identified in many regions across the genome by closed-pig line breeding based on estimated breeding value.

Genome-wide association studies for production, respiratory disease, and immune-related traits in Landrace pigs.

Identification of a quantitative trait locus (QTL) related to a chronic respiratory disease such as Mycoplasmal pneumonia of swine (MPS) and immune-related traits is important for the genetic improvement of disease resistance in pigs. The objective of this study was to detect a novel QTL for a total of 22 production, respiratory disease, and immune-related traits in Landrace pigs. A total of 874 Landrace purebred pigs, which were selected based on MPS resistance, were genotyped using the Illumina PorcineSNP60 BeadChip. We performed single nucleotide polymorphism (SNP)-based and haplotype-based genome-wide association studies (GWAS) to detect a novel QTL and to evaluate the possibility of a pleiotropic QTL for these traits. SNP-based GWAS detected a total of six significant regions in backfat thickness, ratio of granular leucocytes to lymphatic cells, plasma concentration of cortisol at different ages, and complement alternative pathway activity in serum. The significant region detected by haplotype-based GWAS was overlapped across the region detected by SNP-based GWAS. Most of these detected QTL regions were novel regions with some candidate genes located in them. With regard to a pleiotropic QTL among traits, only three of these detected QTL regions overlapped among traits, and many detected regions independently affected the traits.

Differential Metabolomics Profiles Identified by CE-TOFMS between High and Low Intramuscular Fat Amount in Fattening Pigs.

The amount of intramuscular fat (IMF) present in the loin eye area is one of the most important characteristics of high-quality pork. IMF measurements are currently impractical without a labor-intensive process. Metabolomic profiling could be used as an IMF indicator to avoid this process; however, no studies have investigated their use during the fattening period of pigs. This study examined the metabolite profiles in the plasma of two groups of pigs derived from the same Duroc genetic line and fed the same diet. Five plasma samples were collected from each individual the day before slaughter. Capillary electrophoresis-time of flight mass spectrometry (CE-TOFMS) was used to analyze the purified plasma from each sample. Principle component analysis (PCA) and partial least squares (PLS) were used to find the semi-quantitative values of the compounds. The results indicate that branched-chain amino acids are significantly associated with high IMF content, while amino acids are associated with low IMF content. These differences were validated using the quantification analyses by high-performance liquid chromatograph, which supported our results. These results suggest that the concentration of branched-chain amino acids in plasma could be an indicative biomarker for the IMF content in the loin eye area.

Inferring phenotypic causal structure among farrowing and weaning traits in pigs.

Direct selection for litter size or weight at weaning in pigs is often hindered by external interventions such as cross-fostering. The objective of this study was to infer the causal structure among phenotypes of reproductive traits in pigs to enable subsequent direct selection for these traits. Examined traits included: number born alive (NBA), litter size on day 21 (LS21), and litter weight on day 21 (LW21). The study included 6,240 litters from 1,673 Landrace dams and 5,393 litters from 1,484 Large White dams. The inductive causation (IC) algorithm was used to infer the causal structure, which was then fitted to a structural equation model (SEM) to estimate causal coefficients and genetic parameters. Based on the IC algorithm and temporal and biological information, the causal structure among traits was identified as: NBA â†’ LS21 â†’ LW21 and NBA â†’ LW21. Owing to the causal effect of NBA on LS21 and LW21, the genetic, permanent environmental, and residual variances of LS21 and LW21were much lower in the SEM than in the multiple-trait model for both breeds. Given the strong effect of NBA on LS21 and LW21, the SEM and causal information might assist with selective breeding for LS21 and LW21 when cross-fostering occurs.

Correlated response of peripheral blood cytokines with selection for reduced mycoplasma pneumonia of swine lesions in Landrace pigs.

Mycoplasma pneumonia of swine (MPS) is responsible for significant economic losses in the swine industry. We selected Landrace pigs for reduced MPS pulmonary lesions over five generations, and measured concentrations of the following cytokines: interleukin (IL)-10, IL-13, IL-17, tumor necrosis factor (TNF)-α and interferon (IFN)-γ to estimate their correlation with MPS lesions. Sheep red blood cells (SRBC) were injected twice intramuscularly at 70 and 95 kg body weight. Blood serum samples were collected after 1 week of secondary SRBC inoculation and cytokine concentrations were analyzed by ELISA. Genetic parameters and breeding values were estimated. The heritability estimates of IL-10, IL-13, IL-17, TNF-α and IFN-γ were 0.20 ± 0.06, 0.12 ± 0.06, 0.27 ± 0.07, 0.20 ± 0.10 and 0.05 ± 0.03, respectively. Genetic correlations of IL-17 and TNF-α with pulmonary MPS lesions were high (-0.86 ± 0.13 and 0.69 ± 0.29, respectively) and those of IFN-γ and IL-13 with MPS lesions were moderately negative (-0.45). Through selection, the breeding values of IL-17 and IFN-γ increased substantially and those of TNF-α decreased. These results suggest that innate and cellular immunity are more important for the suppression of pulmonary lesions in MPS than humoral-mediated immunity, such as antibody response.

Correlated responses of respiratory disease and immune capacity traits of Landrace pigs selected for Mycoplasmal pneumonia of swine (MPS) lesion.

Five generations of Landrace pigs selected for average daily gain, backfat thickness, Mycoplasmal pneumonia of swine (MPS) lesion score, and plasma cortisol levels, was executed to decrease the MPS lesion score. Genetic parameters and correlated genetic responses for respiratory disease and peripheral blood immune traits were estimated in 1395 Landrace pigs. We estimated the negative genetic correlation of MPS lesion score with phagocytic activity (PA) at 7 weeks of age (-0.67). The breeding values of PA at 7 weeks of age and 105 kg body weight and the correlated selection response of the ratio of granular leukocytes to lymphocytes at 105 kg body weight were significantly increased, and sheep red blood cell-specific antibody production (AP) was significantly decreased in a selection-dependent manner. Increasing of natural immunological indicators (e.g. PA) and decreasing of humoral immunological indicator (e.g. AP) were observed due to genetically decreasing MPS lesion score.

Effect of dietary addition of seaweed and licorice on the immune performance of pigs.

In pig production, dietary additive antibiotics are usually used for growth stimulation and disease prevention, although there is public concern about the increased incidence of resistant antibiotics and food safety. It is possible that such antibiotics might be replaced by naturally derived products such as seaweed and licorice. In this study, we evaluated the effect of dietary addition of seaweed and licorice on enhancing the immune function in swine. The animals of each group (eight animals per group) were sensitized at day 42 and 49, and the immunoglobulin production and the expression of cytokines were detected by the ELISA and real-time PCR. As the results, saliva IgA production of the seaweed-treated group increased around five times compared to that of control (day 56). Delayed hypersensitivity reaction and IgG production of the seaweed-treated group increased around 1.8-2.0 times. In addition, enhanced saliva IgA production was detected at day 50 (around two times) and day 51 (around five times) by the licorice treatment, and lower expression level of tumor necrosis factor-α messenger RNA at day 51 (around 1/25) was observed in the licorice treatment. We conclude that the replacement of antibiotics by naturally derived dietary additives might be feasible for immune system enhancement.

Chitinolytic bacteria in the intestinal tract of Japanese coastal fishes.

Intestinal bacteria from several coastal fish species were screened on 1/20 PYBG medium containing 0.2% colloidal chitin, and 361 bacteria capable of decomposing colloidal chitin were isolated. These isolates were subsequently screened on media containing either 0.5% alpha-chitin or 0.5% beta-chitin resulting in the identification of 31 alpha-chitinolytic and 275 beta-chitinolytic bacterial isolates. Partial 16S rRNA gene sequencing was carried out and homology searches of the resultant sequences against the DDBJ, EMBL, and GenBank databases revealed that the majority (99%) of the chitinolytic bacteria isolated belonged to the Vibrionaceae. Phylogenetic analysis using a Bayesian approach showed that the alpha-chitinolytic bacteria belonging to the Vibrionaceae formed a separate cluster from the non-alpha-chitinolytic bacteria in the Vibrionaceae.