The domestic sheep (Ovis aries) is an important agricultural species raised for meat, wool, and milk across the world. A high-quality reference genome for this species enhances the ability to discover genetic mechanisms influencing biological traits. Furthermore, a high-quality reference genome allows for precise functional annotation of gene regulatory elements. The rapid advances in genome assembly algorithms and emergence of sequencing technologies with increasingly long reads provide the opportunity for an improved de novo assembly of the sheep reference genome. Short-read Illumina (55× coverage), long-read Pacific Biosciences (75× coverage), and Hi-C data from this ewe retrieved from public databases were combined with an additional 50× coverage of Oxford Nanopore data and assembled with canu v1.9. The assembled contigs were scaffolded using Hi-C data with Salsa v2.2, gaps filled with PBsuitev15.8.24, and polished with Nanopolish v0.12.5. After duplicate contig removal with PurgeDups v1.0.1, chromosomes were oriented and polished with 2 rounds of a pipeline that consisted of freebayes v1.3.1 to call variants, Merfin to validate them, and BCFtools to generate the consensus fasta. The ARS-UI_Ramb_v2.0 assembly is 2.63 Gb in length and has improved continuity (contig NG50 of 43.18 Mb), with a 19- and 38-fold decrease in the number of scaffolds compared with Oar_rambouillet_v1.0 and Oar_v4.0. ARS-UI_Ramb_v2.0 has greater per-base accuracy and fewer insertions and deletions identified from mapped RNA sequence than previous assemblies. The ARS-UI_Ramb_v2.0 assembly is a substantial improvement in contiguity that will optimize the functional annotation of the sheep genome and facilitate improved mapping accuracy of genetic variant and expression data for traits in sheep.

Unbalanced metabolic status in the weeks after calving predisposes dairy cows to metabolic and infectious diseases. Blood glucose, IGF-I, non-esterified fatty acids (NEFA) and β-hydroxybutyrate (BHB) are used as indicators of the metabolic status of cows. This work aims to (1) evaluate the potential of milk mid-IR spectra to predict these blood components individually and (2) to evaluate the possibility of predicting the metabolic status of cows based on the clustering of these blood components. Blood samples were collected from 241 Holstein cows on six experimental farms, at days 14 and 35 after calving. Blood samples were analyzed by reference analysis and metabolic status was defined by k-means clustering (k=3) based on the four blood components. Milk mid-IR analyses were undertaken on different instruments and the spectra were harmonized into a common standardized format. Quantitative models predicting blood components were developed using partial least squares regression and discriminant models aiming to differentiate the metabolic status were developed with partial least squares discriminant analysis. Cross-validations were performed for both quantitative and discriminant models using four subsets randomly constituted. Blood glucose, IGF-I, NEFA and BHB were predicted with respective R 2 of calibration of 0.55, 0.69, 0.49 and 0.77, and R 2 of cross-validation of 0.44, 0.61, 0.39 and 0.70. Although these models were not able to provide precise quantitative values, they allow for screening of individual milk samples for high or low values. The clustering methodology led to the sharing out of the data set into three groups of cows representing healthy, moderately impacted and imbalanced metabolic status. The discriminant models allow to fairly classify the three groups, with a global percentage of correct classification up to 74%. When discriminating the cows with imbalanced metabolic status from cows with healthy and moderately impacted metabolic status, the models were able to distinguish imbalanced group with a global percentage of correct classification up to 92%. The performances were satisfactory considering the variables are not present in milk, and consequently predicted indirectly. This work showed the potential of milk mid-IR analysis to provide new metabolic status indicators based on individual blood components or a combination of these variables into a global status. Models have been developed within a standardized spectral format, and although robustness should preferably be improved with additional data integrating different geographic regions, diets and breeds, they constitute rapid, cost-effective and large-scale tools for management and breeding of dairy cows.

The aim of this study was to identify genomic regions associated with 305‐day milk yield and lactation curve parameters on primiparous (n = 9,910) and multiparous (n = 11,158) Holstein cows. The SNP solutions were estimated using a weighted single‐step genomic BLUP approach and imputed high‐density panel (777k) genotypes. The proportion of genetic variance explained by windows of 50 consecutive SNP (with an average of 165 Kb) was calculated, and regions that accounted for more than 0.50% of the variance were used to search for candidate genes. Estimated heritabilities were 0.37, 0.34, 0.17, 0.12, 0.30 and 0.19, respectively, for 305‐day milk yield, peak yield, peak time, ramp, scale and decay for primiparous cows. Genetic correlations of 305‐day milk yield with peak yield, peak time, ramp, scale and decay in primiparous cows were 0.99, 0.63, 0.20, 0.97 and −0.52, respectively. The results identified three windows on BTA14 associated with 305‐day milk yield and the parameters of lactation curve in primi‐ and multiparous cows. Previously proposed candidate genes for milk yield supported by this work include GRINA, CYHR1, FOXH1, TONSL, PPP1R16A, ARHGAP39, MAF1, OPLAH and MROH1, whereas newly identified candidate genes are MIR2308, ZNF7, ZNF34, SLURP1, MAFA and KIFC2 (BTA14). The protein lipidation biological process term, which plays a key role in controlling protein localization and function, was identified as the most important term enriched by the identified genes.

Background: Since their domestication 10,500 years ago, goat populations with distinctive genetic backgrounds have adapted to a broad variety of environments and breeding conditions. The VarGoats project is an international 1000-genome resequencing program designed to understand the consequences of domestication and breeding on the genetic diversity of domestic goats and to elucidate how speciation and hybridization have modeled the genomes of a set of species representative of the genus Capra . Findings: A dataset comprising 652 sequenced goats and 507 public goat sequences, including 35 animals representing eight wild species, has been collected worldwide. We identified 74,274,427 single nucleotide polymorphisms (SNPs) and 13,607,850 insertion-deletions (InDels) by aligning these sequences to the latest version of the goat reference genome (ARS1). A Neighbor-joining tree based on Reynolds genetic distances showed that goats from Africa, Asia and Europe tend to group into independent clusters. Because goat breeds from Oceania and Caribbean (Creole) all derive from imported animals, they are distributed along the tree according to their ancestral geographic origin. Conclusions: We report on an unprecedented international effort to characterize the genome-wide diversity of domestic goats. This large range of sequenced individuals represents a unique opportunity to ascertain how the demographic and selection processes associated with post-domestication history have shaped the diversity of this species. Data generated for the project will also be extremely useful to identify deleterious mutations and polymorphisms with causal effects on complex traits, and thus will contribute to new knowledge that could be used in genomic prediction and genome-wide association studies.

The aim of this study was to describe metabolism of early-lactation dairy cows by clustering cows based on glucose, insulin-like growth factor I (IGF-I), free fatty acid, and β-hydroxybutyrate (BHB) using the k-means method. Predictive models for metabolic clusters were created and validated using 3 sets of milk biomarkers (milk metabolites and enzymes, glycans on the immunogamma globulin fraction of milk, and Fourier-transform mid-infrared spectra of milk). Metabolic clusters are used to identify dairy cows with a balanced or imbalanced metabolic profile. Around 14 and 35 d in milk, serum or plasma concentrations of BHB, free fatty acids, glucose, and IGF-I were determined. Cows with a favorable metabolic profile were grouped together in what was referred to as the "balanced" group (n = 43) and were compared with cows in what was referred to as the "other balanced" group (n = 64). Cows with an unfavorable metabolic profile were grouped in what was referred to as the "imbalanced" group (n = 19) and compared with cows in what was referred to as the "other imbalanced" group (n = 88). Glucose and IGF-I were higher in balanced compared with other balanced cows. Free fatty acids and BHB were lower in balanced compared with other balanced cows. Glucose and IGF-I were lower in imbalanced compared with other imbalanced cows. Free fatty acids and BHB were higher in imbalanced cows. Metabolic clusters were related to production parameters. There was a trend for a higher daily increase in fat- and protein-corrected milk yield in balanced cows, whereas that of imbalanced cows was higher. Dry matter intake and the daily increase in dry matter intake were higher in balanced cows and lower in imbalanced cows. Energy balance was continuously higher in balanced cows and lower in imbalanced cows. Weekly or twice-weekly milk samples were taken and milk metabolites and enzymes (milk glucose, glucose-6-phosphate, BHB, lactate dehydrogenase, N-acetyl-β-d-glucosaminidase, isocitrate), immunogamma globulin glycans (19 peaks), and Fourier-transform mid-infrared spectra (1,060 wavelengths reduced to 15 principal components) were determined. Milk biomarkers with or without additional cow information (days in milk, parity, milk yield features) were used to create predictive models for the metabolic clusters. Accuracy for prediction of balanced (80%) and imbalanced (88%) cows was highest using milk metabolites and enzymes combined with days in milk and parity. The results and models of the present study are part of the GplusE project and identify novel milk-based phenotypes that may be used as predictors for metabolic and performance traits in early-lactation dairy cows.

Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.