Current challenges and future of agricultural genomes to phenomes in the USA.

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research.

Purebred-crossbred genetic parameters for reproductive traits in swine.

For swine breeding programs, testing and selection programs are usually within purebred (PB) populations located in nucleus units that are generally managed differently and tend to have a higher health level than the commercial herds in which the crossbred (CB) descendants of these nucleus animals are expected to perform. This approach assumes that PB animals selected in the nucleus herd will have CB progeny that have superior performance at the commercial level. There is clear evidence that this may not be the case for all traits of economic importance and, thus, including data collected at the commercial herd level may increase the accuracy of selection for commercial CB performance at the nucleus level. The goal for this study was to estimate genetic parameters for five maternal reproductive traits between two PB maternal nucleus populations (Landrace and Yorkshire) and their CB offspring: Total Number Born (TNB), Number Born Alive (NBA), Number Born Alive > 1 kg (NBA > 1 kg), Total Number Weaned (TNW), and Litter Weight at Weaning (LWW). Estimates were based on single-step GBLUP by analyzing any two combinations of a PB and the CB population, and by analyzing all three populations jointly. The genomic relationship matrix between the three populations was generated by using within-population allele frequencies for relationships within a population, and across-population allele frequencies for relationships of the CB with the PB animals. Utilization of metafounders for the two PB populations had no effect on parameter estimates, so the two PB populations were assumed to be genetically unrelated. Joint analysis of two (one PB plus CB) vs. three (both PB and CB) populations did not impact estimates of heritability, additive genetic variance, and genetic correlations. Heritabilities were generally similar between the PB and CB populations, except for LWW and TNW, for which PB populations had about four times larger estimates than CB. Purebred-crossbred genetic correlations (rpc) were larger for Landrace than for Yorkshire, except for NBA > 1 kg. These estimates of rpc indicate that there is potential to improve selection of PB animals for CB performance by including CB information for all traits in the Yorkshire population, but that noticeable additional gains may only occur for NBA > 1 kg and TNW in the Landrace population.

Host Genetics of Response to Porcine Reproductive and Respiratory Syndrome in Sows: Reproductive Performance.

Porcine Reproductive and Respiratory Syndrome (PRRS) is historically the most economically important swine disease worldwide that severely affects the reproductive performance of sows. However, little is still known about the genetic basis of reproductive performance in purebred herds during a PRRS outbreak through the comparison of maternal and terminal breeds. Thus, the objective of this work was to explore the host genetics of response to PRRS in purebred sows from two breeds. Reproductive data included 2546 Duroc and 2522 Landrace litters from 894 and 813 purebred sows, respectively, which had high-density genotype data available (29,799 single nucleotide polymorphisms; SNPs). The data were split into pre-PRRS, PRRS, and post-PRRS phases based on standardized farrow-year-week estimates. Heritability estimates for reproductive traits were low to moderate (≤0.20) for Duroc and Landrace across PRRS phases. On the other hand, genetic correlations of reproductive traits between PRRS phases were overall moderate to high for both breeds. Several associations between MARC0034894, a candidate SNP for response to PRRS, with reproductive performance were identified (P-value < 0.05). Genomic analyses detected few QTL for reproductive performance across all phases, most explaining a small percentage of the additive genetic variance (≤8.2%, averaging 2.1%), indicating that these traits are highly polygenic. None of the identified QTL within a breed and trait overlapped between PRRS phases. Overall, our results indicate that Duroc sows are phenotypically more resilient to PRRS than Landrace sows, with a similar return to PRRS-free performance between breeds for most reproductive traits. Genomic prediction results indicate that genomic selection for improved reproductive performance under a PRRS outbreak is possible, especially in Landrace sows, by training markers using data from PRRS-challenged sows. On the other hand, the high genetic correlations with reproductive traits between PRRS phases suggest that selection for improved reproductive performance in a clean environment could improve performance during PRRS, but with limited efficiency due to their low heritability estimates. Thus, we hypothesize that an indicator trait that could be indirectly selected to increase the response to selection for these traits would be desirable and would also improve the reproductive performance of sows during a PRRS outbreak.

Host Genetics of Response to Porcine Reproductive and Respiratory Syndrome in Sows: Antibody Response as an Indicator Trait for Improved Reproductive Performance.

Antibody response to porcine reproductive and respiratory syndrome (PRRS) virus (PRRSV) infection, measured as sample-to-positive (S/P) ratio, has been proposed as an indicator trait for improved reproductive performance during a PRRS outbreak in Landrace sows. However, this result has not yet been validated in Landrace sows or evaluated in terminal sire lines. The main objectives of this work were to validate the use of S/P ratio as an indicator trait to select pigs during a PRRS outbreak and to explore the genetic basis of antibody response to PRRSV. Farrowing data included 2,546 and 2,522 litters from 894 Duroc and 813 Landrace sows, respectively, split into pre-PRRS, PRRS, and post-PRRS phases. Blood samples were taken from 1,231 purebred sows (541 Landrace and 690 Duroc) following a PRRS outbreak for subsequent PRRSV ELISA analysis for S/P ratio measurement. All animals had high-density genotype data available (29,799 single nucleotide polymorphisms; SNPs). Genetic parameters and genome-wide association studies (GWAS) for S/P ratio were performed for each breed separately. Heritability estimates (± standard error) of S/P ratio during the PRRS outbreak were moderate, with 0.35 ± 0.08 for Duroc and 0.34 ± 0.09 for Landrace. During the PRRS outbreak, favorable genetic correlations of S/P ratio with the number of piglets born alive (0.61 ± 0.34), number of piglets born dead (-0.33 ± 0.32), and number of stillborn piglets (-0.27 ± 0.31) were observed for Landrace sows. For Duroc, the GWAS identified a major quantitative trait locus (QTL) on chromosome (Chr) 7 (24-15 megabases; Mb) explaining 15% of the total genetic variance accounted for by markers (TGVM), and another one on Chr 8 (25 Mb) explaining 2.4% of TGVM. For Landrace, QTL on Chr 7 (24-25 Mb) and Chr 7 (108-109 Mb), explaining 31% and 2.2% of TGVM, respectively, were identified. Some of the SNPs identified in these regions for S/P ratio were associated with reproductive performance but not during the PRRS outbreak. Genomic prediction accuracies for S/P ratio were moderate to high for the within-breed analysis. For the between-breed analysis, these were overall low. These results further support the use of S/P ratio as an indicator trait for improved reproductive performance during a PRRS outbreak in Landrace sows.

Whole Genome Sequence Data Provides Novel Insights Into the Genetic Architecture of Meat Quality Traits in Beef.

Tenderness is a major quality attribute for fresh beef steaks in the United States, and meat quality traits in general are suitable candidates for genomic research. The objectives of the present analysis were to (1) perform genome-wide association (GWA) analysis for marbling, Warner-Bratzler shear force (WBSF), tenderness, and connective tissue using whole-genome data in an Angus population, (2) identify enriched pathways in each GWA analysis; (3) construct a protein-protein interaction network using the associated genes and (4) perform a µ-calpain proteolysis assessment for associated structural proteins. An Angus-sired population of 2,285 individuals was assessed. Animals were transported to a commercial packing plant and harvested at an average age of 457 ± 46 days. After 48 h postmortem, marbling was recorded by graders' visual appraisal. Two 2.54-cm steaks were sampled from each muscle for recording of WBSF, and tenderness, and connective tissue by a sensory panel. The relevance of additive effects on marbling, WBSF, tenderness, and connective tissue was evaluated on a genome-wide scale using a two-step mixed model-based approach in single-trait analysis. A tissue-restricted gene enrichment was performed for each GWA where all polymorphisms with an association p-value lower than 1 × 10-3 were included. The genes identified as associated were included in a protein-protein interaction network and a candidate structural protein assessment of proteolysis analyses. A total of 1,867, 3,181, 3,926, and 3,678 polymorphisms were significantly associated with marbling, WBSF, tenderness, and connective tissue, respectively. The associate region on BTA29 (36,432,655-44,313,046 bp) harbors 13 highly significant markers for meat quality traits. Enrichment for the GO term GO:0005634 (Nucleus), which includes transcription factors, was evident. The final protein-protein network included 431 interations between 349 genes. The 42 most important genes based on significance that encode structural proteins were included in a proteolysis analysis, and 81% of these proteins were potential µ-Calpain substrates. Overall, this comprehensive study unraveled genetic variants, genes and mechanisms of action responsible for the variation in meat quality traits. Our findings can provide opportunities for improving meat quality in beef cattle via marker-assisted selection.

Genetic Basis of Blood-Based Traits and Their Relationship With Performance and Environment in Beef Cattle at Weaning.

The objectives of this study were to explore the usefulness of blood-based traits as indicators of health and performance in beef cattle at weaning and identify the genetic basis underlying the different blood parameters obtained from complete blood counts (CBCs). Disease costs represent one of the main factors determining profitability in animal production. Previous research has observed associations between blood cell counts and an animal's health status in some species. CBC were recorded from approximately 570 Angus based, crossbred beef calves at weaning born between 2015 and 2016 and raised on toxic or novel tall fescue. The calves (N = ∼600) were genotyped at a density of 50k SNPs and the genotypes (N = 1160) were imputed to a density of 270k SNPs. Genetic parameters were estimated for 15 blood and 4 production. Finally, with the objective of identifying the genetic basis underlying the different blood-based traits, genome-wide association studies (GWAS) were performed for all traits. Heritability estimates ranged from 0.11 to 0.60, and generally weak phenotypic correlations and strong genetic correlations were observed among blood-based traits only. Genome-wide association study identified ninety-one 1-Mb windows that accounted for 0.5% or more of the estimated genetic variance for at least 1 trait with 21 windows overlapping across two or more traits (explaining more than 0.5% of estimated genetic variance for two or more traits). Five candidate genes have been identified in the most interesting overlapping regions related to blood-based traits. Overall, this study represents one of the first efforts represented in scientific literature to identify the genetic basis of blood cell traits in beef cattle. The results presented in this study allow us to conclude that: (1) blood-based traits have weak phenotypic correlations but strong genetic correlations among themselves. (2) Blood-based traits have moderate to high heritability. (3) There is evidence of an important overlap of genetic control among similar blood-based traits which will allow for their use in improvement programs in beef cattle.

Haplotype-based genome-wide association studies for carcass and growth traits in chicken.

There have been several genome-wide association study (GWAS) reported for carcass, growth, and meat traits in chickens. Most of these studies have been based on single SNPs GWAS. In contrast, haplotype-based GWAS reports have been limited. In the present study, 2 Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF) and genotyped with the chicken 60K SNP chip were used to perform a haplotype-based GWAS. The lean and fat chicken lines were selected for abdominal fat content for 11 yr. Abdominal fat weight was significantly different between the 2 lines; however, there was no difference for body weight between the lean and fat lines. A total of 132 haplotype windows were significantly associated with abdominal fat weight. These significantly associated haplotype windows were primarily located on chromosomes 2, 4, 8, 10, and 26. Seven candidate genes, including SHH, LMBR1, FGF7, IL16, PLIN1, IGF1R, and SLC16A1, were located within these associated regions. These genes may play important roles in the control of abdominal fat content. Two regions on chromosomes 3 and 10 were significantly associated with testis weight. These 2 regions were previously detected by the single SNP GWAS using this same resource population. TCF21 on chromosome 3 was identified as a potentially important candidate gene for testis growth and development based on gene expression analysis and the reported function of this gene. TCF12, which was previously detected in our SNP by SNP interaction analysis, was located in a region on chromosome 10 that was significantly associated with testis weight. Six candidate genes, including TNFRSF1B, PLOD1, NPPC, MTHFR, EPHB2, and SLC35A3, on chromosome 21 may play important roles in bone development based on the known function of these genes. In addition, several regions were significantly associated with other carcass and growth traits, but no candidate genes were identified. The results of the present study may be helpful in understanding the genetic mechanisms of carcass and growth traits in chickens.

A Vision for Development and Utilization of High-Throughput Phenotyping and Big Data Analytics in Livestock.

Automated high-throughput phenotyping with sensors, imaging, and other on-farm technologies has resulted in a flood of data that are largely under-utilized. Drastic cost reductions in sequencing and other omics technology have also facilitated the ability for deep phenotyping of livestock at the molecular level. These advances have brought the animal sciences to a cross-roads in data science where increased training is needed to manage, record, and analyze data to generate knowledge and advances in Agriscience related disciplines. This paper describes the opportunities and challenges in using high-throughput phenotyping, "big data," analytics, and related technologies in the livestock industry based on discussions at the Livestock High-Throughput Phenotyping and Big Data Analytics meeting, held in November 2017 (see: https://www.animalgenome.org/bioinfo/community/workshops/2017/). Critical needs for investments in infrastructure for people (e.g., "big data" training), data (e.g., data transfer, management, and analytics), and technology (e.g., development of low cost sensors) were defined by this group. Though some subgroups of animal science have extensive experience in predictive modeling, cross-training in computer science, statistics, and related disciplines are needed to use big data for diverse applications in the field. Extensive opportunities exist for public and private entities to harness big data to develop valuable research knowledge and products to the benefit of society under the increased demands for food in a rapidly growing population.

Identification of genome-wide SNP-SNP interactions associated with important traits in chicken.

BACKGROUND: In addition to additive genetic effects, epistatic interactions can play key roles in the control of phenotypic variation of traits of interest. In the current study, 475 male birds from lean and fat chicken lines were utilized as a resource population to detect significant epistatic effects associated with growth and carcass traits. RESULTS: A total of 421 significant epistatic effects were associated with testis weight (TeW), from which 11 sub-networks (Sub-network1 to Sub-network11) were constructed. In Sub-network1, which was the biggest network, there was an interaction between GGA21 and GGAZ. Three genes on GGA21 (SDHB, PARK7 and VAMP3) and nine genes (AGTPBP1, CAMK4, CDC14B, FANCC, FBP1, GNAQ, PTCH1, ROR2 and STARD4) on GGAZ that might be potentially important candidate genes for testis growth and development were detected based on the annotated gene function. In Sub-network2, there was a SNP on GGA19 that interacted with 8 SNPs located on GGA10. The SNP (Gga_rs15834332) on GGA19 was located between C-C motif chemokine ligand 5 (CCL5) and MIR142. There were 32 Refgenes on GGA10, including TCF12 which is predicted to be a target gene of miR-142-5p. We hypothesize that miR-142-5p and TCF12 may interact with one another to regulate testis growth and development. Two genes (CDH12 and WNT8A) in the same cadherin signaling pathway were implicated as potentially important genes in the control of metatarsus circumference (MeC). There were no significant epistatic effects identified for the other carcass and growth traits, e.g. heart weight (HW), liver weight (LW), spleen weight (SW), muscular and glandular stomach weight (MGSW), carcass weight (CW), body weight (BW1, BW3, BW5, BW7), chest width (ChWi), metatarsus length (MeL). CONCLUSIONS: The results of the current study are helpful to better understand the genetic basis of carcass and growth traits, especially for testis growth and development in broilers.