J Gen Virol Try IJSEM Online
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 87 (2006), 3297-3305; DOI 10.1099/vir.0.82235-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Figures
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Gritsun, T. S.
Right arrow Articles by Gould, E. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Gritsun, T. S.
Right arrow Articles by Gould, E. A.
Agricola
Right arrow Articles by Gritsun, T. S.
Right arrow Articles by Gould, E. A.
© 2006 Society for General Microbiology

Direct repeats in the 3' untranslated regions of mosquito-borne flaviviruses: possible implications for virus transmission

T. S. Gritsun and E. A. Gould

Centre for Ecology and Hydrology, Mansfield Road, Oxford OX1 3SR, UK

Correspondence
T. S. Gritsun
tsg{at}ceh.ac.uk


   ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Direct repeats (DRs) of 20–45 nucleotide conserved sequences (CS) and repeated CS (RCS), separated by non-conserved sequences up to 100 nucleotides long, were previously described in the 3' untranslated region (3'UTR) of the three major mosquito-borne flavivirus (MBFV) subgroups, represented by Japanese encephalitis virus, Yellow fever virus and Dengue virus. Each subgroup exhibits a specific pattern of DRs, the biological significance of which has not yet been adequately addressed. The DRs were originally identified using conventional alignment programs based on the assumption that genetic variation is driven primarily by nucleotide substitutions. Since there are no recognized alignment programs that can adequately accommodate very divergent sequences, a method has been devised to construct and analyse a substantially improved 3'UTR alignment between these highly divergent viruses, based on the concept that deletions and/or insertions, in addition to substitutions, are important drivers of 3'UTR evolution. This ‘robust alignment’ approach demonstrated more extensive homologies in the 3'UTR than had been recognized previously and revealed the presence of similar DRs, either intact or as sequence ‘remnants’, in all the MBFV subgroups. The relevance of these observations is discussed in relation to (i) the function of DRs as elements of replication enhancement, (ii) the evolution of RNA secondary structures and (iii) the significance of DRs and secondary structures in MBFV transmissibility between vertebrate and invertebrate hosts.

Supplementary figures are available in JGV Online.


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
The genus Flavivirus (family Flaviviridae) contains four groups, the tick-borne flaviviruses (TBFV), the mosquito-borne flaviviruses (MBFV), the no-known-vector flaviviruses (NKVFV) and the non-classified flaviviruses (NCFV). The MBFV are further subdivided into three subgroups, the Japanese encephalitis virus (JEV), Yellow fever virus (YFV) and Dengue virus (DENV) subgroups, each of which includes human pathogens (Heinz et al., 2000Down; Gould et al., 2001Down). They are enveloped viruses (50 nm in diameter) with single-stranded positive-sense RNA genomes of approximately 11 000 bases. The RNA contains a single-stranded open reading frame (ORF) that encodes a polyprotein (about 3400 amino acids) which is processed co-translationally to produce the structural (envelope, membrane and capsid) and non-structural proteins, the latter providing the proteolytic and replicative functions (Lindenbach & Rice, 2001Down).

The ORF is flanked by the 5' and 3' untranslated regions (UTRs). In common with other positive-stranded RNA viruses, the 5'UTR initiates virus translation and the 3'UTR might also contribute to this process (Li & Brinton, 2001Down; Edgil et al., 2003Down; Holden & Harris, 2004Down; Chiu et al., 2005Down; Tilgner et al., 2005Down). The 3'UTR initiates virus replication, also mediated by interaction with the 5'UTR (You & Padmanabhan, 1999Down; Ackermann & Padmanabhan, 2001Down; Charrel et al., 2001Down; Khromykh et al., 2001Down; You et al., 2001Down; Corver et al., 2003Down; Lo et al., 2003Down). The mode of replication of flavivirus RNA is asymmetrical, with an excess of positive-strand RNA and no free minus strands in infected cells (Westaway et al., 2002Down, 2003Down).

Previous studies of the MBFV 3'UTR identified short direct repeats (DRs) as 20–45 nucleotide regions conserved between the JEV, YFV and DENV subgroups and separated by more divergent sequences. The MBFV DRs were initially specified as CSs and RCSs (conserved/repeated conserved sequences) and also YFV-R (YFV repeats) (Hahn et al., 1987Down). CS1 and CS2 sequences were present in all three MBFV subgroups, whereas the RCS2 was identified only in the JEV and DENV subgroups (Hahn et al., 1987Down). Surprisingly, CS2 and RCS2 were also present in the 3'UTR of the NKVFV group (Charlier et al., 2002Down), which is phylogenetically quite distant from the MBFV. In addition, CS3 and RCS3 were identified only in the JEV subgroups, whereas the three tandem repeats YF-R1, YF-R2 and YF-R3 were specific for the YFV and YFV-related subgroup (Hahn et al., 1987Down; Mutebi et al., 2004Down).

Experiments using infectious clones and replicons revealed that, with the exception of CS1, all DRs could be deleted without the engineered viruses losing their replicative capacity or viability (Men et al., 1996Down; Khromykh & Westaway, 1997Down; Mandl et al., 1998Down; Pletnev, 2001Down; Bredenbeek et al., 2003Down; Lo et al., 2003Down; Tilgner et al., 2005Down). These observations raised a number of questions: why are DRs so well preserved, why are they preserved as repeats and how did they originate?

We have constructed an alignment between the 3'UTRs of MBFV based on the principle that the evolutionary constraints in the UTR could be different from those in the ORF, allowing greater nucleotide perturbations such as large deletions and/or insertions. This approach enabled us to reveal sequence ‘remnants' of CSs/RCSs in all MBFV subgroups, leading us to propose that all CSs/RCSs were present in early MBFV lineages and that their preservation within the individual subgroups was driven by their specific function as enhancers of virus replication. These findings advance our understanding of the molecular organization of flavivirus UTRs and in particular their evolution and relevance to virus transmission and dissemination between vertebrate and invertebrate hosts.


   METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
The MBFV subgroups are highly divergent, and the mean pairwise nucleotide sequence identity of the 3'UTRs for the JEV, DENV and YFV subgroup viruses is too low (<50 %) to enable automated construction of a robust alignment. To overcome this problem, we initially constructed three comparative alignments for available 3'UTR sequences of 55 viruses of the JEV subgroup, 137 viruses of the DENV subgroup and 39 viruses of the YFV subgroup using CLUSTAL X followed by manual editing. The 3'UTRs of the YFV-related viruses Banzi virus (BANZV), Uganda S virus (UGSV), Sepik virus (SEPV) (Mutebi et al., 2004Down) and Yokose virus (YOKV) (Tajima et al., 2005Down) were aligned separately with YFV, and Kokobera virus (KOKV) and Stratford virus (STRV) (Poidinger et al., 1996Down; Nisbet et al., 2005Down) with JEV. Subsequently, identical or nearly identical 3'UTR sequences were excluded from the subgroup alignments by visual inspection. Accordingly, the remaining 3'UTR sequences show maximal virus divergence within each subgroup and they were used for manual construction of a pan-MBFV alignment (see Supplementary Figs S1 and S2 available in JGV Online). Multiple large and small gaps were introduced manually in the alignments using CSs/RCSs as anchors. The viruses used for the alignments are specified by their accession numbers (Supplementary Figs S1 and S2).


   RESULTS AND DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
In published 3'UTR alignments, the homology between all the MBFV was limited to 20–45 nucleotide DRs, i.e. CSs/RCSs and YFV repeats identified as highly conserved regions separated by 50–100 base variable regions. However, available alignment programs do not search for regions of homology that are separated by extended non-homologous regions although, for many RNA viruses, the evolution of UTRs is mediated not only by substitution mutagenesis, but also by deletions and duplications (Copper et al., 1974Down; Faragher & Dalgarno, 1986Down; Bryan et al., 1992Down; Lai, 1992Down; Santagati et al., 1994Down; Pilipenko et al., 1995Down; Hajjou et al., 1996Down; Peerenboom et al., 1997Down; Shi et al., 1997Down; Warren & Murphy, 2003Down).

The results of construction of robust 3'UTR alignments for the MBFV are presented starting with the alignment for YFV and YFV-related viruses that were previously produced using conventional computer programs (Mutebi et al., 2004Down). Our manually corrected alignment revealed a number of important gaps and demonstrated that, even between closely related YFV, 3'UTR evolution involves deletions and insertions rather than substitutions.

3'UTR alignment of YFV subgroup viruses
Supplementary Fig. S1(a) (available in JGV Online) illustrates a 3'UTR alignment between selected strains of the three recognized YFV genotypes, i.e. West African (W), Central/East African (CE) and South American (SA); it also includes three YFV-related viruses, BANZV, UGSV and SEPV (Wang et al., 1996Down; Mutebi et al., 2001Down, 2004Down; Mutebi & Barrett, 2002Down). Although YOKV is also recognized as a YFV-related virus, we included it in the pan-MBFV alignment due to the presence of a poly-AT stretch that aligns only with JEV (see below).

Originally, three 45-nucleotide-long tandem repeats were described for W strains of YFV (Hahn et al., 1987Down). It was subsequently demonstrated that EC and SA strains have lost one (YFV-R2) and two (YFV-R1 and YFV-R2) repeats, respectively (Wang et al., 1996Down; Mutebi et al., 2004Down).

Our alignment (Supplementary Fig. S1a) shows that the process of losing the repeats occurred gradually, with a series of deletions that reduced, but did not completely delete, the repeat elements in different YFV strains. The YFV W strains contain all three intact YFV-specific repeat sequences, whilst the CE strains contain intact YFV-R1/R3 and partial YFV-R2. The SA strains contain intact YFV-R3 and partial YFV-R2/R1 sequences.

The possibility that multiple template switching of RNA polymerase, rather than point mutations, shaped the phenotype of the YFV 3'UTR was supported by comparison with the YFV-related viruses BANZV, UGSV and SEPV (Supplementary Fig. S1a). Previously, it had been reported that they also contained YFV-specific repeats, with two repeats for SEPV, one for UGSV and three for BANZV (Mutebi et al., 2004Down). The alignment in Supplementary Fig. S1(a) illustrates that all three repeats, in a complete or truncated form, were present in all the viruses. However, to identify these repeats, it was necessary to introduce two gaps, a short G1 and a long G2 between YFV-R2 and YFV-R3.

Mutebi et al. (2004)Down suggested that a short sequence of the YFV-repeat, 5'-AACCGGGATACAAC-3', was acquired from the Drosophila melanogaster genome. However, our alignment reveals an alternative possibility, that each repeat was initially formed by multiple duplications and subsequent mutagenesis of an AACCGGG sequence that was initially produced by ‘stuttering’ of the RNA polymerase. Each YFV repeat could have contained four of these initial AACCGGG repeats. The sequences between YFV repeats were also probably part of the longer repeats. For example, there is a TAAAT(T)GG repeat in G1 and G2 only for BANZV, whereas TAAATGG is present in the G2 of other YFV. Thus, TAAAT(T)GG sequences may have been present in the G1 of all YFV, but were subsequently deleted. Another example is the repeat sequence CTGT that is present in all three YFV-related viruses in G2 and also downstream of YFV-R3 (Supplementary Fig. S1a). These examples also support the suggestion that insertional mutagenesis could contribute to the formation of the G2 region.

In order to assess whether or not the YFV repeats are part of longer repeat sequences, they were placed under each other and aligned (Supplementary Fig. S1b available in JGV Online). The TGAAAC sequence that follows YFV-R2 within the G2 region and which is observed only in YFV (not in SEPV, UGSV or BANV) might be relocated from G2 into the G0 region (compare Supplementary Figs S1a and S1b). Less homology is observed in the original position of this sequence within G2, and therefore the ambiguous position of the sequence TGAAAC in either G2 or G0 is indicated (boxed in Supplementary Fig. 1bDown). These examples suggest that YFV repeats were originally longer. Furthermore, the introduction of large gaps enabled us to detect remnants of YFV-R4 that was probably also present in a YFV ancestral lineage.


Figure 1
View larger version (17K):
[in this window]
[in a new window]
 
Fig. 1. Schematic organization of repeats YFV-R1, -R2 and -R3 as deduced from the alignment in Supplementary Fig. S1 within the MBFV 3'UTR; adapted from previous publications (Hahn et al., 1987Down; Mutebi et al., 2004Down) with modifications. Intact YFV repeats are shown as shaded boxes, whereas partial YFV repeats are shown as open boxes. Solid lines depict extended regions of homology whereas dotted lines indicate regions where homology is confined to short nucleotide regions. Gaps between YFV repeats are specified as G0, G1 and G2 and reduced YFV-R4, in correspondence with the alignment in Supplementary Fig. S1.

 
Examination of the alignment downstream of YFV-R3 also indicates that the evolution of YFV and YFV-related viruses occurred primarily through deletions and insertions rather than substitutions. A 216 nucleotide duplication has also been demonstrated in the 3'UTR of numerous South American genotype 1 YFV isolates (Bryant et al., 2005Down). Whilst the origin of CE strains is not known (beyond the phylogenetic prediction), the evolution of SA strains from W strains is assumed to have resulted from the introduction of virus into South America during the slave trade, 100–400 years ago (Gould et al., 2003Down). Therefore, even during this relatively short period, deletions and duplications have substantially changed the 3'UTR of closely related viruses.

Thus, the deduced consensus sequence might represent a longer 3'UTR, of potentially 760 nucleotides, that was present in an ancestral lineage of the YFV and YFV-like viruses. The current longest identified YFV 3'UTR contains 511 nucleotides; the precise length of the 3'UTR of YFV-related viruses has not been determined. On the basis of this alignment and our interpretation of the findings, a schedule for the possible development of the YFV repeats from a common ancestral lineage is presented in Fig. 1Up.

3'UTR alignment of all MBFV subgroups
A comparative 3'UTR alignment between all available MBFV is demonstrated in Supplementary Fig. S2 (available in JGV Online). The CSs and RCSs were used as anchors to search for homology between them. The elements of the predicted RNA secondary structures, superimposed on the alignment, include (i) the 3'LSH with a conserved pentanucleotide CACAG in loop1 (Brinton et al., 1986Down), (ii) a short upstream SL2 that interacts with the 3'LSH in the pseudoknot (Shi et al., 1996Down), (iii) dumbbell-like structures, DB1 (exposing CS2), shared between JEV, DENV and YFV, and DB2 (exposing RCS2), shared between JEV and DENV (Proutski et al., 1997aDown; Olsthoorn & Bol, 2001Down), and (iv) the 3'CYCL sequence, CATATTGA, which interacts with the inverted repeat 5'CYCL sequence in the capsid gene to circularize the virus genome (Hahn et al., 1987Down), required for the initiation of replication. The experimentally established boundaries that define whether or not a virus is viable (Men et al., 1996Down) are also indicated in Supplementary Fig. S2. The observations made from this alignment are described below.

CSs and RCSs are present in all the MBFV.
Although the CSs/RCSs were previously the only regions of recognized homology between all MBFV, the alignment in Supplementary Fig. S2 reveals a number of other shorter regions of homology which we interpret as remnants of longer preserved elements between the CSs and RCSs. YFV sequence was similar in some regions to the DENV or JEV subgroup outside the CSs, but not necessarily with both virus groups in the same regions. A schematic illustration of the presence of CSs/RCSs in all MBFV is presented in Fig. 2Down.


Figure 2
View larger version (33K):
[in this window]
[in a new window]
 
Fig. 2. Scheme for the development of the CSs and RCSs of the MBFV deduced from the alignments in Supplementary Fig. S2. Regions of high nucleotide conservation (CS, RCS and YFV repeats) are boxed and regions of nucleotide variability between them are depicted as solid lines. Shaded boxes specify intact CSs, RCSs, YFV repeats and poly-AT, whereas partial conservation is depicted as open boxes. Schematically presented RNA secondary structures are depicted against the corresponding regions of the 3'UTR and are specified appropriately.

 
Notably, YOKV, which does not have a recognized arthropod vector but diverged from the mosquito-borne YFV-related virus lineages (Kuno et al., 1998Down; Gould et al., 2003Down; Tajima et al., 2005Down; Cook & Holmes, 2006Down), showed significant homology with YFV and the YFV-related viruses which was identified following the introduction of multiple deletions and insertions. Extensive deletions were also apparent in KOKV and STRV (Supplementary Fig. S2), both of which evolved from JEV-related evolutionary lineages (Poidinger et al., 1996Down; Kuno et al., 1998Down; Gould et al., 2001Down; Nisbet et al., 2005Down).

Fragments of the YFV repeat sequences are present in other MBFV.
In previous publications, no homology between YFV repeats and other MBFV was demonstrated. We examined this in more detail by introducing multiple gaps manually into the alignment. However, only short remnants of homologous sequences of each YFV repeat were found in the DENV and JEV subgroups (Supplementary Fig. S2). The region between YFV-R1 and YFV-R2 also contained traces of similarity between the different MBFV subgroups (Supplementary Fig. S2). The alignment in this region is not robust, but the short regions of homology could still represent the remnants of ancient sequences common to all MBFV in accordance with their monophyletic origin (Gould et al., 2001Down). The most likely scenario is that the three YFV repeats were originally present in an early lineage and then subsequently evolved by deletions, insertions and substitutions into the JEV, DENV and YFV subgroups (Fig. 2Up). However, the YFV repeats were not reproduced robustly in the MBFV subgroups, and therefore they may be unique to the YFV subgroup.

Was the poly-AT region originally present in all the MBFV?
The nucleotide region between the YFV stop codon and the YFV-R1 is AT-enriched and is referred to as the poly-AT region. This region appears to have evolved as the result of multiple stuttering by the viral RNA polymerase around the TAA stop codon. West Nile virus (WNV) and Murray Valley encephalitis virus (MVEV) have the longest AT-enriched regions whereas, in the DENV and YFV subgroups, this region is much shorter, suggesting that AT-enrichment occurred only in the JEV subgroup. However, the equivalent region of YOKV aligns quite closely with JEV, although beyond the poly-AT region YOKV aligns closely with YFV, implying that the precursor lineage of YFV may have had a fully developed AT-enriched domain. In order to illustrate this, the YOKV 3'UTR sequence is placed next to both YFV and JEV in the alignment (Supplementary Fig. S2).

The preservation of the poly-AT region in JEV may reflect a selective advantage in the natural habitat of this virus. In structural and functional terms, the AT-enriched region could promote the formation of relaxed, i.e. unstable secondary structures. This is supported by the observation that the region of the 3'UTR proximal to the stop codon of all flaviviruses is extremely variable and in the laboratory is readily altered by spontaneous deletions (Mandl et al., 1998Down; Mutebi et al., 2004Down). Moreover, engineered deletion of this region did not prevent the production of infectious virus (Men et al., 1996Down; Mandl et al., 1998Down). Computer-predicted folding using MFold (available at http://mfold.burnet.edu.au/) demonstrated the formation of only short AT- and GT-paired stem–loop structures (not shown). The advantage for the virus to preserve these new ‘immature’ sequences might be that they do not form long stable stems and might serve as spacers that separate the ORF and folded distal 3'UTR essential for the assembly of the replication complexes (RC). The folding of YFV repeats as two short stem–loop structures (Mutebi et al., 2004Down) demonstrated a similar possibility.

Superimposition of secondary RNA structures on the MBFV alignment.
Another important observation from Supplementary Fig. S2 is the compatibility between the positions of RNA secondary structures predicted for each subgroup of the MBFV. Although these viruses are divergent, the sequences that form the 3'LSH, SL2, DB1 and DB2 for individual groups were located in equivalent regions of the 3'UTR for each of the MBFV subgroups (Supplementary Fig. S2). Firstly, this could be interpreted as evidence (although indirect) of the robustness of the alignment and, secondly, it implies that the secondary RNA structures in the most divergent MBFV lineages could have been formed by similar nucleotide sequences. Each of these RNA structures exposes the CSs/RCSs, the preservation of which might be required for an essential signalling function for RNA replication. For example, the 3'LSH exposes the pentanucleotide CACAG that is conserved between vector-borne flaviviruses (Charlier et al., 2002Down); an engineered mutation in this region abolishes or impairs virus replication (Gritsun et al., 2001Down; Khromykh et al., 2003Down; Tilgner et al., 2005Down). Similarly, DB1 and DB2 expose CS2 and RCS2 (Olsthoorn & Bol, 2001Down), the significance of which is not obvious since, despite their conservation, they are not essential for viability of experimentally engineered viruses, although mutant viruses had reduced infectivity titres (Men et al., 1996Down; Proutski et al., 1999Down; Bredenbeek et al., 2003Down).

The CS3 and RCS3 that are located in close proximity to the stop codon of JEV also form two identical short stem–loop structures (Fig. 2Up). The deletion of CS3 and RCS3 did not lead to the loss of virus infectivity, although reduced virus replication was observed (Khromykh & Westaway, 1997Down; Lo et al., 2003Down; Tilgner et al., 2005Down).

What is the role of the CSs and RCSs?
Previously, it was proposed that the 3'UTR consists of a promoter region that is essential for virus viability and an enhancer region that is not essential for virus viability but increases the RNA replication rate (Proutski et al., 1999Down). The experimentally established boundary that determines whether or not the virus is viable (Men et al., 1996Down) could correspond to the boundary between the promoter and enhancer region (Supplementary Fig. S2). Thus, the minimal promoter function in MBFV could be provided by 3'LSH–SL2–CS1; the engineered virus with a completely deleted enhancer region is infectious, but only in mosquito (not in vertebrate) cells (Men et al., 1996Down). Accordingly, the virus replication enhancer includes all DRs (CSs/RCSs and YFV repeats) that probably facilitate virus replication in vertebrate cells.

Experiments that demonstrate a non-essential role of the DRs in virus viability and RNA synthesis (Men et al., 1996Down; Khromykh & Westaway, 1997Down; Bredenbeek et al., 2003Down; Lo et al., 2003Down; Tilgner et al., 2005Down) appear to contradict the concept of conservation of the CSs/RCSs (Hahn et al., 1987Down; Mutebi et al., 2004Down). However, laboratory experiments may not adequately reflect the essential requirements of DRs for flaviviruses in nature. Since all engineered MBFV exhibit reduced RNA replication rates, compared with wild-type viruses, it is possible that a high RNA replication rate might exert a selection pressure on MBFV. Mosquitoes feed on vertebrate hosts for no more than a few seconds; consequently, the frequency of viral transmission is increased if the virus rapidly reaches a high titre in the blood, before the immune system reduces viraemia. Indeed, the engineered deletions of CSs and RCSs decrease the replication rate of DENV4 to a greater extent in simian than in mosquito cells (Men et al., 1996Down). This corresponds with the proposition that the virus enhancer has evolved to meet the requirements of the short viraemic period in the vertebrate host. In contrast, infected mosquitoes have a period of weeks to amplify the virus before the next feeding period, and therefore the enhancement function is less essential for virus replication in mosquitoes than in vertebrates.

In molecular terms, how could the enhancer accelerate virus replication? It is recognized that a wide variety of host proteins, with functions not directly related to cellular RNA synthesis, are frequently associated with replication of positive-strand RNA viruses (reviewed by Lai, 1998Down; Boguszewska-Chachulska & Haenni, 2005Down), including flaviviruses (Blackwell & Brinton, 1995Down, 1997Down; Ta & Vrati, 2000Down; De Nova Ocampo et al., 2002Down). These are often identified as proteins with RNA-binding and trafficking activities. Indeed, the individual steps of the flavivirus life cycle, i.e. translation, RNA synthesis and assembly, take place in a variety of cellular compartments (convoluted membranes and vesicular packets) that proliferate in cells in response to virus infection (Mackenzie et al., 2001Down; Westaway et al., 2002Down, 2003Down; Uchil & Satchidanandam, 2003Down). The DRs might interact with cellular/viral proteins that direct viral RNA from the sites of translation into the replication sites where assembly of viral RCs is most efficient. Alternatively, the CSs could be direct/indirect signals for cells to express factors essential for efficient viral RNA synthesis. In the absence of DRs, flavivirus replication could be initiated at low efficiency because of the leakage of viral RNA into the appropriate cell compartment or because of the absence of sufficient levels of cellular components that accelerate the assembly of the viral RC.

The particular signalling function of DRs has presumably evolved in different flaviviruses according to their adaptation to different invertebrate and vertebrate hosts. In other words, DRs as enhancers could be the means by which virus adapts to the molecular environment of a new host to maximize replication and therefore increase transmission efficiency. This would explain why some DRs are shared between MBFV subgroups, whilst others are different. It could also explain the apparent discrepancy between conservation of the CSs and RCSs and their apparent non-essential function in laboratory-maintained viruses.

The reason for preservation of the DRs as duplicated signals is less clear. It is possible that host/viral proteins involved in RNA synthesis are active as dimers and therefore the presence of a double signal would accelerate assembly of the RC. The effect of each individual element of the replication enhancer could be relatively moderate and insignificant for laboratory-maintained viruses, but the cumulative action of DRs might accelerate virus replication and this could be a vital requirement for efficient virus transmission and dissemination in the more fastidious natural environment.

If DRs do provide a signalling function, a plausible explanation can now be offered for the preservation of YFV repeats in W strains and their loss in CE and SA strains; their signal function as triple DRs may be essential only in the particular hosts (mosquitoes and vertebrate) in West Africa. Once virus escapes this environment, the need for a triple signal is abolished as postulated for CE and SA strains during their adaptation to new species.

Why is CS1 conserved?
The exposed positions of the DRs in the secondary RNA structure of MBFV are compatible with a proposed signalling function for these elements and also provide a possible reason for why they are conserved. However, the reason for CS1 conservation is less clear. CS1 is recognized as the most conserved region within the MBFV, but the alignment (Supplementary Fig. S2) shows that this region is not absolutely conserved. Nucleotide differences in CS1 were revealed not only between the MBFV subgroups but also within closely related viruses such as MVEV and other viruses from the JEV subgroup (Supplementary Fig. S2). Secondly, the overlapping of two structural models, i.e. the linear (Proutski et al., 1997aDown, bDown, 1999Down) and the cyclization model (Khromykh et al., 2001Down; Thurner et al., 2004Down), within CS1 enabled us to conclude that CS1 may not represent one continuous domain. In fact, linear CS1 could be subdivided into two domains. The first was already described, i.e. 3'CYCL, which is absolutely conserved between all MBFV and has been proven experimentally to interact with 5'CYCL to initiate RNA replication (You & Padmanabhan, 1999Down; Ackermann & Padmanabhan, 2001Down; Charrel et al., 2001Down; Khromykh et al., 2001Down; You et al., 2001Down; Corver et al., 2003Down; Lo et al., 2003Down). However, for RNA–RNA interactions to occur, there is no requirement for absolute preservation of the sequence, and it was suggested earlier that 3'CYCL, as either a single-stranded or double-stranded sequence, might have a signalling function, not only for replication but also for virion assembly (Proutski et al., 1999Down). Recently, it was also suggested that, in addition to interacting with 5'CYCL, 3'CYCL might interact with the 5'UTR (possibly mediated by cellular proteins) to initiate translation (Chiu et al., 2005Down). It is also noteworthy that the pseudoknot PK1, which is conserved in the JEV and DENV subgroups (Olsthoorn & Bol, 2001Down), overlaps with the CS1 sequence (Supplementary Fig. S2).

The second domain of the CS1, downstream of the 3'CYCL, is less conserved and has not been predicted to be involved in the formation of a well-defined secondary structure. However, sequences of several secondary structures overlap significantly within this domain. One of these is SL2, with a loop that exposes a pseudoknot-interacting sequence (Fig. 2Up). Additionally, in the cyclization model, some sequences of the CS1 beyond 3'CYCL are involved in the formation of a long stem (interacting with the 5'UTR and capsid gene) that is essential to initiate virus replication, although these sequences are not identical within the three MBFV subgroups (Supplementary Fig. S2).

Thus, CS1 probably does not represent an entire entity and its conservation might be explained by the overlapping of several essential linear and structural features. However, a signalling function for CS1, which is not revealed by current RNA folding models, cannot be totally excluded.

How did CS and RCS emerge?
What were the molecular events that resulted in the appearance of 20- to 45-nucleotide repeats separated by longer (~50–130 nucleotides) non-repeated sequences? It is possible that the DRs within the 3'UTR of the MBFV were initially longer, but they were subsequently reduced to short DRs (CSs/RCSs and YFV repeats) by extensive mutagenesis. Therefore, these short DRs could represent remnants of longer tandem repeats that have gradually been deleted and mutated. However, our alignment did not reveal similar remnants outside the CSs, probably because the MBFV subgroups are so divergent that no traces remain.

Nevertheless, we recently demonstrated that formation of the 3'UTR of TBFV was mediated by multiple duplication of long regions of about 200 nucleotides (Gritsun & Gould, 2006bDown); the preservation of these long DRs might be explained by the slower rate of evolution of TBFV in comparison with MBFV (Gritsun et al., 1995Down; Zanotto et al., 1996Down). We also demonstrated that formation of the 3'UTR of NCFV resulted from the duplication of a 600 nucleotide-long sequence (Gritsun & Gould, 2006aDown). Duplication of a region of 216 nucleotides in the 3'UTR of South American YFV isolates effectively doubled the secondary RNA structures and recreated RCS2, not present in ‘conventional’ YFV strains (Bryant et al., 2005Down). Finally, we have recently demonstrated the possible common origin of conserved DRs in the genomes of all flaviviruses, including the most divergent species. We believe that this can be explained by development of their 3'UTRs from a common precursor lineage (manuscript in preparation).

Conclusions
In this paper, we have described the construction and analysis of a robust comparative alignment between the 3'UTRs of three divergent MBFV subgroups, the JEV-, DENV- and YFV-related viruses. Conserved DRs (CSs/RCSs and YFV repeats) were revealed either intact or as remnants in each of the three MBFV subgroups, suggesting their common origin. We conclude that, at the molecular level, DRs function as replication enhancers selected under the constraints of the particular transmission and dissemination imposed by the hosts. Appropriate genetic modification of MBFV DRs and studies of their transmission efficiency would clarify whether or not DRs were preserved as essential elements for virus survival in the environment.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS AND DISCUSSION
REFERENCES
 
Ackermann, M. & Padmanabhan, R. (2001). De novo synthesis of RNA by the dengue virus RNA-dependent RNA polymerase exhibits temperature dependence at the initiation but not elongation phase. J Biol Chem 276, 39926–39937.[Abstract/Free Full Text]

Blackwell, J. L. & Brinton, M. A. (1995). BHK cell proteins that bind to the 3' stem-loop structure of the West Nile virus genome RNA. J Virol 69, 5650–5658.[Abstract]

Blackwell, J. L. & Brinton, M. A. (1997). Translation elongation factor-1 alpha interacts with the 3' stem-loop region of West Nile virus genomic RNA. J Virol 71, 6433–6444.[Abstract]

Boguszewska-Chachulska, A. M. & Haenni, A. L. (2005). RNA viruses redirect host factors to better amplify their genome. Adv Virus Res 65, 29–61.[Medline]

Bredenbeek, P. J., Kooi, E. A., Lindenbach, B., Huijkman, N., Rice, C. M. & Spaan, W. J. (2003). A stable full-length yellow fever virus cDNA clone and the role of conserved RNA elements in flavivirus replication. J Gen Virol 84, 1261–1268.[Abstract/Free Full Text]

Brinton, M. A., Fernandez, A. V. & Dispoto, J. H. (1986). The 3'-nucleotides of flavivirus genomic RNA form a conserved secondary structure. Virology 153, 113–121.[CrossRef][Medline]

Bryan, G. T., Gardner, R. C. & Forster, R. L. (1992). Nucleotide sequence of the coat protein gene of a strain of clover yellow vein virus from New Zealand: conservation of a stem-loop structure in the 3' region of potyviruses. Arch Virol 124, 133–146.[CrossRef][Medline]

Bryant, J. E., Vasconcelos, P. F., Rijnbrand, R. C., Mutebi, J. P., Higgs, S. & Barrett, A. D. (2005). Size heterogeneity in the 3' noncoding region of South American isolates of yellow fever virus. J Virol 79, 3807–3821.[Abstract/Free Full Text]

Charlier, N., Leyssen, P., Pleij, C. W. A. & 7 other authors (2002). Complete genome sequence of Montana Myotis leukoencephalitis virus, phylogenetic analysis and comparative study of the 3' untranslated region of flaviviruses with no known vector. J Gen Virol 83, 1875–1885.[Abstract/Free Full Text]

Charrel, R. N., Zaki, A. M., Attoui, H. & 7 other authors (2001). Complete coding sequence of the Alkhurma virus, a tick-borne flavivirus causing severe hemorrhagic fever in humans in Saudi Arabia. Biochem Biophys Res Commun 287, 455–461.[CrossRef][Medline]

Chiu, W. W., Kinney, R. M. & Dreher, T. W. (2005). Control of translation by the 5'- and 3'-terminal regions of the dengue virus genome. J Virol 79, 8303–8315.[Abstract/Free Full Text]

Cook, S. & Holmes, E. C. (2006). A multigene analysis of the phylogenetic relationships among the flaviviruses (Family: Flaviviridae) and the evolution of vector transmission. Arch Virol 151, 309–325.[CrossRef][Medline]

Copper, P. D., Steiner-Pryor, A., Scotti, P. D. & Delong, D. (1974). On the nature of poliovirus genetic recombinants. J Gen Virol 23, 41–49.[Abstract/Free Full Text]

Corver, J., Lenches, E., Smith, K., Robison, R. A., Sando, T., Strauss, E. G. & Strauss, J. H. (2003). Fine mapping of a cis-acting sequence element in Yellow fever virus RNA that is required for RNA replication and cyclization. J Virol 77, 2265–2270.[Abstract/Free Full Text]

De Nova Ocampo, M., Villegas Sepulveda, N. & del Angel, R. M. (2002). Translation elongation factor-1alpha, La, and PTB interact with the 3' untranslated region of dengue 4 virus RNA. Virology 295, 337–347.[CrossRef][Medline]

Edgil, D., Diamond, M. S., Holden, K. L., Paranjape, S. M. & Harris, E. (2003). Translation efficiency determines differences in cellular infection among dengue virus type 2 strains. Virology 317, 275–290.[CrossRef][Medline]

Faragher, S. G. & Dalgarno, L. (1986). Regions of conservation and divergence in the 3' untranslated sequences of genomic RNA from Ross River virus isolates. J Mol Biol 190, 141–148.[CrossRef][Medline]

Gould, E. A., de Lamballerie, X., Zanotto, P. M. & Holmes, E. C. (2001). Evolution, epidemiology, and dispersal of flaviviruses revealed by molecular phylogenies. Adv Virus Res 57, 71–103.[Medline]

Gould, E. A., de Lamballerie, X., Zanotto, P. M. & Holmes, E. C. (2003). Origins, evolution, and vector/host coadaptations within the genus Flavivirus. Adv Virus Res 59, 277–314.[Medline]

Gritsun, T. S. & Gould, E. A. (2006a). The 3' untranslated regions of Kamiti River virus and Cell fusing agent virus originated by self-duplication. J Gen Virol 87, 2615–2619.[Abstract/Free Full Text]

Gritsun, T. S. & Gould, E. A. (2006b). The 3' untranslated region of tick-borne flaviviruses originated by the duplication of long repeat sequences within the open reading frame. Virology 350, 269–275.[CrossRef][Medline]

Gritsun, T. S., Holmes, E. C. & Gould, E. A. (1995). Analysis of flavivirus envelope proteins reveals variable domains that reflect their antigenicity and may determine their pathogenesis. Virus Res 35, 307–321.[CrossRef][Medline]

Gritsun, T. S., Desai, A. & Gould, E. A. (2001). The degree of attenuation of tick-borne encephalitis virus depends on the cumulative effects of point mutations. J Gen Virol 82, 1667–1675.[Abstract/Free Full Text]

Hahn, C. S., Hahn, Y. S., Rice, C. M., Lee, E., Dalgarno, L., Strauss, E. G. & Strauss, J. H. (1987). Conserved elements in the 3' untranslated region of flavivirus RNAs and potential cyclization sequences. J Mol Biol 198, 33–41.[CrossRef][Medline]

Hajjou, M., Hill, K. R., Subramaniam, S. V., Hu, J. Y. & Raju, R. (1996). Nonhomologous RNA-RNA recombination events at the 3' nontranslated region of the Sindbis virus genome: hot spots and utilization of nonviral sequences. J Virol 70, 5153–5164.[Abstract/Free Full Text]

Heinz, F. X., Collett, M. S., Purcell, R. H., Gould, E. A., Howard, C. R., Houghton, M., Moormann, R. J. M., Rice, C. M. & Thiel, H. J. (2000). Family Flaviviridae. In Virus Taxonomy. Seventh Report of the International Committee for the Taxonomy of Viruses, pp. 859–878. Edited by M. H. V. van Regenmortel, C. M. Fauquet, D. H. L. Bishop, E. Carstens, M. K. Estes, S. Lemon, J. Maniloff, M. A. Mayo, D. McGeoch, C. R. Pringle & R. B. Wickner. San Diego: Academic Press.

Holden, K. L. & Harris, E. (2004). Enhancement of dengue virus translation: role of the 3' untranslated region and the terminal 3' stem-loop domain. Virology 329, 119–133.[CrossRef][Medline]

Khromykh, A. A. & Westaway, E. G. (1997). Subgenomic replicons of the flavivirus Kunjin: construction and applications. J Virol 71, 1497–1505.[Abstract]

Khromykh, A. A., Meka, H., Guyatt, K. J. & Westaway, E. G. (2001). Essential role of cyclization sequences in flavivirus RNA replication. J Virol 75, 6719–6728.[Abstract/Free Full Text]

Khromykh, A. A., Kondratieva, N., Sgro, J. Y., Palmenberg, A. & Westaway, E. G. (2003). Significance in replication of the terminal nucleotides of the flavivirus genome. J Virol 77, 10623–10629.[Abstract/Free Full Text]

Kuno, G., Chang, G. J., Tsuchiya, K. R., Karabatsos, N. & Cropp, C. B. (1998). Phylogeny of the genus Flavivirus. J Virol 72, 73–83.[Abstract/Free Full Text]

Lai, M. M. (1992). RNA recombination in animal and plant viruses. Microbiol Rev 56, 61–79.[Abstract/Free Full Text]

Lai, M. M. (1998). Cellular factors in the transcription and replication of viral RNA genomes: a parallel to DNA-dependent RNA transcription. Virology 244, 1–12.[CrossRef][Medline]

Li, W. & Brinton, M. A. (2001). The 3' stem loop of the West Nile virus genomic RNA can suppress translation of chimeric mRNAs. Virology 287, 49–61.[CrossRef][Medline]

Lindenbach, B. D. & Rice, C. M. (2001). Flaviviridae: the viruses and their replication. In Fields Virology, 4th edn, pp. 991–1042. Edited by D. M. Knipe & P. M. Howley. London, New York & Tokyo: Lippincott Williams & Wilkins.

Lo, M. K., Tilgner, M., Bernard, K. A. & Shi, P. Y. (2003). Functional analysis of mosquito-borne flavivirus conserved sequence elements within 3' untranslated region of West Nile virus by use of a reporting replicon that differentiates between viral translation and RNA replication. J Virol 77, 10004–10014.[Abstract/Free Full Text]

Mackenzie, J. M., Khromykh, A. A. & Westaway, E. G. (2001). Stable expression of noncytopathic Kunjin replicons simulates both ultrastructural and biochemical characteristics observed during replication of Kunjin virus. Virology 279, 161–172.[CrossRef][Medline]

Mandl, C. W., Holzmann, H., Meixner, T., Rauscher, S., Stadler, P. F., Allison, S. L. & Heinz, F. X. (1998). Spontaneous and engineered deletions in the 3' noncoding region of tick-borne encephalitis virus: construction of highly attenuated mutants of a flavivirus. J Virol 72, 2132–2140.[Abstract/Free Full Text]

Men, R., Bray, M., Clark, D., Chanock, R. M. & Lai, C. J. (1996). Dengue type 4 virus mutants containing deletions in the 3' noncoding region of the RNA genome: analysis of growth restriction in cell culture and altered viremia pattern and immunogenicity in rhesus monkeys. J Virol 70, 3930–3937.[Abstract]

Mutebi, J. P. & Barrett, A. D. (2002). The epidemiology of yellow fever in Africa. Microbes Infect 4, 1459–1468.[CrossRef][Medline]

Mutebi, J. P., Wang, H., Li, L., Bryant, J. E. & Barrett, A. D. (2001). Phylogenetic and evolutionary relationships among yellow fever virus isolates in Africa. J Virol 75, 6999–7008.[Abstract/Free Full Text]

Mutebi, J. P., Rijnbrand, R. C., Wang, H., Ryman, K. D., Wang, E., Fulop, L. D., Titball, R. & Barrett, A. D. (2004). Genetic relationships and evolution of genotypes of yellow fever virus and other members of the yellow fever virus group within the Flavivirus genus based on the 3' noncoding region. J Virol 78, 9652–9665.[Abstract/Free Full Text]

Nisbet, D. J., Lee, K. J., van den Hurk, A. F., Johansen, C. A., Kuno, G., Chang, G. J., Mackenzie, J. S., Ritchie, S. A. & Hall, R. A. (2005). Identification of new flaviviruses in the Kokobera virus complex. J Gen Virol 86, 121–124.[Abstract/Free Full Text]

Olsthoorn, R. C. & Bol, J. F. (2001). Sequence comparison and secondary structure analysis of the 3' noncoding region of flavivirus genomes reveals multiple pseudoknots. RNA 7, 1370–1377.[Abstract]

Peerenboom, E., Jacobi, V., Cartwright, E. J., Adams, M. J., Steinbiss, H. H. & Antoniw, J. F. (1997). A large duplication in the 3'-untranslated region of a subpopulation of RNA2 of the UK-M isolate of barley mild mosaic bymovirus. Virus Res 47, 1–6.[CrossRef][Medline]

Pilipenko, E. V., Gmyl, A. P. & Agol, V. I. (1995). A model for rearrangements in RNA genomes. Nucleic Acids Res 23, 1870–1875.[Abstract/Free Full Text]

Pletnev, A. G. (2001). Infectious cDNA clone of attenuated Langat tick-borne flavivirus (strain E5) and a 3' deletion mutant constructed from it exhibits decreased neuroinvasiveness in immunodeficient mice. Virology 282, 288–300.[CrossRef][Medline]

Poidinger, M., Hall, R. A. & Mackenzie, J. S. (1996). Molecular characterization of the Japanese encephalitis serocomplex of the flavivirus genus. Virology 218, 417–421.[CrossRef][Medline]

Proutski, V., Gould, E. A. & Holmes, E. C. (1997a). Secondary structure of the 3' untranslated region of flaviviruses: similarities and differences. Nucleic Acids Res 25, 1194–1202.[Abstract/Free Full Text]

Proutski, V., Gaunt, M. W., Gould, E. A. & Holmes, E. C. (1997b). Secondary structure of the 3'-untranslated region of yellow fever virus: implications for virulence, attenuation and vaccine development. J Gen Virol 78, 1543–1549.[Abstract]

Proutski, V., Gritsun, T. S., Gould, E. A. & Holmes, E. C. (1999). Biological consequences of deletions within the 3'-untranslated region of flaviviruses may be due to rearrangements of RNA secondary structure. Virus Res 64, 107–123.[CrossRef][Medline]

Santagati, M. G., Itaranta, P. V., Koskimies, P. R., Maatta, J. A., Salmi, A. A. & Hinkkanen, A. E. (1994). Multiple repeating motifs are found in the 3'-terminal non-translated region of Semliki Forest virus A7 variant genome. J Gen Virol 75, 1499–1504.[Abstract/Free Full Text]

Shi, P. Y., Brinton, M. A., Veal, J. M., Zhong, Y. Y. & Wilson, W. D. (1996). Evidence for the existence of a pseudoknot structure at the 3' terminus of the flavivirus genomic RNA. Biochemistry 35, 4222–4230.[CrossRef][Medline]

Shi, B., Ding, S. & Symons, R. H. (1997). Two novel subgenomic RNAs derived from RNA 3 of tomato aspermy cucumovirus. J Gen Virol 78, 505–510.[Abstract]

Ta, M. & Vrati, S. (2000). Mov34 protein from mouse brain interacts with the 3' noncoding region of Japanese encephalitis virus. J Virol 74, 5108–5115.[Abstract/Free Full Text]

Tajima, S., Takasaki, T., Matsuno, S., Nakayama, M. & Kurane, I. (2005). Genetic characterization of Yokose virus, a flavivirus isolated from the bat in Japan. Virology 332, 38–44.[CrossRef][Medline]

Thurner, C., Witwer, C., Hofacker, I. L. & Stadler, P. F. (2004). Conserved RNA secondary structures in Flaviviridae genomes. J Gen Virol 85, 1113–1124.[Abstract/Free Full Text]

Tilgner, M., Deas, T. S. & Shi, P. Y. (2005). The flavivirus-conserved penta-nucleotide in the 3' stem-loop of the West Nile virus genome requires a specific sequence and structure for RNA synthesis, but not for viral translation. Virology 331, 375–386.[CrossRef][Medline]

Uchil, P. D. & Satchidanandam, V. (2003). Architecture of the flaviviral replication complex. Protease, nuclease, and detergents reveal encasement within double-layered membrane compartments. J Biol Chem 278, 24388–24398.[Abstract/Free Full Text]

Wang, E., Weaver, S. C., Shope, R. E., Tesh, R. B., Watts, D. M. & Barrett, A. D. (1996). Genetic variation in yellow fever virus: duplication in the 3' noncoding region of strains from Africa. Virology 225, 274–281.[CrossRef][Medline]

Warren, C. E. & Murphy, J. F. (2003). The complete nucleotide sequence of Pepper mottle virus-Florida RNA. Arch Virol 148, 189–197.[CrossRef][Medline]

Westaway, E. G., Mackenzie, J. M. & Khromykh, A. A. (2002). Replication and gene function in Kunjin virus. Curr Top Microbiol Immunol 267, 323–351.[Medline]

Westaway, E. G., Mackenzie, J. M. & Khromykh, A. A. (2003). Kunjin RNA replication and applications of Kunjin replicons. Adv Virus Res 59, 99–140.[CrossRef][Medline]

You, S. & Padmanabhan, R. (1999). A novel in vitro replication system for Dengue virus. Initiation of RNA synthesis at the 3'-end of exogenous viral RNA templates requires 5'- and 3'-terminal complementary sequence motifs of the viral RNA. J Biol Chem 274, 33714–33722.[Abstract/Free Full Text]

You, S., Falgout, B., Markoff, L. & Padmanabhan, R. (2001). In vitro RNA synthesis from exogenous dengue viral RNA templates requires long range interactions between 5'- and 3'-terminal regions that influence RNA structure. J Biol Chem 276, 15581–15591.[Abstract/Free Full Text]

Zanotto, P. M., Gould, E. A., Gao, G. F., Harvey, P. H. & Holmes, E. C. (1996). Population dynamics of flaviviruses revealed by molecular phylogenies. Proc Natl Acad Sci U S A 93, 548–553.[Abstract/Free Full Text]

Received 25 May 2006; accepted 14 July 2006.



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Figures
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Gritsun, T. S.
Right arrow Articles by Gould, E. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Gritsun, T. S.
Right arrow Articles by Gould, E. A.
Agricola
Right arrow Articles by Gritsun, T. S.
Right arrow Articles by Gould, E. A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS