<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Cell. Infect. Microbiol.</journal-id>
<journal-title>Frontiers in Cellular and Infection Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Cell. Infect. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">2235-2988</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fcimb.2024.1404431</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Cellular and Infection Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Exploring the role of endogenous retroviruses in seasonal reproductive cycles: a case study of the ERV-V envelope gene in mink</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Yufei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1706041"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Gaofeng</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhu</surname>
<given-names>Yanzhu</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1999199"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Cao</surname>
<given-names>Xiaodong</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Fang</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Huiping</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Liu</surname>
<given-names>Shuying</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
<xref ref-type="aff" rid="aff6">
<sup>6</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Veterinary Medicine, Inner Mongolia Agricultural University</institution>, <addr-line>Hohhot</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Ulanqab Center for Animal Disease Control and Prevention</institution>, <addr-line>Ulanqab</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Institute of Special Animal and Plant Sciences, Chinese Academy of Agricultural Sciences</institution>, <addr-line>Changchun</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>School of pharmacy, Inner Mongolia Medical University</institution>, <addr-line>Hohhot</addr-line>, <country>China</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>College of Veterinary Medicine, Key Laboratory of Basic Veterinary Medicine, Inner Mongolia Agricultural University</institution>, <addr-line>Hohhot</addr-line>, <country>China</country>
</aff>
<aff id="aff6">
<sup>6</sup>
<institution>Key Laboratory of Clinical Diagnosis and Treatment Technology in Animal Disease, Ministry of Agriculture</institution>, <addr-line>Hohhot</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Ye Liu, Van Andel Institute, United States</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Antoinette van der Kuyl, University of Amsterdam, Netherlands</p>
<p>Jian Chen, South China Agricultural University, China</p>
<p>Viol Dhea Kharisma, Airlangga University, Indonesia</p>
<p>Wibke Bayer, Essen University Hospital, Germany</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Shuying Liu, <email xlink:href="mailto:liushuying1968@imau.edu.cn">liushuying1968@imau.edu.cn</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>28</day>
<month>06</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>14</volume>
<elocation-id>1404431</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>03</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>17</day>
<month>06</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Zhang, Wang, Zhu, Cao, Liu, Li and Liu</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Zhang, Wang, Zhu, Cao, Liu, Li and Liu</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>Endogenous retroviruses (ERVs), which originated from exogenous retroviral infections of germline cells millions of years ago and were inherited by subsequent generations as per Mendelian inheritance patterns, predominantly comprise non-protein-coding sequences due to the accumulation of mutations, insertions, deletions, and truncations. Nevertheless, recent studies have revealed that ERVs play a crucial role in diverse biological processes by encoding various proteins.</p>
</sec> <sec>
<title>Methods</title>
<p>In this study, we successfully identified an ERV envelope <italic>(env</italic>) gene in a mink species. A phylogenetic tree of mink ERV-V <italic>env</italic> and reference sequences was constructed using Bayesian methods and maximum-likelihood inference.</p>
</sec> <sec>
<title>Results</title>
<p>Phylogenetic analyses indicated a significant degree of sequence conservation and positive selection within the <italic>env</italic>-surface open reading frame. Additionally, qRT-PCR revealed diverse patterns of mink ERV-V <italic>env</italic> expression in various tissues. The expression of mink ERV-V <italic>env</italic> gene in testicular tissue strongly correlated with the seasonal reproductive cycles of minks.</p>
</sec> <sec>
<title>Discussion</title>
<p>Our study suggests that the ERV-V <italic>env</italic> gene in mink may have been repurposed for host functions.</p>
</sec>
</abstract>
<kwd-group>
<kwd>endogenous retroviruses</kwd>
<kwd>retrovirus</kwd>
<kwd>envelope</kwd>
<kwd>mink</kwd>
<kwd>seasonal breeder</kwd>
<kwd>testis</kwd>
<kwd>testis regression</kwd>
</kwd-group>
<counts>
<fig-count count="12"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="60"/>
<page-count count="18"/>
<word-count count="6362"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Molecular Viral Pathogenesis</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>The retroviruses (RVs) enter host cells by delivering their capsid and viral RNA genomes into the cytoplasm, where viral reverse transcriptase (RT) catalyzes the synthesis of a double-stranded DNA replica of the viral genome. This DNA is then integrated into the host chromosomal DNA. This incorporated DNA becomes an integral part of the genetic material of the infected cell. When integration takes place in germline cells, the virus may pass down through successive generations. Such transmission events can occur in the form of an endogenous viral (or endogenous retroviral; ERV) element (<xref ref-type="bibr" rid="B50">Skirmuntt and Katzourakis, 2019</xref>; <xref ref-type="bibr" rid="B3">Barnes and Price, 2023</xref>; <xref ref-type="bibr" rid="B40">Liu et&#xa0;al., 2023</xref>). ERVs are grouped into the following three categories based on the similarities between their reverse transcriptase (RT) sequences and those of other Retroviridae viruses; class I (gamma- and epsilon-RVs) and class II (alpha-, beta-, delta-, and lenti-RVs) are known as ortho-RVs; class III shows a high similarity with spuma-RVs (<xref ref-type="bibr" rid="B55">Yedavalli et&#xa0;al., 2021</xref>). Despite originating from viruses that undergo rapid evolution, ERVs undergo genetic drift following integration into the host genome and are influenced by the neutral mutation rate of the host species (<xref ref-type="bibr" rid="B49">Simpson et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B37">Kassiotis, 2023</xref>; <xref ref-type="bibr" rid="B15">Dopkins and Nixon, 2024</xref>). Throughout the process of evolution, ERVs were not preserved due to the lack of selective forces, leading to their removal or inactivation. As a result, the accumulation of mutations hindered the production of infectious virions, rendering most ERVs inactive (<xref ref-type="bibr" rid="B49">Simpson et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B37">Kassiotis, 2023</xref>; <xref ref-type="bibr" rid="B15">Dopkins and Nixon, 2024</xref>). Occasionally, certain ERV domains are repurposed by their hosts for cellular functions; such sequences are protected from mutational decay by purifying selection, which prevents the deterioration of genetic sequences that occurs under neutral selection (<xref ref-type="bibr" rid="B37">Kassiotis, 2023</xref>; <xref ref-type="bibr" rid="B15">Dopkins and Nixon, 2024</xref>; <xref ref-type="bibr" rid="B52">Wang et&#xa0;al., 2024</xref>).</p>
<p>ERVs exhibit a dual nature in their relationship with present-day infectious RVs. While some ERVs share similarities with contemporary infectious RVs, others are paleovirological remnants of either extinct or yet-to-be-discovered infectious viruses (<xref ref-type="bibr" rid="B13">Diehl et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B45">Naville and Volff, 2016</xref>; <xref ref-type="bibr" rid="B21">Grandi et&#xa0;al., 2020</xref>). Proliferation of available genome sequences has facilitated the documentation of the widespread distribution of ERVs across vertebrates and enabled genome analyses to compile a comprehensive inventory of viral subtypes in various species and trace instances of cross-species transmission (<xref ref-type="bibr" rid="B24">Hayward et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B60">Zhuo and Feschotte, 2015</xref>; <xref ref-type="bibr" rid="B22">Greenwood et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B55">Yedavalli et&#xa0;al., 2021</xref>). The discovery of ancient paleo-RVs has prompted researchers to reconstruct ancestral genomes (<xref ref-type="bibr" rid="B7">Blanco-Melo et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B55">Yedavalli et&#xa0;al., 2021</xref>), which may not have existing infectious counterparts. Although most mammalian ERVs undergo mutational degradation or truncation, there are occasional instances where an ERV gene is advantageous to the host, leading to its incorporation into a physiological function that has been preserved for millions of years because of evolutionary pressures (<xref ref-type="bibr" rid="B35">Johnson, 2019</xref>; <xref ref-type="bibr" rid="B49">Simpson et&#xa0;al., 2022</xref>). Examples of domesticated genes include viral envelope (<italic>env</italic>) genes that have been repurposed for placental formation and are commonly referred to as syncytins (<xref ref-type="bibr" rid="B14">Dittmar et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B9">Bu et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B34">Jern and Greenwood, 2024</xref>). Additionally, the <italic>env</italic> and <italic>gag</italic> genes have been found to display antiviral capabilities, as exemplified by Fv4 and Fv1 (<xref ref-type="bibr" rid="B31">Ikeda et&#xa0;al., 1985</xref>; <xref ref-type="bibr" rid="B5">Best et&#xa0;al., 1996</xref>; <xref ref-type="bibr" rid="B30">Ikeda et&#xa0;al., 2001</xref>). Furthermore, regulatory sequences that influence host gene expression have been identified (<xref ref-type="bibr" rid="B10">Buzdin et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B43">Meyer et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B59">Zhou et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B48">Plant et&#xa0;al., 2023</xref>).</p>
<p>Minks (<italic>Neovison vison</italic>) are one of the most important fur animals; besides providing fur, they are widely used as experimental models in biomedical research. Minks exhibit a yearly pattern of testicular activity, involving shifts between suppressed and highly active spermatogenesis. Thus, they are valuable subjects for investigating the fundamental mechanisms governing spermatogenesis (<xref ref-type="bibr" rid="B8">Blottner et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B47">Pelletier et&#xa0;al., 2015</xref>). The circannual cycle of minks encompasses three primary stages of testicular function&#x2014;(i) reactivation of testicular activity in autumn (November), (ii) peak sexual activity (February to March), and (iii) sexual inactivity (May to October) (<xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>).</p>
<p>In this study, we observed upregulation of ERV-derived genes (here named mink ERV-V <italic>env</italic>) during testicular regression. However, the origin and biological function of mink ERV-V <italic>env</italic> remain unclear. A phylogenetic tree was constructed, using Bayesian methods and maximum-likelihood inference, for the mink ERV-V <italic>env</italic> sequences, and reference sequences for both endogenous and exogenous viral envelope genes were included in it. A significant degree of sequence conservation and positive selection was seen within the <italic>env</italic>-surface ORFs in their respective orders; qRT-PCR revealed a wide-ranging expression pattern of mink ERV-V <italic>env</italic> in various tissues. Notably, the expression of mink ERV-V <italic>env</italic> gene in the testicular tissue strongly correlated with their seasonal reproductive cycle. This study would make a significant contribution to the literature because it suggested that the ERV-V <italic>env</italic> gene in mink may have been repurposed to perform different host functions.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Ethics statement</title>
<p>The Institutional Animal Care and Use Committee of the Chinese Academy of Agricultural Sciences in Jilin Province, China, reviewed and approved the protocol for animal experiments (Approval number: 2020007) (<xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>). The experiments were conducted in collaboration by the Institute of Special Animal and Plant Sciences, the Chinese Academy of Agricultural Sciences, and Inner Mongolia Agricultural University. Animal testing was done at the Institute of Special Animal and Plant Sciences and sample analysis at Inner Mongolia Agricultural University. Animal testing was carried out following the guidelines of the Animal Research Reporting of <italic>In Vivo</italic> Experiments (ARRIVE).</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Sample collection</title>
<p>Male black minks (<italic>Neovison vison</italic>) that were eleven months old were obtained from the research facility at the Chinese Academy of Agricultural Sciences. Each mink was kept in its own cage, given a daily portion of commercial food, exposed to sunlight, and allowed to freely eat and drink (<xref ref-type="bibr" rid="B57">Zhang and Xu, 2022</xref>; <xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>).</p>
<p>The minks were euthanized by CO<sub>2</sub> inhalation. Black mink testes were collected in January (n = 6), February (n = 6), March (n = 6), April (n = 6), May (n = 6), June (n = 6), July (n = 6), August (n = 6), September (n = 6), October (n = 6), November (n = 6), and December (n = 6). The minks used in this experiment were identical to those in <xref ref-type="bibr" rid="B57">Zhang and Xu, 2022</xref> and <xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Isolation and subcloning of Mink ERV-V env cDNA</title>
<p>To construct plasmids expressing C-terminal His-tagged mink ERV-<italic>env</italic>, we used total RNA from mink testes and subjected it to RT-PCR using a Superscript III one-step RT-PCR kit (Invitrogen, Grand Island, NY). The RT-PCR program was as follows: incubation at 55&#xb0;C for 30 min, denaturation at 95&#xb0;C for 3 min, 35 cycles of denaturation at 95&#xb0;C for 30 s, annealing at 60&#xb0;C for 30 s, and extension at 72&#xb0;C for 2 min, followed by a final extension step at 72&#xb0;C for 10 min. The relevant primers (F: 5&#x2019;- CGCGGATCCCGGCCATGATGAGACTGATGAAGTGGAC -3,&#x2019; R: 5&#x2019;- CGCGAATTCCGTGAACCTCCACCTCCTGAACCTCCACCTCCGGTAGAGAGCCTCTGCCTGCGTC-3&#x2019;) were used for RT-PCR. The resulting PCR products were purified using a QIAquick gel extraction kit (QIAGEN, Shanghai, China), and the mink-ERV-V<italic>-env</italic> ORFs were cloned into the pCDNA4/Myc-His plasmid using EcoRI (Thermo Fisher Scientific, Shanghai, China) and EcoRV (Thermo Fisher Scientific, Shanghai, China) restriction enzymes.</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Quantitative reverse transcription-PCR amplification</title>
<p>In accordance with the manufacturer&#x2019;s instructions, an RNeasy Mini Extraction Kit and RNase-free DNase Set (Qiagen, Shanghai, China) were used to extract total RNA. The first strand of cDNA was synthesized from 2.5 &#x3bc;g of total RNA using SuperScript III (Invitrogen, Grand Island, NY) primed with oligo dT. All qRT-PCR reactions were performed using 2&#xd7;SYBR Green qPCR Master Mix (Life Technologies, Grand Island, NY). All qRT-PCR analyses were performed using Roche LightCyler 480 equipment with a 20-&#x3bc;L reaction volume. The RT-PCR primers used were <italic>env</italic>-F (5&#x2019;- TAGGAGAGCTCTTGGGCGTT-3&#x2019;) and <italic>env</italic>-R (5&#x2019;-TCCTGGCGGCTCTCTAAACT-3&#x2019;). Relative expression of the selected genes was adjusted using <italic>GAPDH</italic> as an endogenous control (<xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>). The study employed the comparative cycle threshold method to determine the amplification fold, with the cycle threshold (Ct) representing the number of PCR cycles needed to reach the fluorescence threshold. The calculation of &#x394;Ct involved subtracting the Ct value of the reference gene GAPDH from that of the target gene. The fold change was calculated as 2<sup>&#x2013;&#x394;&#x394;Ct</sup>, where &#x394;&#x394;Ct is defined as the difference between the &#x394;Ct values of the sample and control groups. Each time point, three testes were collected from three different animals, each providing one RNA sample.</p>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Database search for mink ERV-V env orthologs</title>
<p>ORF sequences from mink ERV-<italic>env</italic> were utilized as probes in a BLAST query of 66 Carnivora genome assemblies archived in the GenBank repository. A BLASTn query was performed with gap costs set at 5 and 2 for existence and extension, match/mismatch scores of +2/&#x2212;3, repeat masking filter turned off, and an expect threshold of 10<sup>&#x2212;75</sup> established.</p>
</sec>
<sec id="s2_6">
<label>2.6</label>
<title>Phylogenetic analysis</title>
<p>We examined the genomic area encompassing BTN1A1 and BTN2A1 as well as the regions extending 150 kb past SPTLC3 to study synteny in Carnivora for ERV-<italic>env</italic> and associated ERV duplicates (<xref ref-type="bibr" rid="B49">Simpson et&#xa0;al., 2022</xref>). The segments were obtained from the latest GenBank assembly for each species (<xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Table S1</bold>
</xref>), and a sequence alignment was performed using MultiPipMaker (<xref ref-type="bibr" rid="B11">Chen et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B29">Horie and Tomonaga, 2019</xref>; <xref ref-type="bibr" rid="B6">Blanco-Melo et&#xa0;al., 2024</xref>).</p>
<p>Phylogenetic trees were built using the neighbor-joining method in MEGA-X (<xref ref-type="bibr" rid="B38">Kumar et&#xa0;al., 2018</xref>). Two trees were constructed using the RT domain of pol and a portion of the transmembrane subunit of the <italic>env</italic> gene. <xref ref-type="supplementary-material" rid="ST2">
<bold>Supplementary Table S2</bold>
</xref> shows the ERV sequences utilized for tree construction. The JTT matrix-based method (<xref ref-type="bibr" rid="B38">Kumar et&#xa0;al., 2018</xref>) was used to compute evolutionary distances, and a gamma distribution (shape parameter = 1) was employed to model the rate of variation across sites. Positions with site coverage below 95% were removed to ensure less than 5% alignment gaps, missing data, or ambiguous bases at any position (<xref ref-type="bibr" rid="B39">Langmead and Salzberg, 2012</xref>; <xref ref-type="bibr" rid="B1">Alser et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B53">Warnow, 2021</xref>). The amino acid sequences of mink ERV-V <italic>env</italic> were aligned using the MUSCLE algorithm within Geneious Prime with default parameters. Maximum-likelihood phylogenetic trees were constructed using the RaxML program with the General Time Reversible + G + I model and 500 bootstrap replicates for branch support.</p>
</sec>
<sec id="s2_7">
<label>2.7</label>
<title>Bioinformatic analysis</title>
<p>The JPred secondary structure prediction server (<xref ref-type="bibr" rid="B16">Drozdetskiy et&#xa0;al., 2015</xref>) was utilized to query the amino acid sequences of mink ERV-V Env and human EnvV2. The resulting JNetPRED and JNetCONF graphs were extracted and analyzed to determine their secondary structures. ProtScale was used to plot the hydropathicity scores for these proteins (<ext-link ext-link-type="uri" xlink:href="https://web.expasy.org/protscale/">https://web.expasy.org/protscale/</ext-link>).</p>
</sec>
<sec id="s2_8">
<label>2.8</label>
<title>Homology modeling</title>
<p>The mink ERV-V Env sequence was analyzed using the I-TASSER (<xref ref-type="bibr" rid="B54">Yang et&#xa0;al., 2015</xref>) program, which employs a multiple-threading approach to identify homologs. The approach involves identifying templates in the Protein Data Bank (PDB), conducting iterative structure-assembly simulations, selecting and refining models, and annotating structure-based functions.</p>
</sec>
<sec id="s2_9">
<label>2.9</label>
<title>Cell culture, transfection, western blotting, and immunofluorescence</title>
<p>Cells (293T) were grown in a 12-well dish with 200,000 cells per well and then transfected with either 1000 ng pcDNA4-mink-ERV-<italic>env</italic>-His or pcDNA4-empty plasmids using Lipo8000&#x2122; Transfection Reagent (Beyotime Biotechnology, Shanghai, China). Cells were lysed with RIPA buffer (Beyotime Biotechnology, Shanghai, China) after three days using the same number of cells. Equal volumes of lysates were loaded onto a 4&#x2013;12% Bis-Tris polyacrylamide gel for SDS-PAGE analysis. Immunoblot analysis was conducted on a PVDF membrane, and the results were captured with a Tanon 5200 imaging system from Tanon (Shanghai, China). The following antibodies were used for protein detection: Anti-6&#xd7; His tag antibody (Cat. No. ab18184; Abcam Plc., Cambridge, UK); anti-GAPDH antibody (Cat. No. 60004&#x2013;1-Ig, Proteintech Group Inc., Wuhan, China); Goat Anti-Rabbit IgG H&amp;L (HRP) (Cat. No. ab205718; Abcam Plc., Cambridge, UK).</p>
<p>Cultured mink testis cells on 25-mm coverslips were transfected with 200 ng of either pCDNA4/Myc-His or pcDNA4-mink-ERV-<italic>env</italic>-His plasmid in 12-well plates. We fixed the cells with 3.7% formaldehyde and permeabilized them with PBS containing 0.1% Triton X-100 after two days. Primary mouse anti-His antibody (Cat. No. ab18184; Abcam Plc., Cambridge, UK) was added to 10 mM PBS containing 10% normal goat serum. After three washes with PBS, the sections were exposed for 30 minutes to secondary antibodies (Goat Anti-Mouse IgG H&amp;L (Alexa Fluor<sup>&#xae;</sup> 647); Cat. No. ab150115; Abcam Plc., Cambridge, UK) diluted in PBS. In negative controls, primary antibody is PBS. Cell nuclei were stained with DAPI (Beyotime Biotechnology, Shanghai, China) and then placed on glass slides using the ProLong antifade kit (Thermo Fisher Scientific) for observation under a Leica fluorescence microscope DM2500.</p>
</sec>
<sec id="s2_10">
<label>2.10</label>
<title>Statistical analysis</title>
<p>A one-way ANOVA analysis of the data was performed with the PRISM software, and then Tukey&#x2019;s multiple comparison test was conducted using GraphPad Software, Inc. The findings are displayed as the average &#xb1; standard error of the mean of values from a minimum of three biological duplicates. Statistical significance was set at a p-value of less than 0.05.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Differential expression of mink ERV-V env associated with total testicular regression and recrudescence in mink</title>
<p>Male minks of seasonal brooding breeds experienced a decrease in seminiferous tubule diameter and a consequent decrease in testis size during testis regression (<xref ref-type="fig" rid="f1">
<bold>Figures&#xa0;1A&#x2013;D</bold>
</xref>). In minks, the testes underwent significant circannual changes in both structure and function, involving the elimination of the germinative epithelium during the non-breeding period (<xref ref-type="supplementary-material" rid="SF1">
<bold>Supplementary Figure S1</bold>
</xref>). This phenomenon resulted in notable alterations in the processes and structures, such as androgen production, apoptosis, cell proliferation, and integrity of the blood-testis barrier (BTB) (<xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>). In a previous transcriptomic study that examined mink testes at different reproductive stages (active, regressing, and inactive), we identified an interesting gene (specifically, an endogenous retrovirus <italic>env</italic> gene) and named it mink ERV-<italic>env</italic>. In the current study, we found that variation in expression of this gene was highly consistent with testis regression. Quantitative PCR results revealed that its expression was elevated during testis regression in minks, remained at high levels during the non-breeding season, decreased during testis recrudescence, and was maintained at low levels during the breeding season (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1E</bold>
</xref>).</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Yearly seasonal breeding cycle in a typical mature mink. <bold>(A)</bold> Cellular changes in the testes during the annual seasonal reproductive cycle. <bold>(B)</bold> Variation in mink testis weights during the annual seasonal breeding cycle. <bold>(C)</bold> Variation in mink testis major axis during the annual seasonal breeding cycle. <bold>(D)</bold> Variation in mink testis transverse diameter during the annual seasonal breeding cycle. <bold>(E)</bold> ERV-<italic>env</italic> expression in normal adult mink testes during the annual seasonal reproductive cycle. A one-way ANOVA was employed to determine significant differences between the different time points (*p &lt; 0.05; **p &lt; 0.01; ***p &lt; 0.001).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g001.tif"/>
</fig>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Conserved env ORF in carnivores related to the mink ERV-V env ORF</title>
<p>Using tBLASTn, we analyzed annotated mammalian genomes in-silico to discover ERVs associated with mink ERV-<italic>env</italic>. As anticipated, the majority of matches corresponded to <italic>env</italic>V orthologs in Carnivora. A more detailed analysis of the genetic location of this ERV in the mink species <italic>Neovison vison</italic> showed that it is situated between the <italic>BTN1A1</italic> and <italic>BTN2A1</italic> genes (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). Furthermore, a coding region was identified at the 3&#x2032; terminus of this ERV, located before the expected 3&#x2032; long terminal repeat (LTR) (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). Further analysis of mink ERV-V in the Dfam repeat database uncovered remnants of the <italic>gag</italic> and <italic>pol</italic> genes that resembled an ERV cluster known as Prima41, whereas the C-terminal portion showed similarities to the ERV cluster MER66 (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Genomic location of mink ERV-<italic>env</italic>. Location of the mink ERV-<italic>env</italic> gene is illustrated on chromosome 1 (NC_058091.1) in the genome of <italic>Neovison vison</italic> (mink). Green indicates RefSeq annotated genes, while yellow represents the mink ERV-<italic>env</italic> ORF. Blue annotations indicate ERV-derived repetitive elements that were identified in the specified region based on data in the Dfam repeat database.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g002.tif"/>
</fig>
<p>Sections of the mink ERV-V <italic>env</italic> ORF were retrieved and compared with homologous sequences from 18 Carnivora species and two closely related species in the superorder Laurasiatheria (<xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Table S1</bold>
</xref>). The mink ERV-V <italic>env</italic> ORF showed significant sequence similarity across Carnivora, with the majority of exons in the associated <italic>BTN1A1</italic> and <italic>BTN2A1</italic> genes in all species analyzed resembling one another, thereby confirming the orthologous nature of this region (see <xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). Furthermore, the region encompassing the mink ERV-V <italic>env</italic> ORF, along with the rest of the provirus, was found to be absent in all mammalian species other than those in Carnivora (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). The mink ERV-V <italic>env</italic> gene, with its long ORF (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>), is preserved in all Carnivora lineages except Hyaenidae, where premature stop codons in the orthologous gene lead to the formation of shortened ORFs (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Genomic synteny of mink ERV-<italic>env</italic>. Genomic segments containing the mink ERV-<italic>env</italic> ORF and <italic>BTN1A1</italic> and <italic>BTN2A1</italic> genes were extracted from the NCBI genome database for each of the indicated species and aligned. The mink ERV-<italic>env</italic> ORF, along with the BTN1A1 and BTN2A1 exons, are highlighted in red and black arrows in the reference assembly.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g003.tif"/>
</fig>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Conservation of mink ERV-<italic>env</italic> ORF in Carnivora lineages. Mink ERV-env-containing ERV orthologs were identified through BLAST searches of individual genome assemblies. MEGA v7.0.26 was used to create this cladogram, which shows the evolutionary connections among species with complete genomes in the order Carnivora. Classifications of suborder and family are displayed on the right side.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g004.tif"/>
</fig>
<p>In a manner akin to host genes, ERV genes, once integrated and stabilized within a lineage, may undergo conservation through purifying/negative selection to maintain advantageous functionality, reflected by a dN/dS (&#x3c9;) ratio below 1, or alternatively, they may experience evolutionary changes through diversifying/positive selection, leading to adaptive modifications and a &#x3c9; value exceeding 1. To discern the indications of selection pressures within mink ERV-V <italic>env</italic> gene, we employed diverse analytical approaches available on the Datamonkey web server. Three residues (133R, 437D, 448R) within the mink ERV-V <italic>env</italic> were identified as undergoing positive selection by the MEME and FEL programs on the Datamonkey webserver. While the orthologs of mink ERV-V <italic>env</italic> displayed high sequence similarity at most sites, the pairwise &#x3c9; values were predominantly below unity across the majority of species. These findings suggest that mink ERV-V <italic>env</italic> have experienced repeated episodes of positive selection within their respective orders.</p>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Identification of the mink ERV-V proviruses as gamma-like ERVs</title>
<p>Segments of retroviral genomes can exhibit diverse phylogenetic histories because of frequent retroviral recombination, which can involve distantly related RVs (<xref ref-type="bibr" rid="B26">Henzy et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B55">Yedavalli et&#xa0;al., 2021</xref>) or occur between endogenous and exogenous viruses (<xref ref-type="bibr" rid="B2">Bamunusinghe et&#xa0;al., 2017</xref>). RVs are categorized into two subfamilies, namely orthoretroviruses and spumaretroviruses (<xref ref-type="bibr" rid="B49">Simpson et&#xa0;al., 2022</xref>). Orthoretroviruses are further classified into six genera, namely <italic>Alpharetrovirus</italic>, <italic>Betaretrovirus</italic>, <italic>Gammaretrovirus</italic>, <italic>Epsilonretrovirus</italic>, <italic>Lentivirus</italic>, and <italic>Deltaretrovirus</italic> (<xref ref-type="bibr" rid="B49">Simpson et&#xa0;al., 2022</xref>). ERVs are categorized into three classes&#x2014;class I (gamma and epsilon), class II (alpha, beta, delta, and lenti), and class III (spuma)&#x2014; depending on how closely they resemble exogenous RVs (<xref ref-type="bibr" rid="B20">Gifford et&#xa0;al., 2018</xref>). Retroviral envelope proteins can be classified into two subtypes, namely gamma-type and beta-type, depending on the linkage type (covalent or non-covalent) between the SU and TM subunits (<xref ref-type="bibr" rid="B27">Henzy and Johnson, 2013</xref>). In the phylogenetic analysis of both endogenous and exogenous retroviral envelope proteins, mink ERV-V Env proteins grouped together with those from the gamma-type envelope subgroup, which has SU and TM subunits linked together covalently. This subgroup encompassed gamma, delta, and alpha RVs (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>). Mink ERV-V <italic>env</italic>, CARenvV, and ARTenvV formed a group with human EnvV2 (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>). By comparing amino acid sequences of RVs from different sources, including those of mink ERV-V ERVs, we created a phylogenetic tree to classify retroviral groups. All seven RV genera had highly conserved RT cores, and the tree displayed three groups corresponding to ERV classes I&#x2013;III (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>). <xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref> shows that the RT from mink ERV-V ERVs grouped together with external gammaretroviruses and early mammalian gamma-like class I ERVs such as FeLV, PERV-C, and AKV.</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Phylogeny of mink ERV-env protein. The predicted amino acid sequence of mink-ERV-env (highlighted in red) was juxtaposed with the amino acid sequences of selected endogenous and exogenous retroviral Env proteins, leading to the construction of a maximum likelihood tree with 500 bootstraps in RaxML.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g005.tif"/>
</fig>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Phylogeny of mink ERV RT protein. Mink-ERV RTs were aligned with those of endogenous and exogenous RVs from different species, and a maximum likelihood tree was constructed with 500 bootstraps using RaxML.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g006.tif"/>
</fig>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Prediction of conserved domains and functional classification of mink ERV-V ENV protein</title>
<p>To ascertain the presence of retroviral Env protein characteristics in the mink ERV-V Env protein, hydrophobicity plots were generated for the protein with human EnvV2. In addition, predicted secondary structures of these proteins were generated. The results of our analysis suggested the presence of a possible N-terminal signal peptide containing both a hydrophobic area and an alpha helix. Additionally, a potential fusion peptide was discovered in the C-terminus that contained both a hydrophobic section and an alpha helix (see <xref ref-type="fig" rid="f7">
<bold>Figures&#xa0;7</bold>
</xref>, <xref ref-type="fig" rid="f8">
<bold>8</bold>
</xref>). A comparison of the amino acid sequences of human EnvV2 and mink ERV-V Env protein showed that no full transmembrane region was present.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Hydrophobicity profiles are shown for human EnvV2 (green), <italic>Felis catus</italic> CARenvV (blue), and mink ERV-V ENV (red).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g007.tif"/>
</fig>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Secondary structure predictions are shown for human EnvV2, cat CARenvV, dog CARenvV, and mink ERV-V Env. In secondary structure prediction, whorl-like structures indicate alpha helices and arrows indicate beta sheets.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g008.tif"/>
</fig>
<p>The structure and characteristics of the canonical retroviral Env protein are shown in <xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9A</bold>
</xref>. Human EnvV2 showed the most resemblance with these Env proteins in the area right after the expected signal peptide (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9B</bold>
</xref>). The mink ERV-V Env protein was truncated immediately before the membrane-spanning domain, which was located close to the Env C-terminus (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9B</bold>
</xref>). The alignment revealed two conserved cysteine-containing motifs in the mink ERV-V Env protein, CxxC in the SU domain, and Cx6CC in the transmembrane domain, which were present at the same location as in EnvV2 (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9B</bold>
</xref>). The mink ERV-V Env protein contained the RxxR furin cleavage motif (RQT/SR) immediately before the anticipated fusion peptide (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9B</bold>
</xref>). The immunosuppressive domain and heptad repeats found in the mink ERV-V Env protein are similar to those found in the previously discovered EnvV2 immunosuppressive domain and heptad repeats (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9B</bold>
</xref>).</p>
<fig id="f9" position="float">
<label>Figure&#xa0;9</label>
<caption>
<p>Structural features of the mink ERV-V Env protein. <bold>(A)</bold> Protein domain pattern map of mink ERV-V Env. <bold>(B)</bold> Comparison of the amino acid sequences of human EnvV2 with those of ERV-V Env from dog, cat, fox, and mink. The purple vertical line indicates the signal peptide cleavage site. The red square indicates the CXXC motifs that are involved in SU-TM interactions. The furin cleavage site (RXXR) is denoted by the green square. The pink square indicates the fusion peptide. The brown square brackets indicate the putative immunosuppressive domain (ISD). The green square indicates a Cx6CC motif. The purple square indicates the transmembrane domain (TM). Black shading indicates the residues that are identical in all proteins, while different shades of grey background indicate the residues that are identical in four, three, or two proteins.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g009.tif"/>
</fig>
</sec>
<sec id="s3_5">
<label>3.5</label>
<title>Structural analysis of the protein coded by the mink ERV-V Env gene</title>
<p>We next modeled the three-dimensional structure of the mink ERV-V Env protein. Five models were created on the I-TASSER web server, and the best model was selected and visualized using PyMOL (<xref ref-type="fig" rid="f10">
<bold>Figure&#xa0;10A</bold>
</xref>). The possible signal peptide and furin cleavage motif of mink ERV-V Env are shown in <xref ref-type="fig" rid="f10">
<bold>Figure&#xa0;10A</bold>
</xref>. Attempts were made to predict the protein modification sites in mink ERV-V Env; results indicated six potential N-glycosylation sites in the disordered region of mink ERV-V Env at amino acid positions 36 (NISS), 46 (NLSD), 108 (NQTG), 178 (NFSC), 314 (NITS), and 437 (NHTA) (<xref ref-type="fig" rid="f10">
<bold>Figure&#xa0;10B</bold>
</xref>). The STRING database v12.0 was used to study the interactions between mink ERV-V Env and other proteins (<xref ref-type="fig" rid="f10">
<bold>Figure&#xa0;10C</bold>
</xref>). Data from the STRING database indicated that interactions between mink ERV-V Env and ubiquitin carboxyl-terminal hydrolase 35 (USP35), USP38, and Ig-like domain-containing protein. In transfected 293T cells, the mink ERV-V Env protein labeled with histidine showed consistent expression when inserted under the constitutive CMV promoter (<xref ref-type="fig" rid="f10">
<bold>Figure&#xa0;10D</bold>
</xref>). The aforementioned expression plasmids were transfected into the primary testicular cells of minks, which resulted in the accumulation of this protein in the cytoplasm (<xref ref-type="fig" rid="f10">
<bold>Figure&#xa0;10E</bold>
</xref>). Importantly, the Env protein from different RVs was found in the cytoplasm during infection, but as the viral life cycle progressed, this protein localized to the plasma membrane for assembly. A major reason for this is that the mink ERV-V Env protein has lost its transmembrane region, which prevents it from being expressed at the cytoplasmic membrane.</p>
<fig id="f10" position="float">
<label>Figure&#xa0;10</label>
<caption>
<p>Comparative homology modeling and functional analysis of the mink ERV-V Env protein. <bold>(A)</bold> Three-dimensional ab-initio models of mink ERV-V Env proteins. <bold>(B)</bold> Post-translationally modified glycation sites in the mink ERV-V Env protein. <bold>(C)</bold> Map of protein-protein interactions with the mink ERV-V Env protein. <bold>(D)</bold> Western blot examination of His-tagged mink ERV-V Env protein in 293T cells. Lane 1, a pcDNA4-mink-ERV-Env&#x2013;His plasmid was transfected into 293T cells. Lane 2, a pCDNA4/Myc-His plasmid was transfected into 293T cells. <bold>(E)</bold> Confocal microscopic examination of His-tagged mink ERV-V Env protein in mink testis cells shows accumulation in the cytoplasm (scale bar, 100 &#x3bc;m).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g010.tif"/>
</fig>
</sec>
<sec id="s3_6">
<label>3.6</label>
<title>Mink ERV-V env gene expressed in various tissues</title>
<p>To ascertain the tissue-specific expression profile of mink ERV-V Env gene, we employed RT-qPCR. GAPDH was utilized as the internal reference gene. Mink ERV-V Env gene exhibited varying expression levels across a diverse range of tissues (Tissues from minks were collected in February), as shown in <xref ref-type="fig" rid="f11">
<bold>Figure&#xa0;11</bold>
</xref>. While reduced expression levels of mink ERV-V <italic>env</italic> gene were observed in the heart, liver, kidney, and spleen, increased expression was observed in the uterus and pancreas.</p>
<fig id="f11" position="float">
<label>Figure&#xa0;11</label>
<caption>
<p>Mink ERV-Env gene is expressed in a wide variety of tissues. A one-way ANOVA was employed to determine significant differences between the different tissue samples (*p &lt; 0.05; **p &lt; 0.01; ***p &lt; 0.001).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g011.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>Seasonal breeding refers to the absence of continuous reproductive activity throughout the year in specific species or populations residing in non-equatorial regions of the Earth that occur as a consequence of cyclical environmental variations caused by climatic seasons. In females, ovulation is suspended during a specific period, while in males, testicular function is diminished or exhausted during the non-breeding season (<xref ref-type="bibr" rid="B4">Beltran-Frutos et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B23">Gumulka et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B51">Tao et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B46">Pandey and Mohanty, 2023</xref>; <xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>). In seasonal breeding animals such as minks, testes undergo cyclical activation and inactivation annually, which is a particularly interesting example of spermatogenesis regulation (<xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>). In this study, we identified a novel gene, mink ERV-V <italic>env</italic>, which showed expression in testes that was accompanied by cyclic changes in spermatogenesis. The expression of mink ERV-V env gradually increased during the regression stage of mink testes and gradually declined during the transition to the reproductive stage.</p>
<p>In this study, we identified ERVs that are related to mink ERV-V env, the majority of which correspond to Carnivora envV orthologs. The sequences likely represent ancient coding sequences that have been preserved. Endogenous viruses evolve at the host&#x2019;s mutation rate rather than at the rate of their viral progenitors once they enter the host genome through vertical transfer (<xref ref-type="bibr" rid="B50">Skirmuntt and Katzourakis, 2019</xref>). This phenomenon gives rise to two potential outcomes for such sequences, namely the accumulation of mutations, which results in partial loss of the coding sequence and gradual deterioration of its function if it proves detrimental or non-beneficial for the host&#x2019;s survival, or fixation and preservation of these sequences within the genomes of the host population if they confer a selective advantage (<xref ref-type="bibr" rid="B17">Feschotte and Gilbert, 2012</xref>; <xref ref-type="bibr" rid="B50">Skirmuntt and Katzourakis, 2019</xref>). The observation of numerous viral sequences in ERVs associated with Carnivora envV orthologs indicates that these integrations occurred in the distant past and have since been deactivated and undergone decay. A small number of species developed premature stop signals in ORFs homologous to mink ERV-V <italic>env</italic>, which are preserved in every lineage across orders. Alignment and phylogenetic analyses of ERV-V <italic>env</italic> sequences will not be affected by this premature stop signal.</p>
<p>In our evolutionary tree, we noticed that mink ERV-V Env proteins grouped together with those from the gamma-type envelope subgroup, known for having surface and transmembrane subunits linked by covalent bonds. Various vertebrate species, including species of mammals, birds, reptiles, and amphibians, harbor ERV loci related to <italic>Gammaretrovirus</italic> (<xref ref-type="bibr" rid="B20">Gifford et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B28">Hogan and Johnson, 2023</xref>). Overall, Gammaretrovirus-related ERV represent a &#x201c;fossil&#x201d; record of the history of retroviruses. This study revealed that mink ERV-V Env contains a possible N-terminal signal peptide that includes both a hydrophobic area and an alpha helix. Additionally, a fusion peptide was predicted in the C-terminus, which contained a hydrophobic segment and an alpha helix. Although mink ERV-V Env lost the transmembrane region and the cytoplasmic tail. But the mink ERV-V Env TM domain has one hydrophobic stretch that flanks an ectodomain. The latter contains two heptad repeats (hr1 and hr2) flanked by two or three cysteine residues (<xref ref-type="bibr" rid="B27">Henzy and Johnson, 2013</xref>). Heptad repeats are of utmost importance in facilitating dynamic rearrangement of the trimer during the fusion process as well as in the formation of the widely conserved coiled-coil structure observed in numerous viral fusion proteins (<xref ref-type="bibr" rid="B12">Chen et&#xa0;al., 2022</xref>). Additionally, the ectodomain sequence of the mink ERV-V Env transmembrane subunit encompasses a segment referred to as the immunosuppressive domain, which consists of 20 amino acids directly preceding the initial cysteine residue and is distinguished by its conserved residues.</p>
<p>Envelope glycoproteins of the gamma type undergo post-translational modifications such as glycosylation and furin-mediated cleavage, similar to the envelope glycoproteins of other RVs. Gamma-type Envs are estimated to have 4&#x2013;13 glycosylation sites, the majority of which are located in the surface subunit of the protein. Our study identified six potential N-glycosylation sites in the disordered region of the mink ERV-V Env surface subunit. While there was significant variation in the glycosylation sites, there was far more conservation in the furin cleavage sites. The mink ERV-V Env protein also contains an RxxR furin cleavage motif (RQT/SR) immediately before the anticipated fusion peptide.</p>
<p>Our study results indicate that the mink ERV-V Env mRNA is present in different tissues, and the ORF of this gene can reliably produce proteins when activated by a constitutively active promoter. However, the underlying reasons for the conservation of these retroviral genes remain unclear. It is important to note that ERV loci in vertebrate genetic material are subject to natural selection, with increasing evidence showing the incorporation of ERV-produced proteins during the evolution of host cells and organisms for various functions (<xref ref-type="bibr" rid="B41">Malfavon-Borja and Feschotte, 2015</xref>; <xref ref-type="bibr" rid="B18">Frank and Feschotte, 2017</xref>; <xref ref-type="bibr" rid="B35">Johnson, 2019</xref>). For instance, <italic>Gammaretrovirus</italic>-related ERVs encompass placental syncytin genes that have been evolutionarily co-opted by various mammalian lineages (<xref ref-type="bibr" rid="B56">Zhang et&#xa0;al., 2015</xref>). In addition, specific endogenous retroviral genes have been utilized throughout host evolution to provide protection against exogenous retroviruses (<xref ref-type="bibr" rid="B25">Helmy and Selvarajoo, 2021</xref>; <xref ref-type="bibr" rid="B12">Chen et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B44">Miyake et&#xa0;al., 2022</xref>). In addition to the presence of syncytin genes and receptor interference, which originated tens of millions of years ago, there are many cases of ERV loci with intact gamma-type Env ORFs of unknown purpose (<xref ref-type="bibr" rid="B21">Grandi et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B42">McCarthy et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B49">Simpson et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B28">Hogan and Johnson, 2023</xref>), with some of them remaining unchanged for more than 100 million years (<xref ref-type="bibr" rid="B36">Kasperek et&#xa0;al., 2022</xref>). The extensive tissue expression of mink ERV-V <italic>env</italic> genes suggests that they may have more than just placenta-specific roles (<xref ref-type="bibr" rid="B19">Funk et&#xa0;al., 2019</xref>). It is possible, however, that these genes interfere with RVs that bind the same receptor, as shown by certain genes found in mice and cats that are derived from ERV env (<xref ref-type="bibr" rid="B33">Ito et&#xa0;al., 2013</xref>, <xref ref-type="bibr" rid="B32">Ito et&#xa0;al., 2016</xref>). To explore these possibilities, further investigation into the tertiary structure of mink ERV-V Env and specific cell surface receptors for each protein would be needed. Furthermore, investigation of mink ERV-V Env ORF expression was accomplished by expressing tagged versions of the protein; however, the use of specific antibodies in primary cells from related species may reveal more about its natural expression levels and its location within the cell. In the future, investigation of the role of mink ERV-V Env in testicular regression is warranted.</p>
<p>Cyclical changes in the growth and regression of testes in minks, including seasonal variations in testis size and sperm production, have been well documented (<xref ref-type="bibr" rid="B57">Zhang and Xu, 2022</xref>; <xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>). Notably, the male gonads of minks exhibit distinct characteristics during different periods, including a sexually active phase in March, followed by a period of reduced spermatogenesis and subsequent reactivation of spermatogenesis in late November (<xref ref-type="bibr" rid="B57">Zhang and Xu, 2022</xref>; <xref ref-type="bibr" rid="B58">Zhang et&#xa0;al., 2023</xref>). In the present study, the changes in sperm production in mature minks across seasons were strongly associated with changes in the cellular makeup of testes, as shown in <xref ref-type="fig" rid="f12">
<bold>Figure&#xa0;12A</bold>
</xref>. Hence, we postulated that substantial alterations in the cellular composition of testes serve as causal determinants of the differences in mink ERV-V <italic>env</italic> gene expression. Additionally, spermatozoan generation was significantly correlated with testosterone synthesis, and these factors displayed synchronized seasonal patterns (<xref ref-type="fig" rid="f12">
<bold>Figure&#xa0;12B</bold>
</xref>). It is plausible that these endocrine modifications directly affected activation and deactivation of the mink ERV-V <italic>env</italic> gene, resulting in fluctuations in its expression. More comprehensive investigations are required to address these questions.</p>
<fig id="f12" position="float">
<label>Figure&#xa0;12</label>
<caption>
<p>Fluctuations in the hypothalamus-pituitary-gonadal axis (HPG) lead to variations in testis size and mink ERV-env expression in minks (Neovison vison) throughout the year. <bold>(A)</bold> Schematic representation of the pattern of regulation of the reproductive cycle and mink ERV-env expression in minks. The expression of mink ERV-env in the testes is accompanied by cyclic changes in spermatogenesis. The expression of this gene gradually increases during the regression stage of mink testes, followed by a gradual decline during the transition to the reproductive stage. The diagram illustrates significant changes in the cellular makeup of the testes as a key factor in differences in mink ERV-V env gene expression. <bold>(B)</bold> Seasonal variations in the HPG axis and associated alterations in testicular activity and mink ERV-env expression are observed. Short photoperiods, which are less than 12 hours of light per day, stimulate the hypothalamus to release GnRH, leading to the production and release of FSH and LH from the pituitary gland. Gonadotrophins are responsible for supporting the production of sperms and hormones in the testes. Meanwhile, the mink ERV-V env gene switches can be fully turned off. Extended exposure to light (&gt; 12 hours per day) suppresses the release of GnRH, ultimately halting the generation of both sperms and testosterone. In addition, the mink ERV-V env gene switches can be fully turned on. Thus, we suggest that changes in hormones can regulate the expression of the mink ERV-V env gene, causing it to transition from an active state to a less active state and then to a repressed state.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-14-1404431-g012.tif"/>
</fig>
<p>We have identified a novel retroviral env gene in <italic>Neovison vison</italic>. Through phylogenetic and other analyses of functional motifs specific to RVs or RV subtypes, we characterized the mink ERV-V <italic>env</italic> gene and found it to align with class I ERV genes. Additionally, we observed a gradual increase in mink ERV-V <italic>env</italic> expression during the regression stage of mink testes, followed by a gradual decline during the transition to their reproductive stage. By further investigating the expression and replication of mink ERV-V <italic>env in vitro</italic>, we can gain further valuable insights. In particular, studies to understand the role of mink ERV-V Env in testicular regression are warranted.</p>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="supplementary-material" rid="SF1">
<bold>Supplementary Material</bold>
</xref>. Further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6" sec-type="ethics-statement">
<title>Ethics statement</title>
<p>The animal study was approved by The Institutional Animal Care and Use Committee of Inner Mongolia Agricultural University, China, approved all the protocols utilized in this work (Approval number: 2020007). The study was conducted in accordance with the local legislation and institutional requirements.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>YZ: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing &#x2013; original draft, Writing &#x2013; review &amp; editing. GW: Conceptualization, Data curation, Formal analysis, Methodology, Resources, Writing &#x2013; review &amp; editing. YaZ: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing &#x2013; review &amp; editing. XC: Data curation, Investigation, Methodology, Software, Validation, Writing &#x2013; review &amp; editing. FL: Conceptualization, Data curation, Formal analysis, Investigation, Project administration, Resources, Supervision, Writing &#x2013; review &amp; editing. HL: Formal analysis, Investigation, Methodology, Resources, Writing &#x2013; review &amp; editing. SL: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing &#x2013; review &amp; editing.</p>
</sec>
</body>
<back>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the National Natural Science Foundation of China (Grant numbers 32072819, and 32360863), Science and Technology Major Special Project of Inner Mongolia (Grant numbers 2021ZD0010), Grassland Talents Innovative Team Project of Inner Mongolia (Grant numbers 20151031), the Innovation team project of Education Department of Inner Mongolia (Grant numbers NMGIRT2412), the Scientific Research Innovation Team Project of Inner Mongolia Colleges and Universities (Grant numbers 2200001130) and High-level Talents Introduction to Scientific Research Start-up Project of Inner Mongolia Agricultural University (NDYB2022-6).</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>We would like to express our gratitude to everyone who took part in this study.</p>
</ack>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec id="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors&#xa0;and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s11" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fcimb.2024.1404431/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fcimb.2024.1404431/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Image_1.jpeg" id="SF1" mimetype="image/jpeg">
<label>Supplementary Figure&#xa0;1</label>
<caption>
<p>Hematoxylin and eosin-stained histological sections of the annual seasonal reproductive cycle in the normal adult mink. There are three main phases of testicular function in the circannual cycle of mink: (i) testis recrudescence begins in November; (ii) Mating is observed from late February to late March; and (iii) testicular regression and a long period of sexual rest (May to October).</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Table_1.doc" id="ST1" mimetype="application/msword"/>
<supplementary-material xlink:href="Table_2.doc" id="ST2" mimetype="application/msword"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alser</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Rotman</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Deshpande</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Taraszka</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Baykal</surname> <given-names>P. I.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Technology dictates algorithms: recent developments in read alignment</article-title>. <source>Genome Biol.</source> <volume>22</volume>, <fpage>249</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s13059-021-02443-7</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bamunusinghe</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Plishka</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Dolan</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Skorski</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Oler</surname> <given-names>A. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Recombinant origins of pathogenic and nonpathogenic mouse gammaretroviruses with polytropic host range</article-title>. <source>J. Virol.</source> <volume>91</volume>, <elocation-id>e00855-17</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/JVI.00855-17</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barnes</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Price</surname> <given-names>D. C.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Endogenous viral elements in ixodid tick genomes</article-title>. <source>Viruses</source> <volume>15</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/v15112201</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Beltran-Frutos</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Seco-Rovira</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Martinez-Hernandez</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Ferrer</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Serrano-Sanchez</surname> <given-names>M. I.</given-names>
</name>
<name>
<surname>Pastor</surname> <given-names>L. M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Cellular modifications in spermatogenesis during seasonal testicular regression: an update review in mammals</article-title>. <source>Anim. (Basel)</source> <volume>12</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/ani12131605</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Best</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Le Tissier</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Towers</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Stoye</surname> <given-names>J. P.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>Positional cloning of the mouse retrovirus restriction gene Fv1</article-title>. <source>Nature</source> <volume>382</volume>, <fpage>826</fpage>&#x2013;<lpage>829</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/382826a0</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Blanco-Melo</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Campbell</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Dennis</surname> <given-names>T. P. W.</given-names>
</name>
<name>
<surname>Modha</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Lytras</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>). <article-title>A novel approach to exploring the dark genome and its application to mapping of the vertebrate virus fossil record</article-title>. <source>Genome Biol.</source> <volume>25</volume>, <fpage>120</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s13059-024-03258-y</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Blanco-Melo</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Gifford</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Bieniasz</surname> <given-names>P. D.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Reconstruction of a replication-competent ancestral murine endogenous retrovirus-L</article-title>. <source>Retrovirology</source> <volume>15</volume>, <fpage>34</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12977-018-0416-3</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Blottner</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Schon</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Jewgenow</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Seasonally activated spermatogenesis is correlated with increased testicular production of testosterone and epidermal growth factor in mink (Mustela vison)</article-title>. <source>Theriogenology</source> <volume>66</volume>, <fpage>1593</fpage>&#x2013;<lpage>1598</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.theriogenology.2006.01.041</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Ren</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>S. W.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Syncytin-1 nonfusogenic activities modulate inflammation and contribute to preeclampsia pathogenesis</article-title>. <source>Cell Mol. Life Sci.</source> <volume>79</volume>, <fpage>290</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00018-022-04294-2</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Buzdin</surname> <given-names>A. A.</given-names>
</name>
<name>
<surname>Prassolov</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Garazha</surname> <given-names>A. V.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Friends-enemies: endogenous retroviruses are major transcriptional regulators of human DNA</article-title>. <source>Front. Chem.</source> <volume>5</volume>, <elocation-id>35</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fchem.2017.00035</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Zheng</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Kishima</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Genomic fossils reveal adaptation of non-autonomous pararetroviruses driven by concerted evolution of noncoding regulatory sequences</article-title>. <source>PloS Pathog.</source> <volume>13</volume>, <elocation-id>e1006413</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.ppat.1006413</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Liao</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Song</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y. Y.</given-names>
</name>
<name>
<surname>Cui</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Evolution and genetic diversity of the retroviral envelope in anamniotes</article-title>. <source>J. Virol.</source> <volume>96</volume>, <elocation-id>e0207221</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/jvi.02072-21</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Diehl</surname> <given-names>W. E.</given-names>
</name>
<name>
<surname>Patel</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Halm</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>W. E.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Tracking interspecies transmission and long-term evolution of an ancient retrovirus using the genomes of modern mammals</article-title>. <source>Elife</source> <volume>5</volume>, <fpage>e12704</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.7554/eLife.12704.026</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dittmar</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Weiler</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Hass</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Cell-cell fusion mediated by viruses and HERV-derived fusogens in cancer initiation and progression</article-title>. <source>Cancers (Basel)</source> <volume>13</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/cancers13215363</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dopkins</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Nixon</surname> <given-names>D. F.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Activation of human endogenous retroviruses and its physiological consequences</article-title>. <source>Nat. Rev. Mol. Cell Biol.</source> <volume>25</volume>, <fpage>212</fpage>&#x2013;<lpage>222</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41580-023-00674-z</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Drozdetskiy</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Cole</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Procter</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Barton</surname> <given-names>G. J.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>JPred4: a protein secondary structure prediction server</article-title>. <source>Nucleic Acids Res.</source> <volume>43</volume>, <fpage>W389</fpage>&#x2013;<lpage>W394</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkv332</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feschotte</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Gilbert</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Endogenous viruses: insights into viral evolution and impact on host biology</article-title>. <source>Nat. Rev. Genet.</source> <volume>13</volume>, <fpage>283</fpage>&#x2013;<lpage>296</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nrg3199</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Frank</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Feschotte</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Co-option of endogenous viral sequences for host cell function</article-title>. <source>Curr. Opin. Virol.</source> <volume>25</volume>, <fpage>81</fpage>&#x2013;<lpage>89</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.coviro.2017.07.021</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Funk</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Cornelis</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Vernochet</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Heidmann</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Dupressoir</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Conley</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Capture of a hyena-specific retroviral envelope gene with placental expression associated in evolution with the unique emergence among carnivorans of hemochorial placentation in hyaenidae</article-title>. <source>J. Virol.</source> <volume>93</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/JVI.01811-18</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gifford</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Blomberg</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Coffin</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Fan</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Heidmann</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Mayer</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Nomenclature for endogenous retrovirus (ERV) loci</article-title>. <source>Retrovirology</source> <volume>15</volume>, <fpage>59</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12977-018-0442-1</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grandi</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Pisano</surname> <given-names>M. P.</given-names>
</name>
<name>
<surname>Demurtas</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Blomberg</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Magiorkinis</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Mayer</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Identification and characterization of ERV-W-like sequences in Platyrrhini species provides new insights into the evolutionary history of ERV-W in primates</article-title>. <source>Mob DNA</source> <volume>11</volume>, <fpage>6</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s13100-020-0203-2</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Greenwood</surname> <given-names>A. D.</given-names>
</name>
<name>
<surname>Ishida</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>O&#x2019;brien</surname> <given-names>S. P.</given-names>
</name>
<name>
<surname>Roca</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Eiden</surname> <given-names>M. V.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Transmission, evolution, and endogenization: lessons learned from recent retroviral invasions</article-title>. <source>Microbiol. Mol. Biol. Rev.</source> <volume>82</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/MMBR.00044-17</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gumulka</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Hrabia</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Rozenboim</surname> <given-names>I.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Annual changes in cell proliferation and apoptosis and expression of connexin 43 in the testes of domestic seasonal breeding ganders</article-title>. <source>Theriogenology</source> <volume>186</volume>, <fpage>27</fpage>&#x2013;<lpage>39</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.theriogenology.2022.03.031</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hayward</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Cornwallis</surname> <given-names>C. K.</given-names>
</name>
<name>
<surname>Jern</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Pan-vertebrate comparative genomics unmasks retrovirus macroevolution</article-title>. <source>Proc. Natl. Acad. Sci. U.S.A.</source> <volume>112</volume>, <fpage>464</fpage>&#x2013;<lpage>469</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.1414980112</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Helmy</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Selvarajoo</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Systems biology to understand and regulate human retroviral proinflammatory response</article-title>. <source>Front. Immunol.</source> <volume>12</volume>, <elocation-id>736349</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fimmu.2021.736349</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Henzy</surname> <given-names>J. E.</given-names>
</name>
<name>
<surname>Gifford</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>W. E.</given-names>
</name>
<name>
<surname>Coffin</surname> <given-names>J. M.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>A novel recombinant retrovirus in the genomes of modern birds combines features of avian and mammalian retroviruses</article-title>. <source>J. Virol.</source> <volume>88</volume>, <fpage>2398</fpage>&#x2013;<lpage>2405</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/JVI.02863-13</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Henzy</surname> <given-names>J. E.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>W. E.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Pushing the endogenous envelope</article-title>. <source>Philos. Trans. R Soc. Lond B Biol. Sci.</source> <volume>368</volume>, <fpage>20120506</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1098/rstb.2012.0506</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hogan</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>W. E.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Unique structure and distinctive properties of the ancient and ubiquitous gamma-type envelope glycoprotein</article-title>. <source>Viruses</source> <volume>15</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/v15020274</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Horie</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Tomonaga</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Paleovirology of bornaviruses: What can be learned from molecular fossils of bornaviruses</article-title>. <source>Virus Res.</source> <volume>262</volume>, <fpage>2</fpage>&#x2013;<lpage>9</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.virusres.2018.04.006</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ikeda</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Kato</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kitani</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Suzuki</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Yoshida</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Inaguma</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2001</year>). <article-title>Virological properties and nucleotide sequences of Cas-E-type endogenous ecotropic murine leukemia viruses in South Asian wild mice, Mus musculus castaneus</article-title>. <source>J. Virol.</source> <volume>75</volume>, <fpage>5049</fpage>&#x2013;<lpage>5058</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/JVI.75.11.5049-5058.2001</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ikeda</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Laigret</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Martin</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Repaske</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>1985</year>). <article-title>Characterization of a molecularly cloned retroviral sequence associated with Fv-4 resistance</article-title>. <source>J. Virol.</source> <volume>55</volume>, <fpage>768</fpage>&#x2013;<lpage>777</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/jvi.55.3.768-777.1985</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ito</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Baba</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Kawasaki</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Nishigaki</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Ancestral Mutations Acquired in Refrex-1, a Restriction Factor against Feline Retroviruses, during its Cooption and Domestication</article-title>. <source>J. Virol.</source> <volume>90</volume>, <fpage>1470</fpage>&#x2013;<lpage>1485</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/JVI.01904-15</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ito</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Watanabe</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Hiratsuka</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Kuse</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Odahara</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Ochi</surname> <given-names>H.</given-names>
</name>
<etal/>
</person-group>. (<year>2013</year>). <article-title>Refrex-1, a soluble restriction factor against feline endogenous and exogenous retroviruses</article-title>. <source>J. Virol.</source> <volume>87</volume>, <fpage>12029</fpage>&#x2013;<lpage>12040</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/JVI.01267-13</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jern</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Greenwood</surname> <given-names>A. D.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Wildlife endogenous retroviruses: colonization, consequences, and cooption</article-title>. <source>Trends Genet.</source> <volume>40</volume>, <fpage>149</fpage>&#x2013;<lpage>159</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tig.2023.10.014</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Johnson</surname> <given-names>W. E.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Origins and evolutionary consequences of ancient endogenous retroviruses</article-title>. <source>Nat. Rev. Microbiol.</source> <volume>17</volume>, <fpage>355</fpage>&#x2013;<lpage>370</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41579-019-0189-2</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kasperek</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Beguin</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Bawa</surname> <given-names>O.</given-names>
</name>
<name>
<surname>De Azevedo</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Job</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Massard</surname> <given-names>C.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Therapeutic potential of the human endogenous retroviral envelope protein HEMO: a pan-cancer analysis</article-title>. <source>Mol. Oncol.</source> <volume>16</volume>, <fpage>1451</fpage>&#x2013;<lpage>1473</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/1878-0261.13069</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kassiotis</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>The immunological conundrum of endogenous retroelements</article-title>. <source>Annu. Rev. Immunol.</source> <volume>41</volume>, <fpage>99</fpage>&#x2013;<lpage>125</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev-immunol-101721-033341</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kumar</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Stecher</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Knyaz</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Tamura</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>MEGA X: molecular evolutionary genetics analysis across computing platforms</article-title>. <source>Mol. Biol. Evol.</source> <volume>35</volume>, <fpage>1547</fpage>&#x2013;<lpage>1549</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/molbev/msy096</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Langmead</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Salzberg</surname> <given-names>S. L.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Fast gapped-read alignment with Bowtie 2</article-title>. <source>Nat. Methods</source> <volume>9</volume>, <fpage>357</fpage>&#x2013;<lpage>359</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nmeth.1923</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Ren</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Fan</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>L.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Resurrection of endogenous retroviruses during aging reinforces senescence</article-title>. <source>Cell</source> <volume>186</volume>, <fpage>287</fpage>&#x2013;<lpage>304 e226</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.cell.2022.12.017</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Malfavon-Borja</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Feschotte</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Fighting fire with fire: endogenous retrovirus envelopes as restriction factors</article-title>. <source>J. Virol.</source> <volume>89</volume>, <fpage>4047</fpage>&#x2013;<lpage>4050</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/JVI.03653-14</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McCarthy</surname> <given-names>K. R.</given-names>
</name>
<name>
<surname>Timpona</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Jenni</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Bloyet</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Brusic</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>W. E.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Structure of the receptor binding domain of envP(b)1, an endogenous retroviral envelope protein expressed in human tissues</article-title>. <source>mBio</source> <volume>11</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/mBio.02772-20</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Meyer</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Rosenkrantz</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Carbone</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Chavez</surname> <given-names>S. L.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Endogenous retroviruses: with us and against us</article-title>. <source>Front. Chem.</source> <volume>5</volume>, <elocation-id>23</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fchem.2017.00023</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Miyake</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Ngo</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>Wulandari</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Shimojima</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Nakagawa</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Kawasaki</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Convergent evolution of antiviral machinery derived from endogenous retrovirus truncated envelope genes in multiple species</article-title>. <source>Proc. Natl. Acad. Sci. U.S.A.</source> <volume>119</volume>, <elocation-id>e2114441119</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.2114441119</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Naville</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Volff</surname> <given-names>J. N.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Endogenous retroviruses in fish genomes: from relics of past infections to evolutionary innovations</article-title>? <source>Front. Microbiol.</source> <volume>7</volume>, <elocation-id>1197</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmicb.2016.01197</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pandey</surname> <given-names>S. P.</given-names>
</name>
<name>
<surname>Mohanty</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Role of the testicular capsule in seasonal modulation of the testis</article-title>. <source>J. Exp. Zool A Ecol. Integr. Physiol.</source> <volume>339</volume>, <fpage>898</fpage>&#x2013;<lpage>910</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/jez.2740</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pelletier</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Akpovi</surname> <given-names>C. D.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Kumar</surname> <given-names>N. M.</given-names>
</name>
<name>
<surname>Vitale</surname> <given-names>M. L.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Complementary expression and phosphorylation of Cx46 and Cx50 during development and following gene deletion in mouse and in normal and orchitic mink testes</article-title>. <source>Am. J. Physiol. Regul. Integr. Comp. Physiol.</source> <volume>309</volume>, <fpage>R255</fpage>&#x2013;<lpage>R276</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1152/ajpregu.00152.2015</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Plant</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Bellefroid</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Van Lint</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A complex network of transcription factors and epigenetic regulators involved in bovine leukemia virus transcriptional regulation</article-title>. <source>Retrovirology</source> <volume>20</volume>, <fpage>11</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12977-023-00623-w</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Simpson</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Kozak</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Boso</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Cross-species transmission of an ancient endogenous retrovirus and convergent co-option of its envelope gene in two mammalian orders</article-title>. <source>PloS Genet.</source> <volume>18</volume>, <elocation-id>e1010458</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pgen.1010458</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Skirmuntt</surname> <given-names>E. C.</given-names>
</name>
<name>
<surname>Katzourakis</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>The evolution of endogenous retroviral envelope genes in bats and their potential contribution to host biology</article-title>. <source>Virus Res.</source> <volume>270</volume>, <fpage>197645</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.virusres.2019.197645</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tao</surname> <given-names>C. Y.</given-names>
</name>
<name>
<surname>Harley</surname> <given-names>J. Z.</given-names>
</name>
<name>
<surname>Spencer</surname> <given-names>S. L.</given-names>
</name>
<name>
<surname>Cohen</surname> <given-names>R. E.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Characterizing seasonal transitions: Breeding-like morphology and behavior during the late non-breeding season in green anole lizards</article-title>. <source>Horm. Behav.</source> <volume>139</volume>, <fpage>105106</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.yhbeh.2021.105106</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>G. H.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Endogenous retroviruses in development and health</article-title>. <source>Trends Microbiol.</source> <volume>32</volume>, <fpage>342</fpage>&#x2013;<lpage>354</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tim.2023.09.006</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Warnow</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Revisiting evaluation of multiple sequence alignment methods</article-title>. <source>Methods Mol. Biol.</source> <volume>2231</volume>, <fpage>299</fpage>&#x2013;<lpage>317</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/978-1-0716-1036-7_17</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Yan</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Roy</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Poisson</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>The I-TASSER Suite: protein structure and function prediction</article-title>. <source>Nat. Methods</source> <volume>12</volume>, <fpage>7</fpage>&#x2013;<lpage>8</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nmeth.3213</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yedavalli</surname> <given-names>V. R. K.</given-names>
</name>
<name>
<surname>Patil</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Parrish</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Kozak</surname> <given-names>C. A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A novel class III endogenous retrovirus with a class I envelope gene in African frogs with an intact genome and developmentally regulated transcripts in Xenopus tropicalis</article-title>. <source>Retrovirology</source> <volume>18</volume>, <fpage>20</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12977-021-00564-2</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Recent advances in the study of active endogenous retrovirus envelope glycoproteins in the mammalian placenta</article-title>. <source>Virol. Sin.</source> <volume>30</volume>, <fpage>239</fpage>&#x2013;<lpage>248</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12250-015-3617-0</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Involvement of testicular N-glycoproteome heterogeneity in seasonal spermatogenesis of the American mink (Neovison vison)</article-title>. <source>Front. Vet. Sci.</source> <volume>9</volume>, <elocation-id>1001431</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fvets.2022.1001431</pub-id>
</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Cao</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Cell adhesion function was altered during the seasonal regression of the seminiferous epithelium in the mink species Neovison vison</article-title>. <source>J. Anim. Sci.</source> <volume>101</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/jas/skad190</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Sam</surname> <given-names>T. W.</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>A. Y.</given-names>
</name>
<name>
<surname>Leung</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Mouse strain-specific polymorphic provirus functions as cis-regulatory element leading to epigenomic and transcriptomic variations</article-title>. <source>Nat. Commun.</source> <volume>12</volume>, <fpage>6462</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41467-021-26630-z</pub-id>
</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhuo</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Feschotte</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Cross-species transmission and differential fate of an endogenous retrovirus in three mammal lineages</article-title>. <source>PloS Pathog.</source> <volume>11</volume>, <elocation-id>e1005279</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.ppat.1005279</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>