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

Transcriptional Regulation of Endogenous Retroviruses and Their Misregulation in Human Diseases

Endogenous retroviruses (ERVs), deriving from exogenous retroviral infections of germ line cells occurred millions of years ago, represent ~8% of human genome. Most ERVs are highly inactivated because of the accumulation of mutations, insertions, deletions, and/or truncations. However, it is becoming increasingly apparent that ERVs influence host biology through genetic and epigenetic mechanisms under particular physiological and pathological conditions, which provide both beneficial and deleterious effects for the host. For instance, certain ERVs expression is essential for human embryonic development. Whereas abnormal activation of ERVs was found to be involved in numbers of human diseases, such as cancer and neurodegenerative diseases. Therefore, understanding the mechanisms of regulation of ERVs would provide insights into the role of ERVs in health and diseases. Here, we provide an overview of mechanisms of transcriptional regulation of ERVs and their dysregulation in human diseases.

Transcriptome Analyses Implicate Endogenous Retroviruses Involved in the Host Antiviral Immune System through the Interferon Pathway

Human endogenous retroviruses (HERVs) are the remains of ancient retroviruses that invaded our ancestors' germline cell and were integrated into the genome. The expression of HERVs has always been a cause for concern because of its association with various cancers and diseases. However, few previous studies have focused on specific activation of HERVs by viral infections. Our previous study has shown that dengue virus type 2 (DENV-2) infection induces the transcription of a large number of abnormal HERVs loci; therefore, the purpose of this study was to explore the relationship between exogenous viral infection and HERV activation further. In this study, we retrieved and reanalyzed published data on 21 transcriptomes of human cells infected with various viruses. We found that infection with different viruses could induce transcriptional activation of HERV loci. Through the comparative analysis of all viral datasets, we identified 43 key HERV loci that were up-regulated by DENV-2, influenza A virus, influenza B virus, Zika virus, measles virus, and West Nile virus infections. Furthermore, the neighboring genes of these HERVs were simultaneously up-regulated, and almost all such neighboring genes were interferon-stimulated genes (ISGs), which are enriched in the host's antiviral immune response pathways. Our data supported the hypothesis that activation of HERVs, probably via an interferon-mediated mechanism, plays an important role in innate immunity against viral infections.

Telescope: Characterization of the retrotranscriptome by accurate estimation of transposable element expression

Characterization of Human Endogenous Retrovirus (HERV) expression within the transcriptomic landscape using RNA-seq is complicated by uncertainty in fragment assignment because of sequence similarity. We present Telescope, a computational software tool that provides accurate estimation of transposable element expression (retrotranscriptome) resolved to specific genomic locations. Telescope directly addresses uncertainty in fragment assignment by reassigning ambiguously mapped fragments to the most probable source transcript as determined within a Bayesian statistical model. We demonstrate the utility of our approach through single locus analysis of HERV expression in 13 ENCODE cell types. When examined at this resolution, we find that the magnitude and breadth of the retrotranscriptome can be vastly different among cell types. Furthermore, our approach is robust to differences in sequencing technology and demonstrates that the retrotranscriptome has potential to be used for cell type identification. We compared our tool with other approaches for quantifying transposable element (TE) expression, and found that Telescope has the greatest resolution, as it estimates expression at specific TE insertions rather than at the TE subfamily level. Telescope performs highly accurate quantification of the retrotranscriptomic landscape in RNA-seq experiments, revealing a differential complexity in the transposable element biology of complex systems not previously observed. Telescope is available at https://github.com/mlbendall/telescope.

Human endogenous retrovirus-K (HML-2): a comprehensive review

The human genome contains a large number of retroviral elements acquired over the process of evolution, some of which are specific to primates. However, as many of these are defective or silenced through epigenetic changes, they were historically considered "junk DNA" and their potential role in human physiology or pathological circumstances have been poorly studied. The most recently acquired, human endogenous retrovirus-K (HERV-K), has multiple copies in the human genome and some of them have complete open reading frames that are transcribed and translated, especially in early embryogenesis. Phylogenetically, HERV-K is considered a supergroup of viruses. One of the subtypes, termed HML-2, seems to be the most active and hence, it is the best studied. Aberrant expression of HML-2 in adult tissues has been associated with certain types of cancer and with neurodegenerative diseases. This review discusses the discovery of these viruses, their classification, structure, regulation and potential for replication, physiological roles, and their involvement in disease pathogenesis. Finally, it presents different therapeutic approaches being considered to target these viruses.

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.

Human Endogenous Retroviruses in Neurological Diseases

The causes of multiple sclerosis and amyotrophic lateral sclerosis have long remained elusive. A new category of pathogenic components, normally dormant within human genomes, has been identified: human endogenous retroviruses (HERVs). These represent ∼8% of the human genome, and environmental factors have reproducibly been shown to trigger their expression. The resulting production of envelope (Env) proteins from HERV-W and HERV-K appears to engage pathophysiological pathways leading to the pathognomonic features of MS and ALS, respectively. Pathogenic HERV elements may thus provide a missing link in understanding these complex diseases. Moreover, their neutralization may represent a promising strategy to establish novel and more powerful therapeutic approaches.

Microarray analysis reveals global modulation of endogenous retroelement transcription by microbes

BACKGROUND

A substantial proportion of both the mouse and human genomes comprise of endogenous retroelements (REs), which include endogenous retroviruses. Over evolutionary time, REs accumulate inactivating mutations or deletions and thus lose the ability to replicate. Additionally, REs can be transcriptionally repressed by dedicated mechanisms of the host. Nevertheless, many of them still possess and express intact open reading frames, and their transcriptional activity has been associated with many physiological and pathological processes of the host. However, this association remains tenuous due to incomplete understanding of the mechanism by which RE transcription is regulated. Here, we use a bioinformatics tool to examine RE transcriptional activity, measured by microarrays, in murine and human immune cells responding to microbial stimulation.

RESULTS

Immune cell activation by microbial signals in vitro caused extensive changes in the transcription not only of the host genes involved in the immune response, but also of numerous REs. Modulated REs were frequently found near or embedded within similarly-modulated host genes. Focusing on probes reporting single-integration, intergenic REs, revealed extensive transcriptional responsiveness of these elements to microbial signals. Microbial stimulation modulated RE expression in a cell-intrinsic manner. In line with these results, the transcriptional activity of numerous REs followed characteristics in different tissues according to exposure to environmental microbes and was further heavily altered during viral infection or imbalances with intestinal microbiota, both in mice and humans.

CONCLUSIONS

Together, these results highlight the utility of improved methodologies in assessing RE transcription profiles in both archived and new microarray data sets. More importantly, application of this methodology suggests that immune activation, as a result of infection with pathogens or dysbiosis with commensal microbes, causes global modulation of RE transcription. RE responsiveness to external stimuli should, therefore, be considered in any association between RE transcription and disease.

Microarray-based identification of individual HERV loci expression: application to biomarker discovery in prostate cancer

The prostate-specific antigen (PSA) is the main diagnostic biomarker for prostate cancer in clinical use, but it lacks specificity and sensitivity, particularly in low dosage values¹​​. ‘How to use PSA' remains a current issue, either for diagnosis as a gray zone corresponding to a concentration in serum of 2.5-10 ng/ml which does not allow a clear differentiation to be made between cancer and noncancer² or for patient follow-up as analysis of post-operative PSA kinetic parameters can pose considerable challenges for their practical application^3,4. Alternatively, noncoding RNAs (ncRNAs) are emerging as key molecules in human cancer, with the potential to serve as novel markers of disease, e.g. PCA3 in prostate cancer^5,6 and to reveal uncharacterized aspects of tumor biology. Moreover, data from the ENCODE project published in 2012 showed that different RNA types cover about 62% of the genome. It also appears that the amount of transcriptional regulatory motifs is at least 4.5x higher than the one corresponding to protein-coding exons. Thus, long terminal repeats (LTRs) of human endogenous retroviruses (HERVs) constitute a wide range of putative/candidate transcriptional regulatory sequences, as it is their primary function in infectious retroviruses. HERVs, which are spread throughout the human genome, originate from ancestral and independent infections within the germ line, followed by copy-paste propagation processes and leading to multicopy families occupying 8% of the human genome (note that exons span 2% of our genome). Some HERV loci still express proteins that have been associated with several pathologies including cancer^7-10. We have designed a high-density microarray, in Affymetrix format, aiming to optimally characterize individual HERV loci expression, in order to better understand whether they can be active, if they drive ncRNA transcription or modulate coding gene expression. This tool has been applied in the prostate cancer field (Figure 1).

Transcriptional profiling of human endogenous retrovirus group HERV-K(HML-2) loci in melanoma

Recent studies suggested a role for the human endogenous retrovirus (HERV) group HERV-K(HML-2) in melanoma because of upregulated transcription and expression of HERV-K(HML-2)-encoded proteins. Very little is known about which HML-2 loci are transcribed in melanoma. We assigned >1,400 HML-2 cDNA sequences generated from various melanoma and related samples to genomic HML-2 loci, identifying a total of 23 loci as transcribed. Transcription profiles of loci differed significantly between samples. One locus was found transcribed only in melanoma-derived samples but not in melanocytes and might represent a marker for melanoma. Several of the transcribed loci harbor ORFs for retroviral Gag and/or Env proteins. Env-encoding loci were transcribed only in melanoma. Specific investigation of rec and np9 transcripts indicated transcription of protein encoding loci in melanoma and melanocytes hinting at the relevance of Rec and Np9 in melanoma. UVB irradiation changed transcription profiles of loci and overall transcript levels decreased in melanoma and melanocytes. We further identified transcribed HML-2 loci formed by reverse transcription of spliced HML-2 transcripts by L1 machinery or in a retroviral fashion, with loci potentially encoding HML-2-like proteins. We reveal complex, sample-specific transcription of HML-2 loci in melanoma and related samples. Identified HML-2 loci and proteins encoded by those loci are particularly relevant for further studying the role of HML-2 in melanoma. Transcription of HERVs appears as a complex mechanism requiring specific studies to elucidate which HERV loci are transcribed and how transcribed HERVs may be involved in disease.

HERV-E-mediated modulation of PLA2G4A transcription in urothelial carcinoma

Human endogenous retroviruses (HERV) and related elements account for more than 8% of the human genome and significantly contribute to the human transcriptome by long terminal repeat (LTR) promoter activity. In this context, HERVs are thought to intervene in the expression of adjacent genes by providing regulatory sequences (cis-effect) or via noncoding RNA including natural antisense transcripts. To address the potential impact of HERV activity in urothelial carcinoma, we comparatively analyzed the HERV transcription profiles in paired samples of non-malignant urothelium and urothelial carcinoma derived from 13 patients with bladder cancer by means of a retrovirus-specific microarray (RetroArray). We established a characteristic HERV signature consisting of six ubiquitously active HERV subgroups (E4-1, HERV-Rb, ERV9, HERV-K-T47D, NMWV3, HERV-KC4). The transcription pattern is largely identical in human urothelial carcinoma, non-malignant urothelial tissue, four tumor-derived cell lines and in a non-malignant urothelial cell line (UROtsa). Quantitative reverse transcriptase PCR (qRT-PCR) of HERV-E4-1, HERV-K(HML-6) and HERV-T(S71-TK1) revealed a bias to lower HERV activity in carcinoma samples compared to non-malignant tissue. Determination of active HERV-E4-1 loci by cloning and sequencing revealed six HERV-E4-1 proviral loci that are differentially regulated in urothelial carcinoma cells and normal tissue. Two full-length HERV-E4-1 proviruses, HERV-Ec1 and HERV-Ec6, are located in antisense orientation in introns of the genes PLA2G4A and RNGTT, respectively. PLA2G4A encodes a cytosolic phospholipase A2 (cPLA2) that is dysregulated in many human tumors. PLA2G4A and HERV-Ec1 displayed reciprocal transcript levels in 7 of 11 urothelial carcinoma patients. Moreover, reciprocal shifts were observed after treatment of UROtsa cells with HERV-Ec1 and PLA2G4A-directed siRNAs or 5-aza-2′-deoxycytidine (aza-dC) pointing to an antagonistic regulation of PLA2G4A and HERV-Ec1 transcription in human urothelial cells. We suggest that transcription of HERV-Ec1 contributes to fine tuning of cPLA2 expression, thereby facilitating tumorigenesis.

Microarray-based sketches of the HERV transcriptome landscape

Human endogenous retroviruses (HERVs) are spread throughout the genome and their long terminal repeats (LTRs) constitute a wide collection of putative regulatory sequences. Phylogenetic similarities and the profusion of integration sites, two inherent characteristics of transposable elements, make it difficult to study individual locus expression in a large-scale approach, and historically apart from some placental and testis-regulated elements, it was generally accepted that HERVs are silent due to epigenetic control. Herein, we have introduced a generic method aiming to optimally characterize individual loci associated with 25-mer probes by minimizing cross-hybridization risks. We therefore set up a microarray dedicated to a collection of 5,573 HERVs that can reasonably be assigned to a unique genomic position. We obtained a first view of the HERV transcriptome by using a composite panel of 40 normal and 39 tumor samples. The experiment showed that almost one third of the HERV repertoire is indeed transcribed. The HERV transcriptome follows tropism rules, is sensitive to the state of differentiation and, unexpectedly, seems not to correlate with the age of the HERV families. The probeset definition within the U3 and U5 regions was used to assign a function to some LTRs (i.e. promoter or polyA) and revealed that (i) autonomous active LTRs are broadly subjected to operational determinism (ii) the cellular gene density is substantially higher in the surrounding environment of active LTRs compared to silent LTRs and (iii) the configuration of neighboring cellular genes differs between active and silent LTRs, showing an approximately 8 kb zone upstream of promoter LTRs characterized by a drastic reduction in sense cellular genes. These gathered observations are discussed in terms of virus/host adaptive strategies, and together with the methods and tools developed for this purpose, this work paves the way for further HERV transcriptome projects.

Smoking increases transcription of human endogenous retroviruses in a newly established in vitro cell model and in normal urothelium

Human endogenous retroviruses (HERVs) accounting for 9% of the human genome are considered as surrogate markers for genetic instability and as a driving force of genetic variation. Moreover, they modulate regular gene activities and give rise to expression of disease-associated peptides that may serve as diagnostic markers or even targets for T cell-based immune responses. To date, no data are available on the potential link between urothelial carcinogenesis, HERV activity, and tobacco smoking, the main risk for bladder cancer. Here, we report on potential alterations in HERV transcription induced by smoking in a newly established in vitro model and in human urothelium. Normal human dermal fibroblasts were cultivated with urine from never (n = 6) and current smokers (n = 6) and transcription levels for the HERV subfamilies HERV-E 4-1, HERV-T S71-TK1, and HERV-K HML-6 were measured by quantitative real-time PCR. Tendencies toward increased mean transcript levels were detected for cells treated with urine from current smokers. Equally, activity measured in human urothelium supported an increase of HERV transcription in current smokers (n = 9) compared to never smokers (n = 4).

Potential mechanisms of endogenous retroviral-mediated genomic instability in human cancer

Malignancy results from a complex combination of genetic and epigenetic changes, the full effects of which are still largely unknown. Here we summarize current knowledge of the origin, retrotranspositional activity, epigenetic state, and transcription of human endogenous retroviruses (HERVs), and then discuss the potential effects of their deregulation in cancer. Evidence suggests that cancer-associated epigenetic changes most likely underlie potential HERV-mediated effects on genome and transcriptome instability and may play a role in malignancy. Despite our currently limited understanding of the importance of HERVs or other transposable elements in cancer development, we believe that the emerging era of high-throughput sequencing of cancer genomes, epigenomes, and transcriptomes will provide unprecedented opportunities to investigate these roles in the future.

Broadly targeted triplex real-time PCR detection of influenza A, B and C viruses based on the nucleoprotein gene and a novel "MegaBeacon" probe strategy

A PCR assay that covers animal and human influenza A, B and C viruses, i.e., most of Orthomyxoviridae, is needed. Influenza types are distinguished based on differences in the nucleoprotein (NP) present in the virus. Conserved NP regions were therefore used to design a TaqMan-based triplex reverse transcription real-time PCR method. Variability of influenza A within the probe target region mandated the development of a novel molecular beacon, the "Mega" molecular beacon (MegaBeacon; MegB), for the detection of influenza A with this method. MegaBeacon is a mismatch-tolerant molecular beacon that is also a TaqMan probe. The triplex method (3QPCR-MegB) was evaluated with influenza A isolates covering 18 HxNx combinations, two influenza B isolates, and five Japanese influenza C isolates, as well as influenza A, B and C synthetic DNA targets. One to ten viral RNA and cDNA genome equivalents were detected per PCR reaction for influenza A, B and C. Seventy-one human nasopharyngeal aspirates from respiratory infections yielded 30 influenza A, 11 influenza B and 0 influenza C with 3QPCR-MegB, where immunofluorescence (IF) found 28 influenza A and 10 influenza B. 3QPCR-MegB was more mismatch-tolerant than a variant PCR with an influenza A TaqMan probe (3QPCR) and is a sensitive and rational method to detect influenza viruses of animal and human origin. MegaBeacon probes hold promise for variable target nucleic acids.

Broadly targeted multiprobe QPCR for detection of coronaviruses: Coronavirus is common among mallard ducks (Anas platyrhynchos)

Coronaviruses (CoVs) can cause trivial or fatal disease in humans and in animals. Detection methods for a wide range of CoVs are needed, to understand viral evolution, host range, transmission and maintenance in reservoirs. A new concept, “Multiprobe QPCR”, which uses a balanced mixture of competing discrete non- or moderately degenerated nuclease degradable (TaqMan^®) probes was employed. It provides a broadly targeted and rational single tube real-time reverse transcription PCR (“NQPCR”) for the generic detection and discovery of CoV. Degenerate primers, previously published, and the new probes, were from a conserved stretch of open reading frame 1b, encoding the replicase. This multiprobe design reduced the degree of probe degeneration, which otherwise decreases the sensitivity, and allowed a preliminary classification of the amplified sequence directly from the QPCR trace. The split probe strategy allowed detection of down to 10 viral nucleic acid equivalents of CoV from all known CoV groups. Evaluation was with reference CoV strains, synthetic targets, human respiratory samples and avian fecal samples. Infectious-Bronchitis-Virus (IBV)-related variants were found in 7 of 35 sample pools, from 100 wild mallards (Anas platyrhynchos). Ducks may spread and harbour CoVs. NQPCR can detect a wide range of CoVs, as illustrated for humans and ducks.

Expression patterns of transcribed human endogenous retrovirus HERV-K(HML-2) loci in human tissues and the need for a HERV Transcriptome Project

Background

A significant proportion of the human genome is comprised of human endogenous retroviruses (HERVs). HERV transcripts are found in every human tissue. Expression of proviruses of the HERV-K(HML-2) family has been associated with development of human tumors, in particular germ cell tumors (GCT). Very little is known about transcriptional activity of individual HML-2 loci in human tissues, though.

Results

By employing private nucleotide differences between loci, we assigned ~1500 HML-2 cDNAs to individual HML-2 loci, identifying, in total, 23 transcriptionally active HML-2 proviruses. Several loci are active in various human tissue types. Transcription levels of some HML-2 loci appear higher than those of other loci. Several HML-2 Rec-encoding loci are expressed in GCT and non-GCT tissues. A provirus on chromosome 22q11.21 appears strongly upregulated in pathologic GCT tissues and may explain high HML-2 Gag protein levels in GCTs. Presence of Gag and Env antibodies in GCT patients is not correlated with activation of individual loci. HML-2 proviruses previously reported capable of forming an infectious HML-2 variant are transcriptionally active in germ cell tissue. Our study furthermore shows that Expressed Sequence Tag (EST) data are insufficient to describe transcriptional activity of HML-2 and other HERV loci in tissues of interest.

Conclusion

Our, to date, largest-scale study reveals in greater detail expression patterns of individual HML-2 loci in human tissues of clinical interest. Moreover, large-scale, specialized studies are indicated to better comprehend transcriptional activity and regulation of HERVs. We thus emphasize the need for a specialized HERV Transcriptome Project.

Variable transcriptional activity of endogenous retroviruses in human breast cancer

Human endogenous retroviruses (HERVs) account for up to 9% of the human genome and include more than 800 elements related to betaretroviruses. While mouse mammary tumor virus (MMTV) is the accepted etiological agent of mammary tumors in mice, the role of retroviral elements in human breast cancer remains elusive. Here, we performed a comprehensive microarray-based analysis of overall retroviral transcriptional activities in 46 mammary gland tissue specimens representing pairs of nonmalignant and tumor samples from 23 patients. An analysis of nonmalignant tissue samples revealed a distinct, mammary gland-specific HERV expression profile that consists of 18 constitutively active HERV taxa. For corresponding tumor samples, a general trend toward lower levels of HERV transcription was observed, suggesting common regulatory mechanisms. In various subsets of patients, however, increased transcript levels of single class I HERV families (HERV-T, HERV-E, and HERV-F) and several class II families, including HML-6, were detected. An analysis of transcribed HML-6 sequences revealed either the activation of some or the increased activity of several proviral loci. No evidence for MMTV or human MMTV-like virus transcripts was found, indicating that transcriptionally active, MMTV analogous, exogenous viruses were not present in the breast cancer samples analyzed.

Quantitative multiplex degenerate PCR for human endogenous retrovirus expression profiling

Expression of human endogenous retroviruses (HERV) has been recurrently observed during cellular differentiation or transformation processes in both cell culture and in vivo. Quantitative approaches that analyze variations in HERV transcription could therefore be valuable for cancer diagnosis. We have developed a quantitative assay combining multiplex degenerate PCR (MD-PCR) and a colorimetric Oligo Sorbent Array (OLISA). Quantification of the expression of these multifamily genes relies on the optimization of the amplification primer mix, that is, the primer degeneracy, the relative concentration of each primer and the total amount of primer. Amplification products of each of the nine studied HERV families are independently and specifically detected and quantified using the OLISA microarray. This method constitutes an improvement over previous pan-retrovirus amplification-based methods, which are mainly qualitative. Furthermore, as MD-PCR/OLISA simultaneously monitors several HERV families, it challenges single-family quantitative RT-PCR. Last, the protocol below provides general rules for the design of MD-PCR applications. Once primers have been designed and optimized, the procedure can be completed in 2 days.

Automated recognition of retroviral sequences in genomic data--RetroTector

Eukaryotic genomes contain many endogenous retroviral sequences (ERVs). ERVs are often severely mutated, therefore difficult to detect. A platform independent (Java) program package, RetroTector^© (ReTe), was constructed. It has three basic modules: (i) detection of candidate long terminal repeats (LTRs), (ii) detection of chains of conserved retroviral motifs fulfilling distance constraints and (iii) attempted reconstruction of original retroviral protein sequences, combining alignment, codon statistics and properties of protein ends. Other features are prediction of additional open reading frames, automated database collection, graphical presentation and automatic classification. ReTe favors elements >1000-bp long due to its dependence on order of and distances between retroviral fragments. It detects single or low-copy-number elements. ReTe assigned a ‘retroviral’ score of 890–2827 to 10 exogenous retroviruses from seven genera, and accurately predicted their genes. In a simulated model, ReTe was robust against mutational decay. The human genome was analyzed in 1–2 days on a LINUX cluster. Retroviral sequences were detected in divergent vertebrate genomes. Most ReTe detected chains were coincident with Repeatmasker output and the HERVd database. ReTe did not report most of the evolutionary old HERV-L related and MalR sequences, and is not yet tailored for single LTR detection. Nevertheless, ReTe rationally detects and annotates many retroviral sequences.

Methods for estimating human endogenous retrovirus activities from EST databases

Background

Human endogenous retroviruses (HERVs) are surviving traces of ancient retrovirus infections and now reside within the human DNA. Recently HERV expression has been detected in both normal tissues and diseased patients. However, the activities (expression levels) of individual HERV sequences are mostly unknown.

Results

We introduce a generative mixture model, based on Hidden Markov Models, for estimating the activities of the individual HERV sequences from EST (expressed sequence tag) databases. We use the model to estimate the relative activities of 181 HERVs. We also empirically justify a faster heuristic method for HERV activity estimation and use it to estimate the activities of 2450 HERVs. The majority of the HERV activities were previously unknown.

Conclusion

(i) Our methods estimate activity accurately based on experiments on simulated data. (ii) Our estimate on real data shows that 7% of the HERVs are active. The active ones are spread unevenly into HERV groups and relatively uniformly in terms of estimated age. HERVs with the retroviral env gene are more often active than HERVs without env. Few of the active HERVs have open reading frames for retroviral proteins.