Impact of inclusion non-additive effects on genome-wide association and variance's components in Scottish black sheep.

CONTEXT: It's well-documented that most economic traits have a complex genetic structure that is controlled by additive and non-additive gene actions. Hence, knowledge of the underlying genetic architecture of such complex traits could aid in understanding how these traits respond to the selection in breeding and mating programs. Computing and having estimates of the non-additive effect for economic traits in sheep using genome-wide information can be important because; non-additive genes play an important role in the prediction accuracy of genomic breeding values and the genetic response to the selection. AIM: This study aimed to assess the impact of non-additive effects (dominance and epistasis) on the estimation of genetic parameters for body weight traits in sheep. METHODS: This study used phenotypic and genotypic belonging to 752 Scottish Blackface lambs. Three live weight traits considered in this study were included in body weight at 16, 20, and 24 weeks). Three genetic models including additive (AM), additive + dominance (ADM), and additive + dominance + epistasis (ADEM), were used. KEY RESULTS: The narrow sense heritability for weight at 16 weeks of age (BW16) were 0.39, 0.35, and 0.23, for 20 weeks of age (BW20) were 0.55, 0.54, and 0.42, and finally for 24 weeks of age (BW24) were 0.16, 0.12, and 0.02, using the AM, ADM, and ADEM models, respectively. The additive genetic model significantly outperformed the non-additive genetic model (p < 0.01). The dominance variance of the BW16, BW20, and BW24 accounted for 38, 6, and 30% of the total phenotypic, respectively. Moreover, the epistatic variance accounted for 39, 0.39, and 47% of the total phenotypic variances of these traits, respectively. In addition, our results indicated that the most important SNPs for live weight traits are on chromosomes 3 (three SNPS including s12606.1, OAR3_221188082.1, and OAR3_4106875.1), 8 (OAR8_16468019.1, OAR8_18067475.1, and OAR8_18043643.1), and 19 (OAR19_18010247.1), according to the genome-wide association analysis using additive and non-additive genetic model. CONCLUSIONS: The results emphasized that the non-additive genetic effects play an important role in controlling body weight variation at the age of 16-24 weeks in Scottish Blackface lambs. IMPLICATIONS: It is expected that using a high-density SNP panel and the joint modeling of both additive and non-additive effects can lead to better estimation and prediction of genetic parameters.

Predictive assessment of single-step BLUP with linear and non-linear similarity RKHS kernels: A case study in chickens.

Single-step GBLUP (ssGBLUP) to obtain genomic prediction was proposed in 2009. Many studies have investigated ssGBLUP in genomic selection in animals and plants using a standard linear kernel (similarity matrix) called genomic relationship matrix (G). More general kernels should allow capturing non-additive effects as well, whereas GBLUP is based on additive gene action. In this study, we generalized ssBLUP to accommodate two non-linear kernels, the averaged Gaussian kernel (AK) and the recently developed arc-cosine deep kernel (DK). We evaluated the methodology using body weight (BW) and hen-housing production (HHP) traits, recorded on a sample of phenotyped and genotyped commercial broiler chickens. There were, thus, different ssGBLUP models corresponding to G, AK and DK. We used random replication of training (TRN) and testing (TST) layouts at different genotyping rates (20%, 40%, 60% and 80% of all birds) in three selective genotyping scenarios. The selections were genotyping the youngest individuals in the pedigree (YS), random genotyping (RS) and genotyping based on parent average (PA). Predictive abilities were measured using rank correlations between the observed and the predictive phenotypic values in TST for each random partition. Prediction accuracy was influenced by the type of kernel when a large proportion of birds was genotyped. An advantage of non-linear kernels (AK and DK) was more apparent when 60 and 80% of birds had been genotyped. For BW, the lowest rank correlations were obtained with G (0.093 ± 0.015 using RS by 20% genotyped individuals) and the highest values with DK (0.320 ± 0.016 in the PA setting with 80% genotyped individuals). For HHP, the lowest and highest rank correlations were obtained by AK with 20% and 80% genotyped individuals, 0.071 ± 0.016 (in RS) and 0.23 ± 0.016 (in PA) respectively. Our results indicated that AK and DK are more effective than G when a large proportion of the target population is genotyped. Our expectation is that ssGBLUP with AK or DK models can perform even better than G when non-additive genetic effects influence the underlying variability of complex traits.

A selection signatures study among Middle Eastern and European sheep breeds.

Selection, both natural and artificial, leaves patterns on the genome during domestication of animals and leads to changes in allele frequencies among populations. Detecting genomic regions influenced by selection in livestock may assist in understanding the processes involved in genome evolution and discovering genomic regions related to traits of economic and ecological interests. In the current study, genetic diversity analyses were conducted on 34,206 quality-filtered SNP positions from 450 individuals in 15 sheep breeds, including six indigenous breeds from the Middle East, namely Iranian Balouchi, Afshari, Moghani, Qezel, Karakas and Norduz, and nine breeds from Europe, namely East Friesian Sheep, Ile de France, Mourerous, Romane, Swiss Mirror, Spaelsau, Suffolk, Comisana and Engadine Red Sheep. The SNP genotype data generated by the Illumina OvineSNP50 Genotyping BeadChip array were used in this analysis. We applied two complementary statistical analyses, FST (fixation index) and xp-EHH (cross-population extended haplotype homozygosity), to detect selection signatures in Middle Eastern and European sheep populations. FST and xp-EHH detected 629 and 256 genes indicating signatures of selection, respectively. Genomic regions identified using FST and xp-EHH contained the CIDEA, HHATL, MGST1, FADS1, RTL1 and DGKG genes, which were reported earlier to influence a number of economic traits. Both FST and xp-EHH approaches identified 60 shared genes as the signatures of selection, including four candidate genes (NT5E, ADA2, C8A and C8B) that were enriched for two significant Gene Ontology (GO) terms associated with the adenosine metabolic procedure. Knowledge about the candidate genomic regions under selective pressure in sheep breeds may facilitate identification of the underlying genes and enhance our understanding on these genes role in local adaptation.

Integration of single nucleotide variants and whole-genome DNA methylation profiles for classification of rheumatoid arthritis cases from controls.

This study evaluated the use of multiomics data for classification accuracy of rheumatoid arthritis (RA). Three approaches were used and compared in terms of prediction accuracy: (1) whole-genome prediction (WGP) using SNP marker information only, (2) whole-methylome prediction (WMP) using methylation profiles only, and (3) whole-genome/methylome prediction (WGMP) with combining both omics layers. The number of SNP and of methylation sites varied in each scenario, with either 1, 10, or 50% of these preselected based on four approaches: randomly, evenly spaced, lowest p value (genome-wide association or epigenome-wide association study), and estimated effect size using a Bayesian ridge regression (BRR) model. To remove effects of high levels of pairwise linkage disequilibrium (LD), SNPs were also preselected with an LD-pruning method. Five Bayesian regression models were studied for classification, including BRR, Bayes-A, Bayes-B, Bayes-C, and the Bayesian LASSO. Adjusting methylation profiles for cellular heterogeneity within whole blood samples had a detrimental effect on the classification ability of the models. Overall, WGMP using Bayes-B model has the best performance. In particular, selecting SNPs based on LD-pruning with 1% of the methylation sites selected based on BRR included in the model, and fitting the most significant SNP as a fixed effect was the best method for predicting disease risk with a classification accuracy of 0.975. Our results showed that multiomics data can be used to effectively predict the risk of RA and identify cases in early stages to prevent or alter disease progression via appropriate interventions.

Performance of pedigree and various forms of marker-derived relationship coefficients in genomic prediction and their correlations.

In recent years, with development and validation of different genotyping panels, several methods have been proposed to build efficient similarity matrices among individuals to be used for genomic selection. Consequently, the estimated genetic parameters from such information may deviate from their counterpart using traditional family information. In this study, we used a pedigree-based numerator relationship matrix (A) and three types of marker-based relationship matrices ( G ) including two identical by descent, that is G K and G M and one identical by state, G V as well as four Gaussian kernel ( GK ) similarity kernels with different smoothing parameters to predict yet to be observed phenotypes. Also, we used different kinship matrices that are a linear combination of marker-derived IBD or IBS matrices with A, constructed as K = λ G + 1 - λ A , where the weight ( λ ) assigned to each source of information varied over a grid of values. A Bayesian multiple-trait Gaussian model was fitted to estimate the genetic parameters and compare the prediction accuracy in terms of predictive correlation, mean square error and unbiasedness. Results show that the estimated genetic parameters (heritability and correlations) are affected by the source of the information used to create kinship or the weight placed on the sources of genomic and pedigree information. The superiority of GK-based model depends on the smoothing parameters (θ) so that with an optimum θ value, the GK-based model statistically yielded better performance (higher predictive correlation, lowest MSE and unbiased estimates) and more stable correlations and heritability than the model with IBD, IBS or A kinship matrices or any of the linear combinations.

Quantifying genomic connectedness and prediction accuracy from additive and non-additive gene actions.

BACKGROUND: Genetic connectedness is classically used as an indication of the risk associated with breeding value comparisons across management units because genetic evaluations based on best linear unbiased prediction rely for their success on sufficient linkage among different units. In the whole-genome prediction era, the concept of genetic connectedness can be extended to measure a connectedness level between reference and validation sets. However, little is known regarding (1) the impact of non-additive gene action on genomic connectedness measures and (2) the relationship between the estimated level of connectedness and prediction accuracy in the presence of non-additive genetic variation. RESULTS: We evaluated the extent to which non-additive kernel relationship matrices increase measures of connectedness and investigated its relationship with prediction accuracy in the cross-validation framework using best linear unbiased prediction and coefficients of determination. Simulated data assuming additive, dominance, and epistatic gene action scenarios and real swine data were analyzed. We found that the joint use of additive and non-additive genomic kernel relationship matrices or non-parametric relationship matrices led to increased capturing of connectedness, up to 25%, and improved prediction accuracies compared to those of baseline additive relationship counterparts in the presence of non-additive gene action. CONCLUSIONS: Our findings showed that connectedness metrics can be extended to incorporate non-additive genetic variation of complex traits. Use of kernel relationship matrices designed to capture non-additive gene action increased measures of connectedness and improved whole-genome prediction accuracy, further broadening the scope of genomic connectedness studies.

A predictive assessment of genetic correlations between traits in chickens using markers.

BACKGROUND: Genomic selection has been successfully implemented in plant and animal breeding programs to shorten generation intervals and accelerate genetic progress per unit of time. In practice, genomic selection can be used to improve several correlated traits simultaneously via multiple-trait prediction, which exploits correlations between traits. However, few studies have explored multiple-trait genomic selection. Our aim was to infer genetic correlations between three traits measured in broiler chickens by exploring kinship matrices based on a linear combination of measures of pedigree and marker-based relatedness. A predictive assessment was used to gauge genetic correlations. METHODS: A multivariate genomic best linear unbiased prediction model was designed to combine information from pedigree and genome-wide markers in order to assess genetic correlations between three complex traits in chickens, i.e. body weight at 35 days of age (BW), ultrasound area of breast meat (BM) and hen-house egg production (HHP). A dataset with 1351 birds that were genotyped with the 600 K Affymetrix platform was used. A kinship kernel (K) was constructed as K = λ G + (1 - λ)A, where A is the numerator relationship matrix, measuring pedigree-based relatedness, and G is a genomic relationship matrix. The weight (λ) assigned to each source of information varied over the grid λ = (0, 0.2, 0.4, 0.6, 0.8, 1). Maximum likelihood estimates of heritability and genetic correlations were obtained at each λ, and the "optimum" λ was determined using cross-validation. RESULTS: Estimates of genetic correlations were affected by the weight placed on the source of information used to build K. For example, the genetic correlation between BW-HHP and BM-HHP changed markedly when λ varied from 0 (only A used for measuring relatedness) to 1 (only genomic information used). As λ increased, predictive correlations (correlation between observed phenotypes and predicted breeding values) increased and mean-squared predictive error decreased. However, the improvement in predictive ability was not monotonic, with an optimum found at some 0 < λ < 1, i.e., when both sources of information were used together. CONCLUSIONS: Our findings indicate that multiple-trait prediction may benefit from combining pedigree and marker information. Also, it appeared that expected correlated responses to selection computed from standard theory may differ from realized responses. The predictive assessment provided a metric for performance evaluation as well as a means for expressing uncertainty of outcomes of multiple-trait selection.

Neurofilament light chain concentration does not correlate with disease status in Labrador Retrievers affected with idiopathic laryngeal paralysis.

OBJECTIVE: The aim of this study was to investigate whether plasma neurofilament light chain (pNfL) concentration was altered in Labrador Retrievers with idiopathic laryngeal paralysis (ILP) compared to a control population. A secondary aim was to investigate relationships between age, height, weight, and body mass index in the populations studied. ANIMALS: 123 dogs: 62 purebred Labrador Retrievers with ILP (ILP Cases) and 61 age-matched healthy medium- to large-breed dogs (Controls). METHODS: Dogs, recruited from August 1, 2016, to March 1, 2022, were categorized as case or control based on a combination of physical exam, neurologic exam, and history. Blood plasma was collected, and pNfL concentration was measured. pNfL concentrations were compared between ILP Cases and Controls. Covariables including age, height, and weight were collected. Relationships between pNfL and covariables were analyzed within and between groups. In dogs where 2 plasma samples were available from differing time points, pNfL concentrations were measured to evaluate alterations over time. RESULTS: No significant difference in pNfL concentration was found between ILP Cases and Control (P = .36). pNfL concentrations had moderate negative correlations with weight and height in the Control group; other variables did not correlate with pNfL concentrations in ILP Case or Control groups. pNfL concentrations do not correlate with ILP disease status or duration in Labrador Retrievers. CLINICAL RELEVANCE: There is no evidence that pNfL levels are altered due to ILP disease duration or progression when compared with healthy controls. When evaluating pNfL concentrations in the dog, weight and height should be considered.

Pathologic classification of a late-onset peripheral neuropathy in a spontaneous Labrador retriever dog model.

Late-onset peripheral neuropathy (LPN) is a heritable canine neuropathy commonly found in Labrador retrievers and is characterized by laryngeal paralysis and pelvic limb paresis. Our objective was to establish canine LPN as a model for human hereditary peripheral neuropathy by classifying it as either an axonopathy or myelinopathy and evaluating length-dependent degeneration. We conducted a motor nerve conduction study of the sciatic and ulnar nerves, electromyography (EMG) of appendicular and epaxial musculature, and histologic analysis of sciatic and recurrent laryngeal nerves in LPN-affected and control dogs. LPN-affected dogs exhibited significant decreases in compound muscle action potential (CMAP) amplitude, CMAP area, and pelvic limb latencies. However, no differences were found in motor nerve conduction velocity, residual latencies, or CMAP duration. Distal limb musculature showed greater EMG changes in LPN-affected dogs. Histologically, LPN-affected dogs exhibited a reduction in the number of large-diameter axons, especially in distal nerve regions. In conclusion, LPN in Labrador retrievers is a common, spontaneous, length-dependent peripheral axonopathy that is a novel animal model of age-related peripheral neuropathy that could be used for fundamental research and clinical trials.

Genetic architecture and polygenic risk score prediction of degenerative suspensory ligament desmitis (DSLD) in the Peruvian Horse.

Introduction: Spontaneous rupture of tendons and ligaments is common in several species including humans. In horses, degenerative suspensory ligament desmitis (DSLD) is an important acquired idiopathic disease of a major energy-storing tendon-like structure. DSLD risk is increased in several breeds, including the Peruvian Horse. Affected horses have often been used for breeding before the disease is apparent. Breed predisposition suggests a substantial genetic contribution, but heritability and genetic architecture of DSLD have not been determined. Methods: To identify genomic regions associated with DSLD, we recruited a reference population of 183 Peruvian Horses, phenotyped as DSLD cases or controls, and undertook a genome-wide association study (GWAS), a regional window variance analysis using local genomic partitioning, a signatures of selection (SOS) analysis, and polygenic risk score (PRS) prediction of DSLD risk. We also estimated trait heritability from pedigrees. Results: Heritability was estimated in a population of 1,927 Peruvian horses at 0.22 ± 0.08. After establishing a permutation-based threshold for genome-wide significance, 151 DSLD risk single nucleotide polymorphisms (SNPs) were identified by GWAS. Multiple regions of enriched local heritability were identified across the genome, with strong enrichment signals on chromosomes 1, 2, 6, 10, 13, 16, 18, 22, and the X chromosome. With SOS analysis, there were 66 genes with a selection signature in DSLD cases that was not present in the control group that included the TGFB3 gene. Pathways enriched in DSLD cases included proteoglycan metabolism, extracellular matrix homeostasis, and signal transduction pathways that included the hedgehog signaling pathway. The best PRS predictive performance was obtained when we fitted 1% of top SNPs using a Bayesian Ridge Regression model which achieved the highest mean of R2 on both the probit and logit liability scales, indicating a strong predictive performance. Discussion: We conclude that within-breed GWAS of DSLD in the Peruvian Horse has further confirmed that moderate heritability and a polygenic architecture underlies the trait and identified multiple DSLD SNP associations in novel tendinopathy candidate genes influencing disease risk. Pathways enriched with DSLD risk variants include ones that influence glycosaminoglycan metabolism, extracellular matrix homeostasis, signal transduction pathways.

Across-breed genetic investigation of canine hip dysplasia, elbow dysplasia, and anterior cruciate ligament rupture using whole-genome sequencing.

Here, we report the use of genome-wide association study (GWAS) for the analysis of canine whole-genome sequencing (WGS) repository data using breed phenotypes. Single-nucleotide polymorphisms (SNPs) were called from WGS data from 648 dogs that included 119 breeds from the Dog10K Genomes Project. Next, we assigned breed phenotypes for hip dysplasia (Orthopedic Foundation for Animals (OFA) HD, n = 230 dogs from 27 breeds; hospital HD, n = 279 dogs from 38 breeds), elbow dysplasia (ED, n = 230 dogs from 27 breeds), and anterior cruciate ligament rupture (ACL rupture, n = 279 dogs from 38 breeds), the three most important canine spontaneous complex orthopedic diseases. Substantial morbidity is common with these diseases. Previous within- and between-breed GWAS for HD, ED, and ACL rupture using array SNPs have identified disease-associated loci. Individual disease phenotypes are lacking in repository data. There is a critical knowledge gap regarding the optimal approach to undertake categorical GWAS without individual phenotypes. We considered four GWAS approaches: a classical linear mixed model, a haplotype-based model, a binary case-control model, and a weighted least squares model using SNP average allelic frequency. We found that categorical GWAS was able to validate HD candidate loci. Additionally, we discovered novel candidate loci and genes for all three diseases, including FBX025, IL1A, IL1B, COL27A1, SPRED2 (HD), UGDH, FAF1 (ED), TGIF2 (ED & ACL rupture), and IL22, IL26, CSMD1, LDHA, and TNS1 (ACL rupture). Therefore, categorical GWAS of ancestral dog populations may contribute to the understanding of any disease for which breed epidemiological risk data are available, including diseases for which GWAS has not been performed and candidate loci remain elusive.

Selection signature analyses and genome-wide association reveal genomic hotspot regions that reflect differences between breeds of horse with contrasting risk of degenerative suspensory ligament desmitis.

Degenerative suspensory ligament desmitis is a progressive idiopathic condition that leads to scarring and rupture of suspensory ligament fibers in multiple limbs in horses. The prevalence of degenerative suspensory ligament desmitis is breed related. Risk is high in the Peruvian Horse, whereas pony and draft breeds have low breed risk. Degenerative suspensory ligament desmitis occurs in families of Peruvian Horses, but its genetic architecture has not been definitively determined. We investigated contrasts between breeds with differing risk of degenerative suspensory ligament desmitis and identified associated risk variants and candidate genes. We analyzed 670k single nucleotide polymorphisms from 10 breeds, each of which was assigned one of the four breed degenerative suspensory ligament desmitis risk categories: control (Belgian, Icelandic Horse, Shetland Pony, and Welsh Pony), low risk (Lusitano, Arabian), medium risk (Standardbred, Thoroughbred, Quarter Horse), and high risk (Peruvian Horse). Single nucleotide polymorphisms were used for genome-wide association and selection signature analysis using breed-assigned risk levels. We found that the Peruvian Horse is a population with low effective population size and our breed contrasts suggest that degenerative suspensory ligament desmitis is a polygenic disease. Variant frequency exhibited signatures of positive selection across degenerative suspensory ligament desmitis breed risk groups on chromosomes 7, 18, and 23. Our results suggest degenerative suspensory ligament desmitis breed risk is associated with disturbances to suspensory ligament homeostasis where matrix responses to mechanical loading are perturbed through disturbances to aging in tendon (PIN1), mechanotransduction (KANK1, KANK2, JUNB, SEMA7A), collagen synthesis (COL4A1, COL5A2, COL5A3, COL6A5), matrix responses to hypoxia (PRDX2), lipid metabolism (LDLR, VLDLR), and BMP signaling (GREM2). Our results do not suggest that suspensory ligament proteoglycan turnover is a primary factor in disease pathogenesis.

Genic and non-genic SNP contributions to additive and dominance genetic effects in purebred and crossbred pig traits.

The present research has estimated the additive and dominance genetic variances of genic and intergenic segments for average daily gain (ADG), backfat thickness (BFT) and pH of the semimembranosus dorsi muscle (PHS). Further, the predictive performance using additive and additive dominance models in a purebred Piétrain (PB) and a crossbred (Piétrain × Large White, CB) pig population was assessed. All genomic regions contributed equally to the additive and dominance genetic variations and lead to the same predictive ability that did not improve with the inclusion of dominance genetic effect and inbreeding in the models. Using all SNPs available, additive genotypic correlations between PB and CB performances for the three traits were high and positive (> 0.83) and dominance genotypic correlation was very inaccurate. Estimates of dominance genotypic correlations between all pairs of traits in both populations were imprecise but positive for ADG-BFT in CB and BFT-PHS in PB and CB with a high probability (> 0.98). Additive and dominance genotypic correlations between BFT and PHS were of different sign in both populations, which could indicate that genes contributing to the additive genetic progress in both traits would have an antagonistic effect when used for exploiting dominance effects in planned matings.

Multi-omic data integration for the study of production, carcass, and meat quality traits in Nellore cattle.

Data integration using hierarchical analysis based on the central dogma or common pathway enrichment analysis may not reveal non-obvious relationships among omic data. Here, we applied factor analysis (FA) and Bayesian network (BN) modeling to integrate different omic data and complex traits by latent variables (production, carcass, and meat quality traits). A total of 14 latent variables were identified: five for phenotype, three for miRNA, four for protein, and two for mRNA data. Pearson correlation coefficients showed negative correlations between latent variables miRNA 1 (mirna1) and miRNA 2 (mirna2) (-0.47), ribeye area (REA) and protein 4 (prot4) (-0.33), REA and protein 2 (prot2) (-0.3), carcass and prot4 (-0.31), carcass and prot2 (-0.28), and backfat thickness (BFT) and miRNA 3 (mirna3) (-0.25). Positive correlations were observed among the four protein factors (0.45-0.83): between meat quality and fat content (0.71), fat content and carcass (0.74), fat content and REA (0.76), and REA and carcass (0.99). BN presented arcs from the carcass, meat quality, prot2, and prot4 latent variables to REA; from meat quality, REA, mirna2, and gene expression mRNA1 to fat content; from protein 1 (prot1) and mirna2 to protein 5 (prot5); and from prot5 and carcass to prot2. The relations of protein latent variables suggest new hypotheses about the impact of these proteins on REA. The network also showed relationships among miRNAs and nebulin proteins. REA seems to be the central node in the network, influencing carcass, prot2, prot4, mRNA1, and meat quality, suggesting that REA is a good indicator of meat quality. The connection among miRNA latent variables, BFT, and fat content relates to the influence of miRNAs on lipid metabolism. The relationship between mirna1 and prot5 composed of isoforms of nebulin needs further investigation. The FA identified latent variables, decreasing the dimensionality and complexity of the data. The BN was capable of generating interrelationships among latent variables from different types of data, allowing the integration of omics and complex traits and identifying conditional independencies. Our framework based on FA and BN is capable of generating new hypotheses for molecular research, by integrating different types of data and exploring non-obvious relationships.

Modeling multiple phenotypes in wheat using data-driven genomic exploratory factor analysis and Bayesian network learning.

Inferring trait networks from a large volume of genetically correlated diverse phenotypes such as yield, architecture, and disease resistance can provide information on the manner in which complex phenotypes are interrelated. However, studies on statistical methods tailored to multidimensional phenotypes are limited, whereas numerous methods are available for evaluating the massive number of genetic markers. Factor analysis operates at the level of latent variables predicted to generate observed responses. The objectives of this study were to illustrate the manner in which data-driven exploratory factor analysis can map observed phenotypes into a smaller number of latent variables and infer a genomic latent factor network using 45 agro-morphological, disease, and grain mineral phenotypes measured in synthetic hexaploid wheat lines (Triticum aestivum L.). In total, eight latent factors including grain yield, architecture, flag leaf-related traits, grain minerals, yellow rust, two types of stem rust, and leaf rust were identified as common sources of the observed phenotypes. The genetic component of the factor scores for each latent variable was fed into a Bayesian network to obtain a trait structure reflecting the genetic interdependency among traits. Three directed paths were consistently identified by two Bayesian network algorithms. Flag leaf-related traits influenced leaf rust, and yellow rust and stem rust influenced grain yield. Additional paths that were identified included flag leaf-related traits to minerals and minerals to architecture. This study shows that data-driven exploratory factor analysis can reveal smaller dimensional common latent phenotypes that are likely to give rise to numerous observed field phenotypes without relying on prior biological knowledge. The inferred genomic latent factor structure from the Bayesian network provides insights for plant breeding to simultaneously improve multiple traits, as an intervention on one trait will affect the values of focal phenotypes in an interrelated complex trait system.

Heritability and genetic variance estimation of Osteosarcoma (OSA) in Irish Wolfhound, using deep pedigree information.

BACKGROUND: Osteosarcoma (OSA) is a devastating disease that is common in the Irish Wolfhound breed. The aim of this study was to use a pedigree-based approach to determine the heritability of OSA in the Irish Wolfhound using data from a large publically available database. RESULTS: The pedigree used for this study included 5110 pure-bred Irish Wolfhounds, including 332 dogs diagnosed with OSA and 360 control dogs; dogs were considered controls if they lived over 10 years of age and were not reported to have developed OSA. The estimated heritability of OSA in the Irish Wolfhound was 0.65. CONCLUSION: The results of this study indicate that OSA in the Irish Wolfhound is highly heritable, and support the need for future research investigating associated genetic mutations.

Osteosarcoma is a devastating condition that is prevalent in the Irish Wolfhound breed. In this study, our aim was to estimate heritability of osteosarcoma in the Irish Wolfhound breed. We undertook a pedigree-based analysis to estimate heritability of osteosarcoma in the Irish Wolfhound. The pedigree used included 5110 pure-bred Irish Wolfhounds, including 332 dogs diagnosed with osteosarcoma and 360 control dogs. We considered dogs to be controls if they were over 10 years of age and were not reported to have developed osteosarcoma. This study found the heritability estimate of osteosarcoma in the Irish Wolfhound to be 0.65. This score means that osteosarcoma in this breed is: 1) highly heritable and 2) a complex trait, which means that both environmental and genetic factors influence disease risk. Overall, our results provide support for further investigation into the genetic variants involved in the development of osteosarcoma in Irish Wolfhounds.

Biologically Enhanced Genome-Wide Association Study Provides Further Evidence for Candidate Loci and Discovers Novel Loci That Influence Risk of Anterior Cruciate Ligament Rupture in a Dog Model.

Anterior cruciate ligament (ACL) rupture is a common condition that disproportionately affects young people, 50% of whom will develop knee osteoarthritis (OA) within 10 years of rupture. ACL rupture exhibits both hereditary and environmental risk factors, but the genetic basis of the disease remains unexplained. Spontaneous ACL rupture in the dog has a similar disease presentation and progression, making it a valuable genomic model for ACL rupture. We leveraged the dog model with Bayesian mixture model (BMM) analysis (BayesRC) to identify novel and relevant genetic variants associated with ACL rupture. We performed RNA sequencing of ACL and synovial tissue and assigned single nucleotide polymorphisms (SNPs) within differentially expressed genes to biological prior classes. SNPs with the largest effects were on chromosomes 3, 5, 7, 9, and 24. Selection signature analysis identified several regions under selection in ACL rupture cases compared to controls. These selection signatures overlapped with genome-wide associations with ACL rupture as well as morphological traits. Notable findings include differentially expressed ACSF3 with MC1R (coat color) and an association on chromosome 7 that overlaps the boundaries of SMAD2 (weight and body size). Smaller effect associations were within or near genes associated with regulation of the actin cytoskeleton and the extracellular matrix, including several collagen genes. The results of the current analysis are consistent with previous work published by our laboratory and others, and also highlight new genes in biological pathways that have not previously been associated with ACL rupture. The genetic associations identified in this study mirror those found in human beings, which lays the groundwork for development of disease-modifying therapies for both species.

Prediction of biological age and evaluation of genome-wide dynamic methylomic changes throughout human aging.

The use of DNA methylation signatures to predict chronological age and aging rate is of interest in many fields, including disease prevention and treatment, forensics, and anti-aging medicine. Although a large number of methylation markers are significantly associated with age, most age-prediction methods use a few markers selected based on either previously published studies or datasets containing methylation information. Here, we implemented reproducing kernel Hilbert spaces (RKHS) regression and a ridge regression model in a Bayesian framework that utilized phenotypic and methylation profiles simultaneously to predict chronological age. We used over 450,000 CpG sites from the whole blood of a large cohort of 4409 human individuals with a range of 10-101 years of age. Models were fitted using adjusted and un-adjusted methylation measurements for cell heterogeneity. Un-adjusted methylation scores delivered a significantly higher prediction accuracy than adjusted methylation data, with a correlation between age and predicted age of 0.98 and a root mean square error (RMSE) of 3.54 years in un-adjusted data, and 0.90 (correlation) and 7.16 (RMSE) years in adjusted data. Reducing the number of predictors (CpG sites) through subset selection improved predictive power with a correlation of 0.98 and an RMSE of 2.98 years in the RKHS model. We found distinct global methylation patterns, with a significant increase in the proportion of methylated cytosines in CpG islands and a decreased proportion in other CpG types, including CpG shore, shelf, and open sea (P < 5e-06). Epigenetic drift seemed to be a widespread phenomenon as more than 97% of the age-associated methylation sites had heteroscedasticity. Apparent methylomic aging rate (AMAR) had a sex-specific pattern, with an increase in AMAR in females with age related to males.

Neurofilament light plasma concentration positively associates with age and negatively associates with weight and height in the dog.

Plasma neurofilament light chain (pNfL) concentration is a biomarker for neuroaxonal injury and degeneration and can be used to monitor response to treatment. Spontaneous canine neurodegenerative diseases are a valuable comparative resource for understanding similar human conditions and as large animal treatment models. The features of pNfL concentration in healthy dogs is not well established. We present data reporting basic pNfL concentration trends in the Labrador Retriever breed. Fifty-five Labrador Retrievers were enrolled. pNfL concentration was measured and correlated to age, sex, neuter status, height, weight, body mass index, and coat color. We found increased pNfL with age (P < 0.0001), shorter stature (P = 0.009) and decreased body weight (P < 0.001). These are similar to findings reported in humans. pNfL concentration did not correlate with sex, BMI or coat color. This data further supports findings that pNfL increase with age in a canine population but highlights a need to consider weight and height when determining normal pNfL concentration in canine populations.

Multi-trait random regression models increase genomic prediction accuracy for a temporal physiological trait derived from high-throughput phenotyping.

Random regression models (RRM) are used extensively for genomic inference and prediction of time-valued traits in animal breeding, but only recently have been used in plant systems. High-throughput phenotyping (HTP) platforms provide a powerful means to collect high-dimensional phenotypes throughout the growing season for large populations. However, to date, selection of an appropriate statistical genomic framework to integrate multiple temporal traits for genomic prediction in plants remains unexplored. Here, we demonstrate the utility of a multi-trait RRM (MT-RRM) for genomic prediction of daily water usage (WU) in rice (Oryza sativa) through joint modeling with shoot biomass (projected shoot area, PSA). Three hundred and fifty-seven accessions were phenotyped daily for WU and PSA over 20 days using a greenhouse-based HTP platform. MT-RRMs that modeled additive genetic and permanent environmental effects for both traits using quadratic Legendre polynomials were used to assess genomic correlations between traits and genomic prediction for WU. Predictive abilities of the MT-RRMs were assessed using two cross-validation (CV) scenarios. The first scenario was designed to predict genetic values for WU at all time points for a set of accessions with unobserved WU. The second scenario was designed to forecast future genetic values for WU for a panel of known accessions with records for WU at earlier time periods. In each scenario we evaluated two MT-RRMs in which PSA records were absent or available for time points in the testing population. Weak to strong genomic correlations between WU and PSA were observed across the days of imaging (0.29-0.870.38-0.80). In both CV scenarios, MT-RRMs showed better predictive abilities compared to single-trait RRM, and prediction accuracies were greatly improved when PSA records were available for the testing population. In summary, these frameworks provide an effective approach to predict temporal physiological traits that are difficult or expensive to quantify in large populations.