Calculation of Evolutionary Correlation between Individual Genes and Full-Length Genome: A Method Useful for Choosing Phylogenetic Markers for Molecular Epidemiology

Individual genes or regions are still commonly used to estimate the phylogenetic relationships among viral isolates. The genomic regions that can faithfully provide assessments consistent with those predicted with full-length genome sequences would be preferable to serve as good candidates of the phylogenetic markers for molecular epidemiological studies of many viruses. Here we employed a statistical method to evaluate the evolutionary relationships between individual viral genes and full-length genomes without tree construction as a way to determine which gene can match the genome well in phylogenetic analyses. This method was performed by calculation of linear correlations between the genetic distance matrices of aligned individual gene sequences and aligned genome sequences. We applied this method to the phylogenetic analyses of porcine circovirus 2 (PCV2), measles virus (MV), hepatitis E virus (HEV) and Japanese encephalitis virus (JEV). Phylogenetic trees were constructed for comparisons and the possible factors affecting the method accuracy were also discussed in the calculations. The results revealed that this method could produce results consistent with those of previous studies about the proper consensus sequences that could be successfully used as phylogenetic markers. And our results also suggested that these evolutionary correlations could provide useful information for identifying genes that could be used effectively to infer the genetic relationships.     

戊型肝炎病毒通用性PCR引物的设计及其基因分型的研究

针对戊型肝炎病毒(HEV)基因分型尚无统一标准的现状,本研究通过分析GenBank中现有的82株HEV基因组全序列,设计一组HEV通用性PCR引物(HEVuPrimer)用于扩增不同基因型HEV可测序长片段,分别基于全基因序列、HEVuPrimer扩增序列和较常用的MXJ引物扩增序列对82株HEV进行基因分型,再以HEVuPrimer扩增1～4型HEV参考株和HEVIg M抗体阳性的临床标本,进行HEV基因分型的研究。研究结果表明HE-VuPrimer扩增区段与全基因序列对82株HEV的基因分型结果完全一致;HEVuPrimer与HEV序列的匹配程度明显高于MXJ引物,且HEVuPrimer区段的基因分型结果较MXJ区段的更为准确。HEVuPrimer可同时扩增出不同基因型HEV基因片段。124份临床标本中,60份测出特异性HEV RNA,阳性率为48.4%,基因分型结果均为4型HEV,但分属于4个不同的基因亚型,核苷酸同源性为80.0%～99.9%,其中有6例近期分离的HEV毒株形成一新的基因亚型。因此,基于HEVuPrimer扩增区段的HEV基因分型方法具有较高的可靠性和可信度,为建立统一、可行的HEV基...     

Parallel phylogenetic analyses using the N, G or Nv gene from a fixed group of VHSV isolates reveal the same overall genetic typing

Different genetic regions representing the viral phospho- (P), nucleocapsid- (N) or glyco-protein (G) gene have been used for phylogenetic studies of viral haemorrhagic septicaemia virus (VHSV). Since these analyses were performed on different virus isolates using various genomic regions, it has been difficult to evaluate how the choice of target region affects the output of the analyses. To address this, we sequenced and performed parallel phylogenetic analysis of an N gene fragment, the entire Nv (non-structural protein) and G genes, and 4 different fragments of the G gene from a fixed virus panel. The overall genotyping of the selected isolates was identical for the 7 target regions, but separation of Genotype I sub-lineages was best when the analysis was performed on the full length G gene (1524 nucleotides, nt). Good resolution was furthermore obtained using smaller sequencing windows represented by a G gene fragment (nt 360 to 720) or the Nv gene (366 nt), although these regions had different characteristics with respect to resolution of Genotype I sublineages and resolution within Sub-lineage Ia. Phylogenetic analysis based on the deduced amino acid sequences was also performed. The phylogenetic relationship between the nucleotide and amino acid sequences of the isolates corresponded best in the case of the N gene/protein. For the 6 other genomic regions, genetically distant isolates occasionally grouped together when compared at protein levels. No clear relationship between the G gene genotyping and serotyping with neutralising (G protein specific) antibodies was observed, stressing that epidemiological analysis based on phenotypic characteristics such as serotype could be misleading.     

Detection of genomic islands via segmental genome heterogeneity

While the recognition of genomic islands can be a powerful mechanism for identifying genes that distinguish related bacteria, few methods have been developed to identify them specifically. Rather, identification of islands often begins with cataloging individual genes likely to have been recently introduced into the genome; regions with many putative alien genes are then examined for other features suggestive of recent acquisition of a large genomic region. When few phylogenetic relatives are available, the identification of alien genes relies on their atypical features relative to the bulk of the genes in the genome. The weakness of these ‘bottom–up’ approaches lies in the difficulty in identifying robustly those genes which are atypical, or phylogenetically restricted, due to recent foreign ancestry. Herein, we apply an alternative ‘top–down’ approach where bacterial genomes are recursively divided into progressively smaller regions, each with uniform composition. In this way, large chromosomal regions with atypical features are identified with high confidence due to the simultaneous analysis of multiple genes. This approach is based on a generalized divergence measure to quantify the compositional difference between segments in a hypothesis-testing framework. We tested the proposed genome island prediction algorithm on both artificial chimeric genomes and genuine bacterial genomes.     

Genetic relationship between cocirculating Human enteroviruses species C

Recombination events between human enteroviruses (HEV) are known to occur frequently and to participate in the evolution of these viruses. In a previous study, we reported the isolation of a panel of viruses belonging to the Human enterovirus species C (HEV-C) that had been cocirculating in a small geographic area of Madagascar in 2002. This panel included type 2 vaccine-derived polioviruses (PV) that had caused several cases of acute flaccid paralysis in humans. Previous partial sequencing of the genome of these HEV-C isolates revealed considerable genetic diversity, mostly due to recombination. In the work presented herein, we carried out a more detailed characterization of the genomes of viruses from this collection. First, we determined the full VP1 sequence of 41 of these isolates of different types. These sequences were compared with those of HEV-C isolates obtained from other countries or in other contexts. The sequences of the Madagascan isolates of a given type formed specific clusters clearly differentiated from those formed by other strains of the same type isolated elsewhere. Second, we sequenced the entire genome of 10 viruses representing most of the lineages present in this panel. All but one of the genomes appeared to be mosaic assemblies of different genomic fragments generated by intra- and intertypic recombination. The location of the breakpoints suggested potential preferred genomic regions for recombination. Our results also suggest that recombination between type HEV-99 and other HEV-C may be quite rare. This first exhaustive genomic analysis of a panel of non-PV HEV-C cocirculating in a small human population highlights the high frequency of inter and intra-typic genetic recombination, constituting a widespread mechanism of genetic plasticity and continually shifting the HEV-C biodiversity.     

Characterization of full-length enterovirus 71 strains from severe and mild disease patients in northeastern China

Human enterovirus 71 (EV71)-associated hand, foot, and mouth disease (HFMD) has been a leading cause of childhood infection in China since 2008. Epidemic and molecular characteristics of HFMD have been examined in many areas of China, including the central and southern regions. However, clinical and genetic characterization of EV71 in the northeastern region of China is scarce. In this study, a series of analyses were performed on seven full-length EV71 sequences from HFMD patients who had either severe or mild disease. We have determined that these seven circulating EV71 viruses from Changchun, China are actually complex recombinant viruses involving multiple type A human enterovirus (HEV). Classified as EV71 subtype C4 (EV71 C4), these Changchun EV71 viruses contain genetic recombination events between the CA4, CA5, EV71B4 and EV71C1 strains. Most of the structural protein region (P1) of these viruses resembled that of the prototype EV71 C1 strains. The non-structural protein domains (P2 and P3) showed a high degree of similarity with CA4, CA5 and EV71 B4 in different regions. The 5'UTR had unclassified recombination,while partial 3D region of these viruses showed a high degree of similarity to CA16. Phylogenetic analysis of full-length or partial sequences of isolates from severe or mild disease patients in Changchun always formed a single cluster in various phylogenetic analyses of different genomic regions, suggesting that all seven strains originated from one single common ancestor. There was no correlation between viral genomic sequence and virulence. Thus, we found that circulating recombinant forms of EV71 are prevalent among HFMD patients in Northeastern China. The existence of a unique cluster of EV71 related viruses in Northeast China has important implications for vaccine development that would address the increasing prevalence of HFMD.     

Genomic Evolution of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Isolates Revealed by Deep Sequencing

Most studies on PRRSV evolution have been limited to a particular region of the viral genome. A thorough genome-wide understanding of the impact of different mechanisms on shaping PRRSV genetic diversity is still lacking. To this end, deep sequencing was used to obtain genomic sequences of a diverse set of 16 isolates from a region of Hong Kong with a complex PRRSV epidemiological record. Genome assemblies and phylogenetic typing indicated the co-circulation of strains of both genotypes (type 1and type 2) with varying Nsp2 deletion patterns and distinct evolutionary lineages (“High Fever”-like and local endemic type). Recombination analyses revealed genomic breakpoints in structural and non-structural regions of genomes of both genotypes with evidence of many recombination events originating from common ancestors. Additionally, the high fold of coverage per nucleotide allowed the characterization of minor variants arising from the quasispecies of each strain. Overall, 0.56–2.83% of sites were found to be polymorphic with respect to cognate consensus genomes. The distribution of minor variants across each genome was not uniform indicating the influence of selective forces. Proportion of variants capable of causing an amino acid change in their respective codons ranged between 25–67% with many predicted to be non-deleterious. Low frequency deletion variants were also detected providing one possible mechanism for their sudden emergence as cited in previous reports.     

Phylogeny and Differentiation of Reptilian and Amphibian Ranaviruses Detected in Europe

Ranaviruses in amphibians and fish are considered emerging pathogens and several isolates have been extensively characterized in different studies. Ranaviruses have also been detected in reptiles with increasing frequency, but the role of reptilian hosts is still unclear and only limited sequence data has been provided. In this study, we characterized a number of ranaviruses detected in wild and captive animals in Europe based on sequence data from six genomic regions (major capsid protein (MCP), DNA polymerase (DNApol), ribonucleoside diphosphate reductase alpha and beta subunit-like proteins (RNR-α and -β), viral homolog of the alpha subunit of eukaryotic initiation factor 2, eIF-2α (vIF-2α) genes and microsatellite region). A total of ten different isolates from reptiles (tortoises, lizards, and a snake) and four ranaviruses from amphibians (anurans, urodeles) were included in the study. Furthermore, the complete genome sequences of three reptilian isolates were determined and a new PCR for rapid classification of the different variants of the genomic arrangement was developed. All ranaviruses showed slight variations on the partial nucleotide sequences from the different genomic regions (92.6–100%). Some very similar isolates could be distinguished by the size of the band from the microsatellite region. Three of the lizard isolates had a truncated vIF-2α gene; the other ranaviruses had full-length genes. In the phylogenetic analyses of concatenated sequences from different genes (3223 nt/10287 aa), the reptilian ranaviruses were often more closely related to amphibian ranaviruses than to each other, and most clustered together with previously detected ranaviruses from the same geographic region of origin. Comparative analyses show that among the closely related amphibian-like ranaviruses (ALRVs) described to date, three recently split and independently evolving distinct genetic groups can be distinguished. These findings underline the wide host range of ranaviruses and the emergence of pathogen pollution via animal trade of ectothermic vertebrates.     

Plant Virology and Next Generation Sequencing: Experiences with a Potyvirus

Next generation sequencing is quickly emerging as the go-to tool for plant virologists when sequencing whole virus genomes, and undertaking plant metagenomic studies for new virus discoveries. This study aims to compare the genomic and biological properties of Bean yellow mosaic virus (BYMV) (genus Potyvirus), isolates from Lupinus angustifolius plants with black pod syndrome (BPS), systemic necrosis or non-necrotic symptoms, and from two other plant species. When one Clover yellow vein virus (ClYVV) (genus Potyvirus) and 22 BYMV isolates were sequenced on the Illumina HiSeq2000, one new ClYVV and 23 new BYMV sequences were obtained. When the 23 new BYMV genomes were compared with 17 other BYMV genomes available on Genbank, phylogenetic analysis provided strong support for existence of nine phylogenetic groupings. Biological studies involving seven isolates of BYMV and one of ClYVV gave no symptoms or reactions that could be used to distinguish BYMV isolates from L. angustifolius plants with black pod syndrome from other isolates. Here, we propose that the current system of nomenclature based on biological properties be replaced by numbered groups (I–IX). This is because use of whole genomes revealed that the previous phylogenetic grouping system based on partial sequences of virus genomes and original isolation hosts was unsustainable. This study also demonstrated that, where next generation sequencing is used to obtain complete plant virus genomes, consideration needs to be given to issues regarding sample preparation, adequate levels of coverage across a genome and methods of assembly. It also provided important lessons that will be helpful to other plant virologists using next generation sequencing in the future.     

Genomic diversity and evolution of the lyssaviruses

Lyssaviruses are RNA viruses with single-strand, negative-sense genomes responsible for rabies-like diseases in mammals. To date, genomic and evolutionary studies have most often utilized partial genome sequences, particularly of the nucleoprotein and glycoprotein genes, with little consideration of genome-scale evolution. Herein, we report the first genomic and evolutionary analysis using complete genome sequences of all recognised lyssavirus genotypes, including 14 new complete genomes of field isolates from 6 genotypes and one genotype that is completely sequenced for the first time. In doing so we significantly increase the extent of genome sequence data available for these important viruses. Our analysis of these genome sequence data reveals that all lyssaviruses have the same genomic organization. A phylogenetic analysis reveals strong geographical structuring, with the greatest genetic diversity in Africa, and an independent origin for the two known genotypes that infect European bats. We also suggest that multiple genotypes may exist within the diversity of viruses currently classified as 'Lagos Bat'. In sum, we show that rigorous phylogenetic techniques based on full length genome sequence provide the best discriminatory power for genotype classification within the lyssaviruses.     

Characterization of a hypervariable region in the genome of Chlamydophila pneumoniae

Chlamydophila pneumoniae displays surprisingly little genomic variation, as seen by comparisons of the published genomes from three different isolates and sequencing of four different genes from different isolates. We have in the present study, however, demonstrated genomic variation between 10 C. pneumoniae isolates in the 11690-bp region between the two outer membrane protein genes pmp1 and pmp2. This region of the C. pneumoniae CWL-029 isolate contains seven C. pneumoniae-specific open reading frames (hb1-7, encoding hydrophobic beta-sheet-containing proteins). We identified additionally 12 open reading frames in the C. pneumoniae CWL-029 genome encoding hypothetical proteins with similarity to the seven hypothetical Hb-proteins. Compared to other isolates, genomic variation is seen to cause frame-shifting of three of the 19 hb-open reading frames, which are proposed to be three full-length genes and eight frame-shifted pseudogenes. The hypothetical proteins encoded by these proposed genes contain an N-terminally located highly hydrophobic stretch of 50-60 residues. A similar motif is found in all identified Chlamydia inclusion membrane proteins and therefore the Hb-proteins are candidate inclusion proteins.     

Tobamovirus evolution: gene overlaps, recombination, and taxonomic implications

Tobamoviruses, mostly isolated from solanaceous plants, may represent ancient virus lineages that have codiverged with their hosts. Recently completed nucleotide sequences of six nonsolanaceous tobamoviruses allowed assessment of the codivergence hypothesis and support a third subgroup within tobamoviruses. The genomic sequences of 12 tobamoviruses and the partial sequences of 11 others have been analyzed. Comparisons of the predicted protein sequences revealed three clusters of tobamoviruses, corresponding to those infecting solanaceous species (subgroup 1), those infecting cucurbits and legumes (subgroup 2), and those infecting crucifers. The orchid-infecting odontoglossum ringspot tobamovirus was associated with subgroup 1 genomes by its coat and movement protein sequences, but with the crucifer-pathogenic tobamoviruses by the remainder of its genome, suggesting that it is the progeny of a recombinant. For four of five genomic regions, subgroup 1 and 3 genomes were equidistant from a subgroup 2 genome chosen for comparison, suggesting uniform rates of evolution. A phylogenetic tree of plant families based on the tobamoviruses they harbor was congruent with that based on rubisco sequences but had a different root, suggesting that codivergence was tempered by rare events of viruses of one family colonizing another family. The proposed subgroup 3 viruses probably have an origin of virion assembly in the movement protein gene, a large (25-codon) overlap of movement and coat protein open reading frames, and a comparably shorter genome. Codon-position- dependent base compositions and codon prevalences suggested that the coat protein frame of the overlap region was ancestral. Bootstrapped parsimony analysis of the nucleotides in the overlap region and of the sequences translated from the -1 frame (the subgroup 3 movement protein frame) of this region produced trees inconsistent with those deduced from other regions. The results are consistent with a model in which a no or short overlap organization was ancestral. Despite encoding of subgroup 2 and 3 movement protein C-termini by nonhomologous nucleotides, weak similarities between their amino acid sequences suggested convergent sequence evolution.      

Genome analysis and CRISPR typing of Salmonella enterica serovar Virchow

Background

Salmonella enterica subsp. enterica serovar Virchow has been recognized as a significant health burden in Asia, Australia and Europe. In addition to its global distribution, S. Virchow is clinically significant due to the frequency at which it causes invasive infections and its association with outbreaks arising from food-borne transmission. Here, we examine the genome of an invasive isolate of S. Virchow SVQ1 (phage type 8) from an outbreak in southeast Queensland, Australia. In addition to identifying new potential genotyping targets that could be used for discriminating between S. Virchow strains in outbreak scenarios, we also aimed to carry out a comprehensive comparative analysis of the S. Virchow genomes.

Results

Genome comparisons between S. Virchow SVQ1 and S. Virchow SL491, a previously published strain, identified a high degree of genomic similarity between the two strains with fewer than 200 single nucleotide differences. Clustered Regularly Interspaced Palindromic Repeats (CRISPR) regions were identified as a highly variable region that could be used to discriminate between S. Virchow isolates. We amplified and sequenced the CRISPR regions of fifteen S. Virchow isolates collected from seven different outbreaks across Australia. We observed three allelic types of the CRISPR region from these isolates based on the presence/absence of the spacers and were able to discriminate S. Virchow phage type 8 isolates originating from different outbreaks. A comparison with 27 published Salmonella genomes found that the S. Virchow SVQ1 genome encodes 11 previously described Salmonella Pathogenicity Islands (SPI), as well as additional genomic islands including a remnant integrative conjugative element that is distinct from SPI-7. In addition, the S. Virchow genome possesses a novel prophage that encodes the Type III secretion system effector protein SopE, a key Salmonella virulence factor. The prophage shares very little similarity to the SopE prophages found in other Salmonella serovars suggesting an independent acquisition of sopE.

Conclusions

The availability of this genome will serve as a genome template and facilitate further studies on understanding the virulence and global distribution of the S. Virchow serovar, as well as the development of genotyping methods for outbreak investigations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-389) contains supplementary material, which is available to authorized users.     

Molecular epidemiology of clinically high-risk Pseudomonas aeruginosa strains: Practical overview

In recent years, numerous outbreaks of multidrug-resistant Pseudomonas aeruginosa have been reported across the world. Once an outbreak occurs, besides routinely testing isolates for susceptibility to antimicrobials, it is required to check their virulence genotypes and clonality profiles. Replacing pulsed-field gel electrophoresis DNA fingerprinting are faster, easier-to-use, and less expensive polymerase chain reaction (PCR)-based methods for characterizing hospital isolates. P. aeruginosa possesses a mosaic genome structure and a highly conserved core genome displaying low sequence diversity and a highly variable accessory genome that communicates with other Pseudomonas species via horizontal gene transfer. Multiple-locus variable-number tandem-repeat analysis and multilocus sequence typing methods allow for phylogenetic analysis of isolates by PCR amplification of target genes with the support of Internet-based services. The target genes located in the core genome regions usually contain low-frequency mutations, allowing the resulting phylogenetic trees to infer evolutionary processes. The multiplex PCR-based open reading frame typing (POT) method, integron PCR, and exoenzyme genotyping can determine a genotype by PCR amplifying a specific insertion gene in the accessory genome region using a single or a multiple primer set. Thus, analyzing P. aeruginosa isolates for their clonality, virulence factors, and resistance characteristics is achievable by combining the clonality evaluation of the core genome based on multiple-locus targeting methods with other methods that can identify specific virulence and antimicrobial genes. Software packages such as eBURST, R, and Dendroscope, which are powerful tools for phylogenetic analyses, enable researchers and clinicians to visualize clonality associations in clinical isolates.     

Eight New Genomes and Synthetic Controls Increase the Accessibility of Rapid Melt-MAMA SNP Typing of Coxiella burnetii

The case rate of Q fever in Europe has increased dramatically in recent years, mainly because of an epidemic in the Netherlands in 2009. Consequently, there is a need for more extensive genetic characterization of the disease agent Coxiella burnetii in order to better understand the epidemiology and spread of this disease. Genome reference data are essential for this purpose, but only thirteen genome sequences are currently available. Current methods for typing C. burnetii are criticized for having problems in comparing results across laboratories, require the use of genomic control DNA, and/or rely on markers in highly variable regions. We developed in this work a method for single nucleotide polymorphism (SNP) typing of C. burnetii isolates and tissue samples based on new assays targeting ten phylogenetically stable synonymous canonical SNPs (canSNPs). These canSNPs represent previously known phylogenetic branches and were here identified from sequence comparisons of twenty-one C. burnetii genomes, eight of which were sequenced in this work. Importantly, synthetic control templates were developed, to make the method useful to laboratories lacking genomic control DNA. An analysis of twenty-one C. burnetii genomes confirmed that the species exhibits high sequence identity. Most of its SNPs (7,493/7,559 shared by >1 genome) follow a clonal inheritance pattern and are therefore stable phylogenetic typing markers. The assays were validated using twenty-six genetically diverse C. burnetii isolates and three tissue samples from small ruminants infected during the epidemic in the Netherlands. Each sample was assigned to a clade. Synthetic controls (vector and PCR amplified) gave identical results compared to the corresponding genomic controls and are viable alternatives to genomic DNA. The results from the described method indicate that it could be useful for cheap and rapid disease source tracking at non-specialized laboratories, which requires accurate genotyping, assay accessibility and inter-laboratory comparisons.     

Patterns of genomic sequence diversity among their simian immunodeficiency viruses suggest that L'Hoest monkeys (Cercopithecus lhoesti) are a natural lentivirus reservoir

Recently, we described a novel simian immunodeficiency virus (SIVlhoest) from a wild-caught L'Hoest monkey (Cercopithecus lhoesti) from a North American zoo. To investigate whether L'Hoest monkeys are the natural host for these viruses, we have screened blood samples from 14 wild animals from the Democratic Republic of Congo. Eight (57%) were found to be seropositive for SIV. Nearly full-length genome sequences were obtained for SIV isolates from three of these monkeys and compared to the original isolate and to other SIVs. The four samples of SIVlhoest formed a distinct cluster in phylogenetic trees. Two of these isolates differed on average at only about 5% of nucleotides, suggesting that they were epidemiologically linked; otherwise, the SIVlhoest isolates differed on average by 18%. Both the level of diversity and the pattern of its variation along the genome were very similar to those seen among isolates of SIVagm from vervet monkeys, pointing to similarities in the nature of, and constraints on, SIV evolution in these two species. Discordant phylogenetic relationships among the SIVlhoest isolates for different genomic regions indicated that mosaic viruses have been generated by recombination, implying that individual monkeys have been coinfected by more than one strain of SIV. Taken together, these observations provide strong evidence that L'Hoest monkeys constitute a natural reservoir for SIV.     

Reference-free SNP calling: improved accuracy by preventing incorrect calls from repetitive genomic regions

ABSTRACT: BACKGROUND: Single nucleotide polymorphisms (SNPs) are the most abundant type of genetic variation in eukaryotic genomes and have recently become the marker of choice in a wide variety of ecological and evolutionary studies. The advent of next-generation sequencing (NGS) technologies has made it possible to efficiently genotype a large number of SNPs in the non-model organisms with no or limited genomic resources. Most NGS-based genotyping methods require a reference genome to perform accurate SNP calling. Little effort, however, has yet been devoted to developing or improving algorithms for accurate SNP calling in the absence of a reference genome. RESULTS: Here we describe an improved maximum likelihood (ML) algorithm called iML, which can achieve high genotyping accuracy for SNP calling in the non-model organisms without a reference genome. The iML algorithm incorporates the mixed Poisson/normal model to detect composite read clusters and can efficiently prevent incorrect SNP calls resulting from repetitive genomic regions. Through analysis of simulation and real sequencing datasets, we demonstrate that in comparison with ML or a threshold approach, iML can remarkably improve the accuracy of de novo SNP genotyping and is especially powerful for the reference-free genotyping in diploid genomes with high repeat contents. CONCLUSIONS: The iML algorithm can efficiently prevent incorrect SNP calls resulting from repetitive genomic regions, and thus outperforms the original ML algorithm by achieving much higher genotyping accuracy. Our algorithm is therefore very useful for accurate de novo SNP genotyping in the non-model organisms without a reference genome.     

2009年云南省2株肠道病毒71型分离株全基因组序列遗传特性分析

对2009年云南省肠道病毒71型分离株KMM09和KM186-09进行全基因组序列测序,并与我国及其它国家流行的EV71基因型进行比较和进化分析。KMM09和KM186-09基因组长为7 409bp,编码2 193个氨基酸,VP1系统进化分析显示2009年云南分离株属于C4基因型的C4a亚型。在结构区,与其它基因型相比较,C基因型之间的核苷酸和氨基酸的同源性高于其它基因型;而在非结构区,C4与B基因型和CA16原型株G10同源性高于其它C基因亚型。通过RDP3重组软件和blast比对分析,发现EV71C4基因型与B3基因型,与CA16原型株G10的基因组在非结构区存在重组。EV71全基因组序列的比较和分析,对了解引起我国手足口病暴发或流行C4基因亚型EV71毒株的遗传特性具有重要意义。     

The ORF2 protein of hepatitis E virus binds the 5' region of viral RNA

Hepatitis E virus (HEV) is a major human pathogen in much of the developing world. It is a plus-strand RNA virus with a 7.2-kb polyadenylated genome consisting of three open reading frames, ORF1, ORF2, and ORF3. Of these, ORF2 encodes the major capsid protein of the virus and ORF3 encodes a small protein of unknown function. Using the yeast three-hybrid system and traditional biochemical techniques, we have studied the RNA binding activities of ORF2 and ORF3, two proteins encoded in the 3' structural part of the genome. Since the genomic RNA from HEV has been postulated to contain secondary structures at the 5' and 3' ends, we used these two terminal regions, besides other regions within the genome, in this study. Experiments were designed to test for interactions between the genomic RNA fusion constructs with ORF2 and ORF3 hybrid proteins in a yeast cellular environment. We show here that the ORF2 protein contains RNA binding activity. The ORF2 protein specifically bound the 5' end of the HEV genome. Deletion analysis of this protein showed that its RNA binding activity was lost when deletions were made beyond the N-terminal 111 amino acids. Finer mapping of the interacting RNA revealed that a 76-nucleotide (nt) region at the 5' end of the HEV genome was responsible for binding the ORF2 protein. This 76-nt region included the 51-nt HEV sequence, conserved across alphaviruses. Our results support the requirement of this conserved sequence for interaction with ORF2 and also indicate an increase in the strength of the RNA-protein interaction when an additional 44 bases downstream of this 76-nt region were included. Secondary-structure predictions and the location of the ORF2 binding region within the HEV genome indicate that this interaction may play a role in viral encapsidation.