Profiling of rotavirus 3′UTR-binding proteins reveals the ATP synthase subunit ATP5B as a host factor that supports late-stage virus replication

Genome replication and virion assembly of segmented RNA viruses are highly coordinated events, tightly regulated by sequence and structural elements in the UTRs of viral RNA. This process is poorly defined and likely requires the participation of host proteins in concert with viral proteins. In this study, we employed a proteomics-based approach, named RNA–protein interaction detection (RaPID), to comprehensively screen for host proteins that bind to a conserved motif within the rotavirus (RV) 3′ terminus. Using this assay, we identified ATP5B, a core subunit of the mitochondrial ATP synthase, as having high affinity to the RV 3′UTR consensus sequences. During RV infection, ATP5B bound to the RV 3′UTR and co-localized with viral RNA and viroplasm. Functionally, siRNA-mediated genetic depletion of ATP5B or other ATP synthase subunits such as ATP5A1 and ATP5O reduced the production of infectious viral progeny without significant alteration of intracellular viral RNA levels or RNA translation. Chemical inhibition of ATP synthase diminished RV yield in both conventional cell culture and in human intestinal enteroids, indicating that ATP5B positively regulates late-stage RV maturation in primary intestinal epithelial cells. Collectively, our results shed light on the role of host proteins in RV genome assembly and particle formation and identify ATP5B as a novel pro-RV RNA-binding protein, contributing to our understanding of how host ATP synthases may galvanize virus growth and pathogenesis.

Genome replication and virion assembly of segmented RNA viruses are highly coordinated events, tightly regulated by sequence and structural elements in the UTRs of viral RNA. This process is poorly defined and likely requires the participation of host proteins in concert with viral proteins. In this study, we employed a proteomics-based approach, named RNA-protein interaction detection (RaPID), to comprehensively screen for host proteins that bind to a conserved motif within the rotavirus (RV) 3 terminus. Using this assay, we identified ATP5B, a core subunit of the mitochondrial ATP synthase, as having high affinity to the RV 3UTR consensus sequences. During RV infection, ATP5B bound to the RV 3UTR and co-localized with viral RNA and viroplasm. Functionally, siRNA-mediated genetic depletion of ATP5B or other ATP synthase subunits such as ATP5A1 and ATP5O reduced the production of infectious viral progeny without significant alteration of intracellular viral RNA levels or RNA translation. Chemical inhibition of ATP synthase diminished RV yield in both conventional cell culture and in human intestinal enteroids, indicating that ATP5B positively regulates late-stage RV maturation in primary intestinal epithelial cells. Collectively, our results shed light on the role of host proteins in RV genome assembly and particle formation and identify ATP5B as a novel pro-RV RNAbinding protein, contributing to our understanding of how host ATP synthases may galvanize virus growth and pathogenesis.
Viruses are obligate intracellular pathogens that usurp host cellular machinery for efficient replication and production of progeny infectious particles. Despite the availability of multiple genetic and biochemical tools that enable the examination of rotavirus-host interactions in a unbiased manner, including genome-wide siRNA screens (1,2), interactome analysis of protein-protein interactions (3), and the recently developed CRISPR-Cas9 screens (4), the molecular interactions between rotaviral RNA and host proteins remain relatively poorly understood. RNA-binding proteins play many important roles in cellular functions (5). With the exception of cytoplasmic RNA sensors, RIG-I-like receptors, including RIG-I, MDA5, and LGP2, and effector molecules such as RNase L and protein kinase R (6 -9), the identities of host proteins that specifically bind to virus RNA genome or other replication intermediates and by-products, have rarely been characterized. Lately, interrogation of stem-loop structures within the untranslated regions (UTRs) of Zika virus and dengue virus has led to the discovery of host proteins Musashi-1 and TRIM25 that play critical roles in the replication cycle of these flaviviruses (10,11). It was also found that, for influenza virus and coxsackievirus, defective RNA secondary structures markedly diminish binding of host exosome complex and PCBP2, respectively, and result in reduced replication (12,13).
Rotaviruses (RVs) 4 are nonenveloped, segmented dsRNA viruses in the Reoviridae family (14). RV infections are the leading cause of infantile diarrhea and severe gastroenteritis, resulting in around 215,000 deaths annually worldwide (15). Similar to those of influenza virus and orthoreovirus, RV gene segments have to be precisely organized and packaged into newly synthesized icosahedral virus particles during genome assembly (16). A previous study of a highly-related bluetongue virus and recent studies of RV suggest that the smaller RV RNA segments may be sorted into the newly forming virion particles first, and then the larger RNA segments enter the particle via RNA-RNA interactions in a sequential and NSP2-dependent manner (17)(18)(19). Hydrodynamic studies that examined the "stiffness" of RV RNA segments in vitro also suggested that packaging of the viral RNA segments into the capsid likely necessitates intimate RNA-protein interactions (20). Given the precision required and the likely high-energy consumption needed for RNA packaging, we hypothesize that host RNA-binding proteins contribute to the correct assembly of RV gene segments.
For most organisms, ATP hydrolysis powered by the ATPase machinery is the single most important "energy currency" (21,22). For bacteriophages and large DNA viruses like herpesvirus and poxviruses, a central component of the packaging motor that drives viral genome assembly is the ATPase subunit, provided by the viruses themselves (23)(24)(25). In contrast, RNA viruses, except bluetongue virus (26), by and large harbor a relative small genome and rarely encode viral ATPases (27), inviting the question whether or not these viruses hijack the host ATPase complex as an alternative strategy to obtain energy for genome packaging. Here, we employed a novel and powerful technique named RNA-Protein Interaction Detection (RaPID), recently developed to study RNA-protein interactions (28), to comprehensively profile the host factors that bind to a stretch of conserved sequences within 3ЈUTR of group A RV genomes. Surprisingly, we identified ATP5B, an integral part of the mitochondrial F 1 -F 0 -ATPase complex (29), as a cellular component that co-precipitated and co-localized with RV dsRNA during infection. Functional dissection using small interfering RNA (siRNA) knockdown and a panel of small-molecule pharmacological inhibitors suggested that ATP5B assists RV genome replication and virion assembly. Thus, our study systematically interrogated the host proteins that interact with the RV 3Ј terminus and revealed a tractable method to rapidly identify host proteins that bind to viral RNA sequences of interest in living cells.

Proteomic analysis reveals novel host proteins that interact with RV 3UTR consensus sequences
For group A human and animal RVs, the last seven nucleotides within the mRNA 3ЈUTR, 5Ј-UGUGACC-3Ј, are highly conserved in all 11 gene segments (30 -32) and distinct from those in group B (5Ј-UAUACCC-3Ј) and group C (5Ј-UGUG-GCU-3Ј) RVs (Fig. 1). This short sequence forms a cis-acting signal that contributes to the efficient synthesis of minus strand RNA and is also conducive to viral gene expression in the host cells through interaction with RV protein NSP3 (33,34). However, the nature and identity of host factors that bind to this important RNA region remain unknown.
To comprehensively identify host proteins that interact with RV 3ЈUTR consensus sequence, we took advantage of a new screening approach, named RaPID, that detects, with high sensitivity, any protein in the vicinity of the target RNA molecule (28). In brief, biotin-protein ligase BirA (36) is N-terminally fused to a 22-amino acid N peptide that recognizes bacteriophage BoxB stem-loops with high affinity (37), flanking both 5Ј and 3Ј ends of an RNA of interest ( Fig. 2A). Thus, BirA ligase is brought in close proximity to and biotinylates all proteins that directly or indirectly associate with target RNA. We then perform immunoprecipitation with streptavidin beads and used MS and bioinformatics analysis to identify interacting protein partners. As a proof of principle, we first tested the interaction between a pair of positive controls: a 15-nucleotide UG-rich EDEN15 motif (UGUUUGUUUGUUUGU) that is reported to bind to the CUG triplet repeat RNA-binding protein 1 (CUG-BP1) (38). As expected, specific immunoprecipitation of endogenous CUG-BP1 was only observed with EDEN15 and not with the negative control scrambled sequences (Fig. 2B). Next, we examined the well-documented interaction between the conserved RV 3ЈUTR and an RV nonstructural protein 3 (NSP3). Consistent with previous reports (34), we detected strong binding of the UGUGACC sequence with either ectopically expressed GFP-tagged NSP3 or endogenous NSP3 expressed during RV infection (Fig. 2, C and D), suggesting that this technique is well-suited to interrogate the host protein-binding partners of the RV 3ЈUTR sequences. In addition to examining binding proteins of 3ЈUTR RV monomers, we further tested whether a pentamer of RV 3ЈUTRs with two adenosine spacers between monomers would increase the sensitivity for detecting host protein binding to NSP3. However, our results indicated that a monomer was equally effective for pulling down NSP3 (data not shown). Therefore, for the rest of the study, we used the monomeric probes to reduce the possibility of creating artificial secondary structures.
By comparing a scrambled sequence (AUAGGCGUC) to an authentic conserved RV 3ЈUTR probe (GAUGUGACC) in the presence or absence of RV infection, we were able to identify host proteins that specifically immunoprecipitated with the 3ЈUTR consensus sequences (Table S1). We then used CRAPome filtering analysis (39) and a Significance Analysis of INTeractome (SAINT) threshold score of 0.9 to identify specific host proteins. Our analysis identified three high-confidence "hits": the mitochondrial F1 complex ␤-polypeptide ATP5B; the RAS oncogene family member RAB1A; and the isoleucyl-tRNA synthetase IARS as interacting partners to the RV 3ЈUTR consensus sequences ( Fig. 2E and Table S2). Of note, ATP5B, an ATP synthase ␤ subunit (29), was found to be enriched at 100-fold in its interaction with 3ЈUTR probe over control (Fig. 2E). ATP5B was previously reported to promote human immunodeficiency virus

ATP5B supports rotavirus infection
(HIV)-1 replication and chikungunya virus infection (40,41), prompting us to hypothesize that it may also act as a pro-RV host factor. Importantly, we also detected two other pivotal components of the ATP synthase complex, ATP5A1 and ATP5O, in the pulldown of 3ЈUTR sequence during RV infection (Fig. 2F). It is worth noting that all three ATP synthase subunits were found to interact with the RV 3ЈUTR probe only in the context of RV infection (Fig. 2F). Taken together, this new assay technique to detect host RNA-binding proteins allowed us to identify several novel proteins, including several ATP synthase subunits that had previously not been known to interact with viral RNA.

ATP5B immunoprecipitates and co-localizes with RV 3UTR during infection
We next carried out a set of experiments to validate the association between ATP5B and the conserved RV 3ЈUTR sequences. Using streptavidin beads to directly pull down biotinylated host proteins in mock-or RV-infected cells, we observed strong bind-ing of endogenous ATP5B with an RV 3ЈUTR probe, but only very weak interaction was detected with the scrambled probe (Fig. 3A). Importantly, the ATP5B-3ЈUTR interaction was observed exclusively during active RV infection (Fig. 3A), suggesting that such interaction might be regulated by other host factors, viral factors, or virus-induced cellular changes.
We hypothesize that to physically interact with RV RNA, ATP5B is likely to co-localize with viral RNA in infected cells. We first determined ATP5B localization in the presence or absence of RV infection. In mock-infected cells, almost all ATP5B was found at the mitochondria (Fig. 3B, upper panel), consistent with previous reports (42) and its role as part of the ATP synthase complex. However, in RV-infected cells, we observed a significant condensation of mitochondrial organization in viral antigen VP6positive cells (Fig. 3B, lower panel), recently noted by Green and Pelkmans (1) as well. To examine the subcellular localization of ATP5B relative to viral RNA, we designed specific fluorescent in situ hybridization (FISH) probes that detect RV RNA with high resolution. Our RV FISH probes were highly specific for viral RNA Figure 2. RaPID identifies host proteins that bind to group A RV 3UTR. A, schematics of a two-plasmid system to screen for RV RNA-binding host proteins. BirA fusion protein expression is driven by a cytomegalovirus promoter, and RNA expression is driven by a U6 promoter. B, HEK293T cells were co-transfected with a BirA expression plasmid and a pMotif-scrambled or EDEN15 plasmid for 24 h and incubated with biotin (50 M) for 24 h. Immunoprecipitation was performed using magnetic streptavidin beads. Eluted lysates and input were examined by Western blotting using primary antibody against CUG-BP1. C, same experiment as B except that HEK293T cells were transfected with a GFP-NSP3 expression plasmid together with a pMotif plasmid encoding scrambled sequences or conserved RV 3ЈUTR. D, HEK293T cells were infected with RRV (m.o.i. ϭ 3), transfected with pMotif-scrambled or 3ЈUTR plasmids at 1 hpi, and cultured in the presence of biotin at 24 hpi for another 24 h. Immunoprecipitation and Western blotting were performed to measure the levels of endogenous NSP3. Quantification was performed using ImageJ on the basis of three blots. E, screening results for RV 3ЈUTR binding host proteins. The x axis represents the fold change of specific host protein binding to 3ЈUTR sequences as compared with a scrambled sequence during RV infection. The y axis represents the SAINT score (ranging from 0 to 1.0), calculated based on the reproducibility between three independent immunoprecipitation-MS experiments. F, spectral counts of several ATP synthase complex subunits from a representative MS and bioinformatics analysis. (Scr, scrambled). For all figures, experiments were repeated at least three times with similar results. Experiment in E was performed three times independently.

ATP5B supports rotavirus infection
and only observed within the viroplasm, the distinct punctate viral factories that are positive for the nonstructural protein NSP2 (Fig.  3C). Importantly, viral RNA signal was also seen, at least partially, to overlap with ATP5B localization during infection ( Fig. 3D and inset). In contrast, minimal co-localization was observed between RV RNA with mitochondrial antiviral-signaling protein (MAVS) (Fig. 3D). Therefore, data from both immunoprecipitation and immunofluorescence experiments are consistent with our screening results indicating that ATP5B is a bona fide RV 3ЈUTRinteracting protein.

ATP5B supports group A RV infection in a strain-independent manner
Having established ATP5B-rotaviral RNA interaction, we next sought to examine the functional role of ATP5B during RV infection. Because genetic depletion of ATP5B induces lethality (43), we utilized siRNA to knock down ATP5B expression levels. All three top candidates identified in our proteomics screen, ATP5B, RAB1A, and IARS, were effectively silenced with their specific siRNAs and examined by their mRNA levels by quantitative PCR (Fig. 4A). Both ATP5B and IARS siRNA also led to

ATP5B supports rotavirus infection
a marked reduction of respective protein levels, although the decrease was more significant with IARS than ATP5B (Fig. 4B). Importantly, in a multiple-round infection assay, at 24 h postinfection, knockdown of ATP5B led to 40% reduction in the mRNA levels of NSP5, an RV gene transcript representative of intracellular RV RNA synthesis (Fig. 4C), and we also observed an ϳ40% decrease in the amount of infectious RVs in the supernatant of the ATP5B knockdown cells (Fig. 4D), suggesting that ATP5B facilitates RV replication in host cells.
Given that all group A RVs share consensus sequences of the ATP5B-binding site at their 3ЈUTR ends (Fig. 1), we next examined whether the pro-RV role of ATP5B is conserved across different RV strains. Consistent with our hypothesis, the infectivity of a panel of cell culture adapted group A human and animal RVs was inhibited by siRNA depletion of ATP5B (Fig.  4E). In contrast, the replication of several other RNA viruses, including vesicular stomatitis virus and two strains of influenza A viruses, which do not possess 3Ј-UGUGACC sequences, was

ATP5B supports rotavirus infection
not affected by ATP5B knockdown (Fig. 4, F and G). ATP5B knockdown also did not influence the replication of Salmonella typhimurium, an intracellular bacterial pathogen or even reovirus, another Reoviridae family member with many properties similar to RVs (Fig. 4, H and I), thereby in support of a highlyspecific pro-RV role of ATP5B.

ATP5B contributes to RV genome assembly
To mechanistically determine at which step in the RV replication cycle in which ATP5B is involved, we undertook a series of experiments to target various molecular processes within a single RV replication cycle. First, initial viral adsorption and subsequent endocytosis, assayed by input viral RNA levels 4 and 37°C post-incubation, respectively, and using EDTA wash as a negative control, were not affected by ATP5B siRNA, when compared with control siRNA (Fig. 5, A and B). These two pieces of experimental evidence strongly suggest that early events in RV infection are not perturbed by the loss of ATP5B.
Next, we measured viral RNA by quantitative PCR and protein levels using polyclonal antibody that recognizes RV double-layered particles (44) at 8 h post-infection within a single virus replication cycle. We found that both levels were comparable between control and ATP5B siRNA-transfected cells (Fig.  5, C and D). Thus, we focused our further attention on the formation of viroplasms, the specialized electron-dense inclusions within the cytoplasm where assembly of immature virus particles is initiated (45). RV viroplasms, detected by a mAb against viral protein NSP2 (46), were also similar with or without ATP5B siRNA silencing (Fig. 5E). Importantly, at 8 h postinfection, we found significant co-localization of ATP5B with NSP2, further supporting our previous immunofluorescence analysis on ATP5B and viral RNA (Fig. 3D). ATP5B was found to co-localize only with NSP2 but not with the two other RV proteins VP4 and VP6 (Fig. 5F). This finding was further validated using co-staining of ATP5B, NSP2, and RV RNA in the same infected cells (Fig. 5G). Using a conformation-dependent mAb that specifically recognizes the trimeric form of VP7 on mature RV particles (47), we correlated and estimated the number of fully assembled RV particles. At 8 h post-infection, there was a significant 4-fold decrease in ATP5B-depleted cells (Fig.  5H). Consistent with a role of ATP5B in the later steps in the RV replication cycle, the amount of infectious RVs and virus progeny released into the supernatant was also reduced by ATP5B knockdown (Fig. 5, I and J). Taken together, these data point to a supportive role of ATP5B for RV genome assembly, which is likely to be a high-energy-consuming process and occurs following viroplasm formation and prior to virus egress.

Inhibition of ATP synthase complex restricts RV infection
Using available crystal structures of both the ATP synthase complex (21) and the RV polymerase VP1 binding to UGU-GACC (48), we extracted the viral RNA structure and attempted to computationally model the interaction between ATP5B and RV 3ЈUTR sequences. We applied HDOCK (49), an on-line bioinformatics tool that combinatorially adopts template-based modeling and free docking, to simulate the biochemical association between ATP5B and the RV 3ЈUTR. This algorithm calculates the possibility and potential affinity of a protein-RNA ligand interaction and assigns a root-meansquare deviation (RMSD) score. The lower the score, the more likely this interaction represents a biologically meaningful binding (50). Of note, compared with the negative control RV protein NSP4, which localizes to the endoplasmic reticulum and does not bind viral RNA (51,52), the validated UGU-GACC-binding protein NSP3 had a much lower RMSD value (Fig. 6A). Importantly, the ATP5B-3ЈUTR value was predicted to be similar to our positive control NSP3 value and much lower than cullin-3, a host protein that interacts with RV NSP1 but not viral RNA that we used as a negative control (Fig. 6A) (3). In addition, we found that the inclusion of ATP5A1, another component of ATP synthase complex and also identified by our proteomics screen (Fig. 2F), to ATP5B significantly favored the interaction with RV 3ЈUTR by 24% (Fig. 6A). It is noteworthy that the ATP synthase-viral RNA interaction is predicted to be weaker than that between viral polymerase VP1 and RV RNA, suggesting that ATP5B-RNA binding may not be direct and possibly require additional host/viral factors. We further predicted the spatial position of UGUGACC, which inserts into the canal of the headgroup of F 0 -F 1 -ATP synthase complex, consisting of three ATP5A1 subunits and three ATP5B subunits (Fig. 6B). Such conformation would allow the Arg-286 of ATP5A1 and Pro-276 of ATP5B to create a docking site for the negatively charged RV 3ЈUTR sequences (Fig. 6C).
Consistent with the hypothesis that the full structure of the ATP motor is needed for RV infection, we found that siRNA knockdown of both ATP5A1 and ATP5O phenocopied ATP5B silencing and suppressed RV infection (Fig. 7A), suggesting that the pro-RV role of ATP5B depends on the integrity of the holoenzyme. In line with this observation, the concurrent knockdown of several ATP synthase complex subunits led to an even more prominent reduction in RV infection (Fig. 7A). In addition to siRNA-mediated genetic depletion, a specific smallmolecule ATP synthase inhibitor, isoapoptolidin (53), potently down-regulated intracellular RV RNA levels and virus yield in the supernatant (Fig. 7, B and C). Several other pharmacological inhibitors of ATP synthase activity also inhibited RV replication in a similar manner (Fig. 7D). These inhibitors did not induce cytotoxicity, just like control and ATP5B siRNA (Fig.  7E). Further corroborated by the fact that the replication of other pathogens was not affected by ATP5B knockdown (Fig. 4, F-I), our data demonstrate that reduction in RV infection was not accounted for by the unhealthy cellular environment created by ATP5B inhibition. Collectively, our results support an important function of the ATP synthase complex in promoting efficient RV replication by enhancing viral genome assembly.

ATP synthase supports RV infection in human enteroids
Finally, we extended our findings and tested the role of ATP synthase activity in RV replication in human intestinal enteroids, a primary human intestinal epithelial cell (IEC) three-dimensional culture system that faithfully recapitulates the IEC population diversity in vivo (54,55). We found that isoapoptolidin treatment of human intestinal enteroids significantly decreased plaque-forming units while not affecting RV RNA and protein levels in a one-step replication cycle (Fig. 8). The intracellular NSP5 levels and percentage of viral antigen-

ATP5B supports rotavirus infection
positive cells were comparable in mock and ATP inhibitor treatment (Fig. 8, A and B, respectively). However, blocking ATP synthase activity significantly reduced the formation of progeny infectious RV focus-forming units (Fig. 8C), further supporting a role for ATP5B in RV assembly.

Discussion
Mechanisms of genome replication and particle assembly of segmented icosahedral RNA viruses are poorly understood. In the case of RVs, a prototypical enteric viral pathogen with 11 discrete dsRNA segments, the majority of what we know derives from studies on viral replication intermediates and cis-acting signals on viral RNAs. The current model states that all 11 RV (ϩ) RNA segments are present in the viroplasm and that the viral VP1 polymerase is recruited and interacts with a highly structured region toward the end of a conserved 3ЈUTR UGUGACC, which forms a panhandle structure (31, 48, 56 -60). The additional binding of VP3 with the 5ЈUTR end facilitates the formation of an RNA(ϩ)-VP1-VP3 complex (61). Subsequently, different segments are assorted via a not well-characterized RNA-RNA interaction and packaged in a VP2-dependent manner (18,62,63). However, whether any and which host factors might participate in these processes remains largely unexplored.
In this study, we utilized a recently developed method (RaPID) (28) to identify RNA-binding proteins to study RV-host interaction. Compared with the traditional tools to study RNAprotein interaction (64 -66), this new method has several advantages including the following: 1) using whole-infected cells rather than lysates to enable RNA-protein binding to take place in its "natural" cellular environment; 2) bypassing the need for crosslinking and reversal; and 3) enhanced sensitivity based on biotin protein ligase BirA that promiscuously biotinylates proteins in a proximity-dependent fashion. Compared with other methods of analyzing essential host factors of virus infection, RaPID offers a different and complementary approach, as ATP5B knockout is lethal according to the recent haploid screens for survival genes (43), which probably explains why ATP5B was not identified as a major hit in our recent CRISPR-Cas9 screen for pro-RV host factors (4). Of note, RaPID technique can be easily adopted to examine the RV gene 11 loop regions and helical junctions (19) or a 6-nucleotide-long pyrimidine-rich tract shared by several RV genes (67). This method can also be used to study the structural elements of other RNA viruses, such as the stem-loop regions and subgenomic viral RNA sequences.
In addition to ATP5B, our screen also identified IARS and RAB1A as rotaviral RNA-binding proteins. It is intriguing that these proteins exhibit distinct subcellular localization: ATP5B at mitochondria where the key adaptor protein for cytosolic RNA-sensing pathway MAVS is localized (68); RAB1A at endoplasmic reticulum; and IARS in the cytosol. Because we did not observe significant effects on RV replication from the siRNA knockdown of ATP5B in a single virus replication cycle, it is possible that their interactions with viral RNA do not have direct effect on virus replication but rather serve as modulators of innate immune signaling in the cytoplasm. Interestingly, upon examination of whether the interferon pathway is potentiated or dampened in the absence of ATP5B during RV infec-

ATP5B supports rotavirus infection
tion, we found that both IFN-3 (IFNL3) expression and secretion into the supernatant were significantly reduced in RRVinfected ATP5B knockdown cells as compared with the control siRNA-transfected cells (Fig. 9). We will explore this new direction further in future studies.
An important aspect of this study that remains to be further examined is whether ATP5B directly or indirectly binds to RV 3ЈUTR during infection. Our bioinformatics modeling approach suggests that the ATP5A1-ATP5B head complex could potentially interact with the conserved UGUGACC sequence at a strength comparable with NSP3. However, considering the even higher affinity between the RV polymerase VP1 and 3ЈUTR (Fig. 6A), one might expect that ATP5B-3ЈUTR binding is also dependent on the presence of selected RV proteins, such as VP1. This hypothesis is consistent with our data that the interaction between 3ЈUTR and all three subunits of ATP synthase complex (5A1, 5B, and 5O) is greatly enhanced during RV infection (Fig. 2F) and that several

ATP5B supports rotavirus infection
attempts to directly test ATP5B-RV RNA interaction using electrophoretic mobility shift assay have not been successful (data not shown). It is tempting to speculate that the on and off status of ATP5B may provide critical energy when the ensemble of VP1-VP3-RNA and VP2 building blocks comes together to encapsulate the RNA segments into the newly-forming capsid.

Cells and plasmids
HEK293 (CRL-1573) and HEK293T (CRL-3216) cells were purchased from American Type Culture Collection (ATCC) and cultured in complete DMEM. MA104 (CRL-2378.1) cells were obtained from ATCC and cultured in complete M199 medium. All three cell lines are not listed as commonly misidentified cell lines maintained by International Cell Line Authentication Committee (ICLAC), and they tested negative for mycoplasma contamination.
The pMotif vector that drives the transcription of target RNA (scrambled sequences, EDEN15 motifx5, and RV-3ЈUTR) was digested with BsmBI and ligated as described previously (28). For each Western blotting experiment, 2 g of pMotif and 0.3 g of pRaPID with optimized N sequence were co-transfected into one well of a 6-well plate pre-seeded with HEK293T cells using 5 l of Lipofectamine 2000 (Invitrogen). At 24 h post-transfection, biotin (Sigma, catalog no. B4501, 50 M) was added to media and incubated for 18 h. For each MS experiment, a single 10-cm dish was used instead of a 6-well plate with 12 g of pMotif RNA and 2 g of pRaPID transfected with 30 l of Lipofectamine 2000. The T7 vector that encodes gene 11 from simian RV SA11 strain was purchased from Addgene (catalog no. 89172).

RaPID
On the day of harvest, cells were washed with cold PBS (one time) and lysed with 200 l of lysis buffer at room temperature. 20 l of 25% Triton X-100 was then added and mixed thoroughly. Subsequently, lysates were supplemented with 220 l of cold 50 mM Tris, pH 7.4, sonicated for 10 s, and centrifuged at 4°C for 10 min. Free biotin was removed by Macrosep Advance Spin Filter 3K MWCO 20 ml (VWR 89131-974). Protein concentration in each sample was determined using spectrophotometer with the Pierce protein quantitation assay (Thermo-Fisher Scientific). Protein concentration across samples was normalized, and 40 l of the sample was aliquoted and marked as IP input. Pulldown of biotinylated proteins was performed with magnetic beads MyOne TM Streptavidin C1 (Life Technol-

ATP5B supports rotavirus infection
ogies, Inc., catalog no. 65002) and sample was rotated with beads for 2 h at room temperature. Samples were washed three times, and beads were boiled in Elution Buffer for 45 min at 98°C and labeled as IP. Both input and IP were loaded and resolved on pre-cast Bio-Rad gels, and SDS-PAGE was conducted as described previously (69) using the following primary antibodies: ATP5B (HPA001520, Sigma), CUG-BP1 (3B1, catalog no. sc-20003, Santa Cruz Biotechnology), GAPDH (catalog no. 631402, BioLegend), and NSP3 (kindly provided by Drs. Susana Lopez and Didier Poncet). Secondary incubation was performed with anti-rabbit (CST, catalog no. 7074) or antimouse (CST, catalog no. 7076) IgG horseradish peroxidaselinked antibodies. Protein bands were visualized with Clarity ECL substrate (Bio-Rad, catalog no. 170-5061), Amersham Biosciences Hyperfilm (GE Healthcare), and STRUCTURIX X-ray film processor (GE Healthcare).
For sample preparation for LC-MS/MS, streptavidin magnetic beads were washed with 500 l of 50 mM ammonium bicarbonate three times. Beads were resuspended in 200 l of 50 mM ammonium bicarbonate supplemented with DTT to a final concentration of 5 mM, incubated on a heat block at 50°C for 5 min followed by head-over-head rocking for 30 min at room temperature. Alkylation was performed by the addition of propionamide to a final concentration of 10 mM and head-overhead rocking for 30 min at room temperature. 250 ng of trypsin/LysC (Promega) was added to each sample and digested overnight at room temperature in the head-over-head shaker followed by the addition of formic acid to 1%. Peptides were removed and washed with 50 l of 0.1% formic acid water. The acidified peptide pools were purified with C18 STAGE tip 37 (NEST group) microspin columns and dried in a speed vac. For LC-MS/MS, peptide pools were reconstituted and injected onto a C18 reversed phase analytical column. All MS/MS data were first analyzed in Preview to provide recalibration criteria and then reformatted to .MGF before full analysis with Byonic version 1.4 (ProteinMetrics). Analyses used Uniprot canonical .fasta files for human, concatenated with common contaminant proteins. Data were searched at 10 ppm mass tolerances for precursors, with 0.4-Da fragment mass tolerances assuming up to two missed cleavages and allowing for N-ragged tryptic digestion, and were validated at a 1% false discovery rate using typical reverse-decoy techniques. The resulting identified peptide spectral matches and assigned proteins were then exported for further analysis using custom tools developed in MatLab (MathWorks) to provide visualization and statistical characterization. SAINT scores were generated with spectral counts from experimental samples and controls using the on-line CRAPome tool (39) to minimize the MS background.

Virus infections
All human and animal RV strains used in this study were propagated in MA104 cells, and RV infection was performed as described previously (3). Recombinant VSV (strain Indiana) expressing GFP was a kind gift from Dr. Jack Rose (Yale University) and used as described previously (70). Influenza A viruses (H1N1 A/California/7/2009, H3N2 A/Victoria/261/ 2011) were used as described previously (71). Reovirus T1L strain was kindly provided by Dr. Carlos Arias (Universidad Nacional Autónoma de México). For all virus infection experiments, the initial inoculum was removed after 1 h of incubation at 37°C and washed with complete or serum-free medium to avoid unsynchronized infection. For virus adsorption assay, cells were incubated with RVs at 4°C for 1 h and washed with ice-cold PBS three times. For virus endocytosis assay, cells were infected with RVs at 37°C for 1 h and washed with serum-free medium containing 5 mM EDTA to remove cell-bound virus at the surface. For multiple rounds of RV infection, 0.05 g/ml trypsin was added to the serum-free medium for outer capsid protein VP4 cleavage, and infection lasted for at least 24 h.

Quantitative RT-PCR
Total RNA was harvested in RLT buffer supplemented with ␤-mercaptoethanol and extracted as described previously (72). Random hexamer was used for reverse transcription reaction. qPCR was performed with the Stratagene Mx3005P (Agilent) with a 25-l reaction consisting of 50 ng of cDNA, 12.5 l of Power SYBR Green master mix (Applied Biosystems), and 200 nM both forward and reverse primers. All SYBR Green primers (Table S3) have been validated with both dissociation curves and DNA electrophoresis of the correct amplicon size.

siRNA transfection
HEK293 cells in 24-well plates were transfected with 6 pmol of control or specific siRNAs against host genes (see information in Table S3) with RNAiMAX (Invitrogen) using a reverse transfection method (3). Cells were cultured for at least 48 h prior to qPCR analysis or virus infections. The transfection was performed in antibiotic-free medium and cells cultured for 48 h prior to virus infection, qPCR or Western blot analysis.

Focus-forming unit assay
HEK293 cells were transfected with control or ATP5B siRNA and infected with RRV at m.o.i. ϭ 0.5 in 24-well plates. At different time points post-infection, supernatants were collected and applied to freshly seeded MA104 cells in a 10ϫ dilution series as described previously (71). In brief, infected cells were fixed overnight and incubated with primary antibodies against rabbit anti-DLP at 37°C for 1 h, and secondary incubation was performed with anti-rabbit IgG at 37°C for 1 h, with the subsequent addition of AEC substrate (SK-4200, Vector Laboratories). Color developing was observed, and the number of fluo-ATP5B supports rotavirus infection rescent colonies was counted. Two independent experiments were performed with triplicate infections.

Salmonella invasion assay
Salmonella enterica serovar Typhimurium strain SL1344 was grown on an LB agar plate at 37°C. A single colony was used to inoculate an LB broth culture and grown overnight with shaking at 37°C. Bacteria were subcultured 1:20 into high-salt LB (300 mM NaCl) and grown statically at 37°C in a 5% CO 2 incubator for 2-3 h. Bacteria were pelleted and resuspended in DMEM, 10% FBS at an m.o.i. of 10 and added to wells of HEK293 cells previously transfected with ATP5B or control siRNA. After 30 min of infection, bacteria were removed, and DMEM, 10% FBS with 100 g/ml gentamicin were added to kill all extracellular bacteria. After 90 min, media were changed to 10 g/ml. Intracellular bacteria were quantified at 2 or 8 h postinfection. Infected cells were lysed with 1% Triton X-100 in PBS for 5 min at room temperature, serially diluted, and then plated on LB agar plates to determine the number of colony-forming units.

Cytotoxicity assay
Cell survival rate of ATP synthase inhibitor-treated or siRNA-transfected HEK293 cells were measured by a luminescence-based CellTiter-Glo 2.0 Assay (Promega, G9241) according to the manufacturer's instructions. In brief, ATP was serially diluted from 100 M to 0.01 nM to formulate a standard curve. Cell lysates were harvested, and intracellular ATP levels were measured and normalized to the standard curve.

RNA fluorescent in situ hybridization
A custom set of 48 Stellaris FISH probes to detect the RRV genome is directly conjugated to Quasar 670 dye (LGC Biosearch Technologies). Mock-or RV-infected cells were fixed, permeabilized, and hybridized with FISH probes to visualize viral RNA molecules according to the manufacturer's instructions. Some samples were also stained with anti-ATP5B and anti-NSP2 antibodies to measure co-localization.

Statistical analysis
The bar graphs are displayed as mean Ϯ S.D. Statistical significance in Figs. 4, A and E-I, 5, C and H-J, 7, B and C, and 8, A and C, was calculated by Student's t test using Prism 7.0c (GraphPad). Statistical significance in Figs. 3D, 4, C and D, 5, A,  B, and F, 7, A, D, and E, and 9, A and B, was calculated by pairwise analysis of variance using Prism 7.0. p values are indicated on each figure (*, p Յ 0.05; **, p Յ 0.01; ***, p Յ 0.001). All experiments, unless otherwise noted, have been repeated at least three times. Mass spectrometry experiments in Fig. 2E were performed precisely three times.