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.

Image analysis to automatically classify anemia based on Famacha score in sheep using ocular conjunctiva images.

Haemonchus contortus is the most pathogenic blood-feeding parasitic in sheep, causing anemia and consequently changes in the color of the ocular conjunctiva, from the deep red of healthy sheep to shades of pink to practically white of non-healthy sheep. In this context, the Famacha method has been created for detecting sheep unable to cope with the infection by H. contortus, through visual assessment of ocular conjunctiva coloration. Thus, the objectives of this study were (1) to extract ocular conjunctiva image features to automatically classify Famacha score and compare two classification models (multinomial logistic regression-MLR and random forest-RF) and (2) to evaluate the applicability of the best classification model on three sheep farms. The dataset consisted of 1,156 ocular conjunctiva images from 422 animals. RF model was used to segment the images, i.e., to select the pixels that belong to the ocular conjunctiva. After segmentation, the quantiles (1%, 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, and 99%) of color intensity in each image channel (red, blue, and green) were determined and used as explanatory variables in the classification models, and the Famacha scores 1 (non-anemic) to 5 (severely anemic) were the target classes to be predicted (scores 1 to 5, with 162, 255, 443, 266, and 30 images, respectively). For objective 1, the performance metrics (precision and sensitivity) were obtained using MLR and RF models considering data from all farms randomly split. For objective 2, a leave-one-farm-out cross-validation technique was used to assess prediction quality across three farms (farms A, B, and C, with 726, 205, and 225 images, respectively). The RF provided the best performances in predicting anemic animals, as indicated by the high values of sensitivity for Famacha score 3 (80.9%), 4 (46.2%), and 5 (60%) compared to the MLR model. The precision of the RF was 72.7% for Famacha score 1 and 62.5% for Famacha score 2. These results indicate that is possible to successfully predict Famacha score, especially for scores 2 to 4, in sheep via image analysis and RF model using ocular conjunctiva images collected in farm conditions. As expected, model validation excluding entire farms in cross-validation presented a lower prediction quality. Nonetheless, this setup is closer to reality because the developed models are supposed to be used across farms, including new ones, and with different environments and management conditions.

Exploring associations among morphometric measurements, genetic group of sire, and performance of beef on dairy calves.

Sire selection for beef on dairy crosses plays an important role in livestock systems as it may affect future performance and carcass traits of growing and finishing crossbred cattle. The phenotypic variation found in beef on dairy crosses has raised concerns from meat packers due to animals with dairy-type carcass characteristics. The use of morphometric measurements may help to understand the phenotypic structures of sire progeny for selecting animals with greater performance. In addition, due to the relationship with growth, these measurements could be used to early predict the performance until the transition from dairy farms to sales. The objectives of this study were 1) to evaluate the effect of different beef sires and breeds on the morphometric measurements of crossbred calves including cannon bone (CB), forearm (FA), hip height (HH), face length (FL), face width (FW) and growth performance; and (2) to predict the weight gain from birth to transition from dairy farms to sale (WG) and the body weight at sale (BW) using such morphometric measurements obtained at first days of animals' life. CB, FA, HH, FL, FW, and weight at 7 ±â€…5 d (BW7) (Table 1) were measured on 206 calves, from four different sire breeds [Angus (AN), SimAngus (SA), Simmental (SI), and Limousin (LI)], from five farms. To evaluate the morphometric measurements at the transition from dairy farms to sale and animal performance 91 out of 206 calves sourced from four farms, and offspring of two different sires (AN and SA) were used. To predict the WG and BW, 97 calves, and offspring of three different sires (AN, SA, and LI) were used. The data were analyzed using a mixed model, considering farm and sire as random effects. To predict WG and BW, two linear models (including or not the morphometric measurements) were used, and a leave-one-out cross-validation strategy was used to evaluate their predictive quality. The HH and BW7 were 7.67% and 10.7% higher (P < 0.05) in SA crossbred calves compared to AN, respectively. However, the ADG and adjusted body weight to 120 d were 14.3% and 9.46% greater (P < 0.05) in AN compared to SA. The morphometric measurements improved the model's predictive performance for WG and BW. In conclusion, morphometric measurements at the first days of calves' life can be used to predict animals' performance in beef on dairy. Such a strategy could lead to optimized management decisions and greater profitability in dairy farms.

Using pseudo-labeling to improve performance of deep neural networks for animal identification.

Contemporary approaches for animal identification use deep learning techniques to recognize coat color patterns and identify individual animals in a herd. However, deep learning algorithms usually require a large number of labeled images to achieve satisfactory performance, which creates the need to manually label all images when automated methods are not available. In this study, we evaluated the potential of a semi-supervised learning technique called pseudo-labeling to improve the predictive performance of deep neural networks trained to identify Holstein cows using labeled training sets of varied sizes and a larger unlabeled dataset. By using such technique to automatically label previously unlabeled images, we observed an increase in accuracy of up to 20.4 percentage points compared to using only manually labeled images for training. Our final best model achieved an accuracy of 92.7% on an independent testing set to correctly identify individuals in a herd of 59 cows. These results indicate that it is possible to achieve better performing deep neural networks by using images that are automatically labeled based on a small dataset of manually labeled images using a relatively simple technique. Such strategy can save time and resources that would otherwise be used for labeling, and leverage well annotated small datasets.

Longitudinal analysis of bovine mammary gland development.

Many studies on bovine mammary glands focus on one stage of development. Often missing in those studies are repeated measures of development from the same animals. As milk production is directly affected by amount of parenchymal tissue within the udder, understanding mammary gland growth along with visualization of its structures during development is essential. Therefore, analysis of ultrasound and histology data from the same animals would result in better understanding of mammary development over time. Thus, this research aimed to describe mammary gland development using non-invasive and invasive tools to delineate growth rate of glandular tissue responsible for potential future milk production. Mammary gland ultrasound images, biopsy samples, and blood samples were collected from 36 heifer dairy calves beginning at 10 weeks of age, and evaluated at 26, 39, and 52 weeks. Parenchyma was quantified at 10 weeks of age using ultrasound imaging and histological evaluation, and average echogenicity was utilized to quantify parenchyma at later stages of development. A significant negative correlation was detected between average echogenicity of parenchyma at 10 weeks and total adipose as a percent of histological whole tissue at 52 weeks. Additionally, a negative correlation between average daily gain at 10 and 26 weeks and maximum echogenicity at 52 weeks was present. These results suggest average daily gain and mammary gland development prior to 39 weeks of age is associated with development of the mammary gland after 39 weeks. These findings could be predictors of future milk production, however this must be further explored.

Predicting ribeye area and circularity in live calves through 3D image analyses of body surface.

The use of sexed semen at dairy farms has improved heifer replacement over the last decade by allowing greater control over the number of retained females and enabling the selection of dams with superior genetics. Alternatively, beef semen can be used in genetically inferior dairy cows to produce crossbred (beef x dairy) animals that can be sold at a higher price. Although crossbreeding became profitable for dairy farmers, meat cuts from beef x dairy crosses often lack quality and shape uniformity. Technologies for quickly predicting carcass traits for animal grouping before harvest may improve meat cut uniformity in crossbred cattle. Our objective was to develop a deep learning approach for predicting ribeye area and circularity of live animals through 3D body surface images using two neural networks: 1) nested Pyramid Scene Parsing Network (nPSPNet) for extracting features and 2) Convolutional Neural Network (CNN) for estimating ribeye area and circularity from these features. A group of 56 calves were imaged using an Intel RealSense D435 camera. A total of 327 depth images were captured from 30 calves and labeled with masks outlining the calf body to train the nPSPNet for feature extraction. Additional 42,536 depth images were taken from the remaining 26 calves along with three ultrasound images collected for each calf from the 12/13th ribs. The ultrasound images (three by calf) were manually segmented to calculate the average ribeye area and circularity and then paired with the depth images for CNN training. We implemented a nested cross-validation approach, in which all images for one calf were removed (leave-one-out, LOO), and the remaining calves were further divided into training (70%) and validation (30%) sets within each LOO iteration. The proposed model predicted ribeye area with an average coefficient of determination (R2) of 0.74% and 7.3% mean absolute error of prediction (MAEP) and the ribeye circularity with an average R2 of 0.87% and 2.4% MAEP. Our results indicate that computer vision systems could be used to predict ribeye area and circularity in live animals, allowing optimal management decisions toward smart animal grouping in beef x dairy crosses and purebred.

This work proposes a method for predicting ribeye specific carcass traits of beef x dairy crossbred calves from 3D images using deep learning. This method completely automates the measurement of carcass traits, providing an efficient way of grouping calves for breeding programs or meat quality control.

Would large dataset sample size unveil the potential of deep neural networks for improved genome-enabled prediction of complex traits? The case for body weight in broilers.

BACKGROUND: Deep neural networks (DNN) are a particular case of artificial neural networks (ANN) composed by multiple hidden layers, and have recently gained attention in genome-enabled prediction of complex traits. Yet, few studies in genome-enabled prediction have assessed the performance of DNN compared to traditional regression models. Strikingly, no clear superiority of DNN has been reported so far, and results seem highly dependent on the species and traits of application. Nevertheless, the relatively small datasets used in previous studies, most with fewer than 5000 observations may have precluded the full potential of DNN. Therefore, the objective of this study was to investigate the impact of the dataset sample size on the performance of DNN compared to Bayesian regression models for genome-enable prediction of body weight in broilers by sub-sampling 63,526 observations of the training set. RESULTS: Predictive performance of DNN improved as sample size increased, reaching a plateau at about 0.32 of prediction correlation when 60% of the entire training set size was used (i.e., 39,510 observations). Interestingly, DNN showed superior prediction correlation using up to 3% of training set, but poorer prediction correlation after that compared to Bayesian Ridge Regression (BRR) and Bayes Cπ. Regardless of the amount of data used to train the predictive machines, DNN displayed the lowest mean square error of prediction compared to all other approaches. The predictive bias was lower for DNN compared to Bayesian models, across all dataset sizes, with estimates close to one with larger sample sizes. CONCLUSIONS: DNN had worse prediction correlation compared to BRR and Bayes Cπ, but improved mean square error of prediction and bias relative to both Bayesian models for genome-enabled prediction of body weight in broilers. Such findings, highlights advantages and disadvantages between predictive approaches depending on the criterion used for comparison. Furthermore, the inclusion of more data per se is not a guarantee for the DNN to outperform the Bayesian regression methods commonly used for genome-enabled prediction. Nonetheless, further analysis is necessary to detect scenarios where DNN can clearly outperform Bayesian benchmark models.

Infrared Spectrometry as a High-Throughput Phenotyping Technology to Predict Complex Traits in Livestock Systems.

High-throughput phenotyping technologies are growing in importance in livestock systems due to their ability to generate real-time, non-invasive, and accurate animal-level information. Collecting such individual-level information can generate novel traits and potentially improve animal selection and management decisions in livestock operations. One of the most relevant tools used in the dairy and beef industry to predict complex traits is infrared spectrometry, which is based on the analysis of the interaction between electromagnetic radiation and matter. The infrared electromagnetic radiation spans an enormous range of wavelengths and frequencies known as the electromagnetic spectrum. The spectrum is divided into different regions, with near- and mid-infrared regions being the main spectral regions used in livestock applications. The advantage of using infrared spectrometry includes speed, non-destructive measurement, and great potential for on-line analysis. This paper aims to review the use of mid- and near-infrared spectrometry techniques as tools to predict complex dairy and beef phenotypes, such as milk composition, feed efficiency, methane emission, fertility, energy balance, health status, and meat quality traits. Although several research studies have used these technologies to predict a wide range of phenotypes, most of them are based on Partial Least Squares (PLS) and did not considered other machine learning (ML) techniques to improve prediction quality. Therefore, we will discuss the role of analytical methods employed on spectral data to improve the predictive ability for complex traits in livestock operations. Furthermore, we will discuss different approaches to reduce data dimensionality and the impact of validation strategies on predictive quality.

Comparison of data analytics strategies in computer vision systems to predict pig body composition traits from 3D images.

Computer vision systems (CVS) have been shown to be a powerful tool for the measurement of live pig body weight (BW) with no animal stress. With advances in precision farming, it is now possible to evaluate the growth performance of individual pigs more accurately. However, important traits such as muscle and fat deposition can still be evaluated only via ultrasound, computed tomography, or dual-energy x-ray absorptiometry. Therefore, the objectives of this study were: 1) to develop a CVS for prediction of live BW, muscle depth (MD), and back fat (BF) from top view 3D images of finishing pigs and 2) to compare the predictive ability of different approaches, such as traditional multiple linear regression, partial least squares, and machine learning techniques, including elastic networks, artificial neural networks, and deep learning (DL). A dataset containing over 12,000 images from 557 finishing pigs (average BW of 120 ± 12 kg) was split into training and testing sets using a 5-fold cross-validation (CV) technique so that 80% and 20% of the dataset were used for training and testing in each fold. Several image features, such as volume, area, length, widths, heights, polar image descriptors, and polar Fourier transforms, were extracted from the images and used as predictor variables in the different approaches evaluated. In addition, DL image encoders that take raw 3D images as input were also tested. This latter method achieved the best overall performance, with the lowest mean absolute scaled error (MASE) and root mean square error for all traits, and the highest predictive squared correlation (R2). The median predicted MASE achieved by this method was 2.69, 5.02, and 13.56, and R2 of 0.86, 0.50, and 0.45, for BW, MD, and BF, respectively. In conclusion, it was demonstrated that it is possible to successfully predict BW, MD, and BF via CVS on a fully automated setting using 3D images collected in farm conditions. Moreover, DL algorithms simplified and optimized the data analytics workflow, with raw 3D images used as direct inputs, without requiring prior image processing.

Beef cattle responses to pre-grazing sward height and low level of energy supplementation on tropical pastures.

The objective of this study was to investigate the effects of energy supplementation and pre-grazing sward height on grazing behavior, nutrient intake, digestion, and metabolism of cattle in tropical pastures managed as a rotational grazing system. Eight rumen-cannulated Nellore steers (24 mo of age; 300 ± 6.0 kg body weight [BW]) were used in a replicated 4 × 4 Latin square design. Treatments consisted of two levels of energy supplementation (0% [none] or 0.3% of BW of ground corn on an as-fed basis) and two pre-grazing sward heights (25 cm [defined by 95% light interception (LI)] or 35 cm [defined by ≥ 97.5% LI]) constituting four treatments. Steers grazed Marandu Palisadegrass [Brachiaria brizantha Stapf. cv. Marandu] and post-grazing sward height was 15 cm for all treatments. Forage dry matter intake (DMI) was increased (P = 0.01) when sward height was 25 cm (1.86% vs. 1.32% BW) and decreased (P = 0.04) when 0.3% BW supplement was fed (1.79% vs. 1.38% BW). Total and digestible DMI were not affected by energy supplementation (P = 0.57) but were increased when sward height was 25 cm (P = 0.01). Steers grazing the 25-cm sward height treatment spent less time grazing and more time resting, took fewer steps between feeding stations, and had a greater bite rate compared with 35-cm height treatment (P < 0.05). Energy supplementation reduced grazing time (P = 0.02) but did not affect any other grazing behavior parameter (P = 0.11). Energy supplementation increased (P < 0.01) diet dry matter digestibility but had no effect on crude protein and neutral detergent fiber digestibilities (P = 0.13). Compared with 35-cm pre-grazing sward height, steers at 25 cm presented lower rumen pH (6.39 vs. 6.52) and greater rumen ammonia nitrogen (11.22 vs. 9.77 mg/dL) and N retention (49.7% vs. 20.8%, P < 0.05). The pre-grazing sward height of 25 cm improved harvesting efficiency and energy intake by cattle, while feeding 0.3% of BW energy supplement did not increase the energy intake of cattle on tropical pasture under rotational grazing.

Feeding the combination of essential oils and exogenous α-amylase increases performance and carcass production of finishing beef cattle.

Two experiments were conducted to evaluate the performance responses of finishing feedlot cattle to dietary addition of essential oils and exogenous enzymes. The treatments in each experiment consisted of (DM basis): MON-sodium monensin (26 mg/kg); BEO-a blend of essential oils (90 mg/kg); BEO+MON-a blend of essential oils plus monensin (90 mg/kg + 26 mg/kg, respectively); BEO+AM-a blend of essential oils plus exogenous α-amylase (90 mg/kg + 560 mg/kg, respectively); and BEO+AM+PRO-a blend of essential oils plus exogenous α-amylase and exogenous protease (90 mg/kg + 560 mg/kg + 840 mg/kg, respectively). Exp. 1 consisted of a 93-d finishing period using 300 Nellore bulls in a randomized complete block design. Animals fed BEO had higher DMI (P < 0.001) but similar feed efficiency to animals fed MON (P ≥ 0.98). Compared with MON, the combination of BEO+AM resulted in 810 g greater DMI (P = 0.001), 190 g greater average daily gain (P = 0.04), 18 kg heavier final body weight (P = 0.04), and 12 kg heavier hot carcass weight (P = 0.02), although feed efficiency was not significantly different between BEO+AM and MON (P = 0.89). Combining BEO+MON tended to decrease hot carcass weight compared with BEO alone (P = 0.08) but not compared with MON (P = 0.98). Treatments did not impact observed dietary net energy values (P ≥ 0.74) or the observed:expected net energy ratio (P ≥ 0.11). In Exp. 2, five ruminally cannulated Nellore steers were used to evaluate intake, apparent total tract digestibility of nutrients, and ruminal parameters in a 5 × 5 Latin square design. Feeding BEO increased the total tract digestibility of CP compared to MON (P = 0.03). Compared to MON, feeding the combination of BEO+MON increased the intake of CP (P = 0.04) and NDF (P = 0.05), with no effects on total tract digestibility of nutrients (P ≥ 0.56), except for a tendency (P = 0.09) to increase CP digestibility. Intakes of all nutrients measured, except for ether extract (P = 0.16) were greater in animals fed BEO+AM when compared with MON (P ≤ 0.03), with no differences on total tract nutrient digestibilities (P ≥ 0.11) between these two treatments. In summary, diets containing the BEO used herein enhanced DMI of growing-finishing feedlot cattle compared with a basal diet containing MON without impair feed efficiency. A synergism between BEO and AM was detected, further increasing cattle performance and carcass production compared to MON.

A novel automated system to acquire biometric and morphological measurements and predict body weight of pigs via 3D computer vision.

Computer vision applications in livestock are appealing since they enable measurement of traits of interest without the need to directly interact with the animals. This allows the possibility of multiple measurements of traits of interest with minimal animal stress. In the current study, an automated computer vision system was devised and evaluated for extraction of features of interest, as body measurements and shape descriptors, and prediction of body weight in pigs. From the 655 pigs that had data collected 580 had more than 5 frames recorded and were used for development of the predictive models. The cross-validation for the models developed with data from nursery and finishing pigs presented an R2 ranging from 0.86 (random selected image) to 0.94 (median of images truncated on the third quartile), whereas with the dataset without nursery pigs, the R2 estimates ranged from 0.70 (random selected image) to 0.84 (median of images truncated on the third quartile). However, overall the mean absolute error was lower for the models fitted without data on nursery animals. From the body measures extracted from the image, body volume, area, and length were the most informative for prediction of body weight. The inclusion of the remaining body measurements (width and heights) or shape descriptors to the model promoted significant improvement of the predictions, whereas the further inclusion of sex and line effects were not significant.