A Pilot Detection and Associate Study of Gene Presence-Absence Variation in Holstein Cattle.

Presence-absence variations (PAVs) are important structural variations, wherein a genomic segment containing one or more genes is present in some individuals but absent in others. While PAVs have been extensively studied in plants, research in cattle remains limited. This study identified PAVs in 173 Holstein bulls using whole-genome sequencing data and assessed their associations with 46 economically important traits. Out of 28,772 cattle genes (from the longest transcripts), a total of 26,979 (93.77%) core genes were identified (present in all individuals), while variable genes included 928 softcore (present in 95-99% of individuals), 494 shell (present in 5-94%), and 371 cloud genes (present in <5%). Cloud genes were enriched in functions associated with hormonal and antimicrobial activities, while shell genes were enriched in immune functions. PAV-based genome-wide association studies identified associations between gene PAVs and 16 traits including milk, fat, and protein yields, as well as traits related to health and reproduction. Associations were found on multiple chromosomes, illustrating important associations on cattle chromosomes 7 and 15, involving olfactory receptor and immune-related genes, respectively. By examining the PAVs at the population level, the results of this research provided crucial insights into the genetic structures underlying the complex traits of Holstein cattle.

Systemic interindividual DNA methylation variants in cattle share major hallmarks with those in humans.

BACKGROUND: We recently identified ~ 10,000 correlated regions of systemic interindividual epigenetic variation (CoRSIVs) in the human genome. These methylation variants are amenable to population studies, as DNA methylation measurements in blood provide information on epigenetic regulation throughout the body. Moreover, establishment of DNA methylation at human CoRSIVs is labile to periconceptional influences such as nutrition. Here, we analyze publicly available whole-genome bisulfite sequencing data on multiple tissues of each of two Holstein cows to determine whether CoRSIVs exist in cattle. RESULTS: Focusing on genomic blocks with ≥ 5 CpGs and a systemic interindividual variation index of at least 20, our approach identifies 217 cattle CoRSIVs, a subset of which we independently validate by bisulfite pyrosequencing. Similar to human CoRSIVs, those in cattle are strongly associated with genetic variation. Also as in humans, we show that establishment of DNA methylation at cattle CoRSIVs is particularly sensitive to early embryonic environment, in the context of embryo culture during assisted reproduction. CONCLUSIONS: Our data indicate that CoRSIVs exist in cattle, as in humans, suggesting these systemic epigenetic variants may be common to mammals in general. To the extent that individual epigenetic variation at cattle CoRSIVs affects phenotypic outcomes, assessment of CoRSIV methylation at birth may become an important tool for optimizing agriculturally important traits. Moreover, adjusting embryo culture conditions during assisted reproduction may provide opportunities to tailor agricultural outcomes by engineering CoRSIV methylation profiles.

Large-Sample Genome-Wide Association Study of Resistance to Retained Placenta in U.S. Holstein Cows.

A genome-wide association study of resistance to retained placenta (RETP) using 632,212 Holstein cows and 74,747 SNPs identified 200 additive effects with p-values < 10-8 on thirteen chromosomes but no dominance effect was statistically significant. The regions of 87.61-88.74 Mb of Chr09 about 1.13 Mb in size had the most significant effect in LOC112448080 and other highly significant effects in CCDC170 and ESR1, and in or near RMND1 and AKAP12. Four non-ESR1 genes in this region were reported to be involved in ESR1 fusions in humans. Chr23 had the largest number of significant effects that peaked in SLC17A1, which was involved in urate metabolism and transport that could contribute to kidney disease. The PKHD1 gene contained seven significant effects and was downstream of another six significant effects. The ACOT13 gene also had a highly significant effect. Both PKHD1 and ACOT13 were associated with kidney disease. Another highly significant effect was upstream of BOLA-DQA2. The KITLG gene of Chr05 that acts in utero in germ cell and neural cell development, and hematopoiesis was upstream of a highly significant effect, contained a significant effect, and was between another two significant effects. The results of this study provided a new understanding of genetic factors underlying RETP in U.S. Holstein cows.

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.

A Million-Cow Validation of a Chromosome 14 Region Interacting with All Chromosomes for Fat Percentage in U.S. Holstein Cows.

A genome-wide association study (GWAS) of fat percentage (FPC) using 1,231,898 first lactation cows and 75,198 SNPs confirmed a previous result that a Chr14 region about 9.38 Mb in size (0.14-9.52 Mb) had significant inter-chromosome additive × additive (A×A) effects with all chromosomes and revealed many new such effects. This study divides this 9.38 Mb region into two sub-regions, Chr14a at 0.14-0.88 Mb (0.74 Mb in size) with 78% and Chr14b at 2.21-9.52 Mb (7.31 Mb in size) with 22% of the 2761 significant A×A effects. These two sub-regions were separated by a 1.3 Mb gap at 0.9-2.2 Mb without significant inter-chromosome A×A effects. The PPP1R16A-FOXH1-CYHR1-TONSL (PFCT) region of Chr14a (29 Kb in size) with four SNPs had the largest number of inter-chromosome A×A effects (1141 pairs) with all chromosomes, including the most significant inter-chromosome A×A effects. The SLC4A4-GC-NPFFR2 (SGN) region of Chr06, known to have highly significant additive effects for some production, fertility and health traits, specifically interacted with the PFCT region and a Chr14a region with CPSF1, ADCK5, SLC52A2, DGAT1, SMPD5 and PARP10 (CASDSP) known to have highly significant additive effects for milk production traits. The most significant effects were between an SNP in SGN and four SNPs in PFCT. The CASDSP region mostly interacted with the SGN region. In the Chr14b region, the 2.28-2.42 Mb region (138.46 Kb in size) lacking coding genes had the largest cluster of A×A effects, interacting with seventeen chromosomes. The results from this study provide high-confidence evidence towards the understanding of the genetic mechanism of FPC in Holstein cows.

Participatory investigation of goat farmers' breeding practices, trait preference, and selection criteria in Burkina Faso.

The community-based breeding program (CBBP) is an innovative approach recommended for genetic improvement and sustainable use of animal genetic resources in extensive farming systems. Successful implementation of this approach requires an understanding of the characteristics of production systems, breeding objectives, and farmers' trait preference. This study aimed to identify the selection criteria of goat farmers in rural areas of Burkina Faso and their potential implications in establishing CBBP. Following focus group discussions, a well-structured questionnaire was designed and administered to 372 randomly selected goat farmers in two different agro-ecological zones. A list of traits obtained during focus group discussions was provided to farmers individually, and they were asked to rank the ones they preferentially use to select breeding animals. Statistical tests were conducted to compare data between the two agro-ecological zones. The results showed that the average goat flock per household was higher (P < 0.05) in the Sudanian (15.68 ± 13.76), compared to the Sudano-Sahelian area (12.93 ± 13.3). Adult females were the dominant age-sex group in both areas. Reasons for culling, keeping breeding bucks, and castration practice were significantly different (P < 0.05) among agro-ecological zones. The most important common criterion for selection in the two zones was body size, coat color, and growth rate for the bucks and does, while fertility (0.06) parameters including twining ability (0.18), kidding frequency (0.11), and mothering ability (0.15) were furthermore considered for breeding does selection. These findings provide valuable insights for developing CBBPs tailored to goat production in the study areas.

Updating test-day milk yield factors for use in genetic evaluations and dairy production systems: a comprehensive review.

Various methods have been proposed to estimate daily yield from partial yields, primarily to deal with unequal milking intervals. This paper offers an exhaustive review of daily milk yields, the foundation of lactation records. Seminal advancements in the late 20th century concentrated on two main adjustment metrics: additive additive correction factors (ACF) and multiplicative correction factors (MCF). An ACF model provides additive adjustments to two times AM or PM milk yield, which then becomes the estimated daily yields, whereas an MCF is a ratio of daily yield to the yield from a single milking. Recent studies highlight the potential of alternative approaches, such as exponential regression and other nonlinear models. Biologically, milk secretion rates are not linear throughout the entire milking interval, influenced by the internal mammary gland pressure. Consequently, nonlinear models are appealing for estimating daily milk yields as well. MCFs and ACFs are typically determined for discrete milking interval classes. Nonetheless, large discrete intervals can introduce systematic biases. A universal solution for deriving continuous correction factors has been proposed, ensuring reduced bias and enhanced daily milk yield estimation accuracy. When leveraging test-day milk yields for genetic evaluations in dairy cattle, two predominant statistical models are employed: lactation and test-day yield models. A lactation model capitalizes on the high heritability of total lactation yields, aligning closely with dairy producers' needs because the total amount of milk production in a lactation directly determines farm revenue. However, a lactation yield model without harnessing all test-day records may ignore vital data about the shapes of lactation curves needed for informed breeding decisions. In contrast, a test-day model emphasizes individual test-day data, accommodating various intervals and recording plans and allowing the estimation of environmental effects on specific test days. In the United States, the patenting of test-day models in 1993 used to restrict the use of test-day models to regional and unofficial evaluations by the patent holders. Estimated test-day milk yields have been used as if they were accurate depictions of actual milk yields, neglecting possible estimation errors. Its potential consequences on subsequent genetic evaluations have not been sufficiently addressed. Moving forward, there are still numerous questions and challenges in this domain.

The African Goat Improvement Network: a scientific group empowering smallholder farmers.

The African Goat Improvement Network (AGIN) is a collaborative group of scientists focused on genetic improvement of goats in small holder communities across the African continent. The group emerged from a series of workshops focused on enhancing goat productivity and sustainability. Discussions began in 2011 at the inaugural workshop held in Nairobi, Kenya. The goals of this diverse group were to: improve indigenous goat production in Africa; characterize existing goat populations and to facilitate germplasm preservation where appropriate; and to genomic approaches to better understand adaptation. The long-term goal was to develop cost-effective strategies to apply genomics to improve productivity of small holder farmers without sacrificing adaptation. Genome-wide information on genetic variation enabled genetic diversity studies, facilitated improved germplasm preservation decisions, and provided information necessary to initiate large scale genetic improvement programs. These improvements were partially implemented through a series of community-based breeding programs that engaged and empowered local small farmers, especially women, to promote sustainability of the production system. As with many international collaborative efforts, the AGIN work serves as a platform for human capacity development. This paper chronicles the evolution of the collaborative approach leading to the current AGIN organization and describes how it builds capacity for sustained research and development long after the initial program funds are gone. It is unique in its effectiveness for simultaneous, multi-level capacity building for researchers, students, farmers and communities, and local and regional government officials. The positive impact of AGIN capacity building has been felt by participants from developing, as well as developed country partners.

Does modeling causal relationships improve the accuracy of predicting lactation milk yields?

This study compared 3 correlational (best prediction, linear regression, and feed-forward neural networks) and 2 causal models (recursive structural equation model and recurrent neural networks) for estimating lactation milk yields. The correlational models assumed associations between test-day milk yields (health conditions), while the casual models postulated unidirectional recursive effects between these test-day variables. Wood lactation curves were used to simulate the data and served as a benchmark model. Individual Wood lactation curves provided an excellent parametric interpretation of lactation dynamics, with their prediction accuracies depending on the coverage of the lactation curve dynamics. Best prediction outperformed other models in the absence of mastitis but was suboptimal when mastitis was present and unaccounted for. Recurrent neural networks yielded the highest accuracy when mastitis was present. Although causal models facilitated the inference about the causality underlying lactation, precisely capturing the causal relationships was challenging because the underlying biology was complex. Misspecification of recursive effects in the recursive structural equation model resulted in a loss of accuracy. Hence, modeling causal relationships does not necessarily guarantee improved accuracies. In practice, a parsimonious model is preferred, balancing model complexity and accuracy. In addition to the choice of statistical models, the proper accounting for factors and covariates affecting milk yields is equally crucial.

Using the community-based breeding program (CBBP) model as a collaborative platform to develop the African Goat Improvement Network-Image collection protocol (AGIN-ICP) with mobile technology for data collection and management of livestock phenotypes.

Introduction: The African Goat Improvement Network Image Collection Protocol (AGIN-ICP) is an accessible, easy to use, low-cost procedure to collect phenotypic data via digital images. The AGIN-ICP collects images to extract several phenotype measures including health status indicators (anemia status, age, and weight), body measurements, shapes, and coat color and pattern, from digital images taken with standard digital cameras or mobile devices. This strategy is to quickly survey, record, assess, analyze, and store these data for use in a wide variety of production and sampling conditions. Methods: The work was accomplished as part of the multinational African Goat Improvement Network (AGIN) collaborative and is presented here as a case study in the AGIN collaboration model and working directly with community-based breeding programs (CBBP). It was iteratively developed and tested over 3 years, in 12 countries with over 12,000 images taken. Results and discussion: The AGIN-ICP development is described, and field implementation and the quality of the resulting images for use in image analysis and phenotypic data extraction are iteratively assessed. Digital body measures were validated using the PreciseEdge Image Segmentation Algorithm (PE-ISA) and software showing strong manual to digital body measure Pearson correlation coefficients of height, length, and girth measures (0.931, 0.943, 0.893) respectively. It is critical to note that while none of the very detailed tasks in the AGIN-ICP described here is difficult, every single one of them is even easier to accidentally omit, and the impact of such a mistake could render a sample image, a sampling day's images, or even an entire sampling trip's images difficult or unusable for extracting digital phenotypes. Coupled with tissue sampling and genomic testing, it may be useful in the effort to identify and conserve important animal genetic resources and in CBBP genetic improvement programs by providing reliably measured phenotypes with modest cost. Potential users include farmers, animal husbandry officials, veterinarians, regional government or other public health officials, researchers, and others. Based on these results, a final AGIN-ICP is presented, optimizing the costs, ease, and speed of field implementation of the collection method without compromising the quality of the image data collection.

Estimating test-day milk yields by modeling proportional daily yields: Going beyond linearity.

In the United States, lactation milk yields are not measured directly but are calculated from the test-day milk yields. Still, test-day milk yields are estimated from partial yields obtained from single milkings. Various methods have been proposed to estimate test-day milk yields, primarily to deal with unequal milking intervals dating back to the 1970s and 1980s. The Wiggans model is a de facto method for estimating test-day milk yields in the United States, which was initially proposed for cows milked 3 times daily, assuming a linear relationship between a proportional test-day milk yield and milking interval. However, the linearity assumption did not hold precisely in Holstein cows milked twice daily because of prolonged and uneven milking intervals. The present study reviewed and evaluated the nonlinear models that extended the Wiggans model for estimating daily or test-day milk yields. These nonlinear models, except step functions, demonstrated smaller errors and greater accuracies for estimated test-day milk yields compared with the conventional methods. The nonlinear models offered additional benefits. For example, the locally weighted regression model (e.g., locally estimated scatterplot smoothing) could utilize data information in scalable neighborhoods and weigh observations according to their distance in milking interval time. General additive models provide a flexible, unified framework to model nonlinear predictor variables additively. Another drawback of the conventional methods is a loss of accuracy caused by discretizing milking interval time into large bins while deriving multiplicative correction factors for estimating test-day milk yields. To overcome this problem, we proposed a general approach that allows milk yield correction factors to be derived for every possible milking interval time, resulting in more accurately estimated test-day milk yields. This approach can be applied to any model, including nonparametric models.

Morphometric Characterization of Local Goat Breeds in Two Agroecological Zones of Burkina Faso, West Africa.

In Burkina Faso, goats are the second most numerous ruminant livestock population, with almost exclusively indigenous breeds being reared in extensive production systems in various agroecological zones. This study was carried out to understand the morphological variation of local goat breeds in the Sudano-Sahelian and Sudanian agroecological zones. A total of 511 adult female animals belonging to two presumed populations (Mossi breed in Sudano-Sahelian zone and Djallonké breed in Sudanian zone) were sampled and body weight as well as a range of linear body measurements, following FAO guidelines, were recorded. The least squares means of body measurements of indicated that Sudano-Sahelian goats have significantly (p < 0.001) larger body measurements than Sudanian goats. Furthermore, relative high variability of the two populations in morphometric traits was observed. Principal Component Analysis (PCA) suggested structure between Mossi breed on one side and Djallonké on the other side, but no strict separation was observed, suggesting that gene flow is occurring among the different populations. A dispersion map with four clusters was built based on the first two factors. The least square means of body measurements ranked the four groups from small to large body size, namely Djallonké, Mossi × Djallonké, Mossi, and Sahelian × Mossi. Gene flow from Sahelian goat into other populations of the country, based on migration of the Fulani ethnic group from the Sahel into areas with Mossi and Djallonké breeds, could explain this configuration and confirms the continuous erosion of genetic identity of these two local breeds. The sustainable use of these adapted local goat genetic resources calls for the promotion of sustainable genetic improvement using participatory breeding approaches.

Daily milk yield correction factors: What are they?

Cows are typically milked 2 or more times on a test-day, but not all these milkings are sampled and weighed. The initial approach estimated a test-day yield with doubled morning (AM) or evening (PM) yield in the AM-PM milking plans, assuming equal AM and PM milking intervals. However, AM and PM milking intervals can vary, and milk secretion rates may be different between day and night. Statistical methods have been proposed to estimate daily yields in dairy cows, focusing on various yield correction factors in 2 broad categories: additive correction factors (ACF) and multiplicative correction factors (MCF). The ACF are evaluated by the average differences between AM and PM milk yield for various milking interval classes, coupled with other categorical variables. We show that an ACF model is equivalent to a regression model of daily yield on categorical regressor variables, and a continuous variable for AM or PM yield with a fixed regression coefficient of 2.0. Similarly, a linear regression model can be implemented as an ACF model with the regression coefficient for AM or PM yield estimated from the data. The linear regression models improved the accuracy of the estimates compared with the ACF models. The MCF are ratios of daily yield to yield from single milkings, but their statistical interpretations vary. Overall, MCF were more accurate for estimating daily milk yield than ACF. The MCF have biological and statistical challenges. Systematic biases occurred when ACF or MCF were computed on discretized milking interval classes, leading to accuracy loss. An exponential regression model was proposed as an alternative model for estimating daily milk yields, which improved the accuracy. Characterization of ACF and MCF showed how they improved the accuracy compared with doubling AM or PM yield as the daily milk yield. All the methods performed similarly with equal AM and PM milkings. The methods were explicitly described to estimate daily milk yield in AM and PM milking plans. Still, the principles generally apply to cows milked more than 2 times a day and apply similarly to the estimation of daily fat and protein yields with some necessary modifications.

Invited review: Reliability computation from the animal model era to the single-step genomic model era.

The calculation of exact reliabilities involving the inversion of mixed model equations poses a heavy computational challenge when the system of equations is large. This has prompted the development of different approximation methods. We give an overview of the various methods and computational approaches in calculating reliability from the era before the animal model to the era of single-step genomic models. The different methods are discussed in terms of modeling, development, and applicability in large dairy cattle populations. The paper also describes the problems faced in reliability computation. Many details dispersed throughout the literature are presented in this paper. It is clear that a universal solution applicable to every model and input data may not be possible, but we point out several efficient and accurate algorithms developed recently for a variety of very large genomic evaluations.

PreciseEdge raster RGB image segmentation algorithm reduces user input for livestock digital body measurements highly correlated to real-world measurements.

Computer vision is a tool that could provide livestock producers with digital body measures and records that are important for animal health and production, namely body height and length, and chest girth. However, to build these tools, the scarcity of labeled training data sets with uniform images (pose, lighting) that also represent real-world livestock can be a challenge. Collecting images in a standard way, with manual image labeling is the gold standard to create such training data, but the time and cost can be prohibitive. We introduce the PreciseEdge image segmentation algorithm to address these issues by employing a standard image collection protocol with a semi-automated image labeling method, and a highly precise image segmentation for automated body measurement extraction directly from each image. These elements, from image collection to extraction are designed to work together to yield values highly correlated to real-world body measurements. PreciseEdge adds a brief preprocessing step inspired by chromakey to a modified GrabCut procedure to generate image masks for data extraction (body measurements) directly from the images. Three hundred RGB (red, green, blue) image samples were collected uniformly per the African Goat Improvement Network Image Collection Protocol (AGIN-ICP), which prescribes camera distance, poses, a blue backdrop, and a custom AGIN-ICP calibration sign. Images were taken in natural settings outdoors and in barns under high and low light, using a Ricoh digital camera producing JPG images (converted to PNG prior to processing). The rear and side AGIN-ICP poses were used for this study. PreciseEdge and GrabCut image segmentation methods were compared for differences in user input required to segment the images. The initial bounding box image output was captured for visual comparison. Automated digital body measurements extracted were compared to manual measures for each method. Both methods allow additional optional refinement (mouse strokes) to aid the segmentation algorithm. These optional mouse strokes were captured automatically and compared. Stroke count distributions for both methods were not normally distributed per Kolmogorov-Smirnov tests. Non-parametric Wilcoxon tests showed the distributions were different (p< 0.001) and the GrabCut stroke count was significantly higher (p = 5.115 e-49), with a mean of 577.08 (std 248.45) versus 221.57 (std 149.45) with PreciseEdge. Digital body measures were highly correlated to manual height, length, and girth measures, (0.931, 0.943, 0.893) for PreciseEdge and (0.936, 0. 944, 0.869) for GrabCut (Pearson correlation coefficient). PreciseEdge image segmentation allowed for masks yielding accurate digital body measurements highly correlated to manual, real-world measurements with over 38% less user input for an efficient, reliable, non-invasive alternative to livestock hand-held direct measuring tools.

Assembly of a pangenome for global cattle reveals missing sequences and novel structural variations, providing new insights into their diversity and evolutionary history.

A cattle pangenome representation was created based on the genome sequences of 898 cattle representing 57 breeds. The pangenome identified 83 Mb of sequence not found in the cattle reference genome, representing 3.1% novel sequence compared with the 2.71-Gb reference. A catalog of structural variants developed from this cattle population identified 3.3 million deletions, 0.12 million inversions, and 0.18 million duplications. Estimates of breed ancestry and hybridization between cattle breeds using insertion/deletions as markers were similar to those produced by single nucleotide polymorphism-based analysis. Hundreds of deletions were observed to have stratification based on subspecies and breed. For example, an insertion of a Bov-tA1 repeat element was identified in the first intron of the APPL2 gene and correlated with cattle breed geographic distribution. This insertion falls within a segment overlapping predicted enhancer and promoter regions of the gene, and could affect important traits such as immune response, olfactory functions, cell proliferation, and glucose metabolism in muscle. The results indicate that pangenomes are a valuable resource for studying diversity and evolutionary history, and help to delineate how domestication, trait-based breeding, and adaptive introgression have shaped the cattle genome.

A multi-tissue atlas of regulatory variants in cattle.

Characterization of genetic regulatory variants acting on livestock gene expression is essential for interpreting the molecular mechanisms underlying traits of economic value and for increasing the rate of genetic gain through artificial selection. Here we build a Cattle Genotype-Tissue Expression atlas (CattleGTEx) as part of the pilot phase of the Farm animal GTEx (FarmGTEx) project for the research community based on 7,180 publicly available RNA-sequencing (RNA-seq) samples. We describe the transcriptomic landscape of more than 100 tissues/cell types and report hundreds of thousands of genetic associations with gene expression and alternative splicing for 23 distinct tissues. We evaluate the tissue-sharing patterns of these genetic regulatory effects, and functionally annotate them using multiomics data. Finally, we link gene expression in different tissues to 43 economically important traits using both transcriptome-wide association and colocalization analyses to decipher the molecular regulatory mechanisms underpinning such agronomic traits in cattle.

Statistical Methods Revisited for Estimating Daily Milk Yields: How Well do They Work?

Cost-effective milking plans have been adapted to supplement the standard supervised twice-daily monthly testing scheme since the 1960s. Various methods have been proposed to estimate daily milk yields (DMY), focusing on yield correction factors. The present study evaluated the performance of existing statistical methods, including a recently proposed exponential regression model, for estimating DMY using 10-fold cross-validation in Holstein and Jersey cows. The initial approach doubled the morning (AM) or evening (PM) yield as estimated DMY in AM-PM plans, assuming equal 12-h AM and PM milking intervals. However, in reality, AM milking intervals tended to be longer than PM milking intervals. Additive correction factors (ACF) provided additive adjustments beyond twice AM or PM yields. Hence, an ACF model equivalently assumed a fixed regression coefficient or a multiplier of "2.0" for AM or PM yields. Similarly, a linear regression model was viewed as an ACF model, yet it estimated the regression coefficient for a single milk yield from the data. Multiplicative correction factors (MCF) represented daily to partial milk yield ratios. Hence, multiplying a yield from single milking by an appropriate MCF gave a DMY estimate. The exponential regression model was analogous to an exponential growth function with the yield from single milking as the initial state and the rate of change tuned by a linear function of milking interval. In the present study, all the methods had high precision in the estimates, but they differed considerably in biases. Overall, the MCF and linear regression models had smaller squared biases and greater accuracies for estimating DMY than the ACF models. The exponential regression model had the greatest accuracies and smallest squared biases. Model parameters were compared. Discretized milking interval categories led to a loss of accuracy of the estimates. Characterization of ACF and MCF revealed their similarities and dissimilarities and biases aroused by unequal milking intervals. The present study focused on estimating DMY in AM-PM milking plans. Yet, the methods and relevant principles are generally applicable to cows milked more than two times a day.

Variants Within Genes EDIL3 and ADGRB3 are Associated With Divergent Fecal Egg Counts in Katahdin Sheep at Weaning.

Gastrointestinal nematodes (GIN) pose a severe threat to sheep production worldwide. Anthelmintic drug resistance coupled with growing concern regarding potential environmental effects of drug use have demonstrated the necessity of implementing other methods of GIN control. The aim of this study was to test for genetic variants associated with resistance or susceptibility to GIN in Katahdin sheep to improve the current understanding of the genetic mechanisms responsible for host response to GIN. Linear regression and case-control genome-wide association studies were conducted with high-density genotype data and cube-root transformed weaning fecal egg counts (tFEC) of 583 Katahdin sheep. The case-control GWAS identified two significant SNPs (P-values 1.49e-08 to 1.01e-08) within introns of the gene adhesion G protein-coupled receptor B3 (ADGRB3) associated with lower fecal egg counts. With linear regression, four significant SNPs (P-values 7.82e-08 to 3.34e-08) were identified within the first intron of the gene EGF-like repeats and discoidin domains 3 (EDIL3). These identified SNPs were in very high linkage disequilibrium (r 2 of 0.996-1), and animals with alternate homozygous genotypes had significantly higher median weaning tFEC phenotypes compared to all other genotypes. Significant SNPs were queried through public databases to identify putative transcription factor binding site (TFBS) and potential lncRNA differences between reference and alternate alleles. Changes in TFBS were predicted at two SNPs, and one significant SNP was found to be within a predicted lncRNA sequence with greater than 90% similarity to a known lncRNA in the bovine genome. The gene EDIL3 has been described in other species for its roles in the inhibition and resolution of inflammation. Potential changes of EDIL3 expression mediated through lncRNA expression and/or transcription factor binding may impact the overall immune response and reduce the ability of Katahdin sheep to control GIN infection. This study lays the foundation for further research of EDIL3 and ADGRB3 towards understanding genetic mechanisms of susceptibility to GIN, and suggests these SNPs may contribute to genetic strategies for improving parasite resistance traits in sheep.

Genome-wide recombination map construction from single sperm sequencing in cattle.

BACKGROUND: Meiotic recombination is one of the important phenomena contributing to gamete genome diversity. However, except for human and a few model organisms, it is not well studied in livestock, including cattle. RESULTS: To investigate their distributions in the cattle sperm genome, we sequenced 143 single sperms from two Holstein bulls. We mapped meiotic recombination events at high resolution based on phased heterozygous single nucleotide polymorphism (SNP). In the absence of evolutionary selection pressure in fertilization and survival, recombination events in sperm are enriched near distal chromosomal ends, revealing that such a pattern is intrinsic to the molecular mechanism of meiosis. Furthermore, we further validated these findings in single sperms with results derived from sequencing its family trio of diploid genomes and our previous studies of recombination in cattle. CONCLUSIONS: To our knowledge, this is the first large-scale single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of recombination, genome instability, and male infertility.