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

Identification of distinct genotypes in circulating RSV A strains based on variants in the virus replication-associated genes

Respiratory syncytial virus (RSV) is a common cause of respiratory infection that often leads to hospitalization of infected younger children and older adults. RSV is classified into two strains, A and B, each with several subgroups or genotypes. One issue with the definition of these subgroups is the lack of a unified method of identification or genotyping. We propose that genotyping strategies based on the genes coding for replication-associated proteins could provide critical information on the replication capacity of the distinct subgroups, while clearly distinguishing genotypes. Here, we analyzed the virus replication-associated genes N, P, M2, and L from de novo assembled RSV A sequences obtained from 31 newly sequenced samples from hospitalized patients in Philadelphia and 78 additional publicly available sequences from different geographic locations within the United States. In-depth analysis and annotation of variants in the replication-associated proteins identified the polymerase protein L as a robust target for genotyping RSV subgroups. Importantly, our analysis revealed non-synonymous variations in L that were consistently accompanied by conserved changes in its co-factor P or the M2-2 protein, suggesting associations and interactions between specific domains of these proteins. Similar associations were seen among sequences of the related human metapneumovirus. These results highlight L as an alternative to other RSV genotyping targets and demonstrate the value of in-depth analyses and annotations of RSV sequences as it can serve as a foundation for subsequent in vitro and clinical studies on the efficiency of the polymerase and fitness of different virus isolates.IMPORTANCEGiven the historical heterogeneity of respiratory syncytial virus (RSV) and the disease it causes, there is a need to understand the properties of the circulating RSV strains each season. This information would benefit from an informative and consensus method of genotyping the virus. Here, we carried out a variant analysis that shows a pattern of specific variations among the replication-associated genes of RSV A across different seasons. Interestingly, these variation patterns, which were also seen in human metapneumovirus sequences, point to previously defined interactions of domains within these genes, suggesting co-variation in the replication-associated genes. Our results also suggest a genotyping strategy that can prove to be particularly important in understanding the genotype-phenotype correlation in the era of RSV vaccination, where selective pressure on the virus to evolve is anticipated. More importantly, the categorization of pneumoviruses based on these patterns may be of prognostic value.

A Respiratory Syncytial Virus Attachment Gene Variant Associated with More Severe Disease in Infants Decreases Fusion Protein Expression, Which May Facilitate Immune Evasion

This study identified a genotype of respiratory syncytial virus (RSV) associated with increased acute respiratory disease severity in a cohort of previously healthy term infants. The genotype (2stop+A4G) consists of two components. The A4G component is a prevalent point mutation in the 4th position of the gene end transcription termination signal of the G gene of currently circulating RSV strains. The 2stop component is two tandem stop codons at the G gene terminus, preceding the gene end transcription termination signal. To investigate the biological role of these RSV G gene mutations, recombinant RSV strains harboring either a wild-type A2 strain G gene (one stop codon preceding a wild-type gene end signal), an A4G gene end signal preceded by one stop codon, or the 2stop+A4G virulence-associated combination were generated and characterized. Infection with the recombinant A4G (rA4G) RSV mutant resulted in transcriptional readthrough and lower G and fusion (F) protein levels than for the wild type. Addition of a second stop codon preceding the A4G point mutation (2stop+A4G) restored G protein expression but retained lower F protein levels. These data suggest that RSV G and F glycoprotein expression is regulated by transcriptional and translational readthrough. Notably, while rA4G and r2stop+A4G RSV were attenuated in cells and in naive BALB/c mice compared to that for wild-type RSV, the r2stop+A4G RSV was better able to infect BALB/c mice in the presence of preexisting immunity than rA4G RSV. Together, these factors may contribute to the maintenance and virulence of the 2stop+A4G genotype in currently circulating RSV-A strains.IMPORTANCE Strain-specific differences in respiratory syncytial virus (RSV) isolates are associated with differential pathogenesis in mice. However, the role of RSV genotypes in human infection is incompletely understood. This work demonstrates that one such genotype, 2stop+A4G, present in the RSV attachment (G) gene terminus is associated with greater infant disease severity. The genotype consists of two tandem stop codons preceding an A-to-G point mutation in the 4th position of the G gene end transcription termination signal. Virologically, the 2stop+A4G RSV genotype results in reduced levels of the RSV fusion (F) glycoprotein. A recombinant 2stop+A4G RSV was better able to establish infection in the presence of existing RSV immunity than a virus harboring the common A4G mutation. These data suggest that regulation of G and F expression has implications for virulence and, potentially, immune evasion.

Variability analysis and inter-genotype comparison of human respiratory syncytial virus small hydrophobic gene

Background

Small hydrophobic (SH) gene is one of the mostly diverse genomic regions of human respiratory syncytial virus (HRSV). Its coding region constitutes less than 50% of the complete gene length, enabling SH gene to be highly variable and the SH protein highly conserved. In standard HRSV molecular epidemiology studies, solely sequences of the second hypervariable region of the glycoprotein gene (HVR2) are analyzed. To what extent do the strains identical in HVR2 differ elsewhere in genomes is rarely investigated. Our goal was to investigate whether diversity and inter-genotype differences observed for HVR2 are also present in the SH gene.

Methods

We sequenced 198 clinical samples collected within a limited area and time frame. In this HRSV collection, rapid and significant changes in HVR2 occurred.

Results

Over 20% of strains from this pool (containing HRSV genotypes NA1, ON1, GA5, BA9 and BA10) would be incorrectly assumed to be identical to another strain if only the HVR2 region was analysed. The majority of differences found in SH gene were located in the 5′ untranslated region (UTR). Seven indels were detected, one was genotype GA5 specific. An in-frame deletion of 9 nucleotides (coding for amino acids 49–51) was observed in one of group A strains. Fifteen different SH protein sequences were detected; 68% of strains possessed the consensus sequence and most of others differed from the consensus in only one amino acid (only 4 strains differed in 2 amino acids). The majority of differing amino acids in group A viruses had the same identity as the corresponding amino acids in group B strains. When analysis was restricted to strains with identical HVR2 nucleotide sequences and differing SH protein sequences, 75% of differences observed in the SH ectodomain were located within region coding for amino acids 49–51.

Conclusions

Basing HRSV molecular epidemiology studies solely on HVR2 largely underestimates the complexity of circulating virus populations. In strain identification, broadening of the genomic target sequence to SH gene would provide a more comprehensive insight into viral pool versatility and its evolutionary processes.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-1020-9) contains supplementary material, which is available to authorized users.

Recovery of a Paramyxovirus, the Human Metapneumovirus, from Cloned cDNA

Human metapneumovirus (HMPV), a single-stranded negative-sense RNA virus belonging to the family Paramyxoviridae, is associated with respiratory tract illness, primarily in young children and persons with underlying disease. Based on genetic and antigenic variation, HMPV strains are classified into two serotypes, with isolates NL/1/00 and NL/1/99 as prototypes for serotypes A and B, respectively. The development of plasmid-based reverse genetics systems for both serotypes has resulted in developments of a wide range of vaccine candidates against HMPV infection. The approach to virus rescue of HMPV is similar to that used for other paramyxoviruses, starting with mini-replicon assays for optimizations of the rescue protocols and subsequent replacement of the mini genome with a plasmid expressing the cDNA of the full-length viral RNA genome. Here, we provide detailed information on the reverse genetics systems for HMPV.

Direct whole-genome deep-sequencing of human respiratory syncytial virus A and B from Vietnamese children identifies distinct patterns of inter- and intra-host evolution

Human respiratory syncytial virus (RSV) is the major cause of lower respiratory tract infections in children ,2 years of age. Little is known about RSV intra-host genetic diversity over the course of infection or about the immune pressures that drive RSV molecular evolution. We performed whole-genome deep-sequencing on 53 RSV-positive samples (37 RSV subgroup A and 16 RSV subgroup B) collected from the upper airways of hospitalized children in southern Vietnam over two consecutive seasons. RSV A NA1 and RSV B BA9 were the predominant genotypes found in our samples, consistent with other reports on global RSV circulation during the same period. For both RSV A and B, the M gene was the most conserved, confirming its potential as a target for novel therapeutics. The G gene was the most variable and was the only gene under detectable positive selection. Further, positively selected sites inG were found in close proximity to and in some cases overlapped with predicted glycosylation motifs, suggesting that selection on amino acid glycosylation may drive viral genetic diversity. We further identified hotspots and coldspots of intra-host genetic diversity in the RSV genome, some of which may highlight previously unknown regions of functional importance.

Whole genome sequencing and evolutionary analysis of human respiratory syncytial virus A and B from Milwaukee, WI 1998-2010

Background

Respiratory Syncytial Virus (RSV) is the leading cause of lower respiratory-tract infections in infants and young children worldwide. Despite this, only six complete genome sequences of original strains have been previously published, the most recent of which dates back 35 and 26 years for RSV group A and group B respectively.

Methodology/Principal Findings

We present a semi-automated sequencing method allowing for the sequencing of four RSV whole genomes simultaneously. We were able to sequence the complete coding sequences of 13 RSV A and 4 RSV B strains from Milwaukee collected from 1998–2010. Another 12 RSV A and 5 RSV B strains sequenced in this study cover the majority of the genome. All RSV A and RSV B sequences were analyzed by neighbor-joining, maximum parsimony and Bayesian phylogeny methods. Genetic diversity was high among RSV A viruses in Milwaukee including the circulation of multiple genotypes (GA1, GA2, GA5, GA7) with GA2 persisting throughout the 13 years of the study. However, RSV B genomes showed little variation with all belonging to the BA genotype. For RSV A, the same evolutionary patterns and clades were seen consistently across the whole genome including all intergenic, coding, and non-coding regions sequences.

Conclusions/Significance

The sequencing strategy presented in this work allows for RSV A and B genomes to be sequenced simultaneously in two working days and with a low cost. We have significantly increased the amount of genomic data that is available for both RSV A and B, providing the basic molecular characteristics of RSV strains circulating in Milwaukee over the last 13 years. This information can be used for comparative analysis with strains circulating in other communities around the world which should also help with the development of new strategies for control of RSV, specifically vaccine development and improvement of RSV diagnostics.

Whole genome characterization of non-tissue culture adapted HRSV strains in severely infected children

BACKGROUND

Human respiratory syncytial virus (HRSV) is the most important virus causing lower respiratory infection in young children. The complete genetic characterization of RSV clinical strains is a prerequisite for understanding HRSV infection in the clinical context. Current information about the genetic structure of the HRSV genome has largely been obtained using tissue culture adapted viruses. During tissue culture adaptation genetic changes can be introduced into the virus genome, which may obscure subtle variations in the genetic structure of different RSV strains.

METHODS

In this study we describe a novel Sanger sequencing strategy which allowed the complete genetic characterisation of 14 clinical HRSV strains. The viruses were sequenced directly in the nasal washes of severely hospitalized children, and without prior passage of the viruses in tissue culture.

RESULTS

The analysis of nucleotide sequences suggested that vRNA length is a variable factor among primary strains, while the phylogenetic analysis suggests selective pressure for change. The G gene showed the greatest sequence variation (2-6.4%), while small hydrophobic protein and matrix genes were completely conserved across all clinical strains studied. A number of sequence changes in the F, L, M2-1 and M2-2 genes were observed that have not been described in laboratory isolates. The gene junction regions showed more sequence variability, and in particular the intergenic regions showed a highest level of sequence variation. Although the clinical strains grew slower than the HRSVA2 virus isolate in tissue culture, the HRSVA2 isolate and clinical strains formed similar virus structures such as virus filaments and inclusion bodies in infected cells; supporting the clinical relevance of these virus structures.

CONCLUSION

This is the first report to describe the complete genetic characterization of HRSV clinical strains that have been sequenced directly from clinical material. The presence of novel substitutions and deletions in the vRNA of clinical strains emphasize the importance of genomic characterization of non-tissue culture adapted primary strains.

Selection for gene junction sequences important for VSV transcription

The heptauridine tract at each gene end and intergenic region (IGR) at the gene junctions of vesicular stomatitis virus (VSV) have effects on synthesis of the downstream mRNA, independent of their respective roles in termination of the upstream mRNA. To investigate the role of the U tract and the IGR in downstream gene transcription, we altered the N/P gene junction of infectious VSV such that transcription levels would be affected and result in altered molar ratios of the N and P proteins, which are critical for optimal viral RNA replication. The changes included extended IGRs between the N and P genes and shortening the length of the heptauridine tract upstream of the P gene start. Viruses having various combinations of these changes were recovered from cDNA and selective pressure for efficient viral replication was applied by sequential passage in cell culture. The replicative ability and sequence at the altered intergenic junctions were monitored throughout the passages to compare the effects of the changes at the IGR and U tract. VSV variants with wild-type U tracts upstream of the P gene replicated to levels similar to wt VSV. Variants with shortened U tracts were reduced in their ability to replicate. With passage, populations emerged that replicated to higher levels. Sequence analysis revealed that mutations had been selected for in these populations that increased the length of the U tract. This correlated with an increase in abundance of P mRNA and protein to provide improved N:P protein molar ratios. Extended IGRs resulted in decreased downstream transcription but the effect was not as extensive as that caused by shortened U tracts. Extended IGRs were not selected against in 5 passages. Our results indicate that the size of the upstream gene end U tract is an important determinant of efficient downstream gene transcription in infectious virus.

Unravelling the complexities of respiratory syncytial virus RNA synthesis

Human respiratory syncytial virus (RSV) is the leading cause of paediatric respiratory disease and is the focus of antiviral- and vaccine-development programmes. These goals have been aided by an understanding of the virus genome architecture and the mechanisms by which it is expressed and replicated. RSV is a member of the order Mononegavirales and, as such, has a genome consisting of a single strand of negative-sense RNA. At first glance, transcription and genome replication appear straightforward, requiring self-contained promoter regions at the 3' ends of the genome and antigenome RNAs, short cis-acting elements flanking each of the genes and one polymerase. However, from these minimal elements, the virus is able to generate an array of capped, methylated and polyadenylated mRNAs and encapsidated antigenome and genome RNAs, all in the appropriate ratios to facilitate virus replication. The apparent simplicity of genome expression and replication is a consequence of considerable complexity in the polymerase structure and its cognate cis-acting sequences; here, our understanding of mechanisms by which the RSV polymerase proteins interact with signals in the RNA template to produce different RNA products is reviewed.

Respiratory syncytial virus and other pneumoviruses: a review of the international symposium--RSV 2003

The Respiratory Syncytial Virus 2003 symposium took place from 8th-11th November 2003 in Stone Mountain, Georgia, and brought together more than 200 international investigators engaged in RSV research. RSV biology, pathogenesis, and clinical data, as well as RSV vaccines and antivirals, were addressed in the meeting, and this review will aim to briefly summarize and discuss the implications of new findings. The meeting also served as the inauguration of the Robert M. Chanock Award for lifetime achievement in RSV research, an award named in honor of the person who started the field of RSV research by recovering the first human RS virus from infants with severe bronchiolitis in 1956.

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.

Antigenic and genetic variation in human respiratory syncytial virus

BACKGROUND

Human respiratory syncytial virus (HRSV) is a leading cause of serious pediatric respiratory disease worldwide. Natural infection provides only partial protection as repeat infections occur throughout life. A brief review of the extent of antigenic and genetic variation observed in HRSV clinical isolates is presented.

METHODS AND RESULTS

Recent experimental research is reviewed, describing key factors that may explain the ability of HRSV to cause multiple infections in the same individual even in the presence of an existing immune response. It is well-appreciated that variability of the G protein, both between and within antigenic subgroups A and B, is partially responsible for repeat HRSV infections. A high level of nucleotide change resulting in amino acid change provides strong evidence for selective pressure for change in G sequences, thus new HRSV variants. Although little variation in gene-coding sequences is observed in the F protein (the second major protective antigen), new evidence of genetic variation has identified alteration of gene expression levels by selection of changes in the gene end termination signal that precedes the gene encoding the F protein. Due to obligatory sequential transcription, these changes affect downstream gene expression levels. These data suggest that modulation of F protein levels may provide a selective advantage in the presence of a preexisting immune response.

CONCLUSIONS

Experimental data in HRSV demonstrate that variation exists not only in gene-coding sequences but also in the signals that control gene expression. Thus alteration in the expression of key proteins provides a second type of antigenic "variation." A better understanding of these differences is critical to the development of an effective vaccine.