Front. Cell Dev. Biol. Frontiers in Cell and Developmental Biology Front. Cell Dev. Biol. 2296-634X Frontiers Media S.A. 10.3389/fcell.2020.00547 Cell and Developmental Biology Original Research Mating to Intact, but Not Vasectomized, Males Elicits Changes in the Endometrial Transcriptome: Insights From the Bovine Model Recuero Sandra 1 Sánchez José María 2 Mateo-Otero Yentel 1 Bagés-Arnal Sandra 2 McDonald Michael 2 Behura Susanta K. 3 Spencer Thomas E. 3 Kenny David A. 4 Yeste Marc 1 Lonergan Pat 2 Fernandez-Fuertes Beatriz 1 * 1Biotechnology of Animal and Human Reproduction (TechnoSperm), Department of Biology, Institute of Food and Agricultural Technology, University of Girona, Girona, Spain 2School of Agriculture and Food Science, University College Dublin, Dublin, Ireland 3Division of Animal Sciences, University of Missouri, Columbia, MO, United States 4Animal and Bioscience Research Centre, Teagasc Grange, Meath, Ireland

Edited by: Tomer Avidor-Reiss, The University of Toledo, United States

Reviewed by: Susan S. Suarez, Cornell University, United States; John J. Bromfield, University of Florida, United States

*Correspondence: Beatriz Fernandez-Fuertes, beatriz.fernandez@udg.edu

This article was submitted to Cell Growth and Division, a section of the journal Frontiers in Cell and Developmental Biology

10 07 2020 2020 8 547 05 05 2020 10 06 2020 Copyright © 2020 Recuero, Sánchez, Mateo-Otero, Bagés-Arnal, McDonald, Behura, Spencer, Kenny, Yeste, Lonergan and Fernandez-Fuertes. 2020 Recuero, Sánchez, Mateo-Otero, Bagés-Arnal, McDonald, Behura, Spencer, Kenny, Yeste, Lonergan and Fernandez-Fuertes

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

An appropriate female reproductive environment is essential for pregnancy success. In several species, including mice, pigs and horses, seminal plasma (SP) components have been shown to modulate this environment, leading to increased embryo viability and implantation. Due to the characteristics of mating in the aforementioned species, SP comes into direct contact with the uterus. However, it is questionable whether any SP reaches the uterus in species that ejaculate inside the vagina, such as humans and cattle. Hence, we hypothesized that sperm, perhaps acting as a vehicle for SP factors, play a more important role in the modulation of the maternal uterine environment in these species. In addition, changes elicited by SP and/or sperm may originate in the vagina and propagate to more distal regions of the female reproductive tract. To test these hypotheses, a bovine model in which heifers were mated to intact or vasectomized bulls or were left unmated was used. RNA-sequencing of endometrial samples collected 24 h after mating with a vasectomized bull did not reveal any differentially expressed genes (DEGs) in comparison with control samples. However, the endometrium of heifers mated with intact bulls exhibited 24 DEGs when compared to heifers mated with vasectomized bulls, and 22 DEGs when compared to unmated control heifers. The expression of a set of cytokines (IL6, IL1A, IL8, and TNFA) and candidate genes identified in the endometrial RNA-sequencing (PLA2G10, CX3CL1, C4BPA, PRSS2, BLA-DQB, and CEBPD) were assessed by RT-qPCR in the vagina and oviductal ampulla. No differences in expression of these genes were observed between treatments in any region. However, mating to both intact and vasectomized bulls induced an increase in IL1A and TNFA expression in the vagina compared to the oviduct. These data indicate that sperm, but not secretions from the accessory glands alone, induce modest changes in endometrial gene expression after natural mating in cattle. However, it is not clear whether this effect is triggered by inherent sperm proteins or SP proteins bound to sperm surface at the time of ejaculation.

seminal plasma cattle endometrium transcriptome maternal environment Horizon 202010.13039/501100007601Science Foundation Ireland10.13039/501100001602

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      Introduction

      Embryonic loss is a major contributor to pregnancy failure in livestock species and humans, ranging from 20 to 40% (Macklon et al., 2002; Wiltbank et al., 2016). Most of these losses occur before implantation, highlighting the importance of this period that encompasses such critical events as the first embryonic cleavage divisions; embryonic genome activation; blastocyst formation and hatching; conceptus development; and the preparation of the endometrium to interact with the embryonic trophectoderm (Diskin and Morris, 2008; Niakan et al., 2012; Sandra et al., 2017). Many factors are involved in implantation failure, but in mice and pigs there is growing evidence of a role for the maternal immune system and its regulation by seminal plasma (SP) (Gangnuss et al., 2004; O’Leary et al., 2004, 2006; Song et al., 2016; Glynn et al., 2017).

      Seminal plasma is a complex fluid resulting from the secretions of the testes, epididymides and accessory sex glands (in the bull: ampullae, seminal vesicles, prostate and bulbourethral glands). Although it is difficult to accurately calculate the precise contribution of each organ and gland to the final composition of this fluid, vasectomy in the bull by removal of a portion of each vas deferens and therefore removing the contribution of the epididymides, does not appear to significantly affect SP volume (Alexander et al., 1971), indicating a more prominent role of the accessory sex glands. However, vasectomy does lead to a reduction in amino acids in the bull ejaculate (Alexander et al., 1971), and to slight modifications in the proteome of human SP (Batruch et al., 2011). Traditionally, SP has been viewed as a mere vehicle for sperm that nourishes and supports these cells in the female reproductive tract. However, mounting evidence demonstrates an emerging role for SP components in the modulation of the endometrial and oviductal environment, which results in improved fertility and embryo survival and development (reviewed in Bromfeld, 2016; Morgan and Watkins, 2020). Exposure to SP in mice (Schjenken et al., 2015; Song et al., 2016; Glynn et al., 2017), as well as in pigs (O’Leary et al., 2004, 2006) and mares (Tunon et al., 2000; Palm et al., 2008; Fedorka et al., 2017), induces the expression of several endometrial cytokines, leading to leukocyte recruitment to the uterus. This migration of immune cells was thought to solely serve the purpose of clearing microorganisms and excess sperm (Pandya and Cohen, 1985; Thompson et al., 1992). However, it is now thought that the endometrial cytokine and chemokine cascade induced by SP is important to facilitate maternal tolerance toward paternal antigens (reviewed in Robertson, 2007). Indeed, in mice, mating drives the expansion of CD4+CD25+ T regulatory cells (Robertson et al., 2009; Shima et al., 2015), which can suppress or modulate the immune response of other cells (Sakaguchi et al., 2001). The increase in CD4+CD25+ is not observed when females are mated to vasectomized or seminal-vesicle-excised males, suggesting that this expansion is driven by secretions from the male accessory glands (Robertson et al., 2009). This effect likely explains why mating increases maternal tolerance toward paternal major histocompatibility complex (MHC) antigens (Robertson et al., 2009), which improves the ability of the semi-allogenic embryo to implant and develop normally in this species (Bromfield et al., 2014; Watkins et al., 2018). In addition to modifying the uterine environment, transcervical infusion of SP in pigs has been shown to modulate ovarian function by increasing corpora lutea (CL) weight and progesterone synthesis (O’Leary et al., 2006), which is essential for creating an appropriate uterine environment for the developing embryo. Moreover, in horses, a pivotal role of SP in protecting spermatozoa from neutrophil phagocytosis in the uterus has been suggested, improving fertility in this environment (Troedsson et al., 2002; Alghamdi et al., 2004).

      Due to characteristics of mating in rodents, pigs and horses, SP reaches the uterus and can therefore interact directly with the endometrium (Hunter, 1981; Dean et al., 2011). It is not clear, however, whether any SP reaches the uterus in species that ejaculate intravaginally and in which the volume of the ejaculate is relatively low, such as cattle or humans. It is possible that in those species, SP has an indirect effect on the endometrial environment and/or that sperm act as vehicles for the transport of SP components to more distal regions of the reproductive tract. In this sense, the bovine model could be more appropriate than rodents or pigs in understanding the regulatory properties of SP in the maternal environment of women.

      In vitro studies in humans have demonstrated the potential of SP to induce expression of cytokines and chemokines in vaginal, cervical and endometrial epithelial cell cultures (Gutsche et al., 2003; Sharkey et al., 2007, 2012a; Remes Lenicov et al., 2012). While gene expression changes and leukocyte recruitment have been described in the human cervix after unprotected, but not condom-protected, coitus (Sharkey et al., 2012b), there is currently no evidence of SP-induced changes in the endometrium in vivo.

      In cattle, the expression of several inflammatory mediators (such as colony-stimulating factor 2 – CSF2, interleukins 1B, 6, 17A and 8 - IL1B, IL6, IL17A, IL8; Prostaglandin-endoperoxide synthase 2 – PTGS2, and transforming growth factor beta 1 – TGF-B1) in uterine horns ipsi- and contralateral to the CL was modified after SP infusion into the uterus, in the absence or presence of sperm (Ibrahim et al., 2019). Despite this, uterine infusion of SP at the time of artificial insemination (AI) does not increase pregnancy rate in heifers or cows (Odhiambo et al., 2009; Ortiz et al., 2019). As mentioned above, it is questionable whether SP reaches the uterus during mating in cattle, so results obtained from the infusion of SP into the uterus may not be representative of physiological conditions. Indeed, recently it has been shown that infusion of SP into the vagina, but not into the uterus, modifies endometrial levels of epidermal growth factor (Badrakh et al., 2020), which highlights the importance of considering the ejaculate deposition site in natural conception in these studies. We recently reported a modest increase in conceptus length in embryos that developed from Day 7 to Day 14 in the uterus of heifers mated to a vasectomized bull in comparison to unmated heifers (Mateo-Otero et al., 2020). However, exposure of heifers to vasectomized bulls prior to AI failed to increase pregnancy rates (Pfeiffer et al., 2012). In addition, although bulls that had their seminal vesicles resected exhibited reduced semen volume, there was no apparent effect on their subsequent fertility (Shah et al., 1968). Together with recent work from our group, demonstrating a deleterious effect of bull SP on endometrial RNA integrity due to the presence of a seminal RNase (Fernandez-fuertes et al., 2019), the literature seems to suggest that SP does not play a significant role in pregnancy establishment in cattle.

      Based on these data, we hypothesized that in species that ejaculate inside the vagina, changes in the female reproductive environment begin in this region and then propagate to more distal regions, such as the uterus and/or the oviduct. Also, because of the lack of direct contact with the seminal fluid, sperm probably play a more important role in the modulation of the uterine environment in these species. In order to test these hypotheses, RNA-sequencing analysis of endometrial samples was carried out following natural mating of heifers with vasectomized (whose ejaculate lack sperm and epididymal and testicular fluid) or intact (that ejaculate sperm and SP) bulls. In addition, the expression of a set of interesting candidate genes was assessed in the vagina and oviductal ampulla, with the aim of determining the effects of sperm and accessory gland secretions from the most proximal region of the female reproductive tract (vagina) to the distal region where gamete interaction takes place (oviductal ampulla).

      Materials and Methods

      Unless otherwise stated, all chemicals and reagents were sourced from Sigma-Aldrich (Arklow, Ireland).

      Animals

      All experimental procedures involving animals were approved by the Animal Research Ethics Committee of University College Dublin and licensed by the Health Products Regulatory Authority (HPRA), Ireland, in accordance with Statutory Instrument No. 543 of 2012 (under Directive 2010/63/EU on the Protection of Animals used for Scientific Purposes). For the duration of the study, all animals were housed in groups of 10–15, independent of treatment, and managed identically in terms of feeding and husbandry routines.

      Vasectomy was carried out by removing approximately 5 cm of both vasa deferentia. This procedure took place 5–6 months prior to the trial. Vasectomized bulls underwent semen evaluation to confirm the lack of sperm and all ran as teasers for oestrus detection with 25 cows/heifers each during the breeding season preceding the study. Intact bulls underwent a breeding soundness evaluation prior to the study. Neither intact and vasectomized bulls had access to females for at least 5 months before the trial, nor during it (apart from the controlled mating to the experimental heifers).

      Experimental Design

      Estrous cycles of crossbreed beef heifers (Angus and Holstein-Friesian cross; n = 28) were synchronized using an 8-day intravaginal device (PRID® Delta, 1.55 g progesterone, Ceva Santé Animale, Libourne, France), together with a 2 mL intramuscular injection of a synthetic gonadotrophin releasing hormone (Ovarelin®, equivalent to 100 μg Gonadorelin, Ceva Santé Animale) administered on the day of PRID insertion. One day prior to PRID removal, all heifers received a 5 mL intramuscular injection of prostaglandin F2 alpha (Enzaprost®, equivalent to 25 mg of Dinoprost, Ceva Santé Animale) to induce luteolysis. Only heifers observed in standing estrus were used (n = 22). Heifers were blocked by weight and randomly allocated to one of three treatments (Figure 1): (1) mated to an intact bull (n = 7), (2) mated to a vasectomized bull (n = 8), or (3) left unmated (control; n = 7). Between 0 to 6 h after estrus detection, heifers were separated from the group and placed in a pen (one at a time) with one of three vasectomized Holstein Friesian bulls, or one of two intact Holstein Friesian bulls (Supplementary Table 1). Once the bull mounted and intromission was confirmed, the heifers were returned to the group. Bulls were allowed to mate no more than twice per day and the experiment was carried out over three consecutive days.

      Summary of the experimental design. Only heifers observed in standing estrus were mated with (1) an intact bull (blue; n = 7), (2) a vasectomized bull (green; n = 8) or (3) left unmated (gray; n = 7). A schematic representation of the male reproductive tract is shown, and the structures that contribute to the ejaculate in each case are highlighted (T, testis; E, epididymis; GP, glans penis; VD, vas deferens; B, bladder; A, ampulla; SV, seminal vesicle; P, prostate; BG, bulbourethral gland). Heifers were slaughtered 24 h after mating and their reproductive tracts were recovered. Endometrial samples were obtained from the base of the ipsilateral uterine horn to perform RNA sequencing (Experiment 1). Samples from ampulla and vagina were obtained to assess the expression of a set of genes (PLA2G10, CX3CL1, C4BPA, PRSS2, BLA-DQB, CEBPD, IL1A, IL6, TNFA, and IL8) by RT-qPCR (Experiment 2).

      Tissue Collection

      All heifers were slaughtered in a commercial abattoir 24 h (±6 h) after mating and their reproductive tracts were recovered. The ovaries were examined to determine the site of the preovulatory or freshly ovulated follicle (Supplementary Figure 1). Endometrial tissue samples were obtained from intercaruncular areas of the base of uterine horn ipsilateral to the preovulatory or freshly ovulated follicle. In cattle, the uterine glands are located in the intercaruncular areas of the endometrium, while the caruncular areas are aglandular. These glands are responsible for the secretion of the histotroph, which will nurture the developing embryo, and drive the maternal signals of implantation (Kelleher et al., 2019; Spencer et al., 2019). In addition, epithelial sections of anterior vagina and whole sections of the ampulla of the ipsilateral oviduct were obtained. Immediately after tissue collection, each sample was snap frozen in liquid nitrogen and stored at −80°C.

      Experiment 1: Seminal Plasma Effects on the Endometrial Transcriptome RNA Extraction

      For total mRNA extraction, samples were first homogenized in Trizol reagent (Invitrogen, Carlsbad, CA, United States) using a steel bead and the Qiagen tissue lyzer (2 × 120 s at maximum speed). On-column RNA purification was performed using the Qiagen RNeasy kit (Qiagen, Crawley, Sussex, United Kingdom) per the manufacturer’s instructions. The quantity of RNA was determined using the Nano Drop 1000 spectrophotometer (Thermo Fisher Scientific, Dublin, Ireland). Prior to endometrial RNA sequencing analysis, the RNA quality was assessed by the Agilent Bioanalyzer (Agilent Technologies, Cork, Ireland). Only samples that exhibited a minimum RNA integrity number (RIN) of 8 were used in this experiment (n = 6 heifers in each experimental group).

      RNA Sequencing Analysis

      RNA library preparation and sequencing were performed by the University of Missouri DNA Core Facility as described previously by Moraes et al. (2018). The raw sequences (fastq) were subjected to quality trimming control using fqtrim1. Then, the quality reads were mapped to the bovine reference genome UMD3.1 using Hisat2 mapper (Kim et al., 2015). Read counts mapping to each gene were determined from the binary alignment map files of the samples using FeatureCounts (Liao et al., 2014). Differential expression analysis between different sample groups was performed using edgeR robust (Zhou et al., 2014).

      For the annotated DEGs, the gene ontology analysis was performed using PANTHER2.

      Experiment 2: Seminal Plasma Effects on Gene Expression in the Vagina and Oviduct RNA Extraction and cDNA Synthesis

      Total mRNA extraction was carried out as described above. For each sample, cDNA was prepared from approximately 100 ng of total mRNA using the High Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific) according to the manufacturer’s instructions. The cDNA obtained was diluted using RNase- and DNase-free water in a final concentration of 5 ng/mL and in a total volume of 100 μL. The diluted cDNA samples were stored at −20°C for subsequent analysis. As some samples were lost during the management of the experiment, the number of samples per group varies for the vagina region: n = 6 heifers in the control group, n = 5 in the intact group and n = 7 in the vasectomized group. For the ampulla region, there were n = 7 individuals in each experimental group (control, intact, and vasectomized). In addition, the RNA quality was assessed by the Agilent Bioanalyzer (Agilent Technologies, Cork, Ireland) (RIN = 7.81 ± 0.29; mean ± standard error of the mean, SEM).

      Genes of Interest Selection for Expression Analysis

      In order to determine SP effects on gene expression in vagina and oviduct regions, a set of genes (IL1A, IL6, TNFA, and IL8) were selected based on the literature. These genes were inflammatory mediators expression of which was reported to be modified by SP exposure in cattle and other species (O’Leary et al., 2004; Sharkey et al., 2012b; Schjenken et al., 2015; Introini et al., 2017; Ibrahim et al., 2019). Moreover, expression of some genes (PLA2G10, CX3CL1, C4BPA, PRSS2, BLA-DQB, and CEBPD) found differentially expressed in RNA-sequencing analysis of endometrium samples was also interrogated. The selection of these genes as interesting targets was based on literature searching (Wedel and Lömsziegler-Heitbrock, 1995; Blom et al., 2004; Hannan et al., 2004; Tapia et al., 2008; Hayashi et al., 2017; Neupane et al., 2017; Pinto De Melo et al., 2017; Tríbulo et al., 2018) and GO term analysis.

      Quantitative Real-Time PCR Analysis

      All primers were designed using Primer Blast software3 (Supplementary Table 2). Briefly, RT-qPCR assays were performed per duplicate in a total volume of 20 μL, containing 10 μL FastStart Universal SYBR Green Master (Roche Diagnostics Ltd., West Sussex, United Kingdom), 1.2 μL forward and reverse primer mix (300 nM final concentration), 5.6 μL nuclease-free water and 2 μL cDNA template on the ABI Prism 7500 Real-Time PCR System (Life Technologies). A total of 40 cycles were performed with the following thermo-cycling conditions for each cycle: 50°C for 2 min, 95°C for 10 min followed by 95°C for 15 s, 60°C for 1 min, 95°C for 30 s and 60°C for 15 s. The melt curve was also included to ensure specificity of amplification. The specificity of all targets was confirmed by the presence of a single sharp peak in the melt curve. A total of eight potential reference genes [Glyceraldehyde 3-Phosphate Dehydrogenase (GAPDH), Actin Cytoplasmic 1 (ACTB), 60S Ribosomal Protein L18 (RPL18), Peptidyl-Prolyl Cis-Trans Isomerase A (PPIA), 14-3-3 Protein Zeta/Delta (YWHAZ), RING Finger Protein 11 (RNF11), Histone H3.3 (H3F3A), Succinate Dehydrogenase Complex Subunit A Flavoprotein Variant (SDHA)] were analyzed using the geNorm function with the qbase + package (Biogazelle, Zwijnaarde, Belgium) to identify the best reference genes. Due to the high variability between samples, a total of four reference genes were selected: RNF11, H3F3A, YWHAZ, and GADPH, which were the most stably expressed (average geNorm M ≤ 0.5).

      Primer efficiency was carried out for the genes of interest, and RT-qPCR of 1:4 dilutions of a cDNA mix from a representative pool of samples were analyzed. The presence of a single sharp peak in the melt curve as well as the standard curve was used to confirm primer specificity. The threshold cycle (Ct) for each sample was automatically calculated using the default settings within the SDS software (SDS 1.4, ABI). In order to obtain the relative expression values of the genes of interest, 2–ΔΔCT method was used (Livak and Schmittgen, 2001). For each individual, the expression of the genes of interest was firstly normalized to the average of housekeeping genes previously selected (RNF11, H3F3A, YWHAZ, and GADPH) with the following formula: ΔCt = Ctgene of interest – Ct(RNF11+H3F3A+YWHAZ+GADPH)/4. The values of ΔΔCt were calculated normalizing the results to the mean across all individuals, including both tissue regions (vagina and ampulla), per each gene of interest. The subsequent statistical analysis was performed using ΔCt values whereas the results are represented as 2–ΔΔCT.

      Results expressed as ΔCT were analyzed with IBM SPSS 25.0 for Windows (Armonk; New York, NY, United States). Data were checked for normal distribution (Shapiro–Wilk test) and homoscedasticity (Levene test) to confirm that parametric assumptions were fulfilled. When these premises were not, data (x) were linearly transformed using the square root (√x) and arcsine of the square root (arcsin √x). Thereafter, data (transformed or not depending on the case) were analyzed by a two-way ANOVA followed by a Sidak post hoc test for pair-wise comparisons. The expression of five genes (CX3CL1, PLA2G10, TNFA, IL6, and CXCL8), even after linear transformation, did not match parametric assumptions. For this reason, Scheirer–Ray–Hare and Mann–Whitney tests were used as non-parametric alternatives. In all cases, the significance level was established at P ≤ 0.05.

      Results Ovary Status

      At the time of sample collection (24 ± 6 h after mating), a total of 13 heifers had a freshly ovulated follicle on their ovary, while the remaining animals exhibited a pre-ovulatory follicle (see Supplementary Table 1 and Supplementary Figure 1). In the control group, three animals had ovulated and four exhibited a pre-ovulatory follicle. Regarding the heifers mated with vasectomized bulls, fresh ovulation was found in six animals and only two presented a pre-ovulatory follicle. While in the intact group, four ovulations and three preovulatory follicles were observed. The proportion of animals that had a fresh ovulation was balanced across treatments for subsequent gene expression analysis.

      Effects of Seminal Plasma on the Endometrial Transcriptome

      Sequencing of endometrial samples of heifers recovered 24 h after mating to intact bulls revealed a total of 22 differentially expressed genes (DEGs) compared with contemporary unmated animals (Table 1 and Supplementary File 1). Of those DEGs, 12 were up-regulated and 10 down-regulated [False discovery rate (FDR) < 0.05]. Some of the genes that exhibited the lowest expression (logFC < -2) were serine protease 2 (PRSS2), complement C9 (C9), oxytocin/neurophysin I prepropeptide (OXT), a novel gene encoding for carbonic anhydrase 1 (ENSBTAG00000036116) and an uncharacterized novel gene (ENSBTAG00000050072). On the other hand, the genes with greater transcript abundance (logFC > 2) levels were coiled-coil domain containing 196 (CCDC196), solute carrier family 24 member 2 (SLC24A2), UDP glucuronosyltransferase family 2-member A1 complex locus (UGT2A1) and interferon gamma inducible protein 47 (IFI47). In contrast, the endometrium of heifers exposed only to SP (by mating with a vasectomized bull) did not exhibit DEGs compared with the control group. Comparison of endometrial transcriptomes of intact and vasectomized groups revealed a total of 24 DEGs, 18 up-regulated and 6 down-regulated (Table 2 and Supplementary File 1) (FDR < 0.05). Amongst these, MHC, class II, DQ beta (BOLA-DQB), GSG1 like (GSG1L), potassium voltage-gated channel subfamily E regulatory subunit 1 (KCNE1) and the novel gene previously mentioned (ENSBTAG00000050072) were those that displayed lower logFC values (logFC < −2). In contrast, higher levels of expression (logFC > 2) were exhibited by interleukin 17F (IL17F), complement component 4 binding protein alpha (C4BPA), the aforementioned IFI47 and UGT2A1, and a novel gene (ENSBTAG00000052851) which has been predicted to code for a protein containing an Ig-like domain.

      List of differentially expressed genes (FDR < 0.05) in endometrial samples of heifers mated with intact bulls compared with unmated heifers.

      Ensembl acc. number Gene name Gene description logFC
      ENSBTAG00000014234 CCDC196 Coiled-coil domain containing 196 8.63
      ENSBTAG00000043972 SLC24A2 Solute carrier family 24 member 2 4.15
      ENSBTAG00000004040 UGT2A1 UDP glucuronosyltransferase family 2 member A1 complex locus 3.11
      ENSBTAG00000003529 ASAH2 N-acylsphingosine amidohydrolase 2 2.82
      ENSBTAG00000015727 IFI47 Interferon gamma inducible protein 47 2.62
      ENSBTAG00000049426 STARD2 Phosphatidylcholine transfer protein 1.72
      ENSBTAG00000037929 ADAM28 ADAM-like, decysin 1 1.42
      ENSBTAG00000019636 SCARA5 Scavenger receptor class A member 5 1.06
      ENSBTAG00000017722 F5 Coagulation factor V 0.99
      ENSBTAG00000019625 EHHADH Enoyl-CoA hydratase and 3-hydroxyacyl CoA dehydrogenase 0.83
      ENSBTAG00000008735 VASH1 Vasoinhibin 1 0.77
      ENSBTAG00000001728 IGSF10 Immunoglobulin superfamily member 10 0.72
      ENSBTAG00000046307 CEBPD CCAAT enhancer binding protein delta −0.86
      ENSBTAG00000011079 C18H19orf48 Chromosome 18 C19orf48 homolog −0.92
      ENSBTAG00000007101 F3 Coagulation factor III, tissue factor −1.04
      ENSBTAG00000051812 CA1L Carbonic anhydrase 1-like −2.45
      ENSBTAG00000008026 OXT Oxytocin/neurophysin I prepropeptide −2.76
      ENSBTAG00000039446 PI3L Elafin-like −2.78
      ENSBTAG00000036116 CA1 Carbonic anhydrase 1 −3.39
      ENSBTAG00000016149 C9 Complement C9 −4.08
      ENSBTAG00000050072 Novel gene −4.99
      ENSBTAG00000021565 PRSS2 Serine protease 2 −6.05
      FDR, false discovery rate; logFC, logarithm of fold change.

      List of differentially expressed genes (FDR < 0.05) in endometrial samples of heifers mated with intact bulls compared with heifers mated with vasectomized bulls.

      Ensembl acc. number Gene name Gene description logFC
      ENSBTAG00000052851 Novel gene 4.57
      ENSBTAG00000016835 IL17F Interleukin 17F 4.50
      ENSBTAG00000032884 TNP2 Transition protein 2 3.78
      ENSBTAG00000019132 DMP1 Dentin matrix acidic phosphoprotein 1 3.34
      ENSBTAG00000015727 IFI47 Interferon gamma inducible protein 47 3.18
      ENSBTAG00000009876 C4BPA Complement component 4 binding protein alpha 2.78
      ENSBTAG00000004040 UGT2A1 UDP glucuronosyltransferase family 2 member A1 complex locus 2.47
      ENSBTAG00000037539 Vascular cell adhesion molecule 1-like 1.81
      ENSBTAG00000021764 GLRB Glycine receptor beta 1.62
      ENSBTAG00000002214 TAT Tyrosine aminotransferase 1.38
      ENSBTAG00000026779 LYZ Lysozyme 1.32
      ENSBTAG00000000601 COL11A2 Collagen type XI alpha 2 chain 1.14
      ENSBTAG00000019588 BLA-DQB MHC class II antigen 1.11
      ENSBTAG00000034338 C15H11orf88 Chromosome 15 C11orf88 homolog 1.10
      ENSBTAG00000021526 RPRM Reprimo, TP53 dependent G2 arrest mediator homolog 1.05
      ENSBTAG00000033429 FAM229B Family with sequence similarity 229 member B 0.94
      ENSBTAG00000024869 CX3CL1 C-X3-C motif chemokine ligand 1 0.78
      ENSBTAG00000021522 PLA2G10 Group 10 secretory phospholipase A2 0.72
      ENSBTAG00000012703 GLO1 Glyoxalase I −0.49
      ENSBTAG00000008147 MICAL1 Microtubule associated monooxygenase, calponin and LIM domain containing 1 −0.50
      ENSBTAG00000001150 KCNE1 Potassium voltage-gated channel subfamily E regulatory subunit 1 −2.51
      ENSBTAG00000004607 GSG1L GSG1 like −4.61
      ENSBTAG00000050072 Novel gene −4.94
      ENSBTAG00000021077 BOLA-DQB Major histocompatibility complex, class II, DQ beta −5.05
      FDR, false discovery rate; logFC: logarithm of fold change.

      Three DEGs (UGT2A1, IFI47 and the novel gene ENSBTAG00000050072) were found to be common of DEGs detected between the intact group and the control and those detected when comparing the intact and vasectomized groups.

      Gene Ontology (GO) Terms of DEGs

      For the annotated genes in each comparison, the GO terms are shown in Figure 2. For the molecular function category, the main represented GO term was “catalytic activity” in the intact group compared with control or vasectomized samples (Figure 2A). “Cellular process” and “metabolic process” were the most represented terms for the biological process category in intact samples compared with the control (Figure 2B). Compared with the vasectomized group in the same category, in addition to “cellular process,” “response to stimulus” was the most represented term (Figure 2B). In regard to pathway category, comparing the intact group with the control, all the terms represented were related to vascular regulation (Figure 2D). In contrast, compared with the vasectomized bull treatment group, the represented terms referred to immunity modulation and amino acid biosynthesis (Figure 2D). “Hydrolase” and “receptor” were the most represented protein class terms among the differentially regulated genes in the endometrium exposed to sperm and SP compared with the control samples (Figure 2C). On the other hand, compared with vasectomized samples, the most represented protein class terms were “cytoskeletal protein” and “immunity protein” (Figure 2C).

      (A–D) Representation of GO terms for annotated DEGs in each comparison. The main categories: (A) molecular function, (B) biological process, (C) protein class, and (D) pathway were represented. IMBCCSP, inflammation mediated by chemokine and cytokine signaling pathway.

      Effects of Seminal Plasma on Gene Expression in the Vagina and Oviduct

      Based on the results obtained from the endometrial RNA-sequencing, we were interested in studying whether gene expression changes are more dramatic at the site of semen deposition (the vagina), and whether these changes can propagate to more distal regions (the oviductal ampulla). Thus, six DEGs from the sequencing analysis were selected (PLA2G10, CX3CL1, C4BPA, PRSS2, BLA-DQB, and CEBPD) based on evidence of their reproductive function found in the literature (Wedel and Lömsziegler-Heitbrock, 1995; Blom et al., 2004; Hannan et al., 2004; Tapia et al., 2008; Hayashi et al., 2017; Neupane et al., 2017; Pinto De Melo et al., 2017; Tríbulo et al., 2018), as well as their GO terms. In addition to these, IL6, IL1A, TNFA, and IL8 expression was also assessed, as these are genes that have been observed to be regulated by SP in several species (O’Leary et al., 2004; Sharkey et al., 2012b; Schjenken et al., 2015; Introini et al., 2017; Ibrahim et al., 2019).

      Differences in relative expression of CEBPD (P < 0.01; Figure 3F) and IL8 (P < 0.05 in the control group and P < 0.01 in the intact and vasectomized groups; Figure 3J) between regions were observed in all groups, being up-regulated in the vagina in comparison with the ampulla, whereas CX3CL1 was down-regulated in the vagina compared with the ampulla, only in the control group (P < 0.01; Figure 3B). In addition, TNFA (P < 0.05 in the intact group and P < 0.01 in the vasectomized group; Figure 3I) and IL1A (P < 0.05; Figure 3H) were up-regulated in the vagina compared with the ampulla in heifers that had been mated to an intact or a vasectomized bull, but not in unmated heifers. The remaining genes did not exhibit region-specific changes (Figures 3A,C–E,G). Conversely, when relative abundance of these genes was compared between treatment groups, no differences were detected (P > 0.05). It is also important to note that for many genes, especially those related to inflammation, there was considerable variability between animals.

      (A–J) Relative expression values of all the genes assessed: (A) PLA2G10, (B) CX3CL1, (C) C4BPA, (D) PRSS2, (E) BLA-DQB, (F) CEBPD, (G) IL6, (H) IL1A, (I) TNFA, (J) IL8 in the vagina (circles) and oviductal ampulla (triangle) in the different experimental groups (control, intact or vasectomized). Each circle/triangle refers to an individual heifer. Different color shades within groups correspond to the bull that was mated to that particular heifer, blue for intact group (n = 2 bulls) and green for vasectomized group (n = 3). Bars represent the mean of relative expression values and asterisks indicate significant (P < 0.05, ∗∗P < 0.01) differences.

      Discussion

      The main findings of this study are: (1) mating to an intact bull induces subtle changes in the endometrial transcriptome; however, (2) these transcriptomic changes are not observed in heifers mated to vasectomized bulls; (3) expression of PLA2G10, CX3CL1, C4BPA, PRSS2, BLA-DQB, CEBPD, IL1A, IL6, TNFA, and IL8 in the vagina and ampulla did not differ between treatments; and (4) TNFA and IL1A exhibited regional differences between vagina and ampulla of heifers mated to intact or vasectomized bulls.

      Seminal plasma is a complex fluid the composition of which is determined by the size, storage capacity, and secretory output of different organs of the male reproductive tract, which in the bull include: testes, epididymides, ampullae of the vasa deferentia, seminal vesicles, prostate and bulbourethral glands. After completion of spermatogenesis, sperm enter the epididymis bathed in fluid produced by the rete testis, which will be absorbed in its majority by the epididymal epithelium (Amann et al., 1974). However, secreted factors of epididymal and/or testicular origin are found in the ejaculate. This is evidenced by studies comparing ejaculates before and after vasectomy, which show lower concentration of amino acids in bulls (Alexander et al., 1971), and proteins in humans (Batruch et al., 2011), after the procedure. Despite this, vasectomy does not lead to a decrease in bull ejaculate volume (Alexander et al., 1971), indicating a more prominent role of the male accessory glands in the production of SP components. It is challenging to precisely calculate the contribution of each gland to the final fluid volume and composition (Seidel and Foote, 1970); however, vesiculectomy (excision of the seminal vesicles) in the bull leads to a more than 50% decrease in ejaculate volume, together with a reduction in total protein and ion concentrations (chloride, potassium, calcium, sodium), lower sperm motility, viability and morphology, and increased pH (Faulkner et al., 1968; Alexander et al., 1971).

      Growing evidence exists for a role of SP in the modulation of cellular and molecular events in the maternal tract of several species during early pregnancy (Gangnuss et al., 2004; O’Leary et al., 2004, 2006; Bromfeld, 2016; Song et al., 2016; Glynn et al., 2017; Morgan and Watkins, 2020). However, most of the in vivo evidence comes from pigs and mice, species in which SP reaches the uterus (Gangnuss et al., 2004; O’Leary et al., 2004, 2006; Bromfield et al., 2014; Song et al., 2016; Glynn et al., 2017). Because of the relatively low volume of ejaculate in men (an average of 3.7 mL, Cooper et al., 2009) and bulls (around 5 mL), and the fact that the ejaculate is deposited in the vagina, it is questionable whether any SP reaches the uterus. Therefore, it is not clear whether this fluid has a critical role in the modulation of the uterine environment in these species. At the time of ejaculation, however, sperm come into contact with SP, leading to proteins binding tightly to the sperm plasma membrane (Pini et al., 2016). For example, seminal vesicle-derived Binder of Sperm Proteins (BSP) −1, −3, and −5 (previously called PDC-109 or BSP-A1/A2, BSP-A3 and BSP-30 kDa respectively), which make up approximately 50% of total protein in SP (Nauc and Manjunath, 2000), bind to sperm and play important roles during capacitation (Manjunath and Thérien, 2002) and formation of the sperm oviductal reservoir (Gwathmey et al., 2003, 2006). Thus, it is possible that in intravaginal ejaculators, sperm can act as a vehicle of SP proteins that interact with the reproductive epithelium to induce changes in the maternal environment. Indeed, bull sperm has been shown to interact with endometrial cells and induce a proinflammatory response in vitro (Elweza et al., 2018; Ezz et al., 2019).

      To address this lack of basic knowledge, this study aimed to determine the effects of bovine SP and sperm exposure during natural mating on the endometrial transcriptome. Additionally, gene expression changes were assessed at the site of semen deposition (vagina) and the distal region where gamete interaction takes place (the oviductal ampulla) to determine whether SP-induced changes can propagate throughout the female reproductive tract.

      In the present study, heifers were mated between 0 and 6 h after seen in standing estrus. Average time from estrus onset to ovulation is 27 h (Walker et al., 1996; Valenza et al., 2012; Randi et al., 2018). Thus, some animals had ovulated by the time of sample collection (24 ± 6 h after mating) whereas the rest exhibited a pre-ovulatory follicle. Ovulation and estrous cycle are orchestrated by an accurate hormonal regulation, and under this regulation, the endometrium experiments functional and morphological changes (Arai et al., 2013). In order to avoid any possible confounding factors due to ovulation having occurred or not, the number of ovulated and non-ovulated heifers that were analyzed by RNA-sequencing and RT-qPCR was balanced between treatments.

      Strikingly, when heifers were exposed to SP in the absence of sperm and testicular and epididymal secretions (i.e., mated to a vasectomized bull) the endometrial RNA-sequencing analysis did not reveal any DEGs in comparison with samples from unmated animals. Conversely, the endometrial transcriptome of heifers mated to intact bulls differed from the control and vasectomized groups, exhibiting differential regulation of a small number of genes that may play a role in bovine fertility. Although, these results could be due to SP reaching the uterus in both treatments, but only testicular of epididymal factors inducing a response; it is more likely that, in cattle, SP does not reach the uterus in the 24 h following natural mating, at least in the absence of a vehicle, such as sperm.

      When compared with control endometrial samples, endometrium obtained from heifers mated to an intact bull exhibited up-regulation of 12 genes and down-regulation of 10 genes. Amongst these genes regulated by mating, some have been shown to participate in tissue remodeling, an important step preparing endometrium to embryo implantation. The gene coding for scavenger receptor class A, member 5 (SCARA5), which participates in innate immunity (Jiang et al., 2006) was up-regulated. This gene has also been observed to be up-regulated in the endometrium of cows at day 20 of pregnancy (Mansouri-Attia et al., 2009), and has been proposed to play a role in the regulation of histotroph secretion and tissue remodeling, two critical processes for embryo implantation (Vitorino Carvalho et al., 2019). Additionally, serine protease 2 (PRSS2), which also participates in tissue remodeling by type 1 collagen degradation, was down-regulated by exposure to SP and sperm at mating. Interestingly, this gene is up-regulated in the endometrium of repeat breeder cows, those that are cycling normally and without clinical abnormalities but that fail to conceive after at least two successive inseminations (Hayashi et al., 2017). Mating to an intact bull also affected genes involved in cell proliferation, such as CCAAT enhancer binding protein delta (CEBPD), which was found to be down-regulated in the endometrium of heifers mated with intact bulls. This gene belongs to the C/EBP leucine-zipper transcription factor family involved in fat and hematopoietic progenitor cells differentiation (Wedel and Lömsziegler-Heitbrock, 1995). Another member of this family, CEBPB, has been identified as a regulator of proliferative events during decidualization in mice (Mantena et al., 2006). In regard to modulation of innate immunity, the gene coding for component 9 of complement system (C9) was down-regulated in the endometrium of heifers mated with intact bulls. The C9 component participates in the final steps of the complement cascade, in the formation of membrane attack complex (MAC), which mediates the formation of channels in the target cell membrane, leading to cell lysis and death (Janeway et al., 2001). Despite its importance, the function of the complement system in the context of reproduction is not well known. For instance, complement regulatory proteins have been found in bull sperm surface, such as CD59 (Byrne et al., 2012), which prevents the formation of MAC (Janeway et al., 2001). These complement regulatory proteins were also identified in human and mouse sperm and they have been proposed to play a role protecting sperm in the female tract (Harris et al., 2006). Moreover, mating with an intact bull resulted in a down-regulation of endometrial oxytocin (OXT). In cattle, high levels of OXT have been reported to impair embryo survival by promoting uterine secretion of prostaglandin F (PGF), which induces luteolysis and consequently a drop in progesterone (Lemaster et al., 1999). Although endometrial OXT production has been observed also in mares (Bae and Watson, 2003), it is not known the locally effect of the OXT secreted by the endometrium during early pregnancy in cattle.

      When comparing endometrial samples from heifers mated to intact bulls to those from heifers mated to vasectomized bulls, 18 genes were found to be up-regulated and four down-regulated. Amongst the up-regulated genes, literature on CX3CL1, VCAM1-like, C4BPA, PLA2G10, IFI47, IL17F and BLA-DQB suggest different roles of these genes in early pregnancy in different species. The chemokine CX3CL1 (C-X3-C motif chemokine ligand 1) has been identified as a potential bovine embryokine (Tríbulo et al., 2018). In addition, CX3CL1 induces recruitment of leukocytes during early pregnancy (Hannan et al., 2004), and also promotes trophoblast migration in women (Hannan et al., 2006). With importance for implantation, VCAM1 has been shown to be involved in the adhesion of the bovine conceptus to the endometrium (Bai et al., 2014). Further, in women with unexplained infertility the endometrial expression of VCAM1 at the peri-implantation stage was significantly lower than control women (Konac et al., 2009). Additionally, the mRNA levels of complement component 4-binding protein alpha (C4BPA), a key inhibitor of the complement system (Blom et al., 2004), were increased during the implantation window in women (Tapia et al., 2011), but decreased in women which suffered repeated implantation failure and unexplained recurrent spontaneous abortion (Lee et al., 2007; Tapia et al., 2008). Although cattle and human implantation differs significantly, in both the time at which it takes place (around day 9 in humans and starting at day 21 in cattle) and the structure of the placenta (hemochorial in humans and epitheliochorial in bovine), it is likely that the up-regulation of C4BPA and VCAM1 induced by exposure to SP and sperm during mating regulates peri-implantation events in cattle. On the other hand, group 10 secretory phospholipase A2 (PLA2G10), which was up-regulated in intact group samples, belongs to phospholipase A2 enzyme family, which is known to participate in inflammatory processes and to catalyze the release of arachidonic acid from phospholipids, needed to prostaglandin production (reviewed in Capper and Marshall, 2001). In addition, this gene has been associated with fertility in beef cattle (Neupane et al., 2017), and it was found highly down-regulated in the uterus of cows with negative energy balance (Wathes et al., 2009) which typically have lower pregnancy rates.

      From an immunological point of view, the expression of interferon gamma inducible protein 47 (IFI47) was also found to be up-regulated in the endometrium of heifers mated with intact bulls in comparison with those mated to vasectomized bulls. Although its function in the uterus has not been defined, IFI47 mRNA was more abundant in the endometrium of high fertility heifers compared to heifers classified as infertile (Minten et al., 2013). In addition, IFI47 was up-regulated in the endometrium of heifers 13 days after embryo transfer (Spencer et al., 2013), probably under stimulation of interferon-tau secreted by the conceptus. Another immune-related gene up-regulated by mating with intact bulls is interleukin 17 F (IL17F), a cytokine produced by T helper 17 lymphocytes (reviewed in Iwakura et al., 2011). A similar response in expression of IL17A was observed by Ibrahim et al. (2019). Both IL17A and IL17F share biological functions; indeed, both are highly homologous, can bind to the same receptor and moreover, can be secreted as heterodimers or homodimers (reviewed in Iwakura et al., 2011). Ibrahim et al. (2019) observed increased levels of IL17A after exposing endometrial cell cultures to SP and sperm or sperm alone, but not SP in the absence of sperm. Also, up-regulation of IL17A took place in vivo, after uterine infusion of semen (a combination of SP and sperm) but not SP alone (Ibrahim et al., 2019). Together with the data obtained in the present study, this suggests that regulation of IL17 expression in the bovine endometrium is mediated by sperm action, not SP. This is interesting to note as, in mice, mating to a vasectomized male induces a similar increase in endometrial IL17A expression than mating with an intact male (Song et al., 2016). This gives weight to our hypothesis that in species that ejaculate intravaginally, sperm play a more important role than SP in regulation of the female reproductive environment. Moreover, the endometrial up-regulation of IL17A in bovine explants has been shown to be exclusively regulated by elongated Day 15 conceptuses, but not by interferon-tau, which is the main signal of pregnancy recognition in cattle, suggesting a role in the embryo – endometrium crosstalk during early pregnancy in cattle (Sánchez et al., 2019). During pregnancy, the maternal immune system must tolerate the presence of an embryo that expresses paternal antigens. The increase in endometrial IL17 expression observed after mating in the aforementioned studies is reflective of an increase in the population of T helper 17 lymphocytes (Song et al., 2016), which probably participate in the establishment of this maternal tolerance toward paternal antigens. Another mechanism of maternal tolerance driven by male factors seems to be the regulation of the MHC. In the present study, mating to an intact bull elicited down-regulation of Bovine Lymphocyte antigen (referred to as MHC in other species) class II, DQ beta (BOLA-DQB), consistent with its down-regulation in pregnant heifers after natural breeding (Dickinson et al., 2018). Conversely, another MHC class II member, BLA-DQB was up-regulated in endometrial samples of heifers mated with intact bulls. Indeed, a list of genes of MHC class II family members, including BLA-DQB and BOLA-DQB, have been associated with reproductive performance in cattle (Pinto De Melo et al., 2017).

      The lack of a response to mating with a vasectomized bull, together with the discovery of genes that were regulated by mating to intact bulls, led to the analysis of tissues that have direct contact with seminal fluid (i.e., the vagina). The ejaculation site during natural conception is an important factor to take into consideration since a recent study has shown that SP infusion into the vagina, but not into the uterus, could influence the levels of endometrial epidermal growth factor (Badrakh et al., 2020), which has been associated with fertility restoration in repeat breeder cows (Katagiri and Takahashi, 2006). In this sense, we were also interested in studying how far into the reproductive tract these changes could be observed, so gene expression analysis in the oviductal ampulla was also undertaken. Based on the sequencing results, genes identified as possible key regulators of uterine environment and pregnancy success (PLA2G10, CX3CL1, C4BPA, PRSS2, BLA-DQB, and CEBPD) were selected. In addition, the expression of inflammatory mediators regulated by SP components in other species (IL1A, IL6, TNFA, and IL8; O’Leary et al., 2004; Sharkey et al., 2012b; Schjenken et al., 2015; Introini et al., 2017; Ibrahim et al., 2019) was also analyzed. However, no difference between treatments was observed in the expression of any gene in the vagina or the oviduct. These results are not consistent with data in other species: in the mouse, Il6 was reduced in the oviduct of females mated to males that had undergone a vesiculectomy with or without a vasectomy (Bromfield et al., 2014); while in the pig, natural mating induces up-regulation of CEBPD in the ampulla (Alvarez-Rodriguez et al., 2019); finally, in human cervix, an increase in IL1A, IL6 and IL8 is observed after coitus (Sharkey et al., 2012b). It is important to highlight, however, that the heifers in the present study had been estrous synchronized (both mated and control animals) with an intravaginal device that was removed 48 h prior to sample collection. Although all the animals were managed under the same conditions, the resulting manipulation could have had an impact on the inflammatory status of the vagina. Indeed, a higher dispersion within group is observed in the expression of inflammatory genes in the vagina, than in the ampulla region. This is especially noticeable in the intact group, were one heifer exhibits very high expression of CX3CL1, IL6, IL1A, TNFA, and IL8. Despite the lack of treatment effect, some genes were shown to be differentially expressed between tissues. A higher expression of CX3CL1 was found in the ampulla compared with the vagina in control samples. In women, CXCL3 is present throughout the oviduct and, interestingly, its receptor was found in ejaculated sperm (Zhang et al., 2004). Conversely, CEBPD and IL8 were found up-regulated in the vagina samples compared with ampulla in all treatments. On the one hand, the expression of CEBPD is a crucial factor during inflammatory acute-phase response, under regulation of a range of cytokines and other inflammatory agents (reviewed in Ramji and Foka, 2002) and, on the other hand, IL8 is well known to be a potent neutrophil chemoattractant (Leonard and Yoshimura, 1990). Therefore, a higher basal expression of both genes in the vagina compared with the ampulla might be expected since this tissue has contact with the outside and, thus, is more prone to environmental/external contaminants.

      Interestingly, IL1A and TNFA were more highly expressed in vagina than in the ampulla tissue of heifers mated either to intact or vasectomized bulls, while control heifers did not exhibit this region-specific difference, suggesting a modulatory role induced by mating. In human exposure of ectocervical explants to SP resulted in the increase of IL1A and TNFA expression levels (Introini et al., 2017). Further, unprotected vaginal coitus, but not condom protected, induced the expression of IL1A in women (Sharkey et al., 2012b). However, in our study, the mechanical stimulus of mating cannot be ruled out.

      Conclusion

      The lack of changes in the endometrial transcriptome and in the expression of selected genes in the vagina and oviduct after mating to a vasectomized male do not support a role of SP (in the absence of sperm nor testicular and epididymal secretions) in regulating early pregnancy and uterine environment in cattle. Rather, the subtle changes in the transcriptome of the endometrium and the vagina seem to be elicited by sperm. These data indicate that, in species that ejaculate intravaginally, sperm play a more critical role in the modulation of the female environment. This is most apparent when looking at regulation of IL17, which is driven by SP in mice (Song et al., 2016) and by sperm in cattle. However, further research is needed to elucidate the role of inherent sperm proteins or SP proteins that attach to sperm at ejaculation.

      Data Availability Statement

      Gene expression data are publicly available at the Dryad Digital Repository (https://doi.org/10.5061/dryad.s7h44j14r).

      Ethics Statement

      The animal study was reviewed and approved by the Animal Research Ethics Committee of University College Dublin, Ireland and licensed by the Health Products Regulatory Authority (HPRA), Ireland, in accordance with Statutory Instrument No. 543 of 2012 (under Directive 2010/63/EU on the Protection of Animals used for Scientific Purposes).

      Author Contributions

      SR carried out the laboratory work, analyzed the results, and wrote the draft. YM-O contributed to the laboratory work and the analysis of results. BF-F, JS, SB-A, MM, and DK carried out the animal work, including handling of bulls and heifers, estrus detection, mating and sample collection. SB and TS performed the RNA library preparation and sequencing analysis. MY performed the statistical analysis of the data. BF-F created the illustrations for Figure 1. PL, JS, MY, and BF-F contributed to the critical revision of the manuscript. BF-F, PL, and JS designed the study. All authors read and approved the final manuscript.

      Conflict of Interest

      The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

      Funding. This work was funded by the EU Horizon 2020 Marie Skłodowska-Curie (No. 792212) and Science Foundation Ireland (Grant No. 16/IA/4474). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

      The authors thank the farm staff at Teagasc Grange for their important collaboration and patience. The authors acknowledge the excellent implication of several colleagues for their assistance during sample collection and processing, especially Mary Wade, Dr. Clio Maicas, and Dr. John A. Browne. The authors also thank Dr. Isabel Barranco, Ariadna Delgado-Bermúdez, and Marc Llavanera for their help during data analysis and elaboration of some figures.

      Supplementary Material

      The Supplementary Material for this article can be found online at: /articles/10.3389/fcell.2020.00547/full#supplementary-material

      References Alexander F. C. M. Zemjanis R. Graham E. F. Schmehl M. L. (1971). Semen characteristics and chemistry from bulls before and after seminal vesiculectomy and after vasectomy. J. Dairy Sci. 54 15301535. 10.3168/jds.S0022-0302(71)86059-9 Alghamdi A. S. Foster D. N. Troedsson M. H. T. (2004). Equine seminal plasma reduces sperm binding to polymorphonuclear neurophils (PMN’s) and improves the fertility of fresh semen inseminated into inflamed uteri. Reproduction 127 593600. 10.1530/rep.1.00096 15129015 Alvarez-Rodriguez M. Atikuzzaman M. Venhoranta H. Wright D. Rodriguez-Martinez H. (2019). Expression of immune regulatory genes in the porcine internal genital tract is differentially triggered by spermatozoa and seminal plasma. Int. J. Mol. Sci. 20 120. 10.3390/ijms20030513 30691059 Amann R. P. Kavanaugh J. F. Griel L. C. Voglmayr J. K. (1974). Sperm production of holstein bulls determined from testicular spermatid reserves, after cannulation of rete testis or vas deferens, and by daily ejaculation. J. Dairy Sci. 57 9399. 10.3168/jds.S0022-0302(74)84836-8 Arai M. Yoshioka S. Tasaki Y. Okuda K. (2013). Remodeling of bovine endometrium throughout the estrous cycle. Anim. Reprod. Sci. 142 19. 10.1016/j.anireprosci.2013.08.003 24051170 Badrakh D. Yanagawa Y. Nagano M. Katagiri S. (2020). Effect of seminal plasma infusion into the vagina on the normalization of endometrial epidermal growth factor concentrations and fertility in repeat breeder dairy cows. J. Reprod. Dev. 66 149154. 10.1262/jrd.2019-148 31996487 Bae S. E. Watson E. D. (2003). A light microscopic and ultrastructural study on the presence and location of oxytocin in the equine endometrium. Theriogenology 60 909921. 10.1016/S0093-691X(02)01362-6 Bai R. Bai H. Kuse M. Ideta A. Aoyagi Y. Fujiwara H. (2014). Involvement of VCAM1 in the bovine conceptus adhesion to the uterine endometrium. Reproduction 148 119127. 10.1530/REP-13-0655 24803492 Batruch I. Lecker I. Kagedan D. Smith C. R. Mullen B. J. Grober E. (2011). Proteomic analysis of seminal plasma from normal volunteers and post-vasectomy patients identifies over 2000 proteins and candidate biomarkers of the urogenital system. J. Proteome Res. 10 941953. 10.1021/pr100745u 21142078 Blom A. M. Villoutreix B. O. Dahlbäck B. (2004). Complement inhibitor C4b-binding protein – Friend or foe in the innate immune system? Mol. Immunol. 40 13331346. 10.1016/j.molimm.2003.12.002 15072852 Bromfeld J. J. (2016). A role for Seminal plasma in modulating pregnancy outcomes in domestic species. Reproduction 152 R223R232. 10.1530/REP-16-0313 27601714 Bromfield J. J. Schjenken J. E. Chin P. Y. Care A. S. Jasper M. J. Robertson S. A. (2014). Maternal tract factors contribute to paternal seminal fluid impact on metabolic phenotype in offspring. Proc. Natl. Acad. Sci. U.S.A. 111 22002205. 10.1073/pnas.1305609111 24469827 Byrne K. Leahy T. Mcculloch R. Colgrave M. L. Holland M. K. (2012). Comprehensive mapping of the bull sperm surface proteome. Proteomics 12 35593579. 10.1002/pmic.201200133 23081703 Capper E. A. Marshall L. A. (2001). Mammalian phospholipases A2: mediators of inflammation, proliferation and apoptosis. Prog. Lipid Res. 40 167197. 10.1016/S0163-7827(01)00002-9 Cooper T. G. Noonan E. von Eckardstein S. Auger J. Baker H. W. G. Behre H. M. (2009). World health organization reference values for human semen characteristics. Hum. Reprod. Update 16 231245. 10.1093/humupd/dmp048 19934213 Dean M. D. Findlay G. D. Hoopmann M. R. Wu C. C. MacCoss M. J. Swanson W. J. (2011). Identification of ejaculated proteins in the house mouse (Mus domesticus) via isotopic labeling. BMC Genomics 12:306. 10.1186/1471-2164-12-306 21663664 Dickinson S. E. Griffin B. A. Elmore M. F. Kriese-Anderson L. Elmore J. B. Dyce P. W. (2018). Transcriptome profiles in peripheral white blood cells at the time of artificial insemination discriminate beef heifers with different fertility potential. BMC Genomics 19:129. 10.1186/s12864-018-4505-4 29426285 Diskin M. G. Morris D. G. (2008). Embryonic and early foetal losses in cattle and other ruminants. Reprod. Domest. Anim. 43 260267. 10.1111/j.1439-0531.2008.01171.x 18638133 Elweza A. E. Ezz M. A. Acosta T. J. Talukder A. K. Shimizu T. Hayakawa H. (2018). A proinflammatory response of bovine endometrial epithelial cells to active sperm in vitro. Mol. Reprod. Dev. 85 215226. 10.1002/mrd.22955 29337420 Ezz M. A. Marey M. A. Elweza A. E. Kawai T. Heppelmann M. Pfarrer C. (2019). TLR2/4 signaling pathway mediates sperm-induced inflammation in bovine endometrial epithelial cells in vitro. PLoS ONE 14:e0214516. 10.1371/journal.pone.0214516 30995239 Faulkner L. C. Hopwood M. L. Wiltbank J. N. (1968). Seminal vesiculectomy in bulls. II. Seminal characteristics and breeding trials. J. Reprod. Fertil. 16 179182. 10.1530/jrf.0.0160179 5660277 Fedorka C. E. Scoggin K. E. Woodward E. M. Squires E. L. Ball B. A. Troedsson M. H. T. (2017). The effect of select seminal plasma proteins on endometrial mRNA cytokine expression in mares susceptible to persistent mating-induced endometritis. Reprod. Domest. Anim. 52 8996. 10.1111/rda.12813 27686063 Fernandez-fuertes B. Sánchez J. M. Bagés-arnal S. Mcdonald M. Yeste M. Lonergan P. (2019). Species-specific and collection method-dependent differences in endometrial susceptibility to seminal plasma-induced RNA degradation. Sci. Rep. 9 112. 10.1038/s41598-019-51413-4 31636362 Gangnuss S. Sutton-McDowall M. L. Robertson S. A. Armstrong D. T. (2004). Seminal plasma regulates corpora lutea macrophage populations during early pregnancy in mice. Biol. Reprod. 71 11351141. 10.1095/biolreprod.104.027425 15175232 Glynn D. J. Heng K. Russell D. L. Sharkey D. J. Robertson S. A. Anand-Ivell R. (2017). Male seminal relaxin contributes to induction of the post-mating cytokine response in the female mouse uterus. Front. Physiol. 8:422. 10.3389/fphys.2017.00422 28674503 Gutsche S. von Wolff M. Strowitzki T. Thaler C. J. (2003). Seminal plasma induces mRNA expression of IL-1β, IL-6 and LIF in endometrial epithelial cells in vitro. Mol. Hum. Reprod. 9 785791. 10.1093/molehr/gag095 14614040 Gwathmey T. M. Ignotz G. G. Mueller J. L. Manjunath P. Suarez S. S. (2006). Bovine seminal plasma proteins PDC-109, BSP-A3, and BSP-30-kDa share functional roles in storing sperm in the oviduct. Biol. Reprod. 75 501507. 10.1095/biolreprod.106.053306 16790686 Gwathmey T. M. Ignotz G. G. Suarez S. S. (2003). PDC-109 (BSP-A1/A2) promotes bull sperm binding to oviductal epithelium in vitro and may be involved in forming the oviductal sperm reservoir. Biol. Reprod. 69 809815. 10.1095/biolreprod.102.010827 12748117 Hannan N. J. Jones R. L. Critchley H. O. D. Kovacs G. J. Rogers P. A. W. Affandi B. (2004). Coexpression of fractalkine and its receptor in normal human endometrium and in endometrium from users of progestin-only contraception supports a role for fractalkine in leukocyte recruitment and endometrial remodeling. J. Clin. Endocrinol. Metab. 89 61196129. 10.1210/jc.2003-031379 15579768 Hannan N. J. Jones R. L. White C. A. Salamonsen L. A. (2006). The chemokines, CX3CL1, CCL14, and CCL4, promote human trophoblast migration at the feto-maternal interface. Biol. Reprod. 74 896904. 10.1095/biolreprod.105.045518 16452465 Harris C. L. Mizuno M. Morgan B. P. (2006). Complement and complement regulators in the male reproductive system. Mol. Immunol. 43 5767. 10.1016/j.molimm.2005.06.026 16023728 Hayashi K. G. Hosoe M. Kizaki K. Fujii S. Kanahara H. Takahashi T. (2017). Differential gene expression profiling of endometrium during the mid-luteal phase of the estrous cycle between a repeat breeder (RB) and non-RB cows. Reprod. Biol. Endocrinol. 15:20. 10.1186/s12958-017-0237-6 28335821 Hunter R. H. F. (1981). Sperm transport and reservoirs in the pig oviduct in relation to the time of ovulation. J. Reprod. Fertil. 63 109117. 10.1530/jrf.0.0630109 6895091 Ibrahim L. A. Rizo J. A. Fontes P. L. P. Lamb G. C. Bromfield J. J. (2019). Seminal plasma modulates expression of endometrial inflammatory meditators in the bovine. Biol. Reprod. 100 660671. 10.1093/biolre/ioy226 30329018 Introini A. Boström S. Bradley F. Gibbs A. Glaessgen A. Tjernlund A. (2017). Seminal plasma induces inflammation and enhances HIV-1 replication in human cervical tissue explants. PLoS Pathog. 13:e1006492. 10.1371/journal.ppat.1006492 28700681 Iwakura Y. Ishigame H. Saijo S. Nakae S. (2011). Functional specialization of interleukin-17 family members. Immunity 34 149162. 10.1016/j.immuni.2011.02.012 21349428 Janeway C. A. Jr. Travers P. Walport M. Shlomchik M. (2001). Immunobiology: The immune System in Health and Disease, 5th Edn. New York, NY: Garland Science. Jiang Y. Oliver P. Davies K. E. Platt N. (2006). Identification and characterization of murine SCARA5, a novel class A scavenger receptor that is expressed by populations of epithelial cells. J. Biol. Chem. 281 1183411845. 10.1074/jbc.M507599200 16407294 Katagiri S. Takahashi Y. (2006). Potential relationship between normalization of endometrial epidermal growth factor profile and restoration of fertility in repeat breeder cows. Anim. Reprod. Sci. 95 5466. 10.1016/j.anireprosci.2005.09.001 16207515 Kelleher A. M. DeMayo F. J. Spencer T. E. (2019). Uterine glands: developmental biology and functional roles in pregnancy. Endocr. Rev. 40 14241445. 10.1210/er.2018-00281 31074826 Kim D. Langmead B. Salzberg S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12 357360. 10.1038/nmeth.3317 25751142 Konac E. Alp E. Onen H. I. Korucuoglu U. Biri A. A. Menevse S. (2009). Endometrial mRNA expression of matrix metalloproteinases, their tissue inhibitors and cell adhesion molecules in unexplained infertility and implantation failure patients. Reprod. Biomed. Online 19 391397. 10.1016/s1472-6483(10)60174-5 Lee J. Oh J. Choi E. Park I. Han C. Kim D. H. (2007). Differentially expressed genes implicated in unexplained recurrent spontaneous abortion. Int. J. Biochem. Cell Biol. 39 22652277. 10.1016/j.biocel.2007.06.012 17681868 Lemaster J. W. Seals R. C. Hopkins F. M. Schrick F. N. (1999). Effects of administration of oxytocin on embryonic survival in progestogen supplemented cattle. Prostaglandins Other Lipid Mediat. 57 259268. 10.1016/S0090-6980(99)00007-6 Leonard E. J. Yoshimura T. (1990). Neutrophil attractant/activation protein-l (NAP-l [Interleukin-8]). Am. J. Respir. Cell Mol. Biol. 2 479486. 10.1165/ajrcmb/2.6.479 2189453 Liao Y. Smyth G. K. Shi W. (2014). FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30 923930. 10.1093/bioinformatics/btt656 24227677 Livak K. J. Schmittgen T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25 402408. 10.1006/meth.2001.1262 11846609 Macklon N. S. Geraedts J. P. M. Fauser B. C. J. M. (2002). Conception to ongoing pregnancy: the “black box” of early pregnancy loss. Hum. Reprod. Update 8 333343. 10.1093/humupd/8.4.333 12206468 Manjunath P. Thérien I. (2002). Role of seminal plasma phospholipid-binding proteins in sperm membrane lipid modification that occurs during capacitation. J. Reprod. Immunol. 53 109119. 10.1016/S0165-0378(01)00098-5 Mansouri-Attia N. Aubert J. Reinaud P. Giraud-Delville C. Taghouti G. Galio L. (2009). Gene expression profiles of bovine caruncular and intercaruncular endometrium at implantation. Physiol. Genomics 39 1427. 10.1152/physiolgenomics.90404.2008 19622795 Mantena S. R. Kannan A. Cheon Y. P. Li Q. Johnson P. F. Bagchi I. C. (2006). C/EBPβ is a critical mediator of steroid hormone-regulated cell proliferation and differentiation in the uterine epithelium and stroma. Proc. Natl. Acad. Sci. U.S.A. 103 18701875. 10.1073/pnas.0507261103 16439483 Mateo-Otero Y. Sánchez J. M. Recuero S. Bagés-Arnal S. Mcdonald M. Kenny D. (2020). Effect of exposure to seminal plasma through natural mating in cattle on conceptus length and gene expression. Front. Cell Dev. Biol. 8:341. 10.3389/fcell.2020.00341 32478076 Minten M. A. Bilby T. R. Bruno R. G. S. Allen C. C. Madsen C. A. Wang Z. (2013). Effects of fertility on gene expression and function of the bovine endometrium. PLoS ONE 8:e067444. 10.1371/journal.pone.0069444 23940519 Moraes J. G. N. Behura S. K. Geary T. W. Hansen P. J. Neibergs H. L. Spencer T. E. (2018). Uterine influences on conceptus development in fertility-classified animals. Proc. Natl. Acad. Sci. U.S.A. 115 E1749E1758. 10.1073/pnas.1721191115 29432175 Morgan H. L. Watkins A. J. (2020). The influence of seminal plasma on offspring development and health. Semin. Cell Dev. Biol. 97 131137. 10.1016/j.semcdb.2019.06.008 31254609 Nauc V. Manjunath P. (2000). Radioimmunoassays for bull seminal plasma proteins (BSP-A1/-A2, BSP-A3, and BSP-30-Kilodaltons), and their quantification in seminal plasma and sperm1. Biol. Reprod. 63 10581066. 10.1095/biolreprod63.4.1058 10993827 Neupane M. Geary T. W. Kiser J. N. Burns G. W. Hansen P. J. Spencer T. E. (2017). Loci and pathways associated with uterine capacity for pregnancy and fertility in beef cattle. PLoS ONE 12:e0188997. 10.1371/journal.pone.0188997 29228019 Niakan K. K. Han J. Pedersen R. A. Simon C. Pera R. A. R. (2012). Human pre-implantation embryo development. Development 139 829841. 10.1242/dev.060426 22318624 Odhiambo J. F. Poole D. H. Hughes L. DeJarnette J. M. Inskeep E. K. Dailey R. A. (2009). Pregnancy outcome in dairy and beef cattle after artificial insemination and treatment with seminal plasma or transforming growth factor beta-1. Theriogenology 72 566571. 10.1016/j.theriogenology.2009.04.013 19501392 O’Leary S. Jasper M. J. Robertson S. A. Armstrong D. T. (2006). Seminal plasma regulates ovarian progesterone production, leukocyte recruitment and follicular cell responses in the pig. Reproduction 132 147158. 10.1530/rep.1.01119 16816340 O’Leary S. Jasper M. J. Warnes G. M. Armstrong D. T. Robertson S. A. (2004). Seminal plasma regulates endometrial cytokine expression, leukocyte recruitment and embryo development in the pig. Reproduction 128 237247. 10.1530/rep.1.00160 15280563 Ortiz W. G. Rizo J. A. Carvalheira L. R. Ahmed B. M. S. Estrada-Cortes E. Harstine B. R. (2019). Effects of intrauterine infusion of seminal plasma at artificial insemination on fertility of lactating Holstein cows. J. Dairy Sci. 102 65876594. 10.3168/jds.2019-16251 31103294 Palm F. Walter I. Budik S. Kolodziejek J. Nowotny N. Aurich C. (2008). Influence of different semen extenders and seminal plasma on PMN migration and on expression of IL-1β, IL-6, TNF-α and COX-2 mRNA in the equine endometrium. Theriogenology 70 843851. 10.1016/j.theriogenology.2008.04.054 18584861 Pandya I. J. Cohen J. (1985). The leukocytic reaction of the human uterine cervix to spermatozoa. Fertil. Steril. 43 417421. 10.1016/S0015-0282(16)48442-6 Pfeiffer K. E. Binversie J. A. Rhinehart J. D. Larson J. E. (2012). Exposure of beef females to the biostimulatory effects of bulls with or without deposition of seminal plasma prior to AI. Anim. Reprod. Sci. 133 2734. 10.1016/j.anireprosci.2012.06.011 22763235 Pini T. Leahy T. Soleilhavoup C. Tsikis G. Labas V. Combes-Soia L. (2016). Proteomic investigation of ram spermatozoa and the proteins conferred by seminal plasma. J. Proteome Res. 15 37003711. 10.1021/acs.jproteome.6b00530 27636150 Pinto, De Melo T. De Camargo G. M. F. De Albuquerque L. G. Carvalheiro R. (2017). Genome-wide association study provides strong evidence of genes affecting the reproductive performance of Nellore beef cows. PLoS ONE 12:e0178551. 10.1371/journal.pone.0178551 28562680 Ramji D. P. Foka P. (2002). CCAAT/enhancer-binding proteins: structure, function and regulation. Biochem. J. 365 561575. 10.1042/BJ20020508 12006103 Randi F. McDonald M. Duffy P. Kelly A. K. Lonergan P. (2018). The relationship between external auditory canal temperature and onset of estrus and ovulation in beef heifers. Theriogenology 110 175181. 10.1016/j.theriogenology.2018.01.001 29407899 Remes Lenicov F. Rodriguez Rodrigues C. Sabatté J. Cabrini M. Jancic C. Ostrowski M. (2012). Semen promotes the differentiation of tolerogenic dendritic cells. J. Immunol. 189 47774786. 10.4049/jimmunol.1202089 23066152 Robertson S. A. (2007). Seminal fluid signaling in the female reproductive tract: lessons from rodents and pigs. J. Anim. Sci. 85 3644. 10.2527/jas.2006-578 17085725 Robertson S. A. Guerin L. R. Bromfield J. J. Branson K. M. Ahlström A. C. Care A. S. (2009). Seminal fluid drives expansion of the CD4+CD25+ T regulatory cell pool and induces tolerance to paternal alloantigens in mice1. Biol. Reprod. 80 10361045. 10.1095/biolreprod.108.074658 19164169 Sakaguchi S. Sakaguchi N. Shimizu J. Yamazaki S. Sakihama T. Itoh M. (2001). Immunologic tolerance maintained by CD25+ CD4+ regulatory T cells: their common role in controlling autoimmunity, tumor immunity, and transplantation tolerance. Immunol. Rev. 182 1832. 10.1034/j.1600-065X.2001.1820102.x 11722621 Sánchez J. M. Mathew D. J. Behura S. K. Passaro C. Charpigny G. Butler S. T. (2019). Bovine endometrium responds differentially to age-matched short and long conceptuses†. Biol. Reprod. 101 2639. 10.1093/biolre/ioz060 30977805 Sandra O. Charpigny G. Galio L. Hue I. (2017). Preattachment embryos of domestic animals: insights into development and paracrine secretions. Annu. Rev. Anim. Biosci. 5 205228. 10.1146/annurev-animal-022516-022900 27959670 Schjenken J. E. Glynn D. J. Sharkey D. J. Robertson S. A. (2015). TLR4 signaling is a major mediator of the female tract response to seminal fluid in mice1. Biol. Reprod. 93 113. 10.1095/biolreprod.114.125740 26157066 Seidel G. Foote R. (1970). Compartmental analysis bovine of sources of the ejaculate. Biol. Reprod. 2 189196. 10.1095/biolreprod2.2.189 5520322 Shah B. A. Hopwood M. L. Faulkner L. C. (1968). Seminal vesiculectomy in bulls. I. Seminal biochemistry. J. Reprod. Fertil. 16 171177. 10.1530/jrf.0.0160171 5660276 Sharkey D. J. Macpherson A. M. Tremellen K. P. Mottershead D. G. Gilchrist R. B. Robertson S. A. (2012a). TGF-β mediates proinflammatory seminal fluid signaling in human cervical epithelial cells. J. Immunol. 189 10241035. 10.4049/jimmunol.1200005 22706080 Sharkey D. J. Tremellen K. P. Jasper M. J. Gemzell-Danielsson K. Robertson S. A. (2012b). Seminal fluid induces leukocyte recruitment and cytokine and chemokine mRNA expression in the human cervix after coitus. J. Immunol. 188 24452454. 10.4049/jimmunol.1102736 22271649 Sharkey D. J. Macpherson A. M. Tremellen K. P. Robertson S. A. (2007). Seminal plasma differentially regulates inflammatory cytokine gene expression in human cervical and vaginal epithelial cells. Mol. Hum. Reprod. 13 491501. 10.1093/molehr/gam028 17483528 Shima T. Inada K. Nakashima A. Ushijima A. Ito M. Yoshino O. (2015). Paternal antigen-specific proliferating regulatory T cells are increased in uterine-draining lymph nodes just before implantation and in pregnant uterus just after implantation by seminal plasma-priming in allogeneic mouse pregnancy. J. Reprod. Immunol. 108 7282. 10.1016/j.jri.2015.02.005 25817463 Song Z. H. Li Z. Y. Li D. D. Fang W. N. Liu H. Y. Yang D. D. (2016). Seminal plasma induces inflammation in the uterus through the γδT/IL-17 pathway. Sci. Rep. 6 18. 10.1038/srep25118 27109934 Spencer T. E. Forde N. Dorniak P. Hansen T. R. Romero J. J. Lonergan P. (2013). Conceptus-derived prostaglandins regulate gene expression In the endometrium prior to pregnancy recognition In ruminants. Reproduction 146 377387. 10.1530/REP-13-0165 23966582 Spencer T. E. Kelleher A. M. Bartol F. F. (2019). Development and function of uterine glands in domestic animals. Annu. Rev. Anim. Biosci. 7 125147. 10.1146/annurev-animal-020518-115321 30183326 Tapia A. Gangi L. M. Zegers-Hochschild F. Balmaceda J. Pommer R. Trejo L. (2008). Differences in the endometrial transcript profile during the receptive period between women who were refractory to implantation and those who achieved pregnancy. Hum. Reprod. 23 340351. 10.1093/humrep/dem319 18077318 Tapia A. Vilos C. Marín J. C. Croxatto H. B. Devoto L. (2011). Bioinformatic detection of E47, E2F1 and SREBP1 transcription factors as potential regulators of genes associated to acquisition of endometrial receptivity. Reprod. Biol. Endocrinol. 9 114. 10.1186/1477-7827-9-14 21272326 Thompson L. A. Barratt C. L. R. Bolton A. E. Cooke I. D. (1992). The leukocytic reaction of the human uterine cervix. Am. J. Reprod. Immunol. 28 8589. 10.1111/j.1600-0897.1992.tb00765.x 1285856 Tríbulo P. Siqueira L. G. B. Oliveira L. J. Scheffler T. Hansen P. J. (2018). Identification of potential embryokines in the bovine reproductive tract. J. Dairy Sci. 101 690704. 10.3168/jds.2017-13221 29128220 Troedsson M. H. T. Alghamdi A. S. Mattisen J. (2002). Equine seminal plasma protects the fertility of spermatozoa in an inflamed uterine environment. Theriogenology 58 453456. 10.1016/S0093-691X(02)00862-2 Tunon A.-M. Katila T. Magnusson U. Nummijärvi A. Rodriguez-Martinez H. (2000). T-cell distribution in two different segments of the equine. Theriogenology 54 835841. 10.1016/s0093-691x(00)00395-2 Valenza A. Giordano J. O. Lopes G. Vincenti L. Amundson M. C. Fricke P. M. (2012). Assessment of an accelerometer system for detection of estrus and treatment with gonadotropin-releasing hormone at the time of insemination in lactating dairy cows. J. Dairy Sci. 95 71157127. 10.3168/jds.2012-5639 23040033 Vitorino Carvalho A. Eozenou C. Richard C. Forde N. Healey G. D. Giraud-Delville C. (2019). Bovine scavenger receptor class A (SR-A) exhibit specific patterns of regulation in the endometrium during the oestrous cycle and early pregnancy. Reprod. Fertil. Dev. 31 10781090. 10.1071/RD18411 30922439 Walker W. L. Nebel R. L. McGilliard M. L. (1996). Time of ovulation relative to mounting activity in dairy cattle. J. Dairy Sci. 79 15551561. 10.3168/jds.S0022-0302(96)76517-7 Wathes D. C. Cheng Z. Chowdhury W. Fenwick M. A. Fitzpatrick R. Morris D. G. (2009). Negative energy balance alters global gene expression and immune responses in the uterus of postpartum dairy cows. Physiol. Genomics 39 113. 10.1152/physiolgenomics.00064.2009 19567787 Watkins A. J. Dias I. Tsuro H. Allen D. Emes R. D. Moreton J. (2018). Paternal diet programs offspring health through sperm- and seminal plasma-specific pathways in mice. Proc. Natl. Acad. Sci. U.S.A. 115 1006410069. 10.1073/pnas.1806333115 30150380 Wedel A. Lömsziegler-Heitbrock H. W. (1995). The C/EBP family of transcription factors. Immunobiology 193 171185. 10.1016/S0171-2985(11)80541-3 Wiltbank M. C. Baez G. M. Garcia-Guerra A. Toledo M. Z. Monteiro P. L. J. Melo L. F. (2016). Pivotal periods for pregnancy loss during the first trimester of gestation in lactating dairy cows. Theriogenology 86 239253. 10.1016/j.theriogenology.2016.04.037 27238438 Zhang Q. Shimoya K. Temma K. Kimura T. Tsujie T. Shioji M. (2004). Expression of fractalkine in the Fallopian tube and of CX3CR1 in sperm. Hum. Reprod. 19 409414. 10.1093/humrep/deh085 14747189 Zhou X. Lindsay H. Robinson M. D. (2014). Robustly detecting differential expression in RNA sequencing data using observation weights. Nucleic Acids Res. 42:e91. 10.1093/nar/gku310 24753412

      https://ccb.jhu.edu/software/fqtrim/

      pantherdb.org

      https://www.ncbi.nlm.nih.gov/tools/primer-blast/

      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016www.kjchain.com.cn
      www.jxmmzz.org.cn
      www.jhwywx.com.cn
      geekkg.com.cn
      www.siqzsq.com.cn
      www.suyin.net.cn
      www.tiyanjia.com.cn
      www.tashout.org.cn
      www.wzszyz.org.cn
      wltgsn.com.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p