Shiliang Wang

Widespread Lateral Gene Transfer from Intracellular Bacteria to Multicellular Eukaryotes

Although common among bacteria, lateral gene transfer—the movement of genes between distantly related organisms—is thought to occur only rarely between bacteria and multicellular eukaryotes. However, the presence of endosymbionts, such as Wolbachia pipientis, within some eukaryotic germlines may facilitate bacterial gene transfers to eukaryotic host genomes. We therefore examined host genomes for evidence of gene transfer events from Wolbachia bacteria to their hosts. We found and confirmed transfers into the genomes of four insect and four nematode species that range from nearly the entire Wolbachia genome (>1 megabase) to short (<500 base pairs) insertions. Potential Wolbachia-to-host transfers were also detected computationally in three additional sequenced insect genomes. We also show that some of these inserted Wolbachia genes are transcribed within eukaryotic cells lacking endosymbionts. Therefore, heritable lateral gene transfer occurs into eukaryotic hosts from their prokaryote symbionts, potentially providing a mechanism for acquisition of new genes and functions.

Draft Genome of the Filarial Nematode Parasite Brugia malayi

Parasitic nematodes that cause elephantiasis and river blindness threaten hundreds of millions of people in the developing world. We have sequenced the ∼90 megabase (Mb) genome of the human filarial parasite Brugia malayi and predict ∼11,500 protein coding genes in 71 Mb of robustly assembled sequence. Comparative analysis with the free-living, model nematode Caenorhabditis elegans revealed that, despite these genes having maintained little conservation of local synteny during ∼350 million years of evolution, they largely remain in linkage on chromosomal units. More than 100 conserved operons were identified. Analysis of the predicted proteome provides evidence for adaptations of B. malayi to niches in its human and vector hosts and insights into the molecular basis of a mutualistic relationship with its Wolbachia endosymbiont. These findings offer a foundation for rational drug design.

The Wolbachia Genome of Brugia malayi: Endosymbiont Evolution within a Human Pathogenic Nematode

Complete genome DNA sequence and analysis is presented for Wolbachia, the obligate alpha-proteobacterial endosymbiont required for fertility and survival of the human filarial parasitic nematode Brugia malayi. Although, quantitatively, the genome is even more degraded than those of closely related Rickettsia species, Wolbachia has retained more intact metabolic pathways. The ability to provide riboflavin, flavin adenine dinucleotide, heme, and nucleotides is likely to be Wolbachias principal contribution to the mutualistic relationship, whereas the host nematode likely supplies amino acids required for Wolbachia growth. Genome comparison of the Wolbachia endosymbiont of B. malayi (wBm) with the Wolbachia endosymbiont of Drosophila melanogaster (wMel) shows that they share similar metabolic trends, although their genomes show a high degree of genome shuffling. In contrast to wMel, wBm contains no prophage and has a reduced level of repeated DNA. Both Wolbachia have lost a considerable number of membrane biogenesis genes that apparently make them unable to synthesize lipid A, the usual component of proteobacterial membranes. However, differences in their peptidoglycan structures may reflect the mutualistic lifestyle of wBm in contrast to the parasitic lifestyle of wMel. The smaller genome size of wBm, relative to wMel, may reflect the loss of genes required for infecting host cells and avoiding host defense systems. Analysis of this first sequenced endosymbiont genome from a filarial nematode provides insight into endosymbiont evolution and additionally provides new potential targets for elimination of cutaneous and lymphatic human filarial disease.

Sequencing and Analyses of All Known Human Rhinovirus Genomes Reveal Structure and Evolution

Infection by human rhinovirus (HRV) is a major cause of upper and lower respiratory tract disease worldwide and displays considerable phenotypic variation. We examined diversity by completing the genome sequences for all known serotypes (n = 99). Superimposition of capsid crystal structure and optimal-energy RNA configurations established alignments and phylogeny. These revealed conserved motifs; clade-specific diversity, including a potential newly identified species (HRV-D); mutations in field isolates; and recombination. In analogy with poliovirus, a hypervariable 5′ untranslated region tract may affect virulence. A configuration consistent with nonscanning internal ribosome entry was found in all HRVs and may account for rapid translation. The data density from complete sequences of the reference HRVs provided high resolution for this degree of modeling and serves as a platform for full genome-based epidemiologic studies and antiviral or vaccine development.

Evolutionary dynamics of human rotaviruses: balancing reassortment with preferred genome constellations.

Group A human rotaviruses (RVs) are a major cause of severe gastroenteritis in infants and young children. Yet, aside from the genes encoding serotype antigens (VP7; G-type and VP4; P-type), little is known about the genetic make-up of emerging and endemic human RV strains. To gain insight into the diversity and evolution of RVs circulating at a single location over a period of time, we sequenced the eleven-segmented, double-stranded RNA genomes of fifty-one G3P[8] strains collected from 1974 to 1991 at Childrens Hospital National Medical Center, Washington, D. C. During this period, G1P[8] strains typically dominated, comprising on average 56% of RV infections each year in hospitalized children. A notable exception was in the 1976 and 1991 winter seasons when the incidence of G1P[8] infections decreased dramatically, a trend that correlated with a significant increase in G3P[8] infections. Our sequence analysis indicates that the 1976 season was characterized by the presence of several genetically distinct, co-circulating clades of G3P[8] viruses, which contained minor but significant differences in their encoded proteins. These 1976 lineages did not readily exchange gene segments with each other, but instead remained stable over the course of the season. In contrast, the 1991 season contained a single major clade, whose genome constellation was similar to one of the 1976 clades. The 1991 clade may have gained a fitness advantage after reassorting with as of yet unidentified RV strain(s). This study reveals for the first time that genetically distinct RV clades of the same G/P-type can co-circulate and cause disease. The findings from this study also suggest that, although gene segment exchange occurs, most reassortant strains are replaced over time by lineages with preferred genome constellations. Elucidation of the selective pressures that favor maintenance of RVs with certain sets of genes may be necessary to anticipate future vaccine needs.

VIGOR, an annotation program for small viral genomes

BackgroundThe decrease in cost for sequencing and improvement in technologies has made it easier and more common for the re-sequencing of large genomes as well as parallel sequencing of small genomes. It is possible to completely sequence a small genome within days and this increases the number of publicly available genomes. Among the types of genomes being rapidly sequenced are those of microbial and viral genomes responsible for infectious diseases. However, accurate gene prediction is a challenge that persists for decoding a newly sequenced genome. Therefore, accurate and efficient gene prediction programs are highly desired for rapid and cost effective surveillance of RNA viruses through full genome sequencing.ResultsWe have developed VIGOR (Viral Genome ORF Reader), a web application tool for gene prediction in influenza virus, rotavirus, rhinovirus and coronavirus subtypes. VIGOR detects protein coding regions based on sequence similarity searches and can accurately detect genome specific features such as frame shifts, overlapping genes, embedded genes, and can predict mature peptides within the context of a single polypeptide open reading frame. Genotyping capability for influenza and rotavirus is built into the program. We compared VIGOR to previously described gene prediction programs, ZCURVE_V, GeneMarkS and FLAN. The specificity and sensitivity of VIGOR are greater than 99% for the RNA viral genomes tested.ConclusionsVIGOR is a user friendly web-based genome annotation program for five different viral agents, influenza, rotavirus, rhinovirus, coronavirus and SARS coronavirus. This is the first gene prediction program for rotavirus and rhinovirus for public access. VIGOR is able to accurately predict protein coding genes for the above five viral types and has the capability to assign function to the predicted open reading frames and genotype influenza virus. The prediction software was designed for performing high throughput annotation and closure validation in a post-sequencing production pipeline.

Biologic, Antigenic, and Full-Length Genomic Characterization of a Bovine-Like Coronavirus Isolated from a Giraffe

ABSTRACT Coronaviruses (CoVs) possess large RNA genomes and exist as quasispecies, which increases the possibility of adaptive mutations and interspecies transmission. Recently, CoVs were recognized as important pathogens in captive wild ruminants. This is the first report of the isolation and detailed genetic, biologic, and antigenic characterization of a bovine-like CoV from a giraffe (Giraffa camelopardalis) in a wild-animal park in the United States. CoV particles were detected by immune electron microscopy in fecal samples from three giraffes with mild-to-severe diarrhea. From one of the three giraffe samples, a CoV (GiCoV-OH3) was isolated and successfully adapted to serial passage in human rectal tumor 18 cell cultures. Hemagglutination assays, receptor-destroying enzyme activity, hemagglutination inhibition, and fluorescence focus neutralization tests revealed close biological and antigenic relationships between the GiCoV-OH3 isolate and selected respiratory and enteric bovine CoV (BCoV) strains. When orally inoculated into a BCoV-seronegative gnotobiotic calf, GiCoV-OH3 caused severe diarrhea and virus shedding within 2 to 3 days. Sequence comparisons and phylogenetic analyses were performed to assess its genetic relatedness to other CoVs. Molecular characterization confirmed that the new isolate belongs to group 2a of the mammalian CoVs and revealed closer genetic relatedness between GiCoV-OH3 and the enteric BCoVs BCoV-ENT and BCoV-DB2, whereas BCoV-Mebus was more distantly related. Detailed sequence analysis of the GiCoV-OH3 spike gene demonstrated the presence of a deletion in the variable region of the S1 subunit (from amino acid 543 to amino acid 547), which is a region associated with pathogenicity and tissue tropism for other CoVs. The point mutations identified in the structural proteins (by comparing GiCoV-OH3, BCoV-ENT, BCoV-DB2, and BCoV-Mebus) were most conserved among GiCoV-OH3, BCoV-ENT, and BCoV-DB2, whereas most of the point mutations in the nonstructural proteins were unique to GiCoV-OH3. Our results confirm the existence of a bovine-like CoV transmissible to cattle from wild ruminants, namely, giraffes, but with certain genetic properties different from those of BCoVs.

Bovine-Like Coronaviruses Isolated from Four Species of Captive Wild Ruminants Are Homologous to Bovine Coronaviruses, Based on Complete Genomic Sequences

ABSTRACT We sequenced and analyzed the full-length genomes of four coronaviruses (CoVs), each from a distinct wild-ruminant species in Ohio: sambar deer (Cervus unicolor), a waterbuck (Kobus ellipsiprymnus), a sable antelope (Hippotragus niger), and a white-tailed deer (Odocoileus virginianus). The fecal samples from the sambar deer, the waterbuck, and the white-tailed deer were collected during winter dysentery outbreaks and sporadic diarrhea cases in 1993 and 1994 (H. Tsunemitsu, Z. R. el-Kanawati, D. R. Smith, H. H. Reed, and L. J. Saif, J. Clin. Microbiol. 33:3264-3269, 1995). A fecal sample from a sable antelope was collected in 2003 from an Ohio wild-animal habitat during the same outbreak when a bovine-like CoV from a giraffe (GiCoV) was isolated (M. Hasoksuz, K. Alekseev, A. Vlasova, X. Zhang, D. Spiro, R. Halpin, S. Wang, E. Ghedin, and L. J. Saif, J. Virol. 81:4981-4990, 2007). For two of the CoVs (sambar deer and waterbuck), complete genomes from both the cell culture-adapted and gnotobiotic-calf-passaged strains were also sequenced and analyzed. Phylogenetically, wild-ruminant CoVs belong to group 2a CoVs, with the closest relatedness to recent bovine CoV (BCoV) strains. High nucleotide identities (99.4 to 99.6%) among the wild-ruminant strains and recent BCoV strains (BCoV-LUN and BCoV-ENT, isolated in 1998) further confirm the close relatedness. Comparative genetic analysis of CoVs of captive wild ruminants with BCoV strains suggests that no specific genomic markers are present that allow discrimination between the bovine strains and bovine-like CoVs from captive wild ruminants; furthermore, no specific genetic markers were identified that defined cell cultured or calf-passaged strains or the host origin of strains. The results of this study confirm prior reports of biologic and antigenic similarities between bovine and wild-ruminant CoVs and suggest that cattle may be reservoirs for CoVs that infect captive wild ruminants or vice versa and that these CoVs may represent host range variants of an ancestral CoV.

Complete genomic sequences, a key residue in the spike protein and deletions in nonstructural protein 3b of US strains of the virulent and attenuated coronaviruses, transmissible gastroenteritis virus and porcine respiratory coronavirus.

Abstract Transmissible gastroenteritis virus (TGEV) isolates that have been adapted to passage in cell culture maintain their infectivity in vitro but may lose their pathogenicity in vivo. To better understand the genomic mechanisms for viral attenuation, we sequenced the complete genomes of two virulent TGEV strains and their attenuated counterparts: virulent TGEV Miller M6 and attenuated TGEV Miller M60 and virulent TGEV Purdue and attenuated TGEV Purdue P115, together with the ISU-1 strain of porcine respiratory coronavirus (PRCV-ISU-1), a naturally occurring TGEV deletion mutant with an altered respiratory tropism and reduced virulence. Pairwise comparison at both the nucleotide (nt) and amino acid (aa) levels between virulent and attenuated TGEV strains identified a common change in nt 1753 of the spike gene, resulting in a serine to alanine mutation at aa position 585 of the spike proteins of the attenuated TGEV strains. Alanine was also present in this protein in PRCV-ISU-1. Particularly noteworthy, the serine to alanine mutation resides in the region of the major antigenic site A/B (aa 506–706) that elicits neutralizing antibodies and within the domain mediating the cell surface receptor aminopeptidase N binding (aa 522–744). Comparison of the predicted polypeptide products of ORF3b showed significant deletions in the naturally attenuated PRCV-ISU-1 and TGEV Miller M60; these deletions occurred at a common break point, suggesting a related mechanism of recombination that may affect viral virulence or tropism. Sequence comparisons at both genomic and protein levels indicated that PRCV-ISU-1 had a closer relationship with TGEV Miller strains than Purdue strains. Phylogenetic analyses showed that virulence is an evolutionarily labile trait in TGEV and that TGEV strains as a group share a common ancestor with PRCV.

VIGOR extended to annotate genomes for additional 12 different viruses

A gene prediction program, VIGOR (Viral Genome ORF Reader), was developed at J. Craig Venter Institute in 2010 and has been successfully performing gene calling in coronavirus, influenza, rhinovirus and rotavirus for projects at the Genome Sequencing Center for Infectious Diseases. VIGOR uses sequence similarity search against custom protein databases to identify protein coding regions, start and stop codons and other gene features. Ribonucleicacid editing and other features are accurately identified based on sequence similarity and signature residues. VIGOR produces four output files: a gene prediction file, a complementary DNA file, an alignment file, and a gene feature table file. The gene feature table can be used to create GenBank submission. VIGOR takes a single input: viral genomic sequences in FASTA format. VIGOR has been extended to predict genes for 12 viruses: measles virus, mumps virus, rubella virus, respiratory syncytial virus, alphavirus and Venezuelan equine encephalitis virus, norovirus, metapneumovirus, yellow fever virus, Japanese encephalitis virus, parainfluenza virus and Sendai virus. VIGOR accurately detects the complex gene features like ribonucleicacid editing, stop codon leakage and ribosomal shunting. Precisely identifying the mat_peptide cleavage for some viruses is a built-in feature of VIGOR. The gene predictions for these viruses have been evaluated by testing from 27 to 240 genomes from GenBank.