Ribosome Biogenesis Modulates Ty1 Copy Number Control in Saccharomyces cerevisiae.

Transposons can impact the host genome by altering gene expression and participating in chromosome rearrangements. Therefore, organisms evolved different ways to minimize the level of transposition. In Saccharomyces cerevisiae and its close relative S. paradoxus, Ty1 copy number control (CNC) is mediated by the self-encoded restriction factor p22, which is derived from the GAG capsid gene and inhibits virus-like particle (VLP) assembly and function. Based on secondary screens of Ty1 cofactors, we identified LOC1, a RNA localization/ribosome biogenesis gene that affects Ty1 mobility predominantly in strains harboring Ty1 elements. Ribosomal protein mutants rps0bΔ and rpl7aΔ displayed similar CNC-specific phenotypes as loc1Δ, suggesting that ribosome biogenesis is critical for CNC. The level of Ty1 mRNA and Ty1 internal (Ty1i) transcripts encoding p22 was altered in these mutants, and displayed a trend where the level of Ty1i RNA increased relative to full-length Ty1 mRNA. The level of p22 increased in these mutants, and the half-life of p22 also increased in a loc1Δ mutant. Transcriptomic analyses revealed small changes in the level of Ty1 transcripts or efficiency of translation initiation in a loc1Δ mutant. Importantly, a loc1Δ mutant had defects in assembly of Gag complexes and packaging Ty1 RNA. Our results indicate that defective ribosome biogenesis enhances CNC by increasing the level of p22, and raise the possibility for versatile links between VLP assembly, its cytoplasmic environment, and a novel stress response.

Ty1 is the most abundant retrotransposon in the Saccharomyces cerevisiae reference strain at 32 copies per haploid genome, and is an effective model for understanding the mechanism and consequences of retrotransposition (Curcio ). Ty1 shares many similarities with retroviruses, except that retrotransposition is not infectious. Ty1 long terminal repeats (LTRs) bracket two overlapping genes: GAG encodes the capsid protein of virus-like particles (VLPs), and POL encodes a polyprotein that contains protease (PR), integrase (IN), and reverse transcriptase (RT). Once transcribed by RNA polymerase II from the 5′ LTR to the 3′ LTR, Ty1 mRNA is exported to the cytoplasm and translated into two protein products, Gag-p49 and Gag-Pol-p199. Gag-Pol is generated by a programmed +1 frameshifting event near the end of GAG. Gag-p49 and Gag-Pol-p199 are present in a ratio of 20:1 (Kawakami ). Ty1 mRNA and translated proteins form cytoplasmic granules called the retrosomes, or T-bodies, which are sites where VLPs assemble (Malagon and Jensen 2008; Checkley ; Dutko ). During or after Ty1 VLP assembly, Gag-p49 and Gag-Pol-p199 undergo maturation catalyzed by PR to form mature Gag-p45, PR, IN, and RT. Reverse transcription takes place using packaged dimeric Ty1 RNA as template, and a heterodimeric complex of IN and RT (Chapman ; Wilhelm ; Feng ). Ty1 cDNA and IN likely form a preintegration complex or intasome, which is imported into the nucleus via a nuclear localization sequence present on IN (Kenna ; Moore ; McLane ). Genomic regions upstream of genes transcribed by RNA polymerase III (such as tRNA genes) are favored sites for Ty1 integration, based on functional interactions between IN and subunits of RNA polymerase III (Devine and Boeke 1996; Baller ; Bridier-Nahmias ; Cheung ).

While retrotransposons and retrovirus-like elements comprise almost 50% of the human genome, Ty1, Ty2, Ty3, Ty4, and Ty5 sequences make up ∼3% of the compact S. cerevisiae reference genome (Kim ; Lander ). Ty1 movement occurs at a low rate of ∼1 × 10−6 transposition events/Ty1 element/generation (Curcio and Garfinkel 1991), yet yeast has the capacity to harbor hundreds of additional Ty1 insertions when cells containing Ty1 fused to the promoter undergo multiple rounds of induction (Scheifele ). These observations suggest that cellular or element-encoded products prevent rampant Ty1 retrotransposition in S. cerevisiae. Interestingly, none of the silencing pathways used by other eukaryotes to minimize transposon movement, such as RNAi or DNA methylation, are present in the S. cerevisiae sensu stricto group (Drinnenberg , 2011; Levin and Moran 2011), implying there may be a novel pathway utilized by S. cerevisiae to modulate Ty1 retrotransposition. Copy-number dependent effects on Ty1 mobility, termed copy number control (CNC) (Garfinkel ) were first addressed when a plasmid containing Ty1 tagged with the retromobility indicator was expressed in a strain of S. paradoxus strain that lacks complete Ty1 elements (Curcio and Garfinkel 1991; Moore ). A dramatically higher level of Ty1 mobility is observed in this Ty1-less strain when compared to the S. cerevisiae laboratory strain S288c, which carries 32 Ty1 elements, or when the Ty1-less strain is repopulated with Ty1. Interestingly, the minimal Ty1 sequence required for CNC is located within GAG and the 5′ LTR. In cells undergoing CNC, defective VLPs are formed that contain lower levels of mature RT and IN, and reverse transcription products (Matsuda and Garfinkel 2009; Purzycka ).

Although early work suggested that Ty1 antisense transcripts confer CNC (Matsuda and Garfinkel 2009), extensive biochemical and genetic analyses of the CNC region revealed a self-encoded Ty1 protein, termed p22, that is necessary and sufficient for CNC (Saha ). p22 is encoded by the 3′-half of GAG, and is produced from subgenomic, internally initiated, Ty1(i) transcripts. p22 is processed to p18 by PR cleavage at the same site utilized by the full-length Gag protein. p22/p18 confers CNC and a trans-dominant negative effect on Ty1 mobility. p22/p18 associates with VLPs, co-immunoprecipitates with Ty1 Gag, and colocalizes with Gag in the cytoplasm. p22/p18 disturbs retrosome formation and VLP assembly, and blocks maturation and reverse transcription within the VLPs. In addition, p18 interferes with the nucleic acid chaperone function of Gag-p45 (Nishida ), and p22/p18 targets Gag and disrupts Gag–Gag interactions during VLP assembly (Tucker ).

Like retroviruses, retrotransposons are capable of making only a few proteins critical for replication; therefore, a variety of host functions modulate their life cycle (Maxwell and Curcio 2007; Goff 2008). Host genes from diverse biological processes affect Ty1 transposition and can stimulate (cofactors) or inhibit (restriction factors) various aspects of the Ty1 life cycle. To date, ∼500 cofactor and 110 restriction factor genes have been identified by genome-wide screens using yeast gene deletion mutants (Griffith ; Nyswaner ; Dakshinamurthy ; Risler ) or transposon mutagenesis (Scholes ). Ty1 cofactors are involved in processes such as transcription, chromatin modification, histone deacetylation, ribosome biogenesis and function, and mRNA turnover and transport (Griffith ; Dakshinamurthy ; Risler ). In contrast, ∼80% of Ty1 restriction factors are involved in nuclear processes including DNA repair, transcription, chromatin structure and function, recombination, and the cell cycle (Scholes ; Nyswaner ; Risler ). We set out to identify host modulators that influence the synthesis or activity of the newly discovered Ty1 restriction factor p22. Such host factors would greatly help us understand the mechanism of CNC.

We reasoned that CNC-specific host factors should stimulate or inhibit CNC in a Ty1 copy number-dependent manner. In other words, a host modulator displaying a CNC-specific phenotype would occur predominantly in CNC+ cells containing Ty1 elements, but not in CNC− cells lacking Ty1. To identify CNC-specific modulators, we analyzed previously identified Ty1 cofactor mutants involved in RNA metabolism, which is central to the process of Ty1 retrotransposition (Curcio ), for those that displayed phenotypes strongly associated with CNC (Matsuda and Garfinkel 2009; Saha ). We then analyzed candidate deletion mutants in congenic CNC− and CNC+ strains that differ only in Ty1 copy number. Surprisingly, we identified several genes involved in ribosome biogenesis that increase the level p22 and show CNC-specificity; Ty1 mobility decreased preferentially in CNC+ cells when compared to CNC− cells. In particular, , a gene required for asymmetric localization of mRNA (Heym and Niessing 2012) and rRNA processing (Urbinati ), facilitates Gag–Gag interactions and retrosome formation by modulating the level of p22.

Materials and Methods

Genetic techniques, media, strains, and plasmids

Strains are listed in Supplemental Material, Table S1. Standard yeast genetic and microbiological procedures were used in this work (Guthrie and Fink 1991). The S. cerevisiae CNC− strain DG3453 was derived from DJ12, which was isolated in Djibouti, Africa, and generously provided by Michel Aigle. A MATα ascospore (DJ12-1B) was subjected to 5-fluoroorotic acid (5-FOA) selection (Boeke ), and a spontaneous mutant was isolated to generate DG3093. The results of Southern blot analysis (Garfinkel ; Moore ) suggest there are 10 Ty1 elements in DJ12-1B and DG3093. To eliminate the Ty1 elements from DG3093, we generated a Ty1: targeting fragment containing Ty1 sequences 1702–3945 (Boeke ). The 2243-bp segment was amplified by PCR from DG3093 using primers 1702XhoBglF (5′-CCCGCTGCAGAGATCTGCTCATCACATACACTC-3′) and 3945BamR (5′-CCGCGGATCCTGCAATCAGGTGAATTCGT-3′) and subcloned into pSP70 (Promega, Madison, WI) digested with XhoI and BglII. bracketed by XhoI restriction sites on pB130 (kindly provided by Gerald Fink) was digested and cloned into the Ty1 SalI site (nucleotide 2173) of pSP70/Ty1. The resulting plasmid pBDG1374 was digested with XhoI and EcoRI to release the Ty1: targeting fragment. Ura+ transformants containing putative Ty1: recombination events with individual Ty1 elements were subjected to 5-FOA selection followed by Southern blot analysis to detect LTR-LTR recombinants (Winston ). Changes in Ty1 hybridization patterns suggest that one Ty1 was eliminated with each round of transformation and 5-FOA selection. After removing nine of 10 Ty1 elements, was deleted as described previously using pBDG652 (Garfinkel ) to generate DG3453. The same Ty1 was present in independent elimination lineages, and could not be lost to generate a Ty1-less strain, raising the possibility that the remaining Ty1 insertion may be required for the function of an essential gene. Crosses between DG3093 containing 10 Ty1 elements and strains with a single Ty1, or crosses between the reference strain BY4741 (Brachmann ) and the single Ty1 strains displayed 70% spore viability in 18–20 tetrads/cross. Ascosporal progeny displayed similar variations in growth and expected segregation of genetic markers. DG3453 was repopulated with ∼19 Ty1 elements by induction of a pGTy1 plasmid to generate DG3648, as described previously (Garfinkel ). MATα deletion strains (Giaever ) derived from BY4742 (Brachmann ) were obtained from Invitrogen (Carlsbad, CA). Gene disruptions in DG3453 and DG3648 were carried out using the KanMX4 cassette (Wach ). The KanMX disruption cassettes were PCR amplified using template DNA extracted from deletion mutants with gene-specific primers A and D from the Saccharomyces Genome Database (https://www.yeastgenome.org/). Correct gene disruptions by KanMX for all of the deletion mutant strains used in this study were confirmed by PCR using locus-specific primers. To analyze Ty1 mobility and protein levels, pBDG606 (pGTy1/Cen/) (Dakshinamurthy ) was introduced to BY4742-, DG3453-, and DG3648-derived strains.

Ty1his3-AI mobility

Quantitative and qualitative Ty1 mobility assays were determined as described previously (Curcio and Garfinkel 1991; Garfinkel ; Sharon ) with minor modifications. For strains expressing pGTy1 (pBDG606), a single colony grown on SC-Ura at 30° was resuspended in 1 ml of SC-Ura + 2% raffinose, and incubated for 24 hr at 30° with aeration; 200 μl of the raffinose culture was centrifuged, and the cell pellet was resuspended in 1 ml of SC-Ura + 2% galactose in quadruplicate. The cultures were grown for 16 hr at 22°, washed, diluted, and spread onto SC-Ura and SC-His-Ura plates, which were incubated for 3–4 days at 30° until colonies formed. The frequency of Ty1 mobility was calculated by dividing the number of His+ Ura+ colonies by the number of Ura+ colonies. For qualitative Ty1 mobility assays, cells patched onto SC-Ura were incubated at 30° for 2 days. Cells were replica plated onto SC-Ura + 2% galactose medium plates and incubated at 22° for 16 hr. To detect Ty1 mobility events, galactose-induced cells were replica plated to SC-Ura-His followed by incubation for 2–3 days at 30°.

RNA isolation and northern blot analysis

For strains carrying pBDG606, 10 ml of SC-Ura + 2% raffinose was inoculated with a single colony, and the culture was grown at 30° to an OD600 of 2, or up to 24 hr, depending on the strain; 2 ml of the raffinose culture was centrifuged and cells were suspended in 10 ml of SC-Ura + 2% galactose and grown at 22° for 16 hr. To detect endogenous Ty1i RNA, a single colony was suspended in 5 ml YEPD and grown at 30° overnight. The overnight culture was diluted different amounts depending on the growth of each strain (1:10 for WT, 1:4 for Δ, 1:5 for Δ, and 1:8 for Δ) in a total of 10 ml fresh YEPD and grown at 22° for 8 hr. Total RNA was extracted using the MasterPure yeast RNA purification kit (Epicentre Biotechnologies, Madison, WI) with modifications as described previously (Saha ). Poly(A)+ RNA was isolated from 250 μg total RNA using the NucleoTrap mRNA purification kit (Clontech, Mountain View, CA) following manufacturer’s protocol. Northern blot analysis using 32P-labeled riboprobes was performed as described previously (Saha ). Hybridization signals were visualized and quantified using a STORM 840 phosphorimager and ImageQuant software (GE Healthcare).

RNA sequencing analyses

RNA sequencing (RNA-seq) and ribosome footprint profiling (Ribo-seq) datasets were created by S. Komili, D. Muzzey, and F.P. Roth, and downloaded from GEO (accession no. GSE34438). Three replicate RNA-seq or Ribo-seq libraries were analyzed from the wild-type parent BY4741 and the isogenic Δ derivative (12 libraries in total). Note that deleting confers a transposition defect in BY4741 or BY4742 (Risler ), which are closely related MATa and MATα strains, respectively (Brachmann ). Sequencing reads were analyzed as follows: reads were trimmed of adaptor sequences. For Ribo-seq, a size filter of 28 bp was imposed, as 28-bp Ribo-seq reads tended to have a high proportion of in-frame reads. As Ty1 is present in multiple copies in the genome, the Ty1-H3 reference element (GenBank M18706.1) was used in the initial round of mapping. Mapping was performed with STAR RNA-seq aligner (version 2.4.2a) allowing zero mismatches. The remaining (non-Ty1) reads were mapped to the S. cerevisiae genome (Release 85; Ensembl). Mapped reads were assigned to genes and plotted with custom scripts. Gene counts are the average of the three replicates for each library, with similar results observed for pairwise library comparisons.

Protein extraction and western blot analysis

Total protein was extracted from strains expressing pGTy1 as previously described (Braiterman ). For detecting p22/p18, cells from 2 OD600 of culture were processed by trichloroacetic acid (TCA) extraction, and immunoblotted as described (Saha ). Samples were separated on 7.5 or 8% (for detecting IN), 10% (for detecting Gag-isoforms) or 15% (for detecting Gag-isoforms and p22/p18) SDS-PAGE gels. Antibody dilutions were as follows: anti-p18 1:5000, anti-VLP 1:7000, anti-IN 1:2500, anti-Pgk1 1:40,000, and anti-TY tag (BB2; UAB Epitope Recognition and Immunoreagent Core, Birmingham, AL) 1:20,000 (Bastin ). Immunoreactive protein signals were quantified using Quantity One software (Bio-Rad, Hercules, CA).

Stability of p22

p22 stability in wild type and Δ strains was determined as described previously (Doh ) with minor modifications; 5 ml of SC-Ura+2% raffinose medium was inoculated with a single colony and grown for ∼24 hr at 30°. The raffinose cultures (3 OD600 for HWA625 and 6 OD600 for HWA626) were added to 50 ml SC-Ura+2% galactose medium and grown at 20° for to an OD600 of 0.6–0.8. Thirty OD600 of cells per strain were harvested and washed three times in 5 ml of SC-Ura-Met+2% Gal. Click-iT L-HPG (Life Technologies) was added to 5 ml cell suspensions in SC-Ura-Met+2% Gal to a final concentration of 80 μM. Cultures were incubated at 20° for 30 min with shaking. Cells were washed twice with 5 ml of SC-Ura-Met+2% Gal and then resuspended in chase medium (SC-Ura-Met+2% Gal +50 mM l-methionine). Equal amounts of cells were harvested at 0, 1, 3, and 6 hr time points. Cell pellets were resuspended in 100 μl HB buffer (25 mM Tris, pH 7.5, 125 mM NaCl, 5 mM EDTA, 0.5% IPEGAL) with protease inhibitors (20 μl of 1 mg/ml Aprotonin, Leupeptin, Pepstatin, and 100 μl of 10 mg/ml PMSF) and vortexed at 4° for 5 min two times with 0.15 g of glass beads. Following addition of 5 μl of 10% SDS, the cell extract was boiled for 5 min and 500 μl ice-cold HB buffer was added. The extract was centrifuged for 20 min at 10,000 × g at 4° and precleared with 50 μl of Protein A agarose (Pierce) slurry (50% beads) washed with HB buffer and incubation at 4° for 1 hr; 10 μg of anti-p18 antibody was added to precleared lysate, and rotated at 4° for 16 hr, and then 50 μl of Protein A agarose slurry was added to each sample and rotated at 4° for 2 hr. Samples were washed three times with 500 μl HB buffer and proteins were eluted from the beads with 50 μl of 50 mM Glycine (pH 3.0) equilibrated in 200 μl of 50 mM Tris-Cl, pH 8.0, 1% SDS. Eluted proteins were precipitated by methanol/chloroform. Click-iT Cell Reaction Buffer Kit (Life Technologies) was used to conjugate TAMRA (Life Technologies) to HPG-labeled p22 according to the manufacturer’s protocol. Protein samples were separated on 15% SDS-PAGE gels. The fluorescent signal of TAMRA-conjugated p22 was detected using a Typhoon trio scanner with 580 BP30 filter, and p22 bands were quantified with ImageQuant Software (GE Healthcare). We assumed that p22 degradation follows first-order kinetics, and determined p22 half-life as described previously (Belle ).

Fluorescence in situ hybridization/Immunofluorescence (FISH/IF)

Retrosomes were analyzed by FISH/IF as described previously (Checkley , 2013; Saha ). To visualize endogenous retrosomes, 5 ml YEPD cultures were inoculated with a single colony and grown for 16 hr at 30°. The overnight cultures were diluted 40- (HWA215) or 500-fold (BY4742) into fresh YEPD, and grown at 20° for ∼24 hr to an OD600 of 0.8–1.0. To visualize pBDG606-induced retrosomes, 5 ml of SC-Ura+2% raffinose medium was inoculated with a single colony and grown for 24 hr at 30°. The raffinose cultures were diluted to 20- (HWA169) or 80-fold (HWA15) in SC-Ura+2% galactose medium and grown at 20° for ∼24 hr to an OD600 of 0.8–1.0. Formaldehyde was added directly to the culture to a final concentration of 4% and allowed to fix for 1 hr. Anti-VLP (rabbit polyclonal, 1:2000) and a cocktail of three GAG digoxigenin (DIG)-labeled antisense oligonucleotides (TyA_380, 5′-GCCTTCTCACATTCTTCTGTTTTGGAAGCTGAAACGTCTAACGGATCTTG-3′; TyA_444, 5′-TTCTCTGGAACAGCTGATGAAGCAGGTGTTGTTGTCTGTTGAGAGTTA-3′; TyA_545, 5′-CAACCAGATGGATTGGCTTGGTTTTGGGTCATCATGCACTGCTGTGGGTA-3′) were used to detect Gag and Ty1 mRNA, respectively. Secondary antibodies used here were anti-rabbit-AF594 (1:200; Life Technologies) and Fluorescein (FITC)-conjugated sheep anti-DIG Fab fragment (1:200; Roche Applied Science). Image acquisition was carried out using a ZeissAxio Observer microscope equipped with an AxioCam HSm camera, and images were analyzed with AxioVision v4.6 software (Carl Zeiss Microscopy, Thornwood, NY). The percentage of cells with retrosomes was calculated by dividing the number of cells with colocalized foci (Ty1 mRNA and Gag) by the total number of cells. At least 200 cells were analyzed for each strain.

Sucrose gradient sedimentation

Sedimentation of Gag complexes was analyzed as described (Tucker ), with the following modifications. For strains BY4742 and HWA215, a 100-ml culture was grown as described for FISH/IF microscopy. Equal amounts of total protein (7–8 mg in 300–450 μl) were applied to a 7–47% continuous sucrose gradient and centrifuged at 25,000 rpm in a SW 41 swinging bucket rotor (Beckman Coulter, Brea, CA) for 3 hr at 4°. Nineteen 0.6-ml fractions were collected, and an equal volume of each fraction was subjected to western blot analysis to detect Ty1 Gag as described above.

Nuclease protection

Strain BY4742 and HWA215 were grown as described above for FISH/IF microscopy. Nuclease protection using benzonase (EMD Millipore, Billerica, MA) was performed as described previously (Nishida ). Total RNA extraction and northern blot analysis was performed as described above.

Crosslinking Gag complexes

HWA15 and HWA169 were grown as described above for FISH/IF microscopy. Formaldehyde (Fisher Scientific, Waltham, MA) was added directly to the culture medium to a final concentration of 1%. Fixed cells were incubated for 1 hr at 20°. Proteins were extracted as previously described (Braiterman and Boeke 1994) and separated on NuPAGE 3–8% Tris-Acetate gels (Invitrogen) for detecting Gag multimers and 10% SDS-PAGE gels for detecting Pgk1. An equal amount of protein extract (7.5 µg) was loaded per lane and electrophoresis was performed at 120 V for 2 hr. Proteins were transferred to PVDF membranes and probed with anti-TY tag at a 1:50,000 dilution. Gag1 quantification was performed as described above.

Data availability

The authors state that all data necessary for confirming the conclusions presented in the article are represented fully within the article. All strains and reagents are available upon request.

Results

Identifying CNC-specific Ty1 cofactors

To identify CNC-specific genes, 498 Ty1 cofactors were analyzed using DAVID (http://david.abcc.ncifcrf.gov/), which enables functional annotation of large gene lists (Huang ). Functional annotation clustering was performed to determine the gene ontology (GO) terms enriched in Ty1 cofactor genes with reduced redundancy (for example, nucleic acid transport and RNA export were grouped in one cluster). Among the top four gene clusters with high enrichment scores (>4) (Table S2) were Ty1 modulators involved in post-transcriptional RNA biogenesis (RNA catabolic processes, RNA transport). To identify CNC-specific cofactors, 51 Ty1 cofactor mutants with defects in RNA metabolism were subjected to two secondary screens (Table 1). First, mutants involved in CNC should contain much less mature IN relative to Gag since CNC is associated with a block in the accumulation of mature Pol proteins (Garfinkel ; Matsuda and Garfinkel 2009). Second, CNC-specific mutations should affect Ty1 mobility to a greater extent in CNC+ strains containing Ty1 elements vs. CNC− strains that do not.

Table 1

Ty1 cofactors with roles in post-transcriptional RNA biogenesis

Classification	Gene name
P-body components	DHH1, LSM1, LSM6, PAT1, PUB1, XRN1
Nonsense-mediated decay	NAM7, NMD2, UPF3
Exosome components	LRP1, RRP6, SKI8
snRNA/snoRNA processing	BUD31, REF2, SNU66
RNA splicing	DBR1, LEA1, MUD2, SNT309, SQS1
CCR4-NOT complex	CCR4, MOT2, POP2
TRAMP complex	PAP2, TRF5
Nuclear cap-binding complex	CBC2, NPL3, STO1
RNA export/transport	APQ12, BUD13, KAP123, LOC1, LOS1, MFT1, NUP120, NUP133, NUP170, NUP188, NUP84, PML39, SCP160, SOL1, SXM1, THO1, THP2, TOM1
Miscellaneous	CTH1, MRT4, RPB4, SPT4, SSN3

Selected Ty1 cofactor mutants were assessed for the level of Ty1 Gag and IN by western blot analysis using VLP and IN antisera (Figure S1A and Table S3). Transposition defects were monitored by a qualitative papillation assay for Ty1 mobility (Curcio and Garfinkel 1991) (Table S3). To reliably detect mature Ty1 IN but not overwhelm CNC by overexpressing Ty1 (Figure S1B) (Garfinkel ), a low-copy centromere-based expression plasmid pGTy1 plasmid (pBDG606) was introduced into the Ty1 cofactor deletion mutants (Giaever ). For example, deletion of , a gene involved in biogenesis of 60S ribosomes and asymmetric localization of RNA (Long ; Urbinati ), resulted in the preferential loss of mature IN relative to Gag following galactose induction of pBDG606 (Figure 1). Also, when wild type and Δ cells were monitored qualitatively for Ty1 mobility after induction on SC-Ura + Galactose plates, the mutant harbored few Ty1 mobility events when compared to the wild-type parent. Two mutant classes were evident from the western blot analysis (Figure S1A) and qualitative mobility analyses (Table S3). We identified nine candidates (, ∆, , ∆, , , , ∆, and ∆) with considerably less mature IN relative to Gag and the Pgk1 loading control, and a low level of Ty1 mobility. Also, we identified mutants such as Δ and Δ that displayed a lower level of both Gag and IN. These mutants were not studied further as they may have less Ty1 mRNA or defects in Ty1 protein synthesis or stability.

Figure 1

Identifying Ty1 cofactor mutants that phenocopy CNC. Selected Ty1 cofactor mutants containing pGTy1his3-AI (pBDG606) the BY4742 background were induced on SC-Uracil (U) + Galactose (Gal) solid and liquid media (Figure S1A). Results from loc1Δ are shown as an example. The level of mature IN and Gag isoforms was assessed by western blot analysis using IN and VLP antisera. IN and Pol precursors, mature Gag-p45 and Gag isoforms (†) that include Gag-p49 are noted alongside the blot. Pgk1 (3-phosphoglycerate kinase) served as a loading control. Cell patches from SC-U+Galactose were replicated to SC-U-Histidine (H) to detect Ty1HIS3 mobility events, as monitored by the formation of His+ papillae. Mutants showing a severe decrease in mature IN but not Gag and a decrease in Ty1 mobility were selected for further analysis.

We deleted nine CNC-specific candidate genes in congenic S. cerevisiae CNC− and CNC+ strains DG3453 (+1 Ty1) and DG3648 (+20 Ty1s), respectively. During strain construction, mutants showing CNC-independent Ty1 mobility defects (), or no defect in Ty1 mobility (, ), clonal variability in Ty1 mobility (), or a severe growth defect () were not pursued further. The remaining four mutants (∆, ∆, ∆, and ∆) were assessed for Ty1 mobility (Figure S2A) and Ty1 protein level (Figure S2B) in the CNC− and CNC+ backgrounds. Mutants showing a decrease () or possible increase () of IN in the CNC− background or a relatively mild decrease of IN in the CNC+ background () compared to wild type were not pursued further. Only displayed both phenotypes expected for a CNC-specific mutant: a severe decrease in mature IN relative to Gag, and a decrease in Ty1 mobility in CNC+ but not in CNC− strains. Therefore, we analyzed how modulated Ty1 CNC in greater detail.

Ribosome biogenesis affects Ty1 CNC

is important for ribosome biogenesis and asymmetric RNA localization (Urbinati ; Heym and Niessing 2012). We hypothesized that, if ribosome biogenesis or asymmetric RNA localization were critical for CNC, disrupting additional genes in those pathways would phenocopy the mutant. For the mRNA localization pathway, and were selected as representatives (Heym and Niessing 2012). For ribosome biogenesis, the Ty1 cofactors and as well as their paralogs and were analyzed for CNC-specificity (Dakshinamurthy ; Risler ; Palumbo ). Deletion mutants were generated in S. cerevisiae CNC− and CNC+ backgrounds, and pBDG606-induced Ty1 mobility was determined (Table 2). The frequency of Ty1 mobility in CNC− and CNC+ strains was assessed individually, and also reported as a CNC−/CNC+ ratio, which estimates the increase in CNC observed in mutant strains when compared with wild type. The Δ mutant displayed a 3.4- and a 117.6-fold decrease in Ty1 mobility in the CNC− and CNC+ backgrounds, respectively. The increase in CNC as estimated by the Δ CNC−/CNC+ ratio was over 300-fold, which was ∼30-fold higher CNC than observed in the wild type. Interestingly, the ribosomal protein paralog mutants Δ and Δ, and Δ and Δ all showed an increase in CNC. However, the mutants identified as Ty1 cofactors (Δ and Δ) had higher levels of CNC than their paralogs (Δ and Δ), which were not recovered as cofactors. , previously identified as a Ty1 cofactor (Risler ), displayed a CNC ratio of 69. In contrast, Ty1 mobility and CNC remained at wild-type levels in cells lacking , which neither modulates Ty1 transposition nor affects ribosome biogenesis. In addition, deleting other genes (Δ, Δ, and Δ) involved in asymmetric localization of mRNA (Heym and Niessing 2012) did not affect Ty1 mobility in CNC− or CNC+ strains (Figure S3). Overall, our results show that defects in ribosome biogenesis trigger an increase in Ty1 CNC while defects in mRNA localization do not.

Table 2

Ty1his3-AI mobility in CNC− and CNC+ S. cerevisiae strains

Genotypea	Ty1his3-AI mobility, ×10−6 (SD)		Fold increase in CNC (CNC−/CNC+)
	CNC−	CNC+
Wild typeb	7481 (1805)	823 (318)	9
loc1Δ	2181 (534)	7 (2)	307
puf6Δ	1781 (748)	26 (4)	69
rps0aΔ	5002 (1585)	59 (48)	85
rps0bΔ	9592 (4285)	67 (22)	143
rpl7aΔ	5436 (1540)	89 (22)	61
rpl7bΔ	4212 (2308)	166 (55)	25
she2Δ	9370 (2038)	750 (231)	12

S. cerevisiae CNC− and isogenic CNC+ strains (WT and mutants) transformed with pBDG606 were induced with galactose. Cells were induced with galactose for 16 hr at 22°. Then cells were diluted, spread onto SC-Ura (for total number of colonies) and SC-Ura-His (for His+ colonies) plates, and incubated for 3–4 days at 30° until colonies formed.

Average of five trials.

Ty1 transcript and protein levels in candidate CNC mutants

Total RNA and protein from wild-type and Δ CNC− and CNC+ strains were subjected to northern (Figure 2A) and western blot (Figure 2B) analyses following pBDG606 induction. The level of Ty1 and Ty1 transcripts decreased modestly in wild type and Δ CNC+ strains when hybridized with a 32P-labeled riboprobe derived from GAG and normalized to actin () mRNA. The level of chromosomal Ty1 transcripts from the single copy CNC− strains was below the limit of detection, and the level of pGTy1-induced Ty1 mRNA decreased slightly in the CNC− Δ mutant. Western blot analysis showed that the levels of IN precursors and mature IN remained about the same in the Δ CNC− mutant. However, mature IN was not detected in the Δ CNC+ mutant, which validates the results obtained in the BY4742 knockout strain (Figure 1 and Figure S1A). The level of Gag isoforms was also similar in the CNC− and CNC+ Δ mutants. Together, our results suggest that the CNC-specific loss of mature IN is enhanced in a Δ mutant while Ty1 mRNA level is comparable to wild type.

Figure 2

Ty1 RNA and protein levels in congenic CNC− and CNC+ WT and loc1Δ strains expressing pBDG606. (A) Northern blot analysis of total RNA hybridized with a 32P-labeled riboprobe containing Ty1 nt 238–1702 to detect Ty1his3-AI and chromosomal Ty1 transcripts. ACT1 served as a loading control. Below are relative changes in RNA level compared to WT as determined by phosphorimage analysis. (B) Western blot analysis of Ty1 IN and Gag detected in total cell extracts (40 μg/lane) of CNC− and CNC+ loc1Δ mutant strains. Refer to Figure 1 legend for a description of Ty1 proteins. Pgk1 served as a loading control. Numbers below each panel indicate relative changes in protein level compared to WT as determined by quantification of protein bands by densitometry. ND, not detected.

To determine if the level of Ty1 Gag and IN proteins was similar in strains defective in ribosome biogenesis or RNA localization, total protein from wild type, Δ, Δ, Δ, Δ, and Δ strains expressing pBDG606 was subjected to western blot analysis (Figure 3). A trend emerged suggesting that mutants showing CNC-specificity (Δ, Δ, and Δ) (Table 2) contain less mature IN only in the CNC+ background, except the Δ mutant showed a decrease of mature IN in both CNC− and CNC+ backgrounds. Normal levels of Gag and mature IN were present in a CNC− and CNC+ Δ mutant. Therefore, likely modulates Ty1 transposition via its role in ribosome biogenesis and not through RNA localization.

Figure 3

Ty1 IN and Gag levels in ribosome biogenesis and RNA localization mutants. Total cell extracts from CNC− and CNC+ wild type and loc1Δ, rps0bΔ, rpl7aΔ, she2Δ, and puf6Δ mutant strains expressing pBDG606 were subjected to western blot analysis using IN and VLP antisera. Pgk1 served as a loading control. Refer to Figure 1 for a description of Ty1 proteins.

Ty1i RNA and p22/p18 levels are altered in CNC-specific mutants

Suresh reported that Ty1i RNA and p22/p18 levels increase in several ribosome biogenesis mutants including Δ. Therefore, endogenous levels of Ty1 mRNA, Ty1i RNA, Gag, and p22/p18 were determined in Δ, Δ, and Δ CNC mutants (Figure 4). To clearly distinguish Ty1i RNA from Ty1 mRNA (Saha ; Suresh ), poly(A)+ RNA was subjected to northern blot analysis using a 32P-labeled riboprobe containing Ty1 nucleotides 1266–1601 that hybridizes with Ty1 mRNA and Ty1i transcripts; mRNA served as a loading control (Figure 4A). Ty1 transcripts were below the limit of detection in wild type or CNC− mutants since there is only one Ty1 element present in the CNC− background. The level of Ty1 mRNA decreased fivefold in the CNC+ Δ background relative to wild type, whereas a similar decrease in Ty1 mRNA was not observed in a BY4742 Δ mutant. We did not examine this strain-specific effect further. Ty1 mRNA levels also varied modestly in the CNC+ Δ and Δ mutants. The level of poly(A)+ Ty1i RNA was altered in the CNC+ Δ, Δ, and BY4742 Δ mutants; slightly increased levels were detected in the CNC+ Δ and BY4742 Δ mutants, and a decreased level was detected in the CNC+ Δ mutant. The level of Ty1i RNA remained unchanged in a Δ mutant. Although there was no clear pattern of changes for Ty1i and Ty1 poly(A)+ transcripts in the ribosome biogenesis mutants, the CNC-specific mutants displayed a trend where the ratio of Ty1i/Ty1 mRNA increased twofold. We also compared the ratio of Ty1i RNA and Ty1 mRNA in each lane without normalizing to the mRNA loading control. In the CNC+ background, Ty1i/Ty1 mRNA ratio for wild type was 0.4, while the ratio for ribosome biogenesis mutants were 0.9, 0.8, and 0.8 for Δ, Δ, and Δ, respectively. In BY4742, the Ty1i/Ty1 mRNA ratio for wild type was 0.5, while the ratio for Δ was 0.7. Taken together, the transcript ratios support the idea that the enhancement of CNC correlates with an increase in Ty1i RNA over Ty1 mRNA.

Figure 4

Endogenous Ty1 transcript and p22 levels from CNC WT and ribosome biogenesis mutants, and BY4742 wild type and a loc1Δ mutant. (A) Northern blot analysis of poly(A)+ RNA hybridized with a 32P-labeled riboprobe containing Ty1 nt 1266–1601. This probe detects Ty1 mRNA, and Ty1 internally initiated (i) RNA that encodes p22/p18. ACT1 served as a loading control. The relative amount of Ty1 transcripts is shown below. (B) Western blot analysis of Ty1 and p22/p18 from total cell extracts using p18 antiserum, which detects Gag isoforms (†), p45, p22, and p18. Pgk1 served as a loading control. The relative amount of p22 and p18 is shown below.

To determine the level of endogenous Gag and p22/p18 in the CNC− and CNC+ strains, total protein was subjected to western blot analysis using an antiserum that detects both Gag and p22/p18 (Figure 4B) (Saha ). In accordance with the levels of Ty1 RNAs (Figure 4A), full-length Gag and a low level of p22/18 were detected in CNC+ wild type but not in CNC− wild type strains. In the CNC+ Δ mutant, we observed a decrease in full-length Gag and a 6.8-fold increase in p22 when compared to the CNC+ wild type. An increase in p22 of 3.5-fold, 4.2-fold, and 19.4-fold was observed in the CNC+ Δ and Δ mutants, and the BY4742 Δ mutant, respectively. Also, there was a 5- to 10-fold decrease in the level of p18 in the CNC+ ribosome biogenesis mutants. Some p18 was detected in the BY4742 Δ mutant, although p22 was present in a much higher amount. Recent work suggests that PR-mediated cleavage of p22 to p18 takes place in VLPs (Saha ; Tucker ). Therefore, the lower level of p18 observed in the ribosome biogenesis mutants is consistent with a defect in VLP assembly or function.

Since Ty1 Gag stability decreased in several Ty1 cofactor mutants defective in ribosome biogenesis and translation (Doh ), we examined p22 stability in wild type, and an isogenic Δ mutant in the BY4742 background by pulse-chase immunoprecipitation (Figure S4). Due to the low level of endogenous p22 and p18 in BY4742, a previously characterized galactose-inducible p22 expression plasmid pBDG1565 (pAUG1p22//2μ/) (Nishida ) was used to assess p22 stability. The half-life of p22 increased from 2.6 hr in wild-type cells to 4.1 hr in the isogenic Δ mutant. Although a ∼60% increase in half-life suggests that p22 may be more stable, this result does not fully account for the increased level of p22 in the CNC+ or BY4742 Δ mutants. We did not detect p18 in strains expressing pBDG1565, presumably because of the low-level of endogenous Ty1 protein processing by PR (Curcio and Garfinkel 1992).

Transcriptomic analyses in a loc1Δ mutant

To determine if the changes in Ty1 transcript levels are part of a global response to a defect in ribosome biogenesis, we analyzed publicly available RNA-seq and ribosome profiling (Ribo-seq) datasets generated in BY4741 and a Δ derivative (GEO accession no. GSE34438). First, RNA-seq datasets were analyzed for transcriptional changes. It is not possible to identify unique Ty1i reads from Ty1 as Ty1i transcripts overlap completely with full-length Ty1 mRNA (Saha ). We split the RNA-seq counts into two regions: (1) Ty1 exclusive and (2) Ty1 or Ty1i-derived reads (Figure 5, A and B). We observed a slight decrease in counts assigned to both regions in the Δ mutant. However, as Ty1 mRNA is more abundant than the Ty1i transcripts (Figure 4A) (Saha ; Suresh ), even a large fold change in Ty1i RNA might be undetectable by RNA-seq. In agreement with this idea, only a modest increase in Ty1i RNA relative to Ty1 mRNA was observed in the BY4742 Δ mutant as determined by northern analysis (Figure 4A).

Figure 5

Genomic analysis of RNA level and ribosome initiation in wild type and a loc1Δ mutant. (A) Transcripts and protein produced from Ty1. The element was divided into regions uniquely assignable to Ty1 (red), or Ty1 and Ty1i transcripts (blue). Three proteins are presented, two overlapping in the zero frame (Gag and p22), and a third produced after the +1 frameshift (Pol). These were assigned to Gag (green), Gag or p22 (light blue), or Pol (orange). (B) RNA-seq read counts for all genes (black dots) between wild type and loc1Δ. Each point is the average of three biological replicas for each strain. Highlighted are RNA read counts uniquely assignable to Ty1 (red), or Ty1 and Ty1i (blue). (C) Ribo-seq read counts for all genes (black dots) between wild type and loc1Δ. Analysis was restricted to reads 28 bp in length, and each point is the average of three biological replicas for each strain. Highlighted are read counts uniquely assignable to GAG in the zero frame (green), to either GAG or p22 in the zero frame (light blue), and reads from the POL region in the +1 frame (orange).

Previous work suggested that Loc1 represses translation of mRNA (Komili ). Therefore, we analyzed a Ribo-seq dataset generated in BY4741 and an isogenic Δ mutant to determine if a specific increase in translation initiation can account for the increase of p22 in the Δ mutant (Figure 5, A and C). The Ribo-seq reads were split according to region and frame: GAG (0-frame), GAG or p22 (0-frame), and POL (+1 frame). Ribo-seq coverage in these regions was not altered in the Δ mutant compared to wild type. Thus, there are no remarkable changes in gene expression profiles to suggest a reason for observed CNC phenotype or increased p22 levels in the Δ mutant.

There were 67 and 102 genes where RNA-seq transcript levels increased or decreased more than twofold in a Δ mutant, respectively, with a false discovery rate (FDR)-adjusted P-value (padj) <0.05 (Table S4). Functional annotation clustering was performed to find which GO terms are enriched among these transcripts (Table S4). In general, expression of genes involved in growth-related processes such as cell wall structure, glycolysis/gluconeogenesis, ion transport, and amino acid biosynthesis were affected by the absence of . It remains to be determined whether these effects represent primary or secondary consequences of deleting .

Gag–Gag interactions and Ty1 mRNA packaging are altered in a loc1Δ mutant

Formaldehyde cross-linking and the sedimentation pattern of Gag complexes was used to investigate the role of in VLP assembly. Formaldehyde treatment of cells captures HIV-1 Gag assembly intermediates at the plasma membrane (Kutluay and Bieniasz 2010). However, our application is simplified since retrosome/VLP assembly takes place in the cytoplasm (Checkley ; Dutko ; Doh ). To help ensure there are comparable levels of Gag prior to crosslinking, we induced pBDG606 expression in wild type BY4742 and an isogenic Δ mutant (Figure 6). Following galactose induction, yeast cells were treated with formaldehyde, and total cell extracts were separated on a 3–8% gradient gel. Gag complexes and a Pgk1 loading control were detected by western blot analysis. In the absence of formaldehyde, Ty1 Gag was mostly present in a monomeric form (Gag1) (Figure 6A). We also detected a small amount of dimeric Gag (Gag2) and the 199 kDa Gag-Pol precursor (denoted by an asterisk). In wild-type cells treated with formaldehyde, there was an ordered accumulation of Gag multimers with Gag2 and trimeric Gag3 as the most prominent forms. When the Δ mutant was treated with formaldehyde, a comparable amount of Gag1 was present when compared to wild type; however, there was a lower level of Gag2 and a marked decrease in higher order multimers, which is also evident in densitometric tracings (Figure 6B). There were also similar levels of Gag1 in wild type and Δ cells when compared with Pgk1. Our results suggest that facilitates early steps in VLP assembly as monitored by Gag multimerization.

Figure 6

Gag multimerization in wild type and loc1Δ strains. (A) Total protein from pBDG606-induced WT and loc1Δ strains were cross-linked by formaldehyde treatment in vivo. Protein extracts from untreated (− formaldehyde) and treated (+ formaldehyde) cells were separated on a 3–8% Tris-Acetate gel and probed with Gag antiserum (TY-tag). Equal amounts of protein separated on a 10% SDS-PAGE gel were probed with Pgk1 antiserum, which served as a loading control. Estimated protein sizes are shown in kilodaltons. Extracts from untreated cells contain monomeric Gag (Gag1) and the Gag-Pol-p199 precursor (*), while cross-linked cells produce a ladder of Gag multimers. The gel conditions allow sufficient separation of Gag1–4, while larger multimers (Gag5+) are not resolved. An additional faster-migrating Pgk1 signal appeared in formaldehyde-treated samples, which likely results from cross-linking. Numbers below the panel indicate relative changes in Gag1 level compared to wild type with and without formaldehyde treatment, as determined by quantification of protein bands by densitometry. Both Pgk1 bands in formaldehyde treated cells were used in the quantification. (B) Densitometry traces of cross-linked Gag from WT and loc1Δ. The pixel intensity for each Gag-multimer per lane is shown along the Y-axis.

To visualize more complex assembly intermediates, total protein extracted from BY4742 wild type and Δ cells was fractionated through 7–47% sucrose gradients, and the sedimentation pattern of endogenous Gag was monitored by western blot analysis (Figure 7A). The majority of Gag from wild type cells reproducibly sedimented as larger complexes (fractions 11–19) likely comprised of VLP assembly intermediates, since few endogenous VLPs are detected in the absence of pGTy1 expression (Checkley ). In contrast, Gag from the Δ mutant reproducibly sedimented as smaller complexes (fractions 7–10), suggesting that higher order Gag complex formation is also defective in a Δ mutant.

Figure 7

Gag sedimentation and nuclease protection of Ty1 mRNA in wild type and loc1Δ strains. (A) Total cell extracts from BY4742 WT or a loc1Δ mutant were centrifuged through a 7–47% continuous sucrose gradient. Ten micrograms of protein extract (input) and equal volumes of each fraction (15 μl/lane) were subjected to western blot analysis using Gag antiserum (TY-tag), which recognizes Gag isoforms and Gag-p45. Fractions containing the highest concentrations of Gag, as determined by densitometry tracing, are underlined. Sucrose gradient analysis was repeated three times and a representative experiment is shown here. (B) Ty1 mRNA packaging as monitored by sensitivity to the nuclease benzonase. Equal aliquots of whole cell extract from BY4742 WT and a loc1Δ mutant were incubated with (+) or without (−) benzonase. To detect Ty1 mRNA, total RNA extracted from these samples was analyzed by northern blot analysis. Nuclease protection assays were repeated three times and a representative experiment is shown. WT protection ranged from 26.5 to 31.4%, and loc1Δ protection ranged from 16.1 to 18.8%. ACT1 was used as a control to confirm RNA degradation in the benzonase treated samples.

To examine if Ty1 mRNA packaging into nuclease-resistant Gag assemblies is defective (Purzycka ; Nishida ), cell extracts from BY4742 wild type and Δ were treated with the nuclease benzonase followed by northern blot analysis (Figure 7B). Hybridization signals from the benzonase treated samples were normalized to untreated controls, and mRNA served as a control for RNA integrity and benzonase activity. Ty1 mRNA level reproducibly decreased about threefold when extracts from wild type cells were treated with benzonase, whereas Ty1 mRNA decreased more than sixfold when extracts from the Δ mutant were treated with the nuclease. These results suggest that less Ty1 mRNA is packaged into nuclease resistant Gag complexes in the absence of .

Retrosomes are not detected in loc1Δ cells

BY4742 wild type and Δ cells were subjected to FISH/IF to visualize endogenous (Figure 8) or pBDG606 induced retrosomes (Figure S5). VLP antiserum was used to detect Gag, and a cocktail of DIG-labeled oligonucleotide GAG probes was used to detect full-length Ty1 mRNA (Saha ). About 34% of wild type cells contained retrosomes, which are defined as distinct foci containing Gag and Ty1 mRNA (Malagon and Jensen 2008; Checkley ). Retrosomes were not detected in the Δ mutant. We also failed to observe nondistinct punctate staining for both Ty1 mRNA and Gag, which is present in other Ty1 cofactor mutants (Checkley ; Dutko ). Instead, diffuse cytoplasmic staining was detected in the Δ mutant. Furthermore, retrosomes were also absent when pBDG606 was expressed in Δ cells (Figure S5). Taken together, our results indicate that steps leading to assembly of functional VLPs are disrupted in a Δ mutant.

Figure 8

Retrosome formation in wild type and loc1Δ cells. BY4742 WT and loc1Δ mutant strains were visualized by DIC and fluorescence microscopy. DNA was stained with DAPI, Ty1 mRNA was detected using a cocktail of GAG-DIG probes, and Gag was detected with VLP antiserum. The percentage of cells containing retrosomes, defined as having at least one colocalized Ty1 mRNA/Gag foci, is shown in the DAPI/Ty1 mRNA/Gag merge. The experiment was repeated twice and >200 cells were analyzed for each strain.

Discussion

Genome-wide screens have revealed cellular modulators of Ty1 retrotransposition, and some are conserved with cellular factors involved in retroviral replication (Maxwell and Curcio 2007; Goff 2008). Since all the genetic screens for Ty1 cofactors (Griffith ; Dakshinamurthy ; Risler ) were performed in CNC+ strains containing genomic Ty1 elements and producing the p22 restriction factor (Saha ), a subset of previously identified Ty1 cofactors may affect the level or activity of p22. RNA metabolism is central to the process of retrotransposition; therefore, cellular processes affecting the Ty1 life cycle between RNA transcription and VLP assembly (Ty1 RNA export, translation, and retrosome formation) may contribute to CNC. Here, we utilized a candidate gene approach to identify CNC-specific modulators based on informative secondary screens that monitor preferential loss of mature IN and Ty1 mobility in a CNC+ strain compared to an isogenic CNC− strain. One weakness of a candidate gene approach is that not all possible causative genes may be queried (Zhu and Zhao 2007). However, the hundreds of known Ty1 modulators (Scholes ; Griffith ; Nyswaner ; Dakshinamurthy ; Risler ) provide a rich source of interactions between Ty1 and its host that can be mined by additional phenotypic analyses. Our analysis of CNC-specific modulators revealed an unexpected connection between Ty1 CNC and a large group of host cofactors involved in ribosome biogenesis and function (Suresh ).

The secondary screens identified as a CNC-specific Ty1 cofactor since the level of mature Ty1 IN and retromobility decreased much more in CNC+ than in CNC− strains. encodes a nucleolar protein that binds RNA (Long ), associates with the pre-60S ribosome complex (Harnpicharnchai ), and is necessary for rRNA processing, 60S assembly and nuclear export of pre-60S subunits (Urbinati ). Loc1 also helps localize mRNA to the distal tip of anaphase cells (Long ) via a direct interaction with She2 (Niedner ). To determine whether modulated Ty1 transposition by altering RNA localization or ribosome biogenesis, we analyzed two other genes implicated in mRNA dynamics, and , for their ability to phenocopy the alteration of CNC observed with . Disrupting leads to complete loss of mRNA localization, and Puf6 interacts with Loc1 in the nucleolus, and has a role in translational repression of mRNA (Heym and Niessing 2012). Our data suggest that the mRNA localization pathway is not involved in CNC because deleting a key component of the pathway, , does not affect Ty1 mobility, confer a CNC-specific transposition defect, or result in loss of IN. Furthermore, several other genes implicated in mRNA localization (, , and ) do not affect Ty1 CNC (Figure S3).

Like , is involved in both mRNA localization as well as 60S ribosome biogenesis (Li ). The Δ mutant showed a decrease of IN in both CNC− and CNC+ strains, but its transposition defect was CNC-specific. Perhaps is required for the process of Ty1 retrotransposition, and is also involved in CNC. Interestingly, Loc1 function may be compromised in the absence of Puf6 since these proteins interact (Shen ; Yang ). While their interaction is important for mRNA localization and 60S ribosome assembly, previous studies suggest that Loc1 and Puf6 carry out distinct functions. For example, Δ, Δ, and a Δ Δ double mutant each display different phenotypes with respect to the localization of ribosomal protein paralogs, and assembly of Rpl43 in the 60S ribosome subunit (Komili ; Yang ).

Analysis of focused our attention on how ribosome biogenesis affects Ty1 transposition and CNC. Remarkably, 71 of 458 host genes identified in multiple genetic screens for Ty1 modulators (Scholes ; Griffith ; Nyswaner ; Dakshinamurthy ; Risler ) encode ribosomal subunits, ribosome biogenesis factors, and translation factors. Recent work suggests that these genes affect multiple steps during transposition, including translation initiation, programmed frameshifting, protein stability, and subcellular protein localization (Suresh ). A rare tRNA-Arg(CCU) gene specifically involved in programmed frameshifting is also required for Ty1 mobility (Kawakami ). Δ, Δ, Δ, and mutants all show transcriptional patterns where Ty1i RNA increases relative to Ty1 mRNA. Additionally, Ty1 mRNA and Gag localization in retrosomes is absent in a Δ mutant, further implicating ribosome biogenesis in Ty1 retrotransposition (Doh ).

Ribosome biogenesis and nucleolar functions influence the replication of retroviruses and other retrotransposons. In addition to its primary function in ribosome biogenesis, the nucleolus is involved in ribonucleoprotein remodeling, cell cycle progression, and stress response (Boisvert ). Cellular stress and viral infection change the organization, composition, and morphology of the nucleolus (Boulon ). The nucleolus is the target for viral infection by RNA viruses (Hiscox 2007; Salvetti and Greco 2014), and nucleolar trafficking of viral proteins occurs with retroviruses such as Rous Sarcoma Virus (RSV), Mouse Mammary Tumor Virus (MMTV), and HIV (Jarboui ; Beyer ; Lochmann ). Ty3 and L1 retrotransposon proteins (Lin ; Goodier ) and Alu RNA (Caudron-Herger ) display nucleolar localization. Moreover, ribosomal protein Rpl7 affects HIV-1 Gag nucleic acid chaperone activity (Mekdad ), and Rpl9 helps traffic MMTV Gag proteins (Beyer ; Lee ; Mekdad ). Taken together, our work provides additional evidence linking ribosome biogenesis and retroelement replication.

We expanded the CNC analysis to include the ribosomal protein gene paralogs and , and and . As is the case for other paralogs, only and were recovered as Ty1 cofactors (Dakshinamurthy ; Risler ). Δ and Δ mutants showed loss of IN and a lower level of Ty1 mobility in CNC+ but not in CNC− strains, demonstrating that these mutants are CNC-specific (Table 2). However, the Δ and Δ paralogs also displayed CNC-specificity but at a lower level than Δ and Δ. Expression of these ribosomal protein genes may account for these results. The level of mRNA is fourfold higher than , and mRNA is 1.5-fold higher when compared to (Komili ), which correlates with CNC−/CNC+ ratios of Ty1 mobility (Table 2). Based on mRNA expression, Rpl7a may be the predominant ribosomal protein in the cell. Therefore, the absence of Rpl7a caused a more severe CNC-specific transposition defect. However, losing Rpl7b does not affect CNC because there is still a substantial amount of Rpl7a present in a cell. For Rps0a and Rps0b, comparable amounts of protein may be expressed from each gene. The absence of one RPS0 paralog could impact CNC, as suggested by the mobility defect observed in both CNC+ Δ and Δ mutants. Our results also agree with recent work regarding dosage-dependent rather than paralog-specific function of Rpl7a and Rpl7b on different cellular phenotypes, including Ty1 life cycle stages (Palumbo ). Taken together, any paralog-specific effects of ribosomal proteins on CNC is probably due to expression levels of the duplicated ribosomal protein genes.

Our results did not identify a step during ribosome biogenesis that impacts Ty1 CNC. Although was analyzed in greatest detail as a Δ mutation conferred robust CNC-specificity, both small and large subunit ribosomal protein mutations as well as a ribosome assembly factor mutation also affect CNC to varying degrees. Instead of responding to specific steps in ribosome biogenesis, Ty1 CNC may be modulated by ribosome availability. Polysome profiles of several mutants analyzed for Ty1 CNC suggest there is an imbalance of 40S and 60S subunits. A decrease in 40S and an increase in 60S subunits is observed in a Δ mutant (Suresh ), a decrease of both 40S and 60S subunits occurs in a Δ mutant (Yang ), and a decrease in 60S subunits, 80S monosomes, and polysomes is detected in a Δ mutant (Palumbo ). Defects in functional ribosome production affect cell growth and proliferation (Woolford and Baserga 2013). Interestingly, there is also a correlation between ribosome availability, cell growth, and Ty1 transposition. Δ and Δ mutants grow slowly, contain decreased amounts of 40S and 60S subunits (Yang ), and display a high CNC−/CNC+ ratio (Table 2). Growth of RPL7 paralog mutants also show a similar relationship with Ty1 CNC. An Δ mutant grows more slowly, and impacts Ty1 CNC more than a Δ mutant, which contains more Rpl7a due to the higher expression of (Palumbo ). In addition, the slow growth of ribosomal protein mutants is related to other cellular phenotypes such as enhanced resistance to endoplasmic reticulum stress (Steffen ). Although the magnitude of CNC may correlate with the level of functional ribosomes, other CNC phenotypes may not be directly related. For example, retrosomes are not detected in Δ (Figure 8) and Δ mutants (Doh ); however, a more severe transposition defect is conferred by a Δ mutant (Table 2). Together, our data support the idea that different steps in the Ty1 life cycle have different thresholds for ribosome levels.

The observed CNC-specific phenotypes in Δ and two ribosomal protein mutants likely results from an increased amount of p22 when compared to the isogenic wild type strain (Figure 4B). We examined steps in Ty1 gene expression and protein stability to determine what might account for the increase in p22 observed in several ribosome biogenesis mutants. Our data permit the following conclusions:

The level of poly(A)+ Ty1 mRNA and Ty1i RNA is altered in the absence of , , and , and displays a trend suggesting that more Ty1i RNA is present relative to Ty1 mRNA (Figure 4A). However, the decrease in the level of poly(A)+ Ty1i RNA in the Δ and Δ mutants or the increase in the level of Ty1i RNA observed in other mutants fails to explain the 3.5 to 19.4-fold increase in the level of p22. Perhaps nonpolyadenylated Ty1i transcripts can also be used to synthesize p22 in the ribosome biogenesis mutants (Saha ; Suresh ).

Loc1 is reported to function as a translational repressor (Komili ); therefore, loss of Loc1 may enhance the initiation of p22 translation. However, Ribo-seq analysis of p22 in a Δ mutant is comparable to other Ty1 proteins as well as cellular proteins (Figure 5C). Although this result raises the possibility that initiation of translation is not enhanced in a Δ mutant, it is unclear why the Ribo-seq analysis does not account for the increased level of p22 observed in Δ as well as the Δ and Δ mutants (Figure 4B). Possible reasons for this discrepancy include: (a) the effects on translation initiation are obscured by the overlapping nature of RNAs and proteins in this region, (b) the experimental conditions or strains used for the libraries differ from those used in our work, or (c) Loc1 may exert its effects through some non-RNA or nontranslational mechanism.

Loc1 enhances p22 turnover. Pulse-chase immunoprecipitation analysis suggests that p22 stability increases in a Δ mutant (Figure S4). However, a ∼60% increase in p22 stability does not completely explain the increase in p22 level (Figure 4B).

What could be causing the transcriptional changes of Ty1 mRNA and Ty1i RNA in the ribosome biogenesis mutants? Perhaps defects in ribosome biogenesis cause a form of stress that leads to changes in Ty1 transcript level. A nucleolar stress response pathway activates p53-dependent as well as p53-independent downstream pathways in mammalian cells (James ). Nucleolar stress may also occur in yeast (Gomez-Herreros ; Thapa ), though the signaling pathway has not been defined. We considered the possibility that the environmental stress response (ESR), which alters expression of a core set of genes under diverse environmental changes (Gasch and Werner-Washburne 2002) is activated in a Δ mutant. However, analysis of the Δ and ESR transcriptomes (Gasch ) reveals little overlap with the ESR (Table S4). Interestingly, stress-responsive transcription factors such as Msn2, Tye7, and Gcn4 affect Ty1 mRNA expression may also be involved in Ty1i RNA transcription. Msn2 binds to Ty1 elements in vivo and is capable of acting as a transcriptional activator as well as a repressor (Schmitt and McEntee 1996; Elfving ). Tye7 is a bHLH transcription factor that activates Ty1-mediated gene expression and regulates Ty1 mRNA and antisense RNAs expression upon adenine starvation (Lohning and Ciriacy 1994; Servant ). Gcn4 is a transcription factor that responds to amino acid starvation via the TOR pathway, and can activate Ty1 transcription by recruiting chromatin remodelers Swi/Snf and SAGA (Hinnebusch and Natarajan 2002; Morillon ). Further work will be required to understand the relationship between factors responsible for Ty1 mRNA and Ty1i RNA transcription.

Several Ty1 defects in cells lacking or ribosomal proteins can be reconciled by an increase in the level of p22 (Saha ). These include a marked decrease in mature IN relative to Gag (Figure 2B and Figure 3), altered Gag–Gag interactions (Figure 6 and Figure 7), and formation of endogenous (Figure 8) or pGTy1-induced retrosomes (Figure S5). Previous work (Saha ), along with data presented here, led us to consider retrosome formation in CNC− and CNC+ cells. In CNC− cells, retrosomes may not form, or are undetectable, due to low levels of Ty1 products, including p22. The Ty1 copy number dependence of CNC, however, suggests the low level of p22 present in cells containing just a few elements may be less inhibitory (Garfinkel , 2016). In CNC+ cells, the increase in retrosome abundance may sensitize Ty1 retrotransposition to small changes in p22 levels and disturbances in the retrosome environment or network of interactions.

Recent work on biomolecular condensates (Banani ; Riback ), also known as membrane-less organelles or quinary structures, have influenced our view of retrosomes and their interactions with other cellular components. We postulate that Ty1 retrosomes are biomolecular condensates as they contain Ty1 RNA and proteins, and likely cellular proteins in an undefined assembly, no fixed membranous structure, and are highly dynamic. Phase separation caused by multivalent macromolecular interactions of nucleic acids and proteins drive the formation of molecular condensates. Indeed, P-bodies and stress granules are biomolecular condensates and required for Ty1 transposition. Retrosomes disassemble in P-body component as well as stress granule component mutants (Checkley ; Dutko ; Doh ). The highly dynamic nature of retrosomes is also revealed by their disassembly under conditions, such as glucose deprivation, that causes P-bodies to form (Checkley ). Our work favors the idea that normal ribosome biogenesis and functional ribosomes contribute to the retrosome environment required to maintain these condensates in a phase-separated or demixed state. The continuum from retrosome formation to VLP assembly may be similar to the formation of condensates such as stress granules (Riback ) or transcriptional super-enhancers (Hnisz ) from ribonucleoprotein complexes or “normal” enhancers, respectively. As various chemical and environmental factors or stress conditions influence the structure and function of biomolecular condensates (Mitrea and Kriwacki 2016), changes in the cellular environment resulting from defects in ribosome biogenesis in combination with increased p22 levels contribute to altered retrosome dynamics and defects in Ty1 mobility.

In summary, our study demonstrates that ribosome biogenesis is critical for Ty1 CNC. The diverse defects in ribosome biogenesis that enhance Ty1 CNC raise the possibility of a cellular response that leads to an increase in the amount of p22 relative to its target Gag. Further work will be required to define a stress response in yeast that affects p22 abundance and retrosome formation, characterize the retrosome interactome both genetically and biochemically, and determine what triggers VLP assembly from retrosomes.

Supplementary Material

Supplemental material is available online at www.genetics.org/lookup/suppl/doi:10.1534/genetics.117.300388/-/DC1.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.