Next Article in Journal
Prevalence of Emergent Dolutegravir Resistance Mutations in People Living with HIV: A Rapid Scoping Review
Next Article in Special Issue
The Accumulation of Phenyllactic Acid Impairs Host Glutamine Metabolism and Inhibits African Swine Fever Virus Replication: A Novel Target for the Development of Anti-ASFV Drugs
Previous Article in Journal
Correction: Cook et al. An Optimized Bioassay for Screening Combined Anticoronaviral Compounds for Efficacy against Feline Infectious Peritonitis Virus with Pharmacokinetic Analyses of GS-441524, Remdesivir, and Molnupiravir in Cats. Viruses 2022, 14, 2429
Previous Article in Special Issue
Immune Responses to Influenza D Virus in Calves Previously Infected with Bovine Viral Diarrhea Virus 2
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 16
Issue 3
10.3390/v16030398
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolution of Endogenous Retroviruses in the Subfamily of Caprinae

by
Ali Shoaib Moawad
1,2,†,
Fengxu Wang
3,†,
Yao Zheng
1,
Cai Chen
1,
Ahmed A. Saleh
1,4,
Jian Hou
3 and
Chengyi Song
1,*
1
College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China
2
Department of Animal Production, Faculty of Agriculture, Kafrelsheikh University, Kafrelsheikh 33516, Egypt
3
State Key Laboratory of Animal Biotech Breeding, College of Biological Sciences, China Agricultural University, Beijing 100193, China
4
Animal and Fish Production Department, Faculty of Agriculture (Alshatby), Alexandria University, Alexandria City 11865, Egypt
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2024, 16(3), 398; https://doi.org/10.3390/v16030398
Submission received: 30 January 2024 / Revised: 28 February 2024 / Accepted: 1 March 2024 / Published: 4 March 2024
(This article belongs to the Special Issue Advances in Endemic and Emerging Viral Diseases in Livestock)
Browse Figures
Review Reports Versions Notes

Abstract

:
The interest in endogenous retroviruses (ERVs) has been fueled by their impact on the evolution of the host genome. In this study, we used multiple pipelines to conduct a de novo exploration and annotation of ERVs in 13 species of the Caprinae subfamily. Through analyses of sequence identity, structural organization, and phylogeny, we defined 28 ERV groups within Caprinae, including 19 gamma retrovirus groups and 9 beta retrovirus groups. Notably, we identified four recent and potentially active groups prevalent in the Caprinae genomes. Additionally, our investigation revealed that most long noncoding genes (lncRNA) and protein-coding genes (PC) contain ERV-derived sequences. Specifically, we observed that ERV-derived sequences were present in approximately 75% of protein-coding genes and 81% of lncRNA genes in sheep. Similarly, in goats, ERV-derived sequences were found in approximately 74% of protein-coding genes and 75% of lncRNA genes. Our findings lead to the conclusion that the majority of ERVs in the Caprinae genomes can be categorized as fossils, representing remnants of past retroviral infections that have become permanently integrated into the genomes. Nevertheless, the identification of the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 groups indicates the presence of relatively recent and potentially active ERVs in these genomes. These particular groups may contribute to the ongoing evolution of the Caprinae genome. The identification of putatively active ERVs in the Caprinae genomes raises the possibility of harnessing them for future genetic marker development.

1. Introduction

Endogenous retroviruses (ERVs) are a type of genomic element that are found across a wide range of vertebrates, spanning from basal vertebrates like sharks and rays to mammals [1]. Throughout the process of evolution, exogenous retroviruses have integrated into the germ line, leading to the formation of stably integrated endogenous retroviruses. These endogenous retroviruses are subsequently inherited by offspring through vertical transmission. ERVs typically consist of three primary genes: Group-specific antigen (gag), which instructs the synthesis of CAPSID proteins; Protease–Polymerase (pro-pol), responsible for producing enzymes involved in maturation, replication, Retroviral integrase (IN), which inserts the viral DNA into the host chromosomal DNA, and insertion; and Envelope (env), where the expression of the env gene enables retroviruses to target and attach to specific cell types and to infiltrate the target cell membrane. These genes are flanked by long terminal repeats (LTRs), which are regulatory regions that contain promoters, enhancers, and polyadenylation signals [1]. The relationship between endogenous and exogenous viral sequences is presumed to be by evolutionary descent, and the presumption is that endogenous viral sequences are derived from exogenous viruses that have infected the host organism’s ancestors in the past and become integrated into their genome through evolutionary processes. ERVs can be classified into three different classes based on their relationships to recognized retroviral genera [2]: Class I ERVs are associated with Epsilon retrovirus and Gamma retrovirus, Class II ERVs with Alpha retrovirus, Beta retrovirus, Delta retrovirus, and Lenti retrovirus, and Class III ERVs with Spuma retrovirus [1].
The interest in these genomic elements has grown significantly due to their potential role in diseases like multiple sclerosis and cancer [3]. Additionally, interest in these genomic elements has been fueled by their impact on the evolution of the host genome, including changes in nearby gene expression, among other factors [4]. In addition, ERVs can domesticate in host genomes and evolve into new functional genes, such as Syncytins (Syncytin-1 and Syncytin-2), which are derived from the Env genes of ERVs, and they play essential roles in placental development [5]. The implications of ERVs’ presence in host genomes have been re-evaluated [6,7]. Since the late 1990s, researchers have investigated the presence of ERVs in certain domestic animals using Southern blot analysis [8] and Polymerase Chain Reaction (PCR) analyses [9]. These studies identified the presence of retroviruses in various domestic animals, including pigs, horses, sheep, goats, cattle, yaks, and cats [8,9]. Nevertheless, the comprehensive investigation of ERVs has typically been focused on humans [10] and murids [11]. In recent years, there has been an increasing trend of analyzing ERVs in a broader range of domestic species.
The Caprinae subfamily is a subset of the Bovidae family of ruminants, primarily consisting of medium-sized bovids. Notable members include sheep, goats, and their relatives. This particular subfamily comprises 12 genera; however, the classification of Caprinae is intricate, resulting in various suggested classifications. To address these ongoing taxonomy issues, the International Union for Conservation of Nature (IUCN) has established a Taxonomy Working Group as part of the Caprinae Specialist Group. Caprinids possess distinct adaptations for montane and alpine environments, which is why this subfamily displays more diversity in Eurasia compared to Africa [12,13,14,15]. Many species have become extinct since the last ice age, predominantly due to human influence. Among the remaining species, five are classified as endangered, eight as vulnerable, and seven as requiring conservation measures due to concerns, but at lower risk [16,17,18,19]. The polymorphisms generated by endogenous Jaagsiekte Sheep Retrovirus (enJSRV) have been used as highly informative genetic markers because the presence of each endogenous retrovirus in the host genome is the result of a single integration event in a single animal and is irreversible, so populations sharing the same provirus in the same genomic location are de facto phylogenetically related [20].
Sheep (Ovis aries) have served as valuable models for studying the coevolution between a host and ERVs, specifically the Jaagsiekte Sheep Retrovirus (JSRV) and its endogenous counterpart, enJSRV. The significance of this model lies in the fact that both forms of the retrovirus are actively present and exhibit strong interactions [21]. The remaining ERVs within the sheep genome can be categorized into nine families from Class I (OERV γ1 to γ9) and three families from Class II (OERV β1 to β3) based on the amplification of their pro/pol complex [22]. Given the significant role of enJSRVs in sheep development and placental morphogenesis, and their ability to block the exogenous JSRV, the expression analyses of sheep ERVs have primarily focused on the enJSRV provirus [21].
Recently, the investigation of endogenous retroviruses (ERVs) has extended to both wild and domestic animal populations. However, our understanding of ERV annotation and evolution in livestock genomes remains limited. To address this knowledge gap, our study focuses on the Caprinae subfamily reference genomes, where we utilize multiple pipelines to de novo mine and annotate ERVs. We systematically characterize their classification, coverage, and evolutionary dynamics within these genomes. Furthermore, we investigated the intersection between ERVs and host genes, including protein-encoding and lncRNA genes. Our comprehensive analysis reveals a genome-wide landscape of ERVs across 13 different species within the Caprinae subfamily. This research not only contributes to our understanding of the coevolution between the mobilome and livestock genome but also enhances our comprehension of genome evolution as a whole.

2. Materials and Methods

2.1. Data Sources of Genomes and ERV Reference Sequences Used in This Study

The reference genomes of 13 species belonging to the Caprinae subfamily were obtained from the National Centre for Biotechnology Information (NCBI) database (https://www.ncbi.nlm.nih.gov, accessed on 14 April 2023) and are listed in Table S1 (Supplementary File S1). To define the classification of Caprinae ERVs, the reference reverse transcriptase (RT) DNA sequences of retroviruses from species other than Caprinae and the previously identified retroviruses in the Caprinae subfamily (containing the reference enJSRVs in sheep) were included (Table S2, Supplementary File S1), in addition to the reference DNA sequences of major homology region (MHR) of gag. These sequences were downloaded from the NCBI nucleotide database. The GenBank accession numbers of the ERVs used for phylogenetic analysis are listed in Table S3 (Supplementary File S1).

2.2. ERVs Mining in the Caprinae Genomes

ERVs were identified from Caprinae genomes using the LTRharvest [23] tool embedded in GenomeTools [24] using the parameters-motif TGCA-minlenltr 100-maxlenltr 7000-similar 90-motifmis 1-mintsd 4-maxtsd 6-overlaps best and RepeatModeler (version 2.0) (http://www.repeatmasker.org/RepeatModeler, accessed on 31 May 2023) using the parameters BuildDatabase-name Genome. LTRharvest is a specialized tool that enables the de novo detection of full-length LTR retrotransposons, including ERVs [23]. To extract the sequences of the genes and their flanking regions from the genomes based on the identified coordinates, the bedtools getfasta tool was employed, bedtools (version 2.27.1) [25]. This tool allows for obtaining sequences from a genome using specified coordinates from a bed file. Sequences that were smaller than 4000 bp or larger than 10,000 bp were excluded from further analysis. Only potential intact sequences within this size range were retained. To obtain full-length copies of the potential intact sequences, about 4 kb upstream and downstream flanking regions were added. This extension was achieved using the slop command (−s, −l 4000, −r 4000) in bedtools.
To further filter and analyze the identified ERVs, only those with potential intact reverse transcriptase (RT) regions were retained. For this purpose, a search for reverse transcriptase sequences was conducted in the sequences detected by LTRharvest and RepeatModeler. The search was performed using the hmmsearch command within the HMMER tool (https://www.ebi.ac.uk/Tools/hmmer/search/hmmscan, accessed on 12 September 2023) [26]. To identify conserved reverse transcriptase sequences, the sequences of well-annotated reverse transcriptase regions in different mammals (specially Caprinae) were downloaded from NCBI and the Pfam [27] and then were converted to Hidden Markov Model (HMM) profiles using the hmmbuild command in HMMER tool. The Pfam database contains Hidden Markov Model (HMM) profiles for protein domains and families. HMM profiles are statistical models that represent multiple sequence alignments and capture conserved motifs and patterns within a protein family. The profiles downloaded from NCBI and the Pfam are specifically designed for reverse transcriptase regions and are used to identify similar sequences in the input data.
To cluster the ERV elements based on sequence similarity and alignment, the Vsearch tool (https://github.com/torognes/vsearch, accessed on 1 October 2023) was utilized. The clustering was performed using the parameters: -id 0.50 to set the similarity threshold to 50%, -clusters to create the clusters, and -msaout to obtain multiple sequence alignments as output. From the clustered sequences, representative sequences were selected for each group based on the integrity of sequences containing potential intact LTRs and the availability of different domains (gag, pro, pol, IN, and env). This process involved analyzing the relationships between the sequences and selecting a representative sequence that best represents the characteristics of each group.

2.3. ERVs Annotation in the Caprinae Genomes

RepeatMasker (version 4.0.9) software [28] was used as a sequence search engine to annotate ERVs in the Caprinae genomes using a custom library, which combined the known repeats in sheep and goat genomes from the Repbase [29]. To identify the functional domains in all predicted ERV groups within the Caprinae subfamily genomes, we translated all six frames of all sequences in each group. This step accounts for the possibility of different reading frames in the sequences. The putatively endogenous retroviruses (ERVs) containing complete retroviral functional domains (gag, pro, pol, IN, and env) were translated by Genscan (http://hollywood.mit.edu/GENSCAN.html, accessed on 11 October 2023). The translated sequences were then searched against the Pfam database using the hmmsearch command [30]. By comparing the translated sequences to the Pfam database, we were able to identify the functional domains present in the predicted ERV groups. To search for retroviral genes in the Caprinae subfamily genomes, we chose to utilize HMMER. We obtained HMM profiles specific to retroviral genera such as Gamma retrovirus, Epsilon retrovirus, Alpha retrovirus, Beta retrovirus, Delta retrovirus, Lenti retrovirus, Spuma retrovirus, and Gypsy virus. Additionally, we retrieved the sequences for specific genes or regions within retroviruses, including gag, pro, pol, IN, and env. These sequences were obtained from the NCBI Database (https://www.ncbi.nlm.nih.gov, accessed on 15 October 2023) and converted to HMM profiles using the hmmbuild command in the HMMER tool.

2.4. Phylogenetic Analysis

We generated phylogenetic trees based on alignments of two regions of the genome: the DNA sequences of reverse transcriptase (RT) of pol and the major homology region (MHR) of gag. To construct multiple alignments of these regions from ERV retrotransposons, we utilized the MAFFT program v.7.310 [31]. We opted to use DNA sequences instead of amino acid sequences for constructing multiple alignments and building the phylogenetic trees due to the presence of frame-shift mutations and stop codons in the ancient retrotransposon elements. These mutations can disrupt the reading frame and introduce premature stops, rendering the translation of these sequences into amino acids unreliable. To build the phylogenetic tree, we used the IQ-TREE program [32] and the parameters-m TEST-bb 1000-nt AUTO, which is a software tool specifically designed for phylogenetic analysis.

2.5. Divergence Analysis of Full ERVs

To identify any “active” elements of ERVs, each, of 13, reference genome of Caprinae species was masked using RepeatMasker (version 4.0.9, -nolow) [28] with the custom repeat library of the full-length ERV elements generated in Section 2.3. Then, the Kimura two-parameter distance (divergence) was determined using RepeatMasker’s calcDivergenceFromAlign.pl package from the RepeatMasker software [28].

2.6. PCR Verification

To identify the potential retrotransposon insertion polymorphisms (RIPs) in Caprinae, we mapped the upstream and downstream flanking sequences (500 bp) in the four full-length ERV groups to the reference genome using Blat [33]. The differential insertions, designated as putative ERV insertion polymorphisms, were obtained using a bedtools window (-w 50, -v) [25]. The ERV insertions that did not fall into the same window (ERV insertion site and 400 bp flanking region) were considered to be putative ERV RIPs.
Six breeds of sheep (Barki, Rahmani, Rahmani x Barki crossbred, Awassi, Ossimi, and Sulffok) were used for PCR verification of full-length ERV polymorphisms. From each sheep breed, five individual DNA samples were isolated from blood using a DNA isolation-kit (Tiangen Biotech, Beijing, China) and then pooled together and subjected to PCR analysis in order to detect the predicted ERV insertion polymorphisms. Considering that most full-length ERVs have a length exceeding 8 kb, double PCR primer pairs were designed and utilized for genotype verification. The strategy employed for primer design is illustrated in Figure 1. The LTR and side primers for ERV-RIPs were designed using Oligo v.7 software and are presented in Supplementary Table S4. The primers were synthesized by Vazyme Biotech. Co., Ltd., Nanjing, China. Amplification was performed using Green-SuperMix from TaKaRa, Japan, with 10 pM of primers and 100 ng of genomic DNA per sample. The PCR consisted of 34 cycles, including denaturation at 95 °C for 1 min, annealing at 58–60 °C for 1 min, extension at 72 °C for 1 min, and a final extension step at 72 °C for 2 min. The amplification process was carried out using T100 Thermal Cycler (BIO RAD, Singapore).

2.7. Genotyping Analysis

We employed double PCR primer pairs (LTR_Primer and Side_Primer) for each potential ERV-RIP sequence. After PCR amplification, if the LTR primer showed a band between 250 and 750 bp but the Side_Primer, which produces a band within a size ranging between 1000 and 1500 bp, was absent, it indicated homozygous ERV insertion and was labeled as ERV+/+. On the contrary, if the longer band within the size range 1000–1500 bp appeared, and no short band (250–750 bp) was present, it suggested an absence of homozygous ERV insertion genotype (labeled as ERV−/−). If both shorter and longer bands were observed, it signified that this RIP genotype must be heterozygous (labeled as ERV+/−). The methodology employed in this part of the study was based on the framework proposed by Du et al. [34].

2.8. Intersection Analysis

Using the intersect program in bedtools (version 2.27.1) [25], we analyzed the distribution bias of retrotransposons in sheep and goat genomes and their intersections with host genes, including protein-encoding genes and lncRNA genes. The bedtools intersect makes it possible to examine overlaps between two sets of genomic features. Moreover, it allows one to have fine control as to how the intersections are reported. The genes were obtained from the annotation files of sheep and goats, downloaded from the NCBI database (https://www.ncbi.nlm.nih.gov, accessed on 19 November 2023), and intersected with the results of the RepeatMasker (version 4.0.9) software [28]. RepeatMasker was used to annotate all ERVs in the sheep and goat genomes using a custom library of the 28 ERV elements identified in Section 2.2.

3. Results

3.1. ERV Mining in the Genomes of Caprinae

ERV-derived elements were extracted from 13 genomes of the Caprinae subfamily using the LTRharvest and RepeatModeler pipelines. The number of ERVs obtained varied across species, with LTRharvest outputs between 79,529 in Capra aegagrus and 215,180 in Ovis aries, and RepeatModeler ranging between 118,052 in Ovis aries and 211,456 in Capra hircus (Table S5, Supplementary File S1). ERV-derived elements that lacked the reverse transcriptase (RT) regions were then excluded. ERVs with potential intact RT regions ranged from 2503 in Capra ibex to 6760 in Capra hircus. When considering the number of full-length ERV copies, including gag, pro, pol, IN, and env domains + LTRs, variations were observed among the different genera of the Caprinae subfamily. A total of 218, 131, 71, 119, and 106 full-length ERVs were detected in the species O. aries, Ovis orientalis, Ovis ammon, Ovis canadensis, and Ovis nivicola of the Ovis genus, respectively. In the Capra genus, 226, 96, 57, and 100 full-length ERVs were identified in Capra hircus, Capra aegagrus, Capra ibex, and Capra sibirica, respectively. The genera of relatives yielded 97 full-length ERVs in Hemitragus hylocrius, 65 in Budorcas taxicolor, 70 in Oreamnos americanus, and 126 in Ovibos moschatus (Table S6, Supplementary File S1).

3.2. Classification of ERVs in the Genomes of Caprinae

The number of ERVs containing reverse transcriptase (RT) regions varied significantly across different types of retroviruses and groups. By utilizing the sequences and constructing a phylogenetic tree, we were able to accurately assess the evolutionary relationships and genetic distances among the ERV retrotransposon elements in the Caprinae genomes. Through phylogenetic analysis, we identified 28 groups of Caprinae ERVs (Cap_ERV_1-28), 22 new groups, and 6 previously reported ERV groups (Table 1). Gamma retroviruses generally had a higher number of ERV elements compared to beta retroviruses. Specifically, 19 groups were classified as Class 1 gamma retroviruses, while 9 groups were classified as Class 2 beta retroviruses, together with several reference sequences of enJSRVs in sheep (Figure 2A). The same classification was detected using the major homology region (MHR) of gag from the representative sequences of groups (Figure 2B). The different genes and their domains located in 28 groups of Caprinae ERVs (Cap_ERV_1-28) are presented in Table S7 (Supplementary File S1).
The lengths of the Cap_ERV elements ranged from 4421 to 9580 bp, with the length of the long terminal repeat (LTR) varying from 143 to 583 bp (Table 1). To identify the functional domains within each Cap_ERV element, we translated all sequences in each group through all six frames and searched against the Pfam database using hmmsearch [30]. We successfully identified multiple retrovirally relevant domains. ERV candidates featuring two LTRs and structural polyproteins common to all retroviruses, namely group-specific antigen (gag), protease (Pro), polymerase (pol), integrase (IN), and envelope (env), were classified as full-length ERVs. Most of the ERV groups in the Caprinae subfamily had undergone decay and were presumed to be inactive.
Only four groups (Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25) of ERV candidates, classified as beta retroviruses, were found to process full-length ERVs, with all retroviral functional proteins (gag, pro, pol, IN, and env) being detectable (Table 1). They may potentially retain the functionality of ERVs and, thus, were designated as putatively active ERVs. A total of 1482 full-length copies of potentially active ERVs were identified in 13 species of Caprinae subfamily (Table S6, Supplementary File S1), distributed in groups as follows: 418 in Cap_ERV_20, 331 in Cap_ERV_21, 479 in Cap_ERV_24, and 254 in Cap_ERV_25 (Table 2). Moreover, a primer binding site (PBSLys, PBSIle, PBSLys, and PBSHis) of 18 nucleotides downstream and adjacent to the 5′-LTR was identified in the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25, respectively. Additionally, a small region of polypurine tract (PPT) was found upstream and adjacent to the 3′-LTR. The structure organization of putatively active ERVs detected in the four full-length ERV elements is summarized in Table S6 and Figure S1 (Supplementary Files S1 and S3).

3.3. Distribution and Genomic Coverages of ERVs in Different Species of Caprinae

In this study, the analysis using LTRharvest and RepeatModeler identified a total of 28 Cap_ERV groups present in each one of all species analyzed, indicating their presence across all members of the Caprinae subfamily (Figure 3A). Among these groups, 19 were classified as Class I ERVs, while the remaining 9 were classified as Class II ERVs. This suggests that the Cap_ERV elements were inserted into the genome of the common ancestor of the Caprinae subfamily around 11 million years ago (MYA) or possibly even earlier. The abundance of detected copies for each Cap_ERV group varied, and this information is summarized in Table S5, (Supplementary File S1). Among all the groups, the ones with the highest abundance, determined by the number of matches or mapped reads, were Cap_ERV_7 in Ovis aries (sheep), Cap_ERV_18 in Capra hircus (goat), Cap_ERV_13 in Budorcas taxicolor, and Cap_ERV_19 in the remaining species. The analysis of a single reference genome for each of the 13 species in the Caprinae subfamily reveals variation in the numbers of endogenous retroviruses (ERVs) annotated by RepeaMasker (Figure 3B and Table S5). Furthermore, there is variability in the proportion of endogenous retroviruses (ERVs) about the total amount of retrotransposons. The Budorcas taxicolor genome has the highest proportion of ERVs at 10.63%. In contrast, Capra ibex has the lowest proportion of ERVs in its genome, with a value of 6.89% (Figure 3C).

3.4. Kimura Divergence and Genotyping Analysis of Recent and Potentially Active ERVs in the Genomes of Caprinae

The analysis of ERV divergence provides insights into the evolutionary history of ERVs. Our findings indicate that out of all the ERV candidates analyzed, only four were identified as full-length ERVs. These four ERVs belong to the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 groups of beta retroviruses (Figure 2A and Table 1). Kimura (K) divergence analysis revealed that the predominant recent activities of the four full-length ERV elements were observed for all 13 genomes of Caprinae (Figure S2, Supplementary File S3), with some ERV copies of these elements representing K divergence less than 5%. However, a few copies of putatively active ERVs display recent activity and represent less than 2% of K divergence, suggesting they are recent invaders of the genome, and may potentially retain functionality.
Four insertion polymorphisms in each group of full-length ERVs in the reference genome of sheep were evaluated by PCR using two pairs of primers (LTR_Primer and Side_Primer) with genomic DNA samples from six sheep breeds, as described in the methods, for proving the activity, and the genotyping results are shown in Figure 4. From the electrophoretogram results, we found that 15, 12, 16, and 6 of ERV genotypes in sheep were ERV+/−, while ERV+/+ represented 8, 12, 3, and 7 genotypes in Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25, respectively. Regarding the ERV−/− genotype, it was observed only in Cap_ERV_25 (10 genotypes); the different genotypes are presented in Table S8 (Supplementary File S1). Ultimately, the four full-length ERV elements were confirmed by PCR, indicating that full-length ERV insertions are highly polymorphic, as they were genotyped in a small population (30 individuals) of studied sheep breeds, and all showed polymorphism (as shown in Figure 4 and Supplementary Table S8).

3.5. Most lncRNAs and Protein-Coding Genes Contain ERV-Derived Sequences in Sheep and Goats

The intersection analysis revealed that most lncRNAs and protein-coding genes contain ERV-derived sequences in sheep and goats. Specifically, approximately 75.26% (15,999) of protein-coding genes and 80.66% (3583) of lncRNA genes in sheep contained ERV-derived sequences. Similarly, in goats, 73.64% (15,150) of protein-coding genes and 75.47% (1471) of lncRNA genes contained ERV-derived sequences (Figure 5A). On the other hand, in sheep, it was found that 24.22% (743,666) of ERV insertions (3,008,458) overlapped with protein-coding genes, while 1.63% (48,955) overlapped with lncRNA genes. Furthermore, 0.186% (5607) of ERV insertions were found to overlap with both protein-coding and lncRNA genes. Similarly, in goats, 22.6% (691,235) of ERV insertions (3,057,947) overlapped with protein-coding genes, and 0.75% (22,957) overlapped with lncRNA genes (Figure 5B). These findings indicate that ERV-derived sequences contribute to the evolution of both lncRNA and protein-coding genes. In fact, ERVs themselves may be lncRNAs through transcription from their 3’LTR.
Concerning of the ERV coverages of various genic features, it revealed noticeable biases in their ERV composition. In sheep, a significant depletion of ERVs was observed in transcribed regions, including PCgene, mRNA, transcript, pseudogene, and exons, with a sequence coverage of approximately 5% or lower. However, no significant bias was observed in the distribution of ERVs in the lnc_RNA gene (19.36%), cDNA_match (53.15%), coding sequence CDS (29.10%), and tRNA (37.68%). In goats, the ERV coverages varied for different genic features, with percentages as follows: PCgene (17.08%), lnc_RNAgene (19.92%), mRNA (8.85%), transcript (13.25%), pseudogene (26.55%), exons (32.72%), cDNA_match (21.30%), coding sequence CDS (40.05%), V_gene_segment (28.81%), C_gene_segment (16.19%), and tRNA (35.18%), as shown in Table S9 (Supplementary File S1).

4. Discussion

The mammalian genomes are mainly composed of three major classes of retrotransposons: long interspersed nuclear elements (LINEs), short interspersed nuclear elements (SINEs), and LTR-retrotransposons. ERVs are the major retrotransposons within mammalian genomes, accounting for 5–10% of the genomic sequences [35,36,37,38]. In this study, we conducted a comprehensive analysis to identify and characterize ERV groups in 13 single reference genomes belonging to 13 different species of the Caprinae subfamily, noting variations in their domain fractions among different species. To the best of our knowledge, this is the first study to provide genome-wide predictions of ERVs at the domain level across a wide range of Caprinae species. ERVs, being the primary long terminal repeat (LTR) retrotransposons found in mammalian genomes, have been integrated into mammalian lineages for over 100 million years [39]. Here, we observed that ERVs exhibit varying genomic coverages within the Caprinae subfamily. The range of genomic coverage for ERVs in the Caprinae subfamily was found to be between 6.89% in the genome of Capra ibex and 10.63% in the genome of Budorcas taxicolor. These findings are consistent with what has been observed in the genomes of other mammalian species, where ERVs also tend to occupy a substantial proportion of the genome [36,40,41,42,43].
Within the species examined in this study, we identified thousands of ERV candidates. It was determined that the majority of these ERVs had undergone decay. However, among these candidates, a subset was found to contain intact reverse transcriptase (RT) regions. The number of ERV candidates with intact RT regions varied from 2503 in Capra ibex to 6760 in Capra hircus. Through phylogenetic analysis, these RT-motif-containing candidates were classified into 19 gamma ERVs and 9 beta ERVs similar to the reference enJSRVs in sheep. The same classification of ERV groups was detected using the major homology region (MHR). The gag gene is poorly conserved among ERVs, but all detected Cap_ERV groups contain an MHR. The retroviral Gag protein exhibits extensive amino acid sequence variation overall; however, one region of Gag, termed the major homology region (MHR), is conserved among all retroviruses [44], which is consistent with our results. Based on divergence analysis of four full-length ERVs among the 13 species of Caprinae (Figure S2, Supplementary File S3), some copies of these elements exhibited an extended period of expansion and indicated increased activity in recent years. It is noteworthy that a few copies of these groups encoded long peptides with intact gag, pro, pol, IN, and env domains. These findings align with similar observations made in the pig genome, further supporting our results [40,45]. Based on our findings, it can be concluded that the majority of ERVs in the genomes of Caprinae are considered to be fossils, meaning they are remnants of ancient retroviral infections that have become permanently integrated into the genomes. However, the presence of the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 elements suggest the existence of potentially active ERVs. This may play a role in the ongoing evolution of the Caprinae genome [46,47].
The presence of full-length ERV insertion polymorphisms, which typically consist of large structural variations exceeding 5 kb, can result in a stronger genetic impact and can be considered as an indicator of activity [34]. In this study, we performed PCR verification and genotyping analysis on six sheep breeds, with five individuals from each breed, using the primers of four RIPs from each group of full-length ERVs. We obtained insertion polymorphic sites. Subsequently, we confirmed that these insertions belong to the four full-length ERVs. This finding indicates that certain full-length ERVs still retain the capability to mobilize and influence the sheep genomes, some of which may contribute to phenotypic variations and population differentiation, despite the majority of ERVs in sheep genomes being considered as fossils.
Early studies on genome sequencing have revealed variations in the activity of retroviruses across mammalian species. Mice contain numerous active ERV elements [41], while several domesticated animals, including pigs [40], rabbits [48], and domestic cats [49], reported previously, and goat and sheep reported here, were found to primarily possess inactive ERV elements, but few of these copies are putative active. These differences in retroviral activity among species contribute to our understanding of the evolutionary dynamics and potential impacts of ERVs in different genomes of mammals [36]. We observed an extensive intersection of ERV-derived sequences with genes in the genomes of sheep and goats. In sheep, ERV-derived sequences were found in approximately 75% of protein-coding genes and 81% of lncRNA genes. Similarly, in goats, ERV-derived sequences were present in about 74% of protein-coding genes and 75% of lncRNA genes, indicating that ERVs substantially contribute to the evolution of both genes and genomes in sheep and goats, exemplifying a characteristic feature of mammals [36,41,50].

5. Conclusions

In the present study, de novo mining of ERVs in 13 species of the Caprinae subfamily was performed, allowing us to define their diversity at a group level and the structural organization, evolutionary dynamics, and distribution within the genome. ERVs were classified into 28 groups (Cap_ERV_1–28) in terms of sequence identity, structural organization, and the phylogenetic analysis of the sequences. Differential evolutionary dynamics across these groups were recorded. Most of the ERVs had decayed, and only Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 showed signs of activity in recent years, with a few copies encoding long peptides with intact gag, pro, pol, IN, and env domains. We also investigated the intersection of ERV insertions with lncRNA and protein-coding genes, and we found that the majority of lncRNA and protein-coding genes contain ERV-derived sequences in both sheep and goats. These findings describe ERVs’ evolution in the Caprinae subfamily and provide a better understanding of genomic evolution in mammals. Furthermore, putatively active ERVs may allow future genetic marker development, which could have considerable potential for applications in quantitative trait locus (QTL) mapping and molecular breeding in the Caprinae, and so are worthy of additional evaluation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v16030398/s1. All data generated or analyzed during this study are included in this manuscript and its supplementary information files. Supplementary File (S1): Table S1. Data sources of Genomes and ERV reference sequences in the subfamily of Caprinae; Table S2. The information of ERVs (RT-DNA sequences) used for phylogenetic analysis; Table S3. The information of ERVs (Gag MHR) used for phylogenetic analysis; Table S4. The PCR Primers of ERV-RIPs in Caprinae; Table S5. No. of ERVs (contain RT domains) obtained in Cap_ERV groups in Caprinae species; Table S6. Detailed information of “active” full-length ERVs detected in Caprinae genomes; Table S7. The different genes and their domains located in 28 groups of Caprinae ERVs; Table S8. Genotypes identification of selected full-length ERV elements, and Table S9. Detailed information on gene features overlapped with ERVs. Supplementary File (S2): Supplemental_Dataset_S1. Supplementary File (S3): Figure S1. Structure organization of “active” ERVs, and Figure S2. Divergence analysis of four full-length ERVs among the 13 species of the Caprinae subfamily.

Author Contributions

Conceptualization: C.S., J.H. and A.S.M.; Data curation: A.S.M., F.W., Y.Z., C.C. and A.A.S.; Formal analysis: A.S.M., C.C. and A.A.S.; Funding acquisition: C.S. and J.H.; Investigation: A.S.M. and F.W. Methodology: A.S.M., Y.Z. and C.C. Project administration: C.S. and J.H. Supervision: C.S.; Writing—original draft: A.S.M. and C.S.; Writing—review and editing: J.H. and F.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key Research and Development Program of China (2022YFD1300201) and the High-end Talent Support Program of Yangzhou University to Chengyi Song.

Institutional Review Board Statement

All animals and sampling procedures in this experiment were supervised and approved by Alexandria University, Faculty of Agriculture (Al-Shatby), Egypt (2022, No. AU082209203103/ Name of researcher: Dr. Ahmed A. Saleh). All procedures and experimental protocols were in accordance with the Guide for the Care and Use of Agricultural Animals in Research and Teaching, Federation of Animal Science Societies (FASS, 2010) https://www.aaalac.org/about/Ag_Guide_3rd_ed.pdf (accessed on 1 March 2024). Also, the study was carried out in compliance with the ARRIVE guidelines (https://arriveguidelines.org, accessed on 1 March 2024).

Data Availability Statement

All data needed to evaluate the conclusions in this paper are presented either in the main text or the Supplementary Materials.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Gifford, R.; Tristem, M. The evolution, distribution and diversity of endogenous retroviruses. Virus Genes 2003, 26, 291–315. [Google Scholar] [CrossRef]
  2. Rasmussen, H.B. Interactions between Exogenous and Endogenous Retroviruses. J. Biomed. Sci. 1997, 4, 1–8. [Google Scholar] [CrossRef]
  3. Ruprecht, K.; Mayer, J.; Sauter, M.; Roemer, K.; Mueller-Lantzsch, N. Endogenous retroviruses: Endogenous retroviruses and cancer. Cell. Mol. Life Sci. 2008, 65, 3366–3382. [Google Scholar] [CrossRef]
  4. Isbel, L.; Whitelaw, E. Endogenous retroviruses in mammals: An emerging picture of how ERVs modify expression of adjacent genes. Bioessays 2012, 34, 734–738. [Google Scholar] [CrossRef]
  5. Lavialle, C.; Cornelis, G.; Dupressoir, A.; Esnault, C.; Heidmann, O.; Vernochet, C.; Heidmann, T. Paleovirology of ‘syncytins’, retroviral env genes exapted for a role in placentation. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 2013, 368, 20120507. [Google Scholar] [CrossRef]
  6. Jern, P.; Coffin, J.M. Effects of retroviruses on host genome function. Annu. Rev. Genet. 2008, 42, 709–732. [Google Scholar] [CrossRef] [PubMed]
  7. Stoye, J.P. Studies of endogenous retroviruses reveal a continuing evolutionary saga. Nat. Rev. Microbiol. 2012, 10, 395–406. [Google Scholar] [CrossRef]
  8. Hecht, S.J.; Stedman, K.E.; Carlson, J.O.; DeMartini, J.C. Distribution of endogenous type B and type D sheep retrovirus sequences in ungulates and other mammals. Proc. Natl. Acad. Sci. USA 1996, 93, 3297–3302. [Google Scholar] [CrossRef] [PubMed]
  9. Herniou, E.; Martin, J.; Miller, K.; Cook, J.; Wilkinson, M.; Tristem, M. Retroviral diversity and distribution in vertebrates. J. Virol. 1998, 72, 5955–5966. [Google Scholar] [CrossRef] [PubMed]
  10. Blikstad, V.; Benachenhou, F.; Sperber, G.O.; Blomberg, J. Endogenous retroviruses: Evolution of human endogenous retroviral sequences: A conceptual account. Cell. Mol. Life Sci. 2008, 65, 3348–3365. [Google Scholar] [CrossRef] [PubMed]
  11. Stocking, C.; Kozak, C. Endogenous retroviruses: Murine endogenous retroviruses. Cell. Mol. Life Sci. 2008, 65, 3383–3398. [Google Scholar] [CrossRef] [PubMed]
  12. Estes, R.D. The Behavior Guide to African Mammals: Including Hoofed Mammals, Carnivores, Primates; University of California Press: Berkeley, CA, USA, 2012. [Google Scholar]
  13. Huffman, B. The Ultimate Ungulate Page. 2003. Available online: https://www.ultimateungulate.com/Images/imageuse.html (accessed on 7 December 2023).
  14. Robinson, M. Grzimek’s Animal Life. Ref. Rev. 2011, 25, 47–48. [Google Scholar] [CrossRef]
  15. Shackleton, D.M. Wild Sheep and Goats and Their Relatives: Status Survey and Conservation Action Plan for Caprinae; IUCN: Gland, Switzerland, 1997. [Google Scholar]
  16. Alden, C.; Estes, R.D.; Schlitter, D.; McBride, B. National Audubon Society field guide to African wildlife. In National Audubon Society Field Guides; Knopf Doubleday Publishing Group: New York, NY, USA, 1995; Volume 1, pp. 118–142. [Google Scholar]
  17. Danell, K.; Bergström, R.; Duncan, P.; Pastor, J. Large Herbivore Ecology, Ecosystem Dynamics and Conservation; Cambridge University Press: Cambridge, UK, 2006; Volume 11. [Google Scholar]
  18. Feldhamer, G.A.; Merritt, J.F.; Krajewski, C.; Rachlow, J.L.; Stewart, K.M. Mammalogy: Adaptation, Diversity, Ecology; Johns Hopkins University Press: Baltimore, MD, USA, 2020. [Google Scholar]
  19. Matters, N. Convention on International Trade in Endangered Species of Wild Fauna and Flora; CITES: Geneva, Switzerland, 2020. [Google Scholar]
  20. Chessa, B.; Pereira, F.; Arnaud, F.; Amorim, A.; Goyache, F.; Mainland, I.; Kao, R.R.; Pemberton, J.M.; Beraldi, D.; Stear, M.J.; et al. Revealing the history of sheep domestication using retrovirus integrations. Science 2009, 324, 532–536. [Google Scholar] [CrossRef] [PubMed]
  21. Arnaud, F.; Varela, M.; Spencer, T.; Palmarini, M. Endogenous retroviruses: Coevolution of endogenous Betaretroviruses of sheep and their host. Cell. Mol. Life Sci. 2008, 65, 3422–3432. [Google Scholar] [CrossRef] [PubMed]
  22. Klymiuk, N.; Müller, M.; Brem, G.; Aigner, B. Characterization of endogenous retroviruses in sheep. J. Virol. 2003, 77, 11268–11273. [Google Scholar] [CrossRef] [PubMed]
  23. Ellinghaus, D.; Kurtz, S.; Willhoeft, U. LTRharvest, an efficient and flexible software for de novo detection of LTR retrotransposons. BMC Bioinform. 2008, 9, 18. [Google Scholar] [CrossRef]
  24. Gremme, G.; Steinbiss, S.; Kurtz, S. GenomeTools: A comprehensive software library for efficient processing of structured genome annotations. IEEE/ACM Trans. Comput. Biol. Bioinform. 2013, 10, 645–656. [Google Scholar] [CrossRef]
  25. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef]
  26. Potter, S.C.; Luciani, A.; Eddy, S.R.; Park, Y.; Lopez, R.; Finn, R.D. HMMER web server: 2018 update. Nucleic Acids Res. 2018, 46, W200–W204. [Google Scholar] [CrossRef]
  27. El-Gebali, S.; Mistry, J.; Bateman, A.; Eddy, S.R.; Luciani, A.; Potter, S.C.; Qureshi, M.; Richardson, L.J.; Salazar, G.A.; Smart, A. The Pfam protein families database in 2019. Nucleic Acids Res. 2019, 47, D427–D432. [Google Scholar] [CrossRef] [PubMed]
  28. Chen, N. Using Repeat Masker to identify repetitive elements in genomic sequences. Curr. Protoc. Bioinform. 2004, 5, 4–10. [Google Scholar] [CrossRef] [PubMed]
  29. Bao, W.; Kojima, K.K.; Kohany, O. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mob. DNA 2015, 6, 11. [Google Scholar] [CrossRef] [PubMed]
  30. Johnson, L.S.; Eddy, S.R.; Portugaly, E. Hidden Markov model speed heuristic and iterative HMM search procedure. BMC Bioinform. 2010, 11, 431. [Google Scholar] [CrossRef] [PubMed]
  31. Yamada, K.D.; Tomii, K.; Katoh, K. Application of the MAFFT sequence alignment program to large data—Reexamination of the usefulness of chained guide trees. Bioinformatics 2016, 32, 3246–3251. [Google Scholar] [CrossRef] [PubMed]
  32. Nguyen, L.-T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
  33. Kent, W.J. BLAT—the BLAST-like alignment tool. Genome Res. 2002, 12, 656–664. [Google Scholar] [CrossRef] [PubMed]
  34. Du, Z.; Chen, C.; Zheng, Y.; Wang, X.; Song, C. Retroviral Insertion Polymorphism (RIP) of Porcine Endogenous Retroviruses (PERVs) in Pig Genomes. Animals 2024, 14, 621. [Google Scholar] [CrossRef] [PubMed]
  35. Belshaw, R.; Pereira, V.; Katzourakis, A.; Talbot, G.; Pačes, J.; Burt, A.; Tristem, M. Long-term reinfection of the human genome by endogenous retroviruses. Proc. Natl. Acad. Sci. USA 2004, 101, 4894–4899. [Google Scholar] [CrossRef]
  36. Lander, E.S.; Linton, L.M.; Birren, B.; Nusbaum, C.; Zody, M.C.; Baldwin, J.; Devon, K.; Dewar, K.; Doyle, M.; FitzHugh, W. Initial sequencing and analysis of the human genome. Nature 2001, 409, 860–921. [Google Scholar]
  37. McCarthy, E.M.; McDonald, J.F. Long terminal repeat retrotransposons of Mus musculus. Genome Biol. 2004, 5, R14. [Google Scholar] [CrossRef]
  38. Nelson, P.N.; Hooley, P.; Roden, D.; Davari Ejtehadi, H.; Rylance, P.; Warren, P.; Martin, J.; Murray, P.G. Human endogenous retroviruses: Transposable elements with potential? Clin. Exp. Immunol. 2004, 138, 1–9. [Google Scholar] [CrossRef] [PubMed]
  39. Katzourakis, A.; Gifford, R.J.; Tristem, M.; Gilbert, M.T.P.; Pybus, O.G. Macroevolution of complex retroviruses. Science 2009, 325, 1512. [Google Scholar] [CrossRef] [PubMed]
  40. Chen, C.; Wang, W.; Wang, X.; Shen, D.; Wang, S.; Wang, Y.; Gao, B.; Wimmers, K.; Mao, J.; Li, K. Retrotransposons evolution and impact on lncRNA and protein coding genes in pigs. Mob. DNA 2019, 10, 19. [Google Scholar] [CrossRef] [PubMed]
  41. Consortium, M.G.S.; Waterston, R.; Lindblad-Toh, K.; Birney, E.; Rogers, J.; Abril, J.; Agarwal, P.; Agarwala, R.; Ainscough, R.; Alexandersson, M. Initial sequencing and comparative analysis of the mouse genome. Nature 2002, 420, 520–562. [Google Scholar] [CrossRef]
  42. Li, R.; Fan, W.; Tian, G.; Zhu, H.; He, L.; Cai, J.; Huang, Q.; Cai, Q.; Li, B.; Bai, Y. The sequence and de novo assembly of the giant panda genome. Nature 2010, 463, 311–317. [Google Scholar] [CrossRef]
  43. Lindblad-Toh, K.; Wade, C.M.; Mikkelsen, T.S.; Karlsson, E.K.; Jaffe, D.B.; Kamal, M.; Clamp, M.; Chang, J.L.; Kulbokas, E.J.; Zody, M.C. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature 2005, 438, 803–819. [Google Scholar] [CrossRef]
  44. Tanaka, M.; Robinson, B.A.; Chutiraka, K.; Geary, C.D.; Reed, J.C.; Lingappa, J.R. Mutations of conserved residues in the major homology region arrest assembling HIV-1 Gag as a membrane-targeted intermediate containing genomic RNA and cellular proteins. J. Virol. 2016, 90, 1944–1963. [Google Scholar] [CrossRef]
  45. Fang, X.; Mou, Y.; Huang, Z.; Li, Y.; Han, L.; Zhang, Y.; Feng, Y.; Chen, Y.; Jiang, X.; Zhao, W. The sequence and analysis of a Chinese pig genome. Gigascience 2012, 1, 16. [Google Scholar] [CrossRef]
  46. Łopata, K.; Wojdas, E.; Nowak, R.; Łopata, P.; Mazurek, U. Porcine endogenous retrovirus (PERV)–molecular structure and replication strategy in the context of retroviral infection risk of human cells. Front. Microbiol. 2018, 9, 730. [Google Scholar] [CrossRef] [PubMed]
  47. Niu, D.; Wei, H.-J.; Lin, L.; George, H.; Wang, T.; Lee, I.-H.; Zhao, H.-Y.; Wang, Y.; Kan, Y.; Shrock, E. Inactivation of porcine endogenous retrovirus in pigs using CRISPR-Cas9. Science 2017, 357, 1303–1307. [Google Scholar] [CrossRef] [PubMed]
  48. Yang, N.; Zhao, B.; Chen, Y.; D’Alessandro, E.; Chen, C.; Ji, T.; Wu, X.; Song, C. Distinct retrotransposon evolution profile in the genome of rabbit (Oryctolagus cuniculus). Genome Biol. Evol. 2021, 13, evab168. [Google Scholar] [CrossRef] [PubMed]
  49. Anai, Y.; Ochi, H.; Watanabe, S.; Nakagawa, S.; Kawamura, M.; Gojobori, T.; Nishigaki, K. Infectious endogenous retroviruses in cats and emergence of recombinant viruses. J. Virol. 2012, 86, 8634–8644. [Google Scholar] [CrossRef] [PubMed]
  50. Gentles, A.J.; Wakefield, M.J.; Kohany, O.; Gu, W.; Batzer, M.A.; Pollock, D.D.; Jurka, J. Evolutionary dynamics of transposable elements in the short-tailed opossum Monodelphis domestica. Genome Res. 2007, 17, 992–1004. [Google Scholar] [CrossRef] [PubMed]
Viruses 16 00398 g001
Figure 1. Design principles for full-length endogenous retrovirus (ERV) polymorphic primers. (1 & 2) LTR_Primer produces a shorter band ranging between 250 and 750 bp. (1 & 3) Side_Primer produces a band within a size ranging between 1000 and 1500 bp.
Figure 1. Design principles for full-length endogenous retrovirus (ERV) polymorphic primers. (1 & 2) LTR_Primer produces a shorter band ranging between 250 and 750 bp. (1 & 3) Side_Primer produces a band within a size ranging between 1000 and 1500 bp.
Viruses 16 00398 g001
Viruses 16 00398 g002
Figure 2. Classification of endogenous retroviruses (ERVs) in the genomes of Caprinae; (A) a phylogenetic tree was constructed using representative retroviral DNA sequences, aligned with MAFFT and analyzed with the IQ-TREE program. Reference reverse transcriptase (RT) sequences from various mammalian species were included for comparison, including reference RT sequences of enJSRV (DQ838493.1, KC407989.1, KC407987.1, KC407986.1, KC407985.1, and KC407991.1). (B) Another phylogenetic tree was generated using MAFFT alignment and IQ-TREE program, focusing on the reference major homology region (MHR) of gag sequences for comparison. The ERVs were categorized into 28 ERV groups (Cap_ERV_1–28, represented by colored dots). (C) Two classes of ERVs, Class I (gamma retroviruses) and Class II (beta retroviruses), in the Caprinae subfamily are depicted in structural schematics for clarification.
Figure 2. Classification of endogenous retroviruses (ERVs) in the genomes of Caprinae; (A) a phylogenetic tree was constructed using representative retroviral DNA sequences, aligned with MAFFT and analyzed with the IQ-TREE program. Reference reverse transcriptase (RT) sequences from various mammalian species were included for comparison, including reference RT sequences of enJSRV (DQ838493.1, KC407989.1, KC407987.1, KC407986.1, KC407985.1, and KC407991.1). (B) Another phylogenetic tree was generated using MAFFT alignment and IQ-TREE program, focusing on the reference major homology region (MHR) of gag sequences for comparison. The ERVs were categorized into 28 ERV groups (Cap_ERV_1–28, represented by colored dots). (C) Two classes of ERVs, Class I (gamma retroviruses) and Class II (beta retroviruses), in the Caprinae subfamily are depicted in structural schematics for clarification.
Viruses 16 00398 g002
Viruses 16 00398 g003
Figure 3. Landscape of Cap_ERV elements in Caprinae subfamily; (A) the distribution of Cap_ERV groups in the ancestral lineage of the Caprinae subfamily, encompassing six genera and various species. The divergence times are determined using Timetree (http://www.timetree.org, accessed on 22 October 2023) and are represented in million years ago (MYA). (B) The number of ERVs identified in thirteen species within the Caprinae subfamily through the utilization of LTRharvest and RepeatModeler tools. (C) The relative proportions of ERVs, calculated as the size of ERVs relative to the total genome size, expressed as a percentage (ERV size/Total genome size × 100), identified within the genomes of diverse species belonging to the Caprinae subfamily.
Figure 3. Landscape of Cap_ERV elements in Caprinae subfamily; (A) the distribution of Cap_ERV groups in the ancestral lineage of the Caprinae subfamily, encompassing six genera and various species. The divergence times are determined using Timetree (http://www.timetree.org, accessed on 22 October 2023) and are represented in million years ago (MYA). (B) The number of ERVs identified in thirteen species within the Caprinae subfamily through the utilization of LTRharvest and RepeatModeler tools. (C) The relative proportions of ERVs, calculated as the size of ERVs relative to the total genome size, expressed as a percentage (ERV size/Total genome size × 100), identified within the genomes of diverse species belonging to the Caprinae subfamily.
Viruses 16 00398 g003
Viruses 16 00398 g004
Figure 4. Polymorphism identification of selected full-length ERV elements; Cap_ERV_20, 21, 24, and 25 in 6 sheep breeds; Barki (B), Rahmani (R), Rahmani x Barki crossbred (C), Awassi (A), Ossimi (O), and Sulffok (S) using Polymerase Chain Reaction (PCR) approach for blood samples of studied sheep breeds. M1: represents DNA marker (2000 bp) while M2: represents DNA marker (5000 bp). The target product length of the bands is written in yellow numbers.
Figure 4. Polymorphism identification of selected full-length ERV elements; Cap_ERV_20, 21, 24, and 25 in 6 sheep breeds; Barki (B), Rahmani (R), Rahmani x Barki crossbred (C), Awassi (A), Ossimi (O), and Sulffok (S) using Polymerase Chain Reaction (PCR) approach for blood samples of studied sheep breeds. M1: represents DNA marker (2000 bp) while M2: represents DNA marker (5000 bp). The target product length of the bands is written in yellow numbers.
Viruses 16 00398 g004
Viruses 16 00398 g005
Figure 5. Contribution of ERV elements to protein coding and lncRNA genes in both sheep and goats. (A) The proportion of genes intersecting with ERVs in each protein-coding gene (PC gene), and long noncoding gene (lncRNA gene). (B) The proportion of ERVs intersecting with PC and lncRNA genes.
Figure 5. Contribution of ERV elements to protein coding and lncRNA genes in both sheep and goats. (A) The proportion of genes intersecting with ERVs in each protein-coding gene (PC gene), and long noncoding gene (lncRNA gene). (B) The proportion of ERVs intersecting with PC and lncRNA genes.
Viruses 16 00398 g005
Table 1. Detailed information of obtained Cap_ERV elements in Caprinae genomes.
Table 1. Detailed information of obtained Cap_ERV elements in Caprinae genomes.
GroupsRepbase NameNo of ERVsLength
Representative
Sequence (bp)		LTR (bp)gag (aa)pro (aa)pol (aa)IN (aa)env (aa)
Class I (Gamma retroviruses)
Cap_ERV_1OviAri-5.324_LTR14258749260301-891-156
Cap_ERV_2-11238668440364-1021-616
Cap_ERV_3-28728530292443-932-471
Cap_ERV_4OviAri-5.2557_int17598123264559-1176-629
Cap_ERV_5-16558294388341-1010-385
Cap_ERV_6-9467492257137-45157-
Cap_ERV_7-17927316274275-514-290
Cap_ERV_8-5257049388488-633-302
Cap_ERV_9-9706894436526-569-301
Cap_ERV_10OviAri-1.272_LTR24697428381126-333-98
Cap_ERV_11-62197316149179-286-103
Cap_ERV_12OviAri-6.205631907762369102-45237-
Cap_ERV_13OviAri-1.30667508198442565-1168-468
Cap_ERV_14OviAri-3.284_LTR7416889412241-52357321
Cap_ERV_15-7576754276553-376-539
Cap_ERV_16-9247459380348-360-145
Cap_ERV_17-14535920156171-306-111
Cap_ERV_18-39866687183168-349-113
Cap_ERV_19-72779580374257-450-103
Class II (Beta retroviruses)
Cap_ERV_20-1985774530353223570244159
Cap_ERV_21-1470729758352723657244338
Cap_ERV_22-274884425062182342--
Cap_ERV_23-547276511585410558--
Cap_ERV_24-1403787041061628982646618
Cap_ERV_25-472769143752419939343160
Cap_ERV_26-43964351438821750843-
Cap_ERV_27-329484392654439466142-
Cap_ERV_28-25514421330125127401--
Note: ERVs: Endogenous retroviruses, LTR: Long terminal repeats, gag: Group-specific antigen, pro: Protease, pol: Polymerase, IN: Integrase, env: Envelope protein, bp: base pair, and aa: amino acid.
Table 2. Number of recent and potentially active full-length ERVs detected in Caprinae genomes.
Table 2. Number of recent and potentially active full-length ERVs detected in Caprinae genomes.
SpeciesRecent and Potentially Active Full-Length Copies of ERVs
Cap_ERV_20Cap_ERV_21Cap_ERV_24Cap_ERV_25
SheepOvis aries60526541
Ovis ammon15182117
Ovis canadensis34313321
Ovis nivicola33262819
Ovis orientalis40323722
GoatsCapra hircus493411627
Capra aegagrus34172718
Capra ibex1814205
Capra sibirica28213417
RelativesHemitragus hylocrius28183615
Budorcas taxicolor18211412
Oreamnos americanus25181611
Ovibos moschatus36293229
Total418331479254
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Moawad, A.S.; Wang, F.; Zheng, Y.; Chen, C.; Saleh, A.A.; Hou, J.; Song, C. Evolution of Endogenous Retroviruses in the Subfamily of Caprinae. Viruses 2024, 16, 398. https://doi.org/10.3390/v16030398

AMA Style

Moawad AS, Wang F, Zheng Y, Chen C, Saleh AA, Hou J, Song C. Evolution of Endogenous Retroviruses in the Subfamily of Caprinae. Viruses. 2024; 16(3):398. https://doi.org/10.3390/v16030398

Chicago/Turabian Style

Moawad, Ali Shoaib, Fengxu Wang, Yao Zheng, Cai Chen, Ahmed A. Saleh, Jian Hou, and Chengyi Song. 2024. "Evolution of Endogenous Retroviruses in the Subfamily of Caprinae" Viruses 16, no. 3: 398. https://doi.org/10.3390/v16030398

APA Style

Moawad, A. S., Wang, F., Zheng, Y., Chen, C., Saleh, A. A., Hou, J., & Song, C. (2024). Evolution of Endogenous Retroviruses in the Subfamily of Caprinae. Viruses, 16(3), 398. https://doi.org/10.3390/v16030398

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop