Discovery of Novel dsRNA Viral Sequences by In Silico Cloning and Implications for Viral Diversity, Host Range and Evolution

Genome sequence of viruses can contribute greatly to the study of viral evolution, diversity and the interaction between viruses and hosts. Traditional molecular cloning methods for obtaining RNA viral genomes are time-consuming and often difficult because many viruses occur in extremely low titers. DsRNA viruses in the families, Partitiviridae, Totiviridae, Endornaviridae, Chrysoviridae, and other related unclassified dsRNA viruses are generally associated with symptomless or persistent infections of their hosts. These characteristics indicate that samples or materials derived from eukaryotic organisms used to construct cDNA libraries and EST sequencing might carry these viruses, which were not easily detected by the researchers. Therefore, the EST databases may include numerous unknown viral sequences. In this study, we performed in silico cloning, a procedure for obtaining full or partial cDNA sequence of a gene by bioinformatics analysis, using known dsRNA viral sequences as queries to search against NCBI Expressed Sequence Tag (EST) database. From this analysis, we obtained 119 novel virus-like sequences related to members of the families, Endornaviridae, Chrysoviridae, Partitiviridae, and Totiviridae. Many of them were identified in cDNA libraries of eukaryotic lineages, which were not known to be hosts for these viruses. Furthermore, comprehensive phylogenetic analysis of these newly discovered virus-like sequences with known dsRNA viruses revealed that these dsRNA viruses may have co-evolved with respective host supergroups over a long evolutionary time while potential horizontal transmissions of viruses between different host supergroups also is possible. We also found that some of the plant partitiviruses may have originated from fungal viruses by horizontal transmissions. These findings extend our knowledge of the diversity and possible host range of dsRNA viruses and offer insight into the origin and evolution of relevant viruses with their hosts.     

A Survey of Overlooked Viral Infections in Biological Experiment Systems

It is commonly accepted that there are many unknown viruses on the planet. For the known viruses, do we know their prevalence, even in our experimental systems? Here we report a virus survey using recently published small (s)RNA sequencing datasets. The sRNA reads were assembled and contigs were screened for virus homologues against the NCBI nucleotide (nt) database using the BLASTn program. To our surprise, approximately 30% (28 out of 94) of publications had highly scored viral sequences in their datasets. Among them, only two publications reported virus infections. Though viral vectors were used in some of the publications, virus sequences without any identifiable source appeared in more than 20 publications. By determining the distributions of viral reads and the antiviral RNA interference (RNAi) pathways using the sRNA profiles, we showed evidence that many of the viruses identified were indeed infecting and generated host RNAi responses. As virus infections affect many aspects of host molecular biology and metabolism, the presence and impact of viruses needs to be actively investigated in experimental systems.     

DING proteins are from Pseudomonas

DING proteins have been described as animal and plant proteins with potential biomineralisation, receptor or signalling roles that have been characterised by an N-terminal DINGGG-sequence. However, these sequences have only ever been identified as either N-terminal peptides or partial cDNA sequences, and have yet to be detected in any of the many genomic animal and plant genomes now available. Microbial relatives of the DING proteins have been described, which appear to be periplasmic phosphate-binding proteins. Recently, full-length Pseudomonas aeruginosa UCBPP-PA14 and Hypericum perforatum genes have been sequenced that show high homology to the published DING protein N-terminal sequences, and small peptides previously identified in conjunction with the peptide sequencing of DING proteins can also be mapped to regions across these full-length sequences. Searching with these sequences identifies other plant and animal cDNA fragments in the public nucleotide databases, and, additionally, an unordered rat genomic contig that contains a DING-like sequence on a small fragment. Analysing the codon usage of these DNA sequences identifies all of these sequences as of Pseudomonas origin, suggesting that DING proteins do not exist in eukaryotes, but instead are potentially due to microbial contamination or infection.     

Previously unknown and highly divergent ssDNA viruses populate the oceans

Single-stranded DNA (ssDNA) viruses are economically important pathogens of plants and animals, and are widespread in oceans; yet, the diversity and evolutionary relationships among marine ssDNA viruses remain largely unknown. Here we present the results from a metagenomic study of composite samples from temperate (Saanich Inlet, 11 samples; Strait of Georgia, 85 samples) and subtropical (46 samples, Gulf of Mexico) seawater. Most sequences (84%) had no evident similarity to sequenced viruses. In total, 608 putative complete genomes of ssDNA viruses were assembled, almost doubling the number of ssDNA viral genomes in databases. These comprised 129 genetically distinct groups, each represented by at least one complete genome that had no recognizable similarity to each other or to other virus sequences. Given that the seven recognized families of ssDNA viruses have considerable sequence homology within them, this suggests that many of these genetic groups may represent new viral families. Moreover, nearly 70% of the sequences were similar to one of these genomes, indicating that most of the sequences could be assigned to a genetically distinct group. Most sequences fell within 11 well-defined gene groups, each sharing a common gene. Some of these encoded putative replication and coat proteins that had similarity to sequences from viruses infecting eukaryotes, suggesting that these were likely from viruses infecting eukaryotic phytoplankton and zooplankton.     

Evolution of the HIV-1 nef gene in HLA-B*57 Positive Elite Suppressors

During pilot studies to investigate the presence of viral RNA of xenotropic murine leukemia virus (MLV)-related virus (XMRV) infection in sera from chronic fatigue syndrome (CFS) patients in Japan, a positive band was frequently detected at the expected product size in negative control samples when detecting a partial gag region of XMRV using a one-step RT-PCR kit. We suspected that the kit itself might have been contaminated with small traces of endogenous MLV genome or XMRV and attempted to evaluate the quality of the kit in two independent laboratories. We purchased four one-step RT-PCR kits from Invitrogen, TaKaRa, Promega and QIAGEN in Japan. To amplify the partial gag gene of XMRV or other MLV-related viruses, primer sets (419F and 1154R, and GAG-I-F and GAG-I-R) which have been widely used in XMRV studies were employed. The nucleotide sequences of the amplicons were determined and compared with deposited sequences of a polytropic endogenous MLV (PmERV), XMRV and endogenous MLV-related viruses derived from CFS patients. We found that the enzyme mixtures of the one-step RT-PCR kit from Invitrogen were contaminated with RNA derived from PmERV. The nucleotide sequence of a partial gag region of the contaminant amplified by RT-PCR was nearly identical (99.4% identity) to a PmERV on chromosome 7 and highly similar (96.9 to 97.6%) to recently identified MLV-like viruses derived from CFS patients. We also determined the nucleotide sequence of a partial env region of the contaminant and found that it was almost identical (99.6%) to the PmERV. In the investigation of XMRV infection in patients of CFS and prostate cancer, researchers should prudently evaluate the test kits for the presence of endogenous MLV as well as XMRV genomes prior to PCR and RT-PCR tests.     

Patterns of variation among distinct alleles of the Flag silk gene from Nephila clavipes

Spider silk proteins and their genes are very attractive to researchers in a wide range of disciplines because they permit linking many levels of organization. However, hypotheses of silk gene evolution have been built primarily upon single sequences of each gene each species, and little is known about allelic variation within a species. Silk genes are known for their repeat structure with high levels of homogenization of nucleotide and amino acid sequence among repeated units. One common explanation for this homogeneity is gene convergence. To test this model, we sequenced multiple alleles of one intron-exon segment from the Flag gene from four populations of the spider Nephila clavipes and compared the new sequences to a published sequence. Our analysis revealed very high levels of heterozygosity in this gene, with no pattern of population differentiation. There was no evidence of gene convergence within any of these alleles, with high levels of nucleotide and amino acid substitution among the repeating motifs. Our data suggest that minimally, there is relaxed selection on mutations in this gene and that there may actually be positive selection for heterozygosity.     

Sequence analysis of members of the human Ig VH4 gene family derived from a single VH locus. Identification of novel germ-line members.

We have generated a mouse x human heterohybridoma that contains a single copy of chromosome 14 and, thus, a haploid set of Ig VH genes. This cell line was used to investigate the germ-line content and nucleotide sequences of members of the VH4 gene family in a polymerase chain reaction-based approach. The analysis of 58 full-length sequences revealed the presence of 12 different germ-line VH4 genes, each of which is potentially functional. These germ-line VH4 genes were compared with the nucleotide sequences of published VH4 genes. Three VH4 genes were 100% identical to previously published sequences and belong to a group of VH4 genes that are strongly conserved and highly prevalent in the human population. Three VH4 genes in our collection displayed greater than 99.3% sequence identity with reported germ-line VH4 sequences and likely represent allelic counterparts of these genes. Six genes displayed less than 97.2% sequence identity with published VH4 genes and were identified as novel members of the human VH4 gene family or more distantly related alleles of known VH4 genes. Collectively, these data suggest that, overall, the human VH4 gene family may be more diverse than hitherto assumed, whereas a number of individual members are nonpolymorphic and extremely well conserved.     

Large-scale search for genes on which positive selection may operate

We conducted a systematic search for the candidate genes on which positive selection may operate, on the premise that for such genes the number of nonsynonymous substitution is expected to be larger than that of synonymous substitutions when the nucleotide sequences of genes under investigation are compared with each other. By obtaining 3,595 groups of homologous sequences from the DDBJ, EMBL, and GenBank DNA sequence databases, we found that 17 gene groups can be the candidates for the genes on which positive selection may operate. Thus, such genes are found to occupy only about 0.5% of the vast number of gene groups so far available. Interestingly enough 9 out of the 17 gene groups were the surface antigens of parasites or viruses.      

Polymerase errors accumulating during natural evolution of the glycoprotein gene of vesicular stomatitis virus Indiana serotype isolates

We report the entire glycoprotein (G) gene nucleotide sequences of 26 vesicular stomatitis virus Indiana serotype (VSV IND) type 1 isolates from North and Central America. These sequences are also compared with partial G gene sequences of VSV IND type 2 (Cocal) and type 3 (Alagoas) viruses and the complete G gene sequences of the more distantly related VSV New Jersey (NJ) and Chandipura viruses. Phylogenetic analysis of the G gene sequences by maximum parsimony revealed four major lineages or subtypes within the classical VSV IND (type 1) viruses, each with a distinct geographic distribution. A high degree of VSV genetic diversity was found in Central America, with several virus subtypes of both VSV IND and NJ serotypes existing in this mainly enzootic disease region. Nineteen percent sequence variation but no deletions or insertions were evident within the 5' noncoding and the coding regions of the VSV IND type 1 G genes. In addition to numerous base substitutions, the 3' noncoding regions of these viruses also contained numerous base insertions and deletions. This resulted in striking variation in G gene sizes, with gene lengths ranging from 1,652 to 1,868 nucleotides. As the VSV IND type 1 subtypes have diverged from the common ancestor with the NJ subtypes, their G mRNAs have accumulated more 3' noncoding sequence inserts, ranging up to 303 nucleotides in length. These primarily consist of an imprecise reiteration of the sequence UUUUUAA, apparently generated by a unique polymerase stuttering error. Analysis of the deduced amino acid sequence differences among VSV IND type 1 viruses revealed numerous substitutions within defined antigenic epitopes, suggesting that immune selection may play a role in the evolution of these viruses.     

Abundant human DNA contamination identified in non-primate genome databases

During routine screens of the NCBI databases using human repetitive elements we discovered an unlikely level of nucleotide identity across a broad range of phyla. To ascertain whether databases containing DNA sequences, genome assemblies and trace archive reads were contaminated with human sequences, we performed an in depth search for sequences of human origin in non-human species. Using a primate specific SINE, AluY, we screened 2,749 non-primate public databases from NCBI, Ensembl, JGI, and UCSC and have found 492 to be contaminated with human sequence. These represent species ranging from bacteria (B. cereus) to plants (Z. mays) to fish (D. rerio) with examples found from most phyla. The identification of such extensive contamination of human sequence across databases and sequence types warrants caution among the sequencing community in future sequencing efforts, such as human re-sequencing. We discuss issues this may raise as well as present data that gives insight as to how this may be occurring.     

The smaller human VH gene families display remarkably little polymorphism.

We report the nucleotide sequence of 30 distinct human VH gene segments from the VHIV, VHV and VHVI gene families. When these sequences were compared to previously published sequences from these smaller human VH families a surprisingly low level of polymorphism was noted. Two VHIV gene segments from unrelated individuals were identical to two previously published VHIV sequences. Five VHV sequences were identical and seven VHVI gene segments were identical. Where differences were found between the sequences, allele specific oligonucleotide probes were used to verify the germline nature of the change and to test for segregation in several large kindreds. These data provide evidence that at least some human VH gene segments are remarkably stable.     

麻疹病毒S191疫苗株N基因稳定性研究及免疫细胞表位分析

研究麻疹病毒(Measles virus,MeV)疫苗株S191毒种和传代病毒核蛋白(Nucleoprotein,N)基因稳定性及其遗传与变异特点;对该序列一些重要位点的氨基酸进行比较,探讨其功能结构及生物学活性变化以及S191疫苗株的保护效果。利用RT-PCR方法扩增S191减毒株23、26、27、29、32、37不同代次N基因,测序进行比对分析。S191传代病毒N基因序列之间核苷酸同源性99.7%~100%,氨基酸同源性为99.6%~100%;S191株与7个疫苗株之间核苷酸序列同源性达99.1%~99.4%;S191和中国流行代表株序列同源性在95.0%~95.4%;S191与世界流行代表株同源性达94.7%~99.4%;S191疫苗株和中国流行代表株CHN93/7(H1a)的4个重要T细胞表位氨基酸保持一致。S191各传代病毒基因具有较高稳定性,该疫苗有一定的保护作用。     

Characterization of the env gene and long terminal repeat of molecularly cloned Friend mink cell focus-inducing virus DNA

The highly oncogenic erythroleukemia-inducing Friend mink cell focus-inducing (MCF) virus was molecularly cloned in phage lambda gtWES.lambda B, and the DNA sequences of the env gene and the long terminal repeat were determined. The nucleotide sequences of Friend MCF virus and Friend spleen focus-forming virus were quite homologous, supporting the hypothesis that Friend spleen focus-forming virus might be generated via Friend MCF virus from an ecotropic Friend virus mainly by some deletions. Despite their different pathogenicity, the nucleotide sequences of the env gene of Friend MCF virus and Moloney MCF virus were quite homologous, suggesting that the putative parent sequence for the generation of both MCF viruses and the recombinational mechanism for their generation might be the same. We compare the amino acid sequences in lymphoid leukemia-inducing ecotropic Moloney virus and Moloney MCF virus, and erythroblastic leukemia-inducing ecotropic Friend virus, Friend-MCF virus, and Friend spleen focus-forming virus. The Friend MCF virus long terminal repeat was found to be 550 base pairs long. This contained two copies of the 39-base-pair tandem repeat, whereas the spleen focus-forming virus genome contained a single copy of the same sequence.     

A complete analysis of HA and NA genes of influenza A viruses

Background

More and more nucleotide sequences of type A influenza virus are available in public databases. Although these sequences have been the focus of many molecular epidemiological and phylogenetic analyses, most studies only deal with a few representative sequences. In this paper, we present a complete analysis of all Haemagglutinin (HA) and Neuraminidase (NA) gene sequences available to allow large scale analyses of the evolution and epidemiology of type A influenza.

Methodology/Principal Findings

This paper describes an analysis and complete classification of all HA and NA gene sequences available in public databases using multivariate and phylogenetic methods.

Conclusions/Significance

We analyzed 18975 HA sequences and divided them into 280 subgroups according to multivariate and phylogenetic analyses. Similarly, we divided 11362 NA sequences into 202 subgroups. Compared to previous analyses, this work is more detailed and comprehensive, especially for the bigger datasets. Therefore, it can be used to show the full and complex phylogenetic diversity and provides a framework for studying the molecular evolution and epidemiology of type A influenza virus. For more than 85% of type A influenza HA and NA sequences into GenBank, they are categorized in one unambiguous and unique group. Therefore, our results are a kind of genetic and phylogenetic annotation for influenza HA and NA sequences. In addition, sequences of swine influenza viruses come from 56 HA and 45 NA subgroups. Most of these subgroups also include viruses from other hosts indicating cross species transmission of the viruses between pigs and other hosts. Furthermore, the phylogenetic diversity of swine influenza viruses from Eurasia is greater than that of North American strains and both of them are becoming more diverse. Apart from viruses from human, pigs, birds and horses, viruses from other species show very low phylogenetic diversity. This might indicate that viruses have not become established in these species. Based on current evidence, there is no simple pattern of inter-hemisphere transmission of avian influenza viruses and it appears to happen sporadically. However, for H6 subtype avian influenza viruses, such transmissions might have happened very frequently and multiple and bidirectional transmission events might exist.     

Rabbit Ig kappa 1b6 gene structure

Previous studies employing Southern blot analyses have detected multiple kappa-homologous sequences within EcoRI-digested DNA isolated from kappa 1b6 homozygous rabbits and kappa 1b6 L chain secreting RMH H158 cell line. These results are very unexpected because the published partial protein sequence for the kappa 1b6 C region is incompatible with an EcoRI restriction endonuclease recognition sequence at the nucleotide level for this allotype. To determine their identity, the kappa-homologous sequences were isolated from DNA extracted from a kappa 1b6 L chain secreting RMH H158 cell line by molecular cloning. Structural analyses demonstrated these sequences to contain genetic information encoding the majority of the kappa 1b6 L chain gene locus. The protein sequence deduced from the kappa 1b6 C region gene was shown to differ from the published partial kappa 1b6 C region protein sequence at five amino acid positions. One of these differences results in a glycine to serine interchange that introduces an EcoRI restriction endonuclease recognition site within the kappa 1b6 C region gene. Subsequent genomic Southern blot analyses confirmed this structural assignment. Based on these data, the EcoRI-sensitive kappa-homologous fragments present within the genomes of the RMH H158 cell line and kappa 1b6 homozygous rabbits represent the nominal kappa 1 gene and not an alternative kappa isotype or kappa pseudogene. Rabbit Ig kappa 1 allelic nucleotide sequence homology comparisons have shown the isolated kappa 1b6 J-C gene locus to display common structural features previously identified in other kappa 1 alleles.     

Structural analysis of a Lotus japonicus genome. III. Sequence features and mapping of sixty-two TAC clones which cover the 6.7 Mb regions of the genome.

A total of sixty-two clones were selected from a TAC (transformation-competent artificial chromosome) genomic library of the Lotus japonicus accession MG-20 based on the sequence information of expressed sequence tags (ESTs), cDNA and gene information, and their nucleotide sequences were determined. The length of the sequenced regions in this study is 6,682,189 bp, and the total length of the regions sequenced so far is 18,711,484 bp together with the nucleotide sequences of 121 TAC clones previously reported. By comparison with the sequences in protein and EST databases and analysis with computer programs for gene modeling, a total of 573 potential protein-coding genes with known or predicted functions, 91 gene segments and 272 pseudogenes were identified in the newly sequenced regions. Each of the sequenced clones was localized onto the linkage map of two accessions of L. japonicus, Gifu B-129 and Miyakojima MG-20, using simple sequence repeat length polymorphism (SSLP) or derived cleaved amplified polymorphic sequence (dCAPS) markers generated based on the nucleotide sequences of the clones. The sequence data, gene information and mapping information are available through the World Wide Web at http://www.kazusa.or.jp/lotus/.     

Polymorphism of the IGHA gene in sheep

Genetic variation in immunoglobulin A, the most abundant immunoglobulin in mammalian cells, has not been reported in ruminants. In this study, variation in the immunoglobulin heavy alpha chain constant gene (IGHA) of sheep was investigated by amplification of a fragment that included the hinge coding sequence, followed by single-strand conformational polymorphism (SSCP) analysis and DNA sequencing. Three novel sequences, each characterized by unique SSCP banding patterns, were identified. One or two sequences were detected in individual sheep and all the sequences identified shared high homology to the published ovine and bovine IGHA sequences, suggesting that these sequences represent allelic variants of the IGHA gene in sheep. Sequence alignment showed that these sequences differed mainly in the 3′ end of exon 1 and in the coding sequence of the hinge region. There was either a deletion or an insertion of two codons in the hinge coding region in these allelic variants. Codon usage in the hinge coding region was quite different from that in the non-hinge coding regions of the gene, suggesting different evolution of the IGHA hinge sequence. Three novel amino acid sequences of ovine IGHA were also predicted, and variation in these sequences might not only affect antigen recognition but also susceptibility to cleavage by bacterial or parasitic proteases. Nucleotide sequence data reported in this paper have been submitted to the NCBI GenBank nucleotide sequence database and have been assigned the accession nos. AY956424–AY956426.