Identification of the avian myeloblastosis virus genome. II. Restriction endonuclease analysis of DNA from lambda proviral recombinants and leukemic myeoblast clones.

Two lambda proviral DNA recombinants were characterized with a number of restriction endonucleases. One recombinant contained a complete presumptive avian myeloblastosis virus (AMV) provirus flanked by cellular sequences on either side, and the second recombinant contained 85% of a myeloblastosis-associated virus type 1 (MAV-1)-like provirus with cellular sequences adjacent to the 5' end of the provirus. Comparing the restriction maps for the proviral DNAs contained in each lambda hybrid showed that the putative AMV and MAV-1-like genomes shared identical enzyme sites for 3.6 megadaltons beginning at the 5' termini of the proviruses with respect to viral RNA. Two enzyme sites near the 3'-end of the MAV-1-like provirus were not present in the putative AMV genome. We also examined a number of leukemic myeloblast clones for proviral content and cell-provirus integration sites. The presumptive AMV provirus was present in all the leukemic myeloblast clones regardless of the endogenous proviral content of the target cells or the AMV pseudotype used for conversion. Multiple cellular sites were suitable for integration of the putative AMV genome and the helper genomes. The proviral genomes were all integrated colinearly with respect to linear viral DNA.     

Discriminating antigen and non-antigen using proteome dissimilarity II: viral and fungal antigens

Immunogenicity arises via many synergistic mechanisms, yet the overall dissimilarity of pathogenic proteins versus the host proteome has been proposed as a key arbiter. We have previously explored this concept in relation to Bacterial antigens; here we extend our analysis to antigens of viral and fungal origin. Sets of known viral and fungal antigenic and non-antigenic protein sequences were compared to human and mouse proteomes. Both antigenic and non-antigenic sequences lacked human or mouse homologues. Observed distributions were compared using the non-parametric Mann-Whitney test. The statistical null hypothesis was accepted, indicating that antigen and non-antigens did not differ significantly. Likewise, we could not determine a threshold able meaningfully to separate non-antigen from antigen. We conclude that viral and fungal antigens cannot be predicted from pathogen genomes based solely on their dissimilarity to mammalian genomes.     

A novel repetitive DNA sequence in the genus Oryza

Summary Repetitive DNA sequences in the genus Oryza (rice) represent a large fraction of the nuclear DNA. The isolation and characterization of major repetitive DNA sequences will lead to a better understanding of rice genome organization and evolution. Here we report the characterization of a novel repetitive sequence, CC-1, from the CC genome. This repetitive sequence is present as long tandem arrays with a repeat unit 194 bp in length in the CC-diploid genome but 172 bp in length in the BBCC and CCDD tetraploid genomes. This repetitive sequence is also present, though at lower copy numbers, in the AA and BB genomes, but is absent in the EE and FF genomes. Hybridization experiments revealed considerable differences both in copy numbers and in restriction fragment patterns of CC-1 both between and within rice species. The results support the hypothesis that the CC genome is more closely related to the AA genome than to the BB genome, and most distantly related to the EE and FF genomes.     

Massive peptide sharing between viral and human proteomes

Thirty viral proteomes were examined for amino acid sequence similarity to the human proteome, and, in parallel, a control of 30 sets of human proteins was analyzed for internal human overlapping. We find that all of the analyzed 30 viral proteomes, independently of their structural or pathogenic characteristics, present a high number of pentapeptide overlaps to the human proteome. Among the examined viruses, human T-lymphotropic virus 1, Rubella virus, and hepatitis C virus present the highest number of viral overlaps to the human proteome. The widespread and ample distribution of viral amino acid sequences through the human proteome indicates that viral and human proteins are formed of common peptide backbone units and suggests a fluid compositional chimerism in phylogenetic entities canonically classified distantly as viruses and Homo sapiens. Importantly, the massive viral to human peptide overlapping calls into question the possibility of a direct causal association between virus-host sharing of amino acid sequences and incitement to autoimmune reactions through molecular recognition of common motifs.     

ZWA: Viral genome assembly and characterization hindrances from virus-host chimeric reads; a refining approach

Viral metagenomics, also known as virome studies, have yielded an unprecedented number of novel sequences, essential in recognizing and characterizing the etiological agent and the origin of emerging infectious diseases. Several tools and pipelines have been developed, to date, for the identification and assembly of viral genomes. Assembly pipelines often result in viral genomes contaminated with host genetic material, some of which are currently deposited into public databases. In the current report, we present a group of deposited sequences that encompass ribosomal RNA (rRNA) contamination. We highlight the detrimental role of chimeric next generation sequencing reads, between host rRNA sequences and viral sequences, in virus genome assembly and we present the hindrances these reads may pose to current methodologies. We have further developed a refining pipeline, the Zero Waste Algorithm (ZWA) that assists in the assembly of low abundance viral genomes. ZWA performs context-depended trimming of chimeric reads, precisely removing their rRNA moiety. These, otherwise discarded, reads were fed to the assembly pipeline and assisted in the construction of larger and cleaner contigs making a substantial impact on current assembly methodologies. ZWA pipeline may significantly enhance virus genome assembly from low abundance samples and virus metagenomics approaches in which a small number of reads determine genome quality and integrity.     

Vertebrate genomes code excess proteins with charge periodicity of 28 residues

All amino acid sequences derived from 248 prokaryotic genomes, 10 invertebrate genomes (plants and fungi) and 10 vertebrate genomes were analysed by the autocorrelation function of charge sequences. The analysis of the total amino acid sequences derived from the 268 biological genomes showed that a significant periodicity of 28 residues is observable for the vertebrate genomes, but not for the other genomes. When proteins with a charge periodicity of 28 residues (PCP28) were selected from the total proteomes, we found that PCP28 in fact exists in all proteomes, but the number of PCP28 is much larger for the vertebrate proteomes than for the other proteomes. Although excess PCP28 in the vertebrate proteomes are only poorly characterized, a detailed inspection of the databases suggests that most excess PCP28 are nuclear proteins.     

Errant processing and structural alterations of genomes present in a varicella-zoster virus vaccine.

Five minority populations of aberrant, varicella-zoster virus (VZV)-derived genomes were identified among the encapsidated DNAs obtained from the nuclear and cytoplasmic fractions of an in vitro infection initiated with a lyophilized sample of the BIKEN VZV vaccine (strain Oka). These were (i) VZV genomes, present within nuclear but not cytoplasmic viral capsids, which had been cleaved at a specific site within the short segment and which were, therefore, 3.15 megadaltons (approximately 4% of the VZV genome length) short of full length; (ii) highly deleted, repetitive VZV genomes which contained the errant cleavage site but not the usual VZV genome terminal sequences; (iii) VZV genomes into which multiples of 1 through 5 defective genome repeat units had been inserted into a homologous site; (iv) VZV genomes with additions of 0.1 or 0.18 megadaltons of DNA at both the terminal and internal ends of the short segment; and (v) VZV DNA which had lost the HindIII restriction site at map position 0.11.     

Recombination-Dependent Oligomerization of Human Papillomavirus Genomes upon Transient DNA Replication

We describe the extensive and progressive oligomerization of human papillomavirus (HPV) genomes after transfection into the U2OS cell line. The HPV genomic oligomers are extrachromosomal concatemeric molecules containing the viral genome in a head-to-tail orientation. The process of oligomerization does not depend on the topology of the input DNA, and it does not require any other viral factors besides replication proteins E1 and E2. We provide evidence that oligomerization of the HPV18 and HPV11 genomes involves homologous recombination. We also demonstrate oligomerization of the HPV18 and HPV11 genomes in SiHa, HeLa, and C-33 A cell lines and provide examples of oligomeric HPV genomes in clinical samples obtained from HPV-infected patients.     

How independent are the appearances of n-mers in different genomes?

MOTIVATION: Analysis of statistical properties of DNA sequences is important for evolutional biology as well as for DNA probe and PCR technologies. These technologies, in turn, can be used for organism identification, which implies applications in the diagnosis of infectious diseases, environmental studies, etc. RESULTS: We present results of the correlation analysis of distributions of the presence/absence of short nucleotide subsequences of different length ('n-mers', n = 5-20) in more than 1500 microbial and virus genomes, together with five genomes of multicellular organisms (including human). We calculate whether a given n-mer is present or absent (frequency of presence) in a given genome, which is not the usually calculated number of appearances of n-mers in one or more genomes (frequency of appearance). For organisms that are not close relatives of each other, the presence/absence of different 7-20mers in their genomes are not correlated. For close biological relatives, some correlation of the presence of n-mers in this range appears, but is not as strong as expected. Suppressed correlations among the n-mers present in different genomes leads to the possibility of using random sets of n-mers (with appropriately chosen n) to discriminate genomes of different organisms and possibly individual genomes of the same species including human with a low probability of error.     

Free and integrated recombinant murine leukemia virus DNAs appear in preleukemic thymuses of AKR/J mice.

We studied the appearance and structure of murine leukemia viral genomes in preleukemic AKR/J mice by Southern hybridization. Up to an average of one to two copies per thymocyte of unintegrated murine leukemia virus DNA appears in the thymuses of preleukemic mice beginning at 4 to 5 months of age and disappears in leukemic thymuses. The free viral genomes are absent in the spleens, livers, and brains of preleukemic mice. Using a series of ecotropic and nonecotropic murine leukemia virus hybridization probes, we showed that the unintegrated viral genomes are structurally analogous to those of recombinant mink cell focus-forming viruses that appear as proviruses in leukemic AKR thymocytes, suggesting that these free viral DNAs are the direct precursors to the leukemia-specific proviruses. The mosaic of ecotropic and nonecotropic sequences within these unintegrated viral DNAs varies from one preleukemic thymus to another but often appears structurally homogeneous within individual thymuses, indicating that often each thymus was being infected by a unique mink cell focus-forming virus. Analysis of high-molecular-weight DNA shows that recombinant proviruses reside in the chromosomal DNA of thymocytes within the preleukemic thymus, with the number rising to an average of several copies per thymocyte, but we do not detect any preferred integration sites. These results suggest that, in general, before the development of thymic leukemias in AKR mice there is a massive infection by a unique mink cell focus-forming virus which then integrates into many different sites of individual thymocytes, one of which grows out to become a tumor.     

Many human endogenous retrovirus K (HERV-K) proviruses are unique to humans.

BACKGROUND: Endogenous retroviruses contribute to the evolution of the host genome and can be associated with disease. Human endogenous retrovirus K (HERV-K) is related to the mouse mammary tumor virus and is present in the genomes of humans, apes and cercopithecoids (Old World monkeys). It is unknown how long ago in primate evolution the full-length HERV-K proviruses that are in the human genome today were formed. RESULTS: Ten full-length HERV-K proviruses were cloned from the human genome. Using provirus-specific probes, eight of the ten were found to be present in a genetically diverse set of humans but not in other extant hominoids. Intact preintegration sites for each of these eight proviruses were present in the apes. A ninth provirus was detected in the human, chimpanzee, bonobo and gorilla genomes, but not in the orang-utan genome. The tenth was found only in humans, chimpanzees and bonobos. Complete sequencing of six of the human-specific proviruses showed that full-length open reading frames for the retroviral protein precursors Gag-Pro-Pol or Env were each present in multiple proviruses. CONCLUSIONS: At least eight full-length HERV-K genomes that are in the human germline today integrated after humans diverged from chimpanzees. All of the viral open reading frames and cis-acting sequences necessary for HERV-K replication must have been intact during the recent time when these proviruses formed. Multiple full-length open reading frames for all HERV-K proteins are present in the human genome today.     

Comprehensive genome analysis of 203 genomes provides structural genomics with new insights into protein family space

We present an analysis of 203 completed genomes in the Gene3D resource (including 17 eukaryotes), which demonstrates that the number of protein families is continually expanding over time and that singleton-sequences appear to be an intrinsic part of the genomes. A significant proportion of the proteomes can be assigned to fewer than 6000 well-characterized domain families with the remaining domain-like regions belonging to a much larger number of small uncharacterized families that are largely species specific. Our comprehensive domain annotation of 203 genomes enables us to provide more accurate estimates of the number of multi-domain proteins found in the three kingdoms of life than previous calculations. We find that 67% of eukaryotic sequences are multi-domain compared with 56% of sequences in prokaryotes. By measuring the domain coverage of genome sequences, we show that the structural genomics initiatives should aim to provide structures for less than a thousand structurally uncharacterized Pfam families to achieve reasonable structural annotation of the genomes. However, in large families, additional structures should be determined as these would reveal more about the evolution of the family and enable a greater understanding of how function evolves.     

Regulation of expression and chromosomal subunit conformation of avian retrovirus genomes.

We have investigated the copy number, chromosomal subunit conformation and regulation of expression of integrated avian retrovirus genomes. Our results indicate that there are approximately two copies of the endogenous viral genomes (RAV-O) per haploid cell genome in uninfected chick embryo fibroblasts (CEF) and red blood cells (RBC). The copy number and subunit conformation (as measured by DNAasel sensitivity) of the RAV-O genomes are independent of the level of expression of these viral DNA sequences. In cells isolated from embryos of the V+, gs-chf- and gs+chf+ phenotypes, approximately one of the two viral genomes is in a DNAase l-sensitive conformation. Upon infection with an exogenous Rous sarcoma virus (PR-RSV-C), one new viral genome is integrated per haploid CEF genome. The newly integrated RSV genome is completely sensitive to DNAase l, and the subunit conformation of the endogenous viral genomes is not altered by the integration of additional exogenous proviruses. Both the endogenous and newly integrated exogenous viral genomes are present in "nu-body" structures, and the selective sensitivity of these proviral DNA sequences to DNAase l is maintained in isolated nucleosomes. Our experiments revealing the DNAase l sensitivity of one of the two RAV-O genomes in gs-chf-CEF led us to reexamine the level of viral specific RNA in CEF of various GS genotypes. We find that GS/GS CEF contain approximately 100 copies of viral RNA per cell, gs/gs CEF contain no detectable viral RNA, and the heterozygote GS/gs CEF contain approximately 50 copies of viral specific RNA per cell. These results suggest that the GS gene controls production of RAV-O RNA sequences in CEF in a "cis" fashion. In RBCs, however, the expression of the RAV-O genome is independent of the GS gene, with both GS/GS and gs/gs RBCs containing roughly equivalent amounts of viral specific RNA. Our results suggest that the chromosomal structure of the endogenous viral genes is independent of the GS gene, and that the GS gene is cis-acting and tissue-specific.     

Eucaryotic transposable genetic elements with inverted terminal repeats

DNA carrying inverted repeats was tested for transposition within the Drosophila genome. Five Bam HI segments containing related inverted repeats were isolated from D. melanogaster and analyzed by electron microscopy and restriction mapping. Southern blot experiments using single-copy flanking sequences as probes allowed the study of DNA arrangements at specific sites in the genomes of five closely related strains. We found that in some genomes the sequences with inverted repeats were present at a particular site, whereas in other genomes they were absent from this site. These results indicated that three of the sequences are transposable genetic elements. In one case we have purified the two corresponding DNA segments, with and without the sequence containing inverted repeats, thereby confirming the mobility of this sequence. These DNA elements were found to be distinct in two ways from copia and others previously described: first, they contain inverted terminal repeats, and second, they have a more heterogeneous construction.     

Higher Rates of Processed Pseudogene Acquisition in Humans and Three Great Apes Revealed by Long-Read Assemblies

LINE-1-mediated retrotransposition of protein-coding mRNAs is an active process in modern humans for both germline and somatic genomes. Prior works that surveyed human data mostly relied on detecting discordant mappings of paired-end short reads, or exon junctions contained in short reads. Moreover, there have been few genome-wide comparisons between gene retrocopies in great apes and humans. In this study, we introduced a more sensitive and accurate method to identify processed pseudogenes. Our method utilizes long-read assemblies, and more importantly, is able to provide full-length retrocopy sequences as well as flanking regions which are missed by short-read based methods. From 22 human individuals, we pinpointed 40 processed pseudogenes that are not present in the human reference genome GRCh38 and identified 17 pseudogenes that are in GRCh38 but absent from some input individuals. This represents a significantly higher discovery rate than previous reports (39 pseudogenes not in the reference genome out of 939 individuals). We also provided an overview of lineage-specific retrocopies in chimpanzee, gorilla, and orangutan genomes.