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.

Combining accelerometers and direct visual observations to detect sickness and pain in cows of different ages submitted to systemic inflammation.

Cattle suffering from inflammatory infection display sickness and pain-related behaviours. As these behaviours may be transient and last only a few hours, one may miss them. The aim of this study was to assess the benefit of combining continuous monitoring of cow behaviour via collar-attached accelerometers with direct visual observations to detect sickness and pain-related behavioural responses after a systemic inflammatory challenge (intravenous lipopolysaccharide injection) in cows of two different ages, proven by clinical, physiological and blood parameters. Twelve cloned Holstein cows (six 'old' cows aged 10-15 years old and six 'young' cows aged 6 years old) were challenged and either directly observed at five time-points from just before the lipopolysaccharide injection up to 24 h post-injection (hpi) or continuously monitored using collar-attached accelerometers in either control or challenge situations. Direct observations identified specific sickness and pain behaviours (apathy, changes in facial expression and body posture, reduced motivation to feed) expressed partially at 3 hpi and fully at 6 hpi. These signs of sickness and pain behaviours then faded, and quicker for the young cows. Accelerometers detected changes in basic activities (low ingesting, low ruminating, high inactivity) and position (high time standing up) earlier and over a longer period of time than direct observations. The combination of sensors and direct observations improved the detection of behavioural signs of sickness and pain earlier on and over the whole study period, even when direct signs were weak especially in young cows. This system could provide great benefit for better earlier animal care.

Prepubertal nutritional modulation in the bull and its impact on sperm DNA methylation.

Enhanced pre-pubertal nutrition in Holstein bulls increased reproductive hormone production and sperm production potential with no negative effects on sperm quality. However, recent trends in human epigenetic research have identified pre-pubertal period to be critical for epigenetic reprogramming in males. Our objective was to evaluate the methylation changes in sperm of bulls exposed to different pre-pubertal diets. One-week-old Holstein bull calves (n = 9), randomly allocated to 3 groups, were fed either a high, medium or low diet (20%, 17% or 12.2% crude protein and 67.9%, 66% or 62.9% total digestible nutrients, respectively) from 2 to 32 weeks of age, followed by medium nutrition. Semen collected from bulls at two specific time points, i.e. 55-59 and 69-71 weeks, was diluted, cryopreserved and used for reduced representation bisulfite sequencing. Differential methylation was detected for dietary treatment, but minimal differences were detected with age. The gene ontology term, "regulation of Rho protein signal transduction", implicated in sperm motility and acrosome reaction, was enriched in both low-vs-high and low-vs-medium datasets. Furthermore, several genes implicated in early embryo and foetal development showed differential methylation for diet. Our results therefore suggest that sperm epigenome keeps the memory of diet during pre-pubertal period in genes important for spermatogenesis, sperm function and early embryo development.

Sperm DNA methylation patterns at discrete CpGs and genes involved in embryonic development are related to bull fertility.

BACKGROUND: Despite a multifactorial approach being taken for the evaluation of bull semen quality in many animal breeding centres worldwide, reliable prediction of bull fertility is still a challenge. Recently, attention has turned to molecular mechanisms, which could uncover potential biomarkers of fertility. One of these mechanisms is DNA methylation, which together with other epigenetic mechanisms is essential for the fertilising sperm to drive normal embryo development and establish a viable pregnancy. In this study, we hypothesised that bull sperm DNA methylation patterns are related to bull fertility. We therefore investigated DNA methylation patterns from bulls used in artificial insemination with contrasting fertility scores. RESULTS: The DNA methylation patterns were obtained by reduced representative bisulphite sequencing from 10 high-fertility bulls and 10 low-fertility bulls, having average fertility scores of - 6.6 and + 6.5%, respectively (mean of the population was zero). Hierarchical clustering analysis did not distinguish bulls based on fertility but did highlight individual differences. Despite this, using stringent criteria (DNA methylation difference ≥ 35% and a q-value < 0.001), we identified 661 differently methylated cytosines (DMCs). DMCs were preferentially located in intergenic regions, introns, gene downstream regions, repetitive elements, open sea, shores and shelves of CpG islands. We also identified 10 differently methylated regions, covered by 7 unique genes (SFRP1, STXBP4, BCR, PSMG4, ARSG, ATP11A, RXRA), which are involved in spermatogenesis and early embryonic development. CONCLUSION: This study demonstrated that at specific CpG sites, sperm DNA methylation status is related to bull fertility, and identified seven differently methylated genes in sperm of subfertile bulls that may lead to altered gene expression and potentially influence embryo development.

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.

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.

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.

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.

Altered DNA methylation associated with an abnormal liver phenotype in a cattle model with a high incidence of perinatal pathologies.

Cloning enables the generation of both clinically normal and pathological individuals from the same donor cells, and may therefore be a DNA sequence-independent driver of phenotypic variability. We took advantage of cattle clones with identical genotypes but different developmental abilities to investigate the role of epigenetic factors in perinatal mortality, a complex trait with increasing prevalence in dairy cattle. We studied livers from pathological clones dying during the perinatal period, clinically normal adult clones with the same genotypes as perinatal clones and conventional age-matched controls. The livers from deceased perinatal clones displayed histological lesions, modifications to quantitative histomorphometric and metabolic parameters such as glycogen storage and fatty acid composition, and an absence of birth-induced maturation. In a genome-wide epigenetic analysis, we identified DNA methylation patterns underlying these phenotypic alterations and targeting genes relevant to liver metabolism, including the type 2 diabetes gene TCF7L2. The adult clones were devoid of major phenotypic and epigenetic abnormalities in the liver, ruling out the effects of genotype on the phenotype observed. These results thus provide the first demonstration of a genome-wide association between DNA methylation and perinatal mortality in cattle, and highlight epigenetics as a driving force for phenotypic variability in farmed animals.

Breeding animals for quality products: not only genetics.

The effect of the Developmental Origins of Health and Disease on the spread of non-communicable diseases is recognised by world agencies such as the United Nations and the World Health Organization. Early environmental effects on offspring phenotype also apply to domestic animals and their production traits. Herein, we show that maternal nutrition not only throughout pregnancy, but also in the periconception period can affect offspring phenotype through modifications of gametes, embryos and placental function. Because epigenetic mechanisms are key processes in mediating these effects, we propose that the study of epigenetic marks in gametes may provide additional information for domestic animal selection.

Genome-wide analysis of methylation in bovine clones by methylated DNA immunoprecipitation (MeDIP).

Methylated DNA immunoprecipitation (MeDIP), when coupled to high-throughput sequencing or microarray hybridization, allows for the identification of methylated loci at a genome-wide scale. Genomic regions affected by incomplete reprogramming after nuclear transfer can potentially be delineated by comparing the MeDIP profiles of bovine clones and non-clones. This chapter presents a MeDIP protocol largely inspired from Mohn and colleagues (Mohn et al., Methods Mol Biol 507:55-64, 2009), with PCR primers specific for cattle, and when possible, overviews of experimental designs adapted to the comparison between clones and non-clones.

Inositol 1,4,5-triphosphate receptor-binding protein released with inositol 1,4,5-triphosphate (IRBIT) associates with components of the mRNA 3' processing machinery in a phosphorylation-dependent manner and inhibits polyadenylation.

IRBIT is a recently identified protein that modulates the activities of both inositol 1,4,5-triphosphate receptor and pancreas-type Na(+)/HCO(3)(-) cotransporter 1, and the multisite phosphorylation of IRBIT is required for achieving this modulatory action. Here, we report the identification of the cleavage and polyadenylation specificity factor (CPSF), which is a multi-protein complex involved in 3' processing of mRNA precursors, as an additional binding partner for IRBIT. We found that IRBIT interacted with CPSF and was recruited to an exogenous polyadenylation signal-containing RNA. The main target for IRBIT in CPSF was Fip1 subunit, and the phosphorylation of the serine-rich region of IRBIT was required both for direct association with Fip1 in vitro and for redistribution of Fip1 into the cytoplasm of intact cells. Furthermore, tert-butylhydroquinone (tBHQ), an agent that induces oxidative stress, increased the phosphorylation level of IRBIT in vivo and in parallel enhanced the interaction between IRBIT and CPSF and promoted the cytoplasmic distribution of endogenous Fip1. In addition to CPSF, IRBIT interacted in vitro with poly(A) polymerase (PAP), which is the enzyme recruited by CPSF to elongate the poly(A) tail, and inhibited PAP activity in a phosphorylation-dependent manner. These findings raise the possibility that IRBIT modulates the polyadenylation state of specific mRNAs, both by controlling the cytoplasmic/nuclear partitioning of Fip1 and by inhibiting PAP activity, in response to a stimulus that alters its phosphorylation state.

IRBIT suppresses IP3 receptor activity by competing with IP3 for the common binding site on the IP3 receptor.

The inositol 1,4,5-trisphosphate (IP3) receptors (IP3Rs) are IP3-gated intracellular Ca2+ channels. We previously identified an IP3R binding protein, IRBIT, which binds to the IP3 binding domain of IP3R and is dissociated from IP3R in the presence of IP3. In the present study, we showed that IRBIT suppresses the activation of IP3R by competing with IP3 by [3H]IP3 binding assays, in vitro Ca2+ release assays, and Ca2+ imaging of intact cells. Multiserine phosphorylation of IRBIT was essential for the binding, and 10 of the 12 key amino acids in IP3R for IP3 recognition participated in binding to IRBIT. We propose a unique mode of IP3R regulation in which IP3 sensitivity is regulated by IRBIT acting as an endogenous "pseudoligand" whose inhibitory activity can be modulated by its phosphorylation status.

ZENON, a novel POZ Kruppel-like DNA binding protein associated with differentiation and/or survival of late postmitotic neurons.

The rat tyrosine hydroxylase gene promoter contains an E-box/dyad motif and an octameric and heptameric element that may be recognized by classes of transcription factors highly expressed during nervous system development. In a one-hybrid genetic screen, we used these sites as targets to isolate cDNAs encoding new transcription factors present in the brain. We identified ZENON, a novel rat POZ protein that contains two clusters of Kruppel-like zinc fingers and that presents several features of a transcription factor. ZENON is found in nuclei following transient transfection with the cDNA. The N-terminal zinc finger cluster contains a DNA binding domain that interacts with the E box. Cotranfection experiments revealed that ZENON induces tyrosine hydroxylase promoter activity. Unlike other POZ proteins, the ZENON POZ domain is not required for either activation of transcription or self-association. In the embryonic neural tube, ZENON expression is restricted to neurons that have already achieved mitosis and are engaged in late stages of neuronal differentiation (late postmitotic neurons). ZENON neuronal expression persists in the adult brain; therefore, ZENON can be considered a marker of mature neurons. We propose that ZENON is involved in the maintenance of panneuronal features and/or in the survival of mature neurons.