Molecular Profiling of DNA Methylation and Alternative Splicing of Genes in Skeletal Muscle of Obese Rabbits.

DNA methylation and the alternative splicing of precursor messenger RNAs (pre-mRNAs) are two important genetic modification mechanisms. However, both are currently uncharacterized in the muscle metabolism of rabbits. Thus, we constructed the Tianfu black rabbit obesity model (obese rabbits fed with a 10% high-fat diet and control rabbits from 35 days to 70 days) and collected the skeletal muscle samples from the two groups for Genome methylation sequencing and RNA sequencing. DNA methylation data showed that the promoter regions of 599 genes and gene body region of 2522 genes had significantly differential methylation rates between the two groups, of which 288 genes had differential methylation rates in promoter and gene body regions. Analysis of alternative splicing showed 555 genes involved in exon skipping (ES) patterns, and 15 genes existed in differential methylation regions. Network analysis showed that 20 hub genes were associated with ubiquitinated protein degradation, muscle development pathways, and skeletal muscle energy metabolism. Our findings suggest that the two types of genetic modification have potential regulatory effects on skeletal muscle development and provide a basis for further mechanistic studies in the rabbit.

Multi-Omics Analysis of Key microRNA-mRNA Metabolic Regulatory Networks in Skeletal Muscle of Obese Rabbits.

microRNAs (miRNAs), small non-coding RNA with a length of about 22 nucleotides, are involved in the energy metabolism of skeletal muscle cells. However, their molecular mechanism of metabolism in rabbit skeletal muscle is still unclear. In this study, 16 rabbits, 8 in the control group (CON-G) and 8 in the experimental group (HFD-G), were chosen to construct an obese model induced by a high-fat diet fed from 35 to 70 days of age. Subsequently, 54 differentially expressed miRNAs, 248 differentially expressed mRNAs, and 108 differentially expressed proteins related to the metabolism of skeletal muscle were detected and analyzed with three sequencing techniques (small RNA sequencing, transcriptome sequencing, and tandem mass tab (TMT) protein technology). It was found that 12 miRNAs and 12 core genes (e.g., CRYL1, VDAC3 and APIP) were significantly different in skeletal muscle from rabbits in the two groups. The network analysis showed that seven miRNA-mRNA pairs were involved in metabolism. Importantly, two miRNAs (miR-92a-3p and miR-30a/c/d-5p) regulated three transcription factors (MYBL2, STAT1 and IKZF1) that may be essential for lipid metabolism. These results enhance our understanding of molecular mechanisms associated with rabbit skeletal muscle metabolism and provide a basis for future studies in the metabolic diseases of human obesity.

Growth, behavioural, serum biochemical and morphological changes in female rabbits fed high-fat diet.

This study aimed to determine whether high-fat diet (HFD) could cause growth, behavioural, biochemical and morphological changes in young female rabbits. Thirty-six female rabbits were randomly divided into two groups fed with either a high-fat diet (HFD) or a standard normal diet (SND) for 5 weeks. Growth and behavioural changes were recorded during the 5-week feeding period. Tissue samples, including blood and adipose tissue, were obtained after slaughter. HFD rabbits weighed more by the end of the feeding period, had a higher percent body weight and adipose tissue weight change and had longer body and bust lengths than SND rabbits. HFD rabbits significantly reduced their feed intake and feeding frequency during the fourth and fifth weeks. HFD rabbits also showed lower frequency of drinking and resting and increased stereotypical behaviour. Besides, HFD rabbits showed significant physiological abnormalities. HFD rabbits had higher serum cholesterol (TC) and triglycerides (TG) levels than SND rabbits at the end of the feeding period, and higher free fatty acid (FFA) levels than rabbits in the SND group after the third week of feeding. Serum thyroxine (T4) increased significantly in week 2 and week 5 and triiodothyronine (T3) increased significantly in week four. However, there was no significant change in serum glucose (GLU) and insulin (INS) levels. Additionally, HFD reduced the area and diameter of perirenal and subcutaneous fat cells and increased their density. Our findings suggest that HFD rabbits had higher weight gains, accumulation of fat, and more behavioural changes than SND rabbits. Although high levels of fat in the diet had a low impact on hyperglycaemia, it could lead to hyperlipidemia and hyperthyroidism. Our results also suggest that sustained HFD may cause the proliferation of adipocytes in young female rabbits.

Identification of eQTLs and sQTLs associated with meat quality in beef.

BACKGROUND: Transcription has a substantial genetic control and genetic dissection of gene expression could help us understand the genetic architecture of complex phenotypes such as meat quality in cattle. The objectives of the present research were: 1) to perform eQTL and sQTL mapping analyses for meat quality traits in longissimus dorsi muscle; 2) to uncover genes whose expression is influenced by local or distant genetic variation; 3) to identify expression and splicing hot spots; and 4) to uncover genomic regions affecting the expression of multiple genes. RESULTS: Eighty steers were selected for phenotyping, genotyping and RNA-seq evaluation. A panel of traits related to meat quality was recorded in longissimus dorsi muscle. Information on 112,042 SNPs and expression data on 8588 autosomal genes and 87,770 exons from 8467 genes were included in an expression and splicing quantitative trait loci (QTL) mapping (eQTL and sQTL, respectively). A gene, exon and isoform differential expression analysis previously carried out in this population identified 1352 genes, referred to as DEG, as explaining part of the variability associated with meat quality traits. The eQTL and sQTL mapping was performed using a linear regression model in the R package Matrix eQTL. Genotype and year of birth were included as fixed effects, and population structure was accounted for by including as a covariate the first PC from a PCA analysis on genotypic data. The identified QTLs were classified as cis or trans using 1 Mb as the maximum distance between the associated SNP and the gene being analyzed. A total of 8377 eQTLs were identified, including 75.6% trans, 10.4% cis, 12.5% DEG trans and 1.5% DEG cis; while 11,929 sQTLs were uncovered: 66.1% trans, 16.9% DEG trans, 14% cis and 3% DEG cis. Twenty-seven expression master regulators and 13 splicing master regulators were identified and were classified as membrane-associated or cytoskeletal proteins, transcription factors or DNA methylases. These genes could control the expression of other genes through cell signaling or by a direct transcriptional activation/repression mechanism. CONCLUSION: In the present analysis, we show that eQTL and sQTL mapping makes possible positional identification of gene and isoform expression regulators.

Genome wide association and gene enrichment analysis reveal membrane anchoring and structural proteins associated with meat quality in beef.

BACKGROUND: Meat quality related phenotypes are difficult and expensive to measure and predict but are ideal candidates for genomic selection if genetic markers that account for a worthwhile proportion of the phenotypic variation can be identified. The objectives of this study were: 1) to perform genome wide association analyses for Warner-Bratzler Shear Force (WBSF), marbling, cooking loss, tenderness, juiciness, connective tissue and flavor; 2) to determine enriched pathways present in each genome wide association analysis; and 3) to identify potential candidate genes with multiple quantitative trait loci (QTL) associated with meat quality. RESULTS: The WBSF, marbling and cooking loss traits were measured in longissimus dorsi muscle from 672 steers. Out of these, 495 animals were used to measure tenderness, juiciness, connective tissue and flavor by a sensory panel. All animals were genotyped for 221,077 markers and included in a genome wide association analysis. A total number of 68 genomic regions covering 52 genes were identified using the whole genome association approach; 48% of these genes encode transmembrane proteins or membrane associated molecules. Two enrichment analysis were performed: a tissue restricted gene enrichment applying a correlation analysis between raw associated single nucleotide polymorphisms (SNPs) by trait, and a functional classification analysis performed using the DAVID Bioinformatic Resources 6.8 server. The tissue restricted gene enrichment approach identified eleven pathways including "Endoplasmic reticulum membrane" that influenced multiple traits simultaneously. The DAVID functional classification analysis uncovered eleven clusters related to transmembrane or structural proteins. A gene network was constructed where the number of raw associated uncorrelated SNPs for each gene across all traits was used as a weight. A multiple SNP association analysis was performed for the top five most connected genes in the gene-trait network. The gene network identified the EVC2, ANXA10 and PKHD1 genes as potentially harboring multiple QTLs. Polymorphisms identified in structural proteins can modulate two different processes with direct effect on meat quality: in vivo myocyte cytoskeletal organization and postmortem proteolysis. CONCLUSION: The main result from the present analysis is the uncovering of several candidate genes associated with meat quality that have structural function in the skeletal muscle.

Pathway enrichment and protein interaction network analysis for milk yield, fat yield and age at first calving in a Thai multibreed dairy population.

OBJECTIVE: This research aimed to determine biological pathways and protein-protein interaction (PPI) networks for 305-d milk yield (MY), 305-d fat yield (FY), and age at first calving (AFC) in the Thai multibreed dairy population. METHODS: Genotypic information contained 75,776 imputed and actual single nucleotide polymorphisms (SNP) from 2,661 animals. Single-step genomic best linear unbiased predictions were utilized to estimate SNP genetic variances for MY, FY, and AFC. Fixed effects included herd-year-season, breed regression and heterosis regression effects. Random effects were animal additive genetic and residual. Individual SNP explaining at least 0.001% of the genetic variance for each trait were used to identify nearby genes in the National Center for Biotechnology Information database. Pathway enrichment analysis was performed. The PPI of genes were identified and visualized of the PPI network. RESULTS: Identified genes were involved in 16 enriched pathways related to MY, FY, and AFC. Most genes had two or more connections with other genes in the PPI network. Genes associated with MY, FY, and AFC based on the biological pathways and PPI were primarily involved in cellular processes. The percent of the genetic variance explained by genes in enriched pathways (303) was 2.63% for MY, 2.59% for FY, and 2.49% for AFC. Genes in the PPI network (265) explained 2.28% of the genetic variance for MY, 2.26% for FY, and 2.12% for AFC. CONCLUSION: These sets of SNP associated with genes in the set enriched pathways and the PPI network could be used as genomic selection targets in the Thai multibreed dairy population. This study should be continued both in this and other populations subject to a variety of environmental conditions because predicted SNP values will likely differ across populations subject to different environmental conditions and changes over time.

Accuracy of genomic-polygenic estimated breeding value for milk yield and fat yield in the Thai multibreed dairy population with five single nucleotide polymorphism sets

Objective: The objectives were to compare variance components, genetic parameters, prediction accuracies, and genomic-polygenic EBV rankings for milk yield (MY) and fat yield (FY) in the Thai multibreed dairy population computed using five SNP sets from GeneSeek GGP80K chip. Methods: The dataset contained monthly MY and FY of 8,361 first-lactation cows from 810 farms. Variance components, genetic parameters, and EBV for five SNP sets from the GeneSeek GGP80K chip were obtained using a 2-trait single-step average-information REML procedure. The SNP sets were the complete SNP set (all available SNP; SNP100), top 75% set (SNP75), top 50% set (SNP50), top 25% set (SNP25) and top 5% set (SNP5). The 2-trait models included herd-year-season, heterozygosity and age at first calving as fixed effects, and animal additive genetic and residual as random effects. Results: The estimates of additive genetic variances for MY and FY from SNP subsets were mostly higher than those of the complete set. The SNP25 MY and FY heritability estimates (0.276 and 0.183) were higher than those from SNP75 (0.265 and 0.168), SNP50 (0.275 and 0.179), SNP5 (0.231 and 0.169) and SNP100 (0.251and 0.159). The SNP25 EBV accuracies for MY and FY (39.76% and 33.82%) were higher than for SNP75 (35.01% and 32.60%), SNP50 (39.64% and 33.38%), SNP5 (38.61% and 29.70%) and SNP100 (34.43% and 31.61%). All rank correlations between SNP100 and SNP subsets were above 0.98 for both traits, except for SNP100 and SNP5 (0.93 for MY; 0.92 for FY). Conclusion: The high SNP25 estimates of genetic variances, heritabilities, EBV accuracies, and rank correlations between SNP100 and SNP25 for MY and FY indicated that genotyping animals with SNP25 dedicated chip would be a suitable alternative to maintain genotyping costs low while speeding up genetic progress for MY and FY in the Thai dairy population.

Modeling correlated marker effects in genome-wide prediction via Gaussian concentration graph models.

In genome-wide prediction, independence of marker allele substitution effects is typically assumed; however, since early stages in the evolution of this technology it has been known that nature points to correlated effects. In statistics, graphical models have been identified as a useful and powerful tool for covariance estimation in high dimensional problems and it is an area that has recently experienced a great expansion. In particular, Gaussian concentration graph models (GCGM) have been widely studied. These are models in which the distribution of a set of random variables, the marker effects in this case, is assumed to be Markov with respect to an undirected graph G. In this paper, Bayesian (Bayes G and Bayes G-D) and frequentist (GML-BLUP) methods adapting the theory of GCGM to genome-wide prediction were developed. Different approaches to define the graph G based on domain-specific knowledge were proposed, and two propositions and a corollary establishing conditions to find decomposable graphs were proven. These methods were implemented in small simulated and real datasets. In our simulations, scenarios where correlations among allelic substitution effects were expected to arise due to various causes were considered, and graphs were defined on the basis of physical marker positions. Results showed improvements in correlation between phenotypes and predicted additive genetic values and accuracies of predicted additive genetic values when accounting for partially correlated allele substitution effects. Extensions to the multiallelic loci case were described and some possible refinements incorporating more flexible priors in the Bayesian setting were discussed. Our models are promising because they allow incorporation of biological information in the prediction process, and because they are more flexible and general than other models accounting for correlated marker effects that have been proposed previously.

Joint genome-wide prediction in several populations accounting for randomness of genotypes: A hierarchical Bayes approach. II: Multivariate spike and slab priors for marker effects and derivation of approximate Bayes and fractional Bayes factors for the complete family of models.

This study corresponds to the second part of a companion paper devoted to the development of Bayesian multiple regression models accounting for randomness of genotypes in across population genome-wide prediction. This family of models considers heterogeneous and correlated marker effects and allelic frequencies across populations, and has the ability of considering records from non-genotyped individuals and individuals with missing genotypes in any subset of loci without the need for previous imputation, taking into account uncertainty about imputed genotypes. This paper extends this family of models by considering multivariate spike and slab conditional priors for marker allele substitution effects and contains derivations of approximate Bayes factors and fractional Bayes factors to compare models from part I and those developed here with their null versions. These null versions correspond to simpler models ignoring heterogeneity of populations, but still accounting for randomness of genotypes. For each marker loci, the spike component of priors corresponded to point mass at 0 in RS, where S is the number of populations, and the slab component was a S-variate Gaussian distribution, independent conditional priors were assumed. For the Gaussian components, covariance matrices were assumed to be either the same for all markers or different for each marker. For null models, the priors were simply univariate versions of these finite mixture distributions. Approximate algebraic expressions for Bayes factors and fractional Bayes factors were found using the Laplace approximation. Using the simulated datasets described in part I, these models were implemented and compared with models derived in part I using measures of predictive performance based on squared Pearson correlations, Deviance Information Criterion, Bayes factors, and fractional Bayes factors. The extensions presented here enlarge our family of genome-wide prediction models making it more flexible in the sense that it now offers more modeling options.

Joint genome-wide prediction in several populations accounting for randomness of genotypes: A hierarchical Bayes approach. I: Multivariate Gaussian priors for marker effects and derivation of the joint probability mass function of genotypes.

It is important to consider heterogeneity of marker effects and allelic frequencies in across population genome-wide prediction studies. Moreover, all regression models used in genome-wide prediction overlook randomness of genotypes. In this study, a family of hierarchical Bayesian models to perform across population genome-wide prediction modeling genotypes as random variables and allowing population-specific effects for each marker was developed. Models shared a common structure and differed in the priors used and the assumption about residual variances (homogeneous or heterogeneous). Randomness of genotypes was accounted for by deriving the joint probability mass function of marker genotypes conditional on allelic frequencies and pedigree information. As a consequence, these models incorporated kinship and genotypic information that not only permitted to account for heterogeneity of allelic frequencies, but also to include individuals with missing genotypes at some or all loci without the need for previous imputation. This was possible because the non-observed fraction of the design matrix was treated as an unknown model parameter. For each model, a simpler version ignoring population structure, but still accounting for randomness of genotypes was proposed. Implementation of these models and computation of some criteria for model comparison were illustrated using two simulated datasets. Theoretical and computational issues along with possible applications, extensions and refinements were discussed. Some features of the models developed in this study make them promising for genome-wide prediction, the use of information contained in the probability distribution of genotypes is perhaps the most appealing. Further studies to assess the performance of the models proposed here and also to compare them with conventional models used in genome-wide prediction are needed.

Effect of daily fluctuations in ambient temperature on reproductive failure traits of Landrace and Yorkshire sows under Thai tropical environmental conditions.

The aim of this study was to determine the effects of daily ranges and maximum ambient temperatures, and other risk factors on reproductive failure of Landrace (L) and Yorkshire (Y) sows under an open-house system in Thailand. Daily ambient temperatures were added to information on 35,579 litters from 5929 L sows and 1057 Y sows from three commercial herds. The average daily temperature ranges (ADT) and the average daily maximum temperatures (PEAK) in three gestation periods from the 35th day of gestation to parturition were classified. The considered reproductive failure traits were the occurrences of mummified fetuses (MM), stillborn piglets (STB), and piglet death losses (PDL) and an indicator trait for number of piglets born alive below the population mean (LBA). A multiple logistic regression model included farrowing herd-year-season (HYS), breed group of sow (BG), parity group (PAR), number of total piglets born (NTB), ADT1, ADT2, ADT3, PEAK1, PEAK2, and PEAK3 as fixed effects, while random effects were animal, repeated observations, and residual. Yorkshire sows had a higher occurrence of LBA than L sows (P = 0.01). The second to fifth parities sows had lower reproductive failures than other parities. The NTB regression coefficients of log-odds were positive (P < 0.01) for all traits. Narrower ranges of ADT3 increased the occurrence of MM, STB, and PDL (P < 0.01), while higher PEAK3 increased the occurrence of MM, STB, PDL, and LBA (P < 0.001). To reduce the risk of reproductive failures, particularly late in gestation, producers would need to closely monitor their temperature management strategies.

Genetic correlations between first parity and accumulated second to last parity reproduction traits as selection aids to improve sow lifetime productivity.

OBJECTIVE: The objective of this research was to estimate genetic correlations between number of piglets born alive in the first parity (NBA1), litter birth weight in the first parity (LTBW1), number of piglets weaned in the first parity (NPW1), litter weaning weight in the first parity (LTWW1), number of piglets born alive from second to last parity (NBA2+), litter birth weight from second to last parity (LTBW2+), number of piglets weaned from second to last parity (NPW2+) and litter weaning weight from second to last parity (LTWW2+), and to identify the percentages of animals (the top 10%, 25%, and 50%) for first parity and sums of second and later parity traits. METHODS: The 9,830 records consisted of 2,124 Landrace (L), 724 Yorkshire (Y), 2,650 LY, and 4,332 YL that had their first farrowing between July 1989 and December 2013. The 8-trait animal model included the fixed effects of first farrowing year-season, additive genetic group, heterosis of the sow and the litter, age at first farrowing, and days to weaning (NPW1, LTWW1, NPW2+, and LTWW2+). Random effects were animal and residual. RESULTS: Heritability estimates ranged from 0.08±0.02 (NBA1 and NPW1) to 0.29±0.02 (NPW2+). Genetic correlations between reproduction traits in the first parity and from second to last parity ranged from 0.17±0.08 (LTBW1 and LTBW2+) to 0.67±0.06 (LTWW1 and LTWW2+). Phenotypic correlations between reproduction traits in the first parity and from second to last parity were close to zero. Rank correlations between LTWW1 and LTWW2+ estimated breeding value tended to be higher than for other pairs of traits across all replacement percentages. CONCLUSION: These rank correlations indicated that selecting boars and sows using genetic predictions for first parity reproduction traits would help improve reproduction traits in the second and later parities as well as lifetime productivity in this swine population.

Imputation Accuracy from Low to Moderate Density Single Nucleotide Polymorphism Chips in a Thai Multibreed Dairy Cattle Population.

The objective of this study was to investigate the accuracy of imputation from low density (LDC) to moderate density SNP chips (MDC) in a Thai Holstein-Other multibreed dairy cattle population. Dairy cattle with complete pedigree information (n = 1,244) from 145 dairy farms were genotyped with GeneSeek GGP20K (n = 570), GGP26K (n = 540) and GGP80K (n = 134) chips. After checking for single nucleotide polymorphism (SNP) quality, 17,779 SNP markers in common between the GGP20K, GGP26K, and GGP80K were used to represent MDC. Animals were divided into two groups, a reference group (n = 912) and a test group (n = 332). The SNP markers chosen for the test group were those located in positions corresponding to GeneSeek GGP9K (n = 7,652). The LDC to MDC genotype imputation was carried out using three different software packages, namely Beagle 3.3 (population-based algorithm), FImpute 2.2 (combined family- and population-based algorithms) and Findhap 4 (combined family- and population-based algorithms). Imputation accuracies within and across chromosomes were calculated as ratios of correctly imputed SNP markers to overall imputed SNP markers. Imputation accuracy for the three software packages ranged from 76.79% to 93.94%. FImpute had higher imputation accuracy (93.94%) than Findhap (84.64%) and Beagle (76.79%). Imputation accuracies were similar and consistent across chromosomes for FImpute, but not for Findhap and Beagle. Most chromosomes that showed either high (73%) or low (80%) imputation accuracies were the same chromosomes that had above and below average linkage disequilibrium (LD; defined here as the correlation between pairs of adjacent SNP within chromosomes less than or equal to 1 Mb apart). Results indicated that FImpute was more suitable than Findhap and Beagle for genotype imputation in this Thai multibreed population. Perhaps additional increments in imputation accuracy could be achieved by increasing the completeness of pedigree information.

Estimation of Genetic Parameters and Trends for Length of Productive Life and Lifetime Production Traits in a Commercial Landrace and Yorkshire Swine Population in Northern Thailand.

The objective of this research was to estimate genetic parameters and trends for length of productive life (LPL), lifetime number of piglets born alive (LBA), lifetime number of piglets weaned (LPW), lifetime litter birth weight (LBW), and lifetime litter weaning weight (LWW) in a commercial swine farm in Northern Thailand. Data were gathered during a 24-year period from July 1989 to August 2013. A total of 3,109 phenotypic records from 2,271 Landrace (L) and 838 Yorkshire sows (Y) were analyzed. Variance and covariance components, heritabilities and correlations were estimated using an Average Information Restricted Maximum Likelihood (AIREML) procedure. The 5-trait animal model contained the fixed effects of first farrowing year-season, breed group, and age at first farrowing. Random effects were sow and residual. Estimates of heritabilities were medium for all five traits (0.17±0.04 for LPL and LBA to 0.20±0.04 for LPW). Genetic correlations among these traits were high, positive, and favorable (p<0.05), ranging from 0.93±0.02 (LPL-LWW) to 0.99±0.02 (LPL-LPW). Sow genetic trends were non-significant for LPL and all lifetime production traits. Sire genetic trends were negative and significant for LPL (-2.54±0.65 d/yr; p = 0.0007), LBA (-0.12±0.04 piglets/yr; p = 0.0073), LPW (-0.14±0.04 piglets/yr; p = 0.0037), LBW (-0.13±0.06 kg/yr; p = 0.0487), and LWW (-0.69±0.31 kg/yr; p = 0.0365). Dam genetic trends were positive, small and significant for all traits (1.04±0.42 d/yr for LPL, p = 0.0217; 0.16±0.03 piglets/yr for LBA, p<0.0001; 0.12±0.03 piglets/yr for LPW, p = 0.0002; 0.29±0.04 kg/yr for LBW, p<0.0001 and 1.23±0.19 kg/yr for LWW, p<0.0001). Thus, the selection program in this commercial herd managed to improve both LPL and lifetime productive traits in sires and dams. It was ineffective to improve LPL and lifetime productive traits in sows.

On the Bayesness, minimaxity and admissibility of point estimators of allelic frequencies.

In this paper, decision theory was used to derive Bayes and minimax decision rules to estimate allelic frequencies and to explore their admissibility. Decision rules with uniformly smallest risk usually do not exist and one approach to solve this problem is to use the Bayes principle and the minimax principle to find decision rules satisfying some general optimality criterion based on their risk functions. Two cases were considered, the simpler case of biallelic loci and the more complex case of multiallelic loci. For each locus, the sampling model was a multinomial distribution and the prior was a Beta (biallelic case) or a Dirichlet (multiallelic case) distribution. Three loss functions were considered: squared error loss (SEL), Kulback-Leibler loss (KLL) and quadratic error loss (QEL). Bayes estimators were derived under these three loss functions and were subsequently used to find minimax estimators using results from decision theory. The Bayes estimators obtained from SEL and KLL turned out to be the same. Under certain conditions, the Bayes estimator derived from QEL led to an admissible minimax estimator (which was also equal to the maximum likelihood estimator). The SEL also allowed finding admissible minimax estimators. Some estimators had uniformly smaller variance than the MLE and under suitable conditions the remaining estimators also satisfied this property. In addition to their statistical properties, the estimators derived here allow variation in allelic frequencies, which is closer to the reality of finite populations exposed to evolutionary forces.

Genetic factors influencing milk and fat yields in tropically adapted dairy cattle: insights from quantitative trait loci analysis and gene associations.

OBJECTIVE: The objective of this study was to identify genes associated with 305-day milk yield (MY) and fat yield (FY) that also influence the adaptability of the Thai multibreed dairy cattle population to tropical conditions. METHODS: A total of 75,776 imputed and actual single nucleotide polymorphisms (SNPs) from 2,661 animals were used to identify genomic regions associated with MY and FY using the single-step genomic best linear unbiased predictions. Fixed effects included herd-yearseason, breed regression, heterosis regression and calving age regression effects. Random effects were animal additive genetic and residual. Individual SNPs with a p-value smaller than 0.05 were selected for gene mapping, function analysis, and quantitative trait loci (QTL) annotation analysis. RESULTS: A substantial number of QTLs associated with MY (9,334) and FY (8,977) were identified by integrating SNP genotypes and QTL annotations. Notably, we discovered 17 annotated QTLs within the health and exterior QTL classes, corresponding to nine unique genes. Among these genes, Rho GTPase activating protein 15 (ARHGAP15) and catenin alpha 2 (CTNNA2) have previously been linked to physiological traits associated with tropical adaptation in various cattle breeds. Interestingly, these two genes also showed signs of positive selection, indicating their potential role in conferring tolerance to trypanosomiasis, a prevalent tropical disease. CONCLUSION: Our findings provide valuable insights into the genetic basis of MY and FY in the Thai multibreed dairy cattle population, shedding light on the underlying mechanisms of tropical adaptation. The identified genes represent promising targets for future breeding strategies aimed at improving milk and fat production while ensuring resilience to tropical challenges. This study significantly contributes to our understanding of the genetic factors influencing milk production and adaptability in dairy cattle, facilitating the development of sustainable genetic selection strategies and breeding programs in tropical environments.

A comparison of five sets of overlapping and non-overlapping sliding windows for semen production traits in the Thai multibreed dairy population.

OBJECTIVE: This study compared five distinct sets of biological pathways and associated genes related to semen volume (VOL), number of sperm (NS), and sperm motility (MOT) in the Thai multibreed dairy population. METHODS: The phenotypic data included 13,533 VOL records, 12,773 NS records, and 12,660 MOT records from 131 bulls. The genotypic data consisted of 76,519 imputed and actual single nucleotide polymorphisms (SNPs) from 72 animals. The SNP additive genetic variances for VOL, NS, and MOT were estimated for SNP windows of one SNP (SW1), ten SNP (SW10), 30 SNP (SW30), 50 SNP (SW50), and 100 SNP (SW100) using a single-step genomic best linear unbiased prediction approach. The fixed effects in the model were contemporary group, ejaculate order, bull age, ambient temperature, and heterosis. The random effects accounted for animal additive genetic effects, permanent environment effects, and residual. The SNPs explaining at least 0.001% of the additive genetic variance in SW1, 0.01% in SW10, 0.03% in SW30, 0.05% in SW50, and 0.1% in SW100 were selected for gene identification through the NCBI database. The pathway analysis utilized genes associated with the identified SNP windows. RESULTS: Comparison of overlapping and non-overlapping SNP windows revealed notable differences among the identified pathways and genes associated with the studied traits. Overlapping windows consistently yielded a larger number of shared biological pathways and genes than non-overlapping windows. In particular, overlapping SW30 and SW50 identified the largest number of shared pathways and genes in the Thai multibreed dairy population. CONCLUSION: This study yielded valuable insights into the genetic architecture of VOL, NS, and MOT. It also highlighted the importance of assessing overlapping and non-overlapping SNP windows of various sizes for their effectiveness to identify shared pathways and genes influencing multiple traits.

Influence of proportion of Brahman genetics on productivity of Brahman-Angus cows at weaning.

This study evaluated the association between the proportion of Brahman genetics and productivity of Brahman-Angus cows at weaning using a 31-yr dataset containing 6,312 cows and 5,405 pregnancies. Cows were contemporaneously reared and enrolled in yearly breeding seasons under subtropical conditions of North-Central Florida. They were evenly distributed in six-breed groups (G) according to the proportion of Brahman genetics: G0% to 19%, G21% to 34%, G38% (Brangus), G41% to 59%, G63% to 78%, and G81% to 100%. The proportion of cows calving (84.9%) did not differ across the six-breed groups. However, cows in the G81% to 100% weaned fewer calves (90.8%) than cows in the G0% to 19% and G21% to 34% (95.7%, each). The weaning rate of cows in the G38% (94.3%), G41% to 59% (94.2%), and G63% to 78% (93.0%) was intermediate between these three breed groups. The preweaning calf mortality was greater for cows in the G81% to 100% (9.2%) than cows in the G0% to 19% and G21% to 34% (4.3%, each), but intermediate for cows in the G38% (5.7%), G41% to 59% (5.8%), and G63% to 78% (7.0%). Cows in the G81% to 100% also weaned lighter calves (220.6 kg) than cows in the G0% to 19% (245.2 kg), G21% to 34% (250.2 kg), G38% (247.9 kg), G41% to 59% (252.5 kg), and G63% to 78% (245.2 kg). Cows in the G0% to 19% weaned lighter calves than cows with 21% to 78% of Brahman genetics. The 205-d adjusted weaning weight evidenced the less productive results of cows in G0% to 19% and G81% to 100% compared with other genetic groups, as they calved at the fastest and slowest rate, respectively. Thus, the 205-d adjusted weaning weight eliminated this bias. Additionally, younger cows weaned lighter calves; and male calves were heavier at weaning than female calves. Both parity order of cow and calf sex altered the magnitude of the described association between breed group of cows and calf weaning weights. Overall, after adjusting for weaning rate and age of calves at weaning, the number of kilograms produced per cow submitted to reproduction was less for cows in the G0% to 19% (191.1 kg) and G81% to 100 (181.8 kg) compared with cows in the G21% to 34 (197.0 kg), G38 (195.9 kg), G41% to 59 (199.7), and G63% to 78 (196.2). Cows in the G81% to 100% were the least productive. Thus, a proportion of Brahman genetics between 21% and 78% ensured superior productivity of Brahman-Angus cows subjected to subtropical conditions.

Variance components and genetic parameters for milk production and lactation pattern in an ethiopian multibreed dairy cattle population.

The objective of this study was to estimate variance components and genetic parameters for lactation milk yield (LY), lactation length (LL), average milk yield per day (YD), initial milk yield (IY), peak milk yield (PY), days to peak (DP) and parameters (ln(a) and c) of the modified incomplete gamma function (MIG) in an Ethiopian multibreed dairy cattle population. The dataset was composed of 5,507 lactation records collected from 1,639 cows in three locations (Bako, Debre Zeit and Holetta) in Ethiopia from 1977 to 2010. Parameters for MIG were obtained from regression analysis of monthly test-day milk data on days in milk. The cows were purebred (Bos indicus) Boran (B) and Horro (H) and their crosses with different fractions of Friesian (F), Jersey (J) and Simmental (S). There were 23 breed groups (B, H, and their crossbreds with F, J, and S) in the population. Fixed and mixed models were used to analyse the data. The fixed model considered herd-year-season, parity and breed group as fixed effects, and residual as random. The single and two-traits mixed animal repeatability models, considered the fixed effects of herd-year-season and parity subclasses, breed as a function of cow H, F, J, and S breed fractions and general heterosis as a function of heterozygosity, and the random additive animal, permanent environment, and residual effects. For the analysis of LY, LL was added as a fixed covariate to all models. Variance components and genetic parameters were estimated using average information restricted maximum likelihood procedures. The results indicated that all traits were affected (p<0.001) by the considered fixed effects. High grade B×F cows (3/16B 13/16F) had the highest least squares means (LSM) for LY (2,490±178.9 kg), IY (10.5±0.8 kg), PY (12.7±0.9 kg), YD (7.6±0.55 kg) and LL (361.4±31.2 d), while B cows had the lowest LSM values for these traits. The LSM of LY, IY, YD, and PY tended to increase from the first to the fifth parity. Single-trait analyses yielded low heritability (0.03±0.03 and 0.08±0.02) and repeatability (0.14±0.01 to 0.24±0.02) estimates for LL, DP and parameter c. Medium heritability (0.21±0.03 to 0.33±0.04) and repeatability (0.27±0.02 to 0.53±0.01) estimates were obtained for LY, IY, PY, YD and ln(a). Genetic correlations between LY, IY, PY, YD, ln(a), and LL ranged from 0.59 to 0.99. Spearman's rank correlations between sire estimated breeding values for LY, LL, IY, PY, YD, ln(a) and c were positive (0.67 to 0.99, p<0.001). These results suggested that selection for IY, PY, YD, or LY would genetically improve lactation milk yield in this Ethiopian dairy cattle population.

Estimation of heritabilities and additive genetic correlations for reproduction traits in swine: insights for tropical commercial production systems using multiple trait animal models.

OBJECTIVE: This study was to estimate heritabilities, additive genetic correlations, and phenotypic correlations between number of piglets born alive (NBA), litter birth weight (LTBW), number of piglets weaned (NPW) and litter weaning weight (LTWW) in different parities of Landrace (L), Yorkshire (Y), Landrace×Yorkshire (LY), and Yorkshire×Landrace (YL) sows in a commercial swine operation in Northern Thailand. METHODS: Two models were utilized, a single trait repeatability model (RM) and a multiple trait animal model (MTM). The RM assumed reproductive records from different parities to be repeated values of the same trait, whereas the MTM assumed these records to be different traits. The two models accounted for the fixed effects of farrowing year-season, genetic group of the sow, heterosis, and age at first farrowing, and the random effects of sow, boar, and residual. RESULTS: Heritability estimates from RM were 0.02±0.01 for NBA, 0.10±0.01 for LTBW, 0.04±0.01 for NPW, and 0.11±0.01 for LTWW. Heritability estimates from MTM fluctuated across parities, ranging from 0.04±0.01 in parity 2 to 0.09±0.02 in parity 4 for NBA, 0.07±0.02 in parity 2 to 0.16±0.02 in parity 3 for LTBW, 0.04±0.02 in parity 4 to 0.08±0.01 in parity 1 for NPW, and 0.16±0.02 in parity 1 to 0.20±0.02 in parity 2 for LTWW. Additive genetic correlation estimates from MTM were also variable, ranging from 0.29±0.24 between NBA in parity 1 and NBA in parity 2 to 0.99±0.05 between LTWW in parity 3 and LTWW in parity 4. CONCLUSION: The findings of this study highlight the advantage of using MTM for the genetic improvement of reproductive traits in swine and contribute to the development of sustainable swine breeding programs in Thailand.