Cet article synthétise une réflexion collective sur les orientations génétiques des bovins laitiers et allaitants de France métropolitaine nécessaires à l’adaptation des animaux aux évolutions des systèmes de production à l’horizon 2040. Il présente les facteurs susceptibles de faire évoluer les productions bovines françaises (démographie des éleveurs, changement climatique, raréfaction des ressources, progrès techniques et évolutions technologiques, réglementations et politiques agricoles, libéralisation du marché, évolution de la demande et organisation des filières). Les systèmes de production bovins devraient néanmoins conserver une diversité importante dans un contexte économique et réglementaire plus contraint, de changement climatique structurant et de raréfaction des ressources non renouvelables. Cet article recommande d’axer la sélection sur la rentabilité économique globale du système plutôt que sur le rendement animal. Pour les bovins laitiers, il s’agira d’augmenter la longévité fonctionnelle, réduire le format des animaux et augmenter la persistance laitière. Pour les bovins à viande, il est souhaitable d’augmenter la précocité sexuelle et la précocité de développement, améliorer la capacité des animaux à valoriser l’herbe et les fourrages, l’aptitude laitière, l’autonomie des animaux et, pour les systèmes d’engraissement herbager, réduire le format. Enfin, pour les bovins laitiers et à viande, il est suggéré d’augmenter la capacité d’ingestion de fourrages grossiers. En regard du changement climatique, il est indispensable de prendre en compte les problèmes sanitaires associés (en particulier émergents), améliorer la tolérance au stress thermique, et réduire les émissions de méthane. La dernière partie de cet article discute des besoins des programmes de sélection induits par ces évolutions.

Among epigenetic modifications, DNA methylation is extensively studied in sperm, and variations in sperm DNA methylation patterns have been associated with fertility, sperm quality, and in vitro fertilization outcomes in several species. Reduced representation bisulfite sequencing (RRBS) is a cost-effective method that has been successfully applied to sperm. However, RRBS library preparation remains sensitive and labor-intensive and can be subjected to diverse sources of technical variation. In order to improve RRBS library preparation performance, as well as its reproducibility, we adapted the reference protocol and implemented it on a Hamilton pipetting automaton.

Background. DNA methylation contributes to the elaboration of phenotypes and is hypothesized to account for interindividual variations in farm animals. Currently, methodologies available to investigate DNA methylation in cattle rely on high throughput sequencing, which cannot be applied to large cohorts. To enable large scale DNA methylation analysis, we developed the RUMIGEN EpiChip, the first DNA methylation array specifically designed for cattle. Results. Manufactured by Illumina, the EpiChip contains 43,317 CpG markers allowing analysis of phenotypes of agronomical interest as well as the study of regulatory processes. The assay design drew on data from numerous studies achieved by the scientific community, and includes CpGs where methylation varies with health status, physiological stage, fertility, and environmental challenges, as well as CpGs located in functional genomic elements (promoters, CTCF binding sites, expression quantitative trait loci). The technical performances of the EpiChip were tested on several semen and blood DNA samples and in two laboratories, showing excellent repeatability, accuracy and interoperability. Methylation values were also concordant with those obtained by reduced representation bisulfite sequencing. The EpiChip has demonstrated a good ability to explore biological processes such as genomic imprinting and differences between cell types, opening the possibility of inferring blood cell composition. Finally, analysis of longitudinal ear biopsies allowed accurate age prediction, suggesting the potential of the array to refine epigenetic clocks. Conclusions. The RUMIGEN EpiChip is a cost-effective and versatile resource that opens new opportunities for the study of regulatory mechanisms underlying phenotypic variation and for large-scale association analyses in cattle. ### Competing Interest Statement The authors have declared no competing interest. European Union’s Horizon 2020, 101000226 Agence Nationale de la Recherche, https://ror.org/00rbzpz17, ANR-21-CE20-0021 APIS-GENE

INTRODUCTION:Although the in vitro production of bovine embryos now exceeds in vivo production, the quality of in vitro-produced embryos remains inferior. Metabolomic analysis of the spent culture medium used for embryos is considered a relevant source of markers for subsequent developmental capacity, genetic status and sex. However, little is known about the sources of variability in this metabolome and their respective significance. OBJECTIVE:We compared bovine embryo spent culture media from in vivo developed (IVD) and in vitro produced (IVP) embryos, and analyzed how they varied with embryo stage, grade and sex. METHODS:Embryos were produced in vitro under two different conditions: SOF medium supplemented with 1% fetal calf serum or IVF Bioscience media, or developed in vivo. They were recovered at Day-6 and cultured individually for 26 hours. The culture medium was analyzed using 1H-NMR spectroscopy, and the embryos were staged, graded and sexed. RESULTS:In vitro production significantly reduced inter-embryo variability of the 24 metabolites measured. The mean variance of metabolite content was 5 to 10 times lower in IVP than in IVD embryos (depending on the IVP conditions). IVD embryo stages and grades at the start or end of the culture period contributed to 8% to 18% of variability in the media content, but only 1% of variability for IVP embryos. IVD embryo stage at the start of culture significantly impacted metabolite content. IVD embryo sex did not globally impact composition of the spent culture medium, but the kinetics of decrease of about half of the metabolites differed between male and female blastocysts. CONCLUSIONS:Standardization of the metabolic properties of IVP embryos invites new questions regarding the impact of embryo biotechnology. Differences in metabolism between male and female IVD embryos are highly transient during the morula-blastocyst stages.

Bull fertility is an important economic trait, and the use of subfertile semen for artificial insemination decreases the global efficiency of the breeding sector. Although the analysis of semen functional parameters can help to identify infertile bulls, no tools are currently available to enable precise predictions and prevent the commercialization of subfertile semen. Because male fertility is a multifactorial phenotype that is dependent on genetic, epigenetic, physiological and environmental factors, we hypothesized that an integrative analysis might help to refine our knowledge and understanding of bull fertility. We combined -omics data (genotypes, sperm DNA methylation at CpGs and sperm small non-coding RNAs) and semen parameters measured on a large cohort of 98 Montbéliarde bulls with contrasting fertility levels. Multiple Factor Analysis was conducted to study the links between the datasets and fertility. Four methodologies were then considered to identify the features linked to bull fertility variation: Logistic Lasso, Random Forest, Gradient Boosting and Neural Networks. Finally, the features selected by these methods were annotated in terms of genes, to conduct functional enrichment analyses. The less relevant features in -omics data were filtered out, and MFA was run on the remaining 12,006 features, including the 11 semen parameters and a balanced proportion of each type of-omics data. The results showed that unlike the semen parameters studied the-omics datasets were related to fertility. Biomarkers related to bull fertility were selected using the four methodologies mentioned above. The most contributory CpGs, SNPs and miRNAs targeted genes were all found to be involved in development. Interestingly, fragments derived from ribosomal RNAs were overrepresented among the selected features, suggesting roles in male fertility. These markers could be used in the future to identify subfertile bulls in order to increase the global efficiency of the breeding sector.

Docosahexaenoic acid (DHA) is an n-3 polyunsaturated fatty acid (PUFA) that improves fertility by increasing membrane fluidity. Moreover, embryos produced by donor females supplied with n-3 PUFA did not show any difference in terms of the lipid profile after 7 days of culture. The present study aimed to investigate the effects of DHA (20 and 100 μM) coupled with carnosine (5 mg/mL), an antioxidant, during oocyte maturation and embryo development on the developmental and cryosurvival rates and the number of pluripotent cells. Free fatty acid receptor-4 (FFAR4), which is able to bind DHA, was visualised by immunostaining. The addition of DHA in the in vitro development (IVD) medium decreased the percentage of pluripotent SOX2 positive cells compared with the control (8.4% vs. 10.9%) without affecting the number of cells (196.7 vs. 191.6 cells) or the developmental (20.9% vs. 23.9% blastocysts rate on D7) and cryosurvival rates (86.3% vs 86.2%). Such a decrease in pluripotent cells, relevant to the differentiation of the first lineage within the inner cell mass, represents an improvement in the embryo quality. On the contrary, embryos without any pluripotent SOX2-positive cells would not be able to achieve gestation. Future studies should follow up these results by carrying out embryo transfers to assess the beneficial effects of DHA supplementation.

Paratuberculosis is a chronic infection of the intestine, mainly the ileum, caused by Mycobacterium avium subsp. paratuberculosis in cattle and other ruminants. This enzootic disease is present worldwide and has a negative impact on the dairy cattle industry. For this subspecies, the current genotyping tools do not provide the needed resolution to investigate the genetic diversity of closely related strains. These limitations can be overcome by the application of whole-genome sequencing (WGS), particularly for clonal populations such as M. avium subsp. paratuberculosis. The purpose of the present study was to undertake a WGS analysis with a panel of 200 animal field M. avium subsp. paratuberculosis strains selected based on a previous large-scale longitudinal study of Prim'Holstein and Normande dairy breeds naturally infected with M. avium subsp. paratuberculosis in the West of France. The pangenome analysis revealed that M. avium subsp. paratuberculosis has a closed pangenome. The phylogeny, based on alignment of 2,786 nonhomoplasic single nucleotide polymorphisms (SNPs), showed that the strain population is structured into three clades independently of the cattle breed or geographic distribution. The increased resolution of phylogeny obtained by WGS confirmed the homoplasic nature of the markers variable-number tandem repeat (VNTR) and short sequence repeat (SSR) used for M. avium subsp. paratuberculosis genotyping. These phylogenetic data also revealed independent introductions of the different genotypes in two main waves since at least 2003. WGS applied to this sampling demonstrated the presence of mixed infections in herds and at the individual animal level. Collectively, the phylogeny results inferred with French isolates compared to M. avium subsp. paratuberculosis isolates from around the world suggest introductions of M. avium subsp. paratuberculosis genotypes through the animal trade. Relationships between genetic traits and epidemiological data can now be investigated to better understand transmission dynamics of the disease. IMPORTANCE Mycobacterium avium subsp. paratuberculosis causes Johne's disease in ruminants, which is present worldwide and has significant negative impacts on the dairy cattle industry and animal welfare. Prevention and control of M. avium subsp. paratuberculosis infection are hampered by knowledge gaps in strain virulence, genotype distribution, and transmission dynamics. This work has revealed new insights into M. avium subsp. paratuberculosis strains currently circulating in western France and how they are related to strains circulating globally. We applied whole-genome sequencing (WGS) to obtain comprehensive information on genome evolution and discrimination of closely related strains. This approach revealed the history of M. avium subsp. paratuberculosis infection in France, refined the pangenomic characteristics of M. avium subsp. paratuberculosis, and demonstrated the existence of mixed infection in animals. Finally, this study identified predominant genotypes, which allow a better understanding of disease transmission dynamics. This information will facilitate tracking of this pathogen on farms and across agricultural regions, thus informing transmission pathways and disease control points.

Background Bovine paratuberculosis, or Johne’s disease (JD), is a contagious and incurable disease caused by Mycobacterium avium subsp. paratuberculosis (MAP). It has adverse effects on animal welfare and is very difficult to control, leading to serious economic consequences. An important line of defense to this disease is host genetic resistance to MAP, which, when it will be more fully understood, could be improved through selective breeding. Using a large dataset of Holstein cows (161,253 animals including 56,766 cows with ELISA serological phenotypes and 12,431 animals with genotypes), we applied a single-step single nucleotide polymorphism (SNP) best linear unbiased prediction approach to investigate the genetic determinism underlying resistance to this disease (heritability estimate and identification of relevant genomic regions) and estimated genetic trends, reliability, and relative risk factors associated with genomic predictions. Results Resistance to JD was moderately heritable (0.14) and 16 genomic regions were detected that accounted for at least 0.05% of the breeding values variance (GV) in resistance to JD, and were located on chromosomes 1, 3, 5, 6, 7, 19, 20, 21, 23, 25, and 27, with the highest percentage of variance explained by regions on chromosomes 23 (0.36% GV), 5 (0.22% GV), 1 (0.14% GV), and 3 (0.13% GV). When estimated for the whole chromosomes, the autosomes with the largest overall contributions were chromosomes 3 (5.3% GV), 10 (4.8%), 23 (4.7%), 1 (3.6%), 7 (3.4%), 5 (2.9%), 12 (2.5%), 11 (2.2%), and 13 (2%). We estimated a slightly favorable genetic trend in resistance to JD over the last two decades, which can be explained by a low positive genetic correlation between resistance to JD and total merit index (+ 0.06). Finally, in a validation population of 907 cows, relatively reliable genomic predictions (reliability = 0.55) were obtained, which allowed the identification of cows at high risk of infection. Conclusions This study provides new insights into the genetic determinism of resistance to JD and shows that this trait can be predicted from SNP genotypes. It has led to the implementation of a single-step genomic evaluation that should rapidly become an effective tool for controlling paratuberculosis on French Holstein farms.

Embryo lipid profile is affected by in vitro culture conditions that lead to an increase in lipids. Efforts have been made to optimize embryo lipid composition as it is associated with their quality. The objective of this study was to evaluate whether the diet supplementation of donor cows (n-3 or n-6 polyunsaturated fatty acids), or the slow freezing protocols (ethylene glycol sucrose vs. glycerol-trehalose), or the physiological stage of the donor (nulliparous heifers vs. primiparous lactating cows) may impact the bovine embryo lipid profile. Lipid extracts of 97 embryos were individually analyzed by liquid chromatography-high resolution mass spectrometry, highlighting 246 lipids, including 85% being overabundant in cow embryos compared to heifer embryos. Among 105 differential lipids, 72 were overabundant after ethylene glycol sucrose protocol, including a single glycerophosphate PA(32:1) representing 27.3% of the significantly modulated lipids, suggesting that it is degraded when glycerol-trehalose protocol is used. No lipids were different according to the n-3 or n-6 supplementation of the donor cows. In conclusion, the embryonic lipid profile was mainly affected by the physiological stage of the donors and the slow freezing protocols. The overabundance of lipids in lactating cow embryos and the resulting lower quality of these embryos are consistent with the lower pregnancy rate observed in cows compared to heifers. Unlike glycerol-trehalose protocol, ethylene glycol sucrose freezing allowed to preserve glycerophospholipids, potentially improving the slow freezing of in vitro-produced embryos. Further studies are required to modulate embryo quality and freezability by modulating the lipidome and by integrating all stages of embryonic production.

Conflicting results regarding alterations to sperm DNA methylation in cases of spermatogenesis defects, male infertility and poor developmental outcomes have been reported in humans. Bulls used for artificial insemination represent a relevant model in this field, as the broad dissemination of bull semen considerably alleviates confounding factors and enables the precise assessment of male fertility. This study was therefore designed to assess the potential for sperm DNA methylation to predict bull fertility. A unique collection of 100 sperm samples was constituted by pooling 2–5 ejaculates per bull from 100 Montbéliarde bulls of comparable ages, assessed as fertile (n = 57) or subfertile (n = 43) based on non-return rates 56 days after insemination. The DNA methylation profiles of these samples were obtained using reduced representation bisulfite sequencing. After excluding putative sequence polymorphisms, 490 fertility-related differentially methylated cytosines (DMCs) were identified, most of which were hypermethylated in subfertile bulls. Interestingly, 46 genes targeted by DMCs are involved in embryonic and fetal development, sperm function and maturation, or have been related to fertility in genome-wide association studies; five of these were further analyzed by pyrosequencing. In order to evaluate the prognostic value of fertility-related DMCs, the sperm samples were split between training (n = 67) and testing (n = 33) sets. Using a Random Forest approach, a predictive model was built from the methylation values obtained on the training set. The predictive accuracy of this model was 72% on the testing set and 72% on individual ejaculates collected from an independent cohort of 20 bulls. This study, conducted on the largest set of bull sperm samples so far examined in epigenetic analyses, demonstrated that the sperm methylome is a valuable source of male fertility biomarkers. The next challenge is to combine these results with other data on the same sperm samples in order to improve the quality of the model and better understand the interplay between DNA methylation and other molecular features in the regulation of fertility. This research may have potential applications in human medicine, where infertility affects the interaction between a male and a female, thus making it difficult to isolate the male factor.

Semen infertility or sub-fertility, whether in humans or livestock species, remains a major concern for clinicians and technicians involved in reproduction. Indeed, they can cause tragedies in human relationships or have a dramatic overall negative impact on the sustainability of livestock breeding. Understanding and predicting semen fertility issues is therefore crucial and quality control procedures as well as biomarkers have been proposed to ensure sperm fertility. However, their predictive values appeared to be too limited and additional relevant biomarkers are still required to diagnose sub-fertility efficiently. During the last decade, the study of molecular mechanisms involved in spermatogenesis and sperm maturation highlighted the regulatory role of a variety of small non-coding RNAs (sncRNAs) and led to the discovery that sperm sncRNAs comprise both remnants from spermatogenesis and post-testicular sncRNAs acquired through interactions with extracellular vesicles along epididymis. This has led to the hypothesis that sncRNAs may be a source of relevant biomarkers, associated either with sperm functionality or embryo development. This review aims at providing a synthetic overview of the current state of knowledge regarding implication of sncRNA in spermatogenesis defects and their putative roles in sperm maturation and embryo development, as well as exploring their use as fertility biomarkers.

A gain-of-function mutation in the fibroblast growth factor receptor 3 gene (FGFR3) results in achondroplasia (ACH), the most frequent form of dwarfism. Constitutive activation of FGFR3 impairs bone formation and elongation and many signal transduction pathways. Identification of new and relevant compounds targeting the FGFR3 signaling pathway is of broad importance for the treatment of ACH, and natural plant compounds are prime drug candidate sources. Here, we found that the phenolic compound (-)-epicatechin, isolated from Theobroma cacao, effectively inhibited FGFR3's downstream signaling pathways. Transcriptomic analysis in an Fgfr3 mouse model showed that ciliary mRNA expression was modified and influenced significantly by the Indian hedgehog and PKA pathways. (-)-Epicatechin is able to rescue mRNA expression impairments that control both the structural organization of the primary cilium and ciliogenesis-related genes. In femurs isolated from a mouse model (Fgfr3Y367C/+) of ACH, we showed that (-)-epicatechin eliminated bone growth impairment during 6 days of ex vivo culture. In vivo, we confirmed that daily subcutaneous injections of (-)-epicatechin to Fgfr3Y367C/+ mice increased bone elongation and rescued the primary cilium defects observed in chondrocytes. This modification to the primary cilia promoted the typical columnar arrangement of flat proliferative chondrocytes and thus enhanced bone elongation. The results of the present proof-of-principle study support (-)-epicatechin as a potential drug for the treatment of ACH.

Despite passing stringent quality control, bulls used in artificial insemination can vary significantly in their fertility, emphasizing the need for reliable markers of sperm quality. This study aimed to identify sperm proteins acting as biomarkers of fertility in 2 different populations of dairy bulls classified based on their field fertility. Semen was collected and cryopreserved from: 54 Holstein bulls located in Ireland, classified according to fertility indexes as low fertility (LF, n = 23), medium fertility (n = 14), or high fertility (HF, n = 17); and 18 Holstein bulls located in Denmark, classified as LF (n = 8) or HF (n = 10). The proteome was measured through liquid chromatography-mass spectrometry and data were analyzed with the R software. Differentially abundant proteins between HF and LF bulls and biomarker proteins were determined through a modified t-test and random forest, respectively, selecting 301 differentially abundant proteins and 34 biomarker proteins. The predictive ability of the 34 biomarkers was evaluated employing support vector machine as the classifier, using their abundance levels in the Irish bulls to train the model and in the Danish bulls for validation. The prediction accuracy was 94.4%, with only one HF bull misclassified, corresponding to the lowest fertility index bull in the HF group. The biomarkers more abundant in sperm of HF bulls enriched axoneme assembly and sperm motility (false discovery rate <0.05), according to functional analysis. In conclusion, a robust model coupled with the application of appropriate bioinformatic tools allowed the identification of functionally relevant sperm proteins predictive of the fertility of Holstein bulls used in artificial insemination.

Bovine paratuberculosis, or Johne’s disease (JD), is a contagious, incurable and very difficult to control disease due to Mycobacterium avium subsp. paratuberculosis (MAP). Using a single step GBLUP, including candidate SNPs for resistance to JD, applied to a pedigree of 161,253 Holstein animals (56,766 with serological phenotypes and 12,431 with genotypes), we estimated a moderate heritability value (0.14) of this trait and detected genomic regions associated to resistance to JD on chromosomes 23, 21, 12, 3, 7, 27, 5, 10, and 20. We also estimated a slightly favorable genetic trend of resistance to JD over the last two decades and relatively reliable genomic predictions (CD=0.55) in a validation population allowing the identification of cows at high risk of infection. Genomic predictions should therefore rapidly become an effective tool for controlling paratuberculosis in French Holstein cow farms.

https://doi.org/10.1093/af/vfab062 This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. © Kiefer, Sellem, Bonnet-Garnier, Pannetier, Costes, Schibler, Jammes Feature Article The epigenome of male germ cells and the programming of phenotypes in cattle

Currently, in vitro embryo production (IVP) is successfully commercially applied in cattle. However, the high sensitivity of embryos to cryopreservation in comparison to in vivo (IVD) embryos slows the dissemination of this biotechnology. Reduced cryotolerance is frequently associated with lipid accumulation in the cytoplasm mainly due to in vitro culture conditions. The objective of this study was to evaluate the lipid composition of biopsied and sexed embryos, produced either in vivo or in vitro from the same Holstein heifers before and after a slow freezing protocol. Lipid extracts were analysed by liquid chromatography-high resolution mass spectrometry, which enabled the detection of 496 features. Our results highlighted a lipid enrichment of IVP embryos in triglycerides and oxidised glycerophospholipids and a reduced abundance in glycerophospholipids. The slow freezing process affected the lipid profiles of IVP and IVD embryos similarly. Lysophosphatidylcholine content was reduced when embryos were frozen/thawed. In conclusion, the embryonic lipid profile is impacted by IVP and slow freezing protocols but not by sex. Lysophosphatidylcholine seemed highly sensitive to cryopreservation and might contribute to explain the lower quality of frozen embryos. Further studies are required to improve embryo freezability by modulating the lipidome.