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

BACKGROUND: Human endogenous retroviruses (HERVs) represent the inheritance of ancient germ-line cell infections by exogenous retroviruses and the subsequent transmission of the integrated proviruses to the descendants. ERVs have the same internal structure as exogenous retroviruses. While no replication-competent HERVs have been recognized, some retain up to three of four intact ORFs. HERVs have been classified before, with varying scope and depth, notably in the RepBase/RepeatMasker system. However, existing classifications are bewildering. There is a need for a systematic, unifying and simple classification. We strived for a classification which is traceable to previous classifications and which encompasses HERV variation within a limited number of clades. RESULTS: The human genome assembly GRCh 37/hg19 was analyzed with RetroTector, which primarily detects relatively complete Class I and II proviruses. A total of 3173 HERV sequences were identified. The structure of and relations between these proviruses was resolved through a multi-step classification procedure that involved a novel type of similarity image analysis ("Simage") which allowed discrimination of heterogeneous (noncanonical) from homogeneous (canonical) HERVs. Of the 3173 HERVs, 1214 were canonical and segregated into 39 canonical clades (groups), belonging to class I (Gamma- and Epsilon-like), II (Beta-like) and III (Spuma-like). The groups were chosen based on (1) sequence (nucleotide and Pol amino acid), similarity, (2) degree of fit to previously published clades, often from RepBase, and (3) taxonomic markers. The groups fell into 11 supergroups. The 1959 noncanonical HERVs contained 31 additional, less well-defined groups. Simage analysis revealed several types of mosaicism, notably recombination and secondary integration. By comparing flanking sequences, LTRs and completeness of gene structure, we deduced that some noncanonical HERVs proliferated after the recombination event. Groups were further divided into envelope subgroups (altogether 94) based on sequence similarity and characteristic "immunosuppressive domain" motifs. Intra and inter(super)group, as well as intraclass, recombination involving envelope genes ("env snatching") was a common event. LTR divergence indicated that HERV-K(HML2) and HERVFC had the most recent integrations, HERVL and HUERSP3 the oldest. CONCLUSIONS: A comprehensive HERV classification and characterization approach was undertaken. It should be applicable for classification of all ERVs. Recombination was common among HERV ancestors.

Development of real-time PCRs for detection and quantitation of human MMTV-like (HML) sequences HML expression in human tissues.

The human genome contains around 1000 betaretrovirus-like copies, human mouse mammary tumour virus (MMTV)-like (HML) groups 1-10, also referred to as human endogenous retrovirus "HERV-K". Despite many efforts, it is not established whether betaretroviruses, exo- or endogenous, are involved in the etiology of breast cancer, or other cancer diseases, in humans. Quantitative real-time PCR (QPCR) TaqMan-based assays for HML groups 1-7, targeting the conserved reverse transcriptase (RT) and integrase (IN) domains of the pol gene were designed. Plasmids containing the entire pol gene of HML1-7 were used as standards. The RT and IN based QPCRs could detect 10(0)-10(3) copies per PCR reaction of the plasmids. However, not all plasmids gave a signal in both RT and IN QPCRs, probably due to mismatches. Furthermore, RT and IN based HML6 specific QPCRs were developed. They were specific for amplification of transcripts for the whole HML6 group. The methods allow the monitoring in body fluids and tissues of expression of a wide range of betaretrovirus-like sequences. Betaretrovirus-like RNA was studied in normal human tissues and of HML6 in brains of multiple sclerosis (MS) patients. Brain, adrenal gland and testis had a high betaretrovirus-like expression. Multiple sclerosis plaques contained the same HML6 RNA concentration as control tissue. These assays are expected to enhance studies on involvement of betaretroviruses in physiology and disease.

The phylogeny of orthoretroviral long terminal repeats (LTRs).

LTRs are sequence elements in retroviruses and retrotransposons which are difficult to align due to their variability. One way of handling such cases is to use Hidden Markov Models (HMMs). In this work HMMs of LTRs were constructed for three groups of orthoretroviruses: the betaretroviruslike human MMTV-like (HML) endogenous retroviruses, the lentiviruses, including HIV, and gammaretroviruslike human endogenous retroviruses (HERVs). The HMM-generated LTR alignments and the phylogenetic trees constructed from them were compared with trees based on alignments of the pol gene at the nucleic acid level. The majority of branches in the LTR and pol based trees had the same order for the three retroviral genera, showing that HMM methods are successful in aligning and constructing phylogenies of LTRs. The HML LTR tree deviated somewhat from the pol tree for the groups HML3, HML7 and HML6. Among the gammaretroviruslike proviruses, the exogenous Mouse Leukemia Virus (MLV) was highly related to HERV-T in the pol based tree, but not in the LTR based tree. Aside from these differences, the similarity between the trees indicates that LTRs and pol coevolved in a largely monophyletic way.

Classification and nomenclature of endogenous retroviral sequences (ERVs): problems and recommendations.

The genomes of many species are crowded with repetitive mobile sequences. In the case of endogenous retroviruses (ERVs) there is, for various reasons, considerable confusion regarding names assigned to families/groups of ERVs as well as individual ERV loci. Human ERVs have been studied in greater detail, and naming of HERVs in the scientific literature is somewhat confusing not just to the outsider. Without guidelines, confusion for ERVs in other species will also probably increase if those ERVs are studied in greater detail. Based on previous experience, this review highlights some of the problems when naming and classifying ERVs, and provides some guidance for detecting and characterizing ERV sequences. Because of the close relationship between ERVs and exogenous retroviruses (XRVs) it is reasonable to reconcile their classification with that of XRVs. We here argue that classification should be based on a combination of similarity, structural features, (inferred) function, and previous nomenclature. Because the RepBase system is widely employed in genome annotation, RepBase designations should be considered in further taxonomic efforts. To lay a foundation for a phylogenetically based taxonomy, further analyses of ERVs in many hosts are needed. A dedicated, permanent, international consortium would best be suited to integrate and communicate our current and future knowledge on repetitive, mobile elements in general to the scientific community.

Evolutionary conservation of orthoretroviral long terminal repeats (LTRs) and ab initio detection of single LTRs in genomic data.

BACKGROUND: Retroviral LTRs, paired or single, influence the transcription of both retroviral and non-retroviral genomic sequences. Vertebrate genomes contain many thousand endogenous retroviruses (ERVs) and their LTRs. Single LTRs are difficult to detect from genomic sequences without recourse to repetitiveness or presence in a proviral structure. Understanding of LTR structure increases understanding of LTR function, and of functional genomics. Here we develop models of orthoretroviral LTRs useful for detection in genomes and for structural analysis. PRINCIPAL FINDINGS: Although mutated, ERV LTRs are more numerous and diverse than exogenous retroviral (XRV) LTRs. Hidden Markov models (HMMs), and alignments based on them, were created for HML- (human MMTV-like), general-beta-, gamma- and lentiretroviruslike LTRs, plus a general-vertebrate LTR model. Training sets were XRV LTRs and RepBase LTR consensuses. The HML HMM was most sensitive and detected 87% of the HML LTRs in human chromosome 19 at 96% specificity. By combining all HMMs with a low cutoff, for screening, 71% of all LTRs found by RepeatMasker in chromosome 19 were found. HMM consensus sequences had a conserved modular LTR structure. Target site duplications (TG-CA), TATA (occasionally absent), an AATAAA box and a T-rich region were prominent features. Most of the conservation was located in, or adjacent to, R and U5, with evidence for stem loops. Several of the long HML LTRs contained long ORFs inserted after the second A rich module. HMM consensus alignment allowed comparison of functional features like transcriptional start sites (sense and antisense) between XRVs and ERVs. CONCLUSION: The modular conserved and redundant orthoretroviral LTR structure with three A-rich regions is reminiscent of structurally relaxed Giardia promoters. The five HMMs provided a novel broad range, repeat-independent, ab initio LTR detection, with prospects for greater generalisation, and insight into LTR structure, which may aid development of LTR-targeted pharmaceuticals.