Generalizability Theory in the Evaluation of Psychological Profile in Track and Field

Generalizability theory (GT) has been used throughout the scientific literature to ensure validity, reliability, and generalizability in different sport contexts. However, there is a small number of studies examining the measurement of psychological profiles in sport from this perspective. Therefore, this study's main goal is the sources of variability and the optimal measurement design estimation for a good assessment of the psychological profile in track and field. The sample consisted of 470 participants (age: Average= 32.1; Standar Desviation = 13.5). The analysis of variance and generalizability component analysis has been performed in order to test the reliability and generalizability of the sample. The profile included the following variables: flow, motivation (from Self-Determination Theory and Achievement Goals), self-confidence, and psychological skills. Results confirm that the sample has a high degree of reliability and generalizability in all the tested models. So, a detailed study on the validity, reliability, and generalizability of samples and measures should be an inherent element in the practice of psychological counseling in sports.

Genetic mapping reveals new loci and alleles for flowering time and plant height using the double round-robin population of barley

Flowering time and plant height are two critical determinants of yield potential in barley (Hordeum vulgare). Despite their role in plant physiological regulation, a complete overview of the genetic complexity of flowering time and plant height regulation in barley is still lacking. Using a double round-robin population originated from the crossings of 23 diverse parental inbred lines, we aimed to determine the variance components in the regulation of flowering time and plant height in barley as well as to identify new genetic variants by single and multi-population QTL analyses and allele mining. Despite similar genotypic variance, we observed higher environmental variance components for plant height than flowering time. Furthermore, we detected new QTLs for flowering time and plant height. Finally, we identified a new functional allelic variant of the main regulatory gene Ppd-H1. Our results show that the genetic architecture of flowering time and plant height might be more complex than reported earlier and that a number of undetected, small effect, or low-frequency genetic variants underlie the control of these two traits.

Genetic parameters for calf feeding traits derived from automated milk feeding machines and number of bovine respiratory disease treatments in North American Holstein calves

Precision livestock farming technologies, such as automatic milk feeding machines, have increased the availability of on-farm data collected from dairy operations. We analyzed feeding records from automatic milk feeding machines to evaluate the genetic background of milk feeding traits and bovine respiratory disease (BRD) in North American Holstein calves. Data from 10,076 preweaning female Holstein calves were collected daily over a period of 6 yr (3 yr included per-visit data), and daily milk consumption (DMC), per-visit milk consumption (PVMC), daily sum of drinking duration (DSDD), drinking duration per-visit, daily number of rewarded visits (DNRV), and total number of visits per day were recorded over a 60-d preweaning period. Additional traits were derived from these variables, including total consumption and duration variance (TCV and TDV), feeding interval, drinking speed (DS), and preweaning stayability. A single BRD-related trait was evaluated, which was the number of times a calf was treated for BRD (NTT). The NTT was determined by counting the number of BRD incidences before 60 d of age. All traits were analyzed using single-step genomic BLUP mixed-model equations and fitting either repeatability or random regression models in the BLUPF90+ suite of programs. A total of 10,076 calves with phenotypic records and genotypic information for 57,019 SNP after the quality control were included in the analyses. Feeding traits had low heritability estimates based on repeatability models (0.006 ± 0.0009 to 0.08 ± 0.004). However, total variance traits using an animal model had greater heritabilities of 0.21 ± 0.023 and 0.23 ± 0.024, for TCV and TDV, respectively. The heritability estimates increased with the repeatability model when using only the first 32 d preweaning (e.g., PVMC = 0.040 ± 0.003, DMC = 0.090 ± 0.009, DSDD = 0.100 ± 0.005, DS = 0.150 ± 0.007, DNRV = 0.020 ± 0.002). When fitting random regression models (RRM) using the full dataset (60-d period), greater heritability estimates were obtained (e.g., PVMC = 0.070 [range: 0.020, 0.110], DMC = 0.460 [range: 0.050, 0.680], DSDD = 0.180 [range: 0.010, 0.340], DS = 0.19 [range: 0.070, 0.430], DNRV = 0.120 [range: 0.030, 0.450]) for the majority of the traits, suggesting that RRM capture more genetic variability than the repeatability model with better fit being found for RRM. Moderate negative genetic correlations of -0.59 between DMC and NTT were observed, suggesting that automatic milk feeding machines records have the potential to be used for genetically improving disease resilience in Holstein calves. The results from this study provide key insights of the genetic background of early in-life traits in dairy cattle, which can be used for selecting animals with improved health outcomes and performance.

Estimation of heritability with genomic information by method R

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.

A boosting method to select the random effects in linear mixed models

This paper proposes a novel likelihood-based boosting method for the selection of the random effects in linear mixed models. The nonconvexity of the objective function to minimize, which is the negative profile log-likelihood, requires the adoption of new solutions. In this respect, our optimization approach also employs the directions of negative curvature besides the usual Newton directions. A simulation study and a real-data application show the good performance of the proposal.

Long-read sequencing-based transcriptomic landscape in longissimus dorsi and transcriptome-wide association studies for growth traits of meat rabbits

Rabbits are an attractive meat livestock species that can efficiently convert human-indigestible plant biomass, and have been commonly used in biological and medical researches. Yet, transcriptomic landscape in muscle tissue and association between gene expression level and growth traits have not been specially studied in meat rabbits. In this study Oxford Nanopore Technologies (ONT) long-read sequencing technology was used for comprehensively exploring transcriptomic landscape in Longissimus dorsi for 115 rabbits at 84 days of age, and transcriptome-wide association studies (TWAS) were performed for growth traits, including body weight at 84 days of age and average daily gain during three growth periods. The statistical analysis of TWAS was performed using a mixed linear model, in which polygenic effect was fitted as a random effect according to gene expression level-based relationships. A total of 18,842 genes and 42,010 transcripts were detected, among which 35% of genes and 47% of transcripts were novel in comparison with the reference genome annotation. Furthermore, 45% of genes were widely expressed among more than 90% of individuals. The proportions (±SE) of phenotype variance explained by genome-wide gene expression level ranged from 0.501 ± 0.216 to 0.956 ± 0.209, and the similar results were obtained when explained by transcript expression level. In contrast, neither gene nor transcript was detected by TWAS to be statistically significantly associated with these growth traits. In conclusion, these novel genes and transcripts that have been extensively profiled in a single muscle tissue using long-read sequencing technology will greatly improve our understanding on transcriptional diversity in rabbits. Our results with a relatively small sample size further revealed the important contribution of global gene expression to phenotypic variation on growth performance, but it seemed that no single gene has an outstanding effect; this knowledge is helpful to include intermediate omics data for implementing genetic evaluation of growth traits in meat rabbits.

Genetic parameters for yearling male reproduction traits in tropical composite cattle population

Fertility is economically important but is hard to quantify and measure in breeding programs which has led extensive breeding programs to ignore fertility in their selection criteria. While female fertility traits have been extensively researched, male fertility traits have been largely ignored. It is estimated that 20% to 40% of bulls have sub-fertility, reducing the number of calves born and profits, highlighting the importance of investigating bull fertility. The most practical measure of male fertility is a bull breeding soundness evaluation (BBSE) which assesses structure as well as semen quality and quantity. Generally, traits recorded in a BBSE are neither genetically evaluated nor used for selection in breeding programs. All traits recorded during a BBSE were analyzed through a series of univariate and bivariate linear mixed models using a genomic relationship matrix to estimate genetic parameters. All genotype and phenotype data were obtained from a tropical composite commercial cattle population and imputed to 27,638 single-nucleotide polymorphisms (SNPs) with a total of 2,613 genotyped animals with BBSE records ranging from 616 to 826 animals depending on the trait. The heritabilities of the 27 traits recorded during a BBSE ranged from 0.02 to 0.49. Seven of the male fertility traits were recommended to be included in a breeding program based on their heritability and their phenotypic and genetic correlations. These traits are scrotal circumference, percent normal sperm, proximal droplets, distal midpiece reflex, knobbed acrosomes, vacuoles/teratoids, and sheath score. Using these seven traits in a breeding program would result in higher calving rates, increasing production and profitability.

Partitioning variance in reproductive success, within years and across lifetimes

Variance in reproductive success (s k 2 , with k = number of offspring) plays a large role in determining the rate of genetic drift and the scope within which selection acts. Various frameworks have been proposed to parse factors that contribute to s k 2 , but none has focused on age-specific values of ϕ = s k 2 / k ¯ , which indicate the degree to which reproductive skew is overdispersed (compared to the random Poisson expectation) among individuals of the same age and sex. Instead, within-age effects are generally lumped with residual variance and treated as "noise." Here, an ANOVA sums-of-squares framework is used to partition variance in annual and lifetime reproductive success into between-group and within-group components. For annual reproduction, the between-age effect depends on age-specific fecundity (b x), but relatively few empirical data are available on the within-age effect, which depends on ϕ x. By defining groups by age-at-death rather than age, the same ANOVA framework can be used to partition variance in lifetime reproductive success (LRS) into between-group and within-group components. Analytical methods are used to develop null-model expectations for random contributions to within-group and between-group components. For analysis of LRS, random variation in longevity appears as part of the between-group variance, and effects (if any) of skip breeding and persistent individual differences contribute to the within-group variance. Simulations are used to show that the methods for variance partitioning are asymptotically unbiased. Practical application is illustrated with empirical data for annual reproduction in American black bears and lifetime reproduction in Dutch great tits. Results show that overdispersed within-age variance (1) dominates annual s k 2 in both male and female black bears, (2) is the primary factor that reduces annual effective size to a fraction of the number of adults, and (3) represents most of the opportunity for selection. In contrast, about a quarter of the variance in LRS in great tits can be attributed to random variation in longevity, and most of the rest is due to modest differences in fecundity with age estimated for a single cohort of females. R code is provided that reads generic input files for annual and lifetime reproductive success and allows users to conduct variance partitioning with their own data.

Genetic Variance Estimation over Time in Broiler Breeding Programmes for Growth and Reproductive Traits

Monitoring the genetic variance of traits is a key priority to ensure the sustainability of breeding programmes in populations under directional selection, since directional selection can decrease genetic variation over time. Studies monitoring changes in genetic variation have typically used long-term data from small experimental populations selected for a handful of traits. Here, we used a large dataset from a commercial breeding line spread over a period of twenty-three years. A total of 2,059,869 records and 2,062,112 animals in the pedigree were used for the estimations of variance components for the traits: body weight (BWT; 2,059,869 records) and hen-housed egg production (HHP; 45,939 records). Data were analysed with three estimation approaches: sliding overlapping windows, under frequentist (restricted maximum likelihood (REML)) and Bayesian (Gibbs sampling) methods; expected variances using coefficients of the full relationship matrix; and a "double trait covariances" analysis by computing correlations and covariances between the same trait in two distinct consecutive windows. The genetic variance showed marginal fluctuations in its estimation over time. Whereas genetic, maternal permanent environmental, and residual variances were similar for BWT in both the REML and Gibbs methods, variance components when using the Gibbs method for HHP were smaller than the variances estimated when using REML. Large data amounts were needed to estimate variance components and detect their changes. For Gibbs (REML), the changes in genetic variance from 1999-2001 to 2020-2022 were 82.29 to 93.75 (82.84 to 93.68) for BWT and 76.68 to 95.67 (98.42 to 109.04) for HHP. Heritability presented a similar pattern as the genetic variance estimation, changing from 0.32 to 0.36 (0.32 to 0.36) for BWT and 0.16 to 0.15 (0.21 to 0.18) for HHP. On the whole, genetic parameters tended slightly to increase over time. The expected variance estimates were lower than the estimates when using overlapping windows. That indicates the low effect of the drift-selection process on the genetic variance, or likely, the presence of genetic variation sources compensating for the loss. Double trait covariance analysis confirmed the maintenance of variances over time, presenting genetic correlations >0.86 for BWT and >0.82 for HHP. Monitoring genetic variance in broiler breeding programmes is important to sustain genetic progress. Although the genetic variances of both traits fluctuated over time, in some windows, particularly between 2003 and 2020, increasing trends were observed, which warrants further research on the impact of other factors, such as novel mutations, operating on the dynamics of genetic variance.

Separating Measurement Error and Signal in Environmental Data: Use of Replicates to Address Uncertainty

Measurement uncertainty has long been a concern in the characterizing and interpreting environmental and toxicological measurements. We compared statistical analysis approaches when there are replicates: a Naı̈ve approach that omits replicates, a Hybrid approach that inappropriately treats replicates as independent samples, and a Measurement Error Model (MEM) approach in a random effects analysis of variance (ANOVA) model that appropriately incorporates replicates. A simulation study assessed the effects of sample size and levels of replication, signal variance, and measurement error on estimates from the three statistical approaches. MEM results were superior overall with confidence intervals for the observed mean narrower on average than those from the Naı̈ve approach, giving improved characterization. The MEM approach also featured an unparalleled advantage in estimating signal and measurement error variance separately, directly addressing measurement uncertainty. These MEM estimates were approximately unbiased on average with more replication and larger sample sizes. Case studies illustrated analyzing normally distributed arsenic and log-normally distributed chromium concentrations in tap water and calculating MEM confidence intervals for the true, latent signal mean and latent signal geometric mean (i.e., with measurement error removed). MEM estimates are valuable for study planning; we used simulation to compare various sample sizes and levels of replication.

Variance of age-specific log incidence decomposition (VALID): a unifying model of measured and unmeasured genetic and non-genetic risks

BACKGROUND

The extent to which known and unknown factors explain how much people of the same age differ in disease risk is fundamental to epidemiology. Risk factors can be correlated in relatives, so familial aspects of risk (genetic and non-genetic) must be considered.

DEVELOPMENT

We present a unifying model (VALID) for variance in risk, with risk defined as log(incidence) or logit(cumulative incidence). Consider a normally distributed risk score with incidence increasing exponentially as the risk increases. VALID's building block is variance in risk, Δ2, where Δ = log(OPERA) is the difference in mean between cases and controls and OPERA is the odds ratio per standard deviation. A risk score correlated r between a pair of relatives generates a familial odds ratio of exp(rΔ2). Familial risk ratios, therefore, can be converted into variance components of risk, extending Fisher's classic decomposition of familial variation to binary traits. Under VALID, there is a natural upper limit to variance in risk caused by genetic factors, determined by the familial odds ratio for genetically identical twin pairs, but not to variation caused by non-genetic factors.

APPLICATION

For female breast cancer, VALID quantified how much variance in risk is explained-at different ages-by known and unknown major genes and polygenes, non-genomic risk factors correlated in relatives, and known individual-specific factors.

CONCLUSION

VALID has shown that, while substantial genetic risk factors have been discovered, much is unknown about genetic and familial aspects of breast cancer risk especially for young women, and little is known about individual-specific variance in risk.

Variance component estimates, phenotypic characterization, and genetic evaluation of bovine congestive heart failure in commercial feeder cattle

The increasing incidence of bovine congestive heart failure (BCHF) in feedlot cattle poses a significant challenge to the beef industry from economic loss, reduced performance, and reduced animal welfare attributed to cardiac insufficiency. Changes to cardiac morphology as well as abnormal pulmonary arterial pressure (PAP) in cattle of mostly Angus ancestry have been recently characterized. However, congestive heart failure affecting cattle late in the feeding period has been an increasing problem and tools are needed for the industry to address the rate of mortality in the feedlot for multiple breeds. At harvest, a population of 32,763 commercial fed cattle were phenotyped for cardiac morphology with associated production data collected from feedlot processing to harvest at a single feedlot and packing plant in the Pacific Northwest. A sub-population of 5,001 individuals were selected for low-pass genotyping to estimate variance components and genetic correlations between heart score and the production traits observed during the feeding period. At harvest, the incidence of a heart score of 4 or 5 in this population was approximately 4.14%, indicating a significant proportion of feeder cattle are at risk of cardiac mortality before harvest. Heart scores were also significantly and positively correlated with the percentage Angus ancestry observed by genomic breed percentage analysis. The heritability of heart score measured as a binary (scores 1 and 2 = 0, scores 4 and 5 = 1) trait was 0.356 in this population, which indicates development of a selection tool to reduce the risk of congestive heart failure as an EPD (expected progeny difference) is feasible. Genetic correlations of heart score with growth traits and feed intake were moderate and positive (0.289-0.460). Genetic correlations between heart score and backfat and marbling score were -0.120 and -0.108, respectively. Significant genetic correlation to traits of high economic importance in existing selection indexes explain the increased rate of congestive heart failure observed over time. These results indicate potential to implement heart score observed at harvest as a phenotype under selection in genetic evaluation in order to reduce feedlot mortality due to cardiac insufficiency and improve overall cardiopulmonary health in feeder cattle.

Analysis of the Genetic Parameters for Dairy Linear Appraisal and Zoometric Traits: A Tool to Enhance the Applicability of Murciano-Granadina Goats Major Areas Evaluation System

Selection for zoometrics defines individuals' productive longevity, endurance, enhanced productive abilities and consequently, their long-term profitability. When zoometric analysis is aimed at large highly selected populations or in those at different levels of selection, linear appraisal systems (LAS) provide a timely response. This study estimates genetic and phenotypic parameters for zoometric/LAS traits in Murciano-Granadina goats, estimating genetic and phenotypic correlations among all traits, and determining whether major area selection would be appropriate or if adaptability strategies may need to be followed. Heritability estimates for the zoometric/LAS traits were low to high, ranging from 0.09 to 0.43, and the accuracy of estimation has improved after decades, rendering standard errors negligible. Scale inversion of specific traits may need to be performed before major areas selection strategies are implemented. Genetic and phenotypic correlations suggests that negative selection against thicker bones and higher rear insertion heights indirectly results in the optimization of selection practices in the rest of the traits, especially those in the structure, capacity and mammary system major areas. The integration and implementation of the strategies proposed within the Murciano-Granadina breeding program maximizes selection opportunities and the sustainable international competitiveness of the Murciano-Granadina goat in the dairy goat breed panorama.

Temporal Trends in Variability of Respirable Dust and Respirable Quartz Concentrations in the European Industrial Minerals Sector

While between- and within-worker variability have been studied quite extensively, hardly any research is available that examines long-term trends in the variability of occupational exposure. In this first study on trends in occupational exposure variability temporal changes in the variability of respirable dust and respirable quartz concentrations within the European industrial minerals sector were demonstrated. Since 2000 the European Industrial Minerals Association's Dust Monitoring Program (IMA-DMP) has systematically collected respirable dust and respirable quartz measurements. The resulting IMA-DMP occupational exposure database contains at present approximately 40 000 personal full-shift measurements, collected at 177 sites owned by 39 companies, located in 23 European countries. Repeated measurements of workers performing their duties within a specific site-job-campaign combination allowed estimation of within- and between-worker variability in exposure concentrations. Overall day-to-day variability predominated the between-worker variability for both respirable dust concentrations and quartz concentrations. The within-worker variability in concentrations by job was two to three times higher for respirable quartz than for respirable dust. The median between-worker variability in respirable dust concentrations was low and further reduced over time. For quartz concentrations the same phenomenon albeit somewhat less strong was observed. In contrast, for the within-worker variability in concentrations downward and upward temporal trends were apparent for both respirable dust and respirable quartz. The study shows that the (relative) size of temporal variability is large and unpredictable and therefore regular measurement campaigns are needed to ascertain compliance to occupational exposure limit values.

Phenotypic Analysis of Fourier-Transform Infrared Milk Spectra in Dairy Goats

The infrared spectrum of bovine milk is used to predict many interesting traits, whereas there have been few studies on goat milk in this regard. The objective of this study was to characterize the major sources of variation in the absorbance of the infrared spectrum in caprine milk samples. A total of 657 goats belonging to 6 breeds and reared on 20 farms under traditional and modern dairy systems were milk-sampled once. Fourier-transform infrared (FTIR) spectra were taken (2 replicates per sample, 1314 spectra), and each spectrum contained absorbance values at 1060 different wavenumbers (5000 to 930 × cm^-1), which were treated as a response variable and analyzed one at a time (i.e., 1060 runs). A mixed model, including the random effects of sample/goat, breed, flock, parity, stage of lactation, and the residual, was used. The pattern and variability of the FTIR spectrum of caprine milk was similar to those of bovine milk. The major sources of variation in the entire spectrum were as follows: sample/goat (33% of the total variance); flock (21%); breed (15%); lactation stage (11%); parity (9%); and the residual unexplained variation (10%). The entire spectrum was segmented into five relatively homogeneous regions. Two of them exhibited very large variations, especially the residual variation. These regions are known to be affected by the absorbance of water, although they also exhibited wide variations in the other sources of variation. The average repeatability of these two regions were 45% and 75%, whereas for the other three regions it was about 99%. The FTIR spectrum of caprine milk could probably be used to predict several traits and to authenticate the origin of goat milk.

Study of genetic parameters for pre-weaning growth traits in inner Mongolia white Arbas cashmere goats

Inner Mongolia Arbas white cashmere goats is a dual-purpose breed for producing cashmere and meat. In recent years, its meat has become more and more popular among consumers because of rich nutrients and delicious flavor. Therefore, it is particularly important to study the genetic and non-genetic factors affecting the early growth traits and estimate variance components of pre-weaning growth traits of Inner Mongolia Albas white cashmere goats. A total of 37487 kidding records such as birth weight (BWT), weaning weight (WWT), average daily gain from birth to weaning (ADG) and Kleiber ratio (KR) from 343 sires and 7296 dams were used in this study. The most appropriate model was chosen on the basis of likelihood ratio test by fitting six models which excluding or including maternal genetic, maternal permanent environmental effects. The parameters were estimated under the most appropriate model using AIREML method by WOMBAT software. With the best model (Model 6), heritability estimates were 0.0435, 0.0911, 0.0932 and 0.2339 for BWT, WWT, ADG and KR traits, respectively. Maternal heritability estimates were 0.0143, 0.0246, 0.0220 and 0.0186 for BWT, WWT, ADG, and KR traits respectively. The correlation between different traits was estimated with the most suitable model by using bivariate analysis method. The direct additive genetic correlation among the traits ranged from -0.026 (BWT~KR) to 0.772 (ADG~KR). The maternal permanent environment correlation is between -0.289 (BWT-KR) ~0.900 (WWT-ADG). Results indicated that maternal effects and direct-maternal genetic covariance should be considered in any program aimed at improving pre-weaning growth traits to have an accurate genetic evaluation. In addition, positive and medium to high genetic correlations generally exist among WWT, ADG and KR due to the existence of genetic variation for early growth traits. The results showed that the genetic progress of these traits could be slowly through selection except for KR.

Genetic parameters of immune traits for Landrace and Large White pig breeds

Improving the immunocompetence towards pathogens represents a desirable objective of breeding strategies to increase resilience. However, the immune system is complex and the genetic foundation of the underlying components is not yet clarified. In the present study, we focused on 22 blood parameters of 1,144 Landrace (LR) and Large White (LW) piglets at the age of 6-7 weeks. The immune profiles covered immune cells, red blood cell characteristics and cytokines. Genetic parameters based on pedigree information along with possible environmental effects were estimated. Litter effects play an important role in the expression of immune parameters of their young progenies. Hence, litter impacts on the piglet's immune profile including the immune parameters of the dam itself were investigated by different models. To incorporate the complexity of the immune network, the data were further investigated with a principal component analysis. Immune traits showed low to high breed-specific heritabilities (h² ). Strong positive r_g were estimated among red blood cell characteristics (0.77-0.99) and among cytokines (0.48-0.99). Neutrophils and lymphocytes illustrated a high negative r_g (-0.96 to -0.98). The litter impact on piglet's immunity was examined and strengthened already observed breed differences. In LR, h² (0.22-0.15) and litter effect (c² ) (0.52-0.44) for IFN-γ decreased after statistical consideration of maternal impact. In LW, a decrease in h² (0.32-0.18) for IFN-γ and an increase in c² (0.54-0.56) were observed. Here, sufficient correlations were detected within various immune traits and functional biological networks of principal components. Most immune traits are heritable and are promising to cover global breed-specific immunocompetence in pigs. The analysis of immune traits has to be extended in order to find an optimal range and to characterize relationships between immunity and performance to gain an improved immune system without accidental losses in productivity.

Is single-step genomic REML with the algorithm for proven and young more computationally efficient when less generations of data are present?

Efficient computing techniques allow the estimation of variance components for virtually any traditional dataset. When genomic information is available, variance components can be estimated using genomic REML (GREML). If only a portion of the animals have genotypes, single-step GREML (ssGREML) is the method of choice. The genomic relationship matrix (G) used in both cases is dense, limiting computations depending on the number of genotyped animals. The algorithm for proven and young (APY) can be used to create a sparse inverse of G (GAPY~-1) with close to linear memory and computing requirements. In ssGREML, the inverse of the realized relationship matrix (H-1) also includes the inverse of the pedigree relationship matrix, which can be dense with a long pedigree, but sparser with short. The main purpose of this study was to investigate whether costs of ssGREML can be reduced using APY with truncated pedigree and phenotypes. We also investigated the impact of truncation on variance components estimation when different numbers of core animals are used in APY. Simulations included 150K animals from 10 generations, with selection. Phenotypes (h2 = 0.3) were available for all animals in generations 1-9. A total of 30K animals in generations 8 and 9, and 15K validation animals in generation 10 were genotyped for 52,890 SNP. Average information REML and ssGREML with G-1 and GAPY~-1 using 1K, 5K, 9K, and 14K core animals were compared. Variance components are impacted when the core group in APY represents the number of eigenvalues explaining a small fraction of the total variation in G. The most time-consuming operation was the inversion of G, with more than 50% of the total time. Next, numerical factorization consumed nearly 30% of the total computing time. On average, a 7% decrease in the computing time for ordering was observed by removing each generation of data. APY can be successfully applied to create the inverse of the genomic relationship matrix used in ssGREML for estimating variance components. To ensure reliable variance component estimation, it is important to use a core size that corresponds to the number of largest eigenvalues explaining around 98% of total variation in G. When APY is used, pedigrees can be truncated to increase the sparsity of H and slightly reduce computing time for ordering and symbolic factorization, with no impact on the estimates.

Concordance correlation coefficients estimated by modified variance components and generalized estimating equations for longitudinal overdispersed Poisson data

For longitudinal overdispersed Poisson data sets, estimators of the intra-, inter-, and total concordance correlation coefficient through variance components have been proposed. However, biased estimators of quadratic forms are used in concordance correlation coefficient estimation. In addition, the generalized estimating equations approach has been used in estimating agreement for longitudinal normal data and not for longitudinal overdispersed Poisson data. Therefore, this paper proposes a modified variance component approach to develop the unbiased estimators of the concordance correlation coefficient for longitudinal overdispersed Poisson data. Further, the indices of intra-, inter-, and total agreement through generalized estimating equations are also developed considering the correlation structure of longitudinal count repeated measurements. Simulation studies are conducted to compare the performance of the modified variance component and generalized estimating equation approaches for longitudinal Poisson and overdispersed Poisson data sets. An application of corticospinal diffusion tensor tractography study is used for illustration. In conclusion, the modified variance component approach performs outstandingly well with small mean square errors and nominal 95% coverage rates. The generalized estimating equation approach provides in model assumption flexibility of correlation structures for repeated measurements to produce satisfactory concordance correlation coefficient estimation results.

Estimation of Genetic Parameters of Type Traits in First Parity Cows of the Autochthonous Cika Cattle in Slovenia

The aim of this study was to estimate genetic parameters of 26 individual and four composite type traits in first parity Cika cows. An analysis of variance was performed with the generalized linear model procedure of the SAS/STAT statistical package, where the fixed effects of year of recording, cow's age at recording and days after calving as a linear regression were included in the model. The variance components for the direct additive genetic effect and the herd effect in all type traits were estimated using the REML method in the VCE-6 software package. The estimated heritabilities ranged from 0.42 to 0.67 for the measured body frame traits, from 0.36 to 0.80 for the scored autochthonous traits, from 0.11 to 0.61 for the scored body frame traits, and from 0.20 to 0.47 for the scored udder traits. The estimated heritabilities for the composite traits called "autochthonous characteristics", "muscularity", "body frame" and "udder" were 0.55, 0.19, 0.19, and 0.26, respectively. The estimated genetic correlations among the measured body frame traits were positive and high, while the majority of them among the scored body frame traits were low to moderate. The estimated proportions of variance explained by the herd effect for the composite traits "autochthonous characteristics," "muscularity," "body frame" and "udder" were 0.09, 0.28, 0.14, and 0.10, respectively. The estimated heritabilities for the type traits of first parity Cika cows were similar to those reported for other breeds where breeding values have been routinely predicted for a long time. All estimated genetic parameters are already used for breeding value prediction in the Cika cattle population.

Using population crossover trials to improve the decision process regarding treatment individualization in N-of-1 trials

Healthcare researchers are showing renewed interest in the utilization of N-of-1 clinical trials for the individualization of pharmacological treatments. Here, we propose a frequentist approach to conducting treatment individualization in N-of-1 trials that we call "partial empirical Bayes." We infer the most beneficial treatment for the patient from combining the information provided by a previously conducted population crossover trial with individual patient data. We propose a method for estimating an optimal number of treatment cycles and investigate the statistical conditions under which N-of-1 trials are more beneficial than traditional clinical approaches. We represent the patient population with a random-coefficients linear model and calculate estimators of posttreatment individual disease severities. We show the estimators' consistency under the most common N-of-1 designs and examine their prediction errors and performance with small numbers of patient's responses. We demonstrate by simulating new patients that our approach is equivalent or superior to both the common clinical practice of recommending the on-average best treatment for all patients and the common individualization method that simply compares average responses to the tested treatments. We conclude that some situations exist in which individualization with N-of-1 trials is highly beneficial while other situations exist in which individualization may be unfruitful.

Genetic selection indices for growth traits in Blanco-Orejinegro cattle

Selection indices are used in genetic improvement programs, with the purpose of selectins simultaneous for several economically important traits. The objective of this study was to construct equations for selection indices in the Blanco-Orejinegro (BON) breed and to determine the index that would generate the greatest genetic progress. The information used included birth weight (BW), body weights adjusted to 120, 240, 480, and 720 days old (W120, W240, W240, 480 and W720, respectively), age at first calving (AFC) and interval between first and second calving (IBC) estimated breeding values. Two Smith and Hazel indices were calculated using variances (I1) and literature (I2), with a part two indices designed using information from experts and breeders (I3 and I4). All the indices favored the reduction of weight at birth. The I1 focused mainly on W120 and I2, I3 and I4 focused on 720. In general, the estimated indices obtained similar reliability and expected genetic differences I1 generated a decrease in direct BW. I2 generated the largest increases in BW and AFC. I3 and I4 generated positive changes in growth and reproductive traits, with I3 generating the greatest genetic gains in the population, especially for W240.

GPS Coordinates for Modelling Correlated Herd Effects in Genomic Prediction Models Applied to Hanwoo Beef Cattle

Simple Summary

It is widely known that the environment influences phenotypic expression and that its effects must be accounted for in genetic evaluation programs. The most used method to account for environmental effects is to add herd and the contemporary group to the model. Although generally informative, the herd effect treats different farms as independent units. However, if two farms are located physically close to each other, they potentially share correlated environmental factors. We introduce a method to model herd effects using physical distances between farms based on GPS coordinates as a proxy for the correlation matrix of these effects, aiming to account for similarities and differences between farms due to environmental factors. A population of beef cattle was used to evaluate the impact on the variance components and on the genomic prediction, of modelling herd effects as correlated, in comparison to assuming the farms as completely independent units. The main result was an increase in the reliabilities of the predicted genomic breeding values compared to reliabilities obtained with traditional models, a finding of practical relevance for genetic evaluation programs.

Abstract

It is widely known that the environment influences phenotypic expression and that its effects must be accounted for in genetic evaluation programs. The most used method to account for environmental effects is to add herd and contemporary group to the model. Although generally informative, the herd effect treats different farms as independent units. However, if two farms are located physically close to each other, they potentially share correlated environmental factors. We introduce a method to model herd effects that uses the physical distances between farms based on the Global Positioning System (GPS) coordinates as a proxy for the correlation matrix of these effects that aims to account for similarities and differences between farms due to environmental factors. A population of Hanwoo Korean cattle was used to evaluate the impact of modelling herd effects as correlated, in comparison to assuming the farms as completely independent units, on the variance components and genomic prediction. The main result was an increase in the reliabilities of the predicted genomic breeding values compared to reliabilities obtained with traditional models (across four traits evaluated, reliabilities of prediction presented increases that ranged from 0.05 ± 0.01 to 0.33 ± 0.03), suggesting that these models may overestimate heritabilities. Although little to no significant gain was obtained in phenotypic prediction, the increased reliability of the predicted genomic breeding values is of practical relevance for genetic evaluation programs.

Quantifying the contribution of dominance deviation effects to complex trait variation in biobank-scale data

The proportion of variation in complex traits that can be attributed to non-additive genetic effects has been a topic of intense debate. The availability of biobank-scale datasets of genotype and trait data from unrelated individuals opens up the possibility of obtaining precise estimates of the contribution of non-additive genetic effects. We present an efficient method to estimate the variation in a complex trait that can be attributed to additive (additive heritability) and dominance deviation (dominance heritability) effects across all genotyped SNPs in a large collection of unrelated individuals. Over a wide range of genetic architectures, our method yields unbiased estimates of additive and dominance heritability. We applied our method, in turn, to array genotypes as well as imputed genotypes (at common SNPs with minor allele frequency [MAF] > 1%) and 50 quantitative traits measured in 291,273 unrelated white British individuals in the UK Biobank. Averaged across these 50 traits, we find that additive heritability on array SNPs is 21.86% while dominance heritability is 0.13% (about 0.48% of the additive heritability) with qualitatively similar results for imputed genotypes. We find no statistically significant evidence for dominance heritability (p<0.05/50 accounting for the number of traits tested) and estimate that dominance heritability is unlikely to exceed 1% for the traits analyzed. Our analyses indicate a limited contribution of dominance heritability to complex trait variation.

Inference for high-dimensional linear mixed-effects models: A quasi-likelihood approach

Linear mixed-effects models are widely used in analyzing clustered or repeated measures data. We propose a quasi-likelihood approach for estimation and inference of the unknown parameters in linear mixed-effects models with high-dimensional fixed effects. The proposed method is applicable to general settings where the dimension of the random effects and the cluster sizes are possibly large. Regarding the fixed effects, we provide rate optimal estimators and valid inference procedures that do not rely on the structural information of the variance components. We also study the estimation of variance components with high-dimensional fixed effects in general settings. The algorithms are easy to implement and computationally fast. The proposed methods are assessed in various simulation settings and are applied to a real study regarding the associations between body mass index and genetic polymorphic markers in a heterogeneous stock mice population.