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

Comprehensive Identification and Characterization of HML-9 Group in Chimpanzee Genome

Endogenous retroviruses (ERVs) are related to long terminal repeat (LTR) retrotransposons, comprising gene sequences of exogenous retroviruses integrated into the host genome and inherited according to Mendelian law. They are considered to have contributed greatly to the evolution of host genome structure and function. We previously characterized HERV-K HML-9 in the human genome. However, the biological function of this type of element in the genome of the chimpanzee, which is the closest living relative of humans, largely remains elusive. Therefore, the current study aims to characterize HML-9 in the chimpanzee genome and to compare the results with those in the human genome. Firstly, we report the distribution and genetic structural characterization of the 26 proviral elements and 38 solo LTR elements of HML-9 in the chimpanzee genome. The results showed that the distribution of these elements displayed a non-random integration pattern, and only six elements maintained a relatively complete structure. Then, we analyze their phylogeny and reveal that the identified elements all cluster together with HML-9 references and with those identified in the human genome. The HML-9 integration time was estimated based on the 2-LTR approach, and the results showed that HML-9 elements were integrated into the chimpanzee genome between 14 and 36 million years ago and into the human genome between 18 and 49 mya. In addition, conserved motifs, cis-regulatory regions, and enriched PBS sequence features in the chimpanzee genome were predicted based on bioinformatics. The results show that pathways significantly enriched for ERV LTR-regulated genes found in the chimpanzee genome are closely associated with disease development, including neurological and neurodevelopmental psychiatric disorders. In summary, the identification, characterization, and genomics of HML-9 presented here not only contribute to our understanding of the role of ERVs in primate evolution but also to our understanding of their biofunctional significance.

Investigating mobile element variations by statistical genetics

The integration of structural variations (SVs) in statistical genetics provides an opportunity to understand the genetic factors influencing complex human traits and disease. Recent advances in long-read technology and variant calling methods for short reads have improved the accurate discovery and genotyping of SVs, enabling their use in expression quantitative trait loci (eQTL) analysis and genome-wide association studies (GWAS). Mobile elements are DNA sequences that insert themselves into various genome locations. Insertional polymorphisms of mobile elements between humans, called mobile element variations (MEVs), contribute to approximately 25% of human SVs. We recently developed a variant caller that can accurately identify and genotype MEVs from biobank-scale short-read whole-genome sequencing (WGS) datasets and integrate them into statistical genetics. The use of MEVs in eQTL analysis and GWAS has a minimal impact on the discovery of genome loci associated with gene expression and disease; most disease-associated haplotypes can be identified by single nucleotide variations (SNVs). On the other hand, it helps make hypotheses about causal variants or effector variants. Focusing on MEVs, we identified multiple MEVs that contribute to differential gene expression and one of them is a potential cause of skin disease, emphasizing the importance of the integration of MEVs in medical genetics. Here, I will provide an overview of MEVs, MEV calling from WGS, and the integration of MEVs in statistical genetics. Finally, I will discuss the unanswered questions about MEVs, such as rare variants.

Comprehensive characterization of ERV-K (HML-8) in the chimpanzee genome revealed less genomic activity than humans

Endogenous retroviruses (ERVs) originate from ancestral germline infections caused by exogenous retroviruses. Throughout evolution, they have become fixed within the genome of the animals into which they were integrated. As ERV elements coevolve with the host, they are normally epigenetically silenced and can become upregulated in a series of physiological and pathological processes. Generally, a detailed ERV profile in the host genome is critical for understanding the evolutionary history and functional performance of the host genome. We previously characterized and cataloged all the ERV-K subtype HML-8 loci in the human genome; however, this has not been done for the chimpanzee, the nearest living relative of humans. In this study, we aimed to catalog and characterize the integration of HML-8 in the chimpanzee genome and compare it with the integration of HML-8 in the human genome. We analyzed the integration of HML-8 and found that HML-8 pervasively invaded the chimpanzee genome. A total of 76 proviral elements were characterized on 23/24 chromosomes, including detailed elements distribution, structure, phylogeny, integration time, and their potential to regulate adjacent genes. The incomplete structure of HML-8 proviral LTRs will undoubtedly affect their activity. Moreover, the results indicated that HML-8 integration occurred before the divergence between humans and chimpanzees. Furthermore, chimpanzees include more HML-8 proviral elements (76 vs. 40) and fewer solo long terminal repeats (LTR) (0 vs. 5) than humans. These results suggested that chimpanzee genome activity is less than the human genome and that humans may have a better ability to shape and screen integrated proviral elements. Our work is informative in both an evolutionary and a functional context for ERVs.

Endogenous retroviruses in the origins and treatment of cancer

Endogenous retroviruses (ERVs) are emerging as promising therapeutic targets in cancer. As remnants of ancient retroviral infections, ERV-derived regulatory elements coordinate expression from gene networks, including those underpinning embryogenesis and immune cell function. ERV activation can promote an interferon response, a phenomenon termed viral mimicry. Although ERV expression is associated with cancer, and provisionally with autoimmune and neurodegenerative diseases, ERV-mediated inflammation is being explored as a way to sensitize tumors to immunotherapy. Here we review ERV co-option in development and innate immunity, the aberrant contribution of ERVs to tumorigenesis, and the wider biomedical potential of therapies directed at ERVs.

Virus-derived variation in diverse human genomes

Acquisition of genetic material from viruses by their hosts can generate inter-host structural genome variation. We developed computational tools enabling us to study virus-derived structural variants (SVs) in population-scale whole genome sequencing (WGS) datasets and applied them to 3,332 humans. Although SVs had already been cataloged in these subjects, we found previously-overlooked virus-derived SVs. We detected non-germline SVs derived from squirrel monkey retrovirus (SMRV), human immunodeficiency virus 1 (HIV-1), and human T lymphotropic virus (HTLV-1); these variants are attributable to infection of the sequenced lymphoblastoid cell lines (LCLs) or their progenitor cells and may impact gene expression results and the biosafety of experiments using these cells. In addition, we detected new heritable SVs derived from human herpesvirus 6 (HHV-6) and human endogenous retrovirus-K (HERV-K). We report the first solo-direct repeat (DR) HHV-6 likely to reflect DR rearrangement of a known full-length endogenous HHV-6. We used linkage disequilibrium between single nucleotide variants (SNVs) and variants in reads that align to HERV-K, which often cannot be mapped uniquely using conventional short-read sequencing analysis methods, to locate previously-unknown polymorphic HERV-K loci. Some of these loci are tightly linked to trait-associated SNVs, some are in complex genome regions inaccessible by prior methods, and some contain novel HERV-K haplotypes likely derived from gene conversion from an unknown source or introgression. These tools and results broaden our perspective on the coevolution between viruses and humans, including ongoing virus-to-human gene transfer contributing to genetic variation between humans.

Contemporary Distribution, Estimated Age, and Prehistoric Migrations of Old World Monkey Retroviruses

Old World monkeys (OWM), simians inhabiting Africa and Asia, are currently affected by at least four infectious retroviruses, namely, simian foamy virus (SFV), simian immunodeficiency virus (SIV), simian T-lymphotropic virus (STLV), and simian type D retrovirus (SRV). OWM also show chromosomal evidence of having been infected in the past with four more retroviral species, baboon endogenous virus (BaEV), Papio cynocephalus endogenous virus (PcEV), simian endogenous retrovirus (SERV), and Rhesus endogenous retrovirus-K (RhERV-K/SERV-K1). For some of the viruses, transmission to other primates still occurs, resulting, for instance, in the HIV pandemic. Retroviruses are intimately connected with their host as they are normally spread by close contact. In this review, an attempt to reconstruct the distribution and history of OWM retroviruses will be made. A literature overview of the species infected by any of the eight retroviruses as well as an age estimation of the pathogens will be given. In addition, primate genomes from databases have been re-analyzed for the presence of endogenous retrovirus integrations. Results suggest that some of the oldest retroviruses, SERV and PcEV, have travelled with their hosts to Asia during the Miocene, when a higher global temperature allowed simian expansions. In contrast, younger viruses, such as SIV and SRV, probably due to the lack of a primate continuum between the continents in later times, have been restricted to Africa and Asia, respectively.

Human Endogenous Retrovirus K Rec forms a Regulatory Loop with MITF that Opposes the Progression of Melanoma to an Invasive Stage

The HML2 subfamily of HERV-K (henceforth HERV-K) represents the most recently endogenized retrovirus in the human genome. While the products of certain HERV-K genomic copies are expressed in normal tissues, they are upregulated in several pathological conditions, including various tumors. It remains unclear whether HERV-K(HML2)-encoded products overexpressed in cancer contribute to disease progression or are merely by-products of tumorigenesis. Here, we focus on the regulatory activities of the Long Terminal Repeats (LTR5_Hs) of HERV-K and the potential role of the HERV-K-encoded Rec in melanoma. Our regulatory genomics analysis of LTR5_Hs loci indicates that Melanocyte Inducing Transcription Factor (MITF) (also known as binds to a canonical E-box motif (CA(C/T)GTG) within these elements in proliferative type of melanoma, and that depletion of MITF results in reduced HERV-K expression. In turn, experimentally depleting Rec in a proliferative melanoma cell line leads to lower mRNA levels of MITF and its predicted target genes. Furthermore, Rec knockdown leads to an upregulation of epithelial-to-mesenchymal associated genes and an enhanced invasion phenotype of proliferative melanoma cells. Together these results suggest the existence of a regulatory loop between MITF and Rec that may modulate the transition from proliferative to invasive stages of melanoma. Because HERV-K(HML2) elements are restricted to hominoid primates, these findings might explain certain species-specific features of melanoma progression and point to some limitations of animal models in melanoma studies.

How Retroviruses and Retrotransposons in Our Genome May Contribute to Autoimmunity in Rheumatological Conditions

More than 200 human disorders include various manifestations of autoimmunity. The molecular events that lead to these diseases are still incompletely understood and their causes remain largely unknown. Numerous potential triggers of autoimmunity have been proposed over the years, but very few of them have been conclusively confirmed or firmly refuted. Viruses have topped the lists of suspects for decades, and it seems that many viruses, including those of the Herpesviridae family, indeed can influence disease initiation and/or promote exacerbations by a number of mechanisms that include prolonged anti-viral immunity, immune subverting factors, and mechanisms, and perhaps “molecular mimicry”. However, no specific virus has yet been established as being truly causative. Here, we discuss a different, but perhaps mechanistically related possibility, namely that retrotransposons or retroviruses that infected us in the past and left a lasting copy of themselves in our genome still can provoke an escalating immune response that leads to autoimmune disease. Many of these loci still encode for retroviral proteins that have retained some, or all, of their original functions. Importantly, these endogenous proviruses cannot be eliminated by the immune system the way it can eliminate exogenous viruses. Hence, if not properly controlled, they may drive a frustrated and escalating chronic, or episodic, immune response to the point of a frank autoimmune disorder. Here, we discuss the evidence and the proposed mechanisms, and assess the therapeutic options that emerge from the current understanding of this field.

Differences between human and chimpanzee genomes and their implications in gene expression, protein functions and biochemical properties of the two species

Chimpanzees are the closest living relatives of humans. The divergence between human and chimpanzee ancestors dates to approximately 6,5-7,5 million years ago. Genetic features distinguishing us from chimpanzees and making us humans are still of a great interest. After divergence of their ancestor lineages, human and chimpanzee genomes underwent multiple changes including single nucleotide substitutions, deletions and duplications of DNA fragments of different size, insertion of transposable elements and chromosomal rearrangements. Human-specific single nucleotide alterations constituted 1.23% of human DNA, whereas more extended deletions and insertions cover ~ 3% of our genome. Moreover, much higher proportion is made by differential chromosomal inversions and translocations comprising several megabase-long regions or even whole chromosomes. However, despite of extensive knowledge of structural genomic changes accompanying human evolution we still cannot identify with certainty the causative genes of human identity. Most structural gene-influential changes happened at the level of expression regulation, which in turn provoked larger alterations of interactome gene regulation networks. In this review, we summarized the available information about genetic differences between humans and chimpanzees and their potential functional impacts on differential molecular, anatomical, physiological and cognitive peculiarities of these species.

ERVcaller: identifying polymorphic endogenous retrovirus and other transposable element insertions using whole-genome sequencing data

MOTIVATION

Approximately 8% of the human genome is derived from endogenous retroviruses (ERVs). In recent years, an increasing number of human diseases have been found to be associated with ERVs. However, it remains challenging to accurately detect the full spectrum of polymorphic (unfixed) ERVs using whole-genome sequencing (WGS) data.

RESULTS

We designed a new tool, ERVcaller, to detect and genotype transposable element (TE) insertions, including ERVs, in the human genome. We evaluated ERVcaller using both simulated and real benchmark WGS datasets. Compared to existing tools, ERVcaller consistently obtained both the highest sensitivity and precision for detecting simulated ERV and other TE insertions derived from real polymorphic TE sequences. For the WGS data from the 1000 Genomes Project, ERVcaller detected the largest number of TE insertions per sample based on consensus TE loci. By analyzing the experimentally verified TE insertions, ERVcaller had 94.0% TE detection sensitivity and 96.6% genotyping accuracy. Polymerase chain reaction and Sanger sequencing in a small sample set verified 86.7% of examined insertion statuses and 100% of examined genotypes. In conclusion, ERVcaller is capable of detecting and genotyping TE insertions using WGS data with both high sensitivity and precision. This tool can be applied broadly to other species.

AVAILABILITY AND IMPLEMENTATION

http://www.uvm.edu/genomics/software/ERVcaller.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of the distribution of human endogenous retroviruses K (HML-2) by PCR-based target enrichment sequencing

BACKGROUND

Human endogenous retroviruses (HERVs), suspected to be transposition-defective, may reshape the transcriptional network of the human genome by regulatory elements distributed in their long terminal repeats (LTRs). HERV-K (HML-2), the most preserved group with the least number of accumulated of mutations, has been associated with aberrant gene expression in tumorigenesis and autoimmune diseases. Because of the high sequence similarity between different HERV-Ks, current methods have limitations in providing genome-wide mapping specific for individual HERV-K (HML-2) members, a major barrier in delineating HERV-K (HML-2) function.

RESULTS

In an attempt to obtain detailed distribution information of HERV-K (HML-2), we utilized a PCR-based target enrichment sequencing protocol for HERV-K (HML-2) (PTESHK) loci, which not only maps the presence of reference loci, but also identifies non-reference loci, enabling determination of the genome-wide distribution of HERV-K (HML-2) loci. Here we report on the genomic data obtained from three individuals. We identified a total of 978 loci using this method, including 30 new reference loci and 5 non-reference loci. Among the 3 individuals in our study, 14 polymorphic HERV-K (HML-2) loci were identified, and solo-LTR330 and N6p21.32 were identified as polymorphic for the first time.

CONCLUSIONS

Interestingly, PTESHK provides an approach for the identification of the genome-wide distribution of HERV-K (HML-2) and can be used for the identification of polymorphic loci. Since polymorphic HERV-K (HML-2) integrations are suspected to be related to various diseases, PTESHK can supplement other emerging techniques in accessing polymorphic HERV-K (HML-2) elements in cancer and autoimmune diseases.

Is the HERV-K HML-2 Xq21.33, an endogenous retrovirus mutated by gene conversion of chromosome X in a subset of African populations, associated with human breast cancer?

The human endogenous retroviruses HERV-K HML-2 have been considered a possible cause of human breast cancer (BrC). A HERV-K HML-2 fully intact provirus Xq21.33 was recently identified in some West African people. We used PCR technology to search for the Xq21.33 provirus in DNA from Nigerian women with BrC and controls. to see if Xq21.33 plays any role in predisposing to BrC. This provirus was detected in 27 of 216 (12.5%) women with BrC and in 22 of 219 (10.0%) controls. These results were not statistically significant. The prevalence of provirus in premenopausal control women 44 years or younger [18/157 (11.46%)} vs women with BrC [12/117 (10.26%)] showed no statistical difference. The prevalence of virus in postmenopausal control women > 45 yrs. was 7.4% (4/54) vs 15.31% (15/98) in postmenopausal women with BrC. These changes were not statistically significant at <.05, but the actual p value of <.0.079, suggests that Xq21.33 might play some role in predisposing to BrC in postmenopausal women. Provirus was present in Ghanaian women (6/87), in 1/6 Pygmy populations and in African American men (4/45) and women (6/68), but not in any Caucasian women (0/109). Two BrC cell lines (HCC 70 and DT22) from African American women had Xq21.33. Env regions of the virus which differed by 2-3 SNPs did not alter the protein sequence of the virus. SNP at 5730 and 8529 were seen in all persons with provirus, while 54% had an additional SNP at 7596.Two Nigerian women and 2 Ghanaian women had additional unusual SNPs. Homozygosity was seen in (5/27) BrC and (2/22) control women. The genetic variation and homozygosity patterns suggested that there was gene conversion of this X chromosome associated virus. The suggestive finding in this preliminary data of possible increased prevalence of Xq21.33 provirus in post-menopausal Nigerian women with BrC should be clarified by a more statistically powered study sample to see if postmenopausal African and/or African American women carriers of Xq21.33 might show increased risk of BrC. The implication of finding such a link would be the development of antiretroviral drugs that might aid in preventing BrC in Xq21.33+ women.

From Ancient to Emerging Infections: The Odyssey of Viruses in the Male Genital Tract

The male genital tract (MGT) is the target of a number of viral infections that can have deleterious consequences at the individual, offspring, and population levels. These consequences include infertility, cancers of male organs, transmission to the embryo/fetal development abnormalities, and sexual dissemination of major viral pathogens such as human immunodeficiency virus (HIV) and hepatitis B virus. Lately, two emerging viruses, Zika and Ebola, have additionally revealed that the human MGT can constitute a reservoir for viruses cleared from peripheral circulation by the immune system, leading to their sexual transmission by cured men. This represents a concern for future epidemics and further underlines the need for a better understanding of the interplay between viruses and the MGT. We review here how viruses, from ancient viruses that integrated the germline during evolution through old viruses (e.g., papillomaviruses originating from Neanderthals) and more modern sexually transmitted infections (e.g., simian zoonotic HIV) to emerging viruses (e.g., Ebola and Zika) take advantage of genital tract colonization for horizontal dissemination, viral persistence, vertical transmission, and endogenization. The MGT immune responses to viruses and the impact of these infections are discussed. We summarize the latest data regarding the sources of viruses in semen and the complex role of this body fluid in sexual transmission. Finally, we introduce key animal findings that are relevant for our understanding of viral infection and persistence in the human MGT and suggest future research directions.

Characterising a human endogenous retrovirus(HERV)-derived tumour-associated antigen: enriched RNA-Seq analysis of HERV-K(HML-2) in mantle cell lymphoma cell lines

BACKGROUND

The cell-surface attachment protein (Env) of the HERV-K(HML-2) lineage of endogenous retroviruses is a potentially attractive tumour-associated antigen for anti-cancer immunotherapy. The human genome contains around 100 integrated copies (called proviruses or loci) of the HERV-K(HML-2) virus and we argue that it is important for therapy development to know which and how many of these contribute to protein expression, and how this varies across tissues. We measured relative provirus expression in HERV-K(HML-2), using enriched RNA-Seq analysis with both short- and long-read sequencing, in three Mantle Cell Lymphoma cell lines (JVM2, Granta519 and REC1). We also confirmed expression of the Env protein in two of our cell lines using Western blotting, and analysed provirus expression data from all other relevant published studies.

RESULTS

Firstly, in both our and other reanalysed studies, approximately 10% of the transcripts mapping to HERV-K(HML-2) came from Env-encoding proviruses. Secondly, in one cell line the majority of the protein expression appears to come from one provirus (12q14.1). Thirdly, we find a strong tissue-specific pattern of provirus expression.

CONCLUSIONS

A possible dependency of Env expression on a single provirus, combined with the earlier observation that this provirus is not present in all individuals and a general pattern of tissue-specific expression among proviruses, has serious implications for future HERV-K(HML-2)-targeted immunotherapy. Further research into HERV-K(HML-2) as a possible tumour-associated antigen in blood cancers requires a more targeted, proteome-based, screening protocol that will consider these polymorphisms within HERV-K(HML-2). We include a plan (and necessary alignments) for such work.

Structural variation of centromeric endogenous retroviruses in human populations and their impact on cutaneous T-cell lymphoma, Sézary syndrome, and HIV infection

Background

Human Endogenous Retroviruses type K HML-2 (HK2) are integrated into 117 or more areas of human chromosomal arms while two newly discovered HK2 proviruses, K111 and K222, spread extensively in pericentromeric regions, are the first retroviruses discovered in these areas of our genome.

Methods

We use PCR and sequencing analysis to characterize pericentromeric K111 proviruses in DNA from individuals of diverse ethnicities and patients with different diseases.

Results

We found that the 5′ LTR-gag region of K111 proviruses is missing in certain individuals, creating pericentromeric instability. K111 deletion (−/− K111) is seen in about 15% of Caucasian, Asian, and Middle Eastern populations; it is missing in 2.36% of African individuals, suggesting that the −/− K111 genotype originated out of Africa. As we identified the −/−K111 genotype in Cutaneous T-cell lymphoma (CTCL) cell lines, we studied whether the −/−K111 genotype is associated with CTCL. We found a significant increase in the frequency of detection of the −/−K111 genotype in Caucasian patients with severe CTCL and/or Sézary syndrome (n = 35, 37.14%), compared to healthy controls (n = 160, 15.6%) [p = 0.011]. The −/−K111 genotype was also found to vary in HIV-1 infection. Although Caucasian healthy individuals have a similar frequency of detection of the −/− K111 genotype, Caucasian HIV Long-Term Non-Progressors (LTNPs) and/or elite controllers, have significantly higher detection of the −/−K111 genotype (30.55%; n = 36) than patients who rapidly progress to AIDS (8.5%; n = 47) [p = 0.0097].

Conclusion

Our data indicate that pericentromeric instability is associated with more severe CTCL and/or Sézary syndrome in Caucasians, and appears to allow T-cells to survive lysis by HIV infection. These findings also provide new understanding of human evolution, as the −/−K111 genotype appears to have arisen out of Africa and is distributed unevenly throughout the world, possibly affecting the severity of HIV in different geographic areas.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0505-8) contains supplementary material, which is available to authorized users.

Gorillas have been infected with the HERV-K (HML-2) endogenous retrovirus much more recently than humans and chimpanzees

Human endogenous retrovirus-K (HERV-K) human mouse mammary tumor virus-like 2 (HML-2) is the most recently active endogenous retrovirus group in humans, and the only group with human-specific proviruses. HML-2 expression is associated with cancer and other diseases, but extensive searches have failed to reveal any replication-competent proviruses in humans. However, HML-2 proviruses are found throughout the catarrhine primates, and it is possible that they continue to infect some species today. To investigate this possibility, we searched for gorilla-specific HML-2 elements using both in silico data mining and targeted deep-sequencing approaches. We identified 150 gorilla-specific integrations, including 31 2-LTR proviruses. Many of these proviruses have identical LTRs, and are insertionally polymorphic, consistent with very recent integration. One identified provirus has full-length ORFs for all genes, and thus could potentially be replication-competent. We suggest that gorillas may still harbor infectious HML-2 virus and could serve as a model for understanding retrovirus evolution and pathogenesis in humans.

Variation in proviral content among human genomes mediated by LTR recombination

Background

Human endogenous retroviruses (HERVs) occupy a substantial fraction of the genome and impact cellular function with both beneficial and deleterious consequences. The vast majority of HERV sequences descend from ancient retroviral families no longer capable of infection or genomic propagation. In fact, most are no longer represented by full-length proviruses but by solitary long terminal repeats (solo LTRs) that arose via non-allelic recombination events between the two LTRs of a proviral insertion. Because LTR-LTR recombination events may occur long after proviral insertion but are challenging to detect in resequencing data, we hypothesize that this mechanism is a source of genomic variation in the human population that remains vastly underestimated.

Results

We developed a computational pipeline specifically designed to capture dimorphic proviral/solo HERV allelic variants from short-read genome sequencing data. When applied to 279 individuals sequenced as part of the Simons Genome Diversity Project, the pipeline retrieves most of the dimorphic loci previously reported for the HERV-K(HML2) subfamily as well as dozens of additional candidates, including members of the HERV-H and HERV-W families previously involved in human development and disease. We experimentally validate several of these newly discovered dimorphisms, including the first reported instance of an unfixed HERV-W provirus and an HERV-H locus driving a transcript (ESRG) implicated in the maintenance of embryonic stem cell pluripotency.

Conclusions

Our findings indicate that human proviral content exhibit more extensive interindividual variation than previously recognized, which has important bearings for deciphering the contribution of HERVs to human physiology and disease. Because LTR retroelements and LTR recombination are ubiquitous in eukaryotes, our computational pipeline should facilitate the mapping of this type of genomic variation for a wide range of organisms.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0142-3) contains supplementary material, which is available to authorized users.

Epigenetic Control of Human Endogenous Retrovirus Expression: Focus on Regulation of Long-Terminal Repeats (LTRs)

Transposable elements, including endogenous retroviruses (ERVs), comprise almost 45% of the human genome. This could represent a significant pathogenic burden but it is becoming more evident that many of these elements have a positive contribution to make to normal human physiology. In particular, the contributions of human ERVs (HERVs) to gene regulation and the expression of noncoding RNAs has been revealed with the help of new and emerging genomic technologies. HERVs have the common provirus structure of coding open reading frames (ORFs) flanked by two long-terminal repeats (LTRs). However, over the course of evolution and as a consequence of host defence mechanisms, most of the sequences contain INDELs, mutations or have been reduced to single LTRs by recombination. These INDELs and mutations reduce HERV activity. However, there is a trade-off for the host cells in that HERVs can provide beneficial sources of genetic variation but with this benefit comes the risk of pathogenic activity and spread within the genome. For example, the LTRs are of critical importance as they contain promoter sequences and can regulate not only HERV expression but that of human genes. This is true even when the LTRs are located in intergenic regions or are in antisense orientation to the rest of the gene. Uncontrolled, this promoter activity could disrupt normal gene expression or transcript processing (e.g., splicing). Thus, control of HERVs and particularly their LTRs is essential for the cell to manage these elements and this control is achieved at multiple levels, including epigenetic regulations that permit HERV expression in the germline but silence it in most somatic tissues. We will discuss some of the common epigenetic mechanisms and how they affect HERV expression, providing detailed discussions of HERVs in stem cell, placenta and cancer biology.

Human endogenous retroviruses: friend or foe?

The integration of proviral DNA into host chromosomal DNA as an obligatory step in the replication cycle of retroviruses is a natural event of genetic recombination between virus and host. When integration occurs in cells of the germ line, it results in mendelian inheritance of viral sequences that we call endogenous retroviruses (ERV) and HERV for humans. HERVs and host often establish a symbiotic relationship, especially in the placenta and in pluripotent embryonic stem cells, but HERVs occasionally have deleterious consequences for the host. This special issue of APMIS features the fascinating relationships between HERV and humans in health and disease.

Roles for retrotransposon insertions in human disease

Over evolutionary time, the dynamic nature of a genome is driven, in part, by the activity of transposable elements (TE) such as retrotransposons. On a shorter time scale it has been established that new TE insertions can result in single-gene disease in an individual. In humans, the non-LTR retrotransposon Long INterspersed Element-1 (LINE-1 or L1) is the only active autonomous TE. In addition to mobilizing its own RNA to new genomic locations via a "copy-and-paste" mechanism, LINE-1 is able to retrotranspose other RNAs including Alu, SVA, and occasionally cellular RNAs. To date in humans, 124 LINE-1-mediated insertions which result in genetic diseases have been reported. Disease causing LINE-1 insertions have provided a wealth of insight and the foundation for valuable tools to study these genomic parasites. In this review, we provide an overview of LINE-1 biology followed by highlights from new reports of LINE-1-mediated genetic disease in humans.

Human endogenous retrovirus HERV-K(HML-2) activity in prostate cancer is dominated by a few loci

BACKGROUND

Increased expression of human endogenous retroviruses, especially HERV-K(HML-2) proviruses, has recently been associated with prostate carcinoma progression. In particular, a HML-2 locus in chromosome 22q11.23 (H22q) is upregulated in many cases. We therefore aimed at delineating the extent and repertoire of HML-2 transcription in prostate cancer tissues and cell lines and to define the transcription pattern and biological effects of H22q.

METHODS

Sanger and high throughput amplicon sequencing was used to define the repertoire of expressed HML-2 in a selected set of samples. qRT-PCR was used to quantify expression of selected proviruses in an extended set of prostate cancer tissues. Transcription factor binding sites (TFBS) were compared bioinformatically using the Transfac database. Expression of H22q was further characterized by siRNA-mediated knockdown, 5' RACE mapping of transcriptional start sites (TSS) and identification of splice sites. Functional effects of H22q knockdown were investigated by viability and apoptosis assays.

RESULTS

In addition to H22q, a limited number of other proviruses were found expressed by sequencing. Of these, provirus ERVK-5 and to a lesser degree ERVK-15 were frequently upregulated in prostate cancer. In contrast, expression of ERVK-24, predominant in germ cell tumors, was not detectable in prostatic tissues. While HML-2 LTRs contain binding sites for the androgen receptor and cofactors, no consistent differences in transcription factor binding sites were found between expressed and non-expressed proviruses. The H22q locus contains two 5'-LTRs of which the upstream LTR is predominantly used in prostatic cells, with an imprecise TSS. Splicing of H22q transcripts is complex, generating, among others, a transcript with an Np9-like ORF. Knockdown of H22q did not significantly affect proliferation or apoptosis of prostate cancer cells.

CONCLUSIONS

Our findings further underline that HML-2 expression is commonly highly tissue-specific. In prostate cancer, a limited number of loci become activated, especially H22q and ERVK-5. As expressed and non-expressed proviruses do not differ significantly in TFBS, tissue- and tumor-specific expression may be governed primarily by chromatin context. Overexpression of HML-2 H22q is more likely consequence than cause of prostate cancer progression.

Characterizing novel endogenous retroviruses from genetic variation inferred from short sequence reads

From Illumina sequencing of DNA from brain and liver tissue from the lion, Panthera leo, and tumor samples from the pike-perch, Sander lucioperca, we obtained two assembled sequence contigs with similarity to known retroviruses. Phylogenetic analyses suggest that the pike-perch retrovirus belongs to the epsilonretroviruses, and the lion retrovirus to the gammaretroviruses. To determine if these novel retroviral sequences originate from an endogenous retrovirus or from a recently integrated exogenous retrovirus, we assessed the genetic diversity of the parental sequences from which the short Illumina reads are derived. First, we showed by simulations that we can robustly infer the level of genetic diversity from short sequence reads. Second, we find that the measures of nucleotide diversity inferred from our retroviral sequences significantly exceed the level observed from Human Immunodeficiency Virus infections, prompting us to conclude that the novel retroviruses are both of endogenous origin. Through further simulations, we rule out the possibility that the observed elevated levels of nucleotide diversity are the result of co-infection with two closely related exogenous retroviruses.