Novel endogenous retrovirus-derived transcript expressed in the bovine placenta is regulated by WNT signaling.

Endogenous retroviruses (ERVs) are involved in placentation; perhaps, the most well-known ERVs are the syncytins, actively transcribed env genes involved in cell-cell fusion and possible morphological variations. However, ERVs other than syncytins that play an important role in placental development have not been well characterized. To identify ERV genes expressed during the onset of placentation in the bovine species, we characterized the expression profiles of bovine conceptus transcripts during the peri-attachment period using RNA-seq analysis, and confirming some candidates through real-time PCR. Using in silico and PCR analyses, we identified a novel ERV proviral sequence derived from a gag region, designated bovine endogenous retroviruses (BERV)-K3, containing Gag_p10 and Gag_p24, zinc finger domain. Initial expression of this ERV in bovine conceptuses was on day 20 (day 0 = day of estrus), soon after conceptus attachment to the endometrial epithelium, and its high placental expression was maintained up to the middle of pregnancy. The BERV-K3 transcript was also found in the uterine luminal and glandular epithelia, liver, kidney, intestine, and skin. BERV-K3 is located on chromosome 7 and integrated within LOC100848658, from which noncoding RNA could be transcribed. Furthermore, the expression of endogenous BERV-K3 in bovine trophoblast cell lines was induced by a WNT agonist, a signaling system common to genes expressed in placentas. These data support the argument that during the evolutionary process, mammals incorporated not only similar ERV sequences, but also ERVs unique to individual species. BERV-K3 is in the latter case, likely providing functions unique to ruminant gestation.

Introduction

The endogenous retroviruses (ERVs) are thought to be derived from ancient viral infections of germ cells or their precursors with a retrovirus in mammals and other vertebrates [1]. Indeed, sequences of retroviral origins make up 8 % of the human genome, 10% of the mouse genome, and 18% of the bovine genome [2-4]. Most ERVs have been inactivated by insertions, deletions, and substitutions of nucleotides, and/or epigenetic modifications such as DNA methylation and histone modifications. However, a few open-reading frames (ORFs) of ERV sequences remain intact and are expressed as virus-derived proteins in the host cells [5].

The placenta is remarkably distinct among mammalian species, suggesting a history of rapid evolutionary diversification, resulting from the genes acquired in individual species. It has become apparent that ERV genes play an important role in the development of the placenta and the trophoblast cell lineage in mammalian species, and that during evolution different species may have utilized ERVs of the same as well as different origins. Indeed, different env genes, syncytins, essential for placental morphogenesis have been independently integrated into the genome of humans [6-10], mice [11], rabbits [12], dogs [13], cats [13], sheep [14-18], and cattle [19-23], sheep [21-23], the Rodentia squirrel-related clade [24], Afrotherian tenrecs [25], and marsupials [26]. All identified syncytin env genes exist in different genome sequences and chromosomal locations among species, but the functions such as cell fusion and immune suppression are all shared in mammals. However, the exact evolutionary pathways and the extent to which ERVs function in placental development are still unclear.

It has been determined that the WNT signaling pathway is an important regulator of embryo/conceptus and maternal interaction such as implantation and placental development in mice, sheep, cow, and humans [27-30]. The WNT can induce two downstream signaling cascades, known as the canonical and noncanonical pathways [31]. The canonical WNT pathway is activated when WNT binds to Frizzled (FZD) family receptors and low-density lipoprotein receptor-related protein (LRP) 5/6 coreceptors. This, in turn, leads to cellular accumulation and nuclear translocation of β-catenin (CTNNB), followed by a complex formation with T-cell fate (TCF) to activate transcription of WNT-target genes.

It was shown in mice that one of the Wnt receptors, Fzd5, was necessary for the expression of glial cells missing homolog 1 (Gcm1), an important transcription factor for placental labyrinth development, followed by induction of syncytin-A expression [32]. Matsuura et al. [33] also showed that Wnt/Ctnnb signaling induces Gcm1-syncytin-A, the cell fusion pathway. Although the WNT signal is an important regulator for the expression of ERV during the period of invasive placental formation, whether the WNT signal induces ERV expression in the noninvasive bovine placenta has not been elucidated.

Unlike primate or murine species, conceptus attachment to the uterine endometrial epithelium and subsequent placentation in most ruminants do not occur soon after blastocyst formation [34]. In fact, the conceptus spends a prolonged period within the uterine lumen before developing a definitive attachment to the endometrial epithelium and subsequent formation of placental structures [35]. In the bovine species, ERVs such as BERV-K1 [20], bERVE-A [21], and syncytin-Rum1 [22] have been identified and their potential functions studied. It should be noted, however, that these ERVs are all from env regions; ERVs from other regions such as gag or pol, which could also function in ruminants' placental development, have not been identified or characterized. We hypothesized that ERVs of regions other than env would exist and function in the trophectoderm during the period of placental formation and functioning.

We searched for nucleotide structures of ERV origin, which were expressed in bovine conceptuses during the peri-attachment periods (the criteria are shown in Supplementary Figure S1). Using RNA-seq data, we found that one candidate gene with gag/pol, located on chromosome 7, was minimally expressed on day 20 (day 0 = day of estrus), when conceptus implantation to the uterine epithelium is initiated. Its expression was found to be elevated on day 22, when the placental formation is initiated, and its high expression continued until the middle of pregnancy. Expression of endogenous ERV gene in bovine trophoblast cells was induced by a WNT agonist, a common intracellular signaling for genes expressed in placentas.

Materials and methods

Animals and sampling

All animal procedures in the present study were approved by the Committee for Experimental Animals at Zen-noh Embryo Transfer (ET) Center, Hokkaido and the University of Tokyo, Tokyo, Japan. Estrous synchronization, super-ovulation, artificial insemination and ET processes were performed as previously described [36]. Day 7 blastocysts were collected from super-ovulated and artificially inseminated Japanese black cattle. Sixteen blastocysts derived from the super-ovulation were transferred nonsurgically into the uterine horn of eight estrous synchronized Holstein heifers (n = 2 blastocysts per transfer), ipsilateral to the CL on day 7. Elongated conceptuses were then collected nonsurgically by uterine flushing (500 ml PBS) on day 17 (n = 2), 20 (n = 2), or 22 (n = 4) [37], corresponding to the period before conceptus attachment to the uterine epithelium, when conceptus attachment is initiated, or 2 days into conceptus attachment, respectively. Days 17 and 20 conceptuses (n = 4 each day) in the uterine flushing media were obtained by centrifugation at 400× for 5 min, snap-frozen and transferred to Animal Resource Science Center at the University of Tokyo. The remaining day 22 pregnant heifers (n = 2) were killed and hysterectomized, from which uteri were excised and frozen immediately using the OCT compound (Sakura Finetek, Tokyo, Japan).

Hysterectomy was performed to collect days 45 and 150 pregnant uteri (n = 3 each) from Japanese black cattle at Azabu University Veterinary Teaching Hospital (AUVTH), Kanagawa, Japan and the National Institute of Agrobiological Sciences (NIAS), Ibaraki, Japan, respectively, from which cotyledons and caruncles were dissected, immediately frozen and stored at −80°C until use. At NIAS, various tissues, heart, liver, kidney, intestine, lung, muscle, skin, lymph node, and spleen were harvested from day 150 pregnant Japanese black cattle, and the uteri were removed from three nonpregnant Japanese black cattle. All tissues dissected were frozen immediately and stored at −80°C until use.

RNA extraction from bovine conceptus and whole-genome sequencing by the RNA-seq (SOLiD3) system

RNA extraction from whole conceptus tissues was performed using Isogen (Nippon Gene, Tokyo, Japan) according to the manufacturer's protocol [38]. For a RNA-seq analysis, total RNA was depleted of rRNAs using the RiboMinus Eukaryote Kit (Life Technologies, Carlsbad, CA, U.S.A.). Details of the sequencing procedure using a SOLiD3 sequencer (Life Technologies), and data analysis have been described previously [37]. In brief, high-throughput sequencing libraries were prepared according to the SOLiD3 whole transcriptome library preparation protocol [39]. ERV-derived sequences in the bovine genome (Btau_4.0.55.dna.toplevel.fa) were identified through the use of RetroTector, which was designed to detect and characterize entire or fragmented ERVs in a given genome [40]. Four types of retroviral-like sequences identified in the genome are gag, pro, pol, and env. Processed nucleotide sequences from the RNA-seq for each of the 3 days examined were aligned against the bovine genome. The Applied Biosystems Whole Transcriptome Analysis Pipeline was used to map short reads. During this mapping phase, up to two mismatches were allowed and short reads that aligned to more than 10 locations were removed. In this analysis, we used the reads whose sequence quality scores were 24 or higher, following the standard parameters of the Applied Biosystems Whole Transcriptome Analysis Pipeline [39]. Matching locations were subsequently used to generate counts for identified ERVs and Ensemble-provided coding sequences (Bos_taurus.Btau_4.0.55.dna.toplevel.fa). To evaluate gene expression level independent of variance in gene lengths and the number of reads among samples, the method of reads per kilobase of exon model per million mapped reads (RPKM), a widely recognized quantification measurement, was applied [41]. The RNA-seq data can be downloaded from the DDBJ Sequence Read Archive [42] with accession number DRA000549 [37].

Cell culture

All cells were maintained at 38.5°C in humidified 5% CO2 and cultured up to confluence. Bovine trophoblast BT-1 [43], CT-1 [44], and F3 [45] cells were kindly provided by Prof. K. Hashizume (Iwate University, Japan), Prof. A. Ealy (Virginia Polytech Institute, U.S.A.), and Prof. C. Pfarrer (University of Veterinary Medicine, Hannover, Germany), respectively. BT-1 cells were cultured on plastic plates coated with type I collagen (Nitta Gelatin, Osaka, Japan) in Dulbecco's modified Eagle medium (DMEM)/Ham's F12 medium (F12) (Invitrogen, Carlsbad, CA, U.S.A.) supplemented with 10% (v/v) fetal bovine serum (FBS; JRH Biosciences, Lenexa, KS, U.S.A.) and antibiotic/antimycotic solution (Invitrogen). CT-1 cells were maintained in DMEM (Invitrogen) containing 10% (v/v) FBS (JRH Biosciences) supplemented with 4.5 g/l d-glucose, nonessential amino acids, 2 mM glutamine, 2 mM sodium pyruvate, 55 mM β-mercaptoethanol, and antibiotic/antimycotic solution (Invitrogen). F3 cells were cultured in DMEM/F12 medium supplemented with 10% FBS (JRH Biosciences) and antibiotic/antimycotic solution (Invitrogen). Bovine primary uterine endometrial epithelial cells (EECs) and stromal cells (STRs) were kindly provided by Prof. K. Okuda (Obihiro University, Japan) and were cultured in DMEM/F12 medium (Invitrogen) supplemented with 5% FBS (JRH Biosciences) and antibiotic/antimycotic solution. Bovine intestinal epithelial cells (Bie) were kindly provided by Prof. H. Aso (Tohoku University, Japan) and were maintained in DMEM/F12 (Invitrogen) supplemented with 10% FBS (JRH Biosciences) and antibiotic/antimycotic solution. Bovine ovarian cumulus-granulosa (oCG) cells were obtained from ovarian follicles that had been collected at a local abattoir, and ear-derived fibroblast (EF) cells were obtained from biopsied ear skin of Japanese black bulls (4 months old) [38], both cells were cultured in DMEM (Invitrogen) containing 5% (v/v) FBS (JRH Biosciences) and antibiotic/antimycotic solution. Bovine kidney cell lines (MDBK and CKT-1) and bovine macrophage cells (BoMAC) were cultured in DMEM (Invitrogen) supplemented with 10% v/v FBS (JRH Biosciences) and antibiotic/antimycotic solution. To study the effects of trophoblast attachment to the uterine endometrial epithelial cells, CT-1 cells were cultured without or with a cell culture insert (Falcon, BD Biosciences, Tokyo, Japan), allowing direct CT-1 cell contact to EECs or indirect cell association with EECs, respectively. To further characterize whether any of the candidate ERV genes could be regulated by Wnt signaling, cultured CT-1 or F3 cells were treated with 1 µM Wnt agonist (sc-222416, Santa Cruz Biotechnology, Dallas, TX, U.S.A.) for 24 h.

RNA isolation from bovine tissues and cultured cells

RNA isolation from bovine tissues and cultured cells was performed using the ISOGEN protocol (Nippon Gene), as described previously [38]. Bovine tissues, heart, liver, kidney, intestine, lung, muscle, skin, lymph node, spleen, and uterus were harvested from three Japanese black cattle at NIAS, Ibaraki, Japan. Excised tissues were submerged in RNAlater (Qiagen, Tokyo, Japan) to prevent RNA degradation, and RNA was then extracted from each tissue. RNA was also isolated from bovine cell lines, including trophoblast cell lines (BT-1, CT-1, and F3), EEC, STR, CKT-1, MDBK, Bie, EF, oCG, and BoMAC. Extracted RNAs were then stored at −30°C until use.

PCR analysis

For PCR and real-time PCR analyses of conceptus RNA, isolated RNA (total 0.5 µg) was reverse-transcribed to cDNA using the ReverTra Ace qPCR RT kit (TOYOBO, Osaka, Japan) in a 10 µl reaction volume, and the resulting cDNA (RT template) was stored at 4°C until use. The cDNA reaction mixture was diluted 1 : 10 using DNase-, RNase-free molecular biology grade water. RT template (cDNA) was subjected to PCR or real-time PCR amplification using specific primers (Table 1). PCR-amplified products were separated on 1.5% (w/v) agarose gels after 32 cycles, from which PCR products were subcloned and verified by DNA sequencing. Quantitative PCRs were performed using the SYBR Green kit (Takara Biomedicals, Tokyo, Japan) and the Applied Biosystems thermal cycle system (7900HT, Applied Biosystems, Tokyo, Japan), as previously described [38]. Real-time PCR was performed under the following thermal cycling conditions: 10 min at 95°C, and 40 cycles of 95°C for 10 s followed by 60°C for 30 s. Average cycle threshold (Ct) values for all mRNAs examined were calculated and normalized to Ct values for ACTB mRNA.

Table 1

Primer sets for BERV-K3 and neighboring gene transcripts

BERV-K3 (P1 in Figure 1B)

For detection of the transcript and preparation of probes for in situ hybridization

P1F: 5′-TCGCCCGAAAGCAGGCTAGTGCTAA-3′

P1R: 5′-CAAGGGCGCAGGCTGTTACCTGTTC-3′

For confirming if BERV-K3 is expressed (in Supplementary Figure S3A)

P1F: 5′-TCGCCCGAAAGCAGGCTAGTGCTAA-3′

P1R′: 5′-GCAAGGTTCCGTTTTATGG-3′

For 5′-RACE primer (in Supplementary Figure S3B)

P1R: 5′-CAAGGGCGCAGGCTGTTACCTGTTC-3′

For isolated predicted full length (in Supplementary Figure S3C)

PFF: 5′-ACGGGTAACAAGGAGTCAAAAG-3′

PFR: 5′-CCCTGATGACAAAGTGACCTCC-3′

The primer set to detect LOC100848658 (P2 in Figure 1B) transcript

F: 5′- GCGTCTACCCCAAAACCAGA-3′

R: 5′- ACAGAGAAAGGTGGTCAGGG-3′

The primer set to detect TCF7 (NM_001099186.2)

F: 5′- CTGTGAGCTGGTTCACCCAT-3′

R: 5′- TCCGCAATGACTTTGGCTCT-3′

The primer set to detect ACTB (NM_173979.3)

F: 5′-CTCTTCCAGCCTTCCTTCCT-3′

R: 5′-GGGCAGTGATCTCTTTCTGC-3′

F, forward primer; R, reverse primer.

RNA isolation from bovine conceptus tissues and 5′-RACE for the characterization of the 5′-side of a full-length BERV-K3 transcript

Total RNA was extracted from day 22 bovine conceptuses using the RNeasy Mini Kit together with the RNase-free DNase Set (Qiagen). To identify a full-length BERV-K3 transcript, 5′-RACE with the primer (P1R prime, Table 1 and Supplementary Figure S2) was used to synthesize a first-strand cDNA using the SMARTer RACE 5′/3′ kit (Takara Bio, Inc., Shiga, Japan) according to the manufacturer's instructions. PCR was performed using PrimeSTAR Max (Takara) under the following thermal cycling conditions: 5 cycles of 98°C for 10 s followed by 72°C for 90 s, 5 cycles of 98°C for 10 s followed by 68°C for 90 s, 30 cycles of 98°C for 10 s, 55°C for 10 s followed by 72°C for 90 s, and 1 cycle of 72°C for 10 min. The PCR product was separated by agarose gel electrophoresis and then visualized with a UV trans-illuminator, from which PCR products were subcloned and verified by DNA sequencing.

In situ hybridization

Frozen day 22 uteri embedded in OCT compound were sectioned (10 µm), fixed in 10% neutral buffered formalin, and then subjected to hybridization. Various bovine tissues were fixed in 10% neutral buffered formalin, embedded in paraffin, sectioned (5 µm), and mounted on silane-coated slides (Zyagen, San Diego, CA, U.S.A.). Slide sections were blocked with Block Ace at room temperature for 1 h, incubated with DIG-labeled antisense or sense riboprobe, and then prepared using the BERV-K3 fragment (nucleotide 150 bases) in the pGEM-T Easy Vector with T7 and SP6 promoters (Promega, Madison, WI, U.S.A.) and an RNA transcription kit (Toyobo, Tokyo, Japan) according to the manufacturer's protocol. As described previously [46], hybridization was performed in a humidified chamber at 42°C for 18 h, and bound probes were visualized using alkaline phosphatase-conjugated anti-DIG antibody with nitro blue tetrazolium and 5-brome-4-chloro-3-indolyl phosphate (Promega) as substrates.

Statistical analysis

Quantitative data were subjected to one-way analysis of variance using the Statistical Analyses System (SAS Institute, Cary, NC, U.S.A.); the SEMs shown in the figures were derived from this. When a significant effect of day of pregnancy was detected (P < 0.05), the data were analyzed by Dunnett's test. In these analyses, day of pregnancy was considered an independent source of variation, and replicate was a dependent source.

Results

Searching for ERV-derived elements expressed during the peri-attachment period

We experimentally validated that the expression of ERV-derived sequences that met all of the following criteria: (1) was identified as an ERV-derived sequence from the bovine genome (Btau4.0) by RetroTector, possessing an ORF at least 100 amino acids or has at least one transmembrane domain, (2) was overlapped or located near or in between transcript regions annotated in the Ensembl database (http://www.ensembl.org), allowing two mismatches and annotated more than 10 multiloci, and (3) was detected on at least 2 days (day 17, 20, or 22) by the SOLiD3 sequencing (Supplementary Figure S1). Using these criteria, 10 putative ERVs that are located in between functional genes of the bovine genome were identified (Table 2 and Supplementary Table S1), after which another criterion, up-regulation of transcripts on day 22, was imposed in the analysis. Although the expression of 7 out of 10 putative ERVs declined, the expression of three candidate transcripts was up-regulated in day 20 or 22 pregnant conceptuses. Two transcripts, that both exist on chromosome 7 (ENSBTAT00000052494 and ENSBTAT00000052446, respectively), increased on day 22 (Table 2). In addition, transcript of a third candidate that exists on chromosome 2 (ENSBTAT00000012578) was detected on day 20 when trophoblast attachment to uterine endometrial epithelial cells was initiated, but not on day 22, and was excluded from the subsequent studies.

Table 2

The listing of candidate 10 retroelements detected

Transcript ID	Origin	Tag numbers			RPKM			Length	Chr.	Strand
		D17	D20	D22	D17	D20	D22
ENSBTAT00000012578	Znf	89	105	54	10.814	8.406	5.539	707	2	−
ENSBTAT00000052900	Gag	9	3	0	0.312	0.069	0	2472	4	−
ENSBTAT00000057217	Gag	407	7	0	18.440	0.209	0	1896	6	−
ENSBTAT00000053802	Gag	435	11	0	19.462	0.324	0	1920	6	−
ENSBTAT00000007724	Env	8	6	0	0.376	0.186	0	1827	6	−
ENSBTAT00000052494	Gag	6	1	44	0.318	0.035	1.973	1617	7	+
ENSBTAT00000052446	Gag	7	1	56	0.244	0.030	1.649	2463	7	+
ENSBTAT00000045982	Env	2	1	0	0.102	0.034	0	1680	11	−
ENSBTAT00000056638	Gag	2	4	0	0.070	0.092	0	2469	24	+
ENSBTAT00000029917	Env	157	14	0	6.801	0.400	0	1983	X	−

RPKM, reads per kilo base-pair per million mapped reads; Chr., chromosome.

Identification of candidate ERV sequence and neighboring genes and search for its LTR-like sequences

Using in silico analysis, we re-annotated the equivalent region on chromosome 7 with the recent database, Bos Taurus UMD3.1 version; however, the candidate ERVs previously identified were not detected. PCR analysis subsequently carried out using the primers (P1, Figure 1B) identified a portion of the ERV transcript, ENSBTAT00000052446 (Btau4.0) predicted from the original RNA-seq analysis, from which the sequences were extended through the use of RACE method, resulting in the identification of one ERV sequence (Supplementary Figure S2) existed on chromosome 7 as shown in Figure 1A. Subsequently, we identified the ERV mRNA on chromosome 7, containing Gag_p10, Gag_p24, and zinc finger (ZnF_CCHC) motifs, and short introns in between ORF 1 and 2 (2 nt) and between ORF 2 and 3 (36 nt) as shown in Figure 1B and Supplementary Figure S2. The analysis was extended to determine whether actual transcripts existed in the bovine conceptuses, and found that the single ERV transcript was, in fact, amplified by RT-PCR (Supplementary Figure S3A). The amino acid sequence of the ERV transcript was then used as a query against an ERV ORF database called gEVE [47], identifying that the ERV sequence was quite similar to the gag sequence of BERV-K2 [20] (Supplementary Figure S4A). Because the presence of LTRs (long terminal repeats) or LTR-like sequences was expected, attempts were made to identify the presence of LTR or LTR-like sequences on the chromosomal regions between 20 kb upstream of and 20 kb downstream from the ERV sequence. Two 988 and 1036 nt alignable sequences, located −333 nt upstream of and 6097 nt downstream from the ERV, respectively, were found, of which sequences were also similar to those of BERV-K2 (Supplementary Figure S4B). We, therefore, named the putative ERV transcript as BERV-K3.

Figure 1.

Identification of BERV-K3 and their neighboring genes on chromosome 7 of the bovine genome.

(A) Genes neighboring BERV-K3. Solid-black bar indicates the location of nucleotide sequences and their identity is given below the solid-black bar. Black bars within the gray box is BERV-K3. (B) Structures of BERV-K3, Pfam domains, and primers set to detect transcripts. P1, indicated by a short-black line, is the primer set 1, which detects BERV-K3. P2, indicated by a short-black line, is the primer to detect LOC100848658. Pfam domains (Gag_p10, Gag_p24, zf-CCHC5, dUTPase and RVP) relative to BERV-K3 structures are given below.

BERV-K3 was located between interleukin enhancer-binding factor 3 (ILF3) and QTRT1 genes, and LOC100848658 was also found as a gene that shares a partial sequence (Figure 1A,B). The BERV-K3 sequence was then expected to be located on equivalent regions of chromosomes in human, murine, and sheep genomes; however, BERV-K3 sequence was not found in these species. These data indicated that the BERV-K3 is unique to the bovine species. In some cases, if these LTRs harbored appropriate binding sites for trophoblast-specific transcription factors, they may function as active promoters or enhancers with the potential to modulate neighboring gene expression [48]. We, therefore, examined the expression of ILF3 and QTRT1 with SOLiD3 data (Table 3); however, these gene expressions were not parallel to that of BERV-K3 expression.

Table 3

The Tag numbers and RPKM values of BERV-K3 neighboring genes.

Transcript ID	Gene name	Tag numbers			RPKM
		D17	D20	D22	D17	D20	D22
ENSBTAT00000003403	SLC44A2 (CTL2)	6650	11 253	2454	175.50	195.68	54.67
ENSBTAT00000040076	ILF3	1360	1220	1036	36.04	21.30	23.17
ENSBTAT00000056840	QTRT1	71	28	47	11.25	2.92	6.28
ENSBTAT00000036121	DNM2	984	853	627	24.52	14.00	13.18
ENSBTAT00000014563	TMED1	268	310	416	15.33	11.68	20.08

Expression of ERV-derived transcripts in bovine conceptuses and other tissues

Using real-time PCR analysis, the expression of candidate transcripts was found in days 17, 20, and 22 bovine conceptuses and days 45 and 150 cotyledons. In particular, the expression of BERV-K3 was detected on day 20, when conceptuses begin to attach to uterine epithelial cells, in agreement with RNA-seq and real-time PCR analyses (Figure 2A). LOC100848658 gene was similarly expressed with BERV-K3 in bovine conceptuses, but not in other tissues examined (Figure 2 and Supplementary Figure S5). BERV-K3 was highly expressed in cotyledon on days 45 and 150 pregnancy, and either expression levels were similar to that in day 22 conceptuses. The expression of BERV-K3 in various bovine tissues was then examined (Figure 2B) and found that BERV-K3 was expressed in muscle, uterus, and skin. The expression of BERV-K3 in various bovine cell lines was further examined (Figure 2C) and found that although some expressions were detected in bovine trophoblast cell lines, BERV-K3 was highly expressed in bovine kidney CKT-1 cells.

Figure 2.

Expression of BERV-K3 transcript in bovine conceptuses and various tissues and cell lines.

(A) Changes in BERV-K3 expression. Using primer set P1 described in Figure 1B, expression levels of transcript were determined in day 17, 20, and 22 (D17, D20, and D22) conceptuses (n = 4, each day) and day 45 and 150 cotyledons (n = 3, each day). RNAs from four or three of conceptuses or cotyledons, respectively, were used to perform the qPCR analysis. Data are expressed as fold difference relative to its expression on D17. *, **Statistically significant differences in mRNA levels (P < 0.05 or P < 0.01, respectively) when compared with that in D17. (B) BERV-K3 expression in various tissues. BERV-K3 transcripts in various tissues of day 150 pregnant animals, except that uterine sample was harvested from nonpregnant Japanese black cattle. RNAs were extracted from the spleen, liver, kidney, intestine, skin, lymph node, lung, heart, muscle, and uterus (cyclic) (n = 3 each). Muscle and skin were taken from the upper left hind leg area. Data are expressed as fold difference relative to the expression in the spleen. (C) BERV-K3 expression in bovine cell lines. BERV-K3 transcripts in cultured bovine primary cells and cell lines. CT-1, BT-1, and F3 are bovine trophoblast cell lines. Primary uterine cells (EPI and STR; indicated passage times below), primary cumulus-granulosa cell (oCG; indicated passage times below), primary ear-fibroblast cell (EF; indicated passage times below), intestinal epithelial cell line (bIE), kidney cell lines (MDBK and CKT-1), lymph node cell line (23CLN), and macrophage cell line (BoMAC) were subjected to the present study. Three samples each were used and three independent experiments were conducted. Data are expressed as fold difference relative to the expression in CT-1.

In situ localization of BERV-K3 transcript in the uterus and other tissues

BERV-K3 transcript was found in day 22 uteri, including conceptus, uterine epithelial, and glandular epithelial cells (Figure 3). In addition, BERV-K3 transcript was found in the skin, liver, kidney, and intestine (Figure 4). In the skin, BERV-K3 was localized in outer root and inner root sheath (ORS and IRS, respectively) areas of hair follicles and the stratum spinosum of epidermis. In the liver, BERV-K3 was detected in hepatic laminas of lobule areas, particularly in hepatocytes surrounding central and interlobular veins (ILVs), and in interlobular bile ducts (ILB) of portal canal areas. In the kidney, BERV-K3 transcript was found in epithelial cells of renal-collecting ducts and distal convoluted tubules, and in infiltrating lymphocytes around renal corpuscle; however, BERV-K3 transcript was not found in the cortex region. In the ileum, BERV-K3 was localized in intestinal epithelia and enterocytes in glandules, and in Peyer's patches region. Positive staining was also detected on lymphocytes located in ileum lymphoid follicles and in infiltrating lymphocytes near renal corpuscles. These results indicated that the expression of BERV-K3 was not restricted to the reproductive tissues.

Figure 3.

In situ localization of BERV-K3 transcript in day 22 pregnant uteri.

(A) In situ localization of BERV-K3 transcript in day 22 pregnant animals. (B) In situ localization with the sense-strand probe (negative control). (C) Higher magnification showing luminal epithelia and elongated trophoblasts. (D) Higher magnification showing glandular epithelial regions. GE, glandular epithelia; LE, luminal epithelia; Tr, trophoblast.

Figure 4.

In situ localization of BERV-K3 transcript in the skin, liver, kidney, and ileum.

Bovine tissue sections were purchased from Zymogen (San Diego, CA, U.S.A.). In the skin, BERV-K3 transcript is found in external root and internal root sheath areas of hair follicles and stratum spinosum of epidermis. ORS, outer root sheath; IRS, inner root sheath. In the liver, BERV-K3 is seen in hepatic laminas of lobule areas, particularly in hepatocytes surrounding central veins and ILV, and in interlobular bile ducts of portal canal areas. ILA, interlobular artery; ILB, interlobular bile duct; ILV, interlobular vein. In the kidney, BERV-K3 transcript is found in epithelial cells of renal-collecting ducts and distal convoluted tubules, and in infiltrating lymphocytes around renal-collecting ducts. CD, collecting duct; DCT, distal convoluted tubule; Rc, renal corpuscle. In the ileum, BERV-K3 is seen in intestinal epithelia and enterocytes in glandules, and in PP regions. PP, Peyer's patch. In each, the far right images represent those with the sense-strand cRNA probe.

Induction of BERV-K3 transcript by cell-to-cell contact or canonical WNT agonist

When CT-1 cells were cultured with cell culture insert, not allowing direct CT-1 contact to EECs, no increase in BERV-K3 transcripts was found. However, when CT-1 cells were cultured without the cell insert, allowing direct cell-to-cell contact, increase in BERV-K3 was found (Figure 5A). Because WNT signal is known to be important for placentation following conceptus attachment to the uterine epithelium in the bovine species [29] and our previous study [37] showed that WNT2B and its receptor FZDs mRNA were detected during the conceptus attachment period, we then treated trophoblast CT-1 and F3 cells with 1 µM of canonical WNT agonist for 24 h. BERV-K3 and TCF7, but not LOC100848658, were induced by the WNT agonist in both cells (Figure 5B and Supplementary Figure S6).

Figure 5.

Induction of BERV-K3 expression in CT-1 by cell-to-cell contact or in CT-1 or F3 cells by treatment with the canonical WNT agonist.

(A) Cell-to-cell contact on BERV-K3 expression. CT-1 cells were co-cultured without or with a cell culture insert on EECs for 24 h. (B) Effect of WNT agonist on BERV-K3 expression. CT-1 (left) or F3 (right) cells were treated with 1 µM of canonical WNT agonist for 24 h. Data are expressed as fold difference relative to the expression in control. **Statistically significant differences in mRNA levels (P < 0.05) when compared with that in control.

Discussion

In the present study, we found ERV, putative gag/pol-derived BERV-K3, transcripts specifically expressed in the bovine placenta from early- to mid-gestation. Trophectoderm and fetal membranes had weak expression of BERV-K3 transcript on day 20 and had greater expression from days 22 to 150. Accordingly, the transcripts were detected in bovine trophoblast CT-1, BT-1 and F3 cells. However, the transcripts were also found in the uterus, skin, liver, kidney, and ileum, indicating that their expression appeared somewhat ubiquitous. Cornelis et al. [22] similarly reported that the Bos-Env2 was expressed in the skin, spleen, and muscle, whereas only limited expression was detected in the bovine placenta. To elucidate molecular mechanisms associated with BERV-K3 transcription, the co-culture system with CT-1 cells and EECs [49] was used to study difference in BERV-K3 expression between days 20 and 22, when bovine conceptus begins its attachment to the uterine epithelium cells on day 20, followed by the tight attachment between two cell types on day 22. When CT-1 cells were cultured without the cell culture insert, increase in BERV-K3 transcripts was observed, suggesting that the cell-to-cell contact of trophoblasts to the uterine epithelial cells is required for BERV-K3 transcription to increase.

In our in silico analysis, LTRs (ERVK-type) were found to be located on both upstream and downstream regions of the BERV-K3 gene. However, while BERV-K3 was expressed in the trophectoderm and uterus, additional expression in the skin, liver, kidney, and ileum indicates that these LTRs are not pregnancy-specific. Transcripts of syncytin-Rum1, a fusogenic protein required for ruminants' placentation, were also detected in maternal (caruncle and intercarunclar endometrium) and fetal membranes (chorioallantois) [50]. Several genes are located in close proximity to BERV-K3 on chromosome 7 of the bovine genome, including ILF3 and CTL2 (also known as SLC44A2), located on the upstream region, and DNM2 and TMED, located on the downstream region of the BERV-K3 locus. Increase in BERV-K3 transcripts in day 22 bovine conceptuses was initially considered trophoblast/placenta-specific; however, the expression of these transcripts other than in the trophectoderm indicates that the insertion of BERV-K3 did not cause the placenta-specific expression of TMED and CTL2. It should be noted that although human and murine genomes lack BERV-K3, TMED, and CTL2 are expressed in placentas of these species. In fact, TMED2/p24b1 protein is required for the morphogenesis of the mouse embryo and placental development [51]. The expression of CTL2 transcripts in bovine conceptuses was high, peaking on day 20 when the conceptus begins its attachment to the uterine endometrial epithelium. Rather than the insertion of BERV-K3 in between TMED and CTL2 genes affected the transcriptional regulation of these genes, it is probable that these genes may play a role in the regulation of BERV-K3 gene expression in bovine trophoblasts. Although viral integration to the host genome could occur in a random manner, the integration would have to be locus-specific, if the integrated gene was to become active in certain cell types and/or physiological conditions.

WNT agonist used in the present study as a WNT signal activator was found to induce endogenous BERV-K3 expression in bovine trophoblast CT-1 and F3 cells. The WNT signaling pathway has been established as an important regulator of embryo/conceptus implantation and placental development in mice, sheep, cow, and humans [27-30]. Our previous study in peri-attachment conceptuses revealed that WNT2B and its receptor 9 FZDs mRNA were detected, while DKK-1, which acts as an antagonist of the canonical WNT signaling pathway, was decreased in day 22 bovine conceptuses [37]. Although an experiment determining the effect and function of DKK-1 is beyond the scope of the present study, it is possible that the decrease in the DKK-1 expression after trophoblast attachment to the uterine epithelial cells could lead to activation of canonical WNT signaling. In response to activation of CTNNB and TCF-dependent transcription by the Wnt agonist, BERV-K3 mRNA was increased in trophoblast CT-1 and F3 cells, suggesting that BERV-K3 could be one of the target genes in the canonical WNT pathway. It was also found that epithelial–mesenchymal transition (EMT) occurs in day 22 bovine conceptuses [52,53]. Because the WNT signal is involved in the EMT process, together with WNT, BERV-K3 may play a role in the bovine conceptus EMT and/or its attachment to the uterine epithelium.

Although envelope proteins of HERVs expressed in human trophoblast cells have fusogenic activity and immunosuppressive properties [9,54], the roles that gag polyproteins play have not been definitively elucidated. Most prominent gag- or gag/pol-derived gene is PEG10, which is thought to have been derived from the Ty3/Gypsy family of retrotransposons and is expressed as two proteins: gag-like protein containing a CCHC-type zinc finger domain and gag/pol-like fusion protein with an additional aspartic protease motif [55,56]. The Peg10 knockout mice showed early embryonic lethality due to defects in placental formation [57]. Although the amino acid sequences differ substantially between PEG10 and BERV-K3 (data not shown), BERV-K3 is also a gag/pol-derived ERV and has similar motifs to those of PEG10. BERV-K3 was found to consist of Gag_p24 and Gag_p10 motifs: Gag_p24 forms inner protein layer of the nucleocapsid, while Gag_p10, a retroviral GAG (core) protein, encompasses the p10 region producing the p10 protein upon proteolytic cleavage of GAG by retroviral protease. The p10 is, in fact, associated with viral envelope glycoproteins in most mammalian retroviruses and may be involved in viral particle assembly, transport, and budding [58]. BERV-K3 also consists of the motifs that potentially encode structural proteins, including zinc finger domain, UTPase, and retroviral aspartyl protease; however, the function of these motifs in BERV-K3 has not been elucidated. Nevertheless, because of its insertion in this specific locus, BERV-K3, unique to the bovine species, could have acquired its capability for expression and/or possibly functioning together with neighboring genes in the bovine placenta.

In the present study, BERV-K3 was also detected in IRS, ORS, and stratum spinosum in skins. It was found that the Wnt/Ctnnb signaling system was detected to function in mouse IRS and ORS [59]. In addition, retroviral-like aspartic protease (SASPase) is expressed in the stratum granulosum of mouse skin and is involved in wrinkle formation and adjustments of the stratum corneum hydration [60,61]. BERV-K3 transcripts found in skin epidermis suggest that BERV-K3 protein could be involved in skin development, although further experimentation is required to characterize its molecular mechanism associated with these events. In the liver, BERV-K3 was detected in hepatocytes surrounding central vein and ILV, in which the Wnt/Ctnnb was detected in pericentral hepatocytes in mice [62]. Involvement of canonical Wnt signaling in intestinal epithelial cell proliferation has also been demonstrated [63]. These and results from the present study suggest that in addition to placental formation and its maintenance, BERV-K3 may play a role in the development and/or maintenance of skin, liver, and intestines.

In conclusion, we found the new gag/pol-derived ERV transcript in RNAs extracted from days 20 to 22 bovine conceptuses, of which ERV was integrated into the chromosome 7. The expression of BERV-K3 gene increased in bovine conceptuses following its attachment to uterine epithelia and its high expression was maintained in the placenta up to 150-day gestation. Although the function of the genes is yet to be elucidated, available data suggest that BERV-K3, unique to the bovine species, could be controlled initially through cell-to-cell contact and then by WNT signaling and plays developmental and/or maintenance roles in placenta, skin, liver, and intestines.