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.

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.

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.

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.

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.

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.

Analytical Methodology for a Metabolome Atlas of Goat's Plasma, Milk and Feces Using 1H-NMR and UHPLC-HRMS.

Metabolomics has been increasingly used in animal and food sciences. Animal health is one of the most important factor that can also alter animal integrity and welfare. Some studies have already investigated the link between health and metabolic profile of dairy animals. These studies in metabolomics often consider a single type of sample using a single analytical platform (nuclear magnetic resonance or mass spectrometry). Only few studies with multi-platform approaches are also used with a single or a multi type of sample, but they mainly consider dairy cows' metabolome although dairy goats present similar diseases, that it could be interesting to detect early to preserve animal health and milk production. This study aims to create a metabolic atlas of goat plasma, milk and feces, based on healthy animals. Our study describes a standard operating procedure for three goat matrices: blood plasma, milk, and feces using multiple platforms (NMR (1H), UHPLC (RP)-MS and UHPLC (HILIC)-MS) that follows a unique sample preparation procedure for each sample type to be analyzed on multi-platforms basis. Our method was evaluated for its robustness and allowed a better characterization of goat metabolic profile in healthy conditions.

Draft Genome Sequences of 142 Mycobacterium avium subsp. paratuberculosis Strains Isolated from Naturally Infected Dairy Cattle.

Mycobacterium avium subsp. paratuberculosis is the etiological agent of Johne's disease in ruminants. Here, we report the annotated draft genome sequences of 142 M. avium subsp. paratuberculosis strains that were isolated from dairy cattle in France between 2014 and 2018. The genomes of these strains were sequenced using Illumina technology.

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.

Lipid profile of bovine grade-1 blastocysts produced either in vivo or in vitro before and after slow freezing process.

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.

Dynamics of cattle sperm sncRNAs during maturation, from testis to ejaculated sperm.

BACKGROUND: During epididymal transit, spermatozoa go through several functional maturation steps, resulting from interactions with epididymal secretomes specific to each region. In particular, the sperm membrane is under constant remodeling, with sequential attachment and shedding of various molecules provided by the epididymal lumen fluid and epididymosomes, which also deliver sncRNA cargo to sperm. As a result, the payload of sperm sncRNAs changes during the transit from the epididymis caput to the cauda. This work was designed to study the dynamics of cattle sperm sncRNAs from spermatogenesis to final maturation. RESULTS: Comprehensive catalogues of sperm sncRNAs were obtained from testicular parenchyma, epididymal caput, corpus and cauda, as well as ejaculated semen from three Holstein bulls. The primary cattle sncRNA sperm content is markedly remodeled as sperm mature along the epididymis. Expression of piRNAs, which are abundant in testis parenchyma, decreases dramatically at epididymis. Conversely, sperm progressively acquires miRNAs, rsRNAs, and tsRNAs along epididymis, with regional specificities. For instance, miRNAs and tsRNAs are enriched in epididymis cauda and ejaculated sperm, while rsRNA expression peaks at epididymis corpus. In addition, epididymis corpus contains mainly 20 nt long piRNAs, instead of 30 nt in all other locations. Beyond the bulk differences in abundance of sncRNAs classes, K-means clustering was performed to study their spatiotemporal expression profile, highlighting differences in specific sncRNAs and providing insights into their putative biological role at each maturation stage. For instance, Gene Ontology analyses using miRNA targets highlighted enriched processes such as cell cycle regulation, response to stress and ubiquitination processes in testicular parenchyma, protein metabolism in epididymal sperm, and embryonic morphogenesis in ejaculated sperm. CONCLUSIONS: Our findings confirm that the sperm sncRNAome does not simply reflect a legacy of spermatogenesis. Instead, sperm sncRNA expression shows a remarkable level of plasticity resulting probably from the combination of multiple factors such as loss of the cytoplasmic droplet, interaction with epididymosomes, and more surprisingly, the putative in situ production and/or modification of sncRNAs by sperm. Given the suggested role of sncRNA in epigenetic trans-generational inheritance, our detailed spatiotemporal analysis may pave the way for a study of sperm sncRNAs role in embryo development.

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.

Accelerating Onset of Puberty Through Modification of Early Life Nutrition Induces Modest but Persistent Changes in Bull Sperm DNA Methylation Profiles Post-puberty.

In humans and model species, alterations of sperm DNA methylation patterns have been reported in cases of spermatogenesis defects, male infertility and exposure to toxins or nutritional challenges, suggesting that a memory of environmental or physiological changes is recorded in the sperm methylome. The objective of this study was to ascertain if early life plane of nutrition could have a latent effect on DNA methylation patterns in sperm produced post-puberty. Holstein-Friesian calves were assigned to either a high (H) or moderate (M) plane of nutrition for the first 24 weeks of age, then reassigned to the M diet until puberty, resulting in HM and MM groups. Sperm DNA methylation patterns from contrasted subgroups of bulls in the HM (ejaculates recovered at 15 months of age; n = 9) and in the MM (15 and 16 months of age; n = 7 and 9, respectively) were obtained using Reduced Representation Bisulfite Sequencing. Both 15 and 16 months were selected in the MM treatment as these bulls reached puberty approximately 1 month after the HM bulls. Hierarchical clustering demonstrated that inter-individual variability unrelated to diet or age dominated DNA methylation profiles. While the comparison between 15 and 16 months of age revealed almost no change, 580 differentially methylated CpGs (DMCs) were identified between the HM and MM groups. Differentially methylated CpGs were mostly hypermethylated in the HM group, and enriched in endogenous retrotransposons, introns, intergenic regions, and shores and shelves of CpG islands. Furthermore, genes involved in spermatogenesis, Sertoli cell function, and the hypothalamic-pituitary-gonadal axis were targeted by differential methylation when HM and MM groups were compared at 15 months of age, reflecting the earlier timing of puberty onset in the HM bulls. In contrast, the genes still differentially methylated in MM bulls at 16 months of age were enriched for ATP-binding molecular function, suggesting that changes to the sperm methylome could persist even after the HM and MM bulls reached a similar level of sexual maturity. Together, results demonstrate that enhanced plane of nutrition in pre-pubertal calves associated with advanced puberty induced modest but persistent changes in sperm DNA methylation profiles after puberty.

A comprehensive overview of bull sperm-borne small non-coding RNAs and their diversity across breeds.

BACKGROUND: Mature sperm carry thousands of RNAs, including mRNAs, lncRNAs, tRNAs, rRNAs and sncRNAs, though their functional significance is still a matter of debate. Growing evidence suggests that sperm RNAs, especially sncRNAs, are selectively retained during spermiogenesis or specifically transferred during epididymis maturation, and are thus delivered to the oocyte at fertilization, providing resources for embryo development. However , a deep characterization of the sncRNA content of bull sperm and its expression profile across breeds is currently lacking. To fill this gap, we optimized a guanidinium-Trizol total RNA extraction protocol to prepare high-quality RNA from frozen bull sperm collected from 40 representative bulls from six breeds. Deep sequencing was performed (40 M single 50-bp reads per sample) to establish a comprehensive repertoire of cattle sperm sncRNA. RESULTS: Our study showed that it comprises mostly piRNAs (26%), rRNA fragments (25%), miRNAs (20%) and tRNA fragments (tsRNA, 14%). We identified 5p-halves as the predominant tsRNA subgroup in bull sperm, originating mostly from Gly and Glu isoacceptors. Our study also increased by ~ 50% the sperm repertoire of known miRNAs and identified 2022 predicted miRNAs. About 20% of sperm miRNAs were located within genomic clusters, expanding the list of known polycistronic pri-miRNA clusters and defining several networks of co-expressed miRNAs. Strikingly, our study highlighted the great diversity of isomiRs, resulting mainly from deletions and non-templated additions (A and U) at the 3p end. Substitutions within miRNA sequence accounted for 40% of isomiRs, with G>A, U>C and C>U substitutions being the most frequent variations. In addition, many sncRNAs were found to be differentially expressed across breeds. CONCLUSIONS: Our study provides a comprehensive overview of cattle sperm sncRNA, and these findings will pave the way for future work on the role of sncRNAs in embryo development and their relevance as biomarkers of semen fertility.

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.