On the impact of mass screening for SARS-CoV-2 through self-testing in Greece.

Background: Screening programs that pre-emptively and routinely test population groups for disease at a massive scale were first implemented during the COVID-19 pandemic in a handful of countries. One of these countries was Greece, which implemented a mass self-testing program during 2021. In contrast to most other non-pharmaceutical interventions (NPIs), mass self-testing programs are particularly attractive for their relatively small financial and social burden, and it is therefore important to understand their effectiveness to inform policy makers and public health officials responding to future pandemics. This study aimed to estimate the number of deaths and hospitalizations averted by the program implemented in Greece and evaluate the impact of several operational decisions. Methods: Granular data from the mass self-testing program deployed by the Greek government between April and December 2021 were obtained. The data were used to fit a novel compartmental model that was developed to describe the dynamics of the COVID-19 pandemic in Greece in the presence of self-testing. The fitted model provided estimates on the effectiveness of the program in averting deaths and hospitalizations. Sensitivity analyses were used to evaluate the impact of operational decisions, including the scale of the program, targeting of sub-populations, and sensitivity (i.e., true positive rate) of tests. Results: Conservative estimates show that the program reduced the reproduction number by 4%, hospitalizations by 25%, and deaths by 20%, translating into approximately 20,000 averted hospitalizations and 2,000 averted deaths in Greece between April and December 2021. Conclusion: Mass self-testing programs are efficient NPIs with minimal social and financial burden; therefore, they are invaluable tools to be considered in pandemic preparedness and response.

Human endogenous retroviruses in cancer: Oncogenesis mechanisms and clinical implications.

Human Endogenous Retroviruses (HERVs) are viral sequences integrated into the human genome, resulting from the infection of human germ-line cells by ancient exogenous retroviruses. Despite losing their replication and retrotransposition abilities, HERVs appear to have been co-opted in human physiological functions while their aberrant expression is linked to human disease. The role of HERVs in multiple malignancies has been demonstrated, however, the extent to which HERV activation and expression participate in the development of cancer is not yet fully comprehended. In this review article, we discuss the presumed role of HERVs in carcinogenesis and their promising diagnostic and prognostic implications. Additionally, we explore recent data on the HERVs in cancer therapeutics, either through the manipulation of their expression, to induce antitumor innate immunity responses or as cancer immunotherapy targets. Finally, more precise and higher resolution high-throughput sequencing approaches will further elucidate HERV participation in human physiological and pathological processes.

Postgenomic Bioinformatic Analysis of Nucleic Acid Sequences Expressed in the Salivary Gland Epithelial Cells of Primary Sjögren's Syndrome Patients in Search of Microorganisms and Endogenous Retroviruses.

Several previous studies from our laboratory have indicated that the salivary gland epithelia of primary Sjögren's syndrome (SS) patients are not only the target of autoimmune immune responses, but also key instigators of the chronic salivary gland inflammatory infiltrates of patients. In particular, the comparative analysis of salivary gland tissue specimens and of in-vitro cultured non-neoplastic salivary gland epithelial cell lines (SGEC, of ductal type) from SS-patients and non-SS disease-controls, have unequivocally highlighted the presence of intrinsic activation in the ductal epithelia of SS-patients and of aberrant expression of inflammagenic molecules thereof, that correlate with the severity of local histopathologic changes, as well as of systemic manifestations of the disease. In the same context, we have recently shown that the ductal epithelia of SS-patients manifest cell-autonomous activation of the AIM2 inflammasome owing to the presence of aberrant cytoplasmic accumulations of damaged DNA. These findings not only provide a mechanistic explanation for the intrinsic activation and inflammatory status of SS ductal epithelia, but may also point towards the putative instigating role of an exogenous or endogenous agent (i.e., a micro-organism or an endogenous retrovirus, respectively). On this basis and to further explore the nature of epithelial cell-intrinsic activation in SS, the present proposal aims to investigate the expression of endogenous retroviral and/or non-human nucleic acid sequences of microbial origin in the ductal salivary gland epithelia of SS-patients, using metagenomic analysis of high throughput DNA and RNA genome sequencing data, which will be obtained from SGEC lines derived from SS-patients and disease-controls.

Immunogenicity and Safety of the BNT162b2 mRNA COVID-19 Vaccine in Patients with Melanoma Treated with Immunotherapy.

The BNT162b2 vaccine against SARS-CoV-2 has a proven efficacy and a favorable safety profile. In cancer patients under immunotherapy in the form of immune-checkpoint inhibitors (ICIs), the efficacy of the vaccine has not been thoroughly studied, while a theoretical concern has also been raised about triggering immune-related adverse events (irAEs) by the vaccine. We conducted a prospective, non-interventional study on the immunogenicity and safety of the BNT162b2 vaccine in patients with advanced or metastatic melanoma treated with ICIs. Blood samples were obtained 0-4 days before the first dose and 12-21 days after the second dose of the vaccine for the quantification of the SARS-CoV-2 anti-spike antibody using an ELISA and immunophenotyping of the T and myeloid cell subpopulations. The active recording of AEs for a two-month period was conducted. Forty patients were included in the study. All but one (97.3%) achieved seroconversion after two doses of the vaccine and no correlations of the antibody titers with any of the studied parameters (age, gender, stage and duration of the disease, type of ICI, previous treatment, etc.) were found. Moreover, no differences in the subpopulations of the T cells (including the T-regulatory cells) or the myeloid cells were found pre- and post-vaccination. All AEs were low-grade, while one case of arthritis exacerbation was noted. The seroconversion rate in the studied population was high and was comparable to that of healthy subjects, while no major safety issues were raised during the safety follow-up. Finally, no derangements in the subpopulations of T cells or myeloid cells were noted. This is the first study focusing on the immunogenicity, safety, and effect of anti-SARS-CoV-2 vaccines on the blood-cell immunophenotype status of patients with melanoma treated with ICIs.

BreakAlign: a Perl program to align chimaeric (split) genomic NGS reads and allow visual confirmation of novel retroviral integrations.

BACKGROUND: Retroviruses replicate by integrating a DNA copy into a host chromosome. Detecting novel retroviral integrations (ones not in the reference genome sequence of the host) from genomic NGS data is bioinformatically challenging and frequently produces many false positives. One common method of confirmation is visual inspection of an alignment of the chimaeric (split) reads that span a putative novel retroviral integration site. We perceived the need for a program that would facilitate this by producing a multiple alignment containing both the viral and host regions that flank an integration. RESULTS: BreakAlign is a Perl program that uses blastn to produce such a multiple alignment. In addition to the NGS dataset and a reference viral sequence, the program requires either (a) the ~ 500nt host genome sequence that spans the putative integration or (b) coordinates of this putative integration in an installed copy of the reference human genome (multiple integrations can be processed automatically). BreakAlign is freely available from https://github.com/marchiem/breakalign and is accompanied by example files allowing a test run. CONCLUSION: BreakAlign will confirm and facilitate characterisation of both (a) germline integrations of endogenous retroviruses and (b) somatic integrations of exogenous retroviruses such as HIV and HTLV. Although developed for use with genomic short-read NGS (second generation) data and retroviruses, it should also be useful for long-read (third generation) data and any mobile element with at least one conserved flanking region.

Systematic review with meta-analysis: COVID-19 outcomes in patients receiving anti-TNF treatments.

BACKGROUND: Accumulating evidence suggests a beneficial effective of tumour necrosis factor-alpha (TNF-α) inhibitors on the outcomes of COVID-19 disease, which, however is not validated by all studies. AIMS: To perform a systematic review and meta-analysis of existing reports to investigate the impact of anti-TNF treatments on the clinical outcomes of COVID-19 patients. METHODS: A systematic search at PubMed and SCOPUS databases using specific keywords was performed. All reports of COVID-19 outcomes for patients receiving anti-TNF therapy by September 2021 were included. Pooled effect measures were calculated using a random-effects model. The Newcastle Ottawa Scale for observational studies was used to assess bias. Studies that were not eligible for meta-analysis were described qualitatively. RESULTS: In total, 84 studies were included in the systematic review, and 35 were included in the meta-analysis. Patients receiving anti-TNF treatment, compared to non-anti-TNF, among COVID-19 cases had a lower probability of hospitalisation (eight studies, 2555 patients, pooled OR = 0.53, 95% CI: 0.42-0.67, I2  = 0) and severe disease defined as intensive care unit admission or death (two studies, 1823 patients, pooled OR = 0.63, 95% CI: 0.41-0.96, I2  = 0), after adjustment for validated predictors of adverse disease outcomes. No difference was found for the risk for hospitalisation due to COVID-19 in populations without COVID-19 for patients receiving anti-TNF treatment compared to non-anti-TNF (three studies, 5 994 958 participants, pooled risk ratio = 0.97, 95% CI: 0.68-1.39, I2  = 20) adjusted for age, sex and comorbidities. CONCLUSIONS: TNF-α inhibitors are associated with a lower probability of hospitalisation and severe COVID-19 when compared to any other treatment for an underlying inflammatory disease.

Viral Causality of Human Cancer and Potential Roles of Human Endogenous Retroviruses in the Multi-Omics Era: An Evolutionary Epidemiology Review.

Being responsible for almost 12% of cancers worldwide, viruses are among the oldest known and most prevalent oncogenic agents. The quality of the evidence for the in vivo tumorigenic potential of microorganisms varies, thus accordingly, viruses were classified in 4 evidence-based categories by the International Agency for Research on Cancer in 2009. Since then, our understanding of the role of viruses in cancer has significantly improved, firstly due to the emergence of high throughput sequencing technologies that allowed the "brute-force" recovery of unknown viral genomes. At the same time, multi-omics approaches unravelled novel virus-host interactions in stem-cell biology. We now know that viral elements, either exogenous or endogenous, have multiple sometimes conflicting roles in human pathophysiology and the development of cancer. Here we integrate emerging evidence on viral causality in human cancer from basic mechanisms to clinical studies. We analyze viral tumorigenesis under the scope of deep-in-time human-virus evolutionary relationships and critically comment on the evidence through the eyes of clinical epidemiology, firstly by reviewing recognized oncoviruses and their mechanisms of inducing tumorigenesis, and then by examining the potential role of integrated viruses in our genome in the process of carcinogenesis.

Upregulation of Human Endogenous Retroviruses in Bronchoalveolar Lavage Fluid of COVID-19 Patients.

Severe COVID-19 pneumonia has been associated with the development of intense inflammatory responses during the course of infections with SARS-CoV-2. Given that human endogenous retroviruses (HERVs) are known to be activated during and participate in inflammatory processes, we examined whether HERV dysregulation signatures are present in COVID-19 patients. By comparing transcriptomes of bronchoalveolar lavage fluid (BALF) of COVID-19 patients and healthy controls, and peripheral blood monocytes (PBMCs) from patients and controls, we have shown that HERVs are intensely dysregulated in BALF of COVID-19 patients compared to those in BALF of healthy control patients but not in PBMCs. In particular, upregulation in the expression of specific HERV families was detected in BALF samples of COVID-19 patients, with HERV-FRD being the most highly upregulated family among the families analyzed. In addition, we compared the expression of HERVs in human bronchial epithelial cells (HBECs) without and after senescence induction in an oncogene-induced senescence model in order to quantitatively measure changes in the expression of HERVs in bronchial cells during the process of cellular senescence. This apparent difference of HERV dysregulation between PBMCs and BALF warrants further studies in the involvement of HERVs in inflammatory pathogenetic mechanisms as well as exploration of HERVs as potential biomarkers for disease progression. Furthermore, the increase in the expression of HERVs in senescent HBECs in comparison to that in noninduced HBECs provides a potential link for increased COVID-19 severity and mortality in aged populations. IMPORTANCE SARS-CoV-2 emerged in late 2019 in China, causing a global pandemic. Severe COVID-19 is characterized by intensive inflammatory responses, and older age is an important risk factor for unfavorable outcomes. HERVs are remnants of ancient infections whose expression is upregulated in multiple conditions, including cancer and inflammation, and their expression is increased with increasing age. The significance of this work is that we were able to recognize dysregulated expression of endogenous retroviral elements in BALF samples but not in PBMCs of COVID-19 patients. At the same time, we were able to identify upregulated expression of multiple HERV families in senescence-induced HBECs in comparison to that in noninduced HBECs, a fact that could possibly explain the differences in disease severity among age groups. These results indicate that HERV expression might play a pathophysiological role in local inflammatory pathways in lungs afflicted by SARS-CoV-2 and their expression could be a potential therapeutic target.

Validation of molecular clock inferred HIV infection ages: Evidence for accurate estimation of infection dates.

BACKGROUND: Improving HIV diagnosis, access to care and effective antiretroviral treatment provides our global strategy to reduce HIV incidence. To reach this goal we need to increase our knowledge about local epidemics. HIV infection dates would be an important information towards this goal, but they are largely unknown. To date, methods to estimate the dates of HIV infection are based mainly on laboratory or molecular methods. Our aim was to validate molecular clock inferred infection dates that were estimated by analysing sequences from 145 people living with HIV (PLHIV) with known transmission dates (clinically estimated infection dates). METHODS: All HIV sequences were obtained by Sanger sequencing and were previously found to belong to well-established molecular transmission clusters (MTCs). RESULTS: Our analysis showed that the molecular clock inferred infection dates were correlated with the clinically estimated ones (Spearman's Correlation coefficient = 0.93, p < 0.001) and that there was an agreement between them (Lin's concordance correlation coefficient = 0.92, p < 0.001). For the 61.4% of cases the molecular clock inferred preceded the clinically estimated infection dates. The median difference between clinically and molecularly estimated dates of infection was of 0.18 (IQR: -0.21, 0.89) years. The lowest differences were identified in people who inject drugs of our study population. CONCLUSIONS: The estimated time to more recent common ancestor (tMRCA) of nodes within clusters provides a reliable approximation of HIV infections for PLHIV infected within MTCs. Next-generation sequencing data and molecular clock estimates based on heterochronous sequences provide, probably, more reliable methods for inferring infection dates. However, since these data are not available in most of the HIV clinical laboratories, our approach, under specific conditions, can provide a reliable estimation of HIV infection dates and can be used for HIV public health interventions.

Human Endogenous Retrovirus-K HML-2 integration within RASGRF2 is associated with intravenous drug abuse and modulates transcription in a cell-line model.

HERV-K HML-2 (HK2) has been proliferating in the germ line of humans at least as recently as 250,000 years ago, with some integrations that remain polymorphic in the modern human population. One of the solitary HK2 LTR polymorphic integrations lies between exons 17 and 18 of RASGRF2, a gene that affects dopaminergic activity and is thus related to addiction. Here we show that this antisense HK2 integration (namely RASGRF2-int) is found more frequently in persons who inject drugs compared with the general population. In a Greek HIV-1-positive population (n = 202), we found RASGRF2-int 2.5 times (14 versus 6%) more frequently in patients infected through i.v. drug use compared with other transmission route controls (P = 0.03). Independently, in a United Kingdom-based hepatitis C virus-positive population (n = 184), we found RASGRF2-int 3.6 times (34 versus 9.5%) more frequently in patients infected during chronic drug abuse compared with controls (P < 0.001). We then tested whether RASGRF2-int could be mechanistically responsible for this association by modulating transcription of RASGRF2 We show that the CRISPR/Cas9-mediated insertion of HK2 in HEK293 cells in the exact RASGRF2 intronic position found in the population resulted in significant transcriptional and phenotypic changes. We also explored mechanistic features of other intronic HK2 integrations and show that HK2 LTRs can be responsible for generation of cis-natural antisense transcripts, which could interfere with the transcription of nearby genes. Our findings suggest that RASGRF2-int is a strong candidate for dopaminergic manipulation, and emphasize the importance of accurate mapping of neglected HERV polymorphisms in human genomic studies.

Impact of Interferon-α Receptor-1 Promoter Polymorphisms on the Transcriptome of the Hepatitis B Virus-Associated Hepatocellular Carcinoma.

Background and aims: Genetic polymorphisms within the promoter of interferon-α receptor type-1 (IFNAR1) have been associated with the susceptibility to and the outcome of chronic hepatitis B virus (HBV) infection. However, the impact of these polymorphisms in the transcriptome of the HBV-associated hepatocellular carcinoma (HCC) remains largely unexplored. Methods: Using whole-genome and exome sequencing data from The Cancer Genome Atlas project, we characterized three single-nucleotide polymorphisms (SNPs: -568G/C, -408C/T, -3C/T) and one variable number tandem repeat [VNTR: -77(GT)n] within the IFNAR1 promoter sequence in 49 HCC patients. RNAseq data from 10 genotyped HCC samples were grouped according to their -77VNTR or -3SNP genotype to evaluate the impact of these polymorphisms on the differential expression on the HCC transcriptome. Results: There is a fourfold higher impact of the -77VNTR on the HCC transcriptome compared to the -3SNP (q < 0.1, p < 0.001). The expression of the primary IFNAR1 transcript is not affected by these polymorphisms but a secondary, HCC-specific transcript is expressed only in homozygous -77VNTR ≤8/≤8(GT)n samples (p < 0.05). At the same time, patients carrying at least one -77VNTR >8(GT) allele, presented a strong upregulation of the fibronectin-1 (FN-1) gene, which has been associated with the development of HCC. Gene Ontology and pathway enrichment analysis of the differentially expressed genes revealed a strong disruption of the PI3K-AKT signaling pathway, which can be partially triggered by the extracellular matrix FN-1. Conclusion: The IFNAR-1 promoter polymorphisms are not involved in the expression levels of the main IFNAR-1 transcript. The -77VNTR has a regulatory role on the expression of a secondary, truncated, HCC-specific transcript, which in turn coincides with disruptions in cancer-associated pathways and in FN-1 expression modifications.

Roles of Endogenous Retroviruses in Early Life Events.

A retrovirus that infected our ancestors 100 million years ago became a human gene that is expressed in embryos and cancers, and can be detected in the blood of pregnant women. Accumulating evidence suggests potential roles for endogenous retroviruses in early life events, which may affect adult health.

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.

Reducing HIV infection in people who inject drugs is impossible without targeting recently-infected subjects.

OBJECTIVE: Although our understanding of viral transmission among people who inject drugs (PWID) has improved, we still know little about when and how many times each injector transmits HIV throughout the duration of infection. We describe HIV dynamics in PWID to evaluate which preventive strategies can be efficient. DESIGN: Due to the notably scarce interventions, HIV-1 spread explosively in Russia and Ukraine in 1990s. By studying this epidemic between 1995 and 2005, we characterized naturally occurring transmission dynamics of HIV among PWID. METHOD: We combined publicly available HIV pol and env sequences with prevalence estimates from Russia and Ukraine under an evolutionary epidemiology framework to characterize HIV transmissibility between PWID. We then constructed compartmental models to simulate HIV spread among PWID. RESULTS: In the absence of interventions, each injector transmits on average to 10 others. Half of the transmissions take place within 1 month after primary infection, suggesting that the epidemic will expand even after blocking all the post-first month transmissions. Primary prevention can realistically target the first month of infection, and we show that it is very efficient to control the spread of HIV-1 in PWID. Treating acutely infected on top of primary prevention is notably effective. CONCLUSION: As a large proportion of transmissions among PWID occur within 1 month after infection, reducing and delaying transmissions through scale-up of harm reduction programmes should always form the backbone of HIV control strategies in PWID. Growing PWID populations in the developing world, where primary prevention is scarce, constitutes a public health time bomb.

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.

Human endogenous retrovirus (HERV) expression is not induced by treatment with the histone deacetylase (HDAC) inhibitors in cellular models of HIV-1 latency.

BACKGROUND: While antiretroviral therapies have improved life expectancy and reduced viral loads in HIV-1-positive individuals, the cessation of treatment results in a rebound of viral replication. This suggests that a reservoir of latently-infected cells remains within these patients, the identity of which is ill-defined and therefore difficult to target therapeutically. Current strategies are aimed at using drugs such as histone deacetylase (HDAC) inhibitors to induce the expression of latent HIV-1 proviruses in order to activate and ultimately eradicate this reservoir of infected cells. One concern with the use of HDAC inhibitors is that they could up-regulate human endogenous retroviruses (HERVs), as well as HIV-1, with potentially pathophysiological consequences. RESULTS: In this study, we analysed the transcription of HERV genes in HIV-1 latency T cell (J-LAT 8.4) and monocyte (U1) models following treatment with the HDAC inhibitors, vorinostat, panobinostat and romidepsin. We examined the expression of HERV-K (HML-2) env and pol, as well as the co-opted genes HERV-W env (syncytin-1), HERV-FRD env (syncytin-2), in these cell lines. Finally, we investigated HERV expression in primary human T cells. CONCLUSIONS: We show that HDAC inhibitors did not substantially increase the transcription of the analysed HERV env or pol genes, suggesting that histone acetylation is not crucial for controlling HERV expression in these experimental models and in ex vivo primary human T cells. Importantly, this indicates that unwanted HERV expression does not appear to be a barrier to the use of HDAC inhibitors in HIV-1 cure strategies.

Larger mammalian body size leads to lower retroviral activity.

Retroviruses have been infecting mammals for at least 100 million years, leaving descendants in host genomes known as endogenous retroviruses (ERVs). The abundance of ERVs is partly determined by their mode of replication, but it has also been suggested that host life history traits could enhance or suppress their activity. We show that larger bodied species have lower levels of ERV activity by reconstructing the rate of ERV integration across 38 mammalian species. Body size explains 37% of the variance in ERV integration rate over the last 10 million years, controlling for the effect of confounding due to other life history traits. Furthermore, 68% of the variance in the mean age of ERVs per genome can also be explained by body size. These results indicate that body size limits the number of recently replicating ERVs due to their detrimental effects on their host. To comprehend the possible mechanistic links between body size and ERV integration we built a mathematical model, which shows that ERV abundance is favored by lower body size and higher horizontal transmission rates. We argue that because retroviral integration is tumorigenic, the negative correlation between body size and ERV numbers results from the necessity to reduce the risk of cancer, under the assumption that this risk scales positively with body size. Our model also fits the empirical observation that the lifetime risk of cancer is relatively invariant among mammals regardless of their body size, known as Peto's paradox, and indicates that larger bodied mammals may have evolved mechanisms to limit ERV activity.

Unfixed endogenous retroviral insertions in the human population.

UNLABELLED: 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.