HIV infection reveals widespread expansion of novel centromeric human endogenous retroviruses

Beyond pathogens: the intriguing genetic legacy of endogenous retroviruses in host physiology

The notion that viruses played a crucial role in the evolution of life is not a new concept. However, more recent insights suggest that this perception might be even more expansive, highlighting the ongoing impact of viruses on host evolution. Endogenous retroviruses (ERVs) are considered genomic remnants of ancient viral infections acquired throughout vertebrate evolution. Their exogenous counterparts once infected the host's germline cells, eventually leading to the permanent endogenization of their respective proviruses. The success of ERV colonization is evident so that it constitutes 8% of the human genome. Emerging genomic studies indicate that endogenous retroviruses are not merely remnants of past infections but rather play a corollary role, despite not fully understood, in host genetic regulation. This review presents some evidence supporting the crucial role of endogenous retroviruses in regulating host genetics. We explore the involvement of human ERVs (HERVs) in key physiological processes, from their precise and orchestrated activities during cellular differentiation and pluripotency to their contributions to aging and cellular senescence. Additionally, we discuss the costs associated with hosting a substantial amount of preserved viral genetic material.

Determination of the Copy Number of Porcine Endogenous Retroviruses (PERV) in Auckland Island Pigs Repeatedly Used for Clinical Xenotransplantation and Elimination of PERV-C

Auckland Island pigs represent an inbred population of feral pigs isolated on the sub-Antarctic island for over 100 years. The animals have been maintained under pathogen-free conditions in New Zealand; they are well characterized virologically and have been used as donor sources in first clinical trials of porcine neonatal islet cell transplantation for the treatment of human diabetes patients. The animals do not carry any of the xenotransplantation-relevant viruses, and in the first clinical trials, no porcine viruses, including porcine endogenous retroviruses (PERVs) were transmitted to the human recipients. PERVs pose a special risk in xenotransplantation, since they are part of the pig genome. When the copy number of PERVs in these animals was analyzed using droplet digital PCR and primers binding to a conserved region of the polymerase gene (PERVpol), a copy number typical for Western pigs was found. This confirms previous phylogenetic analyses of microsatellites as well as mitochondrial analyses showing a closer relationship to European pigs than to Chinese pigs. When kidney cells from very young piglets were analyzed, only around 20 PERVpol copies were detected. Using these cells as donors in somatic cell nuclear transfer (SCNT), animals were born showing PERVpol copy numbers between 35 and 56. These data indicate that Auckland Island pigs have a similar copy number in comparison with other Western pig breeds and that the copy number is higher in adult animals compared with cells from young piglets. Most importantly, PERV-C-free animals were selected and the absence of an additional eight porcine viruses was demonstrated.

Kangaroo endogenous retrovirus (KERV) forms megasatellite DNA with a simple repetition pattern in which the provirus structure is retained

The kangaroo endogenous retrovirus (KERV) was previously reported to have undergone a rapid copy number increase in the red-necked wallaby; however, the mode of amplification was left to be clarified. The present study revealed that the long terminal repeat (LTR) (0.6 kb) and internal region (2.0 kb) of a provirus are repeated alternately, forming megasatellite DNA which we named kervRep. This repetition pattern was the same as that observed for walbRep, megasatellite DNA originating from another endogenous retrovirus. Their formation process can be explained using a simple model: pairing slippage followed by homologous recombination. This model features that the initial step is triggered by the presence of two identical sequences within a short distance; the possession of LTRs by endogenous retroviruses fulfills this condition. The discovery of two cases suggests that formation of this type of satellite DNA is one of non-negligible effects of endogenous retroviruses on their host genomes.

Interplay between endogenous and exogenous human retroviruses

In humans, retroviruses thrive more as symbionts than as parasites. Apart from the only two modern exogenous human retroviruses (human T-cell lymphotropic and immunodeficiency viruses; HTLV and HIV, respectively), ~8% of the human genome is occupied by ancient retroviral DNA [human endogenous retroviruses (HERVs)]. Here, we review the recent discoveries about the interactions between the two groups, the impact of infection by exogenous retroviruses on the expression of HERVs, the effect of HERVs on the pathogenicity of HIV and HTLV and on the severity of the diseases caused by them, and the antiviral protection that HERVs can allegedly provide to the host. Tracing the crosstalk between contemporary retroviruses and their endogenized ancestors will provide better understanding of the retroviral world.

The largest HIV-1-infected T cell clones in children on long-term combination antiretroviral therapy contain solo LTRs

Combination antiretroviral therapy (cART) suppresses viral replication but does not cure HIV infection because a reservoir of infectious (intact) HIV proviruses persists in long-lived CD4+T cells. However, a large majority (>95%) of HIV-infected cells that persist on effective cART carry defective (non-infectious) proviruses. Defective proviruses consisting of only a single LTR (solo long terminal repeat) are commonly found as endogenous retroviruses in many animal species, but the frequency of solo-LTR HIV proviruses has not been well defined. Here we show that, in five pediatric donors whose viremia was suppressed on cART for at least 5 years, the proviruses in the nine largest clones of HIV-infected cells were solo LTRs. The sizes of five of these clones were assayed longitudinally by integration site-specific quantitative PCR. Minor waxing and waning of the clones was observed, suggesting that these clones are generally stable over time. Our findings show that solo LTRs comprise a large fraction of the proviruses in infected cell clones that persist in children on long-term cART. IMPORTANCE This work highlights that severely deleted HIV-1 proviruses comprise a significant proportion of the proviral landscape and are often overlooked.

Identification, comprehensive characterization, and comparative genomics of the HERV-K(HML8) integrations in the human genome

Around 8% of the human genome is composed by Human Endogenous Retroviruses (HERVs), ancient viral sequences inherited from the primate germ line after their infection by now extinct retroviruses. Given the still underexplored physiological and pathological roles of HERVs, it is fundamental to increase our information about the genomic composition of the different groups, to lay reliable foundation for functional studies. Among HERVs, the most characterized elements belong to the beta-like superfamily HERV-K, comprising 10 groups (HML1-10) with HML2 being the most recent and studied one. Among HMLs, the HML8 group is the only one still lacking a comprehensive genomic description. In the present work, we investigated HML8 sequences' distribution in the human genome (GRCh38/hg38), identifying 23 novel proviruses and characterizing the overall 78 HML8 proviruses in terms of genome structure, phylogeny, and integration pattern. HML8 elements were significantly enriched in human chromosomes 8 and X (p<0.005) while chromosomes 17 and 20 showed fewer integrations than expected (p<0.025 and p<0.005, respectively). Phylogenetic analyses classified HML8 members into 3 clusters, corresponding to the three LTR types MER11A, MER11B and MER11C. Besides different LTR types, common signatures in the internal structure suggested the potential existence of three different ancestral HML8 variants. Accordingly, time of integration estimation coupled with comparative genomics revealed that these three clusters have a different time of integration in the primates' genome, with MER11C elements being significantly younger than MER11A- and MER11B associated proviruses (p<0.005 and p<0.05, respectively). Approximately 30% of the HML8 elements were found co-localized within human genes, sometimes in exonic portions and with the same orientation, deserving further studies for their possible effects on gene expression. Overall, we provide the first detailed picture of the HML8 group distribution and variety among the genome, creating the backbone for the specific analysis of their transcriptional activity in healthy and diseased conditions.

Human Endogenous Retroviruses in Diseases

Human endogenous retroviruses (HERVs), which are conserved sequences of ancient retroviruses, are widely distributed in the human genome. Although most HERVs have been rendered inactive by evolution, some have continued to exhibit important cytological functions. HERVs in the human genome perform dual functions: on the one hand, they are involved in important physiological processes such as placental development and immune regulation; on the other hand, their aberrant expression is closely associated with the pathological processes of several diseases, such as cancers, autoimmune diseases, and viral infections. HERVs can also regulate a variety of host cellular functions, including the expression of protein-coding genes and regulatory elements that have evolved from HERVs. Here, we present recent research on the roles of HERVs in viral infections and cancers, including the dysregulation of HERVs in various viral infections, HERV-induced epigenetic modifications of histones (such as methylation and acetylation), and the potential mechanisms of HERV-mediated antiviral immunity. We also describe therapies to improve the efficacy of vaccines and medications either by directly or indirectly targeting HERVs, depending on the HERV.

Centromere defects, chromosome instability, and cGAS-STING activation in systemic sclerosis

Centromere defects in Systemic Sclerosis (SSc) have remained unexplored despite the fact that many centromere proteins were discovered in patients with SSc. Here we report that lesion skin fibroblasts from SSc patients show marked alterations in centromeric DNA. SSc fibroblasts also show DNA damage, abnormal chromosome segregation, aneuploidy (only in diffuse cutaneous (dcSSc)) and micronuclei (in all types of SSc), some of which lose centromere identity while retaining centromere DNA sequences. Strikingly, we find cytoplasmic "leaking" of centromere proteins in limited cutaneous SSc (lcSSc) fibroblasts. Cytoplasmic centromere proteins co-localize with antigen presenting MHC Class II molecules, which correlate precisely with the presence of anti-centromere antibodies. CENPA expression and micronuclei formation correlate highly with activation of the cGAS-STING/IFN-β pathway as well as markers of reactive oxygen species (ROS) and fibrosis, ultimately suggesting a link between centromere alterations, chromosome instability, SSc autoimmunity, and fibrosis.

Modulation of HERV Expression by Four Different Encephalitic Arboviruses during Infection of Human Primary Astrocytes

Human retroelements (HERVs) are retroviral origin sequences fixed in the human genome. HERVs induction is associated with neurogenesis, cellular development, immune activation, and neurological disorders. Arboviruses are often associated with the development of encephalitis. The interplay between these viruses and HERVs has not been fully elucidated. In this work, we analyzed RNAseq data derived from infected human primary astrocytes by Zika (ZikV), Mayaro (MayV), Oropouche (OroV) and Chikungunya (ChikV) viruses, and evaluated the modulation of HERVs and their nearby genes. Our data show common HERVs expression modulation by both alphaviruses, suggesting conserved evolutionary routes of transcription regulation. A total of 15 HERVs were co-modulated by the four arboviruses, including the highly upregulated HERV4_4q22. Data on the upregulation of genes nearby to these elements in ChikV, MayV and OroV infections were also obtained, and interaction networks were built. The upregulation of 14 genes common among all viruses was observed in the networks, and 93 genes between MayV and ChikV. These genes are related to cellular processes such as cellular replication, cytoskeleton, cell vesicle traffic and antiviral response. Together, our results support the role of HERVs induction in the transcription regulation process of genes during arboviral infections.

Characterizing mobile element insertions in 5675 genomes

Mobile element insertions (MEIs) are a major class of structural variants (SVs) and have been linked to many human genetic disorders, including hemophilia, neurofibromatosis, and various cancers. However, human MEI resources from large-scale genome sequencing are still lacking compared to those for SNPs and SVs. Here, we report a comprehensive map of 36 699 non-reference MEIs constructed from 5675 genomes, comprising 2998 Chinese samples (∼26.2×, NyuWa) and 2677 samples from the 1000 Genomes Project (∼7.4×, 1KGP). We discovered that LINE-1 insertions were highly enriched in centromere regions, implying the role of chromosome context in retroelement insertion. After functional annotation, we estimated that MEIs are responsible for about 9.3% of all protein-truncating events per genome. Finally, we built a companion database named HMEID for public use. This resource represents the latest and largest genomewide study on MEIs and will have broad utility for exploration of human MEI findings.

HERV-K(HML7) Integrations in the Human Genome: Comprehensive Characterization and Comparative Analysis in Non-Human Primates

Simple Summary

The human genome is not human at all, but it includes a multitude of sequences inherited from ancient viral infections that affected primates’ germ line. These elements can be seen as the fossils of now-extinct retroviruses, and are called Human Endogenous Retroviruses (HERVs). View as “junk DNA” for a long time, HERVs constitute 4 times the amount of DNA needed to produce all cellular proteins, and growing evidence indicates their crucial role in primate brain evolution, placenta development, and innate immunity shaping. HERVs are also intensively studied for a pathological role, even if the incomplete knowledge about their exact number and genomic position has thus far prevented any causal association. Among possible relevant HERVs, the HERV-K supergroup is of particular interest, including some of the oldest (HML5) as well as youngest (HML2) integrations. Among HERV-Ks, the HML7 group still lack a detailed description, and the present work thus aimed to identify and characterize all HML7 elements in the human genome. Results showed that the HML7 group includes 23 elements and an additional 160 “scars” of past infection that invaded in primates mostly between 20 and 30 million years ago, providing an exhaustive background to study their impact on human pathophysiology.

Abstract

Endogenous Retroviruses (ERVs) are ancient relics of infections that affected the primate germ line and constitute about 8% of our genome. Growing evidence indicates that ERVs had a major role in vertebrate evolution, being occasionally domesticated by the host physiology. In addition, human ERV (HERV) expression is highly investigated for a possible pathological role, even if no clear associations have been reported yet. In fact, on the one side, the study of HERV expression in high-throughput data is a powerful and promising tool to assess their actual dysregulation in diseased conditions; but, on the other side, the poor knowledge about the various HERV group genomic diversity and individual members somehow prevented the association between specific HERV loci and a given molecular mechanism of pathogenesis. The present study is focused on the HERV-K(HML7) group that—differently from the other HERV-K members—still remains poorly characterized. Starting from an initial identification performed with the software RetroTector, we collected 23 HML7 proviral insertions and about 160 HML7 solitary LTRs that were analyzed in terms of genomic distribution, revealing a significant enrichment in chromosome X and the frequent localization within human gene introns as well as in pericentromeric and centromeric regions. Phylogenetic analyses showed that HML7 members form a monophyletic group, which based on age estimation and comparative localization in non-human primates had its major diffusion between 20 and 30 million years ago. Structural characterization revealed that besides 3 complete HML7 proviruses, the other group members shared a highly defective structure that, however, still presents recognizable functional domains, making it worth further investigation in the human population to assess the presence of residual coding potential.

Human Endogenous Retrovirus Expression Is Upregulated in the Breast Cancer Microenvironment of HIV Infected Women: A Pilot Study

In people living with HIV (PLWH), chronic inflammation can lead to cancer initiation and progression, besides driving a dysregulated and diminished immune responsiveness. HIV infection also leads to increased transcription of Human Endogenous Retroviruses (HERVs), which could increase an inflammatory environment and create a tumor growth suppressive environment with high expression of pro-inflammatory cytokines. In order to determine the impact of HIV infection to HERV expression on the breast cancer microenvironment, we sequenced total RNA from formalin-fixed paraffin-embedded (FFPE) breast cancer samples of women HIV-negative and HIV-positive for transcriptome and retrotranscriptome analyses. We performed RNA extraction from FFPE samples, library preparation and total RNA sequencing (RNA-seq). The RNA-seq analysis shows 185 differentially expressed genes: 181 host genes (178 upregulated and three downregulated) and four upregulated HERV transcripts in HIV-positive samples. We also explored the impact of HERV expression in its neighboring breast cancer development genes (BRCA1, CCND1, NBS1/NBN, RAD50, KRAS, PI3K/PIK3CA) and in long non-coding RNA expression (AC060780.1, also known as RP11-242D8.1). We found a significant positive association of HERV expression with RAD50 and with AC060780.1, which suggest a possible role of HERV in regulating breast cancer genes from PLWH with breast cancer. In addition, we found immune system, extracellular matrix organization and metabolic signaling genes upregulated in HIV-positive breast cancer. In conclusion, our findings provide evidence of transcriptional and retrotranscriptional changes in breast cancer from PLWH compared to non-HIV breast cancer, including dysregulation of HERVs, suggesting an indirect effect of the virus on the breast cancer microenvironment.

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.

Centromeric Transcription: A Conserved Swiss-Army Knife

In most species, the centromere is comprised of repetitive DNA sequences, which rapidly evolve. Paradoxically, centromeres fulfill an essential function during mitosis, as they are the chromosomal sites wherein, through the kinetochore, the mitotic spindles bind. It is now generally accepted that centromeres are transcribed, and that such transcription is associated with a broad range of functions. More than a decade of work on this topic has shown that centromeric transcripts are found across the eukaryotic tree and associate with heterochromatin formation, chromatin structure, kinetochore structure, centromeric protein loading, and inner centromere signaling. In this review, we discuss the conservation of small and long non-coding centromeric RNAs, their associations with various centromeric functions, and their potential roles in disease.

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.

Genetic and epigenetic effects on centromere establishment

Endogenous chromosomes contain centromeres to direct equal chromosomal segregation in mitosis and meiosis. The location and function of existing centromeres is usually maintained through cell cycles and generations. Recent studies have investigated how the centromere-specific histone H3 variant CENP-A is assembled and replenished after DNA replication to epigenetically propagate the centromere identity. However, existing centromeres occasionally become inactivated, with or without change in underlying DNA sequences, or lost after chromosomal rearrangements, resulting in acentric chromosomes. New centromeres, known as neocentromeres, may form on ectopic, non-centromeric chromosomal regions to rescue acentric chromosomes from being lost, or form dicentric chromosomes if the original centromere is still active. In addition, de novo centromeres can form after chromatinization of purified DNA that is exogenously introduced into cells. Here, we review the phenomena of naturally occurring and experimentally induced new centromeres and summarize the genetic (DNA sequence) and epigenetic features of these new centromeres. We compare the characteristics of new and native centromeres to understand whether there are different requirements for centromere establishment and propagation. Based on our understanding of the mechanisms of new centromere formation, we discuss the perspectives of developing more stably segregating human artificial chromosomes to facilitate gene delivery in therapeutics and research.

The Genomic Landscape of Centromeres in Cancers

Centromere genomics remain poorly characterized in cancer, due to technologic limitations in sequencing and bioinformatics methodologies that make high-resolution delineation of centromeric loci difficult to achieve. We here leverage a highly specific and targeted rapid PCR methodology to quantitatively assess the genomic landscape of centromeres in cancer cell lines and primary tissue. PCR-based profiling of centromeres revealed widespread heterogeneity of centromeric and pericentromeric sequences in cancer cells and tissues as compared to healthy counterparts. Quantitative reductions in centromeric core and pericentromeric markers (α-satellite units and HERV-K copies) were observed in neoplastic samples as compared to healthy counterparts. Subsequent phylogenetic analysis of a pericentromeric endogenous retrovirus amplified by PCR revealed possible gene conversion events occurring at numerous pericentromeric loci in the setting of malignancy. Our findings collectively represent a more comprehensive evaluation of centromere genetics in the setting of malignancy, providing valuable insight into the evolution and reshuffling of centromeric sequences in cancer development and progression.

Retrovirus insertion site analysis of LGL leukemia patient genomes

Background

Large granular lymphocyte (LGL) leukemia is an uncommon cancer characterized by sustained clonal proliferation of LGL cells. Antibodies reactive to retroviruses have been documented in the serum of patients with LGL leukemia. Culture or molecular approaches have to date not been successful in identifying a retrovirus.

Methods

Because a retrovirus must integrate into the genome of an infected cell, we focused our efforts on detecting a novel retrovirus integration site in the clonally expanded LGL cells. We present a new computational tool that uses long-insert mate pair sequence data to search the genome of LGL leukemia cells for retrovirus integration sites. We also utilize recently published methods to interrogate the status of polymorphic human endogenous retrovirus type K (HERV-K) provirus in patient genomes.

Results

Our data show that there are no new retrovirus insertions in LGL genomes of LGL leukemia patients. However, our insertion call tool did detect four HERV-K provirus integration sites that are polymorphic in the human population but absent from the human reference genome, hg19. To determine if the prevalence of these or other polymorphic proviral HERV-Ks differed between LGL leukemia patients and the general population, we used a recently developed tool that reports sites in the human genome occupied by a known proviral HERV-K. We report that there are significant differences in the number of polymorphic HERV-Ks in the genomes of LGL leukemia patients of European origin compared to individuals with European ancestry in the 1000 genomes (KGP) data.

Conclusions

Our study confirms that the clonal expansion of LGL cells in LGL leukemia is not driven by the integration of a new infectious or endogenous retrovirus, although we do not rule out that these cells are responding to retroviral antigens produced in other cell types. However, our computational analyses revealed that the genomes of LGL leukemia patients carry a higher burden of polymorphic HERV-K proviruses compare to individuals from KGP of European ancestry. Our research emphasizes the merits of comprehensive genomic assessment of HERV-K in cancer samples and suggests that further analyses to determine contributions of HERV-K to LGL leukemia are warranted.

Electronic supplementary material

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

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.

Centromere Repeats: Hidden Gems of the Genome

Satellite DNAs are now regarded as powerful and active contributors to genomic and chromosomal evolution. Paired with mobile transposable elements, these repetitive sequences provide a dynamic mechanism through which novel karyotypic modifications and chromosomal rearrangements may occur. In this review, we discuss the regulatory activity of satellite DNA and their neighboring transposable elements in a chromosomal context with a particular emphasis on the integral role of both in centromere function. In addition, we discuss the varied mechanisms by which centromeric repeats have endured evolutionary processes, producing a novel, species-specific centromeric landscape despite sharing a ubiquitously conserved function. Finally, we highlight the role these repetitive elements play in the establishment and functionality of de novo centromeres and chromosomal breakpoints that underpin karyotypic variation. By emphasizing these unique activities of satellite DNAs and transposable elements, we hope to disparage the conventional exemplification of repetitive DNA in the historically-associated context of ‘junk’.

The HERV-K accessory protein Np9 controls viability and migration of teratocarcinoma cells

Human endogenous retroviruses are remnants of ancient germline infections that make up approximately 8% of the modern human genome. The HERV-K (HML-2) family is one of the most recent entrants into the human germline, these viruses appear to be transcriptionally active, and HERV-K viral like particles (VLPs) are found in cell lines from a number of human malignancies. HERV-K VLPs were first found to be produced in teratocarcinoma cell lines, and since then teratocarcinoma has been thought of as the classical model for HERV-Ks, with the NCCIT teratocarcinoma cell line particularly known to produce VLPs. Treatment for teratocarcinoma has progressed since its discovery, with improved prognosis for patients. Since the introduction of platinum based therapy, first year survival has greatly improved even with disseminated disease; however, it is estimated that 20% to 30% of patients present with metastatic germ cell tumor relapse following initial treatments. Also, the toxicity associated with the use of chemotherapeutic agents used to treat germ cell tumors is still a major concern. In this study, we show that the depletion of the HERV-K accessory protein Np9 increases the sensitivity of NCCIT teratocarcinoma cells to bleomycin and cisplatin. While decreasing the expression of Np9 had only a modest effect on the baseline viability of the cells, the reduced expression of Np9 increased the sensitivity of the teratocarcinoma cells to environmental (serum starvation) and chemical (chemotherapeutic) stresses. Np9 is also essential to the migration of NCCIT teratocarcinoma cells: in a wound closure assay, reduced expression of Np9 resulted in cells migrating into the wound at a slower rate, whereas reintroduction of Np9 resulted in NCCIT cells migrating back into the wound in a manner similar to the control. These findings support the implication that the HERV-K accessory protein Np9 has oncogenic potential.

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.

Transcriptional profiling of HERV-K(HML-2) in amyotrophic lateral sclerosis and potential implications for expression of HML-2 proteins

Background

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder. About 90% of ALS cases are without a known genetic cause. The human endogenous retrovirus multi-copy HERV-K(HML-2) group was recently reported to potentially contribute to neurodegeneration and disease pathogenesis in ALS because of transcriptional upregulation and toxic effects of HML-2 Envelope (Env) protein. Env and other proteins are encoded by some transcriptionally active HML-2 loci. However, more detailed information is required regarding which HML-2 loci are transcribed in ALS, which of their proteins are expressed, and differences between the disease and non-disease states.

Methods

For brain and spinal cord tissue samples from ALS patients and controls, we identified transcribed HML-2 loci by generating and mapping HML-2-specific cDNA sequences. We predicted expression of HML-2 env gene-derived proteins based on the observed cDNA sequences. Furthermore, we determined overall HML-2 transcript levels by RT-qPCR and investigated presence of HML-2 Env protein in ALS and control tissue samples by Western blotting.

Results

We identified 24 different transcribed HML-2 loci. Some of those loci are transcribed at relatively high levels. However, significant differences in HML-2 loci transcriptional activities were not seen when comparing ALS and controls. Likewise, overall HML-2 transcript levels, as determined by RT-qPCR, were not significantly different between ALS and controls. Indeed, we were unable to detect full-length HML-2 Env protein in ALS and control tissue samples despite reasonable sensitivity. Rather our analyses suggest that a number of HML-2 protein variants other than full-length Env may potentially be expressed in ALS patients.

Conclusions

Our results expand and refine recent publications on HERV-K(HML-2) and ALS. Some of our results are in conflict with recent findings and call for further specific analyses. Our profiling of HML-2 transcription in ALS opens up the possibility that HML-2 proteins other than canonical full-length Env may have to be considered when studying the role of HML-2 in ALS disease.

Electronic supplementary material

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

LTR-retrotransposon transcriptome modulation in response to endotoxin-induced stress in PBMCs

Background

Human Endogenous Retroviruses (HERVs) and Mammalian apparent LTR-retrotransposons (MaLRs) represent the 8% of our genome and are distributed among our 46 chromosomes. These LTR-retrotransposons are thought to be essentially silent except in cancer, autoimmunity and placental development. Their Long Terminal Repeats (LTRs) constitute putative promoter or polyA regulatory sequences. In this study, we used a recently described high-density microarray which can be used to study HERV/MaLR transcriptome including 353,994 HERV/MaLR loci and 1559 immunity-related genes.

Results

We described, for the first time, the HERV transcriptome in peripheral blood mononuclear cells (PBMCs) using a cellular model mimicking inflammatory response and monocyte anergy observed after septic shock. About 5.6% of the HERV/MaLR repertoire is transcribed in PBMCs. Roughly one-tenth [5.7–13.1%] of LTRs exhibit a putative constitutive promoter or polyA function while one-quarter [19.5–27.6%] may shift from silent to active. Evidence was given that some HERVs/MaLRs and genes may share similar regulation control under lipopolysaccharide (LPS) stimulation conditions. Stimulus-dependent response confirms that HERV expression is tightly regulated in PBMCs. Altogether, these observations make it possible to integrate 62 HERVs/MaLRs and 26 genes in 11 canonical pathways and suggest a link between HERV expression and immune response. The transcriptional modulation of HERVs located close to genes such as OAS2/3 and IFI44/IFI44L or at a great distance from genes was discussed.

Conclusion

This microarray-based approach revealed the expression of about 47,466 distinct HERV loci and identified 951 putative promoter LTRs and 744 putative polyA LTRs in PBMCs. HERV/MaLR expression was shown to be tightly modulated under several stimuli including high-dose and low-dose LPS and Interferon-γ (IFN-γ). HERV incorporation at the crossroads of immune response pathways paves the way for further functional studies and analyses of the HERV transcriptome in altered immune responses in vivo such as in sepsis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4901-9) contains supplementary material, which is available to authorized users.

Transcriptional Modulation of Human Endogenous Retroviruses in Primary CD4+ T Cells Following Vorinostat Treatment

The greatest obstacle to a cure for HIV is the provirus that integrates into the genome of the infected cell and persists despite antiretroviral therapy. A "shock and kill" approach has been proposed as a strategy for an HIV cure whereby drugs and compounds referred to as latency-reversing agents (LRAs) are used to "shock" the silent provirus into active replication to permit "killing" by virus-induced pathology or immune recognition. The LRA most utilized to date in clinical trials has been the histone deacetylase (HDAC) inhibitor-vorinostat. Potentially, pathological off-target effects of vorinostat may result from the activation of human endogenous retroviruses (HERVs), which share common ancestry with exogenous retroviruses including HIV. To explore the effects of HDAC inhibition on HERV transcription, an unbiased pharmacogenomics approach (total RNA-Seq) was used to evaluate HERV expression following the exposure of primary CD4+ T cells to a high dose of vorinostat. Over 2,000 individual HERV elements were found to be significantly modulated by vorinostat, whereby elements belonging to the ERVL family (e.g., LTR16C and LTR33) were predominantly downregulated, in contrast to LTR12 elements of the HERV-9 family, which exhibited the greatest signal, with the upregulation of 140 distinct elements. The modulation of three different LTR12 elements by vorinostat was confirmed by droplet digital PCR along a dose-response curve. The monitoring of LTR12 expression during clinical trials with vorinostat may be indicated to assess the impact of this HERV on the human genome and host immunity.

Transposable elements: genome innovation, chromosome diversity, and centromere conflict

Although it was nearly 70 years ago when transposable elements (TEs) were first discovered “jumping” from one genomic location to another, TEs are now recognized as contributors to genomic innovations as well as genome instability across a wide variety of species. In this review, we illustrate the ways in which active TEs, specifically retroelements, can create novel chromosome rearrangements and impact gene expression, leading to disease in some cases and species-specific diversity in others. We explore the ways in which eukaryotic genomes have evolved defense mechanisms to temper TE activity and the ways in which TEs continue to influence genome structure despite being rendered transpositionally inactive. Finally, we focus on the role of TEs in the establishment, maintenance, and stabilization of critical, yet rapidly evolving, chromosome features: eukaryotic centromeres. Across centromeres, specific types of TEs participate in genomic conflict, a balancing act wherein they are actively inserting into centromeric domains yet are harnessed for the recruitment of centromeric histones and potentially new centromere formation.

Application of viromics: a new approach to the understanding of viral infections in humans

This review is focused at exploring the strengths of modern technology driven data compiled in the areas of virus gene sequencing, virus protein structures and their implication to viral diagnosis and therapy. The information for virome analysis (viromics) is generated by the study of viral genomes (entire nucleotide sequence) and viral genes (coding for protein). Presently, the study of viral infectious diseases in terms of etiopathogenesis and development of newer therapeutics is undergoing rapid changes. Currently, viromics relies on deep sequencing, next generation sequencing (NGS) data and public domain databases like GenBank and unique virus specific databases. Two commonly used NGS platforms: Illumina and Ion Torrent, recommend maximum fragment lengths of about 300 and 400 nucleotides for analysis respectively. Direct detection of viruses in clinical samples is now evolving using these methods. Presently, there are a considerable number of good treatment options for HBV/HIV/HCV. These viruses however show development of drug resistance. The drug susceptibility regions of the genomes are sequenced and the prediction of drug resistance is now possible from 3 public domains available on the web. This has been made possible through advances in the technology with the advent of high throughput sequencing and meta-analysis through sophisticated and easy to use software and the use of high speed computers for bioinformatics. More recently NGS technology has been improved with single-molecule real-time sequencing. Here complete long reads can be obtained with less error overcoming a limitation of the NGS which is inherently prone to software anomalies that arise in the hands of personnel without adequate training. The development in understanding the viruses in terms of their genome, pathobiology, transcriptomics and molecular epidemiology constitutes viromics. It could be stated that these developments will bring about radical changes and advancement especially in the field of antiviral therapy and diagnostic virology.

Rapid molecular assays to study human centromere genomics

The centromere is the structural unit responsible for the faithful segregation of chromosomes. Although regulation of centromeric function by epigenetic factors has been well-studied, the contributions of the underlying DNA sequences have been much less well defined, and existing methodologies for studying centromere genomics in biology are laborious. We have identified specific markers in the centromere of 23 of the 24 human chromosomes that allow for rapid PCR assays capable of capturing the genomic landscape of human centromeres at a given time. Use of this genetic strategy can also delineate which specific centromere arrays in each chromosome drive the recruitment of epigenetic modulators. We further show that, surprisingly, loss and rearrangement of DNA in centromere 21 is associated with trisomy 21. This new approach can thus be used to rapidly take a snapshot of the genetics and epigenetics of each specific human centromere in nondisjunction disorders and other biological settings.

Susceptibility of Human Endogenous Retrovirus Type K to Reverse Transcriptase Inhibitors

Human endogenous retroviruses (HERVs) make up 8% of the human genome. The HERV type K (HERV-K) HML-2 (HK2) family contains proviruses that are the most recent entrants into the human germ line and are transcriptionally active. In HIV-1 infection and cancer, HK2 genes produce retroviral particles that appear to be infectious, yet the replication capacity of these viruses and potential pathogenicity has been difficult to ascertain. In this report, we screened the efficacy of commercially available reverse transcriptase inhibitors (RTIs) at inhibiting the enzymatic activity of HK2 RT and HK2 genomic replication. Interestingly, only one provirus, K103, was found to encode a functional RT among those examined. Several nucleoside analogue RTIs (NRTIs) blocked K103 RT activity and consistently inhibited the replication of HK2 genomes. The NRTIs zidovudine (AZT), stavudine (d4T), didanosine (ddI), and lamivudine (3TC), and the nucleotide RTI inhibitor tenofovir (TDF), show efficacy in blocking K103 RT. HIV-1-specific nonnucleoside RTIs (NNRTIs), protease inhibitors (PIs), and integrase inhibitors (IIs) did not affect HK2, except for the NNRTI etravirine (ETV). The inhibition of HK2 infectivity by NRTIs appears to take place at either the reverse transcription step of the viral genome prior to HK2 viral particle formation and/or in the infected cells. Inhibition of HK2 by these drugs will be useful in suppressing HK2 infectivity if these viruses prove to be pathogenic in cancer, neurological disorders, or other diseases associated with HK2. The present studies also elucidate a key aspect of the life cycle of HK2, specifically addressing how they do, and/or did, replicate.IMPORTANCE Endogenous retroviruses are relics of ancestral virus infections in the human genome. The most recent of these infections was caused by HK2. While HK2 often remains silent in the genome, this group of viruses is activated in HIV-1-infected and cancer cells. Recent evidence suggests that these viruses are infectious, and the potential exists for HK2 to contribute to disease. We show that HK2, and specifically the enzyme that mediates virus replication, can be inhibited by a panel of drugs that are commercially available. We show that several drugs block HK2 with different efficacies. The inhibition of HK2 replication by antiretroviral drugs appears to occur in the virus itself as well as after infection of cells. Therefore, these drugs might prove to be an effective treatment by suppressing HK2 infectivity in diseases where these viruses have been implicated, such as cancer and neurological syndromes.

Molecular mechanisms by which HERV-K Gag interferes with HIV-1 Gag assembly and particle infectivity

Background

Human endogenous retroviruses (HERVs), the remnants of ancient retroviral infections, constitute approximately 8% of human genomic DNA. Since HERV-K Gag expression is induced by HIV-1 Tat in T cells, induced HERV-K proteins could affect HIV-1 replication. Indeed, previously we showed that HERV-K Gag and HIV-1 Gag coassemble and that this appears to correlate with the effect of HERV-K Gag expression on HIV-1 particle release and its infectivity. We further showed that coassembly requires both MA and NC domains, which presumably serve as scaffolding for Gag via their abilities to bind membrane and RNA, respectively. Notably, however, despite possessing these abilities, MLV Gag failed to coassemble with HIV-1 Gag and did not affect assembly and infectivity of HIV-1 particles. It is unclear how the specificity of coassembly is determined.

Results

Here, we showed that coexpression of HERV-K Gag with HIV-1 Gag changed size and morphology of progeny HIV-1 particles and severely diminished infectivity of such progeny viruses. We further compared HERV-K-MLV chimeric constructs to identify molecular determinants for coassembly specificity and for inhibition of HIV-1 release efficiency and infectivity. We found that the CA N-terminal domain (NTD) of HERV-K Gag is important for the reduction of the HIV-1 release efficiency, whereas both CA-NTD and major homology region of HERV-K Gag contribute to colocalization with HIV-1 Gag. Interestingly, these regions of HERV-K Gag were not required for reduction of progeny HIV-1 infectivity.

Conclusions

Our results showed that HERV-K Gag CA is important for reduction of HIV-1 release and infectivity but the different regions within CA are involved in the effects on the HIV-1 release and infectivity. Altogether, these findings revealed that HERV-K Gag interferes the HIV-1 replication by two distinct molecular mechanisms.

Electronic supplementary material

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

A comprehensive hybridization model allows whole HERV transcriptome profiling using high density microarray

Background

Human endogenous retroviruses (HERVs) have received much attention for their implications in the etiology of many human diseases and their profound effect on evolution. Notably, recent studies have highlighted associations between HERVs expression and cancers (Yu et al., Int J Mol Med 32, 2013), autoimmunity (Balada et al., Int Rev Immunol 29:351–370, 2010) and neurological (Christensen, J Neuroimmune Pharmacol 5:326–335, 2010) conditions. Their repetitive nature makes their study particularly challenging, where expression studies have largely focused on individual loci (De Parseval et al., J Virol 77:10414–10422, 2003) or general trends within families (Forsman et al., J Virol Methods 129:16–30, 2005; Seifarth et al., J Virol 79:341–352, 2005; Pichon et al., Nucleic Acids Res 34:e46, 2006).

Methods

To refine our understanding of HERVs activity, we introduce here a new microarray, HERV-V3. This work was made possible by the careful detection and annotation of genomic HERV/MaLR sequences as well as the development of a new hybridization model, allowing the optimization of probe performances and the control of cross-reactions.

Results

HERV-V3 offers an almost complete coverage of HERVs and their ancestors (mammalian apparent LTR-retrotransposons, MaLRs) at the locus level along with four other repertoires (active LINE-1 elements, lncRNA, a selection of 1559 human genes and common infectious viruses). We demonstrate that HERV-V3 analytical performances are comparable with commercial Affymetrix arrays, and that for a selection of tissue/pathological specific loci, the patterns of expression measured on HERV-V3 is consistent with those reported in the literature.

Conclusions

Given its large HERVs/MaLRs coverage and additional repertoires, HERV-V3 opens the door to multiple applications such as enhancers and alternative promoters identification, biomarkers identification as well as the characterization of genes and HERVs/MaLRs modulation caused by viral infection.

Electronic supplementary material

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

A contaminant-free assessment of Endogenous Retroviral RNA in human plasma

Endogenous retroviruses (ERVs) comprise 6–8% of the human genome. HERVs are silenced in most normal tissues, up-regulated in stem cells and in placenta but also in cancer and HIV-1 infection. Crucially, there are conflicting reports on detecting HERV RNA in non-cellular clinical samples such as plasma that suggest the study of HERV RNA can be daunting. Indeed, we find that the use of real-time PCR in a quality assured clinical laboratory setting can be sensitive to low-level proviral contamination. We developed a mathematical model for low-level contamination that allowed us to design a laboratory protocol and standard operating procedures for robust measurement of HERV RNA. We focus on one family, HERV-K HML-2 (HK2) that has been most recently active even though they invaded our ancestral genomes almost 30 millions ago. We extensively validated our experimental design on a model cell culture system showing high sensitivity and specificity, totally eliminating the proviral contamination. We then tested 236 plasma samples from patients infected with HIV-1, HCV or HBV and found them to be negative. The study of HERV RNA for human translational studies should be performed with extensively validated protocols and standard operating procedures to control the widespread low-level human DNA contamination.

The role of molecular mimicry and other factors in the association of Human Endogenous Retroviruses and autoimmunity

Human Endogenous Retroviruses (HERVs) have been implicated in autoimmune and other diseases. Molecular mimicry has been postulated as a potential mechanism of autoimmunity. Exogenous viruses have also been reported to be associated with the same diseases, as have genetic and environmental factors. If molecular mimicry were to be shown to be an initiating mechanism of some autoimmune diseases, then therapeutic options of blocking antibodies and peptides might be of benefit in halting diseases at the outset. Bioinformatic and molecular modelling techniques have been employed to investigate molecular mimicry and the evidence for the association of HERVs and autoimmunity is reviewed. The most convincing evidence for molecular mimicry is in rheumatoid arthritis, where HERV K-10 shares amino acid sequences with IgG1Fc, a target for rheumatoid factor. Systemic lupus erythematosus is an example of a condition associated with several autoantibodies, and several endogenous and exogenous viruses have been reported to be associated with the disease. The lack of a clear link between one virus and this condition, and the spectrum of clinical manifestations, suggests that genetic, environmental and the inflammatory response to a virus or viruses might also be major factors in the pathogenesis of lupus and other autoimmune conditions. Where there are strong associations between a virus and an autoimmune condition, such as in hepatitis C and cryoglobulinaemia, the use of bioinformatics and molecular modelling can also be utilized to help to understand the role of molecular mimicry in how HERVs might trigger disease.

HERV-K(HML-2), a seemingly silent subtenant - but still waters run deep

A large proportion of the human genome consists of endogenous retroviruses, some of which are well preserved, showing transcriptional activity, and expressing retroviral proteins. The HERV-K(HML-2) family represents the most intact members of these elements, with some having open and intact reading frames for viral proteins and the ability to form virus-like particles. Although generally suppressed in most healthy tissues by a variety of epigenetic processes and antiviral mechanisms, there is evidence that some members of this family are (at least partly) still active - particularly in certain stem cells and various tumors. This raises the possibility of their involvement in tumor induction or in developmental processes. In recent years, many new insights into this fascinating field have been attained, and this review focuses on new discoveries about coevolutionary events and intracellular defense mechanisms against HERV-K(HML-2) activity. We also describe what might occur when these mechanisms fail or become modulated by viral proteins or other viruses and discuss the new vistas opened up by the reconstitution of ancestral viral proteins and even complete HML-2 viruses.

Integration and Fixation Preferences of Human and Mouse Endogenous Retroviruses Uncovered with Functional Data Analysis

Endogenous retroviruses (ERVs), the remnants of retroviral infections in the germ line, occupy ~8% and ~10% of the human and mouse genomes, respectively, and affect their structure, evolution, and function. Yet we still have a limited understanding of how the genomic landscape influences integration and fixation of ERVs. Here we conducted a genome-wide study of the most recently active ERVs in the human and mouse genome. We investigated 826 fixed and 1,065 in vitro HERV-Ks in human, and 1,624 fixed and 242 polymorphic ETns, as well as 3,964 fixed and 1,986 polymorphic IAPs, in mouse. We quantitated >40 human and mouse genomic features (e.g., non-B DNA structure, recombination rates, and histone modifications) in ±32 kb of these ERVs' integration sites and in control regions, and analyzed them using Functional Data Analysis (FDA) methodology. In one of the first applications of FDA in genomics, we identified genomic scales and locations at which these features display their influence, and how they work in concert, to provide signals essential for integration and fixation of ERVs. The investigation of ERVs of different evolutionary ages (young in vitro and polymorphic ERVs, older fixed ERVs) allowed us to disentangle integration vs. fixation preferences. As a result of these analyses, we built a comprehensive model explaining the uneven distribution of ERVs along the genome. We found that ERVs integrate in late-replicating AT-rich regions with abundant microsatellites, mirror repeats, and repressive histone marks. Regions favoring fixation are depleted of genes and evolutionarily conserved elements, and have low recombination rates, reflecting the effects of purifying selection and ectopic recombination removing ERVs from the genome. In addition to providing these biological insights, our study demonstrates the power of exploiting multiple scales and localization with FDA. These powerful techniques are expected to be applicable to many other genomic investigations.

A long terminal repeat retrotransposon of Schizosaccharomyces japonicus integrates upstream of RNA pol III transcribed genes

Background

Transposable elements (TEs) are common constituents of centromeres. However, it is not known what causes this relationship. Schizosaccharomyces japonicus contains 10 families of Long Terminal Repeat (LTR)-retrotransposons and these elements cluster in centromeres and telomeres. In the related yeast, Schizosaccharomyces pombe LTR-retrotransposons Tf1 and Tf2 are distributed in the promoter regions of RNA pol II transcribed genes. Sequence analysis of TEs indicates that Tj1 of S. japonicus is related to Tf1 and Tf2, and uses the same mechanism of self-primed reverse transcription. Thus, we wondered why these related retrotransposons localized in different regions of the genome.

Results

To characterize the integration behavior of Tj1 we expressed it in S. pombe. We found Tj1 was active and capable of generating de novo integration in the chromosomes of S. pombe. The expression of Tj1 is similar to Type C retroviruses in that a stop codon at the end of Gag must be present for efficient integration. 17 inserts were sequenced, 13 occurred within 12 bp upstream of tRNA genes and 3 occurred at other RNA pol III transcribed genes. The link between Tj1 integration and RNA pol III transcription is reminiscent of Ty3, an LTR-retrotransposon of Saccharomyces cerevisiae that interacts with TFIIIB and integrates upstream of tRNA genes.

Conclusion

The integration of Tj1 upstream of tRNA genes and the centromeric clustering of tRNA genes in S. japonicus demonstrate that the clustering of this TE in centromere sequences is due to a unique pattern of integration.

Molecular functions of human endogenous retroviruses in health and disease

Human endogenous retroviruses (HERVs) and related genetic elements form 504 distinct families and occupy ~8% of human genome. Recent success of high-throughput experimental technologies facilitated understanding functional impact of HERVs for molecular machinery of human cells. HERVs encode active retroviral proteins, which may exert important physiological functions in the body, but also may be involved in the progression of cancer and numerous human autoimmune, neurological and infectious diseases. The spectrum of related malignancies includes, but not limits to, multiple sclerosis, psoriasis, lupus, schizophrenia, multiple cancer types and HIV. In addition, HERVs regulate expression of the neighboring host genes and modify genomic regulatory landscape, e.g., by providing regulatory modules like transcription factor binding sites (TFBS). Indeed, recent bioinformatic profiling identified ~110,000 regulatory active HERV elements, which formed at least ~320,000 human TFBS. These and other peculiarities of HERVs might have played an important role in human evolution and speciation. In this paper, we focus on the current progress in understanding of normal and pathological molecular niches of HERVs, on their implications in human evolution, normal physiology and disease. We also review the available databases dealing with various aspects of HERV genetics.

Gene activity in primary T cells infected with HIV89.6: intron retention and induction of genomic repeats

Background

HIV infection has been reported to alter cellular gene activity, but published studies have commonly assayed transformed cell lines and lab-adapted HIV strains, yielding inconsistent results. Here we carried out a deep RNA-Seq analysis of primary human T cells infected with the low passage HIV isolate HIV_89.6.

Results

Seventeen percent of cellular genes showed altered activity 48 h after infection. In a meta-analysis including four other studies, our data differed from studies of HIV infection in cell lines but showed more parallels with infections of primary cells. We found a global trend toward retention of introns after infection, suggestive of a novel cellular response to infection. HIV_89.6 infection was also associated with activation of several human endogenous retroviruses (HERVs) and retrotransposons, of interest as possible novel antigens that could serve as vaccine targets. The most highly activated group of HERVs was a subset of the ERV-9. Analysis showed that activation was associated with a particular variant of ERV-9 long terminal repeats that contains an indel near the U3-R border. These data also allowed quantification of >70 splice forms of the HIV_89.6 RNA and specified the main types of chimeric HIV_89.6-host RNAs. Comparison to over 100,000 integration site sequences from the same infected cell populations allowed quantification of authentic versus artifactual chimeric reads, showing that 5′ read-in, splicing out of HIV_89.6 from the D4 donor and 3′ read-through were the most common HIV_89.6-host cell chimeric RNA forms.

Conclusions

Analysis of RNA abundance after infection of primary T cells with the low passage HIV_89.6 isolate disclosed multiple novel features of HIV-host interactions, notably intron retention and induction of transcription of retrotransposons and endogenous retroviruses.

Electronic supplementary material

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

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.

Expansion of a novel endogenous retrovirus throughout the pericentromeres of modern humans

BACKGROUND

Approximately 8% of the human genome consists of sequences of retroviral origin, a result of ancestral infections of the germ line over millions of years of evolution. The most recent of these infections is attributed to members of the human endogenous retrovirus type-K (HERV-K) (HML-2) family. We recently reported that a previously undetected, large group of HERV-K (HML-2) proviruses, which are descendants of the ancestral K111 infection, are spread throughout human centromeres.

RESULTS

Studying the genomes of certain cell lines and the DNA of healthy individuals that seemingly lack K111, we discover new HERV-K (HML-2) members hidden in pericentromeres of several human chromosomes. All are related through a common ancestor, termed K222, which is a virus that infected the germ line approximately 25 million years ago. K222 exists as a single copy in the genomes of baboons and high order primates, but not New World monkeys, suggesting that progenitor K222 infected the primate germ line after the split between New and Old World monkeys. K222 exists in modern humans at multiple loci spread across the pericentromeres of nine chromosomes, indicating it was amplified during the evolution of modern humans.

CONCLUSIONS

Copying of K222 may have occurred through recombination of the pericentromeres of different chromosomes during human evolution. Evidence of recombination between K111 and K222 suggests that these retroviral sequences have been templates for frequent cross-over events during the process of centromere recombination in humans.

HERV-K(HML-2) rec and np9 transcripts not restricted to disease but present in many normal human tissues

Background

Human endogenous retroviruses of the HERV-K(HML-2) group have been associated with the development of tumor diseases. Various HERV-K(HML-2) loci encode retrovirus-like proteins, and expression of such proteins is upregulated in certain tumor types. HERV-K(HML-2)-encoded Rec and Np9 proteins interact with functionally important cellular proteins and may contribute to tumor development. Though, the biological role of HERV-K(HML-2) transcription and encoded proteins in health and disease is less understood. We therefore investigated transcription specifically of HERV-K(HML-2) rec and np9 mRNAs in a panel of normal human tissues.

Results

We obtained evidence for rec and np9 mRNA being present in all examined 16 normal tissue types. A total of 18 different HERV-K(HML-2) loci were identified as generating rec or np9 mRNA, among them loci not present in the human reference genome and several of the loci harboring open reading frames for Rec or Np9 proteins. Our analysis identified additional alternative splicing events of HERV-K(HML-2) transcripts, some of them encoding variant Rec/Np9 proteins. We also identified a second HERV-K(HML-2) locus formed by L1-mediated retrotransposition that is likewise transcribed in various human tissues.

Conclusions

HERV-K(HML-2) rec and np9 transcripts from different HERV-K(HML-2) loci appear to be present in various normal human tissues. It is conceivable that Rec and Np9 proteins and variants of those proteins are part of the proteome of normal human tissues and thus various cell types. Transcription of HERV-K(HML-2) may thus also have functional relevance in normal human cell physiology.

Electronic supplementary material

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

Regulation of the human endogenous retrovirus K (HML-2) transcriptome by the HIV-1 Tat protein

Approximately 8% of the human genome is made up of endogenous retroviral sequences. As the HIV-1 Tat protein activates the overall expression of the human endogenous retrovirus type K (HERV-K) (HML-2), we used next-generation sequencing to determine which of the 91 currently annotated HERV-K (HML-2) proviruses are regulated by Tat. Transcriptome sequencing of total RNA isolated from Tat- and vehicle-treated peripheral blood lymphocytes from a healthy donor showed that Tat significantly activates expression of 26 unique HERV-K (HML-2) proviruses, silences 12, and does not significantly alter the expression of the remaining proviruses. Quantitative reverse transcription-PCR validation of the sequencing data was performed on Tat-treated PBLs of seven donors using provirus-specific primers and corroborated the results with a substantial degree of quantitative similarity.

IMPORTANCE

The expression of HERV-K (HML-2) is tightly regulated but becomes markedly increased following infection with HIV-1, in part due to the HIV-1 Tat protein. The findings reported here demonstrate the complexity of the genome-wide regulation of HERV-K (HML-2) expression by Tat. This work also demonstrates that although HERV-K (HML-2) proviruses in the human genome are highly similar in terms of DNA sequence, modulation of the expression of specific proviruses in a given biological situation can be ascertained using next-generation sequencing and bioinformatics analysis.

HIV-1 infection leads to increased transcription of human endogenous retrovirus HERV-K (HML-2) proviruses in vivo but not to increased virion production

Recent studies suggest that human endogenous retrovirus group K (HERV-K) provirus expression plays a role in the pathogenesis of HIV-1 infection. In particular, RNA from the HML-2 subgroup of HERV-K proviruses has been reported to be highly expressed at the cellular level and detectable in the plasma of HIV-1-infected patients, suggestive of virion production and, perhaps, replication. In this study, we developed an HML-2-specific quantitative-PCR assay that detects 51 of the 89 known HML-2 proviruses in the human genome. Plasma and peripheral blood mononuclear cells (PBMCs) from HIV-negative controls and HIV-1-infected patients were collected for analysis of HML-2 RNA expression. Contrary to previous reports, we did not detect high levels of HML-2 RNA in the plasma of HIV-1-infected patients, but we did observe a significant increase of HML-2 RNA in total PBMCs compared to HIV-negative controls. The level of HML-2 expression in PBMCs does not appear to be related to patient use of antiretrovirals or to HIV-1 plasma RNA, cellular RNA, or cellular DNA levels. To investigate the source of HML-2 RNA expression, patient PBMCs were sorted into CD3+ CD4+, CD3+ CD8+, CD3- CD14+, and CD3- CD20+ cell subsets and then analyzed for HML-2 RNA levels. No single cell subset was enriched for HML-2 RNA expression in HIV-1-infected patients, but there appears to be substantial variability in the level of HML-2 expression depending on the cell type.

IMPORTANCE

Here, we report that human endogenous retrovirus group K (HERV-K) (HML-2) proviruses are expressed at significantly higher levels in peripheral blood mononuclear cells (PBMCs) from patients with HIV-1 infection than in those from uninfected individuals. However, contrary to previous reports, this expression did not lead to detectable virions in the plasma of these patients. In addition, we found that HML-2 proviruses were expressed in multiple blood cell types from HIV-1-infected individuals, and the magnitude of HML-2 expression was not related to HIV-1 disease markers in this patient cohort. These findings may have implications for HML-2-based therapies targeting HIV-1 infection.

Genomic flexibility of human endogenous retrovirus type K

Human endogenous retrovirus type K (HERV-K) proviruses are scattered throughout the human genome, but as no infectious HERV-K virus has been detected to date, the mechanism by which these viruses replicated and populated the genome remains unresolved. Here, we provide evidence that, in addition to the RNA genomes that canonical retroviruses package, modern HERV-K viruses can contain reverse-transcribed DNA (RT-DNA) genomes. Indeed, reverse transcription of genomic HERV-K RNA into the DNA form is able to occur in three distinct times and locations: (i) in the virus-producing cell prior to viral release, yielding a DNA-containing extracellular virus particle similar to the spumaviruses; (ii) within the extracellular virus particle itself, transitioning from an RNA-containing particle to a DNA-containing particle; and (iii) after entry of the RNA-containing virus into the target cell, similar to canonical retroviruses, such as murine leukemia virus and HIV. Moreover, using a resuscitated HERV-K virus construct, we show that both viruses with RNA genomes and viruses with DNA genomes are capable of infecting target cells. This high level of genomic flexibility historically could have permitted these viruses to replicate in various host cell environments, potentially assisting in their many integration events and resulting in their high prevalence in the human genome. Moreover, the ability of modern HERV-K viruses to proceed through reverse transcription and package RT-DNA genomes suggests a higher level of replication competency than was previously understood, and it may be relevant in HERV-K-associated human diseases.

IMPORTANCE

Retroviral elements comprise at least 8% of the human genome. Of all the endogenous retroviruses, HERV-K viruses are the most intact and biologically active. While a modern infectious HERV-K has yet to be found, HERV-K activation has been associated with cancers, autoimmune diseases, and HIV-1 infection. Thus, determining how this virus family became such a prevalent member of our genome and what it is capable of in its current form are of the utmost importance. Here, we provide evidence that HERV-K viruses currently found in the human genome are able to proceed through reverse transcription and historically utilized a life cycle with a surprising degree of genomic flexibility in which both RNA- and DNA-containing viruses were capable of mediating infection.

Unfixed endogenous retroviral insertions in the human population

One lineage of human endogenous retroviruses (HERVs), HERV-K(HML2), is upregulated in many cancers, some autoimmune/inflammatory diseases, and HIV-infected cells. Despite 3 decades of research, it is not known if these viruses play a causal role in disease, and there has been recent interest in whether they can be used as immunotherapy targets. Resolution of both these questions will be helped by an ability to distinguish between the effects of different integrated copies of the virus (loci). Research so far has concentrated on the 20 or so recently integrated loci that, with one exception, are in the human reference genome sequence. However, this viral lineage has been copying in the human population within the last million years, so some loci will inevitably be present in the human population but absent from the reference sequence. We therefore performed the first detailed search for such loci by mining whole-genome sequences generated by next-generation sequencing. We found a total of 17 loci, and the frequency of their presence ranged from only 2 of the 358 individuals examined to over 95% of them. On average, each individual had six loci that are not in the human reference genome sequence. Comparing the number of loci that we found to an expectation derived from a neutral population genetic model suggests that the lineage was copying until at least ∼250,000 years ago.

IMPORTANCE About 5% of the human genome sequence is composed of the remains of retroviruses that over millions of years have integrated into the chromosomes of egg and/or sperm precursor cells. There are indications that protein expression of these viruses is higher in some diseases, and we need to know (i) whether these viruses have a role in causing disease and (ii) whether they can be used as immunotherapy targets in some of them. Answering both questions requires a better understanding of how individuals differ in the viruses that they carry. We carried out the first careful search for new viruses in some of the many human genome sequences that are now available thanks to advances in sequencing technology. We also compared the number that we found to a theoretical expectation to see if it is likely that these viruses are still replicating in the human population today.