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.

Annotation and visualization of endogenous retroviral sequences using the Distributed Annotation System (DAS) and eBioX.

BACKGROUND: The Distributed Annotation System (DAS) is a widely used network protocol for sharing biological information. The distributed aspects of the protocol enable the use of various reference and annotation servers for connecting biological sequence data to pertinent annotations in order to depict an integrated view of the data for the final user. RESULTS: An annotation server has been devised to provide information about the endogenous retroviruses detected and annotated by a specialized in silico tool called RetroTector. We describe the procedure to implement the DAS 1.5 protocol commands necessary for constructing the DAS annotation server. We use our server to exemplify those steps. Data distribution is kept separated from visualization which is carried out by eBioX, an easy to use open source program incorporating multiple bioinformatics utilities. Some well characterized endogenous retroviruses are shown in two different DAS clients. A rapid analysis of areas free from retroviral insertions could be facilitated by our annotations. CONCLUSION: The DAS protocol has shown to be advantageous in the distribution of endogenous retrovirus data. The distributed nature of the protocol is also found to aid in combining annotation and visualization along a genome in order to enhance the understanding of ERV contribution to its evolution. Reference and annotation servers are conjointly used by eBioX to provide visualization of ERV annotations as well as other data sources. Our DAS data source can be found in the central public DAS service repository, http://www.dasregistry.org, or at http://loka.bmc.uu.se/das/sources.

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.

Use of endogenous retroviral sequences (ERVs) and structural markers for retroviral phylogenetic inference and taxonomy.

BACKGROUND: Endogenous retroviral sequences (ERVs) are integral parts of most eukaryotic genomes and vastly outnumber exogenous retroviruses (XRVs). ERVs with a relatively complete structure were retrieved from the genetic archives of humans and chickens, diametrically opposite representatives of vertebrate retroviruses (over 3300 proviruses), and analyzed, using a bioinformatic program, RetroTector, developed by us. This rich source of proviral information, accumulated in a local database, and a collection of XRV sequences from the literature, allowed the reconstruction of a Pol based phylogenetic tree, more extensive than previously possible. The aim was to find traits useful for classification and evolutionary studies of retroviruses. Some of these traits have been used by others, but they are here tested in a wider context than before. RESULTS: In the ERV collection we found sequences similar to the XRV-based genera: alpha-, beta-, gamma-, epsilon- and spumaretroviruses. However, the occurrence of intermediates between them indicated an evolutionary continuum and suggested that taxonomic changes eventually will be necessary. No delta or lentivirus representatives were found among ERVs. Classification based on Pol similarity is congruent with a number of structural traits. Acquisition of dUTPase occurred three times in retroviral evolution. Loss of one or two NC zinc fingers appears to have occurred several times during evolution. Nucleotide biases have been described earlier for lenti-, delta- and betaretroviruses and were here confirmed in a larger context. CONCLUSION: Pol similarities and other structural traits contribute to a better understanding of retroviral phylogeny. "Global" genomic properties useful in phylogenies are i.) translational strategy, ii.) number of Gag NC zinc finger motifs, iii.) presence of Pro N-terminal dUTPase (dUTPasePro), iv.) presence of Pro C-terminal G-patch and v.) presence of a GPY/F motif in the Pol integrase (IN) C-terminal domain. "Local" retroviral genomic properties useful for delineation of lower level taxa are i.) host species range, ii.) nucleotide compositional bias and iii.) LTR lengths.

Self-organizing map-based discovery and visualization of human endogenous retroviral sequence groups.

About 8 per cent of the human genome consists of human endogenous retroviral sequences (HERVs), which are remains from ancient infections. The HERVs may give rise to transcripts or affect the expression of human genes. The first step in understanding HERV function is to classify HERVs into families. In this work we study the relationships of existing HERV families and detect potentially new HERV families. A Median Self-Organizing Map (SOM), a SOM for non-vectorial data, is used to group and visualize a collection of 3661 HERVs. The SOM-based analysis is complemented with estimates of the reliability of the results. A novel trustworthiness visualization method is used to estimate which parts of the SOM visualization are reliable and which not. The reliability of extracted interesting HERV groups is verified by a bootstrap procedure suitable for SOM visualization-based analysis. The SOM detects a group of epsilonretroviral sequences and a group of ERV9, HERVW, and HUERSP3 sequences which suggests that ERV9 and HERVW sequences may have a common origin.

Functional evaluation of cerebral microembolization in the rat.

This study investigates the effects of cerebral microembolism on motor performance and risk assessment behavior in the rat. Cerebral infarcts were produced in rats by injecting small plastic beads into the left heart ventricle under short-acting anesthesia. The functional outcome was tested 24 h later by subjecting the animals to a series of consecutive behavioral tests. Thereafter, the rats were anesthetized and underwent magnetic resonance imaging. On average about seven infarcts per brain were found. The volume of the individual infarcts was largest in the hippocampus (mean=4.26 mm(3)) and smallest in the white matter (mean=0.83 mm(3)). Embolized animals performed spontaneous and evident locomotion. The activity was, however, significantly decreased compared to rats treated with vehicle. More specific tests for motor ability revealed reduced gait capacity and muscular strength. A significant relationship was found between behaviors reflecting motor ability and the total volume of infarcted tissue in the brain stem, cortex and cerebellum. Also the behavioral profile of risk and benefit assessment was found to be altered by the microembolization. It is concluded that the combination of the microembolization method and behavioral tests provides a valuable tool for further studies of the pathophysiology of, and potential treatment for, cerebral infarction.

Conserved structure and inferred evolutionary history of long terminal repeats (LTRs).

BACKGROUND: Long terminal repeats (LTRs, consisting of U3-R-U5 portions) are important elements of retroviruses and related retrotransposons. They are difficult to analyse due to their variability.The aim was to obtain a more comprehensive view of structure, diversity and phylogeny of LTRs than hitherto possible. RESULTS: Hidden Markov models (HMM) were created for 11 clades of LTRs belonging to Retroviridae (class III retroviruses), animal Metaviridae (Gypsy/Ty3) elements and plant Pseudoviridae (Copia/Ty1) elements, complementing our work with Orthoretrovirus HMMs. The great variation in LTR length of plant Metaviridae and the few divergent animal Pseudoviridae prevented building HMMs from both of these groups.Animal Metaviridae LTRs had the same conserved motifs as retroviral LTRs, confirming that the two groups are closely related. The conserved motifs were the short inverted repeats (SIRs), integrase recognition signals (5´TGTTRNR…YNYAACA 3´); the polyadenylation signal or AATAAA motif; a GT-rich stretch downstream of the polyadenylation signal; and a less conserved AT-rich stretch corresponding to the core promoter element, the TATA box. Plant Pseudoviridae LTRs differed slightly in having a conserved TATA-box, TATATA, but no conserved polyadenylation signal, plus a much shorter R region.The sensitivity of the HMMs for detection in genomic sequences was around 50% for most models, at a relatively high specificity, suitable for genome screening.The HMMs yielded consensus sequences, which were aligned by creating an HMM model (a 'Superviterbi' alignment). This yielded a phylogenetic tree that was compared with a Pol-based tree. Both LTR and Pol trees supported monophyly of retroviruses. In both, Pseudoviridae was ancestral to all other LTR retrotransposons. However, the LTR trees showed the chromovirus portion of Metaviridae clustering together with Pseudoviridae, dividing Metaviridae into two portions with distinct phylogeny. CONCLUSION: The HMMs clearly demonstrated a unitary conserved structure of LTRs, supporting that they arose once during evolution. We attempted to follow the evolution of LTRs by tracing their functional foundations, that is, acquisition of RNAse H, a combined promoter/ polyadenylation site, integrase, hairpin priming and the primer binding site (PBS). Available information did not support a simple evolutionary chain of events.

The first sequenced carnivore genome shows complex host-endogenous retrovirus relationships.

Host-retrovirus interactions influence the genomic landscape and have contributed substantially to mammalian genome evolution. To gain further insights, we analyzed a female boxer (Canis familiaris) genome for complexity and integration pattern of canine endogenous retroviruses (CfERV). Intriguingly, the first such in-depth analysis of a carnivore species identified 407 CfERV proviruses that represent only 0.15% of the dog genome. In comparison, the same detection criteria identified about six times more HERV proviruses in the human genome that has been estimated to contain a total of 8% retroviral DNA including solitary LTRs. These observed differences in man and dog are likely due to different mechanisms to purge, restrict and protect their genomes against retroviruses. A novel group of gammaretrovirus-like CfERV with high similarity to HERV-Fc1 was found to have potential for active retrotransposition and possibly lateral transmissions between dog and human as a result of close interactions during at least 10.000 years. The CfERV integration landscape showed a non-uniform intra- and inter-chromosomal distribution. Like in other species, different densities of ERVs were observed. Some chromosomal regions were essentially devoid of CfERVs whereas other regions had large numbers of integrations in agreement with distinct selective pressures at different loci. Most CfERVs were integrated in antisense orientation within 100 kb from annotated protein-coding genes. This integration pattern provides evidence for selection against CfERVs in sense orientation relative to chromosomal genes. In conclusion, this ERV analysis of the first carnivorous species supports the notion that different mammals interact distinctively with endogenous retroviruses and suggests that retroviral lateral transmissions between dog and human may have occurred.

Divergent patterns of recent retroviral integrations in the human and chimpanzee genomes: probable transmissions between other primates and chimpanzees.

The human genome is littered by endogenous retrovirus sequences (HERVs), which constitute up to 8% of the total genomic sequence. The sequencing of the human (Homo sapiens) and chimpanzee (Pan troglodytes) genomes has facilitated the evolutionary study of ERVs and related sequences. We screened both the human genome (version hg16) and the chimpanzee genome (version PanTro1) for ERVs and conducted a phylogenetic analysis of recent integrations. We found a number of recent integrations within both genomes. They segregated into four groups. Two larger gammaretrovirus-like groups (PtG1 and PtG2) occurred in chimpanzees but not in humans. The PtG sequences were most similar to two baboon ERVs and a macaque sequence but neither to other chimpanzee ERVs nor to any human gammaretrovirus-like ERVs. The pattern was consistent with cross-species transfer via predation. This appears to be an example of horizontal transfer of retroviruses with occasional fixation in the germ line.

Sequence variability, gene structure, and expression of full-length human endogenous retrovirus H.

Recently, we identified and classified 926 human endogenous retrovirus H (HERV-H)-like proviruses in the human genome. In this paper, we used the information to, in silico, reconstruct a putative ancestral HERV-H. A calculated consensus sequence was nearly open in all genes. A few manual adjustments resulted in a putative 9-kb HERV-H provirus with open reading frames (ORFs) in gag, pro, pol, and env. Long terminal repeats (LTRs) differed by 1.1%, indicating proximity to an integration event. The gag ORF was extended upstream of the normal myristylation start site. There was a long leader (including a "pre-gag" ORF) region positioned like the N terminus of murine leukemia virus (MLV) "glyco-Gag," potentially encoding a proline- and serine-rich domain remotely similar to MLV pp12. Another ORF, starting inside the 5' LTR, had no obvious similarity to known protein domains. Unlike other hitherto described gammaretroviruses, the reconstructed Gag had two zinc finger motifs. Alternative splicing of sequences related to the HERV-H consensus was confirmed using dbEST data. env transcripts were most prevalent in colon tumors, but also in normal testis. We found no evidence for full length env transcripts in the dbEST. HERV-H had a markedly skewed nucleotide composition, disfavoring guanine and favoring cytidine. We conclude that the HERV-H consensus shared a gene arrangement common to gammaretroviruses with gag separated by stop codon from pro-pol in the same reading frame, while env resides in another reading frame. There was also alternative splicing. HERV-H consensus yielded new insights in gammaretroviral evolution and will be useful as a model in studies on expression and function.

Effects of oxygen and carbon dioxide on human retinal circulation: an investigation using blue field simulation and scanning laser ophthalmoscopy.

PURPOSE: To investigate the quantitative response of retinal circulation in normal human eyes to oxygen and carbon dioxide provocation and to gain more information about retinal mean transit time (MTT). METHODS: Healthy subjects were studied while breathing: room air (O(2) baseline); pure O(2); room air (CO(2) baseline), and CO(2) in room air. The CO(2) level was brought to 30% above individual baseline level. During each gas exposure, blue field simulation was used to estimate the velocity of macular leucocytes and fluorescein angiography was used to determine MTT, arteriovenous passage time (AVP) and vessel diameters. RESULTS: Breathing pure O(2) caused a reduction in macular leucocyte velocity of about 13%, a non-significant change in MTT and AVP from baseline, and a 10% reduction in the diameter of temporal retinal arteries and veins. Breathing CO(2) caused an increase in macular leucocyte velocity of 21%, a non-significant change in MTT and AVP, and no change in retinal vascular diameters. CONCLUSION: Both O(2) and CO(2) had the expected effects on retinal blood flow, as judged by leucocyte velocity in the macular capillaries, but little effect on retinal MTT, indicating concomitant changes in flow and retinal vascular volume. Measuring the diameter of retinal vessels before and during pure O(2) breathing is a simple and effective way of determining the reactivity of retinal vessels.

ERV3 and related sequences in humans: structure and RNA expression.

The ERV3 locus at chromosome 7q11 is a much studied human endogenous retroviral (HERV) sequence, owing to an env open reading frame (ORF) and placental RNA and protein expression. An analysis of the human genome demonstrated that ERV3 is one of a group of 41 highly related elements (ERV3-like HERVs) which use proline, isoleucine, or arginine tRNA in their primer binding sites. In addition to elements closely related to ERV3, the group included the previously known retinoic acid-inducible element, RRHERVI, also referred to as HERV15, but was separate from the related HERV-E elements. The ERV3-like elements are defective. The only element with an ORF among gag, pro, pol, and env genes was the env ORF of the original ERV3 locus. A search in dbEST revealed ERV3 RNA expression in placenta, skin, carcinoid tumor, and adrenal glands. Expression was also studied with newly developed real-time quantitative PCRs (QPCR) of ERV3 and HERV-E(4-1) env sequences. Results from a novel histone 3.3 RNA QPCR result served as the expression control. QPCR results for ERV3 were compatible with previously published results, with a stronger expression in adrenal gland and placenta than in 15 other human tissues. The expression of the envelope (env) of ERV3 at chromosome 7q11 was also studied by using stringent in situ hybridization. Expression was found in corpus luteum, testis, adrenal gland, Hassal's bodies in thymus, brown fat, pituitary gland, and epithelium of the lung. We conclude that ERV3 env is most strongly expressed in adrenal and sebaceous glands as well as in placenta.

Definition and variation of human endogenous retrovirus H.

We defined the abundant human endogenous retrovirus group HERV-H based on pol similarity. Among 3661 pol-containing elements, 1124 integrations were similar to HERV-H RGH2 pol using translated pol sequences. A clustering procedure lessened these to 234 representatives, amenable to detailed study. Among the 1124, 926 clustered into HERV-H and 106 into adjacent HERV-H-like, the remainder being more distant to HERV-H. The HERV-H group was divided into RTVLH2-like (705 elements) and RGH2-like (77 elements) subgroups. Among 926 HERV-H, LTR differences were 1-33%, 10% had env, 78% had gag, 66% had a histidine primer binding site (PBS), and 3% (both subgroups) had a phenylalanine PBS. Allelic differences in env were studied using a convenient temperature gradient gel electrophoresis (TGGE) method and a genomic single nucleotide polymorphism (SNP) search. A pattern of abundant defective elements and less abundant less defective ones led us to formulate a "midwife" master model where more complete elements help the others in trans to transpose.