{"paper_id":"02ec39c0-bd29-4b97-99c7-8d6962e9e31e","body_text":"Reverse transcribed ssDNA derepresses translation of a 1 \nretron antiviral protein 2 \nKaren Zhang1,2†, Matías Rojas-Montero1†, Darshini Poola1†, Josepha M Klas3, Arturo Carabias4, 3 \nDennis J Zhang3, Mario R Mestre3, Ruiliang Zhao3, Ana Dávila-Hidalgo4,5, Guillermo Montoya4, Rafael 4 \nPinilla-Redondo3, Alejandro González-Delgado1*, Seth L Shipman1,2,6*  5 \n 6 \n 7 \n1 Gladstone Institute of Data Science and Biotechnology, San Francisco, CA, USA 8 \n2 Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, 9 \nCA, USA 10 \n3 Section of Microbiology, Department of Biology, University of Copenhagen, 2100 Copenhagen, 11 \nDenmark 12 \n4 Structural Molecular Biology Group, Novo Nordisk Foundation Centre for Protein Research, 13 \nDepartment of Cellular and Molecular Medicine, Faculty of Health and Medical Sciences University 14 \nof Copenhagen, Blegdamsvej 3B, Copenhagen, 2200, Denmark 15 \n5 Centro de Investigación del Cáncer, Consejo Superior de Investigaciones Científicas (CSIC), 16 \nUniversidad de Salamanca, 37007 Salamanca, Spain 17 \n6Chan Zuckerberg Biohub, San Francisco, CA, USA 18 \n* Corresponding authors. Email: seth.shipman@gladstone.ucsf.edu, 19 \nalejandro.gonzalez@gladstone.ucsf.edu 20 \n† These authors contributed equally to this work  21 \n 22 \n 23 \n 24 \nABSTRACT  25 \nRetrons are bacterial immune systems that prevent the spread of phages by initiating a toxic 26 \nresponse within infected hosts. All previously characterized retrons produce high levels of 27 \nmulticopy single-stranded DNA (msDNA) in the cell by reverse transcription, which acts as an 28 \nantitoxin in the absence of phage infection. However, we describe here a non-canonical mechanism 29 \nfor Type VI retrons, which do not produce detectable msDNA in the absence of phage, yet still 30 \nprovide phage defense. Focusing primarily on Retron-Vpa2, a Type VI retron from Vibrio 31 \nparahaemolyticus, we show broad defense against phages and identify triggers of the system within 32 \nphage recombination systems. Within the Retron-Vpa2 operon, we find a highly enriched, 33 \nstructured transcript that we term a hybrid RNA (hyRNA), which contains both the retron’s reverse 34 \ntranscription template and a translationally repressed toxic effector coding sequence. We find that 35 \nphage infection induces the accumulation of high levels of msDNA and that this msDNA is necessary 36 \nfor derepressing translation of the antiviral toxin. These findings present key biological and 37 \nmechanistic insights into a distinct group of retrons while highlighting the diversity of systems that 38 \nparticipate in bacterial immunity. 39 \n 40 \n 41 \n 42 \n  43 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nINTRODUCTION 44 \nBacteria have evolved an arsenal of antiviral defense mechanisms to protect themselves from 45 \nphages, including a class of defense systems called retrons1–3. Retrons operate via abortive infection 46 \n– upon sensing phage, they release a toxic effector that induces host cell death or dormancy before 47 \nthe phage can propagate, sparing the uninfected bacterial population2–5. Retrons are typically three-48 \npart systems, consisting of a reverse transcriptase (RT), a noncoding RNA (ncRNA), and a toxic 49 \neffector1–3,6. The RT recognizes a highly structured msr region on the ncRNA and reverse 50 \ntranscribes an adjacent msd region into multicopy single-stranded DNA (msDNA)7–13. The 51 \naccumulation of high levels of msDNA by ongoing reverse transcription has been considered a 52 \nhallmark of retrons12–16. In fact, the distinctive bands that the msDNA generates on a 53 \npolyacrylamide gel enabled the discovery of the earliest identified retrons in the 1980s and are 54 \nused as a signature for Retron-containing species to this day13–19.  55 \n 56 \nIn canonical retron systems, the msDNA forms a complex with the remaining ncRNA and RT, which 57 \ntogether interact with the effector protein as an antitoxin to neutralize its toxic effect3,20–27. While 58 \ntoxin-antitoxin motifs are common in phage defense, the incorporation of reverse transcribed DNA 59 \nin the antitoxin is unique to retrons, enabling these systems to respond rapidly to phage-encoded 60 \nDNA interacting proteins, such as nucleases3,28–30, single-stranded binding proteins28, and 61 \nmethyltransferases3,16,22. During infection, these phage proteins degrade, perturb, or modify the 62 \nmsDNA, leading to release of the toxin and activation of the defense response3,16,22,28,29.  63 \n 64 \nRetrons are a particularly diverse class of defense systems. A bioinformatics study encompassing 65 \nover 1900 retrons from multiple bacterial phyla classified retrons into thirteen types based on 66 \ndifferent protein domains found in the highly variable toxic effector component6. Following this 67 \nbioinformatic work, an experimental census was conducted to survey retrons from all thirteen 68 \ntypes, including characterization of msDNA production15. This census revealed an unexpected trait 69 \namong a previously untested group: Type VI retrons failed to produce detectable msDNA.  70 \n 71 \nHere, we present the first mechanistic study of Type VI retrons, finding that this type deviates 72 \nsubstantially from the canonical mechanism described in previously studied retrons. These retrons 73 \ndefend against phage despite the lack of detectable msDNA at baseline. We reveal that the msr-msd 74 \nand effector coding sequence occupy a contiguous transcript which we call the hybrid RNA 75 \n(hyRNA). Further characterization of the operon shows that the hyRNA secondary structure and an 76 \nadditional accessory protein result in translational repression of the toxic effector. Counter to all 77 \nother retrons, reverse transcription of the Type VI msDNA is induced during phage infection, 78 \nspecifically by phage-encoded recombination-associated proteins. Rather than serving as an 79 \nantitoxin, the Type VI msDNA derepresses the translation of the effector. This inverted triggering 80 \nmechanism constitutes an intriguing departure from our current understanding of retrons, 81 \nprompting a reconsideration of the mechanistic range of bacterial immunity.   82 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nRESULTS 83 \nRetron-Vpa2 confers broad phage defense without detectable msDNA 84 \nType VI retrons have a characteristic operon that begins w ith a predicted ncRNA region containing 85 \nthe msr-msd, followed by an ORF encoding a small protein (SP) of unknown function (Fig 1a). 86 \nDownstream of this is another unknown accessory protein containing a helix-turn-helix (HTH) 87 \ndomain, followed by the retron reverse transcriptase (RT).  88 \n 89  \nType VI retrons form a well-supported monophyletic clade (clade 3) within the phylogeny of retron 90 \nRTs6 together with the Type IV (e.g. Retron-Eco62) and Type V retrons (e.g. Retron-Sen131,32). In a 91 \npreviously conducted experimental census of all retron types, we found no msDNA production from 92 \nType VI retrons when expressed in E. coli15. Here, we replicated that result (Fig 1c). Five of the 93 \nretrons were retested exactly as in the retron census, with expression of just the retron RT and 94 \npredicted ncRNA (Retrons-Psp1, -Bfr1, -Tde1, -Psp2, and –Cle1). For two of the retrons, we 95 \nsynthesized and expressed the entire retron operon to test whether omission of a component other 96 \nthan the RT and ncRNA in the previous census might have accounted for the lack of msDNA. Of the 97 \ntwo retrons tested as a full operon, one had been previously included in the census in the more 98 \nminimal form (Retron-Vpa2 from Vibrio parahaemolyticus), while the other was newly included to 99 \ncheck whether matching the host species might yield detectable msDNA (Retron-Eco12 from E. 100 \ncoli). Yet, we again observed no msDNA for any Type VI retron tested in any format. 101 \n 102 \nWe next tested whether these retrons confer defense against phage infection. We expressed each 103 \nunder their native promoters in an E. coli BL21-AI derivative with its endogenous Retron-Eco1 104 \ndeleted (henceforth bSLS.114)33. Despite the lack of msDNA production, we observed robust 105 \ndefense against a virulent strain of phage lambda with each of the Type VI retrons (Fig 1d). For 106 \nRetron-Vpa2, we then tested for defense in an expanded panel of 92 phages, including the BASEL 107 \ncollection34 (Fig 1e). We observed broad defense against the Demerecviridae, Myoviridae, and 108 \nDrexlerviridae families of phages.  109 \n 110 \nWe selected several phages that induced strong defense phenotypes (Bas10, Bas34, Bas53, Muut, 111 \nand Stevie_ev116) and infected liquid cultures of Retron-Vpa2-expressing E. coli MG1655 at a range 112 \nof MOIs from 0 to 10 (Fig 1f). We observed a clear abortive infection phenotype in Retron-Vpa2-113 \nexpressing cells for phages Bas10, Bas34, and Stevie_ev116, where the population persists through 114 \nphage infection in a MOI-dependent manner. For phages Bas53 and Muut, cells expressing the 115 \nretron grow equally well regardless of MOI. By contrast, cultures without the retron exhibit collapse 116 \nat every MOI. Overall, these data indicate that Retron-Vpa2 confers protection to the bacterial 117 \npopulation in the presence of phage.  118 \n 119 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n 120 \nFigure 1 – Retron-Vpa2 confers broad phage defense without detectable msDNA. 121 \nA) Schematic of the Type VI retron operon. B) Location of Type VI retrons and their close relatives in clade 3 of a retron RT 122 \nphylogeny. C) PAGE analysis for msDNA abundance comparing several Type VI retrons with Retron-Eco1 serving as a positive 123 \ncontrol for msDNA. A single-stranded DNA ladder with lengths marked in nucleotides is shown for reference. D) Spot assay 124 \nshowing titration of phage lambdavir on Retron-Eco12- and -Vpa2-expressing strains of E. coli bSLS.114 relative to an empty 125 \nvector control. Efficiency of plating is quantified as pfu/mL in the adjacent bar graph (one-way ANOVA P=0.0152, Dunnett’s 126 \ntest corrected for multiple comparisons versus empty vector: Retron-Eco12 P=0.0181; Retron-Vpa2 P=0.0179). E) Fold 127 \ndefense conferred by Retron-Vpa2 against a panel of 92 phages, depicted by shaded boxes. Colored bars denote phage family, 128 \nsubfamily, and genus. F) Growth curves of E. coli MG1655 expressing either Retron-Vpa2 or an empty vector, infected with 129 \nfive different phages at MOIs of 0, 0.01, 0.1, 1, and 10. Additional statistical details in Supplementary Table 1.  130 \n 131 \n 132 \n 133 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n 134 \n 135 \nRecombination systems trigger Retron-Vpa2 in two phages 136 \n 137 \nTo identify factors responsible for triggering Retron-Vpa2 defense, we examined phage escapees 138 \nfor lambdavir and Stevie_ev116. Whole genome sequencing of the lambdavir escapees revealed 139 \nvarious mutations in the single-stranded annealing protein (beta) and exonuclease (exo) genes of 140 \nthe lambda Red operon, a known recombination system that also includes the RecBCD inhibitor 141 \ngam35–37 (Fig 2a). Many of these are early stop and frameshift mutations that result in loss-of-142 \nfunction of the proteins. Of the mutations that do not disrupt the open reading frame, the point 143 \nmutations in exo reside near its catalytic site and may influence its enzymatic activity (Fig S1a-c). 144 \nThe point mutation in beta resides in its C-terminal domain which interacts with exo38 (Fig S1e-f). 145 \nStevie_ev116 escapee mutations were primarily found in an exonuclease (exo) and recombinase 146 \n(rec) contained in the same operon (Fig 2b, S1d, S1g-h). Due to the presence of an ERF family 147 \ndomain on the recombinase, we predict that this operon is also a recombination system39.  148 \n 149 \nWe experimentally validated our candidate trigger genes from lambdavir by co-expressing them 150 \nwith Retron-Vpa2 in liquid cultures of E. coli bSLS.114 and measuring OD600 over 16 hrs. Co-151 \nexpression of the retron with the full lambda Red operon resulted in strong growth inhibition 152 \ncompared to a control strain lacking the retron, confirming that this operon is sufficient to trigger 153 \nthe abortive infection phenotype (Fig 2c-d, S2a). We also introduced two mutations that we 154 \nobserved in the escapees (beta T132fs and exo W163*) to the operon, which each resulted in partial 155 \nrescue of the growth suppression phenotype. Deleting the third component of the Red operon, gam, 156 \nwhich is a known trigger of Type IV and V retrons2,28, did not rescue growth. Given that the 157 \ncombination of beta and exo is sufficient to fully trigger the system, we next tested them separately, 158 \nfinding that neither affected cell growth when expressed individually (Fig 2e, S2b). Interestingly, 159 \ngam resulted in slight inhibition when expressed in the absence of beta and exo (Fig 2e, S2b). This 160 \ndiffers from the effect of gam on other clade 3 retrons, such as Type IV Retron-Eco6 and Type V 161 \nRetron-Sen1. For both these retrons, gam triggers them to the same degree as the full Red operon, 162 \nwith no requirement for beta or exo (Fig 2f-g, S2c). This points to a clear distinction between the 163 \ntriggering mechanism of Retron-Vpa2 and that of its close relatives.  164 \n 165  \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n 166 \nFigure 2 – Recombination systems trigger Retron-Vpa2 in two phages. 167 \nA) Operon of the genes mutated in lambdavir escapees of Retron-Vpa2. Specific mutations are labeled in red. Mutations used 168 \nin trigger assays are in bold. B) Operon of the genes mutated in Stevie_ev116 escapees of Retron-Vpa2. Specific mutations are 169 \nlabeled in red. C) Trigger assay showing growth curves of liquid cultures expressing Retron-Vpa2 with variants of the lambda 170 \nRed operon, relative to an empty vector with the Red operon, measured over 16 hrs. Each line is the mean of three biological 171 \nreplicates, error bands are plotted in Fig S2a. D) Comparison of OD600 measurements from panel (C) at the 10 hr timepoint, 172 \nshown as a percentage of the empty vector control at the same timepoint. Bars show the mean of three biological replicates 173 \nwith individual replicates plotted as white circles (one-way ANOVA P=0.0023, Dunnett’s test corrected for multiple 174 \ncomparisons versus empty vector+Red: Retron-Vpa2+Red P=0.0017, Retron-Vpa2+Beta+Exo P=0.0029, all other conditions 175 \nns). E) Trigger assay showing growth curves of liquid cultures expressing Retron-Vpa2 with Gam, Beta, or Exo, relative to an 176 \nempty vector with gam, measured over 16 hrs. Growth curve of Retron-Vpa2 with the Red operon from panel (C) shown in a 177 \ndotted red line for comparison. Each line is the mean of three biological replicates, error bands are plotted in Fig S2b. F) 178 \nTrigger assay showing growth curves of liquid cultures expressing either Retron-Eco6 or Retron-Sen1 with either the Red 179 \noperon or gam alone, relative to an empty vector with the Red operon, measured over 16 hrs. Each line is the mean of three or 180 \nmore biological replicates, error bands are plotted in Fig S2c. G) Comparison of OD600 measurements using data from panels 181 \n(C), (E) and (F) at the 10 hr timepoint, shown as a percentage of the empty vector control at the same timepoint. The Retron-182 \nVpa2 + Red and Retron-Vpa2 + Beta + Exo bars are repeated from (D) to permit direct comparison. Bars show the mean of 183 \nthree or more biological replicates with individual replicates plotted as white circles (one-way ANOVA P=0.0002, Šídák's test 184 \ncorrected for multiple comparisons versus empty vector+Red: Retron-Vpa2+Red P=0.0073, Retron-Vpa2+Beta+Exo P=0.0145, 185 \nRetron-Eco6+Red P=0.0119, Retron-Eco6+Gam P=0.0027, Retron-Sen1+Red P=0.0118, Retron-Sen1+Gam P=0.0298, all other 186 \nconditions ns; Šídák's test corrected for multiple comparisons Retron-Eco6+Red versus Retron-Eco6+Gam and Retron-187 \nSen1+Red versus Retron-Sen1+Gam both ns ). Additional statistical details in Supplementary Table 1. 188 \n 189 \nA hybrid RNA encodes a toxic effector and forms a complex with the RT and HTH 190 \nWe next investigated each component of the Retron-Vpa2 operon. First, we sought to 191 \nexperimentally validate the msr-msd region that we predicted based on its secondary structure15. 192 \nTranscripts from the ncRNA region of the retron operon are generally highly enriched in a cell2. We 193 \ntherefore performed RNAseq on Retron-Vpa2-expressing E. coli bSLS.114 and found that 194 \nsequencing coverage was enriched not only in the predicted msr-msd region, but also in the 195 \ndownstream SP coding sequence (Fig 3a). Performing the experiment with a catalytically dead RT 196 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n(YADD -> YAAA catalytic site mutation) had no effect on this enrichment (Fig S3). Because this RNA 197 \ncontains both the msr-msd and the effector coding region in a contiguous transcript, we decided to 198 \nname it the hybrid RNA (hyRNA).  199 \n 200 \nThe Vpa2 hyRNA is approximately 470 nt in length. Its predicted secondary structure, based on 201 \nminimal free energy, resembles three arms extending from a central node (Fig 3b). The first arm 202 \ncontains a highly structured region characteristic of a retron msr-msd7,15. The second arm contains 203 \ntwo stem-loops which incorporate the ribosomal binding site (RBS) and start codon of the SP. The 204 \nthird arm contains the remaining coding sequence of the SP. The 3’ end of the transcript then wraps 205 \nback around and anneals to the first arm via a set of inverted repeats, which resembles the a1/a2 206 \npriming region of a canonical retron9. A covariance model of the hyRNA independently supports the 207 \npredicted secondary structure (Fig S4a) and reveals multiple areas of high sequence conservation, 208 \nincluding the stems in the first arm, the inverted repeats that link each arm to the central node, the 209 \nSP RBS and start/stop codons, and a set of direct repeats in the loop regions of the first and second 210 \narm (Fig S4b).  211 \n 212 \nGiven the SP’s unique position within the hyRNA, we were curious about its role in this system. Its 213 \npredicted structure reveals a C-terminal domain with a hydrophobic core that is conserved across 214 \nmultiple Type VI homologs (Fig S5a-d). When expressed independently from a synthetic operon, we 215 \nobserved that the SP impairs cell growth at a similar level to triggering the retron with the lambda 216 \nRed operon (Fig 3c-d, S6a). Introduction of an early stop mutation (K15*) in the SP eliminated the 217 \nretron’s toxicity upon triggering. Additionally, introduction of three point mutations (I39R, F43R, 218 \nF45R) in the SP’s conserved hydrophobic core significantly reduced toxicity under the same 219 \nconditions. From these results, we conclude that the SP alone is the toxic effector of this retron and 220 \nthat its C-terminal domain is essential for toxicity.  221 \n 222 \nWe then proceeded to mutate other components of the system in the context of the full operon (Fig 223 \n3e-f, S6b). We found that deleting the noncoding region of the hyRNA from the 5’ end up until the 224 \nRBS (Δnc-hyRNA) was highly toxic. Similarly, mutating a conserved residue in the HTH protein 225 \n(R7D) was also toxic. However, mutating the RT catalytic domain (YADD -> YAAA)40 was not toxic, 226 \nindicating that reverse transcription is not necessary to neutralize the SP, contrasting with 227 \ncanonical retron mechanisms.  228 \n 229 \nTo understand which components of the system physically interact with each other, we next 230 \nengineered a version of Retron-Vpa2 with FLAG-tagged RT and HA-tagged HTH in a dSP (K15*) 231 \nbackground. We expressed the system in E. coli MG1655-DE3 and performed a pull-down on either 232 \nthe RT or HTH, followed by a Western blot and TBE-Urea PAGE to detect co-immunoprecipitation of 233 \nthe tagged proteins and hyRNA respectively. Pulling down on the RT led to co-immunoprecipitation 234 \nof the HTH and hyRNA (Fig 3g, S7a). Pulling down on the HTH led to co-immunoprecipitation of the 235 \nRT and hyRNA (Fig 3h, S7b). Thus, the RT, HTH, and hyRNA form a complex in vivo. We also checked 236 \nfor msDNA in the co-immunoprecipitation but did not detect any in either pull-down (Fig S7c). The 237 \nco-immunoprecipitated RNA from the RT pull-down was confirmed to be the hyRNA through 238 \nsequencing (Fig S7d). 239 \n 240  \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n 241 \nFigure 3 – A hybrid RNA encodes a toxic effector and forms a complex with the RT and HTH. 242 \nA) RNAseq coverage (normalized to maximum coverage) of Retron-Vpa2 operon. B) Predicted secondary structure of the 243 \nhyRNA with msr-msd, inverted repeats (IR), ribosome binding site (RBS), SP start codon (AUG) and SP coding region labeled. 244 \nC) Growth curves of liquid cultures expressing SP alone, Retron-Vpa2 triggered with the Red operon, and two Retron-Vpa2 SP 245 \nmutants triggered with the Red operon, relative to an empty vector with the Red operon, measured over 16 hrs. Cultures were 246 \ninduced at 2 hrs (vertical dotted line). Each line is the mean of three biological replicates, error bands are plotted in Fig S6a. 247 \nD) Comparison of OD600 measurements from panel (C) at the 10 hr timepoint, shown as a percentage of the empty vector 248 \ncontrol at the same timepoint. Bars show the mean of three biological replicates with individual replicates plotted as white 249 \ncircles (one-way ANOVA P<0.0001, Dunnett’s test corrected for multiple comparisons versus empty vector+Red: SP P<0.0001, 250 \nRetron-Vpa2+Red P<0.0001, Retron-Vpa2 SP(I39R, F43R, F45R) + Red P=0.0035, all other conditions ns). E) Growth curves of 251 \nliquid cultures expressing Retron-Vpa2 WT, Δnc-hyRNA mutant, HTH(R7D) mutant, and dRT(YADD->YAAA) mutant, relative 252 \nto an empty vector, measured over 16 hrs. Cultures were induced at 2 hrs (vertical dotted line). Each line is the mean of three 253 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nbiological replicates, error bands are plotted in Fig S6b. F) Comparison of OD600 measurements from panel (E) at the 10 hr 254 \ntimepoint, shown as a percentage of the empty vector control at the same timepoint. Bars show the mean of three biological 255 \nreplicates with individual replicates plotted as white circles (one-way ANOVA P<0.0001, Dunnett’s test corrected for multiple 256 \ncomparisons versus empty vector: Retron-Vpa2 Δnc-hyRNA P<0.0001, Retron-Vpa2 HTH(R7D) P<0.0001, all other conditions 257 \nns). G) Anti-FLAG pull-down of cultures expressing Retron-Vpa2, with and without a 3xFLAG-tagged RT. On the left are anti-258 \nHA and anti-FLAG Western blots of the pull-down and lysate fractions, with numbers showing protein weight in kDa. On the 259 \nright is a TBE-Urea PAGE analysis showing detection of hyRNA, with numbers showing length in nt. Complete gels are found 260 \nin Fig S7a. H) Anti-HA pull-down of cultures expressing Vpa2, with and without a HA-tagged HTH. On the left are anti-HA and 261 \nanti-FLAG Western blots of the pull-down and lysate fractions, with numbers showing protein weight in kDa. On the right is a 262 \nTBE-Urea PAGE analysis showing detection of hyRNA, with numbers showing length in nt. Complete gels are found in Fig S7b. 263 \nAdditional statistical details in Supplementary Table 1. 264 \n 265 \nmsDNA production is repressed until Retron-Vpa2 is triggered 266 \nGiven that we detected no abundant msDNA and the catalytic activity of the RT is not necessary for 267 \ntoxin neutralization, we wondered whether the RT could play a non-catalytic role in this retron. We 268 \ntested phage defense against lambdavir with a version of Retron-Vpa2 carrying a catalytically dead 269 \nRT and found that inactivating reverse transcription eliminated phage defense (Fig 4a). Given the 270 \nnecessity of reverse transcription for defense, the absence of msDNA was even more puzzling. Thus, 271 \nwe decided to implement a new strategy for detecting msDNA, inspired by our previous work 272 \ndeveloping a CRISPR-based molecular recorder33. We call this strategy Spacer-seq (Fig 4b). Spacer-273 \nseq harnesses the E. coli Type I-E CRISPR-Cas adaptation machinery, which mediates capture of 274 \nforeign DNA elements into CRISPR arrays41–44. This will “record” the presence of retron msDNA into 275 \nthe host cell’s genome over time. In this CRISPR-Cas system, captured DNA is recorded as 33 nt 276 \nspacers via the integrase proteins Cas1 and Cas243,44.  277 \n 278 \nWe co-expressed Retron-Vpa2 with Cas1-Cas2 in E. coli bSLS.114 (which contains an endogenous 279 \nCRISPR array but no endogenous Cas proteins). After 24 hrs of induction in liquid culture, we 280 \nharvested cells and PCR amplified the CRISPR locus. We sequenced these amplicons, identified 281 \nnewly acquired spacer sequences, and mapped them to the Retron-Vpa2 operon. From this, we 282 \nobserved a 63 nt region of enriched coverage (normalized to total spacers acquired) that is absent 283 \nin a dRT control (Fig 4c). Encouragingly, this region covers a long stem loop in the predicted msr-284 \nmsd portion of the hyRNA. We further validated this method by performing the same experiment on 285 \nRetron-Eco6 and Retron-Sen1, where we similarly saw enriched coverage of their known msd 286 \nregions (Fig S8a-b). No other peaks of this magnitude were found when mapping coverage across 287 \nthe entire plasmid or E. coli genome (Fig S8c). Thus, there is some ongoing msDNA production from 288 \nRetron-Vpa2, which does not accumulate to a level we can visualize on a polyacrylamide gel. 289 \n 290  \nTo identify other factors that might influence msDNA production, we proceeded to perform Spacer-291 \nSeq on dSP (K15*) and dHTH mutants of Retron-Vpa2 (Fig 4d). While mutating the SP did not 292 \nchange the distribution of coverage relative to wild-type, mutating the HTH (in a dSP background to 293 \navoid toxicity) enhanced coverage by nearly five-fold. We then looked at msDNA coverage of Vpa2 294 \ndSP in the presence of different components of the lambda Red operon (Fig 4e-h). Notably, we 295 \nobserved a rightward extension of the spacer distribution peak when Gam is expressed (Fig 4e, h), 296 \nmeaning that this version of the msDNA contains an additional 23 nt at its 5’ end, totaling to a 297 \nlength of 86 nt. This extension also occurs in a ΔrecB background, indicating that it is linked to 298 \nGam’s biological function of RecBCD inhibition36,37 (Fig S7a).  299 \n 300 \nGiven the increase in Spacer-seq coverage yielded by the dHTH mutation, we attempted once again 301 \nto visualize msDNA on a gel under similar conditions. Expressing Retron-Vpa2 with a dHTH 302 \nmutation revealed a series of bands between 60 nt and 90 nt, which are not present in a dRT 303 \ncontrol or the background dSP strain (Fig 4i, S7b). Additionally, expressing Retron-Vpa2 dSP in the 304 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\npresence of the lambda Red operon yielded a single band between 80 nt and 90 nt, which increases 305 \nin intensity with a dHTH mutation. Thus, the HTH appears to repress reverse transcription. 306 \nSequencing the product with the single band confirmed that it is the same 86 nt product identified 307 \nin Spacer-seq when Retron-Vpa2 was expressed in the presence of Gam (Fig 4i). We consider this 308 \nsequence as the full-length msDNA and refer to the 63 nt sequence observed in the absence of Gam 309 \nas the truncated msDNA.  310 \n 311 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n 312 \nFigure 4 – msDNA production is repressed until Retron-Vpa2 is triggered. 313 \nA) Spot assay showing titration of phage lambdavir on Retron-Vpa2 wild-type vs. dRT expressed in E. coli bSLS.114, relative 314 \nto an empty vector control. Efficiency of plating is quantified as pfu/mL in the adjacent bar graph. B) Schematic of Spacer-315 \nseq. C) Per base coverage (normalized to total spacers acquired) of spacers aligning to the Retron-Vpa2 operon for the wild-316 \ntype retron compared to a dRT control. Position of aligned spacers is also depicted as raw per base coverage on the predicted 317 \nsecondary structure of the msr-msd region of the hyRNA. Data is a pool of three biological replicates. D) Normalized per base 318 \ncoverage of spacers aligning to the msr-msd region of Retron-Vpa2 for various mutants, compared to the wild-type retron 319 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nfrom panel (C). Data is a pool of three biological replicates. E-H) Per base coverage (normalized to maximum coverage) of 320 \nspacers aligning to the msr-msd region of Retron-Vpa2 for a Retron-Vpa2 dSP mutant with Gam, Beta, Exo, or the Red 321 \noperon, compared to Retron-Vpa2 dSP alone from panel (D). Data is a pool of three biological replicates. I) PAGE analysis for 322 \nmsDNA production of Retron-Vpa2 dSP with various mutants, as well as with co-expression of the Red operon. Retron-Eco6 323 \nband shown as reference. Ladder on the leftmost lane is single-stranded DNA with increments marked in nucleotides. 324 \nComplete gel is found in Fig S7b. Position of aligned reads from sequencing the Retron-Vpa2 dSP + Red sample is depicted as 325 \nraw per base coverage on the predicted structure of the hyRNA msr-msd. 326 \n 327 \nSP is translated in the presence of triggers and msDNA 328 \nThe final step of Retron-Vpa2 phage defense is release of the toxic SP to arrest growth in the 329 \ninfected host. We next aimed to investigate the mechanism of SP neutralization/activation. Because 330 \nthe RBS and start codon of the SP are sequestered within stem loops of the hyRNA, we hypothesized 331 \nthat the toxin is translationally repressed until the retron is triggered. To measure its production in 332 \nvivo, we designed a split GFP complementation assay45 where SP abundance is detected through a 333 \nfluorescence-based readout. In this assay, an optimized superfolded GFP (sfGFP) is split into two 334 \nfragments which must self-assemble to reconstitute a fluorescent protein. The large fragment 335 \ncontains beta sheets 1-10 (sfGFP1-10), which we express from an independent plasmid. The small 336 \nfragment contains beta sheet 11 (sfGFP11), which we fuse to the C-terminus of a truncated SP 337 \n(containing only the first 31 residues) within the Retron-Vpa2 operon. This way, the total length of 338 \nthe hyRNA is preserved and the SP is rendered nontoxic (Fig S10). In E. coli bSLS.114, we co-339 \nexpressed sfGFP1-10 with SP::sfGFP11 in a Δnc-hyRNA mutant of Vpa2 (which we previously saw 340 \nwas toxic with a wild-type SP, thus we expect SP::sfGFP11 to be produced). However, we did not 341 \nobserve fluorescence signal in this condition, nor in a condition where sfGFP11 is not fused to the 342 \nSP (Fig 5a). We decided to redesign our constructs, this time adding the first 31 residues of the SP 343 \nto the N-terminus of the sfGFP1-10 as well, with the rationale being that if there existed natural 344 \ninteractions between SP monomers, it would improve the chance of self-assembly. Tagging both 345 \nsfGFP fragments with part of the SP resulted in a fluorescence signal well over 10000-fold higher 346 \nthan a negative control without sfGFP11 (Fig 5b).  347 \n 348 \nWe used this assay to measure SP production in the context of a Retron-Vpa2 dHTH mutant and 349 \nobserved high fluorescence (Fig 5c), indicating that the SP is translated when the HTH is 350 \ninactivated. Since this HTH mutant was previously shown to result in the accumulation of 351 \ndetectable msDNA, we next checked whether SP production in the dHTH background requires 352 \nmsDNA by mutating the RT so that no msDNA can be produced. We found that mutating the RT 353 \ndramatically reduced fluorescence under these conditions, demonstrating that msDNA is necessary 354 \nfor translation of the SP. Meanwhile, co-expression of the system with variants of the lambda Red 355 \noperon only resulted in fluorescence in the Red Beta + Exo and wild-type Red conditions, 356 \nrecapitulating our earlier trigger assay results which demonstrated the necessity of beta and exo for 357 \ncomplete retron activation (Fig 5d-e). Interestingly, the maximum fluorescence achieved when 358 \ntriggering the retron with Red was roughly 10-fold less than in the Δnc-hyRNA and dHTH 359 \nconditions, but still 1000-fold higher than a control where the system is co-expressed with an 360 \nempty vector. Thus, the SP is translationally repressed, but this repression is alleviated when co-361 \nexpressing phage triggers, removing the nc-hyRNA, or mutating the HTH, as long as the RT is intact 362 \nand msDNA can accumulate. 363 \n 364  \nSo far, we have observed parallels between the conditions that lead to msDNA production and those 365 \nthat lead to SP production. We wondered whether the two processes are directly linked, where the 366 \nmsDNA molecule itself is derepressing the translation of the SP. To test this, we co-expressed our 367 \nsplit GFP system with an engineered version of Retron-Eco1 that produces either a scrambled 368 \nmDNA sequence, the truncated version of Retron-Vpa2’s msDNA, or the full Retron-Vpa2 msDNA. 369 \nWe observed clear upregulation of SP production in the condition with full Retron-Vpa2 msDNA 370 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\ncompared to the scrambled or truncated version (Fig 5f). This shows that abundant Retron-Vpa2 371 \nmsDNA produced in trans from a different retron is sufficient to induce translation of the SP, and 372 \nthat the stem-loop at the 5’ end of the msDNA is essential. 373 \n 374 \n 375 \nFigure 5 – SP is translated in the presence of triggers and msDNA. 376 \nA) Measurement of relative fluorescence units (RFU) for E. coli bSLS.114 liquid cultures expressing a split GFP system over 16 377 \nhrs, with induction at 0 hrs. No fluorescence is observed for sfGFP1-10 co-expressed with either sfGFP11 or SP::sfGFP11 (in 378 \nthe context of Retron-Vpa2 Δnc-hyRNA), compared to a control without sfGFP11. Alphafold3-predicted structures for sfGFP1-379 \n10 and SP::sfGFP11 are also depicted. B) Measurement of RFU for E. coli bSLS.114 liquid cultures expressing a split GFP 380 \nsystem over 16 hrs, with induction at 0 hrs. Solid lines indicate the mean and error bands indicate the standard deviation of 381 \nthree biological replicates. Strong fluorescence is observed for SP::sfGFP1-10 co-expressed with SP::sfGFP11 (in the context of 382 \nRetron-Vpa2 Δnc-hyRNA), compared to the systems tested in panel (A). Alphafold3-predicted protein structure shows 383 \nreconstitution of SP::sfGFP1-10 and SP::sfGFP11. C) Measurement of RFU for E. coli bSLS.114 liquid cultures expressing a split 384 \nGFP system in the context of Retron-Vpa2 Δnc-hyRNA or dHTH mutants over 18 hrs, with induction at 2 hrs (vertical dotted 385 \nline), compared to a control system lacking sfGFP11. Solid lines indicate the mean and error bands indicate the standard 386 \ndeviation of three biological replicates. D) Measurement of RFU for E. coli bSLS.114 liquid cultures expressing a split GFP 387 \nsystem in the context of Retron-Vpa2 over 18 hrs, with induction at 2 hrs (vertical dotted line), co-expressed with either Gam, 388 \nBeta, or Exo, compared to co-expression with an empty vector. Solid lines indicate the mean and error bands indicate the 389 \nstandard deviation of three biological replicates. E) Measurement of RFU for E. coli bSLS.114 liquid cultures expressing a split 390 \nGFP system in the context of Retron-Vpa2 over 18 hrs, with induction at 2 hrs (vertical dotted line), co-expressed with 391 \nvariants of the Red operon, compared to co-expression with an empty vector. Solid lines indicate the mean and error bands 392 \nindicate the standard deviation of three biological replicates. F) Measurement of RFU for E. coli bSLS.114 liquid cultures 393 \nexpressing a split GFP system in the context of Retron-Vpa2 over 18 hrs, with induction at 2 hrs (vertical dotted line), co-394 \nexpressed with Retron-Eco1 engineered to produce either scrambled msDNA, truncated Retron-Vpa2 msDNA, or full-length 395 \nRetron-Vpa2 msDNA. Solid lines indicate the mean and error bands indicate the standard deviation of three biological 396 \nreplicates.   397 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nDISCUSSION 398 \nIn this work, we uncovered a new mechanism of phage defense by a Type VI retron from Vibrio 399 \nparahaemolyticus, which we named Retron-Vpa2 (Fig 6). This system defends against a broad range 400 \nof phages, from which we characterized escapees of two phages and identified their recombination-401 \nassociated trigger genes. This retron, and likely other members of its type, differs from others 402 \nthrough its inverted triggering mechanism, where the msDNA acts as the agent for toxin 403 \nproliferation, rather than inhibition. We determined that Retron-Vpa2 msDNA, which is 404 \nundetectable at baseline, accumulates 1) in the presence of phage infection, 2) with the expression 405 \nof phage triggers, or 3) if the HTH protein is mutated. The accumulation of msDNA is necessary for 406 \ntranslation of the toxic SP, which is transcribed as part of the hyRNA, but translationally repressed 407 \nin the absence of msDNA. Finally, the msDNA activates translation of the SP even if produced from 408 \nanother retron in trans within the same cell. 409 \n 410 \n 411 \nFigure 6 – Model of Retron-Vpa2 phage defense mechanism 412 \nThe core of this system is a complex comprising the RT, HTH, and hyRNA. The hyRNA contains the template for reverse 413 \ntranscription and the coding sequence for SP translation. In its untriggered state, reverse transcription is inhibited by the 414 \nHTH and the SP is translationally repressed. Phage infection and expression of trigger proteins, such as the Red system of 415 \nlambdavir, induce accumulation of msDNA, which leads to translation of toxic SP. 416 \n 417 \nGiven the complexity of this mechanism, there remain several aspects of it that will prompt future 418 \ninvestigation. For example, we identify in this study that recombination-associated genes in phage 419 \nlambda and phage Stevie_ev116 trigger Retron-Vpa2. Still, it is uncertain exactly how this occurs. 420 \nWe hypothesize that the triggers stabilize the msDNA, such that accumulation of msDNA leads to 421 \ntranslation of SP. This is supported by the fact that lambda beta and Stevie_ev116 rec are single-422 \nstranded annealing proteins (SSAPs) that bind ssDNA, while exo also interacts with DNA35,39. 423 \nFurthermore, it is encouraging that we can identify homologs of these proteins in the majority of 424 \nphages that Retron-Vpa2 defends against. In the case of lambda phage, the presence of gam would 425 \nalso promote msDNA accumulation by inhibiting E. coli RecBCD. We think RecBCD, and potentially 426 \nprophage-associated nucleases, mediate turnover of msDNA from basal levels of reverse 427 \ntranscription. It is also possible that repression of reverse transcription is alleviated by interactions 428 \nof phage triggers with the HTH.  429 \n 430 \nThrough our efforts to understand the role of reverse transcription in this system, we establish 431 \nSpacer-seq as a discovery-based tool for identifying retron msDNA. Here, it provided evidence that 432 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nreverse transcription does occur in the absence of triggers, but due to reduced efficiency and/or 433 \nhigh turnover rate, msDNA does not accumulate in the cell, and thus cannot be detected through 434 \nPAGE analysis. We speculate that much of the msDNA captured by the integrases are degradation 435 \nproducts of RecBCD and host nucleases, a known source of prespacers for CRISPR adaptation46,47. 436 \nFurthermore, the appearance of the 23 nt 5’ region in the msDNA when RecBCD is inhibited could 437 \nindicate that the site is normally shielded from the CRISPR integrases by this complex or some 438 \nother protein, or that it is not synthesized at all. 439 \n 440 \nPerhaps the most intriguing aspect of this retron is its repression of the SP, which requires msDNA 441 \nproduction for translation. We think the underlying mechanism involves a type of riboregulator. 442 \nNotably, the SP RBS and start codon are located within stem loops in the hyRNA, which would 443 \nsequester them from ribosomal access. In addition, we notice a set of 7 nt direct repeats in the 444 \nhyRNA that are highly conserved in sequence and structure across the Type VI retrons (Fig S4a-b). 445 \nOne repeat is in a loop region of the msd and the other is in a loop region of the predicted 446 \nriboregulator. When reverse transcribed, the repeat in the msd would be complementary to that in 447 \nthe riboregulator, potentially leading to hybridization and subsequent release of the SP RBS and 448 \nstart codon. Importantly, this 7 nt sequence is on the 5’ end of the msDNA, which we demonstrated 449 \nis specifically necessary for SP translation (Fig 5f). Future work will seek to further understand the 450 \nspecific sequence and structural elements necessary for the msDNA to activate SP translation.  451 \n 452 \nFinally, although the mechanism of SP toxicity remains elusive, we have made some important 453 \ninsights into its biology. We demonstrated that maintaining hydrophobicity of its C-terminal alpha 454 \nhelix is necessary for toxicity (Fig 3b). Furthermore, our results showing augmented fluorescence of 455 \nsplit GFP when both fragments are tagged with the first 31 residues of the SP suggest that the N-456 \nterminus of the protein mediates oligomerization or colocalization (Fig 5b).  457 \n 458 \nOverall, our study into the biology and mechanism of Retron-Vpa2 reveals a novel retron system 459 \nwhere the toxic effector is translationally repressed within a hyRNA transcript and requires phage-460 \ninduced accumulation of msDNA for expression. These insights will inform future studies on the 461 \npreviously uncharacterized Type VI group of retrons while introducing a new perspective in our 462 \nunderstanding of Retron-mediated phage defense.  463 \n 464 \n  465 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nMETHODS 466 \n 467 \nBacterial strains and growth conditions 468 \nThe E. coli strains used in this study were NEB 5-alpha (NEB C2987) for cloning, MG1655 for the large panel 469 \nphage defense assays, MG1655-DE3 for the pull-down assays (see Pull-down assays for details on strain 470 \nconstruction), and bSLS.114 (derivative of BL21-AI lacking Retron-Eco1)48 for all other experiments. A ΔrecB 471 \nderivative of bSLS.114 (bKAZ020) was constructed using lambda Red recombineering49 for the experiment 472 \nshown in Supplementary Figure 9a.  473 \n 474 \nFor experiments with liquid cultures, bacteria were grown in lysogeny broth (LB) at 37°C, shaking at 250 475 \nrpm. Where appropriate, the LB was supplemented with antibiotics and inducers at the following working 476 \nconcentrations: 200 μg/mL ampicillin, 100 μg/mL carbenicillin (GoldBio C-103), 25 μg/mL chloramphenicol 477 \n(GoldBio C-105), 35 μg/mL kanamycin (GoldBio K-120), 25 μg/mL spectinomycin (GoldBio S-140), 1 mM 478 \nisopropyl β-d-1-thiogalactopyranoside (IPTG, GoldBio I2481C), and 2 mg/mL L-arabinose (GoldBio A-300). 479 \n 480 \nPhage strains and propagation 481 \nExperiments with phage lambda used a strictly lytic strain (lambdavir) that was generously provided by 482 \nLuciano Marraffini. To propagate the phage, a saturated culture of E. coli bSLS.114 was diluted 1:100 in 3 mL 483 \nLB supplemented with 0.1 mM MnCl2 and 5 mM MgCl2 (MMB) and grown to an OD600 of 0.25. At this point, the 484 \nculture was infected with phage at a MOI of ~0.2 and allowed to grow for an additional 16 hrs. After this time, 485 \nthe culture was centrifuged for 10 min at 3,434g and the supernatant containing the phage lysate was filtered 486 \nthrough a 0.2-μm filter. The titer of the lysate was determined via plaque assay (see Plaque assays). Lysates 487 \nwere stored at 4°C for further use. 488 \n 489 \nFor other phages used in the large panel defense screen, phages were isolated by using sterile inoculation 490 \nloop to streak from the original stock (obtained from various sources, see Supplementary Table 3) onto a 491 \nlawn of E. coli MG1655 on LB agar. The plate was incubated overnight at 37°C. Then, singe plaques were 492 \npicked and used to infect a 3 mL LB culture of E. coli MG1655 at an OD600 of ~0.2. The culture was then grown 493 \nfor ~5 hrs at 37°C. After this time, the culture was clarified via centrifugation and the supernatant containing 494 \nthe phage lysate was filtered through a 0.45 μm filter. 495 \n 496 \nPlasmid construction 497 \nPlasmids containing Retron-Vpa2 (pMRM27) and Retron-Eco12 (pMRM38) were designed to incorporate the 498 \nwild-type operon (including the native promoter) into a high-copy pET-21 backbone under a T7/lac inducible 499 \npromoter. These constructs were synthesized by Twist Bioscience. Plasmids containing other retrons tested 500 \nfor msDNA expression were taken from the previous retron census manuscript15.  501 \n 502 \nFor plasmids containing lambdavir trigger genes, the wild-type trigger genes were PCR amplified from phage 503 \nlysate and cloned using Gibson assembly (NEB E2621) into a pBR322 backbone under an arabinose-inducible 504 \npromoter (araBAD). Plasmids containing mutations to the trigger genes were cloned via Q5 site-directed 505 \nmutagenesis (NEB E0552).  506 \n 507 \nPlasmids containing mutations to Retron-Vpa2 (e.g. Δnc-hyRNA, dRT, dHTH, dSP, dSP) or introducing protein 508 \ntags to retron components (e.g. HTH-HA, RT-FLAG) were similarly cloned using site-directed mutagenesis.  509 \n 510 \nFor Spacer-seq, the plasmid expressing the Cas1-Cas2 genes under a lac promoter (pSCL565) in a pCDF 511 \nbackbone was taken directly from a previous manuscript50.  512 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n 513 \nFor the split GFP complementation assay, optimized sfGFP1-10 (sequence taken from pAGM22082 in 514 \nPüllmann et al.51, Addgene #153515) was synthesized via Integrated DNA Technologies (IDT) and cloned into 515 \na p15a backbone under a lac promoter using Gibson assembly. Site-directed mutagenesis was used to add the 516 \nfirst 31 amino acids of the SP to the N-terminus of this plasmid (SP-sfGFP1-10). Site-directed mutagenesis 517 \nwas also used to fuse sfGFP11 (5’-RDHMVLHEYVNAAGIT-3’) with the SP gene (SP::sfGFP11) on a Retron-518 \nVpa2-expressing plasmid. To clone the plasmid used for Retron-Eco1 expression of scrambled msDNA, the 519 \nretron operon under a T7/lac inducible promoter was amplified from a Retron-Eco1-based editor plasmid 520 \n(operon is rearranged with the RT preceding the ncRNA and the effector deleted) from a previous 521 \nmanuscript52 and cloned into a pCDF backbone using Gibson assembly. Site-directed mutagenesis was used 522 \non this plasmid to clone versions that expressed either the truncated or full version of the Retron-Vpa2 523 \nmsDNA. 524 \n 525 \nAll oligos used in this study were synthesized by IDT. Details for all plasmids are listed in Supplementary 526 \nTable 4. 527 \n 528 \nPhylogenetic analysis 529 \nHomologs of Retron-Vpa2 were identified using the RT amino acid sequence (WP_069548583.1) as a BLASTP 530 \nsearch query of the NR protein database (max target sequences = 5000). The resulting sequences were 531 \naligned with MAFFT53, and a Maximum-Likelihood phylogenetic reconstruction was made with 532 \nVeryFastTree54 with default parameters. Classification of the different retrons was done by annotating the 533 \ngenomic neighborhoods with PADLOC55 and used to color the tree.  534 \n 535 \nmsDNA expression and gel analysis 536 \nPlasmids containing retrons were transformed into bSLS.114. Individual colonies were picked and grown in 3 537 \nmL LB + carbenicillin at 37°C until saturation. Each culture was diluted 1:100 in 25 mL LB + carbenicillin, 538 \ngrown for ~2 hrs at 37°C (until OD600 reaches ~0.5), then induced with IPTG and arabinose, after which it was 539 \ngrown for an additional 5 hrs.  540 \n 541 \nDNA was isolated from the cells using the Qiagen Plasmid Plus Midi Kit (Qiagen #12945) and eluted in 150 μL 542 \nof molecular biology grade water. The eluted DNA was mixed with Novex TBE-Urea Sample Buffer (Invitrogen 543 \nLC6876) to reach a sample buffer concentration of 1X. This was then heated to 98 C for >=5 min and 15 µL 544 \nwas loaded onto a Novex 15% TBE-Urea gel (Invitrogen EC6885). ss20 DNA Ladder (Simplex Sciences) was 545 \nalso loaded in a separate well as a marker. The gel was run at 200 V for 45 min in preheated (>75°C) TBE 546 \nrunning buffer, then stained with SYBR Gold Nucleic Acid Gel Stain (Invitrogen S11494) and imaged on a 547 \nChemiDoc Imaging System (Bio-Rad).  548 \n 549 \nPlaque assays 550 \nSmall-drop plaque assays were performed similarly to Mazzocco et al.56, which was used to both titer phage 551 \nstocks and test for defense. For titering, plaque assays were performed with E. coli bSLS.114. For defense 552 \nexperiments with lambdavir, a plasmid containing Retron-Vpa2 (or an empty vector control) was first 553 \ntransformed into E. coli bSLS.114. Individual colonies were picked and grown in 3 mL LB + carbenicillin until 554 \nsaturation. This culture was then diluted 1:100 in MMB + carbenicillin and grown for 2-3 hrs at 37°C. For each 555 \nplaque assay performed, 200 μL of the passaged culture was mixed with 2 mL of MMB top agar (MMB with 556 \n0.75% agar) and poured onto a single well of a rectangular 4-well MMB agar plate (MMB with 1.5% agar). 557 \nAfter the top agar has solidified, tenfold serial dilutions in MMB of a phage lysate was spotted onto the plate 558 \nwith 2 μL spots in technical triplicates. After the spots had completely dried, plates placed in an incubator at 559 \n37°C overnight. Plaque-forming units (pfus) were quantified using the formula: pfu count * dilution factor / 560 \nmL of lysate spotted = pfu/mL. 561 \n 562 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nFor the large panel defense assay, the protocol is similar to above, but plasmids were instead transformed 563 \ninto E. coli MG1655 and colonies were grown in LB + ampicillin. Top agar was spiked with 2% of bacterial 564 \nculture directly from a saturated overnight and phage lysate was spotted onto the plate with 5 μL spots in 565 \ntechnical triplicates.  566 \n 567 \nPhage infection growth curves 568 \nA plasmid containing Retron-Vpa2 (or an empty vector control) was transformed into E. coli MG1655. 569 \nIndividual colonies were picked and used to start an overnight culture in LB + ampicillin. This saturated 570 \nculture was then diluted 1:100 into fresh LB + ampicillin and grown to an OD600 of ~0.2. 180 μL of this culture 571 \nwas transferred to a 96-well microtiter plate and infected with phage. All wells received 20 μL of phage lysate 572 \nthat was pre-diluted to various concentrations to achieve final multiplicities of infection (MOIs) of 0, 0.01, 0.1, 573 \n1, and 10. The plate was incubated in a plate reader (BioTek Synergy H1 Multimode Reader) at 37°C with 574 \nshaking at 200 rpm  for 6 hrs, with OD600 measurements taken every 5 min. 575 \n 576 \nPhage escapee isolation and sequencing 577 \n 578 \nPhage lambdavir 579 \n 580 \nIsolation of lambdavir escapees generally follows the protocol in Millman et al.2. 20 uL of lambdavir lysate (at a 581 \ntiter of ~108 pfu/mL) was used to infect 200 µL of Retron-Vpa2-expressing E. coli bSLS.114 culture at an 582 \nOD600 of 0.3, then incubated for 15 min at room temperature. The infected culture was then mixed into 2 mL 583 \nof MMB top agar, poured onto a single well of a rectangular 4-well MMB agar plate, and grown overnight at 584 \n37°C. The next morning, individual plaques were picked and resuspended in 90 μL of phage buffer (50 mM 585 \nTris pH 7.4, 100 mM MgCl2, 10 mM NaCl). Plaques were left in the buffer for 1 hr at room temperature, with 586 \noccasional vortexing to release the phage from the agar. This lysate was used to infect a 3 mL MMB + 587 \ncarbenicillin culture of Retron-Vpa2-expressing E. coli bSLS.114 at an OD600 of ~0.2, which was grown 588 \novernight at 37°C. The next morning, the culture was centrifuged at 4347g for 10 min and filtered through a 589 \n0.2-μm filter. Phage genomic DNA was isolated from 1 mL of this lysate using Norgen’s Phage DNA Isolation 590 \nKit (Norgen Biotek 46800).  591 \n 592 \nThe purified DNA was prepped for nanopore sequencing on a MinION device using the standard Oxford 593 \nNanopore Technologies (ONT) workflow for the Ligation Sequencing Kit (SQK-LSK109) with the Native 594 \nBarcoding Expansion (EXP-NBD196) on flow cell version R9.4.1 (FLO-MIN106D). The ancestral phage strain 595 \nwas also sequenced alongside the escapees. Sequencing data was base called and aligned to a reference 596 \ngenome (accession J02459.1) using Guppy (v6.1.7 from ONT). Aligned reads were used to generate a 597 \nconsensus sequence for each phage using SAMtools57 (v1.18). Escapee-specific mutations were identified 598 \nthrough alignment of escapee phage genomes with the ancestral phage genome using Geneious.  599 \n 600 \nPhage Stevie_ev116 601 \n 602 \nTo isolate Stevie_ev116 escapees, phage lysate was streaked onto a top agar lawn of 2% Retron-Vpa2-603 \nexpressing E. coli MG1655 supplemented with ampicillin. After overnight incubation at 37°C, resulting 604 \nplaques were picked and used to infect a LB + ampicillin culture of Retron-Vpa2-expressing E. coli MG1655 at 605 \nan OD600 of ~0.2. The culture was incubated overnight at 37°C to propagate the phage, and phage lysates 606 \nwere harvested the next morning. This process was repeated for three additional rounds to enrich the phage. 607 \n 608 \nPhage DNA was purified from filtered lysates using a combination of enzymatic digestion and silica column-609 \nbased extraction. To degrade contaminating bacterial nucleic acids, 450 μL of filtered phage lysate was 610 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nincubated with 50 μL of DNase I 10× buffer, 1 μL of DNase I (1 U/μL), and 1 μL of RNase A. Samples were 611 \nincubated at 37 °C for 1.5 hours without shaking. Enzymes were subsequently inactivated by adding 20 μL of 612 \n0.5 M EDTA. To digest remaining proteins, 1.25 μL of proteinase K was added, and samples were incubated at 613 \n56 °C for 1.5 hours without shaking. DNA was extracted using the DNeasy Blood & Tissue Kit (QIAGEN) 614 \nfollowing a modified protocol. 450 μL of Buffer AL was added to each sample and mixed thoroughly by 615 \nvortexing, followed by incubation at 56 °C for 10 minutes. Next, 450 μL of 96-100% ethanol was added, and 616 \nthe mixture was vortexed to ensure homogeneity. The entire sample was then transferred to a DNeasy Mini 617 \nspin column placed in a 2 mL collection tube and centrifuged at approximately 8,000 rpm for 1 minute. The 618 \nspin column was placed in a new collection tube, 500 μL of Buffer AW1 was added, and centrifugation was 619 \nrepeated at 8,000 rpm for 1 minute. With the column in a new collection tube, another 500 μL of Buffer AW2 620 \nwas added and samples were centrifuged at 8,000 rpm for 3 minutes. The spin column was then transferred 621 \nto a clean microcentrifuge tube for DNA elution. DNA was eluted by adding 30 μL of nuclease-free water 622 \ndirectly to the center of the membrane, incubating at room temperature for 1 minute, and centrifuging at 623 \n8,000 rpm for 1 minute. A second elution was performed by adding an additional 20 μL of nuclease-free 624 \nwater, incubating for 1 minute, and centrifuging as before. Purified DNA was stored at 4 °C until use. 625 \n 626 \nThe extracted DNA was quantified using the Qubit dsDNA High Sensitivity Assay Kit (Thermo Fisher 627 \nScientific). Sequencing libraries were prepared using the Nextera XT DNA Library Preparation Kit (Illumina), 628 \nfollowing the manufacturer's protocol with minor modifications. 1 ng of input DNA per sample was subjected 629 \nto enzymatic tagmentation, which simultaneously fragmented the DNA and tagged it with adapter sequences. 630 \nTagmented DNA was then amplified via limited-cycle PCR using Nextera XT barcoded index primers. Post-631 \nPCR clean-up was performed using AMPure XP magnetic beads (Beckman Coulter). Final library 632 \nconcentrations were measured using Qubit and normalized to equimolar concentrations before pooling. 633 \nLibraries were submitted to Novogene (Beijing, China) for high-throughput sequencing. Sequencing was 634 \nperformed on an Illumina platform using paired-end 2×150 bp chemistry. Adapter trimming and removal of 635 \nlow-quality base calls were performed by Novogene using their in-house quality control pipeline, and 636 \ndemultiplexed reads were returned in FASTQ format for downstream analysis. 637 \n 638 \nTrigger assays 639 \nA plasmid containing the retron being tested (or an empty vector control) was co-transformed with a plasmid 640 \ncontaining the trigger gene into E. coli bSLS.114. Individual colonies were picked and grown until saturation in 641 \n3 mL LB + carbenicillin + chloramphenicol at 37°C. Cultures were then diluted 1:100 in 200 uL of fresh LB + 642 \ncarbenicillin + chloramphenicol in a transparent 96-well microtiter plate ( Corning 3596). The plate was 643 \ncovered with an air-permeable seal (Breathe-Easy, Diversified Biotech BEM-1) and incubated in a microplate 644 \nreader (Molecular Devices SpectraMax i3) with shaking at 37°C for 16 hrs with OD600 measured every 20 min. 645 \nAt the 2 hr timepoint, IPTG and arabinose inducers were added to all cultures. Raw OD600 values were adjusted 646 \nby first multiplying by the empirically determined correction factor of 3.4967 , then subtracting the 0 hr 647 \nmeasurement for each sample.  648 \n 649 \nRNA-seq 650 \nRNA-seq was performed similar to described in Millman et al.2, here using a Retron-Vpa2-expressing strain of 651 \nE. coli bSLS.114. A saturated culture of this strain was diluted 1:100 in 5 mL LB + carbenicillin and grown at 652 \n37°C to an OD600 of 0.6. The culture was centrifuged at 4347g for 10 min at 4°C. The supernatant was discarded, 653 \nand the pellet was retained. The pellet was then treated with 100 μL of 2 mg/mL lysozyme solution (in 10 mM 654 \nTris, 1 mM EDTA, pH 8.0) and incubated for 15 min at 37°C in a 1.5 mL tube. At this time, 1 mL of TRI-reagent 655 \n(Sigma-Aldrich 93289) was added and mixed by pipetting, then 200 μL of chloroform was added and mixed by 656 \npipetting. The sample was incubated at room temperature for 5 min for phase separation. The sample was then 657 \ncentrifuged at 18,213g for 30 min at 4°C, after which the upper phase was carefully removed and transferred 658 \nto a fresh tube. This was mixed with 750 μL of cold isopropanol, sodium acetate (NaAOc) to a final concentration 659 \nof 0.3 M, and linear acrylamide (Invitrogen AM9520) to a final concentration of 10 μg/mL. The sample was 660 \nfrozen either at - 20°C overnight or - 80°C for 1 hr , then centrifuged at 18,213g for 30 min at 4 °C. The 661 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nsupernatant was discarded, and the pellet was washed twice with 500 μL of freshly prepared, cold 70% ethanol. 662 \nFor each wash, the ethanol  was carefully added without disturbing the pellet,  then centrifuged for 15 min at 663 \n18,213g at 4°C, and then the supernatant was discarded. Pellet was air-dried for 5 min, then resuspended in 50 664 \nμL of molecular biology grade water  and incubated at 56 °C for 10 min to elute . RNA concentration  was 665 \nquantified via Nanodrop (Thermo Scientific) and Qubit fluorometer (Invitrogen).  666 \n 667 \n15 μg of the eluted RNA was treated with TURBO DNase (Invitrogen AM2238). Then, TRIzol precipitation and 668 \nwashes from the previous paragraph was performed once again on the treated RNA. This time the sample was 669 \neluted in 15-20 μL of molecular biology grade water and quantified via TapeStation (Agilent) to ensure a RIN > 670 \n5. Ribosomal RNA depletion was performed on  1000 ng of  the sample using the Illumina Ribo -Zero rRNA 671 \nRemoval Kit (Illumina 20040526) following the kit protocol and using RNAClean XP beads (Beckman Coulter 672 \nA63987) for the cleanup steps. Synthesis of cDNA and Illumina TruSeq adapter ligation was performed on 100 673 \nng of depleted RNA using the NEBNext Ultra II Directional RNA Library Prep Kit (NEB E7760) following the kit 674 \nprotocol and using AMPure XP beads (Beckman Coulter A63880) for the cleanup steps . The prepared RNA 675 \nlibrary was sequenced on an Illumina NextSeq 2000 system. Sequencing reads were aligned to the retron 676 \noperon using Geneious.  677 \n 678 \nhyRNA secondary structure prediction 679 \nThe hyRNA sequence from RNA-seq was input into RNAfold58, which predicts secondary structure by 680 \ncomputing minimum free energy, and visualized with forna59. The resulting structure (see Fig 3b) was 681 \nindependently supported by covariance modeling (see below and Fig S3a). 682 \n 683 \nhyRNA covariance modeling 684 \nHomologs of Retron-Vpa2 were identified using the RT amino acid sequence (WP_069548583.1) as a BLASTP 685 \nsearch query of the NR protein database (max target sequences = 100). Flanking nucleotide sequences 1 kb 686 \nupstream and downstream of RT genes were extracted, clustered at 99.9% sequence identity to remove 687 \nreplicates with CD-HIT60 (v4.8.1), and aligned with MAFFT53 (v7.525). The resulting alignment was trimmed 688 \nat the 5′ and 3′ ends to the boundaries of the Retron-Vpa2 hyRNA as determined by RNA-seq. These putative 689 \nhyRNA homologs were clustered at 95% sequence identity with CD-HIT and realigned with mLocARNA61 690 \n(v2.0.1). The resulting structure-based multiple sequence alignment was used to build and calibrate a 691 \ncovariance model (CM) of the hyRNA using the Infernal suite62 (v1.1.5). 692 \n 693 \nTo identify other hyRNA homologs across additional Type VI retron loci using this CM, a database of diverse 694 \nType VI retron RTs was generated as previously described (see Phylogenetic analysis). Briefly, a broad RT 695 \nPfam profile (PF00078) was searched against the ClusteredNR protein database. Hits were retrieved and 696 \naligned against the HMM profile using hmmalign –trim from the HMMER suite63, and all sequences in the 697 \nalignment shorter than 150 amino acids were removed with SeqKit64.  The alignment of remaining RTs (n = 698 \n306,351 members) was used to build a tree with VeryFastTree54. Tree members were annotated with the best 699 \nhit from an MMseqs2 easy-search65 against a custom reference database of known RTs, revealing a clade of 700 \nType VI retron RTs (n = 779 members) that contained Retron-Vpa2.  701 \n 702 \nThe 1-kb nucleotide flanking regions of these Type VI retron RTs were extracted and clustered at 95% 703 \nsequence identity using CD-HIT. The CMsearch function of Infernal was then used to scan the previously built 704 \nCM through these sequences to identify additional hyRNA homologs, and the final hits (n = 81 Type VI retron 705 \nloci, including Retron-Vpa2) were evaluated for statistically significant co-varying base pairs with R-scape66 706 \n(v1.4.0) at a default E-value threshold of 0.05 (Fig S3a). To visualize conserved RNA motifs, a sequence logo 707 \nwas generated from the structure-based multiple sequence alignment of these hyRNA hits using WebLogo67 708 \n(v3.7.9) (Fig S3b). 709 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n 710 \nToxicity assays 711 \nA plasmid containing wild-type or mutant Retron-Vpa2 (or an empty vector control) was  transformed into E. 712 \ncoli bSLS.114. Individual colonies were picked and grown until saturation in 3 mL LB + carbenicillin at 37 °C. 713 \nCultures were then diluted 1:100 in 200 uL of fresh LB + carbenicillin in a transparent 96-well microtiter plate 714 \n(Corning 3596). The plate was covered with an air -permeable seal (Breathe-Easy, Diversified Biotech BEM-1) 715 \nand incubated in a microplate reader (Tecan Infinite 200 PRO) with shaking at 37° C for 16 hrs with OD 600 716 \nmeasured every 10 min. At the 2 hr timepoint, IPTG and arabinose inducers were added to all cultures. Raw 717 \nOD600 values were adjusted by first multiplying by the empirically determined correction factor of 718 \n6.230707455, then subtracting the 0 hr measurement for each sample.  719 \n 720 \nPull-down assays 721 \nThree plasmids were constructed containing the Retron-Vpa2 operon with either a C-terminal Flag tag on the 722 \nRT, a C-terminal HA tag on the HTH, or both tags. For assays using RT as bait, a non-tagged RT plasmid was 723 \nused as a negative control. Similarly, for assays using HTH as bait, a non-tagged HTH plasmid served as the 724 \nnegative control. A modified E. coli MG1655 strain lacking bacterial defense systems and stably integrating 725 \nthe λDE3 prophage (carrying an IPTG-inducible T7 RNA polymerase gene) was generated using the λDE3 726 \nLysogenization Kit (Novagen). This strain is referred to as MG1655-DE3. 727 \n 728 \nPlasmids were transformed into E. coli MG1655-DE3 and plated on LB agar supplemented with 100 µg/mL 729 \nampicillin. Liquid cultures were grown at 37°C with shaking until reaching an OD₆₀₀ of 0.6. Expression of the 730 \noperon was induced with 0.2 mM IPTG for 4 hours at 37°C. After induction, OD₆₀₀ was measured and used to 731 \nnormalize sample input across conditions. According to the OD600, cells were harvested at >3000 × g for 10 732 \nminutes at 4°C and resuspended in 0.8–2 mL of lysis buffer (20 mM Tris-HCl pH 7.5, 150 mM KCl, 1 mM 733 \nMgCl₂, 1 mM TCEP, 0.1% Triton X-100) to ensure that the total cell concentration is similar in all tubes. Lysis 734 \nwas performed by sonication (3 pulses of 10 seconds at 25% amplitude). Soluble proteins were separated by 735 \ncentrifugation at >16,000 × g for 30 minutes at 4°C. The clarified lysate (800 µL) was incubated with 50 µL of 736 \nanti-FLAG M2 affinity gel (Sigma-Aldrich, A2220) or anti-HA agarose beads (Thermo, 26181) prewashed with 737 \nlysis buffer, for 30 minutes at 4°C with gentle agitation. Beads were pelleted at 3000 × g for 3 minutes, and 738 \nthe supernatant was removed. Bound complexes were washed four times with lysis buffer lacking Triton X-739 \n100. Beads were resuspended in ~60 µL of lysis buffer without triton X-100 and divided into three aliquots 740 \n(~20 µL each) for protein, RNA, and DNA analysis. 741 \n 742 \nFor Protein Analysis by SDS-PAGE and Western Blot, the samples were mixed with NuPAGE™ LDS Sample 743 \nBuffer (Invitrogen, NP0007) and denatured by boiling for 5–10 minutes at 95°C. Proteins were separated by 744 \nSDS-PAGE (NuPAGE™ Bis-Tris Mini Protein Gels 4 12–% Acrylamide, Invitrogen, 12030166) and transferred 745 \nto a nitrocellulose membrane (iBlot™ 3 Transfer Stacks, Invitrogen, IB33001). Entire pull-down samples (40 746 \nµL) were loaded, as well as 30 µg of total protein from the input lysate to assess expression levels. The 747 \nmembranes were then blocked for 1 hour at room temperature in TBST (Tris-buffered saline with 0.1% 748 \nTween-20) containing 0.5% powdered milk. Detection was performed using HRP-conjugated anti-HA (Anti-749 \nHA-Peroxidase, High Affinity, Roche) or anti-FLAG M\" antibodies (Sigma, A8592), diluted 1:1000 in TBST with 750 \n0.5% powdered milk. Membranes were washed three times for 10 minutes in TBST without milk, and 751 \nchemiluminescence was visualized using SuperSignal™ West Pico PLUS Chemiluminescent Substrate (Thermo 752 \nScientific, 34577). 753 \n 754 \nFor RNA analysis, samples were resuspended in urea loading dye (4M Urea, 1 mM Tris pH 7.5, 5 mM EDTA), 755 \ndenatured by boiling for 10 minutes at 95°C, and cooled on ice for 2 minutes. DNA samples were pretreated 756 \nwith ~ 100 μg/ml of RNAse A for 10 minutes at 37°C before denaturation. RNA and DNA were resolved on 757 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n10% TBE-Urea gels (Invitrogen, EC68755BOX) run at ~180 V for 70–90 minutes. RNA and DNA was 758 \nvisualized by staining with GelRed Nucleic Acid Stain (10000X DMSO, Millipore, SCT122) for 10 minutes prior 759 \nto imaging. 760 \n 761 \nExtraction of RNA/DNA from pull-down samples for sequencing 762 \nPull-down samples were prepared as described above. The entire sample (20 µL of beads) was resuspended 763 \nin 300 µL of nuclease-free ddH₂O and denatured by boiling at 95 °C for 5–10 minutes to release proteins and 764 \nnucleic acids from the beads. RNA was then purified using the phenol:chloroform extraction method, followed 765 \nby precipitation with sodium acetate and ethanol. Purified RNA was directly prepped for sequencing using 766 \nthe NEBNext Ultra II Directional RNA Library Prep Kit (NEB E7760) and sequenced on an Illumina NextSeq 767 \n2000 system. Sequencing reads were aligned to the reference operon using Geneious.  768 \n 769 \nSpacer-seq 770 \nSpacer-seq experiments were conducted similarly to the naïve CRISPR -Cas adaptation  experiments in a 771 \nprevious manuscript50. A plasmid containing wild -type or mutant Retron-Vpa2 was co -transformed with the 772 \nCas1-Cas2-expressing plasmid into E. coli bSLS.114. For experiments involving trigger genes, a third plasmid 773 \ncontaining the trigger gene was also transformed. Three colonies (biological replicates) for each sample were 774 \npicked and grown until saturation in 0.5 mL LB + carbenicillin + spectinomycin (+ chloramphenicol when the 775 \ntrigger-containing plasmid is present) in a 96-well deep well plate (Thermo Scientific 260251) covered with 776 \nan air -permeable seal (Breathe-Easier, Diversified Biotech BERM-2000) at 37°C with shaking at 1000 rpm. 777 \nCultures were then diluted 1:100 in 0.5 mL LB + antibiotics  + IPTG + arabinose in a new deep well plate and 778 \ngrown for 24 hrs at 37°C. Cells were then harvested by boiling 25 μL of culture with 25 μL of molecular biology 779 \ngrade water at 95°C for 10 min. 0.5 μL of each boiled lysate was then used as template for PCR amplification of 780 \nthe CRISPR array . Specifically, the primers used in this PCR (Supplementary Table 5) amplified the region 781 \nbetween the leader sequence of the array and the second pre -existing spacer. The primers also introduce d 782 \nIllumina TruSeq adapters to the amplicons . Three versions of the  forward primer were used , each with a 783 \ndifferent barcode assigned to a different biological replicate, enabling pooling of replicates for each sample 784 \nprior to sequencing. Amplicons were then prepared for Illumina sequencing using standard workflows and 785 \nsequenced on a NextSeq 2000 system.  786 \n 787 \nTo analyze sequencing data, Illumina FASTQ files containing sequencing reads for each sample were 788 \ndemultiplexed using the custom barcodes introduced to each biological replicate. Demultiplexed files were then 789 \nprocessed using custom code ( https://github.com/Shipman-Lab/Spacer-Seq) developed in  a  previous 790 \nmanuscript33, which handles read trimming and extraction of spacer  sequences. From this, a FASTA file of all 791 \nnewly acquired spacers was obtained for each biological replicate of each sample. These were again pooled into 792 \na single FASTA file for each sample and aligned using Bowtie68 (v1.3.1) to a dictionary of references containing 793 \nthe Retron-expressing plasmid, the Cas1 -Cas2-expressing plasmid, and the E. coli  BL21-AI genome . Reads 794 \naligning to multiple references were discarded. Coverage per nucleotide position was determined from aligned 795 \nreads using SAMtools57 (v1.18) and normalized either to the total number of newly acquired spacers for a given 796 \nsample, or to the maximum coverage observed in the Retron-expressing plasmid. Data analysis was performed 797 \non Jupyter Lab69 (v4.0.7). Visualization of spacer coverage on the Retron-Vpa2 hyRNA secondary structure was 798 \nperformed using VARNA 70 (v3-93). Visualization of spacer coverage across the full plasmid and genome  as 799 \nshown in Fig S6c was performed using Circos71 (v0.69-9).  800 \n 801 \nDirect msDNA sequencing  802 \nFollowing a protocol from a previous publication15, 79 μL of eluted DNA from a bacterial culture midiprep (see 803 \nmsDNA expression and gel analysis) was treated with 10 μL of DBR1 (50 ng/μL in-house prep, see previous 804 \nmanuscript23) and RNaseH (NEB M0297)  in 10 μL of rCutSmart buffer (NEB  B6004) and molecular biology 805 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\ngrade water to reach a total reaction volume of 100 μL. The reaction was incubated at 37 °C for 30 min. Then, 806 \nthe product was cleaned up with Zymo’s ssDNA/RNA Clean & Concentrator kit (Zymo Research D7010 ) and 807 \neluted in 15 μL of molecular biology grade water.  808 \n 809 \nPoly(A)-tailing of cleaned up msDNA was performed by mixing 12.25 μL of molecular biology grade water, 2.5 810 \nμL of terminal transferase reaction  buffer (NEB B0315), 0.25 μL dATP (of a 10 mM stock diluted from NEB 811 \nN0440), 3 μL of terminal transferase (TdT, NEB M0315), and 10 μL of msDNA. Reaction was incubated at room 812 \ntemperature for exactly 60 sec and  then stopped by heating to 70°C for 5 min (~25 adenosines should have 813 \nbeen added in this time).  814 \n 815 \nAfter this, complem entary strand synthesis was performed using a poly (T) primer which also contains an  816 \nIllumina TruSeq adapter ( Supplementary Table 5) . The reaction was setup as follows: 28.9 μL of molecular 817 \nbiology grade water, 8 μL of TdT reaction from the previous step, 5 μL of NEBuffer 2 (NEB B7002), 0.1 μL of the 818 \npoly(T) primer (of a 100 uM stock from IDT), and 5 μL of dNTP mix (of a 10 mM stock from NEB N0447) were 819 \nmixed and heated to 80°C, then allowed to cool to room temperature (facilitates primer annealing). Then, 3 μL 820 \n(15 units)  of DNA Polymerase I, Large (Klenow) Fragment (NEB M0210) was added, and the reaction was 821 \nincubated at 37°C for 30 min. The reaction was purified using the QIAquick PCR Purification Kit (Qiagen 28104) 822 \nand eluted in 15 uL of molecular biology grade water.  823 \n 824 \nFollowing this, oligos containing the second Illumina TruSeq adapter (Supplementary Table 5) were annealed 825 \nby mixing 10 μL of each oligo (top and bottom strand, from a 100 μM stock concentration) with 10 μL of 826 \nmolecular biology grade water and 10 μL of NEBuffer 2, then heated to 95°C for 2 min and allowed to cool back 827 \nto room temperature. 1 μL of the annealed adapters were mixed with 4 μL of the extended msDNA product and 828 \n5 μL of Blunt/TA Ligase Master Mix ( NEB M0367). This was incubated for 15 min at room temperature, after 829 \nwhich the reaction was immediately cleaned up using AMPure XP Beads (Beckman Coulter A63880).  Cleanup 830 \nwas performed with 1.8x beads:DNA volume ratio and the product was eluted in 10 uL of molecular biology 831 \ngrade water. Sample was prepared for Illumina sequencing following standard workflows and sequenced on a 832 \nNextSeq 2000 system. Sequencing reads were trimmed for adapter sequences and aligned to the Retron-Vpa2 833 \noperon using Geneious. Visualization of read coverage on the Retron-Vpa2 hyRNA secondary structure was 834 \nperformed using VARNA70 (v3-93). 835 \n 836 \nSplit GFP complementation assays 837 \nIn preliminary experiments developing this assay, a plasmid containing either sfGFP1-10 or SP-sfGFP1-10 with 838 \na kanamycin-resistance cassette and a plasmid containing sfGFP11 (either on its own or fused to the SP in the 839 \ncontext of Retron-Vpa2 Δnc-hyRNA) with a carbenicillin-resistance cassette were co- transformed into E. coli 840 \nbSLS.114. In all subsequent experiments, a plasmid containing a version of Retron-Vpa2 SP::sfGFP11 with a 841 \ncarbenicillin-resistance cassette and a plasmid containing SP-sfGFP1-10 with a kanamycin-resistance cassette 842 \nwere co-transformed into E. coli bSLS.114. For experiments requiring trigger genes, a third plasmid containing 843 \nthe trigger (or an empty vector control) with a chloramphenicol-resistance cassette was also co-transformed. 844 \nFor experiments using  Retron-Eco1 to produce msDNA  in trans, a third plasmid containing the engineered 845 \nretron with a spectinomycin-resistance cassette was also co-transformed.  846 \n 847 \nFor each sample, three individual colonies were picked (biological replicates) and grown until saturation in 0.5 848 \nmL LB + antibiotics in a 96 -well deep well plate (Thermo Scientific 260251 ) covered with an air -permeable 849 \nseal (Breathe-Easier, Diversified Biotech) at 37°C with shaking at 1000 rpm. Cultures were then diluted 1:100 850 \nin 200 uL of fresh LB + antibiotics in a black, clear- bottom 96-well microtiter plate (Corning 3603). Samples 851 \nwere arranged to have at least one empty well between each set of replicates to avoid bleed -through 852 \nfluorescence across samples. The plate was covered with an air -permeable seal (Bre athe-Easier, Diversified 853 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nBiotech BERM-2000) and incubated in a microplate reader (Tecan Infinite 200 PRO) with shaking at 37° C for 854 \n16 hrs with fluorescence measured from the bottom of the wells using an excitation wavelength of 485 nm and 855 \nemission wavelength of 520 nm . Measurements were taken every 10 min. At the 2 hr timepoint, IPTG and 856 \narabinose inducers were added to all cultures. In preliminary experiments developing the assay, inducers were 857 \ninstead added at the 0 hr time point and the experiment was run for 16 hrs. Fluorescence readings were 858 \nadjusted by subtracting the 0 hr measurement for each sample. 859 \n 860 \nACKNOWLEDGEMENTS 861 \nWork was supported by funding from the National Science Foundation (MCB 2509382), the Robert 862 \nJ Kleeberg Jr. and Helen C. Kleberg Foundation, and the Gordon and Betty Moore Foundation.  863 \nK.Z. was supported by a National Science Foundation Graduate Research Fellowship and a UCSF 864 \nDiscovery Fellowship.  865 \nA.C. was supported by a Lundbeck Foundation grant R380-2021-1448.  866 \nD.J.Z. was supported by the Fulbright U.S. Student Program, which is sponsored by the U.S. 867 \nDepartment of State and the Danish-American Fulbright Commission.  868 \nA.D.H was supported by a Ph.D. fellowship from Junta de Castilla y León and European Social Fund 869 \nPlus (EDU/1868/2022) and an EMBO Scientific Exchange Grant (11561).  870 \nG.M. is part of CPR, which is supported financially by the Novo Nordisk Foundation (NNF14CC0001, 871 \nNNF24SA0098829). This work was also supported by the ERC-AdG 101096548 (INTETOOLS), 872 \nNNF0024386, NNF17SA0030214, and NNF18OC0055061 grants to G.M, who is a member of the 873 \nIntegrative Structural Biology Cluster (ISBUC) at the University of Copenhagen.  874 \nR.P.-R. was supported by a Lundbeck Foundation grant (R347-2020-2346) and a research grant 875 \n(VIL60763) from VILLUM FONDEN.   876 \nS.L.S. is a Chan Zuckerberg Biohub – San Francisco Investigator.  877 \n 878 \nWe thank Mylinh Bernardi and Felicia Miller of the Gladstone Genomics Core for their assistance 879 \nwith sequencing co-immunoprecipitated RNA from our pull-down assays.  880 \n 881 \nAUTHOR CONTRIBUTIONS 882 \nAuthor contributions follow the CRedIT taxonomy 883 \n(https://www.elsevier.com/researcher/author/policies-and-guidelines/credit-author-statement) 884 \n 885 \nK.Z.: Methodology, Investigation, Validation, Formal Analysis, Writing – Original Draft, Visualization 886 \nM.R-M.: Methodology, Investigation, Validation, Writing – Review & Editing 887 \nD.P.: Methodology, Investigation, Validation, Writing – Review & Editing 888 \nJ.M.K.: Methodology, Investigation, Validation, Writing – Review & Editing, Visualization 889 \nA.C.: Methodology, Investigation, Validation, Writing – Review & Editing, Visualization 890 \nD.J.Z.: Methodology, Software, Formal Analysis, Writing – Review & Editing, Visualization 891 \nM.R.M.: Methodology, Software, Formal Analysis, Writing – Review & Editing, Visualization 892 \nR.Z.: Methodology, Investigation, Validation 893 \nA.D-H.: Investigation 894 \nG.M.: Methodology, Writing – Review & Editing, Supervision, Funding Acquisition 895 \nR.P-R.: Methodology, Writing – Review & Editing, Supervision, Funding Acquisition 896 \nA.G-D.: Conceptualization, Writing – Review & Editing, Supervision 897 \nS.L.S.: C onceptualization, Writing  – Review & Editing, Visualization, Supervision, Project 898 \nadministration, Funding acquisition 899 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n 900 \nCOMPETING INTERESTS 901 \nG.M. is a stockholder of Ensoma and has been consultant for Orbis Medicines. The remaining 902 \nauthors declare no competing interests.  903 \n 904 \nDATA AND CODE AVAILABILITY 905 \nSequencing data associated with this study are available in the NCBI SRA (PRJNA1320719) 906 \nhttps://www.ncbi.nlm.nih.gov/bioproject/PRJNA1320719  907 \n 908  \nCustom code used to process or analyze data from this study is available at: 909 \nSpacer-seq analysis was performed using code sourced from https://github.com/Shipman-910 \nLab/Spacer-Seq. 911 \n 912  \nSUPPLEMENTARY MATERIALS 913 \nSee supplementary files for Supplementary Figures S1-S10 and Supplementary Tables 1-5.  914 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\nREFERENCES 915 \n1. Gao, L. et al. Diverse enzymatic activities mediate antiviral immunity in prokaryotes. Science 916 \n369, 1077–1084 (2020). 917 \n2. Millman, A. et al. Bacterial Retrons Function In Anti-Phage Defense. Cell 183, 1551-1561.e12 918 \n(2020). 919 \n3. Bobonis, J. et al. Bacterial retrons encode phage-defending tripartite toxin–antitoxin systems. 920 \nNature 609, 144–150 (2022). 921 \n4. Bernheim, A. & Sorek, R. The pan-immune system of bacteria: antiviral defence as a community 922 \nresource. Nat. Rev. Microbiol. 18, 113–119 (2020). 923 \n5. Lopatina, A., Tal, N. & Sorek, R. Abortive Infection: Bacterial Suicide as an Antiviral Immune 924 \nStrategy. Annu. Rev. Virol. 7, 371–384 (2020). 925 \n6. Mestre, M. R., González-Delgado, A., Gutiérrez-Rus, L. I., Martínez-Abarca, F. & Toro, N. 926 \nSystematic prediction of genes functionally associated with bacterial retrons and classification 927 \nof the encoded tripartite systems. Nucleic Acids Res. 48, 12632–12647 (2020). 928 \n7. Hsu, M. Y., Inouye, S. & Inouye, M. Structural Requirements of the RNA Precursor for the 929 \nBiosynthesis of the Branched RNA-linked Multicopy Single-stranded DNA of Myxococcus 930 \nxanthus. J. Biol. Chem. 264, 6214–6219 (1989). 931 \n8. Inouye, S., Hsu, M.-Y., Eagle, S. & Inouye, M. Reverse transcriptase associated with the 932 \nbiosynthesis of the branched RNA-linked msDNA in Myxococcus xanthus. Cell 56, 709–717 933 \n(1989). 934 \n9. Lampson, B. C., Inouye, M. & Inouye, S. Reverse transcriptase with concomitant ribonuclease H 935 \nactivity in the cell-free synthesis of branched RNA-linked msDNA of Myxococcus xanthus. Cell 936 \n56, 701–707 (1989). 937 \n10. Hsu, M. Y., Eagle, S. G., Inouye, M. & Inouye, S. Cell-free synthesis of the branched RNA-linked 938 \nmsDNA from retron-Ec67 of Escherichia coli. J. Biol. Chem. 267, 13823–13829 (1992). 939 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n11. Shimamoto, T., Inouye, M. & Inouye, S. The Formation of the 2′,5′-Phosphodiester Linkage in the 940 \ncDNA Priming Reaction by Bacterial Reverse Transcriptase in a Cell-free System (∗). J. Biol. 941 \nChem. 270, 581–588 (1995). 942 \n12. Shimamoto, T., Kawanishi, H., Tsuchiya, T., Inouye, S. & Inouye, M. In Vitro Synthesis of 943 \nMulticopy Single-Stranded DNA, Using Separate Primer and Template RNAs, by Escherichia coli 944 \nReverse Transcriptase. J. Bacteriol. 180, 2999–3002 (1998). 945 \n13. Simon, A. J., Ellington, A. D. & Finkelstein, I. J. Retrons and their applications in genome 946 \nengineering. Nucleic Acids Res. 47, 11007–11019 (2019). 947 \n14. Yee, T., Furuichi, T., Inouye, S. & Inouye, M. Multicopy single-stranded DNA isolated from a 948 \ngram-negative bacterium, Myxococcus xanthus. Cell 38, 203–209 (1984). 949 \n15. Khan, A. G. et al. An experimental census of retrons for DNA production and genome editing. 950 \nNat. Biotechnol. 43, 914–922 (2025). 951 \n16. Nakamura, K. L., Zhang, K., Mestre, M. R., Rojas-Montero, M. & Shipman, S. L. Phage defense and 952 \ngenome editing using novel retrons sourced from isolated environmental bacteria. 953 \n2025.01.29.635429 Preprint at https://doi.org/10.1101/2025.01.29.635429 (2025). 954 \n17. Lampson, B. C. et al. Reverse Transcriptase in a Clinical Strain of Escherichia coli: Production of 955 \nBranched RNA-Linked msDNA. Science 243, 1033–1038 (1989). 956 \n18. Lim, D. & Maas, W. K. Reverse transcriptase-dependent synthesis of a covalently linked, 957 \nbranched DNA-RNA compound in E. coli B. Cell 56, 891–904 (1989). 958 \n19. Lima, T. M. O. & Lim, D. A Novel Retron That Produces RNA-less msDNA inEscherichia coliUsing 959 \nReverse Transcriptase. Plasmid 38, 25–33 (1997). 960 \n20. Carabias, A. et al. Retron-Eco1 assembles NAD+-hydrolyzing filaments that provide immunity 961 \nagainst bacteriophages. Mol. Cell 84, 2185-2202.e12 (2024). 962 \n21. Wang, Y. et al. Cryo-EM structures of Escherichia coli Ec86 retron complexes reveal architecture 963 \nand defence mechanism. Nat. Microbiol. 7, 1480–1489 (2022). 964 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n22. Wang, Y. et al. DNA methylation activates retron Ec86 filaments for antiphage defense. Cell Rep. 965 \n43, (2024). 966 \n23. Palka, C., Fishman, C. B., Bhattarai-Kline, S., Myers, S. A. & Shipman, S. L. Retron reverse 967 \ntranscriptase termination and phage defense are dependent on host RNase H1. Nucleic Acids 968 \nRes. 50, 3490–3504 (2022). 969 \n24. Wang, B., Hoffman, R., Hou, Y.-M. & Li, H. Structural Basis for Retron Co-option of Anti-phage 970 \nATPase-nuclease. 2025.05.10.653283 Preprint at https://doi.org/10.1101/2025.05.10.653283 971 \n(2025). 972 \n25. Wang, C. et al. Disassembly activates Retron-Septu for antiphage defense. Science 0, eadv3344 973 \n(2025). 974 \n26. Li, X. et al. Architecture and mechanism of a dual-enzyme retron system in prokaryotic 975 \nimmunity. 2025.05.23.655875 Preprint at https://doi.org/10.1101/2025.05.23.655875 976 \n(2025). 977 \n27. García-Rodríguez, F. M., Martínez-Abarca, F., Wilkinson, M. E. & Toro, N. Similar mechanisms of 978 \nretron-mediated anti-phage defense for different families of tailed phages. 2024.02.09.579579 979 \nPreprint at https://doi.org/10.1101/2024.02.09.579579 (2024). 980 \n28. Stokar- Avihail, A. et al. Discovery of phage determinants that confer sensitivity to bacterial 981 \nimmune systems. Cell 186, 1863-1876.e16 (2023). 982 \n29. Pausch, P. et al. Retron Eco2 breaks tRNAs for antiphage defense. Preprint at 983 \nhttps://doi.org/10.21203/rs.3.rs-7196019/v1 (2025). 984 \n30. Wang, Y. et al. Structural basis of the RNA-mediated Retron-Eco2 oligomerization. Cell Discov. 985 \n11, 73 (2025). 986 \n31. Rychlik, I., Sebkova, A., Gregorova, D. & Karpiskova, R. Low-Molecular-Weight Plasmid of 987 \nSalmonella enterica Serovar Enteritidis Codes for Retron Reverse Transcriptase and Influences 988 \nPhage Resistance. J. Bacteriol. 183, 2852–2858 (2001). 989 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n32. Pilousova, L. & Rychlik, I. Retron Se72 utilizes a unique strategy of the self-priming initiation of 990 \nreverse transcription. Cell. Mol. Life Sci. 68, 3607–3617 (2011). 991 \n33. Bhattarai-Kline, S. et al. Recording gene expression order in DNA by CRISPR addition of retron 992 \nbarcodes. Nature 608, 217–225 (2022). 993 \n34. Maffei, E. et al. Systematic exploration of Escherichia coli phage–host interactions with the 994 \nBASEL phage collection. PLOS Biol. 19, e3001424 (2021). 995 \n35. Caldwell, B. J. & Bell, C. E. Structure and mechanism of the Red recombination system of 996 \nbacteriophage λ. Prog. Biophys. Mol. Biol. 147, 33–46 (2019). 997 \n36. Murphy, K. C. Lambda Gam protein inhibits the helicase and chi-stimulated recombination 998 \nactivities of Escherichia coli RecBCD enzyme. J. Bacteriol. 173, 5808–5821 (1991). 999 \n37. Wilkinson, M. et al. Structural basis for the inhibition of RecBCD by Gam and its synergistic 1000 \nantibacterial effect with quinolones. eLife 5, e22963 (2016). 1001 \n38. Smith, C. E. & Bell, C. E. Domain Structure of the Redβ Single-Strand Annealing Protein: the C-1002 \nterminal Domain is Required for Fine-Tuning DNA-binding Properties, Interaction with the 1003 \nExonuclease Partner, and Recombination in vivo. J. Mol. Biol. 428, 561–578 (2016). 1004 \n39. Iyer, L. M., Koonin, E. V. & Aravind, L. Classification and evolutionary history of the single-strand 1005 \nannealing proteins, RecT, Redβ, ERF and RAD52. BMC Genomics 3, 8 (2002). 1006 \n40. Lampson, B. C., Inouye, M. & Inouye, S. Retrons, msDNA, and the bacterial genome. Cytogenet. 1007 \nGenome Res. 110, 491–499 (2005). 1008 \n41. Datsenko, K. A. et al. Molecular memory of prior infections activates the CRISPR/Cas adaptive 1009 \nbacterial immunity system. Nat. Commun. 3, 945 (2012). 1010 \n42. Yosef, I., Goren, M. G. & Qimron, U. Proteins and DNA elements essential for the CRISPR 1011 \nadaptation process in Escherichia coli. Nucleic Acids Res. 40, 5569–5576 (2012). 1012 \n43. Nuñez, J. K. et al. Cas1–Cas2 complex formation mediates spacer acquisition during CRISPR– Cas 1013 \nadaptive immunity. Nat. Struct. Mol. Biol. 21, 528–534 (2014). 1014 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n44. Nuñez, J. K., Harrington, L. B., Kranzusch, P. J., Engelman, A. N. & Doudna, J. A. Foreign DNA 1015 \ncapture during CRISPR–Cas adaptive immunity. Nature 534, S13–S14 (2016). 1016 \n45. Cabantous, S., Terwilliger, T. C. & Waldo, G. S. Protein tagging and detection with engineered 1017 \nself-assembling fragments of green fluorescent protein. Nat. Biotechnol. 23, 102–107 (2005). 1018 \n46. Levy, A. et al. CRISPR adaptation biases explain preference for acquisition of foreign DNA. 1019 \nNature 520, 505–510 (2015). 1020 \n47. Shiriaeva, A. A. et al. Host nucleases generate prespacers for primed adaptation in the E. coli 1021 \ntype I-E CRISPR-Cas system. Sci. Adv. 8, eabn8650 (2022). 1022 \n48. Lopez, S. C., Crawford, K. D., Lear, S. K., Bhattarai-Kline, S. & Shipman, S. L. Precise genome 1023 \nediting across kingdoms of life using retron-derived DNA. Nat. Chem. Biol. 18, 199–206 (2022). 1024 \n49. Datsenko, K. A. & Wanner, B. L. One-step inactivation of chromosomal genes in Escherichia coli 1025 \nK-12 using PCR products. Proc. Natl. Acad. Sci. 97, 6640–6645 (2000). 1026 \n50. Lopez, S. C., Lee, Y., Zhang, K. & Shipman, S. L. SspA is a transcriptional regulator of CRISPR 1027 \nadaptation in E. coli. Nucleic Acids Res. 53, gkae1244 (2025). 1028 \n51. Püllmann, P. et al. A modular two yeast species secretion system for the production and 1029 \npreparative application of unspecific peroxygenases. Commun. Biol. 4, 562 (2021). 1030 \n52. Fishman, C. B. et al. Continuous multiplexed phage genome editing using recombitrons. Nat. 1031 \nBiotechnol. 43, 1299–1310 (2025). 1032 \n53. Katoh, K. & Standley, D. M. MAFFT Multiple Sequence Alignment Software Version 7: 1033 \nImprovements in Performance and Usability. Mol. Biol. Evol. 30, 772–780 (2013). 1034 \n54. Piñeiro, C., Abuín, J. M. & Pichel, J. C. Very Fast Tree: speeding up the estimation of phylogenies 1035 \nfor large alignments through parallelization and vectorization strategies. Bioinformatics 36, 1036 \n4658–4659 (2020). 1037 \n55. Payne, L. J. et al. Identification and classification of antiviral defence systems in bacteria and 1038 \narchaea with PADLOC reveals new system types. Nucleic Acids Res. 49, 10868–10878 (2021). 1039 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n56. Mazzocco, A., Waddell, T. E., Lingohr, E. & Johnson, R. P. Enumeration of Bacteriophages Using 1040 \nthe Small Drop Plaque Assay System. in Bacteriophages 81–85 (Humana Press, 2009). 1041 \ndoi:10.1007/978-1-60327-164-6_9. 1042 \n57. Danecek, P. et al. Twelve years of SAMtools and BCFtools. GigaScience 10, giab008 (2021). 1043 \n58. Lorenz, R. et al. ViennaRNA Package 2.0. Algorithms Mol. Biol. 6, 26 (2011). 1044 \n59. Kerpedjiev, P., Hammer, S. & Hofacker, I. L. Forna (force-directed RNA): Simple and effective 1045 \nonline RNA secondary structure diagrams. Bioinformatics 31, 3377–3379 (2015). 1046 \n60. Li, W. & Godzik, A. Cd-hit: a fast program for clustering and comparing large sets of protein or 1047 \nnucleotide sequences. Bioinformatics 22, 1658–1659 (2006). 1048 \n61. Will, S., Joshi, T., Hofacker, I. L., Stadler, P. F. & Backofen, R. LocARNA-P: Accurate boundary 1049 \nprediction and improved detection of structural RNAs. RNA 18, 900–914 (2012). 1050 \n62. Nawrocki, E. P. & Eddy, S. R. Infernal 1.1: 100-fold faster RNA homology searches. 1051 \nBioinformatics 29, 2933–2935 (2013). 1052 \n63. Eddy, S. R. Accelerated Profile HMM Searches. PLOS Comput. Biol. 7, e1002195 (2011). 1053 \n64. Shen, W., Le, S., Li, Y. & Hu, F. SeqKit: A Cross-Platform and Ultrafast Toolkit for FASTA/Q File 1054 \nManipulation. PLOS ONE 11, e0163962 (2016). 1055 \n65. Steinegger, M. & Söding, J. MMseqs2 enables sensitive protein sequence searching for the 1056 \nanalysis of massive data sets. Nat. Biotechnol. 35, 1026–1028 (2017). 1057 \n66. Rivas, E., Clements, J. & Eddy, S. R. A statistical test for conserved RNA structure shows lack of 1058 \nevidence for structure in lncRNAs. Nat. Methods 14, 45–48 (2017). 1059 \n67. Crooks, G. E., Hon, G., Chandonia, J.-M. & Brenner, S. E. WebLogo: A Sequence Logo Generator. 1060 \nGenome Res. 14, 1188–1190 (2004). 1061 \n68. Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efficient alignment of 1062 \nshort DNA sequences to the human genome. Genome Biol. 10, R25 (2009). 1063 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint \n\n69. Granger, B. E. & Pérez, F. Jupyter: Thinking and Storytelling With Code and Data. Comput. Sci. 1064 \nEng. 23, 7–14 (2021). 1065 \n70. Darty, K., Denise, A. & Ponty, Y. VARNA: Interactive drawing and editing of the RNA secondary 1066 \nstructure. Bioinformatics 25, 1974–1975 (2009). 1067 \n71. Krzywinski, M. I. et al. Circos: An information aesthetic for comparative genomics. Genome Res. 1068 \n(2009) doi:10.1101/gr.092759.109. 1069 \n 1070 \n.CC-BY 4.0 International licenseavailable under a \n(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made \nThe copyright holder for this preprintthis version posted October 22, 2025. ; https://doi.org/10.1101/2025.10.22.683967doi: bioRxiv preprint","source_license":"CC-BY-4.0","license_restricted":false}