Open access
Editor's Pick
Virology
Research Article
23 May 2024

Indiscriminate activities of different henipavirus polymerase complex proteins allow for efficient minigenome replication in hybrid systems

ABSTRACT

The henipaviruses, including Nipah virus (NiV) and Hendra virus (HeV), are biosafety level 4 (BSL-4) zoonotic pathogens that cause severe neurological and respiratory disease in humans. To study the replication machinery of these viruses, we developed robust minigenome systems that can be safely used in BSL-2 conditions. The nucleocapsid (N), phosphoprotein (P), and large protein (L) of henipaviruses are critical elements of their replication machinery and thus essential support components of the minigenome systems. Here, we tested the effects of diverse combinations of the replication support proteins on the replication capacity of the NiV and HeV minigenomes by exchanging the helper plasmids coding for these proteins among the two viruses. We demonstrate that all combinations including one or more heterologous proteins were capable of replicating both the NiV and HeV minigenomes. Sequence alignment showed identities of 92% for the N protein, 67% for P, and 87% for L. Notably, variations in amino acid residues were not concentrated in the N-P and P-L interacting regions implying that dissimilarities in amino acid composition among NiV and HeV polymerase complex proteins may not impact their interactions. The observed indiscriminate activity of NiV and HeV polymerase complex proteins is different from related viruses, which can support the replication of heterologous genomes only when the whole polymerase complex belongs to the same virus. This newly observed promiscuous property of the henipavirus polymerase complex proteins likely attributed to their conserved interaction regions could potentially be harnessed to develop universal anti-henipavirus antivirals.

IMPORTANCE

Given the severity of disease induced by Hendra and Nipah viruses in humans and the continuous emergence of new henipaviruses as well as henipa-like viruses, it is necessary to conduct a more comprehensive investigation of the biology of henipaviruses and their interaction with the host. The replication of henipaviruses and the development of antiviral agents can be studied in systems that allow experiments to be performed under biosafety level 2 conditions. Here, we developed robust minigenome systems for the Nipah virus (NiV) and Hendra virus (HeV) that provide a convenient alternative for studying NiV and HeV replication. Using these systems, we demonstrate that any combination of the three polymerase complex proteins of NiV and HeV could effectively initiate the replication of both viral minigenomes, which suggests that the interaction regions of the polymerase complex proteins could be effective targets for universal and effective anti-henipavirus interventions.

INTRODUCTION

Henipaviruses, including Nipah virus (NiV) and Hendra virus (HeV), are members of the Henipavirus genus within the Paramyxoviridae family. NiV and HeV represent threatening zoonotic pathogens classified as BSL-4 (biosafety level 4) agents due to their high pathogenicity and lack of available vaccines and antivirals (1, 2). Outbreaks of NiV and HeV have occurred frequently in Southeast Asia and Australia (36). NiV was first isolated in 1999 during an outbreak in pigs that led to subsequent cases of encephalitis among pig farmers in Malaysia and Singapore (7, 8). NiV infection can cause severe respiratory symptoms as well as fatal neurological symptoms, and the virus can spread between humans (9). HeV was identified in Australia in 1994 and is associated with severe respiratory and neurological disease in horses (10). The case fatality rates of NiV and HeV in humans are 60%–100% and there are no efficacious antiviral therapeutics or licensed vaccines for human use (11, 12). To date, a vaccine to protect horses from HeV has been commercialized in Australia. The lack of equivalent prophylactics for human populations remains a critical gap in public health. Furthermore, the emergence of novel Henipaviruses such as Langya (LayV), Gamak (GAKV), and Mojiang (MojV) accentuates the dynamic landscape of this viral family, warranting heightened surveillance and the need for effective intervention strategies (1315).
As BSL-4 pathogens, research with live NiV and HeV needs to be carried out in high containment labs imposing significant limitations to the scientific research and the development of therapeutics against these viruses. The establishment of NiV and HeV minigenome systems that can be used in the BSL-2 conditions is an effective strategy to facilitate broader research aimed at understanding the molecular mechanisms involved in virus replication and serve as a platform for testing antivirals that target these processes (1618).
The genomes of NiV and HeV comprise a single-stranded negative-sense RNA molecule of approximately 18.2 kilobases (kb) in length. These genomes encode a repertoire of structural and non-structural proteins pivotal for virus replication and transcription (19). The nucleocapsid (N), phosphoprotein (P), and large (L) proteins are central to this machinery and collectively orchestrate the assembly of ribonucleoprotein (RNP) complexes essential for viral RNA synthesis (18, 20). In general, minigenome systems for paramyxoviruses, including henipaviruses, consist of a minigenome plasmid in which a reporter gene is flanked by the viral leader and trailer promoter sequences and three helper plasmids each expressing the N, P, and L support proteins under the control of an inducible promoter. After the four plasmids are transfected into cells, the N protein coats the minigenome and minigenome replication is carried out by the viral polymerase L aided by its co-factors N and P.
Minigenome systems have been created for several Mononegavirales. In most cases, the helper plasmids need to be homologous to the minigenome parent virus or all helper plasmids need to be used as a set from the same virus that can then function with a closely related heterologous virus. For example, minigenomes for the rhabdovirus infectious hematopoietic necrosis virus can be replicated efficiently by the complete set of helper plasmids (N, P, and L) from the related hemorrhagic septicemia virus, and vice versa. However, replication is highly inefficient and does not occur when the helper plasmids from these viruses are mixed (21). Similarly, sets of heterologous proteins worked in trans to replicate minigenomes of closely related strains of vesicular stomatitis virus but replication did not occur when the helper plasmids came from mixed strains (22). Cross-activity of polymerase and its co-factors has been also reported for the pneumoviruses human, bovine, and ovine respiratory syncytial virus (RSV) (23); however, filovirus helper plasmids do not seem to work in trans even if present as a set (24). Among paramyxoviruses, it has been shown that a full-length infectious clone of Sendai virus (SeV) could be successfully rescued after co-transfection with the helper plasmid set from human parainfluenza virus 1 (HPIV1) and human parainfluenza virus 3 (HPIV3) strains, but mixing helper plasmids or using helper plasmids from a more distant morbillivirus or pneumovirus was ineffective (2527). Also, replication of heterologous minigenomes of three morbilliviruses only happens when full sets of helper plasmids are used (28).
Here, we have developed an efficient minigenome system for studying NiV and HeV replication under BSL-2 conditions as an alternative system to study virus replication. Using these systems, we found unexpected remarkable promiscuity among henipavirus polymerase complex proteins that allows efficient replication of the NiV and HeV minigenomes in hybrid systems without the need for homologous components within these viruses.

RESULTS

Construction of NiV and HeV minigenome systems in BSR-T7/5 cells

The henipavirus viral genomes consist of a leader sequence, six protein-encoding genes, and a trailer sequence (Fig. 1A). We developed minigenome systems to utilize the T7 polymerase that is constitutively expressed in BSR-T7/5 cells (29). The minigenome plasmid (MG) consists of two separate functional units (Fig. 1B). The control unit encodes an internal ribosome entry site (IRES) and an enhanced green fluorescent protein (eGFP) reporter gene under the T7 promoter that acts as a control for the minigenome transfection and the T7 polymerase activity. The viral replicon unit includes the virus 5′ end of the L gene fragment, a reporter mCherry gene flanked by the leader and trailer sequences, and the non-coding region (NCR) of N and L genes, which are cis elements essential for viral replication and transcription. In addition, the N gene 3′ NCR sequence was added between the L gene fragment and the mCherry gene to ensure that mCherry can be transcribed. The reverse complementary sequences of each viral fragment and the mCherry gene were used for the construction of the viral replicon unit. The viral replicon unit is flanked by the self-cleaving hammerhead ribozyme (Hh-Rbz) before the trailer sequence and the hepatitis delta virus ribozyme (HDV-Rbz) (30) after the leader sequence to ensure transcription products have a nucleotide length divisible by six, which is necessary for efficient replication as described by the “rule of six.”
Fig 1
Fig 1 Henipavirus minigenome strategy. (A) Schematic representation of a generic Henipavirus genome. The six structural genes are represented by orange boxes. 3′ NCR and 5′ NCR of each gene are shown by white boxes. (B) Schematic of NiV and HeV minigenome system plasmids. RNA synthesis and translation in the viral replicon unit of the minigenome system are depicted. Hh-Rbz, hammerhead ribozyme. HDV-Rbz, hepatitis delta virus ribozyme. (C) Schematic of the three support/helper plasmids.
After transfecting the MG and helper plasmids into BSR-T7/5 cells, both the viral replicon unit and the control unit are transcribed by the T7 RNA polymerase produced by the BSR-T7/5 cells. The eGFP RNA is then translated by host translation machinery in an IRES-dependent translation manner. For the viral replicon unit, the minigenomic RNA generated by the T7 RNA polymerase is coated by N proteins expressed from a helper plasmid (Fig. 1C) to form a ribonucleoprotein and act as replication template (Fig. 1B). The L protein, expressed from a separate helper plasmid (Fig. 1C), will then recognize the promoter located in the leader sequence of the minigenomic RNA and produce antiminigenomic RNA. Antiminigenomic RNA can also act as a template for minigenomic RNA synthesis. Viral polymerase complex proteins also recognize the gene start (GS) signal in the N gene 5′ NCR and gene end (GE) signal in the N gene 3′ NCR of minigenomic RNA and initiate mCherry mRNA transcription. The L protein then modifies mCherry mRNAs to add the 5′ cap and 3′ polyA tail, and finally host translation machinery produces the mCherry protein. Consequently, the green fluorescence signal of eGFP and the red fluorescence signal of mCherry protein are readouts of this system.
To initially test and validate the minigenomes, all four plasmids (MG, N, P, and L) were transfected into BSR-T7/5 cells (MG: 875 ng; N: 312 ng; P: 200 ng; L: 100 ng). The microscopic fluorescence signal was captured at indicated time points and analyzed using ImageJ. mCherry signal was observed in about 9.1% of cells for NiV and 2.5% for HeV transfection group at 72 h post-transfection (hpt) (Fig. 2A and B). No mCherry signal was observed in the control group which was transfected with the MG, N, and P plasmids but not with the L plasmid. These data demonstrate that the dual reporter genes minigenome systems for NiV and HeV are functional in BSR-T7/5 cells.
Fig 2
Fig 2 Functional validation of henipavirus minigenome systems. (A) Verification of the NiV Bangladesh strain minigenome system. BSR-T7/5 cells were transfected with NiV MG and three helper plasmids (MG: 875 ng; N: 312 ng; P: 200 ng; L: 100 ng). The control group was transfected with NiV MG, pTM-N, pTM-P of NiV, and pTM-1 vector (MG: 875 ng; N: 312 ng; P: 200 ng; pTM1: 100 ng). (B) Verification of the HeV Redlands strain minigenome system. BSR-T7/5 cells were transfected with HeV MG and three helper plasmids (MG: 875 ng; N: 312 ng; P: 200 ng; L: 100 ng). The control group was transfected with HeV MG, pTM-N, pTM-P of HeV, and pTM-1 vector (MG: 875 ng; N: 312 ng; P: 200 ng; pTM1: 100 ng). Expression of eGFP (Green) and mCherry (Red) was observed at 72 hpt by widefield microscopy at 5× magnification. Digital zoomed images are shown in the panel on the right. Cells with both green and red fluorescence signals are indicated by white arrows. The fluorescent images shown are representative of three independent experiments. Scale bar lengths are indicated.

Optimization of transfection efficiency of henipavirus minigenome systems

To optimize the transfection efficiency of both minigenome systems, we tested four different ratios of MG: N: P: L, including ratios that have been previously reported for henipaviruses minigenome systems or virus rescue (3135) (Table 1). “Ratio 2” (MG: 500 ng; N: 150 ng; P: 50 ng; L: 60 ng) resulted in drastically enhanced efficiency at 72 hpt, with 35.6% and 11.0% of cells positive for reporter gene expression for the NiV and HeV minigenomes, respectively, compared to only 2.4% and 6.4% for NiV and HeV, respectively, in “Ratio 1,” the next best tested (Fig. 3A).
Fig 3
Fig 3 Optimization of henipavirus minigenome systems. (A) BSR-T7/5 cells were transfected with NiV MG alongside the three helper plasmids of NiV (left panel) or HeV MG and the three helper plasmids of HeV (right panel) at four different ratios shown in Table 1. Expression of eGFP (Green) and mCherry (red) was observed at 72 hpt. The fluorescent images shown are representative of three independent experiments. Transfection efficiency is shown as mean ± SD below the images. (B) Timelapse microscopy images of BSR-T7/5 cells transfected with MG and three helper plasmids of NiV (MG: 500 ng; N: 150 ng; P: 50 ng; L: 60 ng) at 30–47 hpt and eGFP (green) and mCherry (Red) signal are shown at the indicated timepoints. Scale bar lengths are indicated.
TABLE 1
TABLE 1 Transfection ratios of NiV and HeV minigenome system plasmids
 MG/ngpTM-N/ngpTM-P/ngpTM-L/ng
Ratio 1875312200100
Ratio 25001505060
Ratio 3600600400200
Ratio 41000410260410
Although the “Ratio 2” was identified to have the highest transfection efficiency, we noticed that BSR-T7/5 cells showed several different fluorescence signals among the transfected cells. As expected, most cells showed either only eGFP fluorescence (green) or dual eGFP and mCherry fluorescence (orange), but a few cells showed mCherry only (red), which we did not expect as all mCherry signal should theoretically be accompanied by eGFP in our system (Fig. 1A). We postulated that we were only seeing a snapshot of the reporter protein dynamics and missing the eGFP signal in some cells when the pictures were taken. To test this hypothesis, we used live-cell imaging of cells transfected with the four NiV minigenome system plasmids to investigate the fluorescence signals through time in single cells. Several cells showed a bright green fluorescence signal at 33 hpt and then started showing a yellow fluorescence signal at 36 hpt gradually increasing the expression level of mCherry (Fig. 3B, View 1). In other cases, cells showed a single red fluorescence signal at the beginning and turned to an orange/yellow at 44 hpt when eGFP was expressed (Fig. 3B, View 2). A different subset of cells showed orange/yellow signal during most of the live-cell imaging time course because of similar expression levels of mCherry and eGFP (Fig. 3B, View 3). We did not see any cell that maintained only mCherry signal throughout the entire time course. These results proved that the fluorescence of single cells was dynamic over time, and we concluded that our minigenome system is working as we expected in its optimized conditions.

Cross-activity of NiV and HeV polymerase complex proteins

The polymerase complex proteins N, P, and L are essential components for paramyxovirus transcription and replication. As previously discussed, in most cases, homologous support proteins or the full set of heterologous support proteins can replicate a closely related virus. During the optimization of our minigenome systems for NiV and HeV, we assessed the impact of having heterologous components of the polymerase complex on the minigenome efficiency. As previously reported (20), homologous sets of NiV and HeV polymerase components work well in trans to replicate and transcribe the heterologous genome (Fig. 4A and B, Com. 7). However, to our surprise, efficient trans polymerase activity could be seen with all combinations of different polymerase components. As shown by mCherry expression, the NiV MG successfully replicates when HeV N, P, and/or L helper plasmids are used in any combination with the NiV helper plasmids (Fig. 4A). Similarly, all combinations of heterologous NiV helper plasmids work to replicate the HeV MG (Fig. 4B). To test whether the observed promiscuity of the henipaviruses polymerase activities was limited to closely related viruses, we used helper plasmids from the paramyxovirus SeV with the NiV minigenome. Transfection of the NiV MG with SeV helper plasmids did not result in mCherry signal, nor could NiV support proteins initiate replication of the SeV minigenome (Fig. 4C). These data suggest that polymerase-associated proteins of NiV and HeV work well with each other but not with other paramyxoviruses.
Fig 4
Fig 4 Cross-activity of NiV, HeV, and SeV polymerase complex proteins. (A) Cross-activity of different combinations NiV and HeV support proteins in the NiV minigenome system. NiV MG (500 ng) was transfected into BSR-T7/5 cells with the indicated combinations of NiV and HeV helper plasmids (N: 150 ng; P: 50 ng; L: 60 ng). eGFP (Green) and mCherry (red) were observed at 72 hpt. (B) Cross-activity of different sets NiV and HeV support proteins in the HeV minigenome system. HeV MG (500 ng) was transfected into BSR-T7/5 cells with different combinations of NiV and HeV helper plasmids (N: 150 ng; P: 50 ng; L: 60 ng). eGFP (Green) and mCherry (red) were observed at 72 hpt. (C) Cross-activity verification of NiV and SeV minigenome systems. NiV or SeV MG (500 ng) plasmids were transfected into BSR-T7/5 cells with different combinations of NiV and HeV helper plasmids (N: 150 ng; P: 50 ng; L: 60 ng). eGFP (Green) and mCherry (red) were observed at 48 hpt. The fluorescent images shown are representative of three independent experiments. Com., Combination. Scale bar shows 200 µm.

Conserved domains for protein-protein interactions likely allow efficient replication, irrespective of support protein combinations

To understand why different combinations of helper plasmids did not affect the replication of Henipavirus minigenomes, as is seen for other Mononegavirales, we compared the amino acid sequences of the three polymerase complex proteins of the NiV Bangladesh and HeV Redlands strains used as bases for our minigenomes. We first focused on the N protein since the C-terminal intrinsically disordered domain (Ntail) of the N protein has four defined functional boxes including Box 3, which binds to the C-terminal X domain of viral phosphoprotein (PXD) to tether P onto the nucleocapsid template (Fig. 5A) (36, 37). Sequence alignment showed that NiV and HeV N proteins share 91.7% identity (defined as the percentage of the same amino acids) and 96.8% similarity (defined as the percentage of the same amino acids plus conservatively replaced amino acids) (Fig. S1A; Table 2). Only two amino acids (A488I and A492T) within Box 3 of NiV and HeV showed non-conservative replacement while Box 2 amino acid sequences of SeV N protein differed greatly from Box 3 amino acid sequences of NiV and HeV (Fig. 5B). We also compared all amino acid sequences available in the NCBI virus database for NiV and HeV Box 3 and found that most of the Box 3 amino acid sequences of NiV and HeV were conserved (Fig. S2). Differences were mainly shown in five novel HeV-g2 variants (UCY33663, UCY33672, UCY33681, UCY33690, and QYC64598) that were isolated between 2013 and 2020 in Australia (5, 38).
Fig 5
Fig 5 Polymerase complex proteins sequence analysis of NiV Bangladesh and HeV Redlands strains. (A) Schematic representation of the interaction between the N and P proteins via key domains. Ncore, structured N-terminal domain. Ntail, disordered C-terminal domain. PMD, P multimerization domain; PNT, N-terminal region of P; PCT, C-terminal region of P. PXD, X domain of P. CR, conserved region. (B) Sequence alignment of the Ntail of NiV Bangladesh strain, HeV Redlands strain, and SeV Cantell strain. The first and last amino acid positions for each region based on the NiV Bangladesh strain sequence are displayed above or below. (C) Sequence alignment of NiV Bangladesh strain, HeV Redlands strain, and SeV Cantell strain P protein PMD and PXD. PMD and PXD are indicated by orange boxes. The first and last amino acid positions for PMD and PXD are displayed above or below based on the NiV Bangladesh strain sequence. (D) Sequence differences analysis of each region of Henipavirus L protein. Schematic representation of Henipavirus L protein, including six conserved regions (CR I-CR VI), which are presented in magenta boxes. Similarities of each region are indicated below the boxes. Sequence alignment of the key CR I of NiV Bangladesh strain, HeV Redlands strain, and SeV Cantell strain and the least conserved region (607-711) between CR II and CR III of NiV Bangladesh strain and HeV Redlands strain are shown individually. The amino acid positions are shown based on the NiV Bangladesh strain sequence. Red text indicates non-similar residues. Blue text represents similar but non-identical residues. The dots represent identical amino acids.
TABLE 2
TABLE 2 Identities and similarities of NiV and HeV polymerase complex proteins
 N proteinP proteinL protein
Identity91.7%66.6%87.1%
Similarity96.8%77.1%94.3%
The P proteins shared 66.6% identity and 77.1% similarity (Table 2). Although the identity of the P protein is lower, sequence alignment results showed that most differences were distributed in the N-terminal region of P (PNT) (Fig. S1B). The PXD that binds to Ntail was conserved with only two significant amino acid differences (G683D and N690T) (Fig. 5C). Similarly, the P multimerization domain (PMD) that interacts with the conserved region I(CRI) of the L protein to recruit L onto the nucleocapsid template (39) was conserved between the two henipaviruses. PXD and PMD of SeV displayed less identical amino acids (Fig. 5C). Alignment of the PMD and PXD sequences of all NiV and HeV strains available from NCBI indicated those two domains were highly conserved, except for one non-similar amino acid substitution in the PMD (G504S) and in PXD (G683D) (Fig. S3 and S4).
Finally, we examined NiV and HeV L protein, a multifunctional enzyme that is conserved among Mononegavirales. L has six conserved regions (CRs) named CR I-CR VI (40, 41). CR I, II, and III are in the RNA-dependent RNA polymerization (RdRp) domain, CR IV and V are in the cap addition (Cap) domain and CR VI is in the cap methylation (MT) domain (42). The L proteins of NiV Bangladesh and HeV Redlands shared 87.1% identity and 94.3% similarity (Table 2). The amino acid difference distribution analysis revealed a similarity of over 94.9% for all six CRs. The primary distinction was observed in the region separating CR II to CR III (Fig. 5D). No significant function has been identified for this area of the protein. We also found that the CR I region of L, which binds primarily to the PMD of the P protein to facilitate P-L interactions, was conserved in all available sequences of henipaviruses, except for four of the novel HeV-g2 strains (Fig. S5). CR I of SeV L protein was less conserved than NiV and HeV (Fig. 5D).
These data suggest that dissimilarities within the polymerase complex proteins of NiV and HeV are not concentrated in regions crucial for N, P, and L interaction, which may account for the ability of helper plasmids with diverse arrangement combinations to facilitate replication of both NiV and HeV. Furthermore, sequence analysis revealed the conservation of these crucial interacting regions across a wide range of NiV and HeV isolates, supporting potential cross-activity between viral proteins. The differences in these key interacting domains between SeV and henipavirus replication proteins may explain the failure of replication in henipa-SeV hybrid systems.

DISCUSSION

Henipaviruses are highly pathogenic and global cases are on the rise (5, 43, 44), making the lack of effective therapeutics a pressing concern for human health. The requirement for BSL-4 containment poses significant challenges to the study of these viruses. Minigenome systems are powerful tools to safely circumvent the need for BSL4 conditions and have been widely used in virus research (4547), especially for highly pathogenetic agents including Ebola virus, Zika virus, Marburg virus (MARV), and henipaviruses (16, 18, 48, 49). Here, we describe the establishment of a T7 RNA polymerase-based minigenome system for both NiV and HeV and the use of these platforms for the analysis of viral polymerase complex proteins cross-activity.
Previously reported minigenome systems for henipaviruses have used chloramphenicol acetyltransferase (CAT), luciferase, or fluorescent proteins such as RFP as the reporter genes to replace all viral structural genes (20, 50, 51). We constructed a bi-cistronic minigenome system including two separate units. The viral replicon unit contained 2 kb of the L gene and was originally designed for further research of henipavirus copy-back viral genome production during the replication, and the downstream mCherry gene. The control unit expresses eGFP if the transfection and the T7 polymerase are working well. Unexpectedly, we found a single red fluorescent signal in some cells, which has been demonstrated in another bi-cistronic minigenome system for NiV signal based on the relative expression level of eGFP and mCherry protein (50). Using live imaging, we validated the dynamic fluorescent signal based on the relative expression level of eGFP and mCherry protein. We confirmed that all cells eventually express both reporters, although not all cells expressed both reporters at the same time after transfection (Fig. 3B). After optimizing the ratio of four plasmids, we established a robust minigenome system for both NiV and HeV that can be used for further study of henipaviruses. A similar dual reporter strategy could be applied to the research of pro-viral and anti-viral replication factors of henipaviruses and other viruses.
Replication and transcription of paramyxoviruses require homotypic support proteins, including N, P, and L, although in some conditions heterotypic sets are functional among closely related viruses (20, 25, 28, 52). Interestingly, all previous reports on paramyxoviruses show that only support proteins from the same virus can initiate effective replication of a heterologous viral genome or minigenome, highlighting the importance of interactions between N-P and P-L in viral replication (5355). The viral genome of paramyxoviruses is a negative-sense RNA, which is coated by the N protein to form a ribonucleoprotein (RNP). Functional regions in the Ntail bind to PXD to recruit the P protein to the RNP. At the same time, the PMD of the P protein interacts with CR I of the L protein, attaching L to the template. The precise interaction between the polymerase complex proteins is critical for the successful replication of the viral genome and for this reason, the polymerase complex is a key therapeutic target for paramyxoviruses (52, 53, 56). Consistent with previous reports, we show in our system that the three support proteins of HeV initiate the replication of NiV (20). We also show that N, P, and L of NiV exhibited the ability to replicate HeV minigenome. Strikingly, we also found that various combinations of support proteins enabled replication for both NiV and HeV minigenomes, even when the support proteins were not originating from the same virus (Fig. 4A and B). Interestingly, this phenotype that mixed support proteins from different viruses could initiate the viral replication effectively has not been identified in other paramyxoviruses.
The L protein is conserved among paramyxoviruses, while N and P proteins vary. Four boxes exist in Ntail of henipavirus N protein and Box 3 binds to PXD. This contrasts with both measles virus (MeV) and SeV Ntail, which have only three boxes, and Box 2 interacts with PXD (5759). The length of the P gene varies greatly and is less conserved. The P gene of two henipaviruses, NiV and HeV, showed a lower identity than N or L (Table 2). Of note, the glycoproteins of NiV and HeV demonstrated robust cross-functional compatibility, displaying efficient heterotypic activity with each other. Conversely, no such heterotypic activity was detected when comparing them with the envelope glycoproteins of the morbilliviruses MeV and Canine distemper virus (CDV) (60). The amino acid differences in proteins of the replication complex partly explain the inability of the NiV helper plasmids to initiate SeV minigenome replication. The ability of diverse combinations of helper plasmids to facilitate replication of both NiV and HeV minigenomes may be attributed, at least in part, to the conservation of critical protein interaction regions (Box 3 of N, PMD and PXD of P, and CR I of L). Comprehensive sequence alignment suggests potential cross-activity between proteins across a wide range of henipavirus isolates (Fig. S2 to S5). The increased amino acid variations in emerging HeV-g2 strains underscore the imperative for timely detection and analysis of novel strains. Moreover, these findings suggest that cross-interaction patterns may evolve alongside the emergence of new strains. While henipavirus outbreaks are currently restricted to Southeast Asia and Australia, the emergence of more henipaviruses and henipa-like viruses raises a serious public health concern of a global pandemic (61, 62). Our findings demonstrate the cross-activity between NiV and HeV, suggesting the possibility of recombinant variants. It is possible that cross-activity of polymerase complex proteins may not only exist between NiV and HeV, but also happen among different henipaviruses and henipa-like viruses. While this cross-activity may facilitate virus evolution during coinfection of two or more viruses, representing a threat to public health, it also suggests a therapeutic potential. Subunit vaccines and peptides targeting disruptions in key regions of polymerase complex proteins analog to henipaviruses N Box 3, P PXD or PMD, or L CR I domains could competitively bind to viral N, P, or L proteins. This approach could lead to the development of a novel broad-spectrum antiviral target against henipaviruses that have the indiscriminate activities of polymerase complex proteins. In addition, the identification of critical amino acids for polymerase complex protein interaction could facilitate the development of small-molecule drugs to block interactions between N, P, and L proteins. These analogs and small-molecule drugs could be used universally against paramyxoviruses. Furthermore, recombinant viruses produced through co-infection with other viruses, harboring fragments from other viral N, P, or L proteins, may also be susceptible to the effects of analogs or small-molecule drugs. Nevertheless, a notable concern arises regarding the emergence of novel variants, such as HeV-g2 strains, which may exhibit enhanced divergence within key interaction regions, potentially compromising the efficacy of such therapeutic targeting strategies.
By utilizing two effective minigenome systems of NiV and HeV, we discovered that different henipavirus polymerase complex proteins have indiscriminate activities and can facilitate heterologous replication of NiV and HeV minigenomes. These data pave the road for future studies on henipaviruses and shed light on our understanding of cross-activity between paramyxovirus polymerase complex proteins, raising novel considerations for viral surveillance and therapeutic development.

MATERIALS AND METHODS

Cell culture

BSR-T7/5 cells (hamster kidney cells expressing bacteriophage T7 RNA polymerase, kindly provided by K. Conzelmann (29)) were maintained in Dulbecco’s modified Eagle medium (DMEM) (ThermoFisher) supplemented with 10% fetal bovine serum (FBS), L-glutamine 2 mM (Invitrogen), gentamicin 50 ng/mL (ThermoFisher), sodium pyruvate 1 mM (Invitrogen), and 400 μg G418 Sulfate (Invitrogen) at 37°C with 5% CO2. Cells were treated with mycoplasma removal agent (MP Biomedicals) before use and screened monthly for mycoplasma contamination with MycoAlert Plus mycoplasma testing kit (Lonza).

Construction of minigenome system plasmids

To generate the helper plasmids, HeV N, P, L gene sequences were synthesized by IDT (Integrated DNA Technologies) based on the sequence of HeV Redlands strain (GenBank accession No. HM044317) (63) . NiV N and L gene sequences were were synthesized by IDT based on the NiV Bangladesh strain sequence (GenBank accession No. AY988601). The NiV P gene was based on the NiV Bangladesh genome and internally codon-optimized by IDT to reduce its complexity before synthesis. N, P, and L gene open reading frames (ORF) were amplified with PrimeSTAR Max DNA Polymerase (TAKARA) and ligated between the NcoI and BamHI sites of the pTM1 vector with In-Fusion Snap Assembly Master Mix (TAKARA) according to the manufacturer’s protocols. SeV Cantell (SeV-C) strain (GenBank accession No. OR764764) (64) helper plasmids were also constructed by cloning and inserting N, P, and L gene ORFs into pTM1 vector. All plasmids were confirmed by full-length sequencing. For the generation of the dual reporter minigenome plasmids, the two separate functional units were constructed separately. Five nucleotides were inserted before the mCherry start codon to ensure the sequence between the T7 promoter and T7 terminator of the viral replicon unit conformed the “rule of six.” The minigenome plasmids were constructed by inserting these two units’ sequences into the pSL1180 vector. SeV-C minigenome plasmid was constructed in the same way but containing all the viral elements of the SeV-Cantell strain.

Plasmid transfection

BSR-T7/5 cells were seeded in a 24-well plate the day before transfection. Cells were washed twice with PBS (Invitrogen) and incubated in 500 mL Opti-MEM (ThermoFisher) before transfection. BSR-T7/5 cells were transfected with MG, pTM-N, pTM-P, and 1 pTM-L using TransIT-LT1 (Mirus) and incubated at 37°C. The negative control group was transfected with MG, pTM-N, pTM-P, and pTM1 vector plasmid. Plasmids were transfected into BSR-T7/5 cells according to the ratios shown in Table 1. Fluorescence signal was observed daily after transfection. Images were captured with the 5× and 20× objectives of a Zeiss Axio observer Widefield microscope. The transfection efficiency of the viral replicon unit was calculated as mCherry+cells/eGFP+ cells.

Time-lapse microscopy

BST-T7/5 cells were observed 30–47 h after transfection with NiV MG and the three helper plasmids of NiV at the “Ratio 2” (Table 1). During observation, cells were maintained in Opti-MEM media at 37°C. eGFP and mCherry fluorescence were visualized using the 5× objectives of a Zeiss Axio observer Widefield fluorescence microscope every 15 min.

Sequence analysis

The complete genomes of NiV Bangladesh strain and HeV Redlands strain were obtained from GenBank of the National Center for Biotechnology Information (NCBI). Amino acid identity and similarity of the N, P, and L protein amino acids of the Bangladesh strain and Redlands strain were analyzed by SnapGene 6.1.1. Amino acid divergence was represented by the exact match of amino acids (identity) or similarity in amino acid structure (similarity). The complete N, P, and L protein sequences of all available Henipavirus isolates submitted on or before 9 February 2024 were downloaded from NCBI Virus. Metadata about submitters, counties of isolation, year, and host for all sequences were collected from NCBI Virus. Only sequences that represented >99% of the full-length amino acid sequences of each protein were selected and aligned using MegAlign Pro.

ACKNOWLEDGMENTS

The authors acknowledge Roisin Relly and Matthew Hackbart for carefully editing the manuscript.
This project was funded by a grant from the Novo Nordisk Foundation and the Pandemic Antiviral Discovery (PAD) Initiative (NNF22SA0082041), NIH grant A137062, and the Washington University BJC Investigator program.
Part of Fig. 5A was created with BioRender.

SUPPLEMENTAL MATERIAL

Fig. S1 to S5 - jvi.00503-24-s0001.pdf
Sequence analysis of henipavirus N, P, and L proteins.
ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.

REFERENCES

1.
Devnath P, Wajed S, Chandra Das R, Kar S, Islam I, Masud HMAA. 2022. The pathogenesis of Nipah virus: a review. Microb Pathog 170:105693.
2.
Kaza B, Aguilar HC. 2023. Pathogenicity and virulence of henipaviruses. Virulence 14:2273684.
3.
Ching PKG, de los Reyes VC, Sucaldito MN, Tayag E, Columna-Vingno AB, Malbas FF, Bolo GC, Sejvar JJ, Eagles D, Playford G, Dueger E, Kaku Y, Morikawa S, Kuroda M, Marsh GA, McCullough S, Foxwell AR. 2015. Outbreak of henipavirus infection, Philippines, 2014. Emerg Infect Dis 21:328–331.
4.
Lo MK, Lowe L, Hummel KB, Sazzad HMS, Gurley ES, Hossain MJ, Luby SP, Miller DM, Comer JA, Rollin PE, Bellini WJ, Rota PA. 2012. Characterization of Nipah virus from outbreaks in Bangladesh, 2008-2010. Emerg Infect Dis 18:248–255.
5.
Annand EJ, Horsburgh BA, Xu K, Reid PA, Poole B, de Kantzow MC, Brown N, Tweedie A, Michie M, Grewar JD, et al. 2022. Novel Hendra virus variant detected by sentinel surveillance of horses in Australia. Emerg Infect Dis 28:693–704.
6.
Li H, Kim J-YV, Pickering BS. 2023. Henipavirus zoonosis: outbreaks, animal hosts and potential new emergence. Front. Microbiol 14:1167085.
7.
Centers for Disease Control and Prevention (CDC1999. Update: outbreak of Nipah virus—Malaysia and Singapore,1999. MMWR Morb Mortal Wkly 48:335–337. https://pubmed.ncbi.nlm.nih.gov/10366143/.
8.
Sahani M, Parashar UD, Ali R, Das P, Lye MS, Isa MM, Arif MT, Ksiazek TG, Sivamoorthy M, Nipah Encephalitis Outbreak Investigation Group. 2001. Nipah virus infection among abattoir workers in Malaysia, 1998-1999. Int J Epidemiol 30:1017–1020.
9.
Clayton BA. 2017. Nipah virus: transmission of a zoonotic paramyxovirus. Curr Opin Virol 22:97–104.
10.
Murray K, Selleck P, Hooper P, Hyatt A, Gould A, Gleeson L, Westbury H, Hiley L, Selvey L, Rodwell B, Ketterer P. 1995. A morbillivirus that caused fatal fisease in horses and humans. Science 268:94–97.
11.
Bruno L, Nappo MA, Ferrari L, Di Lecce R, Guarnieri C, Cantoni AM, Corradi A. 2022. Nipah virus disease: epidemiological, clinical, diagnostic and legislative aspects of this unpredictable emerging zoonosis. Animals (Basel) 13:159.
12.
Yuen KY, Fraser NS, Henning J, Halpin K, Gibson JS, Betzien L, Stewart AJ. 2021. Hendra virus: epidemiology dynamics in relation to climate change, diagnostic tests and control measures. One Health 12:100207.
13.
Fernandez-Aguilar X, Pujols J, Velarde R, Rosell R, Lopez-Olvera JR, Marco I, Pumarola M, Segales J, Lavin S, Cabezon O. 2014. Novel Henipa-like virus, Mojiang paramyxovirus, in rats, China, 2012 20:1062–1064.
14.
Mallapaty S. 2022. New 'Langya' virus identified in China: what scientists know so far. Nature 608:656–657.
15.
Lee SH, Kim K, Kim J, No JS, Park K, Budhathoki S, Lee SH, Lee J, Cho SH, Cho S, Lee GY, Hwang J, Kim HC, Klein TA, Uhm CS, Kim WK, Song JW. 2021. Discovery and genetic characterization of novel paramyxoviruses related to the genus Henipavirus in Crocidura species in the Republic of Korea. Viruses 13:2020.
16.
Tao W, Gan T, Guo M, Xu Y, Zhong J. 2017. Novel stable Ebola virus minigenome replicon reveals remarkable stability of the viral genome. J Virol 91:e01316-17.
17.
Ke X, Ye C, Liu R, Liu F, Chen Q. 2024. Establishment of a novel minigenome system for the identification of drugs targeting Nipah virus replication. J Gen Virol 105.
18.
Sleeman K, Bankamp B, Hummel KB, Lo MK, Bellini WJ, Rota PA. 2008. The C, V and W proteins of Nipah virus inhibit minigenome replication. J Gen Virol 89:1300–1308.
19.
Sun B, Jia L, Liang B, Chen Q, Liu D. 2018. Phylogeography, transmission, and viral proteins of Nipah virus. Virol Sin 33:385–393.
20.
Halpin K, Bankamp B, Harcourt BH, Bellini WJ, Rota PA. 2004. Nipah virus conforms to the rule of six in a minigenome replication assay. J Gen Virol 85:701–707.
21.
Hoffmann B, Schütze H, Mettenleiter TC. 2003. Recognition of cis-acting elements of infectious haematopoietic necrosis virus and viral hemorrhagic septicemia virus by homologous and heterologous helper proteins. Virus Res 93:79–89.
22.
Kim GN, Kang CY. 2005. Utilization of homotypic and heterotypic proteins of vesicular stomatitis virus by defective interfering particle genomes for RNA replication and virion assembly: implications for the mechanism of homologous viral interference. J Virol 79:9588–9596.
23.
Yunus AS, Krishnamurthy S, Pastey MK, Huang Z, Khattar SK, Collins PL, Samal SK. 1999. Rescue of a bovine respiratory syncytial virus genomic RNA analog by bovine, human and ovine respiratory syncytial viruses confirms the “functional integrity” and “cross-recognition” of BRSV cis-acting elements by HRSV and ORSV. Arch Virol 144:1977–1990.
24.
Mühlberger E, Weik M, Volchkov VE, Klenk HD, Becker S. 1999. Comparison of the transcription and replication strategies of marburg virus and Ebola virus by using artificial replication systems. J Virol 73:2333–2342.
25.
Pelet T, Marq JB, Sakai Y, Wakao S, Gotoh H, Curran J. 1996. Rescue of Sendai virus cDNA templates with cDNA clones expressing parainfluenza virus type 3 N, P and L proteins. J Gen Virol 77:2465–2469.
26.
Curran JA, Kolakofsky D. 1991. Rescue of a Sendai virus DI genome by other parainfluenza viruses- implications for genome replication. Virology 182:168–176.
27.
Dimock K, Collins PL. 1993. Rescue of synthetic analogs of genomic RNA and replicative-intermediate RNA of human parainfluenza virus type 3. J Virol 67:2772–2778.
28.
Brown DD, Collins FM, Duprex WP, Baron MD, Barrett T, Rima BK. 2005. Rescue' of mini-genomic constructs and viruses by combinations of morbillivirus N, P and L proteins. J Gen Virol 86:1077–1081.
29.
Buchholz UJ, Finke S, Conzelmann KK. 1999. Generation of bovine respiratory syncytial virus (BRSV) from cDNA- BRSV NS2 is not essential for virus replication in tissue culture, and the human RSV leader region acts as a functional BRSV genome promoter. J Virol 73:251–259.
30.
Beaty SM, Park A, Won ST, Hong P, Lyons M, Vigant F, Freiberg AN, Duprex WP, Lee B, tenOever BR. 2017. Efficient and robust paramyxoviridae reverse genetics systems. mSphere 2:e00376-16.
31.
Huang M, Sato H, Hagiwara K, Watanabe A, Sugai A, Ikeda F, Kozuka-Hata H, Oyama M, Yoneda M, Kai C. 2011. Determination of a phosphorylation site in Nipah virus nucleoprotein and its involvement in virus transcription. J Gen Virol 92:2133–2141.
32.
Laing ED, Amaya M, Navaratnarajah CK, Feng YR, Cattaneo R, Wang LF, Broder CC. 2018. Rescue and characterization of recombinant cedar virus, a non-pathogenic henipavirus species. Virol J 15:56.
33.
Welch SR, Tilston NL, Lo MK, Whitmer SLM, Harmon JR, Scholte FEM, Spengler JR, Duprex WP, Nichol ST, Spiropoulou CF. 2020. Inhibition of Nipah virus by defective interfering particles. J Infect Dis 221:S460–S470.
34.
Griffin BD, Leung A, Chan M, Warner BM, Ranadheera C, Tierney K, Audet J, Frost KL, Safronetz D, Embury-Hyatt C, Booth SA, Kobasa D. 2019. Establishment of an RNA polymerase II-driven reverse genetics system for Nipah virus strains from Malaysia and Bangladesh. Sci Rep 9:11171.
35.
Freiberg A, Dolores LK, Enterlein S, Flick R. 2008. Establishment and characterization of plasmid-driven minigenome rescue systems for Nipah virus: RNA polymerase I- and T7-catalyzed generation of functional paramyxoviral RNA. Virology 370:33–44.
36.
Habchi J, Blangy S, Mamelli L, Jensen MR, Blackledge M, Darbon H, Oglesbee M, Shu Y, Longhi S. 2011. Characterization of the interactions between the nucleoprotein and the phosphoprotein of henipavirus. J Biol Chem 286:13583–13602.
37.
Habchi J, Longhi S. 2015. Structural disorder within paramyxoviral nucleoproteins and phosphoproteins in their free and bound forms: from predictions to experimental assessment. Int J Mol Sci 16:15688–15726.
38.
Taylor J, Thompson K, Annand EJ, Massey PD, Bennett J, Eden JS, Horsburgh BA, Hodgson E, Wood K, Kerr J, Kirkland P, Finlaison D, Peel AJ, Eby P, Durrheim DN. 2022. Novel variant Hendra virus genotype 2 infection in a horse in the greater Newcastle region, New South Wales, Australia. One Health 15:100423.
39.
Alam MJ, Sultana MS, Ahmed J. 2015. Structure analysis of interacting domains of RNA dependent RNA polymerase (Rdrp) complex in Nipah virus. Biochem Physiol :4.
40.
Liang B. 2020. Structures of the mononegavirales polymerases. J Virol 94:e00175–20.
41.
Harcourt BH, Tamin A, Halpin K, Ksiazek TG, Rollin PE, Bellini WJ, Rota PA. 2001. Molecular characterization of the polymerase gene and genomic termini of Nipah virus. Virology 287:192–201.
42.
Fearns R, Plemper RK. 2017. Polymerases of paramyxoviruses and pneumoviruses. Virus Res 234:87–102.
43.
Chavda VP, Apostolopoulos V, Sah R. 2023. Langya henipavirus outbreak. Expert Rev Anti Infect Ther 21:1021–1024.
44.
Spiropoulou CF. 2019. Nipah virus outbreaks: still small but extremely lethal. J Infect Dis 219:1855–1857.
45.
Zhang H, Fischer DK, Shuda M, Moore PS, Gao SJ, Ambrose Z, Guo H. 2022. Construction and characterization of two SARS-CoV-2 minigenome replicon systems. J Med Virol 94:2438–2452.
46.
Yamada H, Taniguchi S, Shimojima M, Tan L, Kimura M, Morinaga Y, Fukuhara T, Matsuura Y, Komeno T, Furuta Y, Saijo M, Tani H. 2021. M segment-based minigenome system of severe fever with thrombocytopenia syndrome virus as a tool for antiviral drug screening. Viruses 13:1061.
47.
Su J, Dou Y, You Y, Cai X. 2015. Application of minigenome technology in virology research of the paramyxoviridae family. J Microbiol Immunol Infect 48:123–129.
48.
Vanmechelen B, Stroobants J, Vermeire K, Maes P. 2021. Advancing Marburg virus antiviral screening: optimization of a novel T7 polymerase-independent minigenome system. Antiviral Res 185:104977.
49.
Xie X, Zou J, Shan C, Yang Y, Kum DB, Dallmeier K, Neyts J, Shi PY. 2016. Zika virus replicons for drug discovery. EBioMedicine 12:156–160.
50.
Bruhn JF, Hotard AL, Spiropoulou CF, Lo MK, Saphire EO. 2019. A conserved basic patch and central kink in the Nipah virus phosphoprotein multimerization domain are essential for polymerase function. Structure 27:660–668.
51.
Sugai A, Sato H, Yoneda M, Kai C. 2017. Gene end-like sequences within the 3' non-coding region of the Nipah virus genome attenuate viral gene transcription. Virology 508:36–44.
52.
Sourimant J, Plemper RK. 2016. Organization, function, and therapeutic targeting of the Morbillivirus RNA-dependent RNA polymerase complex. Viruses 8:251.
53.
Cox R, Plemper RK. 2015. The paramyxovirus polymerase complex as a target for next-generation anti-paramyxovirus therapeutics. Front Microbiol 6:459.
54.
Bruhn JF, Barnett KC, Bibby J, Thomas JMH, Keegan RM, Rigden DJ, Bornholdt ZA, Saphire EO. 2014. Crystal structure of the Nipah virus phosphoprotein tetramerization domain. J Virol 88:758–762.
55.
Jordan PC, Liu C, Raynaud P, Lo MK, Spiropoulou CF, Symons JA, Beigelman L, Deval J. 2018. Initiation, extension, and termination of RNA synthesis by a paramyxovirus polymerase. PLoS Pathog 14:e1006889.
56.
Sourimant J, Thakkar VD, Cox RM, Plemper RK. 2020. Viral evolution identifies a regulatory interface between paramyxovirus polymerase complex and nucleocapsid that controls replication dynamics. Sci Adv 6:eaaz1590.
57.
Longhi S. 2015. Structural disorder within paramyxoviral nucleoproteins. FEBS Lett 589:2649–2659.
58.
Longhi S. 2009. Nucleocapsid structure and function. Curr Top Microbiol and Immunol 329:103–128.
59.
Houben K, Marion D, Tarbouriech N, Ruigrok RWH, Blanchard L. 2007. Interaction of the C-terminal domains of Sendai virus N and P proteins: comparison of polymerase-nucleocapsid interactions within the paramyxovirus family. J Virol 81:6807–6816.
60.
Bossart KN, Wang LF, Flora MN, Chua KB, Lam SK, Eaton BT, Broder CC. 2002. Membrane fusion tropism and heterotypic functional activities of the Nipah virus and Hendra virus envelope glycoproteins. J Virol 76:11186–11198.
61.
Caruso S, Edwards SJ. 2023. Recently emerged novel henipa-like viruses: shining a spotlight on the shrew. Viruses 15:2407.
62.
Hernández LHA, da Paz TYB, Silva SP da, Silva FS da, Barros BCV de, Nunes BTD, Casseb LMN, Medeiros DBA, Vasconcelos PF da C, Cruz ACR. 2022. First genomic evidence of a henipa-like virus in Brazil. Viruses 14:2167.
63.
Marsh GA, Todd S, Foord A, Hansson E, Davies K, Wright L, Morrissy C, Halpin K, Middleton D, Field HE, Daniels P, Wang LF. 2010. Genome sequence conservation of Hendra virus isolates during spillover to horses, Australia. Emerg Infect Dis 16:1767–1769.
64.
Castro ÍA, Yang Y, Gnazzo V, Kim D-H, Van Dyken SJ, López CB. 2023. Murine parainfluenza virus persists in lung innate immune cells sustaining chronic lung pathology. bioRxiv:2023.11.07.566103.

Information & Contributors

Information

Published In

cover image Journal of Virology
Journal of Virology
Volume 98Number 613 June 2024
eLocator: e00503-24
Editor: Martin Schwemmle, University Medical Center Freiburg, Freiburg, Germany
PubMed: 38780245

History

Received: 16 March 2024
Accepted: 24 April 2024
Published online: 23 May 2024

Keywords

  1. henipavirus
  2. minigenome system
  3. polymerase complex proteins

Data Availability

All data pertaining to this study are available in the paper and the supplemental material.

Contributors

Authors

Department of Molecular Microbiology and Center for Women's Infectious Diseases Research, Washington University in St Louis, St. Louis, Missouri, USA
Department of Molecular Microbiology and Center for Women's Infectious Diseases Research, Washington University in St Louis, St. Louis, Missouri, USA
Department of Molecular Microbiology and Center for Women's Infectious Diseases Research, Washington University in St Louis, St. Louis, Missouri, USA

Editor

Martin Schwemmle
Editor
University Medical Center Freiburg, Freiburg, Germany

Notes

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note:

  • For recently published articles, the TOTAL download count will appear as zero until a new month starts.
  • There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
  • Citation counts come from the Crossref Cited by service.

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. For an editable text file, please select Medlars format which will download as a .txt file. Simply select your manager software from the list below and click Download.

View Options

Figures

Tables

Media

Share

Share

Share the article link

Share with email

Email a colleague

Share on social media

American Society for Microbiology ("ASM") is committed to maintaining your confidence and trust with respect to the information we collect from you on websites owned and operated by ASM ("ASM Web Sites") and other sources. This Privacy Policy sets forth the information we collect about you, how we use this information and the choices you have about how we use such information.
FIND OUT MORE about the privacy policy