Endogenous retroviruses: still active after all these years?

Functional Bidirectionality of ERV-Derived Long Non-Coding RNAs in Humans

Human endogenous retroviruses (HERVs) are widely recognized as the result of exogenous retroviruses infecting the ancestral germline, stabilizing integration and vertical transmission during human genetic evolution. To date, endogenous retroviruses (ERVs) appear to have been selected for human physiological functions with the loss of retrotransposable capabilities. ERV elements were previously regarded as junk DNA for a long time. Since then, the aberrant activation and expression of ERVs have been observed in the development of many kinds of human diseases, and their role has been explored in a variety of human disorders such as cancer. The results show that specific ERV elements play respective crucial roles. Among them, long non-coding RNAs (lncRNAs) transcribed from specific long-terminal repeat regions of ERVs are often key factors. lncRNAs are over 200 nucleotides in size and typically bind to DNA, RNA, and proteins to perform biological functions. Dysregulated lncRNAs have been implicated in a variety of diseases. In particular, studies have shown that the aberrant expression of some ERV-derived lncRNAs has a tumor-suppressive or oncogenic effect, displaying significant functional bidirectionality. Therefore, theses lncRNAs have a promising future as novel biomarkers and therapeutic targets to explore the concise relationship between ERVs and cancers. In this review, we first summarize the role of ERV-derived lncRNAs in physiological regulation, mainly including immunomodulation, the maintenance of pluripotency, and erythropoiesis. In addition, pathological regulation examples of their aberrant activation and expression leading to carcinogenesis are highlighted, and specific mechanisms of occurrence are discussed.

Different Host-Endogenous Retrovirus Relationships between Mammals and Birds Reflected in Genome-Wide Evolutionary Interaction Patterns

Mammals and birds differ largely in their average endogenous retrovirus loads, namely the proportion of endogenous retrovirus in the genome. The host-endogenous retrovirus relationships, including conflict and co-option, have been hypothesized among the causes of this difference. However, there has not been studies about the genomic evolutionary signal of constant host-endogenous retrovirus interactions in a long-term scale and how such interactions could lead to the endogenous retrovirus load difference. Through a phylogeny-controlled correlation analysis on ∼5,000 genes between the dN/dS ratio of each gene and the load of endogenous retrovirus in 12 mammals and 21 birds, separately, we detected genes that may have evolved in association with endogenous retrovirus loads. Birds have a higher proportion of genes with strong correlation between dN/dS and the endogenous retrovirus load than mammals. Strong evidence of association is found between the dN/dS of the coding gene for leucine-rich repeat-containing protein 23 and endogenous retrovirus load in birds. Gene set enrichment analysis shows that gene silencing rather than immunity and DNA recombination may have a larger contribution to the association between dN/dS and the endogenous retrovirus load for both mammals and birds. The above results together showing different evolutionary patterns between bird and mammal genes can partially explain the apparently lower endogenous retrovirus loads of birds, while gene silencing may be a universal mechanism that plays a remarkable role in the evolutionary interaction between the host and endogenous retrovirus. In summary, our study presents signals that the host genes might have driven or responded to endogenous retrovirus load changes in long-term evolution.

Roles of Human Endogenous Retroviruses and Endogenous Virus-Like Elements in Cancer Development and Innate Immunity

Human endogenous retroviruses (HERVs) are remnants of ancient retroviral infections in the host genome. Although mutations and silencing mechanisms impair their original role in viral replication, HERVs are believed to play roles in various biological processes. Long interspersed nuclear elements (LINEs) are non-LTR retrotransposons that have a lifecycle resembling that of retroviruses. Although LINE expression is typically silenced in somatic cells, it also contributes to various biological processes. The aberrant expression of HERVs and LINEs is closely associated with the development of cancer and/or immunological diseases, suggesting that they are integrated into various pathways related to the diseases. HERVs/LINEs control gene expression depending on the context as promoter/enhancer elements. Some RNAs and proteins derived from HERVs/LINEs have oncogenic potential, whereas others stimulate innate immunity. Non-retroviral endogenous viral elements (nrEVEs) are a novel type of virus-like element in the genome. nrEVEs may also be involved in host immunity. This article provides a current understanding of how these elements impact cellular physiology in cancer development and innate immunity, and provides perspectives for future studies.

Comparative genomics and genome-wide SNPs of endangered Eld's deer provide breeder selection for inbreeding avoidance

Eld's deer, a conserved wildlife species of Thailand, is facing inbreeding depression, particularly in the captive Siamese Eld's deer (SED) subspecies. In this study, we constructed genomes of a male SED and a male Burmese Eld's deer (BED), and used genome-wide single nucleotide polymorphisms to evaluate the genetic purity and the inbreeding status of 35 SED and 49 BED with limited pedigree information. The results show that these subspecies diverged approximately 1.26 million years ago. All SED were found to be purebred. A low proportion of admixed SED genetic material was observed in some BED individuals. Six potential breeders from male SED with no genetic relation to any female SED and three purebred male BED with no relation to more than 10 purebred female BED were identified. This study provides valuable insights about Eld's deer populations and appropriate breeder selection in efforts to repopulate this endangered species while avoiding inbreeding.

Identification and molecular characterization of simian endogenous retrovirus in Macaca fascicularis and Macaca nemestrina from captive breeding facilities in Bogor, Indonesia

Background and Aim

Endogenous retroviruses (ERVs) found in all vertebrates, including non-human primates (NHPs), are known to be genetically inherited. Thus, recent studies have explored ERVs for human immunodeficiency virus vaccine development using human ERV (HERV) due to the hypervariability of exogenous retroviruses which cause conventional vaccines to be ineffective. HERV was also found to be able to induce an immune response in cancer patients. This study aimed to identify and molecularly characterize ERVs from Indonesian NHPs: Macaca fascicularis and Macaca nemestrina. Then, we described the phylogenetic relationship of these isolates with those of the simian ERVs (SERVs) characterized in other species and countries.

Materials and Methods

First, 5 mL of whole blood samples was taken from 131 long-tailed macaques and 58 pig-tailed macaques in captive breeding facilities at Bogor, Indonesia, for DNA extraction. Next, the DNA samples were screened using the SYBR Green real-time polymerase chain reaction (PCR) technique with specific primers for env (simian retroviruses [SRV]1-5 7585U19 and SRV1-5 7695L21). Positive SERV results were those with cycle threshold (CT) values < 24 (CT < 24) and melting temperature (TM) ranges of 80°C-82°C. Then, whole-genome nucleotide sequences from two pig-tailed macaques samples detected as positive SERV were generated using a nucleic acid sequencing technique which utilized the walking primer method. Subsequently, the sequences were analyzed using bioinformatics programs, such as 4Peaks, Clustal Omega, and BLAST (NCBI). Subsequently, a phylogenetic tree was constructed using the neighbor-joining method in MEGA X.

Results

SYBR Green real-time PCR amplification results indicated that SERV (Mn B1 and Mn B140910)-positive samples had CT values of 22.37-22.54 and TM of 82°C. Moreover, whole-genome sequences resulted in 7991 nucleotide sequences, comprising long terminal repeat, gag, pro, pol, and env genes identical between the sequenced samples. Furthermore, the phylogenetic tree results indicated that both samples from M. nemestrina had 99%-100% nucleotide identities to the Mn 92227 sample identified at the National Primate Center University of Washington (NaPRC UW) which was imported from Indonesia in 1998, confirmed as a novel SERV strain. The phylogenetic tree results also indicated that although SERV whole-genome nucleotide and env amino acid sequences were clustered with SRV-2 (identity values of 82% and 79%, respectively), they had a 99%-100% nucleotide identity to Mn 92227. Meanwhile, the gag, pro, and pol amino acids were clustered with SRV-1, SRV-3, SRV-4, SRV-5, SRV-8, and SERV/1997, with 82% and 88% identity values.

Conclusion

Based on the SYBR Green real-time PCR profiles generated, similarities with Mn 92227 were observed. Subsequent phylogenetic analysis confirmed that both samples (Mn B1 and Mn B140919) from pig-tailed macaques in the country of origin were novel SERV strains at NaPRC UW. Therefore, it could be used in biomedical research on ERVs.

Retroviral Elements in Pathophysiology and as Therapeutic Targets for Amyotrophic Lateral Sclerosis

The study of the role of retroviruses in amyotrophic lateral sclerosis (ALS) dates back to the 1960s shortly after transposable elements themselves were first discovered. It was quickly realized that in wild mice both horizontal and vertical transmissions of retroviral elements were key to the development of an ALS-like syndrome leading to the postulate that endogenous retroviruses (ERVs) contribute significantly to the pathogenicity of this disease. Subsequent studies identified retroviral reverse transcriptase activity in brains of individuals with ALS from Guam. However, except for a single study from the former Soviet Union, ALS could not be transmitted to rhesus macaques. The discovery of an ALS-like syndrome in human immunodeficiency virus (HIV) and human T cell leukemia virus infected individuals led to renewed interest in the field and reverse transcriptase activity was found in the blood and cerebrospinal fluid of individuals with sporadic ALS. However, exogenous retroviruses could not be found in individuals with ALS which further reinforced the possibility of involvement of a human ERV (HERV). The first demonstration of the involvement of a HERV was the discovery of the activation of human endogenous retrovirus-K subtype HML-2 in the brains of individuals with ALS. The envelope protein of HML-2 is neurotoxic and transgenic animals expressing the envelope protein develop an ALS-like syndrome. Activation of HML-2 occurs in the context of generalized transposable element activation and is not specific for ALS. Individuals with HIV-associated ALS show a remarkable response to antiretroviral therapy; however, antiretroviral trials in ALS down-regulate HML-2 without ameliorating the disease. This highlights the need for specific drugs to be developed against HML-2 as a novel therapeutic target for ALS. Other approaches might include antisense oligonucleotides, shRNA targeted against the envelope gene or antibodies that can target the extracellular envelope protein. Future clinical trials in ALS should consider combination therapies to control these ERVs.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.

A beginner's guide to manual curation of transposable elements

BACKGROUND

In the study of transposable elements (TEs), the generation of a high confidence set of consensus sequences that represent the diversity of TEs found in a given genome is a key step in the path to investigate these fascinating genomic elements. Many algorithms and pipelines are available to automatically identify putative TE families present in a genome. Despite the availability of these valuable resources, producing a library of high-quality full-length TE consensus sequences largely remains a process of manual curation. This know-how is often passed on from mentor-to-mentee within research groups, making it difficult for those outside the field to access this highly specialised skill.

RESULTS

Our manuscript attempts to fill this gap by providing a set of detailed computer protocols, software recommendations and video tutorials for those aiming to manually curate TEs. Detailed step-by-step protocols, aimed at the complete beginner, are presented in the Supplementary Methods.

CONCLUSIONS

The proposed set of programs and tools presented here will make the process of manual curation achievable and amenable to all researchers and in special to those new to the field of TEs.

Detection and Characterisation of an Endogenous Betaretrovirus in Australian Wild Deer

Endogenous retroviruses (ERVs) are the remnants of past retroviral infections that once invaded the host’s germline and were vertically transmitted. ERV sequences have been reported in mammals, but their distribution and diversity in cervids are unclear. Using next-generation sequencing, we identified a nearly complete genome of an endogenous betaretrovirus in fallow deer (Dama dama). Further genomic analysis showed that this provirus, tentatively named cervid endogenous betaretrovirus 1 (CERV β1), has typical betaretroviral genome features (gag-pro-pol-env) and the betaretrovirus-specific dUTPase domain. In addition, CERV β1 pol sequences were detected by PCR in the six non-native deer species with wild populations in Australia. Phylogenetic analyses demonstrated that CERV β1 sequences from subfamily Cervinae clustered as sister taxa to ERV-like sequences in species of subfamily Muntiacinae. These findings, therefore, suggest that CERV β1 endogenisation occurred after the split of these two subfamilies (between 3.3 and 5 million years ago). Our results provide important insights into the evolution of betaretroviruses in cervids.

The diversity and evolution of retroviruses: perspectives from viral “fossils”

Retroviruses exclusively infect vertebrates, causing a variety of diseases. The replication of retroviruses requires reverse transcription and integration into host genomes. When infecting germline cells, retroviruses become inherited vertically, forming endogenous retroviruses (ERVs). ERVs document past viral infections, providing molecular fossils for studying the evolutionary history of retroviruses. In this review, we summarize the recent advances in understanding the diversity and evolution of retroviruses from the perspectives of viral fossils, and discuss the effects of ERVs on the evolution of host biology.

•

Recent advances in understanding the diversity and evolution of retroviruses.

•

Methods to analyze ERVs.

•

The effects of ERVs on the evolution of host biology.

Sequence features of retrotransposons allow for epigenetic variability

Transposable elements (TEs) are mobile genetic elements that make up a large fraction of mammalian genomes. While select TEs have been co-opted in host genomes to have function, the majority of these elements are epigenetically silenced by DNA methylation in somatic cells. However, some TEs in mice, including the Intracisternal A-particle (IAP) subfamily of retrotransposons, have been shown to display interindividual variation in DNA methylation. Recent work has revealed that IAP sequence differences and strain-specific KRAB zinc finger proteins (KZFPs) may influence the methylation state of these IAPs. However, the mechanisms underlying the establishment and maintenance of interindividual variability in DNA methylation still remain unclear. Here, we report that sequence content and genomic context influence the likelihood that IAPs become variably methylated. IAPs that differ from consensus IAP sequences have altered KZFP recruitment that can lead to decreased KAP1 recruitment when in proximity of constitutively expressed genes. These variably methylated loci have a high CpG density, similar to CpG islands, and can be bound by ZF-CxxC proteins, providing a potential mechanism to maintain this permissive chromatin environment and protect from DNA methylation. These observations indicate that variably methylated IAPs escape silencing through both attenuation of KZFP binding and recognition by ZF-CxxC proteins to maintain a hypomethylated state.

Adeno-Associated Vector-Delivered CRISPR/SaCas9 System Reduces Feline Leukemia Virus Production In Vitro

Feline leukemia virus (FeLV) is a retrovirus of cats worldwide. High viral loads are associated with progressive infection and the death of the host, due to FeLV-associated disease. In contrast, low viral loads, an effective immune response, and a better clinical outcome can be observed in cats with regressive infection. We hypothesize that by lowering viral loads in progressively infected cats, using CRISPR/SaCas9-assisted gene therapy, the cat’s immune system may be permitted to direct the infection towards a regressive outcome. In a step towards this goal, the present study evaluates different adeno-associated vectors (AAVs) for their competence in delivering a gene editing system into feline cells, followed by investigations of the CRISPR/SaCas9 targeting efficiency for different sites within the FeLV provirus. Nine natural AAV serotypes, two AAV hybrid strains, and Anc80L65, an in silico predicted AAV ancestor, were tested for their potential to infect different feline cell lines and feline primary cells. AAV-DJ revealed superior infection efficiency and was thus employed in subsequent transduction experiments. The introduction of double-strand breaks, using the CRISPR/SaCas9 system targeting 12 selected FeLV provirus sites, was confirmed by T7 endonuclease 1 (T7E1), as well as Tracking of Indels by Decomposition (TIDE) analysis. The highest percentage (up to 80%) of nonhomologous end-joining (NHEJ) was found in the highly conserved gag and pol regions. Subsequent transduction experiments, using AAV-DJ, confirmed indel formation and showed a significant reduction in FeLV p27 antigen for some targets. The targeting of the FeLV provirus was efficient when using the CRISPR/SaCas9 approach in vitro. Whether the observed extent of provirus targeting will be sufficient to provide progressively FeLV-infected cats with the means to overcome the infection needs to be further investigated in vivo.

The impact of endogenous Avian Leukosis Viruses (ALVE) on production traits in elite layer lines

Avian Leukosis Virus subgroup E (ALVE) integrations are endogenous retroviral elements found in the chicken genome. The presence of ALVE has been reported to have negative impacts on multiple traits, including egg production and body weight. The recent development of rapid, inexpensive and specific ALVE detection methods has facilitated their characterization in elite commercial egg production lines across multiple generations. The presence of 20 ALVE was examined in 8 elite lines, from 3 different breeds. Seventeen of these ALVE (85%) were informative and found to be segregating in at least one of the lines. To test for an association between specific ALVE inserts and traits, a large genotype by phenotype study was undertaken. Genotypes were obtained for 500 to 1500 males per line, and the phenotypes used were sire-daughter averages. Phenotype data were analyzed by line with a linear model that included the effects of generation, ALVE genotype and their interaction. If genotype effect was significant, the number of ALVE copies was fitted as a regression to estimate additive ALVE gene substitution effect. Significant associations between the presence of specific ALVE inserts and 18 commercially relevant performance and egg quality traits, including egg production, egg weight and albumen height, were observed. When an ALVE was segregating in more than one line, these associations did not always have the same impact (negative, positive or none) in each line. It is hypothesized that the presence of ALVE in the chicken genome may influence production traits by 3 mechanisms: viral protein production may modulate the immune system and impact overall production performance (virus effect); insertional mutagenesis caused by viral integration may cause direct gene alterations or affect gene regulation (gene effect); or the integration site may be within or adjacent to a quantitative trait region which impacts a performance trait (linkage disequilibrium, marker effect).

Endogenization of a Prosimian Retrovirus during Lemur Evolution

Studies of viruses that coevolved with lemurs provide an opportunity to understand the basal traits of primate viruses and provide an evolutionary context for host-virus interactions. Germline integration of endogenous retroviruses (ERVs) are fossil evidence of past infections. Hence, characterization of novel ERVs provides insight into the ancient precursors of extant viruses and the evolutionary history of their hosts. Here, we report the discovery of a novel endogenous retrovirus present in the genome of a lemur, Coquerel's sifaka (Propithecus coquereli). Using next-generation sequencing, we identified and characterized the complete genome sequence of a retrovirus, named prosimian retrovirus 1 (PSRV1). Phylogenetic analyses indicate that PSRV1 is a gamma-type betaretrovirus basal to the other primate betaretroviruses and most closely related to simian retroviruses. Molecular clock analysis of PSRV1 long terminal repeat (LTR) sequences estimated the time of endogenization within 4.56 MYA (± 2.4 MYA), placing it after the divergence of Propithecus species. These results indicate that PSRV1 is an important milestone of lemur evolution during the radiation of the Propithecus genus. These findings may have implications for both human and animal health in that the acquisition of a gamma-type env gene within an endogenized betaretrovirus could facilitate a cross-species jump between vertebrate class hosts.

Increased copy number of syncytin-1 in the trophectoderm is associated with implantation of the blastocyst

Background

A key step in embryo implantation is the adhesion to and invasion of the endometrium by the blastocyst trophectoderm. The envelope proteins of HERV-W and -FRD (human endogenous retrovirus-W and -FRD), syncytin-1 and syncytin-2, are mainly distributed in the placenta, and play important roles in the development of the placenta. The placenta originates from the trophectoderm of the blastocyst. It is unclear whether the envelope proteins of HERV-W and -FRD have an effect on the development of the trophectoderm and whether they have any association with the implantation of the blastocyst.

Methods

The whole-genome amplification products of the human blastocyst trophectoderm were used to measure the copy number of syncytin-1 and syncytin-2 using real time qPCR. In addition, clinical data associated with the outcome of pregnancies was collected, and included age, body mass index (BMI), basic follicle stimulating hormone(bFSH), rate of primary infertility and oligo-astheno-teratospermia, the thickness of the endometrium on the day of endometrial transformation, the levels of estrogen and progestin on the transfer day, the days and the morphological scores of the blastocysts. The expression of mRNA and the copy numbers of syncytin-1 and syncytin-2 in H1 stem cells, and in differentiated H1 cells, induced by BMP4, were measured using real time qPCR.

Results

The relative copy number of syncytin-1 in the pregnant group (median: 424%, quartile: 232%-463%, p < 0.05) was significantly higher than in the non-pregnant group (median: 100%, quartile: 81%-163%). There was a correlation (r _s  = 0.681, p < 0.001) between the copy number of syncytin-1 and blastocyst implantation after embryo transfer. As the stem cells differentiated, the expression of NANOG mRNA decreased, and the expression of caudal type homeobox 2(CDX2) and β-human chorionic gonadotropin (β-hCG) mRNAs increased. Compared to the undifferentiated cells, the relative expression of the syncytin-1 mRNA was 1.63 (quartile: 0.59-6.37, p > 0.05), 3.36 (quartile: 0.85-14.80, p > 0.05), 10.85 (quartile: 3.39-24.46, p < 0.05) and 67.81 (quartile: 54.07-85.48, p < 0.05) on day 1, 3, 5 and 7, respectively, after the differentiation. The relative expression of syncytin-2 was 5.34 (quartile: 4.50-10.30), 7.90 (quartile: 2.46-14.01), 57.44 (quartile: 38.35-103.87) and 344.76 (quartile: 267.72-440.10) on day 1, 3, 5 and 7, respectively, after the differentiation (p < 0.05). The copy number of syncytin-1 increased significantly during differentiation.

Conclusion

Preceding the transfer of frozen embryos, the increased copy number of syncytin-1 in the blastocyst trophectoderm was associated with good outcomes of pregnancies.

HNSCC: Tumour Antigens and Their Targeting by Immunotherapy

Head and neck squamous cell carcinomas (HNSCC) are a heterogeneous group of malignant tumours typically caused by alcohol and tobacco consumption, although an increasing number of HNSCC arise due to persistent infection with high-risk human papilloma virus (HPV). The treatment of HNSCC remains challenging, and the first-line setting is focused on surgery and chemoradiotherapy. A substantial proportion of HNSCC patients die from their disease, especially those with recurrent and metastatic disease. Among factors linked with good outcome, immune cell infiltration appears to have a major role. HPV-driven HNSCC are often T-cell rich, reflecting the presence of HPV antigens that are immunogenic. Tumour-associated antigens that are shared between patients or that are unique to an individual person may also induce varying degrees of immune response; studying these is important for the understanding of the interaction between the host immune system and the cancer. The resulting knowledge is critical for the design of better immunotherapies. Key questions are: Which antigens lead to an adaptive immune response in the tumour? Which of these are exploitable for immunotherapy? Here, we review the current thinking regarding tumour antigens in HNSCC and what has been learned from early phase clinical trials.

Identification and characterisation of endogenous Avian Leukosis Virus subgroup E (ALVE) insertions in chicken whole genome sequencing data

Background

Endogenous retroviruses (ERVs) are the remnants of retroviral infections which can elicit prolonged genomic and immunological stress on their host organism. In chickens, endogenous Avian Leukosis Virus subgroup E (ALVE) expression has been associated with reductions in muscle growth rate and egg production, as well as providing the potential for novel recombinant viruses. However, ALVEs can remain in commercial stock due to their incomplete identification and association with desirable traits, such as ALVE21 and slow feathering. The availability of whole genome sequencing (WGS) data facilitates high-throughput identification and characterisation of these retroviral remnants.

Results

We have developed obsERVer, a new bioinformatic ERV identification pipeline which can identify ALVEs in WGS data without further sequencing. With this pipeline, 20 ALVEs were identified across eight elite layer lines from Hy-Line International, including four novel integrations and characterisation of a fast feathered phenotypic revertant that still contained ALVE21. These bioinformatically detected sites were subsequently validated using new high-throughput KASP assays, which showed that obsERVer was highly precise and exhibited a 0% false discovery rate. A further fifty-seven diverse chicken WGS datasets were analysed for their ALVE content, identifying a total of 322 integration sites, over 80% of which were novel. Like exogenous ALV, ALVEs show site preference for proximity to protein-coding genes, but also exhibit signs of selection against deleterious integrations within genes.

Conclusions

obsERVer is a highly precise and broadly applicable pipeline for identifying retroviral integrations in WGS data. ALVE identification in commercial layers has aided development of high-throughput diagnostic assays which will aid ALVE management, with the aim to eventually eradicate ALVEs from high performance lines. Analysis of non-commercial chicken datasets with obsERVer has revealed broad ALVE diversity and facilitates the study of the biological effects of these ERVs in wild and domesticated populations.

Diversity of endogenous avian leukosis virus subgroup E (ALVE) insertions in indigenous chickens

Background

Avian leukosis virus subgroup E (ALVE) insertions are endogenous retroviruses (ERV) that are restricted to the domestic chicken and its wild progenitor. In commercial chickens, ALVE are known to have a detrimental effect on productivity and provide a source for recombination with exogenous retroviruses. The wider diversity of ALVE in non-commercial chickens and the role of these elements in ERV-derived immunity (EDI) are yet to be investigated.

Results

In total, 974 different ALVE were identified from 407 chickens sampled from village populations in Ethiopia, Iraq, and Nigeria, using the recently developed obsERVer bioinformatics identification pipeline. Eighty-eight percent of all identified ALVE were novel, bringing the known number of ALVE integrations to more than 1300 across all analysed chickens. ALVE content was highly lineage-specific and populations generally exhibited a large diversity of ALVE at low frequencies, which is typical for ERV involved in EDI. A significantly larger number of ALVE was found within or near coding regions than expected by chance, although a relative depletion of ALVE was observed within coding regions, which likely reflects selection against deleterious integrations. These effects were less pronounced than in previous analyses of chickens from commercial lines.

Conclusions

Identification of more than 850 novel ALVE has trebled the known diversity of these retroviral elements. This work provides the basis for future studies to fully quantify the role of ALVE in immunity against exogenous ALV, and development of programmes to improve the productivity and welfare of chickens in developing economies.

The dual lifestyle of genome-integrating virophages in protists

DNA viruses with efficient host genome integration capability were unknown in eukaryotes until recently. The discovery of virophages, satellite-like DNA viruses that depend on lytic giant viruses that infect protists, revealed a genetically diverse group of viruses with high genome mobility. Virophages can act as strong inhibitors of their associated giant viruses, and the resulting beneficial effects on their unicellular hosts resemble a population-based antiviral defense mechanism. By comparing various aspects of genome-integrating virophages, in particular the virophage mavirus, with other mobile genetic elements and parasite-derived defense mechanisms in eukaryotes and prokaryotes, we show that virophages share many features with other host-parasite systems. Yet, the dual lifestyle exhibited by mavirus remains unprecedented among eukaryotic DNA viruses, with potentially far-reaching ecological and evolutionary consequences for the host.

Tracing the history of LINE and SINE extinction in sigmodontine rodents

Background

L1 retrotransposons have co-evolved with their mammalian hosts for the entire history of mammals and currently compose ~ 20% of a mammalian genome. B1 retrotransposons are dependent on L1 for retrotransposition and span the evolutionary history of rodents since their radiation. L1s were found to have lost their activity in a group of South American rodents, the Sigmodontinae, and B1 inactivation preceded the extinction of L1 in the same group. Consequently, a basal group of sigmodontines have active L1s but inactive B1s and a derived clade have both inactive L1s and B1s. It has been suggested that B1s became extinct during a long period of L1 quiescence and that L1s subsequently reemerged in the basal group.

Results

Here we investigate the evolutionary histories of L1 and B1 in the sigmodontine rodents and show that L1 activity continued until after the L1-extinct clade and the basal group diverged. After the split, L1 had a small burst of activity in the former group, followed by extinction. In the basal group, activity was initially low but was followed by a dramatic increase in L1 activity. We found the last wave of B1 retrotransposition was large and probably preceded the split between the two rodent clades.

Conclusions

Given that L1s had been steadily retrotransposing during the time corresponding to B1 extinction and that the burst of B1 activity preceding B1 extinction was large, we conclude that B1 extinction was not a result of L1 quiescence. Rather, the burst of B1 activity may have contributed to L1 extinction both by competition with L1 and by putting strong selective pressure on the host to control retrotransposition.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0164-5) contains supplementary material, which is available to authorized users.

Whole-Genome Analysis of Domestic Chicken Selection Lines Suggests Segregating Variation in ERV Makeups

Retroviruses have invaded vertebrate hosts for millions of years and left an extensive endogenous retrovirus (ERV) record in the host genomes, which provides a remarkable source for an evolutionary perspective on retrovirus-host associations. Here we identified ERV variation across whole-genomes from two chicken lines, derived from a common founder population subjected to 50 years of bi-directional selection on body weight, and a distantly related domestic chicken line as a comparison outgroup. Candidate ERV loci, where at least one of the chicken lines indicated distinct differences, were analyzed for adjacent host genomic landscapes, selective sweeps, and compared by sequence associations to reference assembly ERVs in phylogenetic analyses. Current data does not support selection acting on specific ERV loci in the domestic chicken lines, as determined by presence inside selective sweeps or composition of adjacent host genes. The varying ERV records among the domestic chicken lines associated broadly across the assembly ERV phylogeny, indicating that the observed insertion differences result from pre-existing and segregating ERV loci in the host populations. Thus, data suggest that the observed differences between the host lineages are best explained by substantial standing ERV variation within host populations, and indicates that even truncated, presumably old, ERVs have not yet become fixed in the host population.

Endogenous Retroviruses Function as Gene Expression Regulatory Elements During Mammalian Pre-implantation Embryo Development

Pre-implantation embryo development encompasses several key developmental events, especially the activation of zygotic genome activation (ZGA)-related genes. Endogenous retroviruses (ERVs), which are regarded as “deleterious genomic parasites”, were previously considered to be “junk DNA”. However, it is now known that ERVs, with limited conservatism across species, mediate conserved developmental processes (e.g., ZGA). Transcriptional activation of ERVs occurs during the transition from maternal control to zygotic genome control, signifying ZGA. ERVs are versatile participants in rewiring gene expression networks during epigenetic reprogramming. Particularly, a subtle balance exists between ERV activation and ERV repression in host–virus interplay, which leads to stage-specific ERV expression during pre-implantation embryo development. A large portion of somatic cell nuclear transfer (SCNT) embryos display developmental arrest and ZGA failure during pre-implantation embryo development. Furthermore, because of the close relationship between ERV activation and ZGA, exploring the regulatory mechanism underlying ERV activation may also shed more light on the enigma of SCNT embryo development in model animals.

Identification, Characterization, and Heritability of Murine Metastable Epialleles: Implications for Non-genetic Inheritance

Generally repressed by epigenetic mechanisms, retrotransposons represent around 40% of the murine genome. At the Agouti viable yellow (Avy) locus, an endogenous retrovirus (ERV) of the intracisternal A particle (IAP) class retrotransposed upstream of the agouti coat-color locus, providing an alternative promoter that is variably DNA methylated in genetically identical individuals. This results in variable expressivity of coat color that is inherited transgenerationally. Here, a systematic genome-wide screen identifies multiple C57BL/6J murine IAPs with Avy epigenetic properties. Each exhibits a stable methylation state within an individual but varies between individuals. Only in rare instances do they act as promoters controlling adjacent gene expression. Their methylation state is locus-specific within an individual, and their flanking regions are enriched for CTCF. Variably methylated IAPs are reprogrammed after fertilization and re-established as variable loci in the next generation, indicating reconstruction of metastable epigenetic states and challenging the generalizability of non-genetic inheritance at these regions.

Nomenclature for endogenous retrovirus (ERV) loci

Retroviral integration into germline DNA can result in the formation of a vertically inherited proviral sequence called an endogenous retrovirus (ERV). Over the course of their evolution, vertebrate genomes have accumulated many thousands of ERV loci. These sequences provide useful retrospective information about ancient retroviruses, and have also played an important role in shaping the evolution of vertebrate genomes. There is an immediate need for a unified system of nomenclature for ERV loci, not only to assist genome annotation, but also to facilitate research on ERVs and their impact on genome biology and evolution. In this review, we examine how ERV nomenclatures have developed, and consider the possibilities for the implementation of a systematic approach for naming ERV loci. We propose that such a nomenclature should not only provide unique identifiers for individual loci, but also denote orthologous relationships between ERVs in different species. In addition, we propose that-where possible-mnemonic links to previous, well-established names for ERV loci and groups should be retained. We show how this approach can be applied and integrated into existing taxonomic and nomenclature schemes for retroviruses, ERVs and transposable elements.

LTR12 promoter activation in a broad range of human tumor cells by HDAC inhibition

A considerable proportion of the human genome consists of transposable elements, including the long terminal repeats (LTRs) of endogenous retroviruses. During evolution, such LTRs were occasionally inserted upstream of protein-coding genes, contributing to their regulation. We previously identified the LTR12 from endogenous retrovirus 9 (ERV9) as a regulator of proapoptotic genes such as TP63 or TNFRSF10B. The promoter activity of LTR12 is largely confined to the testes, silenced in testicular carcinoma, but reactivated in testicular cancer cells by broad-range histone deacetylase (HDAC) inhibitors. Here we show that inhibition of HDAC1-3 is sufficient for LTR12 activation. Importantly, HDAC inhibitors induce LTR12 activity not only in testicular cancer cells, but also in cells derived from many additional tumor species. Finally, we characterize the transcription factor NF-Y as a mediator of LTR12 promoter activity and HDAC inhibitor-induced apoptosis, in the context of widespread genomic binding of NF-Y to specific LTR12 sequences. Thus, HDAC inhibitor-driven LTR12 activation represents a generally applicable means to induce proapoptotic genes in human cancer cells.

Acute Myeloid Leukemia With a Rare t(7;14)(q21;q32) and Trisomy 4 With Poor Clinical Outcome: A Case Report

Objectives

Recurrent cytogenetic abnormalities and/or molecular aberrations play an important role in the diagnosis and prognostification of acute myeloid leukemia (AML). We describe a case of a 40 year old woman diagnosed with de novo AML with a novel t(7;14)(q21,q32) and trisomy 4 with poor clinical outcome. Methods: Morphologic, flow cytometry and cytogenetic results of the patient's peripheral blood and bone marrow samples were analyzed.

Results

The diagnostic bone marrow was hypercellular for age (>95%) with increased blasts (62%) that by flow cytometry exhibited myeloid differentiation with a few T/NK lineage markers. Cytogenetics showed a t(7;14)(q21,q32) and trisomy 4. The patient had extremely poor response to two rounds of induction chemotherapy with persistent leukemia following therapy.

Conclusion

To the best of our knowledge, the t(7;14) is a novel cytogenetic abnormality that has not been reported previously in acute myeloid leukemia, and is important to report as it appears to be associated with poor prognosis.

Identification of a novel HERV-K(HML10): comprehensive characterization and comparative analysis in non-human primates provide insights about HML10 proviruses structure and diffusion

Background

About half of the human genome is constituted of transposable elements, including human endogenous retroviruses (HERV). HERV sequences represent the 8% of our genetic material, deriving from exogenous infections occurred millions of years ago in the germ line cells and being inherited by the offspring in a Mendelian fashion. HERV-K elements (classified as HML1–10) are among the most studied HERV groups, especially due to their possible correlation with human diseases. In particular, the HML10 group was reported to be upregulated in persistent HIV-1 infected cells as well as in tumor cells and samples, and proposed to have a role in the control of host genes expression. An individual HERV-K(HML10) member within the major histocompatibility complex C4 gene has even been studied for its possible contribution to type 1 diabetes susceptibility. Following a first characterization of the HML10 group at the genomic level, performed with the innovative software RetroTector, we have characterized in detail the 8 previously identified HML10 sequences present in the human genome, and an additional HML10 partial provirus in chromosome 1p22.2 that is reported here for the first time.

Results

Using a combined approach based on RetroTector software and a traditional Genome Browser Blat search, we identified a novel HERV-K(HML10) sequence in addition to the eight previously reported in the human genome GRCh37/hg19 assembly. We fully characterized the nine HML10 sequences at the genomic level, including their classification in two types based on both structural and phylogenetic characteristics, a detailed analysis of each HML10 nucleotide sequence, the first description of the presence of an Env Rec domain in the type II HML10, the estimated time of integration of individual members and the comparative map of the HML10 proviruses in non-human primates.

Conclusions

We performed an unambiguous and exhaustive analysis of the nine HML10 sequences present in GRCh37/hg19 assembly, useful to increase the knowledge of the group’s contribution to the human genome and laying the foundation for a better understanding of the potential physiological effects and the tentative correlation of these sequences with human pathogenesis.

Electronic supplementary material

The online version of this article (10.1186/s13100-017-0099-7) contains supplementary material, which is available to authorized users.

Chicken (Gallus gallus) endogenous retrovirus generates genomic variations in the chicken genome

Background

Transposable elements (TEs) comprise ~10% of the chicken (Gallus gallus) genome. The content of TEs is much lower than that of mammalian genomes, where TEs comprise around half of the genome. Endogenous retroviruses are responsible for ~1.3% of the chicken genome. Among them is Gallus gallus endogenous retrovirus 10 (GGERV10), one of the youngest endogenous retrovirus families, which emerged in the chicken genome around 3 million years ago.

Results

We identified a total of 593 GGERV10 elements in the chicken reference genome using UCSC genome database and RepeatMasker. While most of the elements were truncated, 49 GGERV10 elements were full-length retaining 5′ and 3′ LTRs. We examined in detail their structural features, chromosomal distribution, genomic environment, and phylogenetic relationships. We compared LTR sequence among five different GGERV10 subfamilies and found sequence variations among the LTRs. Using a traditional PCR assay, we examined a polymorphism rate of the 49 full-length GGERV10 elements in three different chicken populations of the Korean domestic chicken, Leghorn, and Araucana. The result found a breed-specific GGERV10B insertion locus in the Korean domestic chicken, which could be used as a Korean domestic chicken-specific marker.

Conclusions

GGERV10 family is the youngest ERV family and thus might have contributed to recent genomic variations in different chicken populations. The result of this study showed that one of GGERV10 elements integrated into the chicken genome after the divergence of Korean domestic chicken from other closely related chicken populations, suggesting that GGERV10 could be served as a molecular marker for chicken breed identification.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0085-5) contains supplementary material, which is available to authorized users.

Are human endogenous retroviruses triggers of autoimmune diseases? Unveiling associations of three diseases and viral loci

Autoimmune diseases encompass a plethora of conditions in which the immune system attacks its own tissue, identifying them as foreign. Multiple factors are thought to contribute to the development of immune response to self, including differences in genotypes, hormonal milieu, and environmental factors. Viruses including human endogenous retroviruses have long been linked to the occurrence of autoimmunity, but never proven to be causative factors. Endogenous viruses are retroviral sequences embedded in the host germline DNA and transmitted vertically through successive generations in a Mendelian manner. In this study by means of genetic epidemiology, we have searched for the involvement of endogenous retroviruses in three selected autoimmune diseases: multiple sclerosis, type 1 diabetes mellitus, and rheumatoid arthritis. We found that at least one human endogenous retroviral locus was associated with each of the three diseases. Although there was a significant overlap, most loci only occurred in one of the studied disease. Remarkably, within each disease, there was a statistical interaction (synergy) between two loci. Additional synergy between retroviral loci and human lymphocyte antigens is reported for multiple sclerosis. We speculate the possibility that recombinants or mixed viral particles are formed and that the resulting viruses stimulate the innate immune system, thereby initiating the autoimmune response.

Computational Evaluation of the Strict Master and Random Template Models of Endogenous Retrovirus Evolution

Transposable elements (TEs) are DNA sequences that are able to replicate and move within and between host genomes. Their mechanism of replication is also shared with endogenous retroviruses (ERVs), which are also a type of TE that represent an ancient retroviral infection within animal genomes. Two models have been proposed to explain TE proliferation in host genomes: the strict master model (SMM), and the random template (or transposon) model (TM). In SMM only a single copy of a given TE lineage is able to replicate, and all other genomic copies of TEs are derived from that master copy. In TM, any element of a given family is able to replicate in the host genome. In this paper, we simulated ERV phylogenetic trees under variations of SMM and TM. To test whether current phylogenetic programs can recover the simulated ERV phylogenies, DNA sequence alignments were simulated and maximum likelihood trees were reconstructed and compared to the simulated phylogenies. Results indicate that visual inspection of phylogenetic trees alone can be misleading. However, if a set of statistical summaries is calculated, we are able to distinguish between models with high accuracy by using a data mining algorithm that we introduce here. We also demonstrate the use of our data mining algorithm with empirical data for the porcine endogenous retrovirus (PERV), an ERV that is able to replicate in human and pig cells in vitro.

Transposable elements in the mammalian embryo: pioneers surviving through stealth and service

Transposable elements (TEs) are notable drivers of genetic innovation. Over evolutionary time, TE insertions can supply new promoter, enhancer, and insulator elements to protein-coding genes and establish novel, species-specific gene regulatory networks. Conversely, ongoing TE-driven insertional mutagenesis, nonhomologous recombination, and other potentially deleterious processes can cause sporadic disease by disrupting genome integrity or inducing abrupt gene expression changes. Here, we discuss recent evidence suggesting that TEs may contribute regulatory innovation to mammalian embryonic and pluripotent states as a means to ward off complete repression by their host genome.

Unusual RNA plant virus integration in the soybean genome leads to the production of small RNAs

Horizontal gene transfer (HGT) is known to be a major force in genome evolution. The acquisition of genes from viruses by eukaryotic genomes is a well-studied example of HGT, including rare cases of non-retroviral RNA virus integration. The present study describes the integration of cucumber mosaic virus RNA-1 into soybean genome. After an initial metatranscriptomic analysis of small RNAs derived from soybean, the de novo assembly resulted a 3029-nt contig homologous to RNA-1. The integration of this sequence in the soybean genome was confirmed by DNA deep sequencing. The locus where the integration occurred harbors the full RNA-1 sequence followed by the partial sequence of an endogenous mRNA and another sequence of RNA-1 as an inverted repeat and allowing the formation of a hairpin structure. This region recombined into a retrotransposon located inside an exon of a soybean gene. The nucleotide similarity of the integrated sequence compared to other Cucumber mosaic virus sequences indicates that the integration event occurred recently. We described a rare event of non-retroviral RNA virus integration in soybean that leads to the production of a double-stranded RNA in a similar fashion to virus resistance RNAi plants.

Evaluation of Vero-cell-derived simian endogenous retrovirus infection in humans by detection of viral genome in clinicopathological samples and commercialized vaccines and by serology of Japanese general population

BACKGROUND

Vero cells are used in laboratories for the isolation of viruses and the production of vaccines. Recently, the sequence of simian endogenous retrovirus (SERV) was identified in Vero cells (SERVagm-Vero), with homology to exogenously transmitted, pathogenic simian retroviruses (SRVs). Although SERVagm-Vero was shown to be noninfectious to human cells in vitro, SERV infection in humans is controversial. In this study, we evaluated the status of SERV infection in humans by detecting the viral genome in clinicopathological samples and commercialized vaccines, and with a serological survey of the Japanese general population.

METHODS

Real-time polymerase chain reaction (PCR) and reverse transcription-PCR were used to detect the SERVagm-Vero genome. We also examined the seroprevalence of SERV in 1000 individuals in the Japanese population with an enzyme-linked immunoabsorbent assay (ELISA) using mixed SERVagm-Vero gag, pol, and env proteins as antigens.

RESULTS

Real-time PCR failed to detect SERVagm-Vero genomic fragments in 783 human clinicopathological samples, and all 13 human cell lines tested were negative for the SERVagm-Vero genome. Thirteen commercialized vaccines, including five Vero-based vaccines, were also negative for the SERVagm-Vero genome on real-time PCR and reverse transcription-PCR. Eight (0.8%) were seropositive on ELISA, and western blotting showed that all eight sera contained anti-pol antibodies. All SERV-seropositive individuals were born before 1965, suggesting that SERV infection in Japan might not be associated with vaccine, because more than 90% of Japanese children born from 1964 to 2012 have received live poliovirus vaccines containing virus produced in Vero cells since the 1980s.

CONCLUSION

We have confirmed that the vaccines we use today are free of SERVagm-Vero. Moreover, SERV infection in humans is very rare and unlikely to be associated with vaccines in the Japanese general population.

Classification and characterization of human endogenous retroviruses; mosaic forms are common

Background

Human endogenous retroviruses (HERVs) represent the inheritance of ancient germ-line cell infections by exogenous retroviruses and the subsequent transmission of the integrated proviruses to the descendants. ERVs have the same internal structure as exogenous retroviruses. While no replication-competent HERVs have been recognized, some retain up to three of four intact ORFs. HERVs have been classified before, with varying scope and depth, notably in the RepBase/RepeatMasker system. However, existing classifications are bewildering. There is a need for a systematic, unifying and simple classification. We strived for a classification which is traceable to previous classifications and which encompasses HERV variation within a limited number of clades.

Results

The human genome assembly GRCh 37/hg19 was analyzed with RetroTector, which primarily detects relatively complete Class I and II proviruses. A total of 3173 HERV sequences were identified. The structure of and relations between these proviruses was resolved through a multi-step classification procedure that involved a novel type of similarity image analysis (“Simage”) which allowed discrimination of heterogeneous (noncanonical) from homogeneous (canonical) HERVs. Of the 3173 HERVs, 1214 were canonical and segregated into 39 canonical clades (groups), belonging to class I (Gamma- and Epsilon-like), II (Beta-like) and III (Spuma-like). The groups were chosen based on (1) sequence (nucleotide and Pol amino acid), similarity, (2) degree of fit to previously published clades, often from RepBase, and (3) taxonomic markers. The groups fell into 11 supergroups. The 1959 noncanonical HERVs contained 31 additional, less well-defined groups. Simage analysis revealed several types of mosaicism, notably recombination and secondary integration. By comparing flanking sequences, LTRs and completeness of gene structure, we deduced that some noncanonical HERVs proliferated after the recombination event. Groups were further divided into envelope subgroups (altogether 94) based on sequence similarity and characteristic “immunosuppressive domain” motifs. Intra and inter(super)group, as well as intraclass, recombination involving envelope genes (“env snatching”) was a common event. LTR divergence indicated that HERV-K(HML2) and HERVFC had the most recent integrations, HERVL and HUERSP3 the oldest.

Conclusions

A comprehensive HERV classification and characterization approach was undertaken. It should be applicable for classification of all ERVs. Recombination was common among HERV ancestors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0232-y) contains supplementary material, which is available to authorized users.

An Evolutionarily Young Polar Bear (Ursus maritimus) Endogenous Retrovirus Identified from Next Generation Sequence Data

Transcriptome analysis of polar bear (Ursus maritimus) tissues identified sequences with similarity to Porcine Endogenous Retroviruses (PERV). Based on these sequences, four proviral copies and 15 solo long terminal repeats (LTRs) of a newly described endogenous retrovirus were characterized from the polar bear draft genome sequence. Closely related sequences were identified by PCR analysis of brown bear (Ursus arctos) and black bear (Ursus americanus) but were absent in non-Ursinae bear species. The virus was therefore designated UrsusERV. Two distinct groups of LTRs were observed including a recombinant ERV that contained one LTR belonging to each group indicating that genomic invasions by at least two UrsusERV variants have recently occurred. Age estimates based on proviral LTR divergence and conservation of integration sites among ursids suggest the viral group is only a few million years old. The youngest provirus was polar bear specific, had intact open reading frames (ORFs) and could potentially encode functional proteins. Phylogenetic analyses of UrsusERV consensus protein sequences suggest that it is part of a pig, gibbon and koala retrovirus clade. The young age estimates and lineage specificity of the virus suggests UrsusERV is a recent cross species transmission from an unknown reservoir and places the viral group among the youngest of ERVs identified in mammals.

Human Endogenous Retrovirus Type K (HERV-K) Particles Package and Transmit HERV-K-Related Sequences

Human endogenous retroviruses (HERV) make up 8% of the human genome. While the youngest of these retroviruses, HERV-K(HML-2), termed HK2, is able to code for all viral proteins and produce virus-like particles, it is not known if these virus particles package and transmit HK2-related sequences. Here, we analyzed the capacity of HK2 for packaging and transmitting HK2 sequences. We created an HK2 probe, termed Bogota, which can be packaged into HK2 viruses, and transfected it into cells that make HK2 particles. Supernatants of the transfected cells, which contained HK2 viral particles, then were added to target cells, and the transmissibility of the HK2 Bogota reporter was tracked by G418 resistance. Our studies revealed that contemporary HK2 virions produced by some teratocarcinoma and breast cancer cell lines, as well as by peripheral blood lymphocytes from lymphoma patients, can package HK2 Bogota probes, and these viruses transmitted these probes to other cells. After transmission, HK2 Bogota transcripts undergo reverse transcription, a step impaired by antiretroviral agents or by introduction of mutations into the probe sequences required for reverse transcription. HK2 viruses were more efficiently transmitted in the presence of HK2 Rec or HIV-1 Tat and Vif. Transmitted Bogota probes formed episomes but did not integrate into the cellular genome. Resistance to integration might explain the relatively low number of HK2 insertions that were acquired during the last 25 million years of evolution. Whether transient transmission of modern HK2 sequences, which encode two putative oncoproteins, can lead to disease remains to be studied.

IMPORTANCE

Retroviruses invaded the genome of human ancestors over the course of millions of years, yet these viruses generally have been inactivated during evolution, with only remnants of these infectious sequences remaining in the human genome. One of these viruses, termed HK2, still is capable of producing virus particles, although these particles have been regarded as being noninfectious. Using a genetic probe derived from HK2, we have discovered that HK2 viruses produced in modern humans can package HK2 sequences and transmit them to various other cells. Furthermore, the genetic sequences packaged in HK2 undergo reverse transcription. The transmitted probe circularized in the cell and failed to integrate into the cellular genome. These findings suggest that modern HK2 viruses can package viral RNA and transmit it to other cells. Contrary to previous views, we provide evidence of an extracellular viral phase of modern HK2 viruses. We have no evidence of sustained, spreading infection.

Genome-wide amplification of proviral sequences reveals new polymorphic HERV-K(HML-2) proviruses in humans and chimpanzees that are absent from genome assemblies

BACKGROUND

To date, the human population census of proviruses of the Betaretrovirus-like human endogenous retroviral (HERV-K) (HML-2) family has been compiled from a limited number of complete genomes, making it certain that rare polymorphic loci are under-represented and are yet to be described.

RESULTS

Here we describe a suppression PCR-based method called genome-wide amplification of proviral sequences (GAPS) that selectively amplifies DNA fragments containing the termini of HERV-K(HML-2) proviral sequences and their flanking genomic sequences. We analysed the HERV-K(HML-2) proviral content of 101 unrelated humans, 4 common chimpanzees and three centre d'etude du polymorphisme humain (CEPH) pedigrees (44 individuals). The technique isolated HERV-K(HML-2) proviruses that had integrated in the genomes of the great apes throughout their divergence and included evolutionarily young elements still unfixed for presence/absence.

CONCLUSIONS

By examining the HERV-K(HML-2) proviral content of 145 humans we detected a new insertionally polymorphic Type I HERV-K(HML-2) provirus. We also observed provirus versus solo long terminal repeat (LTR) polymorphism within humans at a previously unreported, but ancient, locus. Finally, we report two novel chimpanzee specific proviruses, one of which is dimorphic for a provirus versus solo LTR. Thus GAPS enables the isolation of uncharacterised HERV-K(HML-2) proviral sequences and provides a direct means to assess inter-individual genetic variation associated with HERV-K(HML-2) proviruses.

The effect of life history on retroviral genome invasions

Endogenous retroviruses (ERV), or the remnants of past retroviral infections that are no longer active, are found in the genomes of most vertebrates, typically constituting approximately 10% of the genome. In some vertebrates, particularly in shorter-lived species like rodents, it is not unusual to find active endogenous retroviruses. In longer-lived species, including humans where substantial effort has been invested in searching for active ERVs, it is unusual to find them; to date none have been found in humans. Presumably the chance of detecting an active ERV infection is a function of the length of an ERV epidemic. Intuitively, given that ERVs or signatures of past ERV infections are passed from parents to offspring, we might expect to detect more active ERVs in species with longer generation times, as it should take more years for an infection to run its course in longer than in shorter lived species. This means the observation of more active ERV infections in shorter compared to longer-lived species is paradoxical. We explore this paradox using a modeling approach to investigate factors that influence ERV epidemic length. Our simple epidemiological model may explain why we find evidence of active ERV infections in shorter rather than longer-lived species.

CMV-promoter driven codon-optimized expression alters the assembly type and morphology of a reconstituted HERV-K(HML-2)

The HERV-K(HML-2) family contains the most recently integrated and best preserved endogenized proviral sequences in the human genome. All known elements have nevertheless been subjected to mutations or deletions that render expressed particles non-infectious. Moreover, these post-insertional mutations hamper the analysis of the general biological properties of this ancient virus family. The expression of consensus sequences and sequences of elements with reverted post-insertional mutations has therefore been very instrumental in overcoming this limitation. We investigated the particle morphology of a recently reconstituted HERV-K113 element termed oriHERV-K113 using thin-section electron microscopy (EM) and could demonstrate that strong overexpression by substitution of the 5'LTR for a CMV promoter and partial codon optimization altered the virus assembly type and morphology. This included a conversion from the regular C-type to an A-type morphology with a mass of cytoplasmic immature cores tethered to the cell membrane and the membranes of vesicles. Overexpression permitted the release and maturation of virions but reduced the envelope content. A weaker boost of virus expression by Staufen-1 was not sufficient to induce these morphological alterations.

Chimpanzee-specific endogenous retrovirus generates genomic variations in the chimpanzee genome

Endogenous retroviruses (ERVs), eukaryotic transposable elements, exist as proviruses in vertebrates including primates and contribute to genomic changes during the evolution of their host genomes. Many studies about ERVs have focused on the elements residing in the human genome but only a few studies have focused on the elements which exist in non-human primate genomes. In this study, we identified 256 chimpanzee-specific endogenous retrovirus copies (PtERVs: Pan troglodyte endogenous retroviruses) from the chimpanzee reference genome sequence through comparative genomics. Among the chimpanzee-specific ERV copies, 121 were full-length chimpanzee-specific ERV elements while 110 were chimpanzee-specific solitary LTR copies. In addition, we found eight potential retrotransposition-competent full-length chimpanzee-specific ERV copies containing an intact env gene, and two of them were polymorphic in chimpanzee individuals. Through computational analysis and manual inspection, we found that some of the chimpanzee-specific ERVs have propagated via non-classical PtERV insertion (NCPI), and at least one of the PtERVs may have played a role in creating an alternative transcript of a chimpanzee gene. Based on our findings in this study, we state that the chimpanzee-specific ERV element is one of the sources of chimpanzee genomic variations, some of which might be related to the alternative transcripts in the chimpanzee population.

Hybridization capture reveals evolution and conservation across the entire Koala retrovirus genome

The koala retrovirus (KoRV) is the only retrovirus known to be in the midst of invading the germ line of its host species. Hybridization capture and next generation sequencing were used on modern and museum DNA samples of koala (Phascolarctos cinereus) to examine ca. 130 years of evolution across the full KoRV genome. Overall, the entire proviral genome appeared to be conserved across time in sequence, protein structure and transcriptional binding sites. A total of 138 polymorphisms were detected, of which 72 were found in more than one individual. At every polymorphic site in the museum koalas, one of the character states matched that of modern KoRV. Among non-synonymous polymorphisms, radical substitutions involving large physiochemical differences between amino acids were elevated in env, potentially reflecting anti-viral immune pressure or avoidance of receptor interference. Polymorphisms were not detected within two functional regions believed to affect infectivity. Host sequences flanking proviral integration sites were also captured; with few proviral loci shared among koalas. Recently described variants of KoRV, designated KoRV-B and KoRV-J, were not detected in museum samples, suggesting that these variants may be of recent origin.

Evolutionary dynamics of endogenous Jaagsiekte sheep retroviruses proliferation in the domestic sheep, mouflon and Pyrenean chamois

The oncogenic exogenous Jaagsiekte sheep retrovirus (JSRV), responsible for ovine pulmonary adenocarcinoma, has several endogenous counterparts termed enJSRVs. Although many of these elements have been inactivated over time by the accumulation of deleterious mutations or internal recombination leading to solo long terminal repeat (LTR) formation, several members of enJSRVs have been identified as nearly intact and probably represent recent integration events. To determine the level of enJSRV polymorphism in the sheep population and related species, we have undertaken a study by characterizing enJSRVs copies and independent integration sites in six domestic sheep and two wild species of the sheep lineage. enJSRVs copies were detected by amplifying the env-LTR region by PCR, and for the detection of the insertion sites, we used two approaches: (1) an in silico approach based on the recently published Sheep Reference Genome Assembly (OARv3.0) and (2) an experimental approach based on PCR suppression and inverse PCR techniques. In total, 103 enJSRV sequences were generated across 10 individuals and enJSRV integrations were found on 11 of the 28 sheep chromosomes. These findings suggest that there are still uncharacterized enJSRVs, and that some of the integration sites are variable among the different species, breeds of the same species, subspecies and geographic locations.

A novel recombinant retrovirus in the genomes of modern birds combines features of avian and mammalian retroviruses

Endogenous retroviruses (ERVs) represent ancestral sequences of modern retroviruses or their extinct relatives. The majority of ERVs cluster alongside exogenous retroviruses into two main groups based on phylogenetic analyses of the reverse transcriptase (RT) enzyme. Class I includes gammaretroviruses, and class II includes lentiviruses and alpha-, beta-, and deltaretroviruses. However, analyses of the transmembrane subunit (TM) of the envelope glycoprotein (env) gene result in a different topology for some retroviruses, suggesting recombination events in which heterologous env sequences have been acquired. We previously demonstrated that the TM sequences of five of the six genera of orthoretroviruses can be divided into three types, each of which infects a distinct set of vertebrate classes. Moreover, these classes do not always overlap the host range of the associated RT classes. Thus, recombination resulting in acquisition of a heterologous env gene could in theory facilitate cross-species transmissions across vertebrate classes, for example, from mammals to reptiles. Here we characterized a family of class II avian ERVs, "TgERV-F," that acquired a mammalian gammaretroviral env sequence. Although TgERV-F clusters near a sister clade to alpharetroviruses, its genome also has some features of betaretroviruses. We offer evidence that this unusual recombinant has circulated among several avian orders and may still have infectious members. In addition to documenting the infection of a nongalliform avian species by a mammalian retrovirus, TgERV-F also underscores the importance of env sequences in reconstructing phylogenies and supports a possible role for env swapping in allowing cross-species transmissions across wide taxonomic distances.

IMPORTANCE

Retroviruses can sometimes acquire an envelope gene (env) from a distantly related retrovirus. Since env is a key determinant of host range, such an event affects the host range of the recombinant virus and can lead to the creation of novel retroviral lineages. Retroviruses insert viral DNA into the host DNA during infection, and therefore vertebrate genomes contain a "fossil record" of endogenous retroviral sequences thought to represent past infections of germ cells. We examined endogenous retroviral sequences in avian genomes for evidence of recombination events involving env. Although cross-species transmissions of retroviruses between vertebrate classes (from mammals to birds, for example) are thought to be rare, we here characterized a group of avian retroviruses that acquired an env sequence from a mammalian retrovirus. We offer evidence that this unusual recombinant circulated among songbirds 2 to 4 million years ago and has remained active into the recent past.

Comprehensive analysis of human endogenous retrovirus group HERV-W locus transcription in multiple sclerosis brain lesions by high-throughput amplicon sequencing

Human endogenous retroviruses (HERVs) of the HERV-W group comprise hundreds of loci in the human genome. Deregulated HERV-W expression and HERV-W locus ERVWE1-encoded Syncytin-1 protein have been implicated in the pathogenesis of multiple sclerosis (MS). However, the actual transcription of HERV-W loci in the MS context has not been comprehensively analyzed. We investigated transcription of HERV-W in MS brain lesions and white matter brain tissue from healthy controls by employing next-generation amplicon sequencing of HERV-W env-specific reverse transcriptase (RT) PCR products, thus revealing transcribed HERV-W loci and the relative transcript levels of those loci. We identified more than 100 HERV-W loci that were transcribed in the human brain, with a limited number of loci being predominantly transcribed. Importantly, relative transcript levels of HERV-W loci were very similar between MS and healthy brain tissue samples, refuting deregulated transcription of HERV-W env in MS brain lesions, including the high-level-transcribed ERVWE1 locus encoding Syncytin-1. Quantitative RT-PCR likewise did not reveal differences in MS regarding HERV-W env general transcript or ERVWE1- and ERVWE2-specific transcript levels. However, we obtained evidence for interindividual differences in HERV-W transcript levels. Reporter gene assays indicated promoter activity of many HERV-W long terminal repeats (LTRs), including structurally incomplete LTRs. Our comprehensive analysis of HERV-W transcription in the human brain thus provides important information on the biology of HERV-W in MS lesions and normal human brain, implications for study design, and mechanisms by which HERV-W may (or may not) be involved in MS.