Effects of the donor factors and freezing protocols on the bovine embryonic lipid profile†.

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.

New insights into the genetic resistance to paratuberculosis in Holstein cattle via single-step genomic evaluation.

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.

Sperm-borne sncRNAs: potential biomarkers for semen fertility?

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.

Optimization of Sample Preparation for Metabolomics Exploration of Urine, Feces, Blood and Saliva in Humans Using Combined NMR and UHPLC-HRMS Platforms.

Currently, most clinical studies in metabolomics only consider a single type of sample such as urine, plasma, or feces and use a single analytical platform, either NMR or MS. Although some studies have already investigated metabolomics data from multiple fluids, the information is limited to a unique analytical platform. On the other hand, clinical studies investigating the human metabolome that combine multi-analytical platforms have focused on a single biofluid. Combining data from multiple sample types for one patient using a multimodal analytical approach (NMR and MS) should extend the metabolome coverage. Pre-analytical and analytical phases are time consuming. These steps need to be improved in order to move into clinical studies that deal with a large number of patient samples. Our study describes a standard operating procedure for biological specimens (urine, blood, saliva, and feces) using multiple platforms (1H-NMR, RP-UHPLC-MS, and HILIC-UHPLC-MS). Each sample type follows a unique sample preparation procedure for analysis on a multi-platform basis. Our method was evaluated for its robustness and was able to generate a representative metabolic map.

Identification of the ABCC4, IER3, and CBFA2T2 candidate genes for resistance to paratuberculosis from sequence-based GWAS in Holstein and Normande dairy cattle.

BACKGROUND: Bovine paratuberculosis is a contagious disease, caused by Mycobacterium avium subsp. paratuberculosis (MAP), with adverse effects on animal welfare and serious economic consequences. Published results on host genetic resistance to MAP are inconsistent, mainly because of difficulties in characterizing the infection status of cows. The objectives of this study were to identify quantitative trait loci (QTL) for resistance to MAP in Holstein and Normande cows with an accurately defined status for MAP. RESULTS: From MAP-infected herds, cows without clinical signs of disease were subjected to at least four repeated serum ELISA and fecal PCR tests over time to determine both infected and non-infected statuses. Clinical cases were confirmed using PCR. Only cows that had concordant results for all tests were included in further analyses. Positive and control cows were matched within herd according to their birth date to ensure a same level of exposure to MAP. Cows with accurate phenotypes, i.e. unaffected (control) or affected (clinical or non-clinical cases), were genotyped with the Illumina BovineSNP50 BeadChip. Genotypes were imputed to whole-genome sequences using the 1000 Bull Genomes reference population (run6). A genome-wide association study (GWAS) of MAP status of 1644 Holstein and 649 Normande cows, using either two (controls versus cases) or three classes of phenotype (controls, non-clinical and clinical cases), revealed three regions, on Bos taurus (BTA) chromosomes 12, 13, and 23, presenting significant effects in Holstein cows, while only one of those was identified in Normande cows (BTA23). The most significant effect was found on BTA13, in a short 8.5-kb region. Conditional analyses revealed that only one causal variant may be responsible for the effects observed on each chromosome with the ABCC4 (BTA12), CBFA2T2 (BTA13), and IER3 (BTA23) genes as good functional candidates. CONCLUSIONS: A sequence-based GWAS on cows for which resistance to MAP was accurately defined, was able to identify candidate variants located in genes that were functionally related to resistance to MAP; these explained up to 28% of the genetic variance of the trait. These results are very encouraging for efforts towards implementation of a breeding strategy aimed at improving resistance to paratuberculosis in Holstein cows.

Longitudinal study of Mycobacterium avium ssp. paratuberculosis fecal shedding patterns and concurrent serological patterns in naturally infected dairy cattle.

Mycobacterium avium ssp. paratuberculosis (MAP) is the etiological agent of paratuberculosis, a disease that affects ruminants worldwide. Despite global interest in the control of this disease, gaps exist in our knowledge of fecal shedding patterns and concurrent serological patterns. This longitudinal study in dairy cattle herds with high MAP seroprevalence in France aimed at accurately describing fecal shedding patterns over 1 year; relating those shedding patterns to individual animal characteristics (age, breed, parity); and exploring the association between fecal shedding patterns and serological patterns. To describe temporal fecal shedding patterns and continuity of shedding, along with the standard quantitative PCR (qPCR) threshold cycle we used a cutoff value that related to low or nonculturable fecal shedding. We also defined a threshold cycle indicative of shedding in high quantities to describe infection progression patterns. Twenty-one herds completed the study, and 782 cows were tested 4 times each. We obtained 4 sets of paired fecal qPCR and serum ELISA results from 757 cows. Although we targeted highly likely infectious animals, we found a large diversity of shedding patterns, as well as high variability between herds in the proportion of animals showing a given pattern. The fecal qPCR results of almost 20% of the final study sample were positioned at least once in the range that indicated low or nonculturable fecal shedding (between the adjusted and the standard cutoff value). Although these animals would typically be classified as non-shedders, they could be important to infection dynamics on the farm. Animals that shed at least twice consecutively and animals that shed in high quantities rarely reverted to negativity. Repeated fecal qPCR can be used to detect temporal fecal shedding traits, and the decision to cull an animal could practically be based on temporal, semiquantitative results. Overall, we found a mismatch between fecal shedding and ELISA seropositivity (637 animals were ELISA-negative 4 times, but only 13% of those animals were qPCR-negative 4 times). We found that having more than 2 ELISA-positive samples was strongly related to persistent and continuous shedding. We suggest that although serological testing is much less sensitive than qPCR, it can also be used, particularly over the course of multiple testing events, to identify animals that are most likely to contribute to the contamination of the farm environment.

A multi-scale analysis of bull sperm methylome revealed both species peculiarities and conserved tissue-specific features.

BACKGROUND: Spermatozoa have a remarkable epigenome in line with their degree of specialization, their unique nature and different requirements for successful fertilization. Accordingly, perturbations in the establishment of DNA methylation patterns during male germ cell differentiation have been associated with infertility in several species. While bull semen is widely used in artificial insemination, the literature describing DNA methylation in bull spermatozoa is still scarce. The purpose of this study was therefore to characterize the bull sperm methylome relative to both bovine somatic cells and the sperm of other mammals through a multiscale analysis. RESULTS: The quantification of DNA methylation at CCGG sites using luminometric methylation assay (LUMA) highlighted the undermethylation of bull sperm compared to the sperm of rams, stallions, mice, goats and men. Total blood cells displayed a similarly high level of methylation in bulls and rams, suggesting that undermethylation of the bovine genome was specific to sperm. Annotation of CCGG sites in different species revealed no striking bias in the distribution of genome features targeted by LUMA that could explain undermethylation of bull sperm. To map DNA methylation at a genome-wide scale, bull sperm was compared with bovine liver, fibroblasts and monocytes using reduced representation bisulfite sequencing (RRBS) and immunoprecipitation of methylated DNA followed by microarray hybridization (MeDIP-chip). These two methods exhibited differences in terms of genome coverage, and consistently, two independent sets of sequences differentially methylated in sperm and somatic cells were identified for RRBS and MeDIP-chip. Remarkably, in the two sets most of the differentially methylated sequences were hypomethylated in sperm. In agreement with previous studies in other species, the sequences that were specifically hypomethylated in bull sperm targeted processes relevant to the germline differentiation program (piRNA metabolism, meiosis, spermatogenesis) and sperm functions (cell adhesion, fertilization), as well as satellites and rDNA repeats. CONCLUSIONS: These results highlight the undermethylation of bull spermatozoa when compared with both bovine somatic cells and the sperm of other mammals, and raise questions regarding the dynamics of DNA methylation in bovine male germline. Whether sperm undermethylation has potential interactions with structural variation in the cattle genome may deserve further attention.

A missense variant in the coil1A domain of the keratin 25 gene is associated with the dominant curly hair coat trait (Crd) in horse.

BACKGROUND: Curly horses present a variety of curl phenotypes that are associated with various degrees of curliness of coat, mane, tail and ear hairs. Their origin is still a matter of debate and several genetic hypotheses have been formulated to explain the diversity in phenotype, including the combination of autosomal dominant and recessive alleles. Our purpose was to map the autosomal dominant curly hair locus and identify the causal variant using genome-wide association study (GWAS) and whole-genome sequencing approaches. RESULTS: A GWAS was performed using a Bayesian sparse linear mixed model, based on 51 curly and 19 straight-haired French and North American horses from 13 paternal families genotyped on the Illumina EquineSNP50 BeadChip. A single strong signal was observed on equine chromosome 11, in a region that encompasses the type I keratin gene cluster. This region was refined by haplotype analysis to a segment including 36 genes, among which are 10 keratin genes (KRT-10, -12, -20, -23, -24, -25, -26, -27, -28, -222). To comprehensively identify candidate causal variants within all these genes, whole-genome sequences were obtained for one heterozygous curly stallion and its straight-haired son. Among the four non-synonymous candidate variants identified and validated in the curly region, only variant g.21891160G>A in the KRT25 gene (KRT25:p.R89H) was in perfect agreement with haplotype status in the whole pedigree. Genetic association was then confirmed by genotyping a larger population consisting of 353 horses. However, five discordant curly horses were observed, which carried neither the variant nor the main haplotype associated with curliness. Sequencing of KRT25 for two discordant horses did not identify any other deleterious variant, which suggests locus rather than allelic heterogeneity for the curly phenotype. CONCLUSIONS: We identified the KRT25:p.R89H variant as responsible for the dominant curly trait, but a second dominant locus may also be involved in the shape of hairs within North American Curly horses.

A missense mutation in PFAS (phosphoribosylformylglycinamidine synthase) is likely causal for embryonic lethality associated with the MH1 haplotype in Montbéliarde dairy cattle.

A candidate mutation in the sex hormone binding globulin gene was proposed in 2013 to be responsible for the MH1 recessive embryonic lethal locus segregating in the Montbéliarde breed. In this follow-up study, we excluded this candidate variant because healthy homozygous carriers were observed in large-scale genotyping data generated in the framework of the genomic selection program. We fine mapped the MH1 locus in a 702-kb interval and analyzed genome sequence data from the 1,000 bull genomes project and 54 Montbéliarde bulls (including 14 carriers and 40 noncarriers). We report the identification of a strong candidate mutation in the gene encoding phosphoribosylformylglycinamidine synthase (PFAS), a protein involved in de novo purine synthesis. This mutation, located in a class I glutamine amidotransferase-like domain, results in the substitution of an arginine residue that is entirely conserved among eukaryotes by a cysteine (p.R1205C). No homozygote for the cysteine-encoding allele was observed in a large population of more than 25,000 individuals despite a 6.7% allelic frequency and 122 expected homozygotes under neutrality assumption. Genotyping of 18 embryos collected from heterozygous parents as well as analysis on nonreturn rates suggested that most homozygous carriers died between 7 and 35 d postinsemination. The identification of this strong candidate mutation will enable the accurate testing of the reproducers and the efficient selection against this lethal recessive embryonic defect in the Montbéliarde breed.

Novel equine tissue miRNAs and breed-related miRNA expressed in serum.

BACKGROUND: MiRNAs regulate multiple genes at the post-transcriptional level and therefore play an important role in many biological processes. It has been suggested that miRNA exported outside the cells contribute to inter-cellular communication. Consequently, circulating miRNAs are of particular interest and are promising biomarkers for many diseases. The number of miRNAs annotated in the horse genome is much lower compared to model organisms like human and mouse. We therefore aimed to identify novel equine miRNAs for tissue types and breed in serum. RESULTS: We analysed 71 small RNA-seq libraries derived from nine tissues (gluteus medius, platysma, masseter muscle, heart, liver, cartilage, bone, total blood and serum) using miRDeep2 and miRdentify tools. Known miRNAs represented between 2.3 and 62.9 % of the reads in 71 libraries. A total of 683 novel miRNAs were identified. Breed and tissue type affected the number of miRNAs detected and interestingly, affected its average intensity. A total of 50 miRNAs in serum proved to be potential biomarkers to differentiate specific breed types, of which miR-122, miR-200, miR-483 were over-expressed and miR-328 was under-expressed in ponies compared to Warmbloods. The different miRNAs profiles, as well as the differences in their expression levels provide a foundation for more hypotheses based on the novel miRNAs discovered. CONCLUSIONS: We identified 683 novel equine miRNAs expressed in seven solid tissues, blood and serum. Additionally, our approach evidenced that such data supported identification of specific miRNAs as markers of functions related to breeds or disease tissues.

Next-generation sequencing identifies equine cartilage and subchondral bone miRNAs and suggests their involvement in osteochondrosis physiopathology.

BACKGROUND: MicroRNAs (miRNAs) are an abundant class of small single-stranded non-coding RNA molecules ranging from 18 to 24 nucleotides. They negatively regulate gene expression at the post-transcriptional level and play key roles in many biological processes, including skeletal development and cartilage maturation. In addition, miRNAs involvement in osteoarticular diseases has been proved and some of them were identified as suitable biomarkers for pathological conditions. Equine osteochondrosis (OC) is one of the most prevalent juvenile osteoarticular disorders in horses and represents a major concern for animal welfare and economic reasons. Its etiology and pathology remain controversial and biological pathways as well as molecular mechanisms involved in the physiopathology are still unclear. This study aims to investigate the potential role of miRNAs in equine osteochondrosis (OC) physiopathology.Short-read NGS technology (SOLID™, Life Technologies) was used to establish a comprehensive repertoire of miRNA expressed in either equine cartilage or subchondral bone. Undamaged cartilage and subchondral bone samples from healthy (healthy samples) and OC-affected (predisposed samples) 10-month Anglo-Arabian foals were analysed. Samples were also subjected or not to an experimental mechanical loading to evaluate the role of miRNAs in the regulation of mechano-transduction pathways. Predicted targets of annotated miRNAs were identified using miRmap. RESULTS: Epiphyseal cartilage and subchondral bone miRNome were defined, including about 300 new miRNAs. Differentially expressed miRNAs were identified between bone and cartilage from healthy and OC foals, as well as after an experimental mechanical loading. In cartilage, functional annotation of their predicted targets suggests a role in the maintenance of cartilage integrity through the control of cell cycle and differentiation, energy production and metabolism as well as extracellular matrix structure and dynamics. In bone, miRNA predicited targets were associated with osteoblasts and osteoclasts differentiation, though the regulation of energy production, vesicle transport and some growth factor signaling pathways. CONCLUSION: Taken together, our results suggest a role of miRNAs in equine OC physiopathology and in the cellular response to biomechanical stress in cartilage and bone. In silico target prediction and functional enrichment analysis provides new insight into OC molecular physiopathology.

Omics technologies provide new insights into the molecular physiopathology of equine osteochondrosis.

BACKGROUND: Osteochondrosis (OC(D)) is a juvenile osteo-articular disorder affecting several mammalian species. In horses, OC(D) is considered as a multifactorial disease and has been described as a focal disruption of endochondral ossification leading to the development of osteoarticular lesions. Nevertheless, OC(D) physiopathology is poorly understood. Affected horses may present joint swelling, stiffness and lameness. Thus, OC(D) is a major concern for the equine industry. Our study was designed as an integrative approach using omics technologies for the identification of constitutive defects in epiphyseal cartilage and/or subchondral bone associated with the development of primary lesions to further understand OC(D) pathology. This study compared samples from non-affected joints (hence lesion-free) from OC(D)-affected foals (n = 5, considered predisposed samples) with samples from OC-free foals (n = 5) considered as control samples. Consequently, results are not confounded by changes associated with the evolution of the lesion, but focus on altered constitutive molecular mechanisms. Comparative proteomics and micro computed tomography analyses were performed on predisposed and OC-free bone and cartilage samples. Metabolomics was also performed on synovial fluid from OC-free, OC(D)-affected and predisposed joints. RESULTS: Two lesion subtypes were identified: OCD (lesion with fragment) and OC (osteochondral defects). Modulated proteins were identified using omics technologies (2-DE proteomics) in cartilage and bone from affected foals compare to OC-free foals. These were associated with cellular processes including cell cycle, energy production, cell signaling and adhesion as well as tissue-specific processes such as chondrocyte maturation, extracellular matrix and mineral metabolism. Of these, five had already been identified in synovial fluid of OC-affected foals: ACTG1 (actin, gamma 1), albumin, haptoglobin, FBG (fibrinogen beta chain) and C4BPA (complement component 4 binding protein, alpha). CONCLUSION: This study suggests that OCD lesions may result from a cartilage defect whereas OC lesions may be triggered by both bone and cartilage defects, suggesting that different molecular mechanisms responsible for the equine osteochondrosis lesion subtypes and predisposition could be due to a defect in both bone and cartilage. This study will contribute to refining the definition of OC(D) lesions and may improve diagnosis and development of therapies for horses and other species, including humans.

A novel tyrosine kinase inhibitor restores chondrocyte differentiation and promotes bone growth in a gain-of-function Fgfr3 mouse model.

Activating germline fibroblast growth factor receptor 3 (FGFR3) mutations cause achondroplasia (ACH), the most common form of human dwarfism and a spectrum of skeletal dysplasias. FGFR3 is a tyrosine kinase receptor and constitutive FGFR3 activation impairs endochondral ossification and triggers severe disorganization of the cartilage with shortening of long bones. To decipher the role of FGFR3 in endochondral ossification, we analyzed the impact of a novel tyrosine kinase inhibitor (TKI), A31, on both human and mouse mutant FGFR3-expressing cells and on the skeleton of Fgfr3(Y367C/+) dwarf mice. We found that A31 inhibited constitutive FGFR3 phosphorylation and restored the size of embryonic dwarf femurs using an ex vivo culture system. The increase in length of the treated mutant femurs was 2.6 times more than for the wild-type. Premature cell cycle exit and defective chondrocyte differentiation were observed in the Fgfr3(Y367C/+) growth plate. A31 restored normal expression of cell cycle regulators (proliferating cell nuclear antigen, KI67, cyclin D1 and p57) and allowed pre-hypertrophic chondrocytes to properly differentiate into hypertrophic chondocytes. Our data reveal a specific role for FGFR3 in the cell cycle and chondrocyte differentiation and support the development of TKIs for the treatment of FGFR3-related chondrodysplasias.

Involvement of mitochondrial dysfunction and ER-stress in the physiopathology of equine osteochondritis dissecans (OCD).

Osteochondrosis (OC) is a developmental bone disorder affecting several mammalian species including the horse. Equine OC is described as a focal disruption of endochondral ossification, leading to osteochondral lesions (osteochondritis dissecans, OCD) that may release free bodies within the joint. OCD lesions trigger joint swelling, stiffness and lameness and affects about 30% of the equine population. OCD is considered as multifactorial but its physiopathology is still poorly understood and genes involved in genetic predisposition are still unknown. Our study compared two healthy and two OC-affected 18-month-old French Trotters diagnosed with OCD lesions at the intermediate ridge of the distal tibia. A comparative shot-gun proteomic analysis of non-wounded cartilage and sub-chondral bone from healthy (healthy samples) and OC-affected foals (predisposed samples) identified 83 and 53 modulated proteins, respectively. These proteins are involved in various biological pathways including matrix structure and maintenance, protein biosynthesis, folding and transport, mitochondrial activity, energy and calcium metabolism. Transmission electron microscopy revealed typical features of mitochondrial swelling and ER-stress, such as large, empty mitochondria, and hyper-dilated rough endoplasmic reticulum, in the deep zone of both OC lesions and predisposed cartilage. Abnormal fibril organization surrounding chondrocytes and abnormal features at the ossification front were also observed. Combining these findings with quantitative trait loci and whole genome sequencing results identified about 140 functional candidate genes carrying putative damaging mutations in 30 QTL regions. In summary, our study suggests that OCD lesions may result from defective hypertrophic terminal differentiation associated with mitochondrial dysfunction and ER-stress, leading to impaired cartilage and bone biomechanical properties, making them prone to fractures. In addition, 11 modulated proteins and several candidate mutations located in QTL regions were identified, bringing new insight into the molecular physiopathology and genetic basis of OCD.

Multi-omics data integration for the identification of biomarkers for bull fertility.

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.

A method for proteomic analysis of equine subchondral bone and epiphyseal cartilage.

Proteomic analyses of cartilage and, to a lesser extent, of bone have long been impaired because of technical challenges related to their structure and biochemical properties. We have developed a unified method based on phenol extraction, 2DE, silver staining, and subsequent LC-MS/MS. This method proved to be efficient to characterize the proteome of equine cartilage and bone samples collected in vivo. Since proteins from several cellular compartments could be recovered, our procedure is mainly suitable for in situ molecular physiology studies focused on the cellular content of chondrocytes, osteoblasts, and osteoclasts as well as that of the extracellular matrix, with the exception of proteoglycans. Our method alleviates some drawbacks of cell culture that can mask physiological differences, as well as reduced reproducibility due to fractionation. Proteomic comparative studies between cartilage and bone samples from healthy and affected animals were thus achieved successfully. This achievement will contribute to increasing knowledge on the molecular mechanisms underlying the physiopathology of numerous osteoarticular diseases in horses and in humans.

Effect of DHA on the quality of In vitro produced bovine embryos.

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.

Genetic Features of Mycobacterium avium subsp. paratuberculosis Strains Circulating in the West of France Deciphered by Whole-Genome Sequencing.

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.

Predicting male fertility from the sperm methylome: application to 120 bulls with hundreds of artificial insemination records.

BACKGROUND: 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. RESULTS: 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. CONCLUSION: 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.

Theobroma cacao improves bone growth by modulating defective ciliogenesis in a mouse model of achondroplasia.

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.