Sequencing and analysis of globally obtained human parainfluenza viruses 1 and 3 genomes

Molecular epidemiology of a Parainfluenza Type 3 virus outbreak: Informing infection control measures on adult hematology wards

OBJECTIVES

Parainfluenza virus type 3 (PIV3) outbreaks among hematology patients are associated with high morbidity and mortality. Prompt implementation of infection prevention (IP) measures has proven to be the most efficacious approach for controlling PIV3 outbreaks within this patient population. The most suitable IP measures can vary depending on the mode of virus transmission, which remains unidentified in most outbreaks. We describe the molecular epidemiology of an outbreak of PIV3 among hematology patients and the development of a new method that allows for the differentiation of outbreak and community strains, from which a closed outbreak could be inferred.

METHODS

Patients were screened for respiratory viruses using multiplex-PCR. PIV3 positive samples with a cycle threshold (Ct)-value of <31 underwent a retrospective characterization via an in-house developed sequence analysis of the hemagglutinin-neuraminidase (HN) gene.

RESULTS

Between July and September 2022, 31 hematology patients were identified with PIV3. Although infection control measures were implemented, the outbreak persisted for nine weeks. Sequencing the HN gene of 27 PIV3 strains from 27 patients revealed that all outbreak strains formed a distinct cluster separate from the control strains, suggestive of a nosocomial transmission route.

CONCLUSIONS

Sequencing the HN gene of PIV3 strains in an outbreak setting enables outbreak strains to be distinguished from community strains. Early molecular characterization of PIV3 strains during an outbreak can serve as a tool in determining potential transmission routes. This, in turn, enables rapid implementation of targeted infection prevention measures, with the goal of minimizing the outbreak's duration and reducing associated morbidity and mortality.

Phylogenetic lineage dynamics of global parainfluenza virus type 3 post-COVID-19 pandemic

During the coronavirus disease 2019 (COVID-19) pandemic, outbreaks of parainfluenza virus type 3 (PIV-3) decreased due to infection control measures. However, a post-pandemic resurgence of PIV-3 has recently been observed. Nonetheless, the role of viral genetic epidemiology, possibly influenced by a genetic bottleneck effect, remains unexplored. We investigated the phylogenetic structure of the publicly available PIV-3 whole-genome and hemagglutinin-neuraminidase (HN) gene sequences spanning the last 65 years, including the COVID-19 pandemic. Sequences were retrieved from the nucleotide database of the National Center for Biotechnology Information using the search term "Human respirovirus 3." Sequence subsets covering all six genes of PIV-3 or the HN gene were designated as the whole-genome and HN surveillance data sets, respectively. Using these data sets, we constructed maximum-likelihood phylogenetic trees and performed a time-scaled analysis using a Bayesian SkyGrid coalescent prior. A total of 455 whole-genome and 1,139 HN gene sequences were extracted, revealing 10 and 11 distinct lineages, respectively, with >98% concurrence in lineage assignments. During the 2020 COVID-19 pandemic, only three single-lineage clusters were identified in Japan, Korea, and the USA. The inferred year of origin for PIV-3 was 1938 (1903-1963) for the whole-genome data set and 1955 (1930-1963) for the HN gene data set. Our study suggests that PIV-3 epidemics in the post-COVID era are likely influenced by a pandemic-driven bottleneck phenomenon and supports previous hypotheses suggesting s that PIV-3 originated during the early half of the 20th century.IMPORTANCEUsing publicly available parainfluenza virus type 3 (PIV-3) whole-genome sequences, we estimated that PIV-3 originated during the 1930s, consistent with previous hypotheses. Lineage typing and time-scaled phylogenetic analysis revealed that PIV-3 experienced a bottleneck phenomenon in Korea and the USA during the coronavirus disease 2019 pandemic. We identified the conservative hemagglutinin-neuraminidase gene as a viable alternative marker in long-term epidemiological studies of PIV-3 when whole-genome analysis is limited.

Molecular evolutionary analyses of the fusion protein gene in human respirovirus 1

Few evolutionary studies of the human respiratory virus (HRV) have been conducted, but most of them have focused on HRV3. In this study, the full-length fusion (F) genes in HRV1 strains collected from various countries were subjected to time-scaled phylogenetic, genome population size, and selective pressure analyses. Antigenicity analysis was performed on the F protein. The time-scaled phylogenetic tree using the Bayesian Markov Chain Monte Carlo method estimated that the common ancestor of the HRV1 F gene diverged in 1957 and eventually formed three lineages. Phylodynamic analyses showed that the genome population size of the F gene has doubled over approximately 80 years. Phylogenetic distances between the strains were short (< 0.02). No positive selection sites were detected for the F protein, whereas many negative selection sites were identified. Almost all conformational epitopes of the F protein, except one in each monomer, did not correspond to the neutralising antibody (NT-Ab) binding sites. These results suggest that the HRV1 F gene has constantly evolved over many years, infecting humans, while the gene may be relatively conserved. Mismatches between computationally predicted epitopes and NT-Ab binding sites may be partially responsible for HRV1 reinfection and other viruses such as HRV3 and respiratory syncytial virus.

Molecular Evolution of Human Parainfluenza Virus Type 2 Based on Hemagglutinin-Neuraminidase Gene

To understand the molecular evolution of human parainfluenza virus type 2 (HPIV2), 21 Hemagglutinin-Neuraminidase (HN) gene sequences covering seven Chinese provinces in 2011 and 2017 to 2021 were combined with 90 published HN sequences worldwide for phylogenetic analysis. The result showed that global HPIV2 could be classified into two distinct clusters (I and II), five lineages (IA to IIE), and four sublineages (IB1 and 2, and IIE1 and 2). The minimum genetic distances between different clusters and lineages were 0.049 and 0.014, respectively. In the last decade, one lineage (IID) and three sublineages (IB1, IB2, and IIE1) have been cocirculating in China, with the sublineages IB2 and IIE1 dominating, while sublineages IB1 and IIE1 are dominant globally. In addition, the spread of HPIV2 had relative spatial clustering, and sublineage IB2 has only been detected in China thus far. The overall evolution rate of HPIV2 was relatively low, on the order of 10-4 substitutions/site/year, except for sublineage IB2 at 10-3 substitutions/site/year. Furthermore, human-animal transmission was observed, suggesting that the HPIV2 might have jumped out of animal reservoirs in approximately 1922, the predicted time of a common ancestor. The entire HN protein was under purifying/negative selection, and the specific amino acid changes and two novel N-glycosylation sites (N316 and N517) in sublineages IB1, IB2, and IIE1 were mostly located in the globular head region of the HN protein. In this study, preliminary evolutionary characteristics of HPIV2 based on the HN gene were obtained, increasing the recognition of the evolution and adaptation of HPIV2. IMPORTANCE The phylogenetic analysis showed that global HPIV2 could be classified into two distinct clusters (I and II) and five lineages (IA to IIE) with at least 0.049 and 0.014 genetic distances between clusters and lineages, respectively. Furthermore, lineages IB and IIE could be further divided into two sublineages (IB1-2 and IIE1-2). All China sequences belong to one lineage and three sublineages (IB1, IB2, IID, and IIE1), among which sublineages IB2 and IIE1 are predominant and cocirculating in China, while sublineages IB1 and IIE1 are dominant globally. The overall evolution rate of HPIV2 is on the order of 10-4 substitutions/site/year, with the highest rate of 2.18 × 10-3 for sublineage IB2. The entire HN protein is under purifying/negative selection, and the specific amino acid substitutions and two novel N-glycosylation sites (N316 and N517) in sublineages IB1, IB2, and IIE1 are mostly located in the globular head region of the HN protein.

Seasonal epidemiological and clinical characteristics of pediatric patients with human parainfluenza virus infection by serotype: a retrospective study

BACKGROUND

The development of the polymerase chain reaction (PCR) test promoted the evaluation of the epidemiological and clinical characteristics of human parainfluenza virus (HPIV) type 4, which has been rarely studied using conventional diagnostic methods. This study aimed to determine the seasonal epidemiological and clinical characteristics of all four HPIV serotypes (HPIV-1, HPIV-2, HPIV-3, and HPIV-4) during the era of PCR testing.

METHODS

The medical records of hospitalized pediatric patients diagnosed with HPIV infections by a multiplex PCR test between 2015 and 2021 were retrospectively reviewed to determine the seasonal distributions of each HPIV serotype. For patients with a single HPIV infection, the clinical characteristics of each HPIV serotype were evaluated and compared with one another.

RESULTS

Among the 514 cases of HPIV infection, HPIV-1, HPIV-2, HPIV-3, and HPIV-4 were identified in 27.2%, 11.9%, 42.6%, and 18.3% of cases, respectively. HPIV-3 was most prevalent in spring, and the other three serotypes were most prevalent in autumn. For patients with a single HPIV infection, those infected by HPIV-1 and HPIV-3 were younger than those infected by HPIV-2 and HPIV-4 (P < 0.001). Croup and lower respiratory tract infection (LRI) were most frequently diagnosed in patients infected by HPIV-1 (P < 0.001) and HPIV-4 (P = 0.002), respectively. During 2020-2021, HPIV-3 was most prevalent in autumn and caused fewer LRIs (P = 0.009) and more seizures (P < 0.001) than during 2015-2019.

CONCLUSIONS

Each HPIV serotype exhibited a distinct seasonal predominance, and some differences in the clinical characteristics of the HPIV serotypes were observed. HPIV-4 acted as an important cause of LRI. Considering the recent changes in the epidemiological and clinical characteristics of HPIV-3, more time-series analyses should be conducted.

Divergent Viruses Discovered in Swine Alter the Understanding of Evolutionary History and Genetic Diversity of the Respirovirus Genus and Related Porcine Parainfluenza Viruses

Paramyxoviridae is a rapidly growing family of viruses, whose potential for cross-species transmission makes it difficult to predict the harm of newly emerging viruses to humans and animals. To better understand their diversity, evolutionary history, and co-evolution with their hosts, we analyzed a collection of porcine parainfluenza virus (PPIV) genomes to reconstruct the species classification basis and evolutionary history of the Respirovirus genus. We sequenced 17 complete genomes of porcine respirovirus 1 (also known as porcine parainfluenza virus 1; PPIV-1), thereby nearly tripling the number of currently available PPIV-1 genomes. We found that PPIV-1 was widely prevalent in China with two divergent lineages, PPIV-1a and PPIV-1b. We further provided evidence that a new species, porcine parainfluenza virus 2 (PPIV-2), had recently emerged in China. Our results pointed to a need for revising the current species demarcation criteria of the Respirovirus genus. In addition, we used PPIV-1 as an example to explore recombination and diversity of the Respirovirus genus. Interestingly, we only detected heterosubtypic recombination events between PPIV-1a and PPIV-1b with no intrasubtypic recombination events. The recombination hotspots highlighted a diverse geography-dependent genome structure of paramyxovirus infecting swine in China. Furthermore, we found no evidence of co-evolution between respirovirus and its host, indicating frequent cross-species transmission. In summary, our analyses showed that swine can be infected with a broad range of respiroviruses and recombination may serve as an important evolutionary mechanism for the Respirovirus genus’ greater diversity in genome structure than previously anticipated.

IMPORTANCE Livestock have emerged as critically underrecognized sources of paramyxovirus diversity, including pigs serving as the source of Nipah virus (NiV) and swine parainfluenza virus type 3, and goats and bovines harboring highly divergent viral lineages. Here, we identified a new species of Respirovirus genus named PPIV-2 in swine and proposed to revise the species demarcation criteria of the Respirovirus genus. We found heterosubtypic recombination events and high genetic diversity in PPIV-1. Further, we showed that genetic recombination may have occurred in the Respirovirus genus which may be associated with host range expansion. The continued expansion of Respirovirus genus diversity in livestock with relatively high human contact rates requires enhanced surveillance and ongoing evaluation of emerging cross-species transmission threats.

Virological Surveillance and Molecular Characterization of Human Parainfluenzavirus Infection in Children with Acute Respiratory Illness: Germany, 2015-2019

Human parainfluenza viruses (HPIVs) are important causes of respiratory illness, especially in young children. However, surveillance for HPIV is rarely performed continuously, and national-level epidemiologic and genetic data are scarce. Within the German sentinel system, to monitor acute respiratory infections (ARI), 4463 respiratory specimens collected from outpatients < 5 years of age between October 2015 and September 2019 were retrospectively screened for HPIV 1-4 using real-time PCR. HPIV was identified in 459 (10%) samples. HPIV-3 was the most common HPIV-type, with 234 detections, followed by HPIV-1 (113), HPIV-4 (61), and HPIV-2 (49). HPIV-3 was more frequently associated with age < 2 years, and HPIV-4 was more frequently associated with pneumonia compared to other HPIV types. HPIV circulation displayed distinct seasonal patterns, which appeared to vary by type. Phylogenetic characterization clustered HPIV-1 in Clades 2 and 3. Reclassification was performed for HPIV-2, provisionally assigning two distinct HPIV-2 groups and six clades, with German HPIV-2s clustering in Clade 2.4. HPIV-3 clustered in C1, C3, C5, and, interestingly, in A. HPIV-4 clustered in Clades 2.1 and 2.2. The results of this study may serve to inform future approaches to diagnose and prevent HPIV infections, which contribute substantially to ARI in young children in Germany.

Epidemiology and seasonality of human parainfluenza serotypes 1-3 in Australian children

Background

Parainfluenza viruses are significant contributors to childhood respiratory illness worldwide, although detailed epidemiological studies are lacking. Few recent Australian studies have investigated serotype‐specific PIV epidemiology, and there is a paucity of southern hemisphere PIV reports. We report age‐stratified PIV hospitalisation rates and a mathematical model of PIV seasonality and dynamics in Western Australia (WA).

Methods

We used linked perinatal, hospital admission and laboratory diagnostic data of 469 589 children born in WA between 1996 and 2012. Age‐specific rates of viral testing and PIV detection in hospitalised children were determined using person time‐at‐risk analysis. PIV seasonality was modelled using a compartmental SEIRS model and complex demodulation methods.

Results

From 2000 to 2012, 9% (n = 43 627) of hospitalised children underwent PIV testing, of which 5% (n = 2218) were positive for PIV‐1, 2 or 3. The highest incidence was in children aged 1‐5 months (PIV‐1:62.6 per 100 000 child‐years, PIV‐2:26.3/100 000, PIV‐3:256/100 000), and hospitalisation rates were three times higher for Aboriginal children compared with non‐Aboriginal children overall (IRR: 2.93). PIV‐1 peaked in the autumn of even‐numbered years, and PIV‐3 annually in the spring, whereas PIV‐2 had inconsistent peak timing. Fitting models to the higher incidence serotypes estimated reproduction numbers of 1.24 (PIV‐1) and 1.72 (PIV‐3).

Conclusion

PIV‐1 and 3 are significant contributors towards infant respiratory hospitalisations. Interventions should prioritise children in the first 6 months of life, with respect to the observed autumn PIV‐1 and spring PIV‐3 activity peaks. Continued surveillance of all serotypes and investigation into PIV‐1 and 3 interventions should be prioritised.

Peripheral blood T cells response in human parainfluenza virus-associated lower respiratory tract infection in children

Human Parainfluenza virus (HPIV) causes lower respiratory tract infections (LRTI) mostly in young children. Respiratory viral infections may decline T cells in circulation and display enhanced pathogenicity. This study is aimed to analyze T cells alterations due to HPIV in children with LRTIs. Children (N = 152) with bronchitis or pneumonia, admitted in tertiary care hospitals were included in the study. Respiratory samples (throat or nasopharyngeal swabs) were taken and HPIV genotypes (1-4) were analyzed through RT-PCR. Peripheral blood T cells, CD3+, CD4+, CD8+, and CD19+, were analyzed in confirmed HPIV positive and healthy control group children through flow cytometry. The positivity rate of HPIV was 24.34% and the most prevalent genotype was HPIV-3 (20.40%). HPIV-1 and HPIV-2 were detected in 0.66% and 02% children respectively. The T lymphocyte counts were observed significantly reduced in children infected with HPIV-3. CD4+ cell (1580 ± 97.87) counts did not change significantly but the lowest CD8+ T cell counts (518.5 ± 74.00) were recorded. Similarly, CD3+ and CD19 cell ratios were also reduced. The CD4/CD8 ratio was significantly higher (3.12 ± 0.59) in the study population as compared to the control group (2.18 ± 0.654). Changes in the count of CD8+ T cells were more pronounced in patients with bronchiolitis and pneumonia. It is concluded that CD8+ T cells show a reduced response to HPIV-3 in children with severe LRTIs suggesting a strong association of these cells with disease severity.

Genetic diversity and evolutionary analysis of human respirovirus type 3 strains isolated in Kenya using complete hemagglutinin-neuraminidase (HN) gene

Human respirovirus type 3 (HRV3) is a leading etiology of lower respiratory tract infections in young children and ranks only second to the human respiratory syncytial virus (HRSV). Despite the public health importance of HRV3, there is limited information about the genetic characteristics and diversity of these viruses in Kenya. To begin to address this gap, we analyzed 35 complete hemagglutinin-neuraminidase (HN) sequences of HRV3 strains isolated in Kenya between 2010 and 2013. Viral RNA was extracted from the isolates, and the entire HN gene amplified by RT-PCR followed by nucleotide sequencing. Phylogenetic analyses of the sequences revealed that all the Kenyan isolates grouped into genetic Cluster C; sub-clusters C1a, C2, and C3a. The majority (54%) of isolates belonged to sub-cluster C3a, followed by C2 (43%) and C1a (2.9%). Sequence analysis revealed high identities between the Kenyan isolates and the HRV3 prototype strain both at the amino acid (96.5-97.9%) and nucleotide (94.3-95.6%) levels. No amino acid variations affecting the catalytic/active sites of the HN glycoprotein were observed among the Kenyan isolates. Selection pressure analyses showed that the HN glycoprotein was evolving under positive selection. Evolutionary analyses revealed that the mean TMRCA for the HN sequence dataset was 1942 (95% HPD: 1928-1957), while the mean evolutionary rate was 4.65x10-4 nucleotide substitutions/site/year (95% HPD: 2.99x10-4 to 6.35x10-4). Overall, our results demonstrate the co-circulation of strains of cluster C HRV3 variants in Kenya during the study period. This is the first study to describe the genetic and molecular evolutionary aspects of HRV3 in Kenya using the complete HN gene.