Molecular basis for naturally occurring elevated readthrough transcription across the M-F junction of the paramyxovirus SV5

Analysis of Paramyxovirus Transcription and Replication by High-Throughput Sequencing

High-throughput sequencing (HTS) of virus-infected cells can be used to study in great detail the patterns of virus transcription and replication. For paramyxoviruses, and by analogy for all other negative-strand RNA viruses, we show that directional sequencing must be used to distinguish between genomic RNA and mRNA/antigenomic RNA because significant amounts of genomic RNA copurify with poly(A)-selected mRNA. We found that the best method is directional sequencing of total cell RNA, after the physical removal of rRNA (and mitochondrial RNA), because quantitative information on the abundance of both genomic RNA and mRNA/antigenomes can be simultaneously derived. Using this approach, we revealed new details of the kinetics of virus transcription and replication for parainfluenza virus (PIV) type 2, PIV3, PIV5, and mumps virus, as well as on the relative abundance of the individual viral mRNAs.

We have developed a high-throughput sequencing (HTS) workflow for investigating paramyxovirus transcription and replication. We show that sequencing of oligo(dT)-selected polyadenylated mRNAs, without considering the orientation of the RNAs from which they had been generated, cannot accurately be used to analyze the abundance of viral mRNAs because genomic RNA copurifies with the viral mRNAs. The best method is directional sequencing of infected cell RNA that has physically been depleted of ribosomal and mitochondrial RNA followed by bioinformatic steps to differentiate data originating from genomes from viral mRNAs and antigenomes. This approach has the advantage that the abundance of viral mRNA (and antigenomes) and genomes can be analyzed and quantified from the same data. We investigated the kinetics of viral transcription and replication during infection of A549 cells with parainfluenza virus type 2 (PIV2), PIV3, PIV5, or mumps virus and determined the abundances of individual viral mRNAs and readthrough mRNAs. We found that the mRNA abundance gradients differed significantly between all four viruses but that for each virus the pattern remained relatively stable throughout infection. We suggest that rapid degradation of non-poly(A) mRNAs may be primarily responsible for the shape of the mRNA abundance gradient in parainfluenza virus 3, whereas a combination of this factor and disengagement of RNA polymerase at intergenic sequences, particularly those at the NP:P and P:M gene boundaries, may be responsible in the other viruses.

IMPORTANCE High-throughput sequencing (HTS) of virus-infected cells can be used to study in great detail the patterns of virus transcription and replication. For paramyxoviruses, and by analogy for all other negative-strand RNA viruses, we show that directional sequencing must be used to distinguish between genomic RNA and mRNA/antigenomic RNA because significant amounts of genomic RNA copurify with poly(A)-selected mRNA. We found that the best method is directional sequencing of total cell RNA, after the physical removal of rRNA (and mitochondrial RNA), because quantitative information on the abundance of both genomic RNA and mRNA/antigenomes can be simultaneously derived. Using this approach, we revealed new details of the kinetics of virus transcription and replication for parainfluenza virus (PIV) type 2, PIV3, PIV5, and mumps virus, as well as on the relative abundance of the individual viral mRNAs.

The aberrant gene-end transcription signal of the matrix M gene of human parainfluenza virus type 3 downregulates fusion F protein expression and the F-specific antibody response in vivo

Human parainfluenza virus type 3 (HPIV3), a paramyxovirus, is a major viral cause of severe lower respiratory tract disease in infants and children. The gene-end (GE) transcription signal of the HPIV3 matrix (M) protein gene is identical to those of the nucleoprotein and phosphoprotein genes except that it contains an apparent 8-nucleotide insert. This was associated with an increased synthesis of a readthrough transcript of the M gene and the downstream fusion (F) protein gene. We hypothesized that this insert may function to downregulate the expression of F protein by interfering with termination/reinitiation at the M-F gene junction, thus promoting the production of M-F readthrough mRNA at the expense of monocistronic F mRNA. To test this hypothesis, two similar recombinant HPIV3 viruses from which this insert in the M-GE signal was removed were generated. The M-GE mutants exhibited a reduction in M-F readthrough mRNA and an increase in monocistronic F mRNA. This resulted in a substantial increase in F protein synthesis in infected cells as well as enhanced incorporation of F protein into virions. The efficiency of mutant virus replication was similar to that of wild-type (wt) HPIV3 both in vitro and in vivo. However, the F-protein-specific serum antibody response in hamsters was increased for the mutants compared to wt HPIV3. This study identifies a previously undescribed viral mechanism for attenuating the host adaptive immune response. Repairing the M-GE signal should provide a means to increase the antibody response to a live attenuated HPIV3 vaccine without affecting viral replication and attenuation.

IMPORTANCE

The HPIV3 M-GE signal was previously shown to contain an apparent 8-nucleotide insert that was associated with increased synthesis of a readthrough mRNA of the M gene and the downstream F gene. However, whether this had any significant effect on the synthesis of monocistronic F mRNA or F protein, virus replication, virion morphogenesis, and immunogenicity was unknown. Here, we show that the removal of this insert shifts F gene transcription from readthrough M-F mRNA to monocistronic F mRNA. This resulted in a substantial increase in the amount of F protein expressed in the cell and packaged in the virus particle. This did not affect virus replication but increased the F-specific antibody response in hamsters. Thus, in wild-type HPIV3, the aberrant M-GE signal operates a previously undescribed mechanism that reduces the expression of a major neutralization and protective antigen, resulting in reduced immunogenicity. This has implications for the design of live attenuated HPIV3 vaccines; specifically, the antibody response against F can be elevated by "repairing" the M-GE signal to achieve higher-level F antigen expression, with no effect on attenuation.

Analysis of the highly diverse gene borders in Ebola virus reveals a distinct mechanism of transcriptional regulation

Ebola virus (EBOV) belongs to the group of nonsegmented negative-sense RNA viruses. The seven EBOV genes are separated by variable gene borders, including short (4- or 5-nucleotide) intergenic regions (IRs), a single long (144-nucleotide) IR, and gene overlaps, where the neighboring gene end and start signals share five conserved nucleotides. The unique structure of the gene overlaps and the presence of a single long IR are conserved among all filoviruses. Here, we sought to determine the impact of the EBOV gene borders during viral transcription. We show that readthrough mRNA synthesis occurs in EBOV-infected cells irrespective of the structure of the gene border, indicating that the gene overlaps do not promote recognition of the gene end signal. However, two consecutive gene end signals at the VP24 gene might improve termination at the VP24-L gene border, ensuring efficient L gene expression. We further demonstrate that the long IR is not essential for but regulates transcription reinitiation in a length-dependent but sequence-independent manner. Mutational analysis of bicistronic minigenomes and recombinant EBOVs showed no direct correlation between IR length and reinitiation rates but demonstrated that specific IR lengths not found naturally in filoviruses profoundly inhibit downstream gene expression. Intriguingly, although truncation of the 144-nucleotide-long IR to 5 nucleotides did not substantially affect EBOV transcription, it led to a significant reduction of viral growth.

IMPORTANCE

Our current understanding of EBOV transcription regulation is limited due to the requirement for high-containment conditions to study this highly pathogenic virus. EBOV is thought to share many mechanistic features with well-analyzed prototype nonsegmented negative-sense RNA viruses. A single polymerase entry site at the 3' end of the genome determines that transcription of the genes is mainly controlled by gene order and cis-acting signals found at the gene borders. Here, we examined the regulatory role of the structurally unique EBOV gene borders during viral transcription. Our data suggest that transcriptional regulation in EBOV is highly complex and differs from that in prototype viruses and further the understanding of this most fundamental process in the filovirus replication cycle. Moreover, our results with recombinant EBOVs suggest a novel role of the long IR found in all filovirus genomes during the viral replication cycle.

Genome sequencing and phylogenetic analysis of 39 human parainfluenza virus type 1 strains isolated from 1997-2010

Thirty-nine human parainfluenza type 1 (HPIV-1) genomes were sequenced from samples collected in Milwaukee, Wisconsin from 1997–2010. Following sequencing, phylogenetic analyses of these sequences plus any publicly available HPIV-1 sequences (from GenBank) were performed. Phylogenetic analysis of the whole genomes, as well as individual genes, revealed that the current HPIV-1 viruses group into three different clades. Previous evolutionary studies of HPIV-1 in Milwaukee revealed that there were two genotypes of HPIV-1 co-circulating in 1991 (previously described as HPIV-1 genotypes C and D). The current study reveals that there are still two different HPIV-1 viruses co-circulating in Milwaukee; however, both groups of HPIV-1 viruses are derived from genotype C indicating that genotype D may no longer be in circulation in Milwaukee. Analyses of genetic diversity indicate that while most of the genome is under purifying selection some regions of the genome are more tolerant of mutation. In the 40 HPIV-1 genomes sequenced in this study, the nucleotide sequence of the L gene is the most conserved while the sequence of the P gene is the most variable. Over the entire protein coding region of the genome, 81 variable amino acid residues were observed and as with nucleotide diversity, the P protein seemed to be the most tolerant of mutation (and contains the greatest proportion of non-synonymous to synonymous substitutions) while the M protein appears to be the least tolerant of amino acid substitution.

Isolation and molecular identification of Sunshine virus, a novel paramyxovirus found in Australian snakes

This paper describes the isolation and molecular identification of a novel paramyxovirus found during an investigation of an outbreak of neurorespiratory disease in a collection of Australian pythons. Using Illumina® high-throughput sequencing, a 17,187 nucleotide sequence was assembled from RNA extracts from infected viper heart cells (VH2) displaying widespread cytopathic effects in the form of multinucleate giant cells. The sequence appears to contain all the coding regions of the genome, including the following predicted paramyxoviral open reading frames (ORFs): 3'--Nucleocapsid (N)--putative Phosphoprotein (P)--Matrix (M)--Fusion (F)--putative attachment protein--Polymerase (L)--5'. There is also a 540 nucleotide ORF between the N and putative P genes that may be an additional coding region. Phylogenetic analyses of the complete N, M, F and L genes support the clustering of this virus within the family Paramyxoviridae but outside both of the current subfamilies: Paramyxovirinae and Pneumovirinae. We propose to name this new virus, Sunshine virus, after the geographic origin of the first isolate--the Sunshine Coast of Queensland, Australia.

Coordinate deletion of N-glycans from the heptad repeats of the fusion F protein of Newcastle disease virus yields a hyperfusogenic virus with increased replication, virulence, and immunogenicity

The role of N-linked glycosylation of the Newcastle disease virus (NDV) fusion (F) protein in viral replication and pathogenesis was examined by eliminating potential acceptor sites using a reverse genetics system for the moderately pathogenic strain Beaudette C (BC). The NDV-BC F protein contains six potential acceptor sites for N-linked glycosylation at residues 85, 191, 366, 447, 471, and 541 (sites Ng1 to Ng6, respectively). The sites at Ng2 and Ng5 are present in heptad repeat (HR) domains HR1 and HR2, respectively, and thus might affect fusion. Each N-glycosylation site was eliminated individually by replacing asparagine (N) with glutamine (Q), and a double mutant (Ng2 + 5) involving the two HR domains was also made. Each mutant was successfully recovered by reverse genetics except for the one involving Ng6, which is present in the cytoplasmic domain. All of the F proteins expressed by the recovered mutant viruses were efficiently cleaved and transported to the infected-cell surface. None of the individual mutations affected viral fusogenicity, but the double mutation at Ng2 and Ng5 in HR1 and HR2 increased fusogenicity >12-fold. The single mutations at sites Ng1, Ng2, and Ng5 resulted in modestly reduced multicycle growth in vitro. These three single mutations were also the most attenuating in eggs and 1-day-old chicks and were associated with decreased replication and spread in 2-week-old chickens. In contrast, the combination of the mutations at Ng2 and Ng5 yielded a virus that, compared to the BC parent, replicated >100-fold more efficiently in vitro, was more virulent in eggs and chicks, replicated more efficiently in chickens with enhanced tropism for the brain and gut, and elicited stronger humoral cell responses. These results illustrate the effects of N-glycosylation of the F protein on NDV pathobiology and suggest that the N-glycans in HR1 and HR2 coordinately downregulate viral fusion and virulence.

Illumination of parainfluenza virus infection and transmission in living animals reveals a tissue-specific dichotomy

The parainfluenza viruses (PIVs) are highly contagious respiratory paramyxoviruses and a leading cause of lower respiratory tract (LRT) disease. Since no vaccines or antivirals exist, non-pharmaceutical interventions are the only means of control for these pathogens. Here we used bioluminescence imaging to visualize the spatial and temporal progression of murine PIV1 (Sendai virus) infection in living mice after intranasal inoculation or exposure by contact. A non-attenuated luciferase reporter virus (rSeV-luc(M-F*)) that expressed high levels of luciferase yet was phenotypically similar to wild-type Sendai virus in vitro and in vivo was generated to allow visualization. After direct intranasal inoculation, we unexpectedly observed that the upper respiratory tract (URT) and trachea supported robust infection under conditions that result in little infection or pathology in the lungs including a low inoculum of virus, an attenuated virus, and strains of mice genetically resistant to lung infection. The high permissivity of the URT and trachea to infection resulted in 100% transmission to naïve contact recipients, even after low-dose (70 PFU) inoculation of genetically resistant BALB/c donor mice. The timing of transmission was consistent with the timing of high viral titers in the URT and trachea of donor animals but was independent of the levels of infection in the lungs of donors. The data therefore reveals a disconnect between transmissibility, which is associated with infection in the URT, and pathogenesis, which arises from infection in the lungs and the immune response. Natural infection after transmission was universally robust in the URT and trachea yet limited in the lungs, inducing protective immunity without weight loss even in genetically susceptible 129/SvJ mice. Overall, these results reveal a dichotomy between PIV infection in the URT and trachea versus the lungs and define a new model for studies of pathogenesis, development of live virus vaccines, and testing of antiviral therapies.

Human parainfluenza viruses (HPIVs) are a leading cause of pediatric hospitalization for lower respiratory tract infection, yet it is unknown why primary infection typically induces immunity without causing severe pathology. To study the determinants of PIV spread within the respiratory tracts of living animals, we developed a model for non-invasive imaging of living mice infected with Sendai virus, the murine counterpart of HPIV1. This system allowed us to measure the temporal and spatial dynamics of paramyxovirus infection throughout the respiratory tracts of living animals after direct inoculation or transmission. We found that the upper respiratory tract and trachea were highly permissive to infection, even under conditions that limit lower respiratory infection and pathogenesis. The timing of transmission coincided with high virus growth in the upper respiratory tracts and trachea of donor mice independent of the extent of infection in the lungs. After transmission, infection spread preferentially in the upper respiratory tract and trachea, inducing protective immunity without weight loss. Our work reveals a disconnect between Sendai virus transmissibility and pathogenicity, and the experimental model developed here will be instrumental in studying PIV pathogenesis.

Residues in the heptad repeat a region of the fusion protein modulate the virulence of Sendai virus in mice

While the molecular basis of fusion (F) protein refolding during membrane fusion has been studied extensively in vitro, little is known about the biological significance of membrane fusion activity in parainfluenza virus replication and pathogenesis in vivo. Two recombinant Sendai viruses, F-L179V and F-K180Q, were generated that contain F protein mutations in the heptad repeat A region of the ectodomain, a region of the protein known to regulate F protein activation. In vitro, the F-L179V virus caused increased syncytium formation (cell-cell membrane fusion) yet had a rate of replication and levels of F protein expression and cleavage similar to wild-type virus. The F-K180Q virus had a reduced replication rate along with reduced levels of F protein expression, cleavage, and fusogenicity. In DBA/2 mice, the hyperfusogenic F-L179V virus induced greater morbidity and mortality than wild-type virus, while the attenuated F-K180Q virus was much less pathogenic. During the first week of infection, virus replication and inflammation in the lungs were similar for wild-type and F-L179V viruses. After approximately 1 week of infection, the clearance of F-L179V virus was delayed, and more extensive interstitial inflammation and necrosis were observed in the lungs, affecting entire lobes of the lungs and having significantly greater numbers of syncytial cell masses in alveolar spaces on day 10. On the other hand, the slower-growing F-K180Q virus caused much less extensive inflammation than wild-type virus, presumably due to its reduced replication rate, and did not cause observable syncytium formation in the lungs. Overall, the results show that residues in the heptad repeat A region of the F protein modulate the virulence of Sendai virus in mice by influencing both the spread and clearance of the virus and the extent and severity of inflammation. An understanding of how the F protein contributes to infection and inflammation in vivo may assist in the development of antiviral therapies against respiratory paramyxoviruses.

Characterization of a recombinant Newcastle disease virus vaccine strain

A recombinant La Sota strain (KBNP-C4152R2L) in which fusion (F) and hemagglutinin-neuraminidase (HN) genes were replaced with those of a contemporary genotype VIId virus, KBNP-4152, has been developed. To attenuate the virulence of the recombinant strain, the F cleavage motif was mutated from (112)RRQKR(116) to (112)GRQAR(116), and to reduce pathogenic instability, a codon which does not allow changes to basic amino acids by single point mutation was inserted at codon 115. In addition a six-nucleotide sequence was inserted into the intergenic region between matrix protein and F genes for attenuation without breaking the "rule-of-six." The HN protein length was increased from 571 to 577 as a marker. Serological tests revealed that the antigenicity of KBNP-C4152R2L was similar to that of KBNP-4152 but distinct from that of the La Sota strain. KBNP-C4152R2L was avirulent (intracerebral pathogenicity index, 0.0; mean death time, >168 h) and stable in pathogenicity through in vivo passages. The killed oil emulsion of and live KBNP-C4152R2L were completely protective against mortality and egg drop caused by virulent strains, and KBNP-C4152R2L was applicable to in ovo vaccination. Therefore, KBNP-C4152R2L is a promising vaccine strain and viral vector in terms of antigenicity, productivity, safety, and pathogenic stability.

Region between the canine distemper virus M and F genes modulates virulence by controlling fusion protein expression

Morbilliviruses, including measles and canine distemper virus (CDV), are nonsegmented, negative-stranded RNA viruses that cause severe diseases in humans and animals. The transcriptional units in their genomes are separated by untranslated regions (UTRs), which contain essential transcription and translation signals. Due to its increased length, the region between the matrix (M) protein and fusion (F) protein open reading frames is of particular interest. In measles virus, the entire F 5' region is untranslated, while several start codons are found in most other morbilliviruses, resulting in a long F protein signal peptide (Fsp). To characterize the role of this region in morbillivirus pathogenesis, we constructed recombinant CDVs, in which either the M-F UTR was replaced with that between the nucleocapsid (N) and phosphoprotein (P) genes, or 106 Fsp residues were deleted. The Fsp deletion alone had no effect in vitro and in vivo. In contrast, substitution of the UTR was associated with a slight increase in F gene and protein expression. Animals infected with this virus either recovered completely or experienced prolonged disease and death due to neuroinvasion. The combination of both changes resulted in a virus with strongly increased F gene and protein expression and complete attenuation. Taken together, our results provide evidence that the region between the morbillivirus M and F genes modulates virulence through transcriptional control of the F gene expression.

Polymerase read-through at the first transcription termination site contributes to regulation of borna disease virus gene expression

An unusually long noncoding sequence is located between the N gene of Borna disease virus (BDV) and the genes for regulatory factor X and polymerase cofactor P. Most of these nucleotides are transcribed and seem to control translation of the bicistronic X/P mRNA. We report here that Vero cells persistently infected with mutant viruses containing minor alterations in this control region showed almost normal levels of N, X, and P proteins but exhibited greatly reduced levels of mRNAs coding for these viral gene products. Surprisingly, cells infected with these BDV mutants accumulated a viral transcript 1.9 kb in length that represents a capped and polyadenylated mRNA containing the coding regions of the N, X, and P genes. Cells infected with wild-type BDV also contained substantial amounts of this read-through mRNA, which yielded both N and P protein when translated in vitro. Viruses carrying mutations that promoted read-through transcription at the first gene junction failed to replicate in the brain of adult rats. In the brains of newborn rats, these mutant viruses were able to replicate after acquiring second-site mutations in or near the termination signal located downstream of the N gene. Thus, sequence elements adjacent to the core termination signal seem to regulate the frequency by which the polymerase terminates transcription after the N gene. We conclude from these observations that BDV uses read-through transcription for fine-tuning the expression of the N, X, and P genes which, in turn, influence viral polymerase activity.

Expression of novel genes encoded by the paramyxovirus J virus

Characterization of the J virus or, in keeping with recent nomenclature recommendations, J paramyxovirus (JPV) genome revealed a unique genome structure, consisting of eight genes in the order 3'-N-P/V/C-M-F-SH-TM-G-L-5'. The small hydrophobic (SH) protein and the transmembrane (TM) protein genes are predicted to encode proteins 69 and 258 aa in size, respectively. The 4401 nt attachment (G) protein gene, much larger than most other paramyxovirus attachment protein genes sequenced to date, encodes a putative 709 aa attachment protein and contains distally a second open reading frame (ORF-X) 2115 nt long. Experiments undertaken in this study were intended to confirm the sequence-based gene allocation of JPV and to determine if proteins encoded by the SH gene, the novel TM gene and ORF-X are expressed. Northern blot analyses carried out on mRNA purified from JPV-infected cells indicated that the putative transcription initiation and termination sequences flanking the SH and TM genes are functional, consistent with their allocation as discrete genes, although a high level of read-through was observed across almost all transcriptional boundaries. Probes specific to the G protein coding region and ORF-X both identified an mRNA species corresponding to the predicted length of the G gene, confirming sequence-based predictions. While the SH and TM proteins were both detected in infected cells, no evidence was found for the expression of ORF-X. Preliminary studies indicate that the novel TM protein is a type II glycosylated integral membrane protein, orientated with its C terminus exposed at the cell surface.

Role for the phosphoprotein P subunit of the paramyxovirus polymerase in limiting induction of host cell antiviral responses

Six amino acid substitutions in the shared N-terminal region of the P subunit of the viral polymerase and the accessory V protein convert the noncytopathic paramyxovirus simian virus 5 (SV5), which is a poor inducer of host cell responses, into a P/V mutant (P/V-CPI-) that induces high levels of apoptosis, interferon-beta (IFN-beta), and proinflammatory cytokines. In this study, we addressed the question of whether these new mutant phenotypes are due to the presence of an altered P protein or of an altered V protein or of both proteins. By the use of the P/V-CPI- mutant as a backbone, new mutant viruses were engineered to express the wild-type (WT) V protein (+V-wt) or WT P protein (+P-wt) from an additional gene inserted between the HN and L genes. In human epithelial cell lines, the +V-wt virus showed reduced activation of apoptosis and lower secretion of IFN-beta and proinflammatory cytokines compared to the parental P/V-CPI- virus. The presence of a V protein lacking the C-terminal cysteine-rich domain (corresponding to the SV5 I protein) did not reduce these host cell responses to P/V-CPI- infection. Unexpectedly, the +P-wt virus, which expressed a WT P subunit of the viral polymerase, also induced much lower levels of host cell responses than the parental P/V-CPI- mutant. For both +V-wt and +P-wt viruses, reduced levels of IFN-beta synthesis correlated with reduced IRF-3 dimerization and nuclear localization of IRF-3 and NF-kappaB, suggesting that the WT P and V proteins acted at an early stage in antiviral pathways. Host cell responses induced by the various P/V mutants directly correlated with levels of viral mRNA accumulation but not with steady-state levels of genomic RNA. Our results support the hypothesis that WT P and V proteins limit induction of antiviral responses by controlling the production of key viral inducers. A model is presented for the mechanism by which both the P subunit of the viral polymerase and the V accessory protein contribute to the ability of a paramyxovirus to limit activation of antiviral responses.

Variations in intergenic region sequences of Human respiratory syncytial virus clinical isolates: analysis of effects on transcriptional regulation

Sequences at the beginnings and ends of Human respiratory syncytial virus (HRSV) genes are necessary for efficient initiation and termination of transcription. The gene start sequences are well conserved and contain signals required for initiation, while the semi-conserved sequences at the gene ends direct transcriptional termination with varying efficiencies. The intergenic regions, which lie between the gene ends and the downstream gene start sequences, are not conserved in length or sequence, and certain positions have been reported to play a role in transcriptional regulation. We have previously shown that the gene end sequences in HRSV subgroup A clinical isolates are variable and that variations found at certain gene ends decreased transcriptional termination and downstream mRNA expression. Here, we have extended this work to examine variation in the intergenic regions between the genes of clinical isolates. We determined the sequences of the eight intergenic regions and the M2/L overlap from clinical isolates from the US and UK and found that all of these regions contained variations from the prototype A2 strain. The amount of variation observed was disparate among the different intergenic regions and did not correlate with length. The effects of selected variant sequences on transcription were examined in the context of subgenomic replicons. While some changes in the intergenic regions had minor effects, certain sequence variations significantly altered transcription termination or initiation. A single nucleotide deletion in the M/SH intergenic region decreased initiation at the SH gene start seven-fold, while changes in the F/M2 intergenic region were found that in some cases increased and in others decreased termination at the F gene end. The P/M intergenic region was the most variable, but none of the changes examined affected either termination at the P gene end or initiation of the downstream M gene start. These results show that in HRSV clinical isolates the intergenic region sequences are variable and that changes in these regions have the potential to affect transcriptional control at the gene junctions.

Variations in transcription termination signals of human respiratory syncytial virus clinical isolates affect gene expression

Human respiratory syncytial virus (HRSV) has a single-stranded, negative-sense RNA genome with 10 genes encoding 11 proteins. Sequences at the beginning of the HRSV genes are highly conserved; however, the gene end sequences vary around a semiconserved consensus sequence, and the nontranscribed intergenic regions vary in both length and sequence. The regions at the junctions between HRSV genes (the gene end sequence of an upstream gene, intergenic region, and the gene start sequence of a downstream gene) contain elements required for efficient termination of the upstream gene and transcription of the downstream gene. Previous studies have examined variation in the HRSV coding sequences, but none have systematically analyzed the noncoding transcriptional control regions for variability. We determined the gene start and gene end sequences of each of the 10 HRSV genes from 14 clinical isolates for variations from the sequence of the prototype A2 strain. No changes were found in any of the gene start sequences. Eight of the 10 gene end sequences, however, contained variations. Several of these, a U(4)-tract instead of a U(6)- or U(5)-tract at the M and SH gene ends, respectively, (U(4)A) and an A-to-G change at position four in the G gene end (A4G), were predicted to affect termination and were examined for their effects on transcription. The changes were found to inhibit transcriptional termination, resulting in increased polycistronic readthrough and correspondingly reduced initiation of the downstream monocistronic mRNA. Viruses with the A4G variant G gene end sequence produced less F protein than those with A2-like G gene end sequences. Examination of additional G gene end sequences available in GenBank revealed that the observed A4G variation was restricted to one phylogenetic lineage of HRSV. All viruses examined within this lineage possessed this variant G gene end sequence. The data presented show that the gene end sequences of naturally occurring HRSV clinical isolates vary from those of the prototypic A2 strain and that certain of these changes inhibit efficient transcriptional termination and downstream gene expression.

Mutations affecting transcriptional termination in the p gene end of subacute sclerosing panencephalitis viruses

Numerous mutations are found in subacute sclerosing panencephalitis (SSPE) viruses, and the M gene is the gene most commonly affected. In some SSPE viruses, such as the MF, Osaka-1, Osaka-2, and Yamagata-1 strains, translation of the M protein is complicated by a transcriptional defect that leads to an almost exclusive synthesis of dicistronic P-M mRNA. To understand the molecular mechanisms of this defect, we sequenced the P gene at the P-M gene junction for several virus strains and probed the involvement of several mutations in the readthrough region via their expression in measles virus minigenomes containing different sequences of the P-M gene junction and flanking reporter genes. The deletion of a single U residue in the U tract of the Osaka-1 strain (3'-UAAUAUUUUU-5') compared with the consensus sequence resulted in a marked reduction of the expression of the downstream reporter gene. In addition, the expression of the downstream gene was markedly decreased by (i) the substitution of a C residue in the U tract of the P gene end of the OSA-2/Fr/B strain of the Osaka-2 virus (3'-UGAUAUUCUU-5' compared with the sequence 3'-UGAUAUUUUU-5' from a sibling virus of the same strain, OSA-2/Fr/V), and (ii) the substitution of a G in the sequence of the P gene end of the Yamagata-1 strain at a variable site immediately upstream from the six-U tract (3'-UGAUGUUUUUU-5' instead of 3'-UGAUUUUUUUU-5'). Mutations at the P gene end can account for the readthrough transcription variation at the P-M gene junction, which directly affects M protein expression.

Naturally occurring substitutions in the P/V gene convert the noncytopathic paramyxovirus simian virus 5 into a virus that induces alpha/beta interferon synthesis and cell death

The V protein of the paramyxovirus simian virus 5 (SV5) is responsible for targeted degradation of STAT1 and the block in alpha/beta interferon (IFN-alpha/beta) signaling that occurs after SV5 infection of human cells. We have analyzed the growth properties of a recombinant SV5 that was engineered to be defective in targeting STAT1 degradation. A recombinant SV5 (rSV5-P/V-CPI-) was engineered to contain six naturally occurring P/V protein mutations, three of which have been shown in previous transfection experiments to disrupt the V-mediated block in IFN-alpha/beta signaling. In contrast to wild-type (WT) SV5, human cells infected with rSV5-P/V-CPI- had STAT1 levels similar to those in mock-infected cells. Cells infected with rSV5-P/V-CPI- were found to express higher-than-WT levels of viral proteins and mRNA, suggesting that the P/V mutations had disrupted the regulation of viral RNA synthesis. Despite the inability to target STAT1 for degradation, single-step growth assays showed that the rSV5-P/V-CPI- mutant virus grew better than WT SV5 in all cell lines tested. Unexpectedly, cells infected with rSV5-P/V-CPI- but not WT SV5 showed an activation of a reporter gene that was under control of the IFN-beta promoter. The secretion of IFN from cells infected with rSV5-P/V-CPI- but not WT SV5 was confirmed by a bioassay for IFN. The rSV5-P/V-CPI- mutant grew to higher titers than did WT rSV5 at early times in multistep growth assays. However, rSV5-P/V-CPI- growth quickly reached a final plateau while WT rSV5 continued to grow and produced a final titer higher than that of rSV5-P/V-CPI- by late times postinfection. In contrast to WT rSV5, infection of a variety of cell lines with rSV5-P/V-CPI- induced cell death pathways with characteristics of apoptosis. Our results confirm a role for the SV5 V protein in blocking IFN signaling but also suggest new roles for the P/V gene products in controlling viral gene expression, the induction of IFN-alpha/beta synthesis, and virus-induced apoptosis.

Transcriptional control of the RNA-dependent RNA polymerase of vesicular stomatitis virus

The nonsegmented negative strand (NNS) RNA viruses include some of the mosr problematic human, animal and plant pathogens extant: for example, rabies virus, Ebola virus, respiratory syncytial virus, the parainfluenza viruses, measles and infectious hemapoietic necrosis virus. The key feature of transcriptional control in the NNS RNA viruses is polymerase entry at a single 3' proximal site followed by obligatory sequential transcription of the linear array of genes. The levels of gene expression are primarily regulated by their position on the genome. The promoter proximal gene is transcribed in greatest abundance and each successive downstream gene is synthesized in progressively lower amounts due to attenuation of transcription at each successive gene junction. In addition, NNS RNA virus gene expression is regulated by cis-acting sequences that reside at the beginning and end of each gene and the intergenic junctions. Using vesicular stomatitis virus (VSV), the prototypic NNS, many of these control elements have been identified.The signals for transcription initiation and 5' end modification and for 3' end polyadenylation and termination have been elucidated. The sequences that determine the ability of the polymerase to slip on the template to generate polyadenylate have been identified and polyadenylation has been shown to be template dependent and integral to the termination process. Transcriptional termination is a key element in control of gene expression of the negative strand RNA viruses and a means by which expression of individual genes may be silenced or regulated within the framework of a single transcriptional promoter. In addition, the fundamental question of the site of entry of the polymerase during transcription has been reexamined and our understanding of the process altered and updated. The ability to engineer changes into infectious viruses has confirmed the action of these elements and as a consequence, it has been shown that transcriptional control is key to controlling the outcome of a viral infection. Finally, the principles of transcriptional regulation have been utilized to develop a new paradigm for systematic attenuation of virulence to develop live attenuated viral vaccines.

Transcriptional termination modulated by nucleotides outside the characterized gene end sequence of respiratory syncytial virus

The genes of respiratory syncytial (RS) virus are transcribed sequentially by the viral RNA polymerase from a single 3'-proximal promoter. Polyadenylation and termination are directed by a sequence at the end of each gene, after which the polymerase crosses an intergenic region and reinitiates at the start sequence of the next gene. The 10 viral genes have different gene end sequences and different termination efficiencies, which allow for regulation of gene expression, since termination of each gene is required for initiation of the downstream gene. RNA sequences within the previously characterized 13 nucleotide gene end, including a conserved sequence 3'-UCAAU-5' and a tract of U residues, are important for termination. In this study, two additional sequence elements outside of the 13 nucleotide gene end were found to modulate termination efficiency: the A residue upstream of the 3'-UCAAU-5' sequence, and the first nucleotide of the intergenic region when it follows a U(4) tract.

The long noncoding region of the human parainfluenza virus type 1 f gene contributes to the read-through transcription at the m-f gene junction

Sendai virus (SV) and human parainfluenza virus type 1 (hPIV1) have genomes consisting of nonsegmented negative-sense RNA in which the six genes are separated by well-conserved intergenic (IG) sequences and transcriptional start (S) and end signals. In hPIV1-infected cells, transcriptional termination at the M-F gene junction is ineffective; a large number of M-F read-through transcripts are produced (T. Bousse, T. Takimoto, K. G. Murti, and A. Portner, Virology 232:44-52, 1997). In contrast, few M-F read-through transcripts are detected in SV-infected cells. Sequence analysis indicated that the hPIV1 IG and S sequences in the M-F junction differ from those of SV. Furthermore, the hPIV1 F gene contains an unusually long noncoding sequence. To identify the cis-acting elements that prevent transcriptional termination at the M-F junction, we rescued recombinant SV (rSVhMFjCG) in which its M-F gene junction was replaced by that of hPIV1. Cells infected with rSVhMFjCG produced an abundance of M-F read-through transcripts; this result indicated that the hPIV1 M-F junction is responsible for inefficient termination. When one or both of the IG and S sites in rSVhMFjCG were replaced by those of SV, the efficiency of transcriptional termination increased but not to the level observed in wild-type SV-infected cells. Deletion of most of the long noncoding region of the hPIV1 F gene in rSVhMFjCG in addition to the mutations in IG and S signals resulted in efficient termination that was equivalent to the level observed in wild-type virus-infected cells. Therefore, the long noncoding sequence of the hPIV1 F gene contains cis-acting element(s) that affects transcriptional termination. Our evaluation of the effect of inefficient transcriptional termination on viral replication in culture revealed that cells infected with rSVhMFjCG produced less F protein than cells infected with wild-type SV and that assembly of the recombinant SV in culture was less efficient. These phenotypes seem to be responsible for the extended survival of mice infected with rSVhMFjCG.

Increased readthrough transcription across the simian virus 5 M-F gene junction leads to growth defects and a global inhibition of viral mRNA synthesis

Recombinant simian virus 5 (rSV5) mutants containing substitutions in the M-F intergenic region were generated to determine the effect of increased readthrough transcription on the paramyxovirus growth cycle. We have previously shown, using an SV5 dicistronic minigenome, that replacement of the 22-base M-F intergenic region with a foreign sequence results in a template (Rep22) that directs very high levels of M-F readthrough transcription. An rSV5 containing the Rep22 substitution grew slower and to final titers that were 50- to 80-fold lower than those of wild-type (WT) rSV5. Cells infected with the Rep22 virus produced very low levels of monocistronic M and F mRNA, consistent with the M-F readthrough phenotype. Surprisingly, Rep22 virus-infected cells also displayed a global decrease in the accumulation of viral mRNA from genes located upstream and downstream of the M-F junction, and overall viral protein synthesis was reduced. Second-site revertants of the Rep22 virus that had regained WT transcription and growth properties contained a single base substitution that increased the M gene end U tract from four to eight residues, suggesting that the growth defects originated from higher-than-normal M-F readthrough transcription. Thus, the primary growth defect for the Rep22 virus appears to be in viral RNA synthesis and not in morphogenesis. A second rSV5 virus (G14), which contained a different foreign M-F intergenic sequence, grew to similar or slightly higher titers than WT rSV5 in some cell types and produced ~1.5- to 2-fold more mRNA and viral protein. The data support the hypothesis that inhibition of Rep22 virus growth is due to increased access by the polymerase to the 5' end of the genome and to the resulting overexpression of L protein. We propose that the elevated naturally occurring M-F readthrough which is characteristic of many paramyxoviruses serves as a mechanism to fine-tune the level of polymerase that is optimal for virus growth.

Analysis of the noncoding regions of measles virus strains in the Edmonston vaccine lineage

The noncoding sequence of five Edmonston vaccine viruses (AIK-C, Moraten, Rubeovax, Schwarz, and Zagreb) and those of a low-passage Edmonston wild-type (wt) measles virus have been determined and compared. Twenty-one nucleotide positions were identified at which Edmonston wt and one or more vaccine strains differed. The location of some of these nucleotide substitutions suggests that they may influence the efficiency of mRNA synthesis, processing, and translation, as well as genome replication and encapsidation. Five nucleotide substitutions were conserved in all of the vaccine strains. Two of these were in the genomic 3'-terminal transcriptional control region and could affect RNA synthesis or encapsidation. Three were found within the 5'-untranslated region of the F mRNA, potentially altering translation control sequences. The remaining vaccine virus base changes were found in one to four vaccine strains. Their genomic localization suggests that some may modify cis-acting regulatory domains, including the Kozak consensus element of the P and M genes, the F gene-end signal, and the F mRNA 5'-untranslated sequence.

Spacing constraints on reinitiation of paramyxovirus transcription: the gene end U tract acts as a spacer to separate gene end from gene start sites

The paramyxovirus gene end U tracts are thought to serve as templates for the addition of a 3' polyA tail to viral mRNAs. The goal of the work described here was to determine the function in transcription of the naturally occurring variability in length of the gene end U tracts of the paramyxovirus simian virus 5 (SV5). An anchored RT-PCR assay was developed to test the hypothesis that the variable U tracts template the addition of variable lengths of polyA tails to mRNAs. The results showed that although the SV5 NP, M, and SH genes encode U tracts of seven, four, and six U residues, respectively, their mRNAs contain similar polyA tails of approximately 250-290 bases. These results indicate that the variable gene end U tracts are functionally equivalent in directing polyadenylation. A reverse genetics system based on a dicistronic minigenome containing the SH-HN gene junction was used to test the hypothesis that the variable U tracks affect the efficiency of transcription termination. Minigenome templates containing an SH gene end with a long U tract of six residues (U6) directed efficient transcription termination and reinitiation at the downstream HN start site with no nucleotide preference for the downstream intergenic region. Surprisingly, truncating the SH gene end U tract to four residues (U4) did not affect SH termination but, rather, reduced downstream HN reinitiation to 20-30% of wild-type levels. Efficient HN reinitiation could be restored to mutant U4 templates in either of two ways: by increasing the U-tract length from four to six residues or by increasing the length of the intergenic region. Efficient HN reinitiation required a minimum of six bases between the last nucleotide in SH and the first nucleotide in HN. We propose that for some paramyxoviruses, the gene end U tract serves a previously unrecognized role as a spacer region between the gene end and gene start sites.

Sendai virus gene start signals are not equivalent in reinitiation capacity: moderation at the fusion protein gene

In paramyxovirus transcription, viral RNA polymerase synthesizes each monocistronic mRNA by recognizing the gene start (S) and end (E) signals flanking each gene. These signal sequences are well conserved in the virus family; nevertheless, they do exhibit some variations even within a virus species. In Sendai virus (SeV) Z strain, the E signals are identical for all six genes but there are four (N, P/M/HN, F, and L) different S signals with one or two nucleotide variations. The significance of these variations for in vitro and in vivo replication has been unknown. We addressed this issue by SeV reverse genetics. The luciferase gene was placed between the N and P gene so that recombinant SeVs expressed luciferase under the control of each of the four different S signals. The S signal for the F gene was found to drive a lower level of transcription than that of the other three, which exhibited comparable reinitiation capacities. The polar attenuation of SeV transcription thus appeared to be not linear but biphasic. Then, a mutant SeV whose F gene S signal was replaced with that used for the P, M, and HN genes was created, and its replication capability was examined. The mutant produced a larger amount of F protein and downstream gene-encoded proteins and replicated faster than wild-type SeV in cultured cells and in embryonated eggs. Compared with the wild type, the mutant virus also replicated faster in mice and was more virulent, requiring a dose 20 times lower to kill 50% of mice. On the other hand, the unique F start sequence as well as the other start sequences are perfectly conserved in all SeV isolates sequenced to date, including highly virulent fresh isolates as well as egg-adapted strains, with a virulence several magnitudes lower than that of the fresh isolates. This moderation of transcription at the F gene may therefore be relevant to viral fitness in nature.

Effect of inserting paramyxovirus simian virus 5 gene junctions at the HN/L gene junction: analysis of accumulation of mRNAs transcribed from rescued viable viruses

Simian parainfluenza virus 5 (SV5) is a prototype of the Paramyxoviridae family of nonsegmented negative-sense RNA viruses. The single-stranded RNA genomes of these viruses contain a series of tandemly linked genes separated by intergenic (IG) sequences flanked by gene-end (GE) and gene-start (GS) sequences. The viral RNA polymerase (vRNAP) complex is thought to enter the genome at its 3' end, and synthesis of mRNAs is thought to occur by a stop-start mechanism in a sequential and polar manner, with transcriptional attenuation occurring primarily at the intergenic regions. As a result, multiple nonoverlapping mRNA species are generated for each single entry of the vRNAP. To investigate the functions of GE, IG, and GS sequences in transcription, we constructed plasmids containing cDNAs of the full-length SV5 genome in which the gene junction sequences (GE, IG, and GS sequences) located between the hemagglutinin-neuraminidase (HN) and the polymerase (L) genes were replaced with the counterpart sequences from other gene junctions. By using reverse genetics, we recovered viable viruses from each cDNA construct, although their growth characteristics varied. Analysis of the HN and L mRNAs by quantitative RNase protection assay indicated that the ratios of HN to L mRNAs varied over a fourfold range. The alteration of the gene junction sequences also permitted examination of the hypothesized requirement for hexamer nucleotide position of the GS sites. The recovery of infectious viruses with transcription initiation sites that occurred at nucleotide positions 1, 2, 3, 5, and 6 of the hexamer suggest that the requirement is nonstringent.

Highly diverse intergenic regions of the paramyxovirus simian virus 5 cooperate with the gene end U tract in viral transcription termination and can influence reinitiation at a downstream gene

A dicistronic minigenome containing the M-F gene junction was used to determine the role of the simian virus 5 (SV5) intergenic regions in transcription. The M-F junction differs from the other SV5 junctions by having a short M gene end U tract of only four residues (U4 tract) and a 22-base M-F intergenic sequence between the M gene end and F gene start site. Replacing the 22-base M-F intergenic region with nonviral sequences resulted in a minigenome template (Rep 22) that was defective in termination at the end of the M gene. Efficient M gene termination could be restored to the mutant Rep 22 template in either of two ways: by increasing the U tract length from four to six residues or by restoring a G residue immediately downstream of the wild-type (WT) U4 tract. In a dicistronic SH-HN minigenome, a U4-G combination was functionally equivalent to the naturally occurring SH U6-A gene end in directing SH transcription termination. In addition to affecting termination, the M-F intergenic region also influenced polymerase reinitiation. In the context of the WT U4-G M gene end, substituting nonviral sequences into the M-F intergenic region had a differential effect on F gene reinitiation, where some but not all nonviral sequences inhibited reinitiation. The inhibition of F gene reinitiation correlated with foreign sequences having a high C content. Deleting 6 bases or inserting 18 additional nucleotides into the middle of the 22-base M-F intergenic segment did not influence M gene termination or F gene reinitiation, indicating that M-F intergenic length per se is not a important factor modulating the SV5 polymerase activity. Our results suggest that the sequence diversity at an SV5 gene junction reflects specific combinations which may differentially affect SV5 gene expression and provide an additional level of transcriptional control beyond that which results from the distance of a gene from the 3' end promoter.