key: cord-0002555-nulqcpz5 authors: Hao, Sujuan; Zhang, Junmei; Chen, Zhen; Xu, Huanzhou; Wang, Hanzhong; Guan, Wuxiang title: Alternative Polyadenylation of Human Bocavirus at Its 3′ End Is Regulated by Multiple Elements and Affects Capsid Expression date: 2017-01-18 journal: J Virol DOI: 10.1128/jvi.02026-16 sha: d474315fbb8b2c284039f0caec40aaa334db068b doc_id: 2555 cord_uid: nulqcpz5 Alternative processing of human bocavirus (HBoV) P5 promoter-transcribed RNA is critical for generating the structural and nonstructural protein-encoding mRNA transcripts. The regulatory mechanism by which HBoV RNA transcripts are polyadenylated at proximal [(pA)p] or distal [(pA)d] polyadenylation sites is still unclear. We constructed a recombinant HBoV infectious clone to study the alternative polyadenylation regulation of HBoV. Surprisingly, in addition to the reported distal polyadenylation site, (pA)d, a novel distal polyadenylation site, (pA)d2, which is located in the right-end hairpin (REH), was identified during infectious clone transfection or recombinant virus infection. (pA)d2 does not contain typical hexanucleotide polyadenylation signal, upstream elements (USE), or downstream elements (DSE) according to sequence analysis. Further study showed that HBoV nonstructural protein NS1, REH, and cis elements of (pA)d were necessary and sufficient for efficient polyadenylation at (pA)d2. The distance and sequences between (pA)d and (pA)d2 also played a key role in the regulation of polyadenylation at (pA)d2. Finally, we demonstrated that efficient polyadenylation at (pA)d2 resulted in increased HBoV capsid mRNA transcripts and protein translation. Thus, our study revealed that all the bocaviruses have distal poly(A) signals on the right-end palindromic terminus, and alternative polyadenylation at the HBoV 3′ end regulates its capsid expression. IMPORTANCE The distal polyadenylation site, (pA)d, of HBoV is located about 400 nucleotides (nt) from the right-end palindromic terminus, which is different from those of bovine parvovirus (BPV) and canine minute virus (MVC) in the same genus whose distal polyadenylation is located in the right-end stem-loop structure. A novel polyadenylation site, (pA)d2, was identified in the right-end hairpin of HBoV during infectious clone transfection or recombinant virus infection. Sequence analysis showed that (pA)d2 does not contain typical polyadenylation signals, and the last 42 nt form a stem-loop structure which is almost identical to that of MVC. Further study showed that NS1, REH, and cis elements of (pA)d are required for efficient polyadenylation at (pA)d2. Polyadenylation at (pA)d2 enhances capsid expression. Our study demonstrates alternative polyadenylation at the 3′ end of HBoV and suggests an additional mechanism by which capsid expression is regulated. polyadenylation of BPV and MVC, (pA)d2 is located on the right-end terminus ( Fig. 2A) . Interestingly, the last 42 nt of the right-end stem of HBoV and MVC are almost identical, with only two different nucleotides in the loop ( Fig. 2A) . To confirm our 3= RACE results, total RNA of pHBoV1-WH-transfected HEK293T cells were harvested at different times (12, 24, 36 , and 48 h) and hybridized with oligo(dT) and specific primers upstream of (pA)d. RNA fragments were detected with probes spanning from nt 4745 to nt 5171 after RNase H cleavage. Polyadenylation at (pA)d and (pA)d2 were both detected by Northern blot analysis (Fig. 2D) , which is consistent with 3= RACE data and suggests that a portion of VP1/VP2 encoding mRNA transcripts are polyadenylated at (pA)d2. NS1 and REH are required for polyadenylation at (pA)d2. Since (pA)d2 does not contain a typical polyadenylation signal, we determined whether the viral proteins affect polyadenylation at (pA)d2. All of the structural and nonstructural protein open reading frames were disrupted by single-nucleotide mutation, which resulted in early translation termination (Fig. 3A) . Transfection of the NS1 knockout plasmid (pHBoV1-NS1KO) resulted in a 7-fold reduction in polyadenylation at (pA)d2 (Fig. 3B, lane 3 , and C, lane 2). Cotransfection of pHBoV1-NS1KO and NS1 expression plasmid (pcDNA-NS1) restored polyadenylation at (pA)d2 (Fig. 3B , lane 4, and C, lanes 8 and 13), which suggested that NS1 is necessary for efficient polyadenylation at (pA)d2. Transfection of the NP1 knockout plasmid (pHBoV1-NP1KO) to HEK293T cells resulted in decreased polyadenylation at both (pA)d and (pA)d2 as measured by using 3= RACE (Fig. 3B , lanes . PCR was performed with specific primers and products were resolved in 1.5% agarose gels. (C) 3= RACE. pHBoV1-WH was transfected to A549, HeLa, HEK293, or HEK293T cells, RNA was isolated, and 3= RACE was performed as described for panel B. Lane M, DNA ladder. (D) Northern blot analysis. Total RNA was isolated from HBoV infectious clone-transfected HEK293T cells at the indicated times and hybridized with oligo(dT) and HBoV-specific primers, followed by RNase H cleavage. RNAs were resolved in 1.5% agarose gels, transferred to Hybond-N ϩ membranes, and hybridized with probes spanning nt 4745 to nt 5171. 5 ) and Northern blot analysis (Fig. 3C, lanes 3 and 9) . However, the ratio of polyadenylation at (pA)d2 to (pA)d was not changed. Cotransfection of pHBoV1-NP1KO and NP1 expression plasmid (pXJ40-NP1) did not change the ratio of polyadenylation at (pA)d2 to (pA)d (Fig. 3B, lane 6 , and C, lane 10). However, the abundance of RNA transcripts polyadenylated at both (pA)d and (pA)d2 (Fig. 3B, lane 6 , and C, lane 10) increased more than 10-fold, which is consistent with the previous report that NP1 facilitates capsid protein expression (33) . These results suggested that NP1 is not required for polyadenylation at (pA)d2. Parvovirus VP1 is required for viral infection, and VP2 is the major capsid protein which is involved in viral particle assembly. Transfecting VP1 and VP2 knockout plasmids (pHBoV1-VP1KO and pHBoV1-VP2KO) to HEK293T cells did not result in a change of polyadenylation at (pA)d2 (Fig. 3B , lanes 7 and 8, and C, lanes 4 and 5), which suggested that VP1/VP2 is not involved in the regulation of its 3=-end alternative polyadenylation. (pA)d2 is located in the loop of the right-end palindromic terminus. We next checked whether the hairpins at 5= and/or 3= ends of HBoV are necessary for polyadenylation at (pA)d2. Even deleting all of the left-end hairpin did not affect polyadenylation at (pA)d2 (Fig. 4B , lanes 3 to 6, and C, lanes 2 and 3). Transfecting pHBoV1REH-98, in which the sequences after the (pA)d2 site (nt 5445) were truncated, to HEK293T cells resulted in decreased polyadenylation at (pA)d2, and the (pA)d2/(pA)d ratio was reduced 2-fold (Fig. 4B, lane 7) . Polyadenylation only at (pA)d was observed when pHBoVREH-161, in which the sequences before (pA)d2 were trimmed, was transfected to HEK293T cells (Fig. 4B, lane 8) . The result indicated that the right hairpin plays an important role in polyadenylation at (pA)d2. Taken together, these data show that NS1 and REH are required for polyadenylation at (pA)d2, while NP1 and capsid proteins do not regulate polyadenylation at (pA)d2. cis elements of (pA)d are required for polyadenylation at (pA)d2. To test further which elements are required for efficient polyadenylation at (pA)d2, we constructed an To investigate whether the cis elements of (pA)d are required for polyadenylation at (pA)d2, (pA)d and downstream sequences were replaced with heterologous sequences, , which suggested that cis elements of (pA)d are required for polyadenylation at (pA)d2. We then mutated the hexanucleotide of the (pA)d (pEGFP-mPAS). No polyadenylation at (pA)d2 was observed when pEGFP-mPAS was transfected into cells in the absence or presence of NS1 cotransfection (Fig. 5B , lanes 11 and 12, and C, lanes 11 and 12). Replacing (pA)d and downstream sequences with a synthetic polyadenylation signal (sPA) resulted in strong polyadenylation at sPA without polyadenylation at (pA)d2, even in the presence of NS1 cotransfection (Fig. 5B , lanes 9 and 10, and C, lanes 9 and 10). Those results indicated that the addition of the sPA created a strong enough polyadenylation signal that the majority of RNA transcripts were polyadenylated at (pA)d, and few RNA transcripts were polyadenylated at (pA)d2. Collectively, these experiments showed that NS1, REH, and cis elements of (pA)d are sufficient viral elements for efficient polyadenylation at (pA)d2. DSE of (pA)d regulate alternative polyadenylation at (pA)d and (pA)d2. We then determined whether downstream elements (DSE) of (pA)d affect polyadenylation at (pA)d2 in the context of an infectious clone. Sequence analysis showed that a U-rich stretch is 22 nt downstream of hexanucleotide AAUAAA. Transfecting the U-rich mutation plasmids (pHBoV1mDSE1 and pHBoV1mDSE2) into HEK293T cells resulted in reduced polyadenylation at both (pA)d and (pA)d2, as determined by using Northern blotting (Fig. 6B, lanes 3 and 4) , suggesting that this motif is a major component of downstream elements for efficient polyadenylation at (pA)d and (pA)d2. Mutating the sequences between the U-stretch and AAUAAA resulted in reduced polyadenylation at (pA)d, but polyadenylation at (pA)d2 was not changed (Fig. 6B, lane 5 ). This result shows that polyadenylation at (pA)d is not required for polyadenylation at (pA)d2. There is also a difference in the requirement of cis elements for polyadenylation at the two different sites. The sequences between nt 5159 and nt 5196 are required for efficient polyadenylation at (pA)d. Polyadenylation at (pA)d2 only requires sequences between nt 5178 and nt 5196. Replacing the (pA)d and downstream sequences with a synthetic polyadenylation signal resulted in strong polyadenylation at (pA)d but a loss of polyadenylation at (pA)d2 (Fig. 6B, lane 6) , which further confirmed that a strong polyadenylation at (pA)d inhibited polyadenylation at (pA)d2. The sequences and distance between (pA)d and (pA)d2 affect polyadenylation at (pA)d2. We next investigated whether the distance between (pA)d and (pA)d2 plays a role in the regulation of alternative polyadenylation at the 3= end of HBoV. pHBoV1Del1 and pHBoV1Del2, in which the sequences from nt 5223 or nt 5197 to nt 5356 were deleted, were transfected into HEK293T cells. 3= RACE and Northern blot analysis showed that neither polyadenylation at (pA)d2 nor the (pA)d2/(pA)d ratio changed (Fig. 7B, lanes 3 and 4, and C, lanes 2 and 3) , which indicates shortening the distance does not affect polyadenylation at (pA)d2. As the distance between the two Polyadenylation at (pA)d2 enhances capsid mRNA and protein expression. The outcome of alternative polyadenylation is to regulate the expression level of encoded proteins. To determine whether polyadenylation at (pA)d2 affects protein expression, we first measured whether (pA)d2 affects green fluorescent protein (GFP) expression in the reporter system. GFP expression did not change obviously ( Fig. 8B and C, lanes 3 and 4) when pEGFP-(pA)d was transfected into HEK293T cells in the presence of NS1 cotransfection. However, in the presence of NS1 expression, transfection of pEGFP-(pA)d2 resulted in polyadenylation at both (pA)d and (pA)d2, and the expression of GFP increased as determined using immunofluorescence and Western blot analysis ( Fig. 8B and C, lanes 5 and 6). Transfecting pEGFP-HS resulted in the loss of polyadenylation at (pA)d2 and a decrease in GFP expression ( Fig. 8B and C, lanes 7 and 8) , suggesting that polyadenylation at (pA)d2 regulates mRNA and protein expression levels. To analyze whether polyadenylation affects capsid protein expression, a VP2 expression cassette and the downstream sequence of (pA)d, with or without (pA)d2, were cloned into expression vector pXJ40-FLAG. Cotransfecting pXJ40-VP2(pA)d and an NS1 expression plasmid into HEK293T cells resulted in polyadenylation at (pA)d and the polyadenylation site of pXJ40-FLAG vector (Fig. 9B, lanes 2 and 3) . No polyadenylation at (pA)d2 was detected because of the lack of REH (Fig. 9B , lanes 2 and 3, and C, lanes 1 and 2). Similar levels of VP2 mRNA and protein were observed (Fig. 9C, lanes 1 and 2, and D, lanes 1 and 2) . However, cotransfecting pXJ40-VP2(pA)d2 with NS1 into HEK293T cells resulted in polyadenylation at both (pA)d and (pA)d2 (Fig. 9B, lane 5, and C, lane 4) . The level of mRNA transcripts increased more than 20-fold when NS1 was cotransfected with pXJ40-VP2 (pA)d2 (Fig. 9C, lanes 3 and 4) . The expression level of VP2 protein increased at least 10-fold (Fig. 9D, lanes 3 and 4) , suggesting that polyadenylation at (pA)d2 increased the mRNA and protein expression level of capsid. Taken together, these data suggest that alternative polyadenylation at the 3= end of HBoV regulates the capsid mRNA transcript levels and results in increased capsid protein expression. We discovered a new polyadenylation site, (pA)d2, in the right-end hairpin of human bocavirus. The nonstructural protein NS1, the right-end hairpin, and cis elements of (pA)d are all required for efficient polyadenylation at (pA)d2. The sequences and distance between (pA)d and (pA)d2 also affect polyadenylation efficiency at (pA)d2. Increased polyadenylation at (pA)d2 resulted in elevated capsid expression, indicating that alternative polyadenylation at the 3= end of HBoV is important for virus structural protein expression. The discovery of new distal polyadenylation site (pA)d2. Alternative polyadenylation plays an important role in the parvovirus life cycle. The blockage of RNA transcript maturation by the proximal polyadenylation site is a limiting step of parvovirus infection (29, 30, 34, 35) . When B19V infects permissive cells, most RNA transcripts read through proximal polyadenylation sites, and full-length mRNA transcripts are generated to encode capsid proteins for the production of progeny virus. However, when B19V infects nonpermissive cells, the majority of RNA transcripts are polyadenylated at (pA)p, causing a lack of capsid mRNA, resulting in aborted infection. Internal polyadenylation is also a limiting step of Aleutian mink disease virus (AMDV) genome replication and progeny virus production (30) . We found that the recombinant HBoV infectious clone contains three proximal polyadenylation sites, which is consistent with previous reports (22, 28) . Interestingly, in addition to the reported distal polyadenylation site, (pA)d, a new polyadenylation site, (pA)d2, was discovered when recombinant infectious clones were transfected into HEK293T cells. In contrast to the distal polyadenylation sites of BPV and MVC, which are located in the right-end palindromic terminus, the (pA)d of HBoV is located about 400 nt upstream of the right-end hairpin. VP1/VP2-encoding mRNAs are efficiently polyadenylated downstream of the AAUAAA site. Polyadenylation signals include the hexanucleotide AAUAAA or its variants, the cleavage site CA dinucleotide, downstream U-rich elements, and G/U-rich upstream elements. However, (pA)d2 of HBoV is located in the loop of REH, and no classical hexanucleotide polyadenylation signal was found around (pA)d2. The last 42 nucleotides of the stem-loop structure of HBoV and MVC are almost identical. It is possible that HBoV and MVC use similar mechanisms for polyadenylation at the stem-loop region. The only puzzling thing is that there is no polyadenylation site immediately upstream of (pA)d2. To our knowledge, this is the first report of alternative polyadenylation at a distal polyadenylation site in parvoviruses. Regulation of polyadenylation at (pA)d2 by cis and trans elements. HBoV encodes at least three nonstructural proteins that include NS1, NS1-70K, and NP1. NS1 is a multifunctional protein expressed from the mRNA transcripts from the P5 promoter of the left half of the genome. NS1 contains DNA binding, endonuclease, ATPase, and helicase activities. Disrupting the NS1 open reading frame by point mutation resulted in the loss of polyadenylation at (pA)d2, suggesting that NS1 is indispensable for polyadenylation at (pA)d2. NS2, NS3, and NS4 are dispensable for viral replication (36) . However, NS2 is necessary for efficient replication in primary human airway epithelium cultured at an air-liquid interface (HAE-ALI). Polyadenylation at (pA)d2 was not detected after cotransfection of NS2, NS3, NS4, and NS1-70K expression plasmids with (pA)d2 polyadenylation reporter plasmid pEGFP-(pA)d2 in HEK293T cells (data not shown). These results suggested that all of the functional NS1 domains are required for efficient polyadenylation at (pA)d2. NP1 is unique among parvoviruses and is essential for viral genome DNA replication (21, 32) . NP1 facilitates (pA)p read-through of structural protein-encoding mRNA transcripts (24, 25) and increases capsid protein expression. Disrupting the NP1 open reading frame did not affect polyadenylation at (pA)d and (pA)d2, while the RNA transcript abundance decreased. Cotransfection of NP1-expressing plasmid pXJ40-NP1 with pHBoV1-NP1KO did not change the ratio of polyadenylation at (pA)d2 to (pA)d. However, the total RNA level increased, which indicates NP1 regulates the abundance of mRNA transcripts but not polyadenylation at (pA)d2. Structural proteins VP1/VP2 have no effect on polyadenylation at (pA)d2. Since no typical polyadenylation sequence was found before the (pA)d2 site, whether the cis elements of (pA)d are important for efficient polyadenylation at (pA)d2 was determined by using mutagenesis. Mutating either the hexanucleotide AAUAAA or the DSE resulted in decreased polyadenylation at both (pA)d and (pA)d2, indicating that AAUAAA and DSE of (pA)d are indispensable for the efficient polyadenylation at (pA)d2. Interestingly, we found that polyadenylation at (pA)d2 was detected when polyadenylation at (pA)d was abolished by mutation of sequences from nt 5159 to nt 5177, which suggested that polyadenylation at (pA)d is not required for polyadenylation at (pA)d2. The downstream elements required for polyadenylation at (pA)d and (pA)d2 are also different. The distance and sequence per se between (pA)d and (pA)d2 also affected polyadenylation efficiency at (pA)d2. Polyadenylation at (pA)d2 was decreased when the distance was increased or the sequence was replaced with heterologous sequence. Thus, polyadenylation at (pA)d2 depends on the cis elements of (pA)d. Mutating the right-end hairpin structure destroyed the stem-loop structure where (pA)d2 is located and resulted in the loss of polyadenylation at (pA)d2, indicating that the hairpin structure also plays an important role in polyadenylation at (pA)d2. (pA)d2 function. Alternative polyadenylation is widely spread in metazoan proteincoding transcripts and leads to variable 3= untranslated regions (UTRs), which has been shown to regulate gene expression (37, 38) . In the presence of the NS1 expression plasmid, transfecting pXJ40-VP2(pA)d2 resulted in polyadenylation at (pA)d and (pA)d2. We found that polyadenylation at (pA)d2 increased not only the abundance of VP1/VP2 transcripts but also the expression level of VP1/VP2 proteins. However, our results do not immediately suggest a mechanism to explain the increase in capsid protein expression, a topic we are currently investigating. The distance between (pA)d and (pA)d2 is 274 nt and affects polyadenylation efficiency at (pA)d2. This 3= noncoding region (NCR) of HBoV has been predicted to form two conserved hairpin structures and plays an important part in the replication of bocaviruses (39) (40) (41) . The 3= NCR could function as a binding site for miRNAs or be involved in the production of noncoding RNAs to regulate protein expression (42, 43) . The sequences of this region enhanced polyadenylation at (pA)d2, but the mechanism is currently unknown. Taken together, our study has demonstrated that alternative polyadenylation of HBoV capsid-encoding mRNA transcripts is regulated both by the nonstructural protein NS1 and by multiple cis elements. Cell lines and cell culture. HEK293T (ATCC, CRL-11268), HEK293 (ATCC, CRL-1573), A549 (ATCC, CCL-185), and HeLa cells (ATCC, CCL-2) were cultured in Dulbecco's modified Eagle's medium (DMEM; GIBCO) containing 10% fetal bovine serum (FBS; GIBCO) at 37°C with 5% CO 2 . Human bronchial epithelial Calu-3 cells were cultured in DMEM-nutrient mixture F-12 Hams (1:1) supplemented with 1% MEM nonessential amino acids solution (GIBCO), 1% GlutaMAX (GIBCO), and 10% FBS (GIBCO) at 37°C with 5% CO 2 . Calu-3 cells were seeded at a density of 5 ϫ 10 5 cells/cm 2 on semipermeable membrane inserts (0.6 cm 2 ; Millicell-PCF; Millipore). After 3 to 4 weeks of culture for differentiation, cells with transepithelial electrical resistance (TEER) values over 600 V/cm 2 were used for HBoV infection. Transfection. Plasmids (2 g) were transfected into cells plated on 60-mm dishes with Lipofectamine 2000 reagent (Invitrogen, Life Technologies) according to the manufacturer's instructions. Virus production and purification. The infectious clone pHBoV1-WH (10 g) was transfected into HEK293T cells seeded on 100-mm plates with Lipofectamine 2000. Forty-eight hours posttransfection, the cells were collected and lysed by three rounds of freezing and thawing. The cell lysate was then spun at 10,000 rpm for 30 min after Benzonase (Sigma) treatment for 30 min at 37°C. The supernatant was collected and further purified through discontinuous step gradients of iodixanol, prepared using a 60% (wt/vol) sterile solution of OptiPrep (Axis-Shield) as described previously (44) . Viral DNA was extracted using a QIAamp blood minikit (Qiagen) and quantified using quantitative PCR as described previously (21) . Virus infection. Differentiated Calu-3 cells in Millicell inserts were incubated with HBoV1 virus purified from transfected HEK293T cells at a multiplicity of infection (MOI) of 100 g/cell at 37°C for 2 h, followed by three washes with phosphate-buffered saline (PBS). The cells were then cultured for HBoV RNA analysis. Plasmid construction. (i) HBoV1 infectious clone construction. The HBoV genome was prepared from the feces of a 2-year-old infant using a viral DNA extraction kit (Qiagen) according to the manufacturer's instructions. The nearly full-length genome (lacking the left-end hairpin [LEH] and right-end hairpin [REH] sequences) was amplified using primers (XhoI-F, 5=-GCGCTCGAGCCGGCAGACA TATTGGATTC-3=; XbaI-R, 5=-GCGTCTAGATGTACAACAACAACACATTAAAAGATATAGAG-3=). LEH and REH were synthesized according to the reference sequence JQ923422. The full-length recombinant HBoV clone was ligated into the pBBSmal vector and designated pHBoV1-WH (21, 45) . (ii) Construction of pHBoV1-WH mutants. All of the nonstructural and structural protein knockout plasmids were constructed based on the recombinant HBoV clone pHBoV1-WH by single-nucleotide mutation that resulted in early translation termination of the open reading frame. NS1, NP1, VP1, and VP2 open reading frames were disrupted by mutating nt 542 from T to A, nt 2588 from G to A, nt 3205 from T to A, and nt 3540 from T to G, respectively (Fig. 3A) . (iii) Inverted terminal repeat (ITR) mutants. LEH truncation plasmids pHBoV1LEH-39, pHBoV1LEH-70, pHBoV1LEH-96, and pHBoV1LEH-140 were constructed by deleting the first 39, 70, 96, and 140 nt, respectively, of HBoV in pHBoV1-WH. Similarly, the right-end mutant plasmids pHBoV1REH-98 and pHBoV1REH-161 were constructed by truncating the last 98 and 161 nt of HBoV in pHBoV1LEH-140, which does not contain the LEH. (iv) Reporter constructs to test cis elements required for polyadenylation at (pA)d2. The parent polyadenylation reporter plasmid pEGFP-C1(-SV40 PA) was generated by deleting the SV40 polyadenylation signal of pEGFP-C1. The (pA)d2 reporter plasmids pEGFP-(pA)d and pEGFP-(pA)d2 were constructed by inserting the HBoV sequences of nt 5121 to 5271 and nt 5121 to 5543, respectively, into pEGFP-C1(-SV40 PA). pEGFP-HS was constructed by replacing nt 5121 to 5355 with the heterologous sequences from the open reading frame of kanamycin spanning from nt 57 to nt 291 based on plasmid pEGFP-(pA)d2. pEGFP-sPA was generated by changing nt 5153 to 5201 with a synthesized polyadenylation signal (5=-AATAAAAGATCTTTATTTTCATTAGATCTGTGTGTTGGTTTTTTGTGTG-3=) based on plasmid pEGFP-(pA)d2. pEGFP-mPAS was constructed by replacing (pA)d hexanucleotide AAUAAA with CCTAGG based on pEGFP-(pA)d2. (v) Constructs to identify the DSE of polyadenylation at (pA)d site. pHBoV1-WH(NheI-SmaI) was constructed by mutating nt 4918 to 4921 from CTCA to TAGC and nt 5350 to 5355 from ATGCAA to CCCGGG, resulting in NheI and SmaI restriction enzyme sites in pHBoV1-WH. pHBoV1mDSE1, pHBoV1mDSE2, pHBoV1mDSE3, and pHBoV1mDSE4 were constructed by mutating nt 5159 to 5196, nt 5177 to 5196, nt 5159 to 5177, and nt 5159 to 5201 to 5=-CACTGACTGTCAGTCTGAGTCAGAGTCACTGA CTCAGC-3=, 5=-TCAGAGTCACTGACTCAGC-3=, 5=-CACTGACTGTCAGTCTGAG-3=, and synthetic polyadenylation signal (5=-AGATCTTTATTTTCATTAGATCTGTGTGTTGGTTTTTTGTGTG-3=), respectively, in pHBoV1-WH(NheI-SmaI) (46, 47) . (vi) Constructs to analyze effect of the distance between (pA)d and (pA)d2 on polyadenylation at (pA)d2 site. pHBoV1Del1 and pHBoV1Del2 plasmids were made by deleting nt 5223 to nt 5356 and nt 5197 to nt 5356 on pHBoV1-WH(NheI-SmaI) to shorten the distance between (pA)d and (pA)d2. pHBoV1mut1 was mutated from nt 5223 to nt 5356 with a kanamycin open reading frame from nt 4 to nt 117. pHBoV1HS1 and pHBoV1HS2 were constructed by replacing the sequences from nt 5223 to nt 5356 with a kanamycin open reading frame from nt 4 to nt 217 and nt 4 to nt 617, respectively. (vii) Construction of NS1, NP1, and VP2 expression plasmids. NS1 expression plasmid pcDNA-NS1 was constructed by inserting the NS1 ORF into pcDNA3.1 vector with a C-terminal hemagglutinin (HA) tag. NP1 ORF, VP2 ORF, and downstream sequences were inserted into pXJ40-FLAG to generate NP1 (pXJ40-NP1) and VP2 expression plasmids [pXJ40-VP2(pA)d and pXJ40-VP2(pA)d2]. The plasmid pXJ40-VP2(pA)d contains only the (pA)d sequence. The plasmid pXJ40-VP2(pA)d2 contains both (pA)d and (pA)d2 sequences (Fig. 9A) . Southern blot analysis. Isolation of low-molecular-weight (Hirt) DNA and Southern analysis were performed as previously described (29, 48) with digoxigenin (DIG)-labeled DNA probes (Roche). RNA isolation. Total RNA from transfected cells was harvested using TRIzol reagent (Ambion) according to the manufacturer's instructions. 3= RACE. Reverse transcription was performed with total RNA from transfected or infected cells with oligo(dT) primers (5=-GCTATCATCACAATGGACTTTTTTTTTTTTTTTTTTV-3=). PCR amplification was performed using an adaptor primer (5=-GCTATCATCACAATGGAC-3=) and a gene-specific primer (GSP1, 5=-CCAACTAACCTAGAATACAAACTTC-3=). The second round of PCR was performed with the adaptor primer and another gene-specific primer (GSP2, 5=-CAGATTTCCTATCACCAGAG-3=), and the PCR products were gel purified with the gel extraction kit (Omega) according to the manufacturer's protocol and then sequenced. RNase H cleavage and Northern blotting. Total RNA (20 g) was combined with 100 pmol oligo(dT) (5=-TTTTTTTTTTTTTTTTTTTT-3=) and a gene-specific primer (5=-CATCCATATGTCCCCCACTA-3=). The mixture was incubated at 65°C for 5 min, followed by slow cooling to room temperature. RNase H buffer and enzyme (5 U) were added and samples were incubated at 37°C for 1 h. RNase H-treated samples were ethanol precipitated for at least 30 min at Ϫ20°C and then centrifuged. The RNA pellet was dissolved and run on 1.5% agarose gel containing 2.2 M formaldehyde for 12 h at 28 V. RNA was transferred to Hybond-N ϩ membrane by semidry transferring and then cross-linked using UV. Probe detection was performed by using the DIG luminescence detection kit II (Roche) according to the manufacturer's protocol. Signals were detected with the ChemiDoc MP imaging system (Bio-Rad). Western blotting. Cell lysates were prepared 48 h after transfection and were separated by SDS-12% PAGE, followed by transfer to nitrocellulose. Protein detection was carried out using standard protocols with anti-HA antibody (66006; Proteintech), anti-GFP (66002; Proteintech), anti-FLAG (F1804; Sigma), and anti-actin (sc-47778; Santa Cruz Biotechnology). Luminescent signals were detected with the ChemiDoc MP imaging system (Bio-Rad). Statistical analysis. 3= RACE was repeated at least three times. The ratio of polyadenylation at (pA)d2 to (pA)d was quantified. The averages and standard deviations are presented. We thank Zhi Ning for providing the B19V infectious clone PM20 and all of the members of the laboratory of W.G. for discussions and critical reading. The study was supported by the National Natural Science Foundation of China (31270208 to W.G.). The funders had no role in the design, interpretation, or submitting of this work for publishing. Cloning of a human parvovirus by molecular screening of respiratory tract samples The family Parvoviridae Human bocavirus types 1, 2 and 3 in acute gastroenteritis of childhood Detection of human bocaviruses 1 to 4 from nasopharyngeal swab samples collected from patients with respiratory tract infections Detection of human bocavirus 1 and 2 from children with acute gastroenteritis in Japan Circulating of human bocavirus 1, 2, 3, and 4 in pediatric patients with acute gastroenteritis in Thailand Differential seroprevalence of human bocavirus species 1-4 in Beijing Detection of human bocaviruses 1, 2 and 3 in Irish children presenting with gastroenteritis Detection of human bocavirus 3 in China Seroepidemiology of human bocaviruses 1-4 Human bocavirus species 2 and 3 in Brazil Detection of human bocavirus-2 in children with acute gastroenteritis in South Korea A novel bocavirus associated with acute gastroenteritis in Australian children Human bocavirus in Italian patients with respiratory diseases Frequent detection of bocavirus DNA in German children with respiratory tract infections Evidence of human coronavirus HKU1 and human bocavirus in Australian children Epidemiological profile and clinical associations of human bocavirus and other human parvoviruses Human bocavirus in French children Human bocavirus: a novel parvovirus epidemiologically associated with pneumonia requiring hospitalization in Thailand Establishment of a reverse genetics system for studying human bocavirus in human airway epithelia Characterization of the gene expression profile of human bocavirus Structure of the NS1 protein N-terminal origin recognition/nickase domain from the emerging human bocavirus Characterization of the nonstructural proteins of the bocavirus minute virus of canines NP1 protein of the bocaparvovirus minute virus of canines controls access to the viral capsid genes via its role in RNA processing Phospholipase A2-like activity of human bocavirus VP1 unique region VP1u phospholipase activity is critical for infectivity of full-length parvovirus B19 genomic clones Human bocavirus can be cultured in differentiated human airway epithelial cells Block to the production of full-length B19 virus transcripts by internal polyadenylation is overcome by replication of the viral genome Internal polyadenylation of parvoviral precursor mRNA limits progeny virus production The transcription profile of the bocavirus bovine parvovirus is unlike those of previously characterized parvoviruses Molecular characterization of infectious clones of the minute virus of canines reveals unique features of bocaviruses Nonstructural protein NP1 of human bocavirus 1 plays a critical role in the expression of viral capsid proteins A block in full-length transcript maturation in cells nonpermissive for B19 parvovirus Identification and characterization of two internal cleavage and polyadenylation sites of parvovirus B19 RNA Identification and functional analysis of novel nonstructural proteins of human bocavirus 1 Alternative cleavage and polyadenylation: the long and short of it Proliferating cells express mRNAs with shortened 3= untranslated regions and fewer microRNA target sites Analysis of the cis and trans requirements for DNA replication at the right end hairpin of the human bocavirus 1 genome A study of the human bocavirus replicative genome structures Detection of a bocavirus circular genome in fecal specimens from children with acute diarrhea in Beijing, China Functions of the 3= and 5= genome RNA regions of members of the genus Flavivirus A highly structured, nuclease-resistant, noncoding RNA produced by flaviviruses is required for pathogenicity Recombinant adenoassociated virus purification using novel methods improves infectious titer and yield Molecular and functional analyses of a human parvovirus B19 infectious clone demonstrates essential roles for NS1, VP1, and the 11-kilodalton protein in virus replication and infectivity Definition of an efficient synthetic poly(A) site Internal polyadenylation of the parvovirus B19 precursor mRNA is regulated by alternative splicing The genome of human parvovirus b19 can replicate in nonpermissive cells with the help of adenovirus genes and produces infectious virus