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

Murine Gammaretrovirus Group G3 Was Not Found in Swedish Patients with Myalgic Encephalomyelitis/Chronic Fatigue Syndrome and Fibromyalgia

BACKGROUND: The recent report of gammaretroviruses of probable murine origin in humans, called xenotropic murine retrovirus related virus (XMRV) and human murine leukemia virus related virus (HMRV), necessitated a bioinformatic search for this virus in genomes of the mouse and other vertebrates, and by PCR in humans. RESULTS: Three major groups of murine endogenous gammaretroviruses were identified. The third group encompassed both exogenous and endogenous Murine Leukemia Viruses (MLVs), and most XMRV/HMRV sequences reported from patients suffering from myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS). Two sensitive real-time PCRs for this group were developed. The predicted and observed amplification range for these and three published XMRV/HMRV PCRs demonstrated conspicuous differences between some of them, partly explainable by a recombinatorial origin of XMRV. Three reverse transcription real-time PCRs (RTQPCRs), directed against conserved and not overlapping stretches of env, gag and integrase (INT) sequences of XMRV/HMRV were used on human samples. White blood cells from 78 patients suffering from ME/CFS, of which 30 patients also fulfilled the diagnostic criteria for fibromyalgia (ME/CFS/FM) and in 7 patients with fibromyalgia (FM) only, all from the Gothenburg area of Sweden. As controls we analyzed 168 sera from Uppsala blood donors. We controlled for presence and amplifiability of nucleic acid and for mouse DNA contamination. To score as positive, a sample had to react with several of the XMRV/HMRV PCRs. None of the samples gave PCR reactions which fulfilled the positivity criteria. CONCLUSIONS: XMRV/HMRV like proviruses occur in the third murine gammaretrovirus group, characterized here. PCRs developed by us, and others, approximately cover this group, except for the INT RTQPCR, which is rather strictly XMRV specific. Using such PCRs, XMRV/HMRV could not be detected in PBMC and plasma samples from Swedish patients suffering from ME/CFS/FM, and in sera from Swedish blood donors

Genome-wide analysis reveals the extent of EAV-HP integration in domestic chicken

Background: EAV-HP is an ancient retrovirus pre-dating Gallus speciation, which continues to circulate in modern chicken populations, and led to the emergence of avian leukosis virus subgroup J causing significant economic losses to the poultry industry. We mapped EAV-HP integration sites in Ethiopian village chickens, a Silkie, Taiwan Country chicken, red junglefowl Gallusgallus and several inbred experimental lines using whole-genome sequence data. Results: An average of 75.22 ± 9.52 integration sites per bird were identified, which collectively group into 279 intervals of which 5% are common to 90% of the genomes analysed and are suggestive of pre-domestication integration events. More than a third of intervals are specific to individual genomes, supporting active circulation of EAV-HP in modern chickens. Interval density is correlated with chromosome length (P<2.31−6), and 27 % of intervals are located within 5 kb of a transcript. Functional annotation clustering of genes reveals enrichment for immune-related functions (P<0.05). Conclusions: Our results illustrate a non-random distribution of EAV-HP in the genome, emphasising the importance it may have played in the adaptation of the species, and provide a platform from which to extend investigations on the co-evolutionary significance of endogenous retroviral genera with their hosts

A new look at the LTR retrotransposon content of the chicken genome

BACKGROUND: LTR retrotransposons contribute approximately 10 % of the mammalian genome, but it has been previously reported that there is a deficit of these elements in the chicken relative to both mammals and other birds. A novel LTR retrotransposon classification pipeline, LocaTR, was developed and subsequently utilised to re-examine the chicken LTR retrotransposon annotation, and determine if the proposed chicken deficit is biologically accurate or simply a technical artefact. RESULTS: Using LocaTR 3.01 % of the chicken galGal4 genome assembly was annotated as LTR retrotransposon-derived elements (nearly double the previous annotation), including 1,073 that were structurally intact. Element distribution is significantly correlated with chromosome size and is non-random within each chromosome. Elements are significantly depleted within coding regions and enriched in gene sparse areas of the genome. Over 40 % of intact elements are found in clusters, unrelated by age or genera, generally in poorly recombining regions. The transcription of most LTR retrotransposons were suppressed or incomplete, but individual domain and full length retroviral transcripts were produced in some cases, although mostly with regularly interspersed stop codons in all reading frames. Furthermore, RNAseq data from 23 diverse tissues enabled greater characterisation of the co-opted endogenous retrovirus Ovex1. This gene was shown to be expressed ubiquitously but at variable levels across different tissues. LTR retrotransposon content was found to be very variable across the avian lineage and did not correlate with either genome size or phylogenetic position. However, the extent of previous, species-specific LTR retrotransposon annotation appears to be a confounding factor. CONCLUSIONS: Use of the novel LocaTR pipeline has nearly doubled the annotated LTR retrotransposon content of the chicken genome compared to previous estimates. Further analysis has described element distribution, clustering patterns and degree of expression in a variety of adult tissues, as well as in three embryonic stages. This study also enabled better characterisation of the co-opted gamma retroviral envelope gene Ovex1. Additionally, this work suggests that there is no deficit of LTR retrotransposons within the Galliformes relative to other birds, or to mammalian genomes when scaled for the three-fold difference in genome size

Identification and analysis of HML2 sequences in human genome assembly GRCh37/hg19

Background Human endogenous retroviruses (HERVs) originated from exogenous retroviral infections of the human germ line cells and spread in the human population through vertical transmission over millions of years. Among HERVs, the HML2 proviruses [1] are the most recently integrated and show the most intact proviral genomes. HML-2 expression has tentatively been associated with different pathological conditions, including Hodgkin’s lymphoma, melanoma, breast and testicular cancer. A comprehensive recent study identified 91 HML2 proviruses [2]. Material and methods Human genome (assembly GRCh 37/hg19) was analyzed with RetroTector (ReTe) version 1.01 [3]. ReTe was run on a machine with 4 6-core Xeon processors, 2.66Ghz each, 256 Gb of RAM and 4 Tb of disks, with an estimated execution time of 1-2 days. BLASTN, using HML consensuses (Blikstad et al, unpublished) and the May 2013 Repeatmasker library, ENSEMBL and MEGA5 were used, in successive steps, for classification and identification of locus position and phylogenetic inference. Time since integration was inferred using a neutral substitution rate between cognate LTRs of 0.2 mutations per million years. Results ReTe [2] identified more than 120 HML2 proviruses, many of which were not previously reported, accounting for roughly 0.01% of the total human genome. Among the identified HML2 proviruses more than 50% are ≥ 8000 bp in length and more than 50% have both LTRs. HML2 proviruses bordering to HML1, HML3, HML9 and HML10, as well as recombinant proviruses containing HML2 sequences were detected. HML2 proviruses were present in all chromosomes and showed to form clusters, particularly in chromosomes 1, 4, 8 and 19. Open reading frames (ORFs) predicted by ReTe revealed that 21 proviruses have at least 1 ORFs in gag, pro, pol and env genes, while 6 had ORFs in 3 genes. Age analysis versus reductions of ORFs and proviral length was performed. Phylogenetic analyses were performed with whole element DNA, concatenated Gag, Pro and Pol amino acid sequences, and Pol amino acid sequences. Conclusions In an attempt to establish a comprehensive catalog of HML2 proviruses that could set the basis for further research, we detected over 120 HML2 proviruses and performed a first characterization of them