PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2011; 6(5): e20234.
Published online May 25, 2011. doi:  10.1371/journal.pone.0020234
PMCID: PMC3102101
Human Endogenous Retrovirus K106 (HERV-K106) Was Infectious after the Emergence of Anatomically Modern Humans
Aashish R. Jha,1* Douglas F. Nixon,2 Michael G. Rosenberg,3 Jeffrey N. Martin,4 Steven G. Deeks,4 Richard R. Hudson,5 Keith E. Garrison,#6 and Satish K. Pillai#7*
1Department of Human Genetics, University of Chicago, Chicago, Illinois, United States of America
2Division of Experimental Medicine, University of California San Francisco, San Francisco, California, United States of America
3Albert Einstein College of Medicine, Bronx, New York, United States of America
4HIV/AIDS Division, Department of Medicine, San Francisco General Hospital, University of California San Francisco, San Francisco, California, United States of America
5Department of Evolution and Ecology, The University of Chicago, Chicago, Illinois, United States of America
6Department of Biology, Saint Mary's College of California, Moraga, California, United States of America
7Division of Infectious Diseases, Department of Medicine, San Francisco Veterans Affairs Medical Center, University of California San Francisco, San Francisco, California, United States of America
Michael Hendricks, Editor
Harvard University, United States of America
#Contributed equally.
* E-mail: nepaliaashish/at/gmail.com (ARJ); satish.pillai/at/ucsf.edu (SKP)
Conceived and designed the experiments: ARJ SKP KEG DFN. Performed the experiments: ARJ. Analyzed the data: ARJ SKP RRH. Contributed reagents/materials/analysis tools: RRH KEG DFN MGR JNM SGD. Wrote the paper: ARJ DFN MGR SGD KEG SKP.
Received January 25, 2011; Accepted April 19, 2011.
HERV-K113 and HERV-K115 have been considered to be among the youngest HERVs because they are the only known full-length proviruses that are insertionally polymorphic and maintain the open reading frames of their coding genes. However, recent data suggest that HERV-K113 is at least 800,000 years old, and HERV-K115 even older. A systematic study of HERV-K HML2 members to identify HERVs that may have infected the human genome in the more recent evolutionary past is lacking. Therefore, we sought to determine how recently HERVs were exogenous and infectious by examining sequence variation in the long terminal repeat (LTR) regions of all full-length HERV-K loci. We used the traditional method of inter-LTR comparison to analyze all full length HERV-Ks and determined that two insertions, HERV-K106 and HERV-K116 have no differences between their 5′ and 3′ LTR sequences, suggesting that these insertions were endogenized in the recent evolutionary past. Among these insertions with no sequence differences between their LTR regions, HERV-K106 had the most intact viral sequence structure. Coalescent analysis of HERV-K106 3′ LTR sequences representing 51 ethnically diverse individuals suggests that HERV-K106 integrated into the human germ line approximately 150,000 years ago, after the emergence of anatomically modern humans.
Endogenous retroviruses (ERVs) are the fossilized, germ line-integrated remnants of exogenous retroviruses and comprise approximately 8% of the human genome [1], [2]. Upon integration into germ cells, ERVs are transmitted in a Mendelian fashion, accumulating mutations over time that erode and ultimately eliminate their original viral sequence structure. Just as ERV genomes are subjected to modification within the host, host-ERV interactions may also induce changes in the host genome. For example, homologous recombination between distant HERV loci have likely induced chromosomal rearrangements in the human host, and coevolutionary conflicts with the exogenous precursors of HERVs may have resulted in sequence modifications in human antiretroviral genes such as APOBEC3H and TRIM5α [3], [4], [5]. Thus, identifying HERVs that integrated into the human genome in the recent evolutionary past may allow us to catalog proviruses that influenced the evolution of genes involved in antiviral responses. In addition, HERVs have been implicated in playing a role in a number of diseases [6], [7], [8], [9], [10], [11]. It has been suggested that cross-talk between HIV-1 and HERV transactivation pathways may compromise HIV-1 eradication strategies aimed at purging latently-infected cellular reservoirs [12]. Thus, detailed investigation of the genomic architecture and evolutionary history of the most recent HERV insertions may provide vital insights into viral inactivation and the present-day conflict between humans and infectious retroviruses such as HIV-1 and HTLV-1.
The most ancient HERVs are 60–70 million years old [1], [13]. However, some members of the HERV-K family integrated into the human genome after the hominid-chimpanzee split 6 million years ago (Mya) [14], [15], [16], [17], and these human-specific insertions are considered the youngest of all ERVs in the human genome. In this study, we analyzed all full-length HERV-K (HML-2) loci in the human genome to determine how recently in human evolutionary history the infectious progenitors of endogenous retroviruses were active.
The human genome consists of numerous HERV insertions that contain 5′ and 3′ LTRs and varying remnants of the viral coding genes in between. Many more insertions currently exist as solo-LTRs, which are the remnants of full-length insertions that have been truncated into a single LTR by virtue of inter-LTR recombination events. Inter-LTR comparison is the standard method used to estimate the age of full-length HERV insertions. At the time of insertion, the LTR sequences are presumed to be identical and the differences between the 5′ and 3′ LTRs of the same provirus are believed to arise due to substitutions accumulating in a clocklike manner post-insertion [18]. To elucidate the insertion timeline of the human-specific endogenous retroviruses into the human genome, we obtained the complete genome sequences of all the human-specific full-length HERV-K insertion previously reported [15], [16] using the UCSC genome browser (hg18).
It has previously been hypothesized that gene conversion may be rampant amongst HERV loci [19]. Pervasive gene conversion would prevent us from applying a molecular clock to these data, since observed sequence variation between LTR regions would be attributed to recombination rather than the stepwise accumulation of mutations predicted by the neutral theory of evolution [20]. We reconstructed a maximum likelihood phylogeny of all HERV-K LTRs to determine the prevalence of gene conversion between insertions in the HERV-K family (Figure 1). The clustering of the 5′ and 3′ LTR sequences associated with each insertion locus in our phylogeny suggests that gene conversion is rare among HERV-Ks, with the single exception of HERV-K115. Evidence of a gene conversion event in HERV-K115 has been reported previously [17], [19].
Figure 1
Figure 1
Phylogenetic tree of full-length HERV-K (HML-2) LTR sequences.
To estimate the insertion times of each human-specific endogenous retrovirus, we calculated the age of each human-specific full-length HERV-K insertion using the traditional inter-LTR comparison method [18]. Inter-LTR divergence measurements of the fourteen insertions that did not undergo gene conversion in the LTR region were converted to insertion age estimates by applying an established HERV-K LTR-specific divergence rate of 0.13% per million years (Myr) to these data [21] (Table 1). Based on this method, HERV-K1p31.1 (referred to as HERV-K116 from here onwards) and HERV-K106 have the highest probability of being the youngest full-length endogenous retroviruses in the human genome. Due to the perfect identity between their 5′ and 3′ LTR sequences, an inter-LTR comparison is only informative to calculate a maximum, upper-bound age estimate. A single mutation is predicted to arise within an LTR every 0.8 Myr based on the divergence rate of 0.13%/Myr. Having no sequence differences between their LTRs, inter-LTR analysis of HERV-K116 and HERV-K106 suggests that both insertions must be younger in age than the time required for one sequence polymorphism to arise between the two LTRs, or 0.8 Myr.
Table 1
Table 1
Human specific complete HERV-K (HML-2) proviruses within the human genome.
We used Multalin [22] to align HERV-K113, K115, K106, and K116 with experimentally reconstituted infectious HERVs KCON [23] and Phoenix [24] that are based on artificial consensus sequences of full-length human-specific HERV-K (HML-2) insertions. We observed that while all four HERV insertions (K113, K115, K106, and K116) exhibited similarities to the reconstituted viruses, all four contained mutations that were unique to each insertion (Figure S1). We observed that both HERV-K106 and HERV-K116 are members of the type I HERV-K family as evidenced by the presence of a 292 bp ‘deletion’ in env which is the signature of all type I HERV members. The presence of this 292 bp env deletion in multiple HERV-K type I members suggests that this deletion may have been present in the infectious ancestral precursors of these viruses and probably does not render a HERV insertion dysfunctional on its own. However, HERV-K116 also has a 2846 bp deletion in its pol sequence [16]. In contrast, HERV-K106 exhibits relatively intact retroviral genome architecture (Figure 2A). Thus, these data suggest that HERV-K106 is the youngest endogenous retrovirus that survives largely intact in the human genome.
Figure 2
Figure 2
Genome organization and haplotypes of HERV-K106.
While inter-LTR comparison is a useful tool to estimate the age of HERV insertions, it becomes less informative with fewer mutations between the LTR regions and provides only an upper bound age estimate in the absence of sequence differences between the LTR regions. We recently developed an alternative dating method to infer insertion age when the inter-LTR method is inapplicable (to estimate the insertion dates of solo-LTRs and HERV loci such as K115 that show evidence of gene conversion) [25]. HERV insertions with identical 5′ and 3′ LTR sequences represent another scenario in which our alternative method is useful. We applied this method to the HERV-K106 insertion to derive a more precise estimate of its age. Our approach involves the application of coalescent inference to inter-host sequence variation in one of the proviral LTR sequences to estimate HERV insertion age. We generated complete HERV-K106 3′LTR sequences from 51 individuals representing various ethnicities and three different geographical locations within the United States (Table 2). PCR amplification and sequencing of the HERV-K106 LTR revealed three single nucleotide polymorphic positions (SNPs): 133, 403, and 835 (numbered according to their position in the GenBank reference sequence AF164620) (Table 2). Four HERV-K106 haplotypes could be constructed based on the SNPs identified in the 3′LTR region (Figure 2B).
Table 2
Table 2
Base frequencies and haplotypes of HERV-K106 3'LTR with haplotype frequencies in various ethnic groups within the United States.
The coalescent estimation of insertion age rests on the assumption that the insertion site is evolving neutrally. Therefore, we conducted tests on the HERV-K106 insertion site to determine whether the assumption of neutrality was maintained. We performed 10,000 coalescent simulations using MS software [26] and Schaffner's calibrated model of human genome evolution [27] to calculate the probability of observing exactly 3 SNPs in a neutrally evolving, 960 bp stretch of human DNA. We assumed no recombination, and an inferred human mutation rate of 9.0×10−9 subsitutions per site per generation. These simulations yielded a Gaussian mutational probability distribution with a mean of 8 SNPs, and revealed that the presence of only 3 SNPs in a 960 bp neutrally evolving region of the human genome deviates significantly from expectations based on Schaffner's model (p<0.05). These findings suggest one of two possibilities about HERV-K106. Either it is evolving under selection, or it is evolving neutrally but has not evolved in tandem with the human genome for a sufficient amount of time to conform to the predictions of Schaffner's human genome-based model. We explored the selection hypothesis by performing standard tests of neutrality on the K106 locus, including Tajima's D [28], and Fu and Li's D* and H [29] (performed using DNASP [30]). In all cases, we could not reject the null hypothesis that HERV-K106 is evolving neutrally (p>0.10). HERV-K106 itself may have been evolving neutrally but could have been driven to fixation due to hitchhiking effects if the region flanking the insertion had been under positive selection. We used the HGDP Selection Browser [31] to test whether the genomic region containing the HERV-K106 insertion site is under selection by calculating the iHS and XP-EHH statistics on genotypes in the Human Genome Diversity Panel. The iHS and XP-EHH statistics are haplotype homozygosity-based tests used to detect signatures of recent selection on variants that have not yet reached fixation [32] and can be applied to detect selective sweeps in alleles that have approached fixation in one population but are polymorphic in the overall human population [33]. Even though HERV-K106 is fixed in all humans, the genetic region flanking the K106 insertion may contain SNPs that could reveal if a selective sweep has occurred in this region. We found that HERV-K106 is incorporated into a genomic location with only a few genes nearby (Figure S2) and the region in chromosome 3 containing HERV-K106 exhibited no signatures of selection, as both iHS and XP-EHH did not yield extreme values (Figure S3). These data collectively support that K106 is evolving neutrally and has only shared its evolutionary history with the human genome for a relatively short period of time.
We constructed a maximum likelihood phylogeny of all 94 observed HERV-K106 3′ LTR haplotype sequences to estimate the age of the K106 insertion (Figure S4). According to coalescent theory, the genetic distance to the inferred most recent common ancestor (MRCA) should reflect the time that has elapsed since the establishment of the ancestral sequence, and in this particular case, the age of the proviral insertion itself. We used two previously reported evolutionary rates to translate genetic distance into coalescence time. The upper-bound for the coalescence-based age estimate was inferred using the HERV-K LTR specific mutation rate of 1.3×10−9 mutations/site/year [21], and the inferred mammalian genome mutation rate of 2.2×10−9 mutations/site/year [34] was used to calculate a lower-bound estimate. Based on the two divergence rates, we estimate that HERV-K106 was integrated into the human genome between 91,000 and 154,000 years ago, after the emergence of anatomically modern humans [35], [36].
Our study suggests that HERV-K106 is a recent acquisition in the human genome and is possibly specific to Homo sapiens sapiens. HERV-K106 is fixed in humans, but HERV-K113 and HERV-K115 remain insertionally polymorphic [17], despite being inserted into the human genome several hundreds of thousands of years before K106 [25]. Features associated with the three regions of the human genome harboring these HERV insertions may help to explain the apparent discordance between insertion age and fixation. The fixation probability of an ERV may be inversely correlated with local chromosomal recombination rate and local gene density [37], [38]. We employed the high resolution recombination map of Kong et al. [39] to examine local recombination rates in the genomic neighborhoods of HERV-K106, K113 and K115 insertions. Analysis of the map revealed that HERV-K106 lies in one of only nineteen recognized “recombinaton deserts” in the human genome (crossover rates less than 0.3 c/Mb), which may have accelerated its fixation. Conversely, K113 and K115 reside in regions that are considerably more recombinogenic (1.18 cM/Mb and 2.14 cM/Mb, respectively) [39]. The low rate of recombination in the region surrounding K106 may have accelerated its fixation, or the high rates of recombination in the regions surrounding K113 and K115 may have decelerated the fixation of these older insertions. We next employed the NCBI Human Genome Map Viewer [40] to investigate the potential role of gene density in the neighborhood of these three HERV insertions, since a highly significant inverse correlation between local gene density and HERV fixation rate has been demonstrated [37]. The insertionally polymorphic HERV-K113 resides in chromosome 19, which has the highest gene density of all human chromosomes. We surveyed regions within a window of 1 Mb upstream to 1 Mb downstream of each HERV insertion for the presence of characterized genes to obtain a more localized estimate of gene density. HERV-K106, K113 and K115 are flanked by 32, 64, and 115 known genes, respectively. Again, the low gene density in the region surrounding K106 may have accelerated its fixation, or the higher gene densities in the regions surrounding K113 and K115 may have hindered the fixation of these older insertions. Taken together, the patterns of local recombination rate and gene density support genomic context as a primary determinant of fixation that explains the discordance between insertion age and insertional polymorphism in the human population observed within the HERV-K family.
There are many fixed and insertionally polymorphic HERV-K(HML-2) solo-LTRs in the human genome [16], [41], [42] and it is possible that some of these unfixed solo-LTRs may have been derived from proviral insertions that are younger than HERV-K106. However, the absence of a second LTR from these insertions prevents us from including them in our current investigation because we relied on the inter-LTR comparisons to identify HERVs that recently integrated in the human genome. A second caveat stems from the reliance of our age estimates on previously established mutation rates; changes in local mutation or recombination rate may influence the apparent age of this insertion. However, our analyses of the local recombination rates surrounding the K106 and K113 insertions suggest that this is not likely to be a confounding factor. Furthermore, we previously hervotyped a small subset of individuals in our samples (n = 16) to study haplotype diversity in HERV-K113 [25]. Haplotype diversity for K113 in this smaller subset of samples was higher than the haplotype diversity that we observed for HERV-K106 in the same subset, as would be expected when comparing an older to a younger insertion. The frequencies of the K113 minor haplotypes in individuals of both African and European descent were much higher than the frequencies of minor haplotypes that we observed for HERV-K106 (Figure 3).
Figure 3
Figure 3
Comparison of haplotype frequencies between HERV K106 and HERV-K113.
To the best of our knowledge, this is the first rigorous evidence that the infectious progenitors of human endogenous retroviruses were active after the origin of anatomically modern man, Homo sapiens, in Africa [35], [36].
Samples and DNA extraction
Peripheral blood samples were obtained with written consent from a total of 51 individuals drawn from the pediatric HIV cohort at Jacobi Medical Center (Bronx, NY), the SCOPE cohort of chronic HIV-1 infection (San Francisco, CA), and from blood donors to the Stanford Blood Bank (Palo Alto, CA). Twelve participants were perinatally HIV-1 infected children followed at Jacobi Medical Center, Bronx, NY. Twenty-five samples of HIV-1 infected adults were obtained from the SCOPE cohort in San Francisco, California and fourteen individuals from Stanford Blood Bank were anonymous, healthy HIV-1 negative individuals. The study was approved by the Institutional Review Boards at Jacobi Medical Center and the University of California at San Francisco (UCSF). All participants self-identified their ethnicities. HERV-K106 hervotyping [25] and haplotype construction for this study was carried out without prior knowledge of the ethnicity of the participants. Genomic DNA was isolated according to the manufacturer's protocol using QIAmp DNA Blood Mini Kit (Qiagen, Valencia, California, USA).
PCR and gel electrophoresis
Insertion sites, solo-LTR, and full length HERV-K106 were detected using the HERV-K106 specific primers 5′-GGTGTTGCTGTGGAAGGTATTC-3′ and 5′-TCCATGGCTATCCACGAGA-3′ and the two step PCR protocol previously described [15], [16]. Neither insertional polymorphism nor solo-LTR formation was detected in our sample. We designed two primers HERV-K1063LTR1 5′-ATTTGGTGCCAGGAACTGAG-3′ and HERV-K1063LTR2 5′-AAGAAAAGGGGGAAATGTGG-3′ that were specific to the HERV-K106 3′LTR for sequencing purposes. HERV-K1063LTR1 annealed to the host flanking DNA upstream of HERV-K106 insertion in chromosome 3 and HERV-K1063LTR2 bound downstream of the 3′LTR to the env gene of HERV-K106. PCR was performed in a volume of 50 uL with 25 ng of genomic DNA in 1.0 uM of each primer, 200 uM of each dNTP, 25 mM of MgCl2, and 0.5 units of Taq DNA polymerase obtained from Applied Biosystems, Foster City, California, USA. PCR was conducted with 10 minutes of initial denaturing at 95°C, 30 cycles of 10 seconds at 94°C, 30 seconds at annealing temperature (61.5°C), and 59 seconds of elongation at 72°C followed by a 10 min final extension at 72°C. Amplified PCR products were stored at 4°C. The amplified PCR products were detected as ~1200 bp bands in 2% agarose gel. A low-mass DNA ladder (Invitrogen, Carlsbad, California, USA) was used to verify the product size (Figure S5).
Sequencing and sequence analysis
Sequencing was performed at MCLAB, South San Francisco using an ABI 3730XL sequencer. To determine an individual's hervotype (HERV genotype) for HERV-K106, sequences were aligned and single nucleotide polymorphisms (SNPs) were identified using Sequencher version 4.9 [25]. To rule out errors introduced during the PCR or sequencing process, PCR was performed twice on each genomic DNA sample and SNPs were verified by bidirectional sequencing of the amplicons produced by the two independent PCR on the genomic DNA. There were ten individuals who were heterozygous at the SNP sites. Six were heterozygous at only one site and the haplotypes of these six individuals could be inferred. The remaining four individuals were heterozygous at more than one position and were excluded from the analysis since their haplotypes could not be reconstructed. After excluding the ambiguous heterozygotes, we analyzed the base frequencies at the three SNP sites of the HERV-K106 3′LTR in 102 HERV-K106 sequences. The base frequencies and SNP positions are summarized in table 2. We also compared HERV-K106 diversity between various ethnicities in the USA represented in our sample. While African Americans (n = 13) and European Americans (n = 22) were abundant in our sample, other ethnic groups such as Asians, Hispanics, East Indians, and Native Americans were represented in very low numbers (n<5) in our sample. Thus, we compared HERV-K106 haplotype frequencies in African Americans and European Americans and the rest of the ethnicities were grouped together as ‘others’.
Coalescence and phylogenetic inference
Coalescence time was calculated as described previously [25]. The sequence of the most recent common ancestor (MRCA) of the HERV-K106 insertion was inferred using maximum likelihood inference and K110 (GenBank Accession number AF164617) as an outgroup taxon (Figure S4). Coalescence (insertion) times were approximated using the formula An external file that holds a picture, illustration, etc.
Object name is pone.0020234.e001.jpg where T is the time (My) passed since the insertion event, n is the number of observed taxa, D is the divergence associated with a descendant taxon (cumulative branch lengths between observed haplotype and inferred MRCA), and R represents one of two previously reported divergence (evolutionary) rates (2.2×10−9 and 1.3×10−9 mutations/site/year) used to generate upper and lower bound estimates of insertion times [21], [34]. We used measurements of genetic divergence from inferred ancestral 5' LTR sequences (Figure S4) to estimate the HERV-K106 insertion time.
Figure S1
Comparison of HERV-K106 genome to K113, K115, K116, KCON and Phoenix. Multalin [22] was used to align K106 genome to that of two full length human specific HERV insertions (K113 and K115). We also included in the alignment two experimentally reconstituted HERVs KCON and Phoenix that are infectious and K116 that has a 2846 bp deletion in its pol gene. We observed that while all four HERV insertions (K113, K115, K106, and K116) exhibited similarities to the reconstituted viruses, all four contained mutations that were unique to each insertion. The identical regions in the genomes of all these HERV insertions are shown in blue. Genomic regions in red indicate regions that vary between at least one of these HERV insertions and ‘–’ indicates deleted regions.
(TIF)
Figure S2
Genomic location of HERV-K106. Top: High resolution genomic location of HERV-K106 showing insertion in a gene desert. No genes can be seen within ~1 kb upstream or downstream of HERV-K106. Bottom: In a lower resolution genomic location map of the K106 insertion, the immune genes such as CD200R1 and CD200RL1 are visible far upstream of HERV-K106.
(TIF)
Figure S3
Lack of selective sweep on or near HERV-K106 region. The fixation of HERV-K106 in the human population could be due to drift, selection on HERV-K106, or a hitchhiking effect driven by other gene(s). If HERV-K106 reached fixation because it is under selection or due to a hitchhiking effect (i.e., if selection drove a nearby mutation to fixation and that mutation happened to be on a haplotype containing HERV-K106), then the XP-EHH statistic which measures complete selective sweeps should indicate selection acting on or nearby HERV-K106. We searched for signals of selection as measured by iHS and XP-EHH statistics in the samples from Human Genome Diversity Panels using tools on Prof. Johnathan Pritchard's webpage (http://hgdp.uchicago.edu/). We were unable to detect any signals of selection on or near HERV-K106. We were able to detect signals of selection in regions 500 kb farther away from HERV-K106 insertion region, indicating that selection can be detected in that region but HERV-K106 is not under positive selection.
(TIF)
Figure S4
HERV-K106 3′LTR haplotypes. ML tree of four haplotypesof HERV-K106 3′LTR. The most prevalent haplotypewas CGC (green) which is also the haplotypeof the HERV-K106 3′LTR in Genbank(AF164620). Minor haplotypes CGT (red), TGC (blue), and TAC (orange) are also shown. All haplotypescluster together. The HERV-K110 3′LTR was used as the outgroup.
(TIF)
Figure S5
HERV-K106 3′LTR amplification. The two primers specific to HERV-K106 3′LTR that we designed (HERV-K1063LTR1 5′-ATTTGGTGCCAGGAACTGAG-3′ and HERV-K1063LTR2 5′-AAGAAAAGGGGGAAATGTGG-3′) were used to amplify the complete HERV-K106 3′ LTR as detected by ~1200 bp bands in 2% agarosegel. Low DNA Mass ladder (Invitrogen, Carlsbad, California, USA) was used to verify the product size.
(TIF)
Acknowledgments
We thank the academic editor Dr. Michael Hendricks and the anonymous reviewer for helpful comments.
Footnotes
Competing Interests: The authors have declared that no competing interests exist. Dr. Douglas Nixon is a section editor for PLoS ONE. This manuscript was handled by a different section editor.
Funding: This work was supported by National Institutes of Health grants AI076059, AI060379, AI084113, UCSF CFAR (P-30 AI27763), UCSF CTSI (UL1 RR024131), CNICS (R24 AI067039), NIAID (K24AI069994), and 1K01DA024654 grants. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
1. Bannert N, Kurth R. The evolutionary dynamics of human endogenous retroviral families. Annu Rev Genomics Hum Genet. 2006;7:149–173. [PubMed]
2. Paces J, Pavlicek A, Paces V. HERVd: database of human endogenous retroviruses. Nucleic Acids Res. 2002;30:205–206. [PMC free article] [PubMed]
3. Hughes JF, Coffin JM. Evidence for genomic rearrangements mediated by human endogenous retroviruses during primate evolution. Nat Genet. 2001;29:487–489. [PubMed]
4. OhAinle M, Kerns JA, Li MM, Malik HS, Emerman M. Antiretroelement activity of APOBEC3H was lost twice in recent human evolution. Cell Host Microbe. 2008;4:249–259. [PMC free article] [PubMed]
5. Kaiser SM, Malik HS, Emerman M. Restriction of an extinct retrovirus by the human TRIM5alpha antiviral protein. Science. 2007;316:1756–1758. [PubMed]
6. Christensen T. HERVs in Neuropathogenesis. J Neuroimmune Pharmacol. 2010. [PubMed]
7. Moyes DL, Martin A, Sawcer S, Temperton N, Worthington J, et al. The distribution of the endogenous retroviruses HERV-K113 and HERV-K115 in health and disease. Genomics. 2005;86:337–341. [PubMed]
8. Contreras-Galindo R, Kaplan MH, Leissner P, Verjat T, Ferlenghi I, et al. Human endogenous retrovirus K (HML-2) elements in the plasma of people with lymphoma and breast cancer. J Virol. 2008;82:9329–9336. [PMC free article] [PubMed]
9. Contreras-Galindo R, Lopez P, Velez R, Yamamura Y. HIV-1 infection increases the expression of human endogenous retroviruses type K (HERV-K) in vitro. AIDS Res Hum Retroviruses. 2007;23:116–122. [PubMed]
10. Garrison KE, Jones RB, Meiklejohn DA, Anwar N, Ndhlovu LC, et al. T cell responses to human endogenous retroviruses in HIV-1 infection. PLoS Pathog. 2007;3:e165. [PMC free article] [PubMed]
11. Moyes D, Griffiths DJ, Venables PJ. Insertional polymorphisms: a new lease of life for endogenous retroviruses in human disease. Trends Genet. 2007;23:326–333. [PubMed]
12. Trono D, Van Lint C, Rouzioux C, Verdin E, Barre-Sinoussi F, et al. HIV persistence and the prospect of long-term drug-free remissions for HIV-infected individuals. Science. 2010;329:174–180. [PubMed]
13. Tristem M. Identification and characterization of novel human endogenous retrovirus families by phylogenetic screening of the human genome mapping project database. J Virol. 2000;74:3715–3730. [PMC free article] [PubMed]
14. Glazko GV, Nei M. Estimation of divergence times for major lineages of primate species. Mol Biol Evol. 2003;20:424–434. [PubMed]
15. Barbulescu M, Turner G, Seaman MI, Deinard AS, Kidd KK, et al. Many human endogenous retrovirus K (HERV-K) proviruses are unique to humans. Curr Biol. 1999;9:861–868. [PubMed]
16. Macfarlane C, Simmonds P. Allelic variation of HERV-K(HML-2) endogenous retroviral elements in human populations. J Mol Evol. 2004;59:642–656. [PubMed]
17. Turner G, Barbulescu M, Su M, Jensen-Seaman MI, Kidd KK, et al. Insertional polymorphisms of full-length endogenous retroviruses in humans. Curr Biol. 2001;11:1531–1535. [PubMed]
18. Dangel AW, Baker BJ, Mendoza AR, Yu CY. Complement component C4 gene intron 9 as a phylogenetic marker for primates: long terminal repeats of the endogenous retrovirus ERV-K(C4) are a molecular clock of evolution. Immunogenetics. 1995;42:41–52. [PubMed]
19. Belshaw R, Pereira V, Katzourakis A, Talbot G, Paces J, et al. Long-term reinfection of the human genome by endogenous retroviruses. Proc Natl Acad Sci U S A. 2004;101:4894–4899. [PubMed]
20. Kimura M. DNA and the neutral theory. Philos Trans R Soc Lond B Biol Sci. 1986;312:343–354. [PubMed]
21. Lebedev YB, Belonovitch OS, Zybrova NV, Khil PP, Kurdyukov SG, et al. Differences in HERV-K LTR insertions in orthologous loci of humans and great apes. Gene. 2000;247:265–277. [PubMed]
22. Corpet F. Multiple sequence alignment with hierarchical clustering. Nucleic Acids Res. 1988;16:10881–10890. [PMC free article] [PubMed]
23. Lee YN, Bieniasz PD. Reconstitution of an infectious human endogenous retrovirus. PLoS Pathog. 2007;3:e10. [PMC free article] [PubMed]
24. Dewannieux M, Harper F, Richaud A, Letzelter C, Ribet D, et al. Identification of an infectious progenitor for the multiple-copy HERV-K human endogenous retroelements. Genome Res. 2006;16:1548–1556. [PubMed]
25. Jha AR, Pillai SK, York VA, Sharp ER, Storm EC, et al. Cross-sectional dating of novel haplotypes of HERV-K 113 and HERV-K 115 indicate these proviruses originated in Africa before Homo sapiens. Mol Biol Evol. 2009;26:2617–2626. [PMC free article] [PubMed]
26. Hudson RR. Generating samples under a Wright-Fisher neutral model of genetic variation. Bioinformatics. 2002;18:337–338. [PubMed]
27. Schaffner SF, Foo C, Gabriel S, Reich D, Daly MJ, et al. Calibrating a coalescent simulation of human genome sequence variation. Genome Res. 2005;15:1576–1583. [PubMed]
28. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–595. [PubMed]
29. Fu YX, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993;133:693–709. [PubMed]
30. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–1452. [PubMed]
31. Browser THS The HGDP Selection Browser.
32. Voight BF, Kudaravalli S, Wen X, Pritchard JK. A map of recent positive selection in the human genome. PLoS Biol. 2006;4:e72. [PMC free article] [PubMed]
33. Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, et al. Genome-wide detection and characterization of positive selection in human populations. Nature. 2007;449:913–918. [PMC free article] [PubMed]
34. Kumar S, Subramanian S. Mutation rates in mammalian genomes. Proc Natl Acad Sci U S A. 2002;99:803–808. [PubMed]
35. McDougall I, Brown FH, Fleagle JG. Stratigraphic placement and age of modern humans from Kibish, Ethiopia. Nature. 2005;433:733–736. [PubMed]
36. White TD, Asfaw B, DeGusta D, Gilbert H, Richards GD, et al. Pleistocene Homo sapiens from Middle Awash, Ethiopia. Nature. 2003;423:742–747. [PubMed]
37. Katzourakis A, Pereira V, Tristem M. Effects of recombination rate on human endogenous retrovirus fixation and persistence. J Virol. 2007;81:10712–10717. [PMC free article] [PubMed]
38. Dolgin ES, Charlesworth B. The effects of recombination rate on the distribution and abundance of transposable elements. Genetics. 2008;178:2169–2177. [PubMed]
39. Kong A, Gudbjartsson DF, Sainz J, Jonsdottir GM, Gudjonsson SA, et al. A high-resolution recombination map of the human genome. Nat Genet. 2002;31:241–247. [PubMed]
40. Viewer NHGM NCBI Human Genome Map Viewer.
41. Belshaw R, Dawson AL, Woolven-Allen J, Redding J, Burt A, et al. Genomewide screening reveals high levels of insertional polymorphism in the human endogenous retrovirus family HERV-K(HML2): implications for present-day activity. J Virol. 2005;79:12507–12514. [PMC free article] [PubMed]
42. Hughes JF, Coffin JM. Human endogenous retrovirus K solo-LTR formation and insertional polymorphisms: implications for human and viral evolution. Proc Natl Acad Sci U S A. 2004;101:1668–1672. [PubMed]
43. Heslin DJ, Murcia P, Arnaud F, Van Doorslaer K, Palmarini M, et al. A single amino acid substitution in a segment of the CA protein within Gag that has similarity to human immunodeficiency virus type 1 blocks infectivity of a human endogenous retrovirus K provirus in the human genome. J Virol. 2009;83:1105–1114. [PMC free article] [PubMed]
Articles from PLoS ONE are provided here courtesy of
Public Library of Science