Complete sequence and comparative analysis of the genome of herpes B virus (Cercopithecine herpesvirus 1) from a rhesus monkey

The complete DNA sequence of herpes B virus (Cercopithecine herpesvirus 1) strain E2490, isolated from a rhesus macaque, was determined. The total genome length is 156,789 bp, with 74.5% G+C composition and overall genome organization characteristic of alphaherpesviruses. The first and last residues of the genome were defined by sequencing the cloned genomic termini. There were six origins of DNA replication in the genome due to tandem duplication of both oriL and oriS regions. Seventy-four genes were identified, and sequence homology to proteins known in herpes simplex viruses (HSVs) was observed in all cases but one. The degree of amino acid identity between B virus and HSV proteins ranged from 26.6% (US5) to 87.7% (US15). Unexpectedly, B virus lacked a homolog of the HSV gamma(1)34.5 gene, which encodes a neurovirulence factor. Absence of this gene was verified in two low-passage clinical isolates derived from a rhesus macaque and a zoonotically infected human. This finding suggests that B virus most likely utilizes mechanisms distinct from those of HSV to sustain efficient replication in neuronal cells. Despite the considerable differences in G+C content of the macaque and B virus genes (51% and 74.2%, respectively), codons used by B virus are optimal for the tRNA population of macaque cells. Complete sequence of the B virus genome will certainly facilitate identification of the genetic basis and possible molecular mechanisms of enhanced B virus neurovirulence in humans, which results in an 80% mortality rate following zoonotic infection.     

Proteins Specified by bovine herpesvirus 1 (infectious bovine rhinotracheitis virus)

An electrophoretic analysis of radioactively labeled, purified, "empty" and DNA-containing infectious bovine rhinotracheitis virions revealed the presence of 25 to 33 structural (virion) polypeptides. A total of 11 of these polypeptides could be labeled with [3H]glucosamine and were identified as glycoproteins. In addition to the 25 structural polypeptides, infectious bovine rhinotracheitis virus infected cells also contained at least 15 nonstructural (nonvirion) polypeptides that were not present in purified virions. Expression of the viral polypeptides in infected cells was controlled temporally. Thus, most viral polypeptides could be categorized as "alpha" (immediate early), "beta" (early), or "gamma" (late) on the basis of their order of appearance in infected cells and whether their syntheses were dependent upon prior viral protein or DNA synthesis. None of the glycoproteins belongs to the alpha class, although at least one (GVP11) was synthesized in the absence of viral DNA synthesis. Serum from a cow in which infectious bovine rhinotracheitis virus lesions were reactivated by dexamethasone precipitated both structural and nonstructural polypeptides.     

Features of DNA fragments obtained by random amplified polymorphic DNA (RAPD) assays

Random amplified polymorphic DNA (RAPD) fragments were prepared from samples of Calonectris diomedea (Cory's shearwater, Aves) and Haemonchus contortus (Nematoda) DNA by polymerase chain reaction (PCR) using decamers containing two restriction enzyme sites as primers. Six of 19 studied RAPD fragments probably originated from traces of commensal microorganisms. Many rearranged fragments, absent in the original genomic DNA, were synthesized and amplified during the processing of all the DNA samples, indicating that interactions occur within and between strands during the annealing step of PCR. The model of interactions between molecular species during DNA amplification with a single arbitrary oligonucleotide primer was modified to include nested primer annealing and interactions within and between strands. The presence of these artefacts in the final RAPD have a major effect on the interpretation of polymorphism studies.     

Characteristics of the tomato nuclear genome as determined by sequencing undermethylated EcoRI digested fragments

A collection of 9,990 single-pass nuclear genomic sequences, corresponding to 5 Mb of tomato DNA, were obtained using methylation filtration (MF) strategy and reduced to 7,053 unique undermethylated genomic islands (UGIs) distributed as follows: (1) 59% non-coding sequences, (2) 28% coding sequences, (3) 12% transposons—96% of which are class I retroelements, and (4) 1% organellar sequences integrated into the nuclear genome over the past approximately 100 million years. A more detailed analysis of coding UGIs indicates that the unmethylated portion of tomato genes extends as far as 676 bp upstream and 766 bp downstream of coding regions with an average of 174 and 171 bp, respectively. Based on the analysis of the UGI copy distribution, the undermethylated portion of the tomato genome is determined to account for the majority of the unmethylated genes in the genome and is estimated to constitute 61±15 Mb of DNA (~5% of the entire genome)—which is significantly less than the 220 Mb estimated for gene-rich euchromatic arms of the tomato genome. This result indicates that, while most genes reside in the euchromatin, a significant portion of euchromatin is methylated in the intergenic spacer regions. Implications of the results for sequencing the genome of tomato and other solanaceous species are discussed.     

Potential of genomic selection for improvement of resistance to ostreid herpesvirus in Pacific oyster (Crassostrea gigas)

In genomic selection (GS), genome-wide SNP markers are used to generate genomic estimated breeding values for selection candidates. The application of GS in shellfish looks promising and has the potential to help in dealing with one of the main issues currently affecting Pacific oyster production worldwide, which is the ‘summer mortality syndrome’. This causes periodic mass mortality in farms worldwide and has mainly been attributed to a specific variant of the ostreid herpesvirus (OsHV-1). In the current study, we evaluated the potential of genomic selection for host resistance to OsHV-1 in Pacific oysters, and compared it with pedigree-based approaches. An OsHV-1 disease challenge was performed using an immersion-based virus exposure treatment for oysters for 7 days. A total of 768 samples were genotyped using the medium-density SNP array for oysters. A GWAS was performed for the survival trait using a GBLUP approach in blupf 90 software. Heritability ranged from 0.25 ± 0.05 to 0.37 ± 0.05 (mean ± SE) based on pedigree and genomic information respectively. Genomic prediction was more accurate than pedigree prediction, and SNP density reduction had little impact on prediction accuracy until marker densities dropped below approximately 500 SNPs. This demonstrates the potential for GS in Pacific oyster breeding programmes, and importantly, demonstrates that a low number of SNPs might suffice to obtain accurate genomic estimated breeding values, thus potentially making the implementation of GS more cost effective.     

Analysis and characterization of the complete genome of tupaia (tree shrew) herpesvirus

The tupaia herpesvirus (THV) was isolated from spontaneously degenerating tissue cultures of malignant lymphoma, lung, and spleen cell cultures of tree shrews (Tupaia spp.). The determination of the complete nucleotide sequence of the THV strain 2 genome resulted in a 195,857-bp-long, linear DNA molecule with a G+C content of 66.5%. The terminal regions of the THV genome and the loci of conserved viral genes were found to be G+C richer. Furthermore, no large repetitive DNA sequences could be identified. This is in agreement with the previous classification of THV as the prototype species of herpesvirus genome group F. The search for potential coding regions resulted in the identification of 158 open reading frames (ORFs) regularly distributed on both DNA strands. Seventy-six out of the 158 ORFs code for proteins that are significantly homologous to known herpesvirus proteins. The highest homologies found were to primate and rodent cytomegaloviruses. Biological properties, protein homologies, the arrangement of conserved viral genes, and phylogenetic analysis revealed that THV is a member of the subfamily Betaherpesvirinae. The evolutionary lineages of THV and the cytomegaloviruses seem to have branched off from a common ancestor. In addition, it was found that the arrangements of conserved genes of THV and murine cytomegalovirus strain Smith, both of which are not able to form genomic isomers, are colinear with two different human cytomegalovirus (HCMV) strain AD169 genomic isomers that differ from each other in the orientation of the long unique region. The biological properties and the high degree of relatedness of THV to the mammalian cytomegaloviruses allow the consideration of THV as a model system for investigation of HCMV pathogenicity.     

Systemic herpes-like virus in catfish Ictalurus melas (Italy) differs from Ictalurid herpesvirus 1 (North America)

A herpesvirus was isolated during 2 occurrences of mass mortality among adult catfish Ictalurus melas raised in different farms in northern Italy. The agent replicated in the channel catfish ovary (CCO) cell line from channel catfish I. punctatus, inducing a cytopathic effect similar to that caused by Ictalurid herpesvirus 1 (also referred to as channel catfish herpesvirus, CCV). The new herpesvirus, designated I. melas herpesvirus (IcmHV) did not react with polyclonal rabbit or monoclonal antibodies directed to CCV in either neutralization or indirect immunofluorescence assays. The virions of IcmHV possessed a hexagonal nucleocapsid of 107 nm in diameter surrounded by an envelope with a diameter of 227 nm (n = 20) typical for members of the family Herpesviridae. Virions of IcmHV purified from infected CCO cells contained 17 polypeptides ranging in size from 17.5 to 175 kDa and most differed in molecular weight from those found for CCV. The IcmHV was also distinct from CCV when compared by restriction fragment length polymorphisms (RFLP) of genomic DNA following digestions with the endonucleases Kpn I and Sac I. Lastly, the virulence of IcmHV for channel catfish fry and juveniles, respectively, was demonstrated by experimental infections induced by bath exposure or intraperitoneal injection that resulted in 78 to 96% cumulative mortality in groups of exposed fish. Preventing the introduction of this agent into geographic regions where significant channel catfish production occurs should be a high priority.     

Detection of ostreid herpesvirus-1 (OsHV-1) by PCR using a rapid and simple method of DNA extraction from oyster larvae

A DNA extraction procedure was developed for the detection of ostreid herpesvirus-1 (OsHV-1) using the polymerase chain reaction (PCR) in oyster larvae. The DNA extraction procedure developed was tested on 8 larval samples. Abnormal nuclei with characteristic features associated with OsHV-1 infections were only observed in samples in which the viral DNA was detected by PCR. A previously described competitive PCR method was applied to detect inhibition during PCR reactions. The results show that the method can be used on small amounts of oyster larvae (3 mg) for the detection of OsHV-1 DNA by PCR.     

Biological characterization and next-generation genome sequencing of the unclassified Cotia virus SPAn232 (Poxviridae)

Cotia virus (COTV) SPAn232 was isolated in 1961 from sentinel mice at Cotia field station, São Paulo, Brazil. Attempts to classify COTV within a recognized genus of the Poxviridae have generated contradictory findings. Studies by different researchers suggested some similarity to myxoma virus and swinepox virus, whereas another investigation characterized COTV SPAn232 as a vaccinia virus strain. Because of the lack of consensus, we have conducted an independent biological and molecular characterization of COTV. Virus growth curves reached maximum yields at approximately 24 to 48 h and were accompanied by virus DNA replication and a characteristic early/late pattern of viral protein synthesis. Interestingly, COTV did not induce detectable cytopathic effects in BSC-40 cells until 4 days postinfection and generated viral plaques only after 8 days. We determined the complete genomic sequence of COTV by using a combination of the next-generation DNA sequencing technologies 454 and Illumina. A unique contiguous sequence of 185,139 bp containing 185 genes, including the 90 genes conserved in all chordopoxviruses, was obtained. COTV has an interesting panel of open reading frames (ORFs) related to the evasion of host defense, including two novel genes encoding C-C chemokine-like proteins, each present in duplicate copies. Phylogenetic analysis revealed the highest amino acid identity scores with Cervidpoxvirus, Capripoxvirus, Suipoxvirus, Leporipoxvirus, and Yatapoxvirus. However, COTV grouped as an independent branch within this clade, which clearly excluded its classification as an Orthopoxvirus. Therefore, our data suggest that COTV could represent a new poxvirus genus.     

Inhibition of human immunodeficiency virus 1 (HIV-1) by variant B of human herpesvirus 6 (HHV-6).

Four HHV-6 strains were initially isolated during attempts to observe HIV-1 replication in cultured primary lymphocytes from 48 patients with AIDS. HHV-6 DNA from each strain was extracted from primary cell cultures and amplified using specific primers in a nested polymerase chain reaction (PCR) assay. All HHV-6 strains were classified as B variants by submitting the PCR products to the digestion of two restriction enzymes (Hind III and Bgl II). Since in primary cultures, the appearance of HHV-6 cytopathic effect was followed by a progressive reduction of HIV-1 replication, we tried to reproduce the observed inhibition in vitro. Two HHV-6 strains, used throughout the experiments, showed their ability to suppress HIV-1 replication when the viruses co-infected CD4+T lymphocyte cultures. While the intrinsic mechanism of this finding still remains unclear, the inhibition of HIV-1 replication was observed only when a high multiplicity of infection (m.o.i.) of HHV-6 and a low m.o.i. of HIV-1 were used in dually infected cell cultures. By using a semiquantitative determination of HIV-1 cDNA by PCR, it appears that the inhibition begins in infected cell cultures and, once established, does not allow any further HIV-1 replication.     

Substrate specificity of three viral thymidine kinases (TK): vaccinia virus TK, feline herpesvirus TK, and canine herpesvirus TK

In search of novel suicide gene candidates we have cloned and characterized thymidine kinases from three viruses; vaccinia virus TK (VVTK), feline herpesvirus TK (FHV-TK), and canine herpesvirus TK (CHV-TK). Our studies showed that VVTK primarily is a thymidine kinase, with a substrate specificity mainly restricted to dThd and only minor affinity for dCyd. VVTK also is related closely to mammalian thymidine kinase 1 (TK1), with 66% identity and 75% general homology. Although CHV-TK and FHV-TK are sequence related to herpes simplex virus types 1 thymidine kinase (HSV1-TK), with 31% and 35% identity and a general similarity of 54%, the substrate specificity of these enzymes was restricted to dThd and thymidine analogs.     

Examination of host genome for the presence of integrated fragments of Solenopsis invicta virus 1

A series of oligonucleotide primer pairs covering the entire genome of Solenopsis invicta virus 1 (SINV-1) were used to probe the genome of its host, S. invicta, for integrated fragments of the viral genome. All of the oligonucleotide primer sets yielded amplicons of anticipated size from cDNA created from an RNA template from SINV-1. However, no corresponding amplification was observed when genomic DNA (from 32 colonies of S. invicta) was used as template for the PCR amplifications. Host DNA integrity was verified by amplification of an ant-specific gene, SiGSTS1. The representation of fire ant colonies included both social forms, monogyne and polygyne, and those infected and uninfected with SINV-1. Furthermore, no amplification was observed from genomic DNA from ant samples collected from Argentina or the US. Thus, it appears that SINV-1 genome integration, or a portion therein, has not likely occurred within the S. invicta host genome.     

Psittacid herpesvirus 1 and infectious laryngotracheitis virus: Comparative genome sequence analysis of two avian alphaherpesviruses

Psittacid herpesvirus 1 (PsHV-1) is the causative agent of Pacheco's disease, an acute, highly contagious, and potentially lethal respiratory herpesvirus infection in psittacine birds, while infectious laryngotracheitis virus (ILTV) is a highly contagious and economically significant avian herpesvirus which is responsible for an acute respiratory disease limited to galliform birds. The complete genome sequence of PsHV-1 has been determined and compared to the ILTV sequence, assembled from published data. The PsHV-1 and ILTV genomes exhibit similar structural characteristics and are 163,025 bp and 148,665 bp in length, respectively. The PsHV-1 genome contains 73 predicted open reading frames (ORFs), while the ILTV genome contains 77 predicted ORFs. Both genomes contain an inversion in the unique long region similar to that observed in pseudorabies virus. PsHV-1 is closely related to ILTV, and it is proposed that it be assigned to the Iltovirus genus. These two avian herpesviruses represent a phylogenetically unique clade of alphaherpesviruses that are distinct from the Marek's disease-like viruses (Mardivirus). The determination of the complete genomic nucleotide sequences of PsHV-1 and ILTV provides a tool for further comparative and functional analysis of this unique class of avian alphaherpesviruses.     

Whole genome sequencing of an ethnic Pathan (Pakhtun) from the north-west of Pakistan

Background

Pakistan covers a key geographic area in human history, being both part of the Indus River region that acted as one of the cradles of civilization and as a link between Western Eurasia and Eastern Asia. This region is inhabited by a number of distinct ethnic groups, the largest being the Punjabi, Pathan (Pakhtuns), Sindhi, and Baloch.

Results

We analyzed the first ethnic male Pathan genome by sequencing it to 29.7-fold coverage using the Illumina HiSeq2000 platform. A total of 3.8 million single nucleotide variations (SNVs) and 0.5 million small indels were identified by comparing with the human reference genome. Among the SNVs, 129,441 were novel, and 10,315 nonsynonymous SNVs were found in 5,344 genes. SNVs were annotated for health consequences and high risk diseases, as well as possible influences on drug efficacy. We confirmed that the Pathan genome presented here is representative of this ethnic group by comparing it to a panel of Central Asians from the HGDP-CEPH panels typed for ~650 k SNPs. The mtDNA (H2) and Y haplogroup (L1) of this individual were also typical of his geographic region of origin. Finally, we reconstruct the demographic history by PSMC, which highlights a recent increase in effective population size compatible with admixture between European and Asian lineages expected in this geographic region.

Conclusions

We present a whole-genome sequence and analyses of an ethnic Pathan from the north-west province of Pakistan. It is a useful resource to understand genetic variation and human migration across the whole Asian continent.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1290-1) contains supplementary material, which is available to authorized users.     

High-throughput sequencing of three Lemnoideae (duckweeds) chloroplast genomes from total DNA

Background

Chloroplast genomes provide a wealth of information for evolutionary and population genetic studies. Chloroplasts play a particularly important role in the adaption for aquatic plants because they float on water and their major surface is exposed continuously to sunlight. The subfamily of Lemnoideae represents such a collection of aquatic species that because of photosynthesis represents one of the fastest growing plant species on earth.

Methods

We sequenced the chloroplast genomes from three different genera of Lemnoideae, Spirodela polyrhiza, Wolffiella lingulata and Wolffia australiana by high-throughput DNA sequencing of genomic DNA using the SOLiD platform. Unfractionated total DNA contains high copies of plastid DNA so that sequences from the nucleus and mitochondria can easily be filtered computationally. Remaining sequence reads were assembled into contiguous sequences (contigs) using SOLiD software tools. Contigs were mapped to a reference genome of Lemna minor and gaps, selected by PCR, were sequenced on the ABI3730xl platform.

Conclusions

This combinatorial approach yielded whole genomic contiguous sequences in a cost-effective manner. Over 1,000-time coverage of chloroplast from total DNA were reached by the SOLiD platform in a single spot on a quadrant slide without purification. Comparative analysis indicated that the chloroplast genome was conserved in gene number and organization with respect to the reference genome of L. minor. However, higher nucleotide substitution, abundant deletions and insertions occurred in non-coding regions of these genomes, indicating a greater genomic dynamics than expected from the comparison of other related species in the Pooideae. Noticeably, there was no transition bias over transversion in Lemnoideae. The data should have immediate applications in evolutionary biology and plant taxonomy with increased resolution and statistical power.