Rotation-induced macromolecular spooling of DNA

Genetic information is stored in a linear sequence of base-pairs; however, thermal fluctuations and complex DNA conformations such as folds and loops make it challenging to order genomic material for in vitro analysis. In this work, we discover that rotation-induced macromolecular spooling of DNA around a rotating microwire can monotonically order genomic bases, overcoming this challenge. We use single-molecule fluorescence microscopy to directly visualize long DNA strands deforming and elongating in shear flow near a rotating microwire, in agreement with numerical simulations. While untethered DNA is observed to elongate substantially, in agreement with our theory and numerical simulations, strong extension of DNA becomes possible by introducing tethering. For the case of tethered polymers, we show that increasing the rotation rate can deterministically spool a substantial portion of the chain into a fully stretched, single-file conformation. When applied to DNA, the fraction of genetic information sequentially ordered on the microwire surface will increase with the contour length, despite the increased entropy. This ability to handle long strands of DNA is in contrast to modern DNA sample preparation technologies for sequencing and mapping, which are typically restricted to comparatively short strands resulting in challenges in reconstructing the genome. Thus, in addition to discovering new rotation-induced macromolecular dynamics, this work inspires new approaches to handling genomic-length DNA strands.


INTRODUCTION
Although the double-helix structure of DNA stores genetic information linearly as a sequence of bases at the molecular level, entropy randomizes the threedimensional conformations of DNA polymers in free solution, making accessing genomic information from long polymers exceedingly challenging. To manipulate and study DNA, physical approaches have been developed to spatially arrange DNA in a controlled format, including magnetic and optical traps [1,2], microfluidic approaches [3][4][5][6], and nanoconfinement approaches [7], such as Convex Lens-induced Confinement (CLiC) [8] or nanopore confinement [9]. However, in vivo manipulation of genomic-length biopolymers remains an outstanding challenge that, with existing approaches, becomes more difficult with increasing strand length.
In this work, we untangle and order individual DNA molecules in long robust sections. Using fluorescence imaging, we demonstrate that DNA strands deform and elongate in shear flow near a rotating microwire when the rotation rate exceeds a critical value. While single strands of untethered DNA are observed to elongate substantially in agreement with scaling theory and simulations, full DNA extension becomes obtainable by combining tethering with this DNA spooling approach. We find that tethering one end of the DNA to the rotating microwire enhances extension, and produces novel "shofar" conformations. By slowly increasing the experimentally realizable rotation rate above an additional critical value, our simulations show a substantial portion of the tethered DNA is spooled into a strongly stretched, single-file curvilinear conformation. Since the conformation consists of a "shofar"-type tail and a base-ordered stem, we refer to this as a "French-horn" conformation. In this rotation-induced macromolecular spooling, the fraction of the polymer in the fully ordered stem of the French-horn conformation is found to increase with strand length.

Rotation-induced Spooling
The physical mechanism for compactly wrapping stretched DNA around a spinning microwire is presented in Fig. 1a for untethered DNA. The rotation drives an azimuthal flow, which shears the deformable DNA and draws it towards the surface of the microwire. When the DNA is untethered, the biopolymer deforms and elongates. While the rotation rates required to strongly stretch untethered DNA are very high moderate rotation rates are sufficient to strongly stretch significant portions of the DNA when one end is tethered to the rotating microwire (Figs. 1c and 1d). The stretched DNA is tautly wrapped and unraveled into a coordinated single-file conformation. By gradually increasing the rotation rate, the Simulation snapshot portraying a typical steady-state conformation of an untethered DNA chain. c. The tethered system with the French-horn DNA conformation: The blob size ξ is a function of angle φ from the untethered end. For sufficiently large Weissenberg numbers Wi, a portion of the DNA is fully stretched into a stem section, while the rest is in a horn conformation. For smaller Wi there is no stem and the polymer is in the shofar conformation. d. Simulation snapshots portraying typical steady-state tethered conformations. Increasing Wi changes the state from relaxed, to shofar, to French-horn conformations. Sufficiently long, rapidly rotated chains wrap around the circumference forming a rim of strongly stretched DNA.
fraction of segments in the strongly stretched stem is increased, while chain-breaking events can be avoided and the excluded volume interactions between the blobs ensure that overlaps are rare. In this way, substantial portions of single molecules of sufficiently long DNA can be organized into single-file conformations, deterministically ordering the genetic information on a cylindrical surface.
In prior work using planar geometries, DNA and other long macromolecules have been demonstrated to deform when subjected to experimentally achievable flows [10,11]. For example, a wall-tethered flexible polymer subject to sufficient shear rates deforms into a string of blobs [12,13]. The physical mechanisms leveraged in this study can be understood by considering the deformation of untethered DNA. We thus consider untethered strands, after first introducing our methods, before subsequently analyzing tethered chains.

METHODS
When a cylindrical microwire of radius a is rotated with rate Ω, the no-slip boundary condition generates a flow profile v φ (r) = Ωa 2 /r and a rapidly decaying shear rateγ (r) = −Ω (a/r) 2 at a distance r from the centre of the cylinder (Fig. 1). Each strand of DNA is composed of N Kuhn segments of length b ≈ 100 nm and characteristic relaxation time τ b = ηb 3 /k B T ≈ 2.4 × 10 −4 s in a solvent with dynamic viscosity η and thermal energy k B T . For a chain in a good solvent with a Flory exponent ν ≈ 3/5, unperturbed DNA has a relaxation time τ τ b N 3ν and a corresponding undeformed radius of gyration R g0 bN ν .
By suspending a microwire directly above a coverslip (Fig. 2a), submerging it in a drop of solution containing DNA, and gradually increasing its rotation rate, we experimentally observe the dynamics of DNA in a rotationally-induced shear flow. Far from the slowly rotating wire, the shear is small enough that untethered DNA appears relaxed and diffusive. Nearer to the microwire, the flow dominates and the DNA is advected. Some polymers, which are initially distributed evenly throughout the solution, are drawn towards the rotating microwire [14]. This radial migration across streamlines towards the microwire is expected to arise from the combined effects of the hydrodynamic interactions with the wire, the nonhomogeneous flow, and the decrease of diffusivity with stretching [15]. Within the advection regime, many strands are experimentally observed to be elongated by the shear flow (Fig. 2b). The Weissenberg number Wi =γτ expresses the competition between the characteristic shear forces and the polymer's intrinsic relaxation time τ . The shear deforms the DNA when Wi 1 (Fig. 2b and c). By substituting in the relaxation time and shear rate, we predict that deformation occurs for r/a (Ωτ b ) 1/2 N 3ν/2 . This hints at the benefit of our technique when dealing with long or potentially even genomic-length DNA, since the minimum rotation rate required to achieve deformation is Ω * Wi τ b N −3ν , which drops rapidly with DNA length.

Experimental Apparatus
The apparatus used to perform the experiments is shown in Fig. 3. A 50 µm diameter tungsten wire (Malin co.) is held and spun at both ends by two stepper motors (Pololu Robotics & Electronics part 1204). Each end of the wire is held within a short piece of fluorinated ethylene propylene (FEP) microfludic tubing with an inner diameter of 50 µm (IDEX Health & Science), each of which is mounted to one of the stepper motors using couplers machined from aluminum. The small inner diameter of the tubing serves to reduce microwire wobble. The mounting points of the wires do not exactly align, which results in minor wobbling about the rotation axis.
Each stepper motor is mounted on a z-axis-stage (Thorlabs part MS1S) which allows its height to be varied. This allows the vertical position of the wire above the sample coverslip, as well as the angle of the wire to the horizontal, to be set precisely. One of the stepper motors is also mounted on a horizontal stage which allows the distance between the motors to be varied, so that a controlled tension can be applied to the microwire. The motors and stages are mounted to a base machined from aluminum. The rotation of the motors is controlled by an Arduino microcontroller board. The current maximum rotation rate is not limited by the motors but rather by vibrations and wobbling that hinder accurate imaging of the DNA.
The apparatus mounts on an inverted fluorescence microscope (Nikon Ti-E), and experiments are performed with a 60× NA 1.0 water immersion objective (Nikon CFI Plan Apo VC 60XWI) and an Andor iXon3 EMCCD camera at an EM gain of 270. We estimate the depthof-field of our optical system to be about 1.2 µm. This means that, due to the curvature of the microwire, the entire length of an extended DNA molecule cannot be in perfect focus at once. The longest segment of DNA which can lie within the focal plane at once is about 16 µm for a 50 µm-diameter wire, so some sections of the measured T4-DNA molecules are inevitably slightly out of focus. In order to observe the DNA, the YOYO-1 dye (Section ) is excited using 1.5 mW from a 488 nm Coherent Sapphire laser.

DNA Preparation
For untethered experiments, T4-phage DNA (166 kbp, New England Biosciences) is used at a concentration of 0.1 µg/mL. The DNA are stained with YOYO-1 fluorescent dye (Life Technologies) at a ratio of 1 dye molecule per 10 base-pairs. YOYO-1 increases the DNA contour length at this staining ratio from 56.4 µm to 65(2) µm for T4-DNA [16]. Experiments are carried out in 0.5× tris-borate-EDTA (TBE) buffer, with 3 % β-mercaptoethanol (BME) added as an anti- c. Image of the microwire 1 s after it has stopped rotating: DNA strands have begun to relax to more entropic states. d. Conformational relaxation of a single DNA strand in the period of time immediately after the microwire has stopped rotating (see Movie S1). e. Quantitative measurements of the extension of a single DNA strand during relaxation, along with the fit and associated studentized residuals, used to determine the fully extended length of the polymer (t = 0). photobleaching agent. Before beginning untethered experiments, the microwire is first passivated using a solution of 10 % 55 kDa polyvinylpyrrolidone (PVP) in 0.5X TBE to reduce sticking of DNA to the glass coverslip and tungsten microwire.
For a solvent viscosity similar to that of water, the bulk relaxation time and radius of gyration of T4-DNA are τ = 0.74 s and R g0 = 1.46 µm [17] with approximately N ≈ 500 Kuhn segments. The Flory exponent for T4-DNA in 0.5× TBE (a good solvent) is ν = 0.558, which is intermediate between a Gaussian chain and a swollen excluded volume chain to account for the finite length and width of the double-stranded DNA [18].

Measuring Extension
To accurately measure the end-to-end lengths L of the untethered DNA in the experimental apparatus, the relaxation from the extended state is fit. The microwire rotates, extending the near-wire polymers, and then the rotation is suddenly halted. The DNA strands relax to a more highly entropic state with decreased extension and the initial length where t is the time since the rotation of the microwire has stopped and A is a proportionality constant that is allowed to vary when fitting. To estimate the extension when the wire is rotating, this function is fit to the relaxation of each polymer, and the initialextended length when the wire is still rotating L(t = 0) = 4R 2 g0 + A N b is extracted from the fit (Fig. 2e). The fit to Eq. 1 is found to be accurate for the range of rotation rates utilized in this study.
Each individual measurement of the elongation length of the DNA strand at time t requires locating the positions of its two ends. The uncertainty in the experimental measurement of each of these positions is dominated by the pixel size of the camera, corresponding to 267 nm. The uncertainty in each measurement is then δL = 533 nm, while the uncertainty in the fully extended length of the polymer is extracted from the fit to its relaxation. There is large variance in the extent of elongation as expected from results on tumbling DNA in simple shears [20,21].

DPD Algorithm
The rotation-induced macromolecular spooling does not rely on chemical details of DNA-it is a non-specific physical phenomenon. For this reason, a coarse-grain Dissipative Particle Dynamics (DPD) numerical model is used to simulate generic DNA molecules to extend our understanding of DNA dynamics in three dimensional azimuthal flows at various rotation rates with hydrodynamic interactions. The apparatus used to collect data on the rotationinduced macromolecular spooling of DNA. The microwire is situated directly above the cover slip using the z stages of the two motor mounts and the tension stage is adjusted so that the microwire is taut. Solution containing DNA is dropped onto the cover slip such that a portion of the wire is submerged. The motors are driven simultaneously by an Arduino microcontroller.
Both DNA monomers and fluid particles are represented by soft particles of mass m and size σ. The i th particle is described by its coordinates in continuousspace, which evolve in response to inter-particle forces with other particles (labeled j). In DPD, these are commonly pair-wise conservative forces F C ij with a soft repulsion with a cut off r c and strength a ij , dissipative forces F D ij with inter-particle frictional drag γ DPD , and random forces F R ij due to the thermal energy k B T . An extended discussion of these forces can be found in Ref. [22].
The DPD beads that form the DNA strand have an additional conservative bonding force to sequentially polymerize N = 100 beads. For this, we opt for the widely used FENE spring potential U FENE , in conjuction with the purely repulsive truncated Lennard-Jones potential U WCA [22]. We note that the both U FENE and U WCA are only applied between consecutive coarse-grained beads. The effective coarse-grained monomer size is determined from the U WCA and chosen to be the same as the DPD beads r c = 1. The DPD repulsion suffices to model longranged excluded-volume interactions, which yield the expected size exponent for a real polymer in a good solvent [23]. This allows us to obtain a monomer relaxation time 1 < τ b < 0.1 [23], which we approximate as τ b ≈ 0.3. The presence of the DPD fluid ensures that hydrodynamic interactions between monomers are reproduced, leading to Zimm dynamics. The bulk relaxation time of DPD chains kept in the vicinity of the surface of the microwire by tethering in quiescent solvent is measured to be τ = 920, corresponding to τ = 0.74 s when modelling T4-DNA. Likewise, each DPD bead represents roughly ten Kuhn segments of T4-DNA, since the DPD chain consists of N = 100 beads, while T4-DNA is composed of N ∼ 10 3 Kuhn segments. Such long DNA can be rescaled and coarse-grained as a freely-jointed chain on scales above the DPD bead size, while dynamics on smaller length scales are not modelled.
We numerically model the microwire as a semi-infinite cylindrical surface of radius a = 5 composed of three layers of DPD beads, which rotate with a fixed angular velocity. The rotation of these wall beads produces the shear rateγ that deforms the polymer, which hydrodynamically interacts with the cylinder surface through the DPD solvent. Relative to R g0 , the simulated cylinder is smaller than the experimentally employed tungsten wire, which may impact cross-streamline migration. A concentric cylinder of radius a out = 35 provides a finite size to the system and is constituted by a single layer of DPD beads that rotates with the theoretically expected speed v φ (r = a out ) = Ωa 2 /a out . The resulting shears match the theoretical expectations except very near the microwire surface, where minor packing and slip-effects cause the flow profile to be over predicted. We apply periodic boundary conditions along the microwire's axial direction with a system thickness of 15. We observe that the soft nature of the DPD beads can lead to slight polymer penetration into the inner rod when the rotation rate is sufficiently high. To reduce this effect, we add a smooth mathematical cylindrical surface that is transparent to the solvent beads (hence does not explicitly affect the flow) but unto which we apply U WCA for every polymer bead. With the surface characterized by a radius a − 1/2, the polymer feels a steric repulsion and does not penetrate the surface of the nominal rod.
The DNA model is initialized tightly wrapped around the cylinder and left to equilibrate. By noting the monomer position in a cylindrical coordinate system r i (r i , φ i , z i ), it is possible to find the instantaneous span L s ≡ r ∆φ = r (φ max − φ min ). A tightly wrapped polymer is thus initialized at a high value of L s , which exponentially decreases to a steady state value (characterized by τ exp ). An exponential fit to Eq. 1 is conducted to determine the steady-state curvilinear extension L. As in experiments, only DNA strands in the elongated phase of the tumbling cycle [21] are analyzed.
Although our simulations demonstrate that the monomers are driven towards the cylinder surface as expected, when the rotation rate is low, untethered polymers can sometimes escape into the bulk. This may be due to the fact that the simulated cylinder is effectively smaller than than in experiments, which may affect radial migration. Although re-capture events are observed, escaped polymers are not analyzed. As in simple shear, tumbling dynamics [24] and U-turn conformations [25] are observed in rotation-induced macromolecular spooling. This sets a lower bound to the rotation rate since the model DNA commonly escapes from the surface of the rotating wire at slow rotation rates.

UNTETHERED DNA
Weakly Extended Untethered DNA The curvilinear extension L of elongated polymers are experimentally measured at a number of rotation rates (Fig. 4). Below Ω ≈ 7 rpm, some strands are drawn towards the rotating microwire; however, their elongation is weak. In this weak elongation regime (Wi 1), the DNA behaves as a thermal spring subjected to hydrodynamic stretching forces. The incident shear rate around a rotating cylinder is expected to bė γ = Ωa 2 / (a + ξ) 2 Ω, which enacts a characteristic dimensionless drag F b/k B T Ωτ b N 2ν on the weakly deformed DNA. The drag force on T4-DNA at the microwire surface is roughly F ∼ ηR g0 ×(γR g0 ) ∼ 10 −2 pN. Since the dimensionless thermal spring force scales as L/ (N b), we predict that in the weak elongation regime which can also be written L/R g0 WiN 1−2ν and is independent of wire diameter. The DNA extension is thus predicted to increase rapidly with rotation rate as L ∼ Wi in this regime, which is consistent with the low rotation rate data in Fig. 4. However, a crossover from an initially rapid rise for Wi 1 to a much slower rise when Wi > 1 is expected; the experimental data shows a marked crossover from the weak-extension regime to a much lower scaling at moderate rotation rates.

Moderately Extended Untethered DNA
When the rotation rate is moderate (Wi 1), the DNA deforms and aligns along the azimuthal streamlines at the surface of the rotating cylinder (Fig. 2c-d). For a limited regime, this moderate elongation can be roughly regarded as a string of blobs, in which the statistical breadth of the deformed chain (transverse to the flow direction) is described by the characteristic size ξ = k B T /F (Fig. 1b). This moderate elongation regime is bordered by the weakly extended regime discussed in Section for forces that do not deform the conformation (k B T /R g0 ∼ 3 × 10 −3 pN) and by the upper limit of strong deformation, in which the blob framework breaks down as ξ approaches the Kuhn length (k B T /b ∼ 6 × 10 −2 pN). While the upper limit is a mechanical property of double stranded DNA, the lower limit decreases as R g0 ∼ N ν increases. In this limited moderate extension regime, genomic-length DNA of many Kuhn segments can be considered a procession of blobs, in which tension is equally distributed throughout the chain and each blob has a theoretical relaxation time τ ξ τ b (ξ/b) 3 . Since shear and relaxation balance when the Weissenberg number for each is unity, the scaling prediction for the untethered characteristic size scale is This result can be written as ξ/R g0 Wi −1/3 . The conformational breadth of the DNA transverse to the flow direction is experimentally inaccessible due to the diffraction limited resolution. Therefore, we consult numerical experiments to directly investigate the conformational breadth. DPD simulations explicitly estimate size as the average distance of the monomers from the surface of the wire ξ = r − a (Fig. 4 inset). For sufficient rotation rates, ξ scales with rotation rate in the manner predicted by Eq. 3. At low rotation rates, the DNA often diffuses into bulk solution from the wrapped state, placing a lower bound on the rotation rates that can be numerically investigated (see Section ).
The curvilinear end-to-end distance of the highly deformed DNA on the wire surface is predicted to be by considering the number of Kuhn segments in a single blob and the total number of blobs. This result is equivalent to L/R g0 Wi (1−ν)/(3ν) . The extension increases linearly with contour length N b such that that longer chains stretch more easily than shorter chains in this moderate extension regime, but the extension is only weakly dependent on rotation rate in this moderately extended regime. Inserting ν = 0.558 for T4-DNA predicts a weak dependence on rotation rate (L/R g0 ∼ Wi 0.26 ) for this regime, which stands in contrast to the initial linear dependence in Eq. 2 for weakly stretched chains. Though the effective Flory exponent may vary slightly for blobs consisting of fewer Kuhn lengths, since double stranded DNA is neither sufficiently thin nor stiff to be a Gaussian chain (ν = 1/2) the double stranded minimum value (ν = 0.535 [18]) is only 4% less than the full chain value for T4-DNA (ν = 0.558 [18]).
At the highest experimentally achievable rotation rates in Fig. 4, the extension of the untethered DNA scales in accordance with Eq. 4. Due to limitations on the current apparatus, the microwire cannot be steadily driven above Ω ≈ 130 rpm without inducing significant wire wobble. DPD simulations further support this scaling at moderate rotation rates. The coarse-grained simulations replicate the experimental scaling in the moderate extension regime described by Eq. 4.
Although the moderate extension regime is evident in both the experimental and numerical untethered results (Fig. 4), relatively large experimental rotation rates are required to deform the T4-DNA substantially more. The DPD data extend to higher Wi and exhibit evidence of a strong stretching limit by deviating from Eq. 4 (Fig. 4). In this highly extended regime, the dimensionless restoring force goes as which is consistent with the high Wi DPD data. In this strong force limit, the coarse-grained DPD simulations are not expected to quantitatively reproduce the semiflexible nature of double stranded DNA; however, we expect the qualitative conclusion that the curvilinear extension slows at high Weissenberg number to hold for double stranded DNA. To experimentally obtain such strongly stretched chains would require ξ → b and, therefore, Ω * full τ −1 b ≈ 4000 s −1 , which would be experimentally challenging. However, by anchoring one end of the DNA to the rotating microwire, a substantial portion of sufficiently long chains can be strongly stretched.

TETHERED DNA
Untethered strands of DNA readily deform in the vicinity of the rotating microwire, but reaching the fully stretched limit is infeasible and fluctuations due to tumbling dynamics are significant. However, by anchoring one end of the DNA to the rotating microwire, a substantial portion of sufficiently long chains can be strongly stretched.
For the tethered experiments, λ-phage DNA (48.5 kbp, New England Biosciences) is used, also at a concentration of 0.1 µg/mL and stained with YOYO-1 fluorescent dye as in the untethered case. Staining increases the λ-DNA contour length from 16.5 µm to 19.0(7) µm [16]. Custom Cy5-labeled oligos complementary to the 12-bp singlestranded overhangs at the ends of λ-phage DNA are added (IDT). The hydrophobic character of the Cy5 endlabels results in a significant increase in surface-ahesion. The glass coverslip is passivated with PVP as in the untethered case, but the microwire is not. Tethered DNA molecules, for which one end has attached to the wire, are used for quantitative analysis of the tethered case. Since the hydrophobic Cy5-labels are positioned at the DNA ends, immobilization occurs at the ends. In DPD simulations, the polymer is tethered at one end via a FENE potential.
Similar to the classical differences between a chain dragged by one end through a quiescent fluid and a wall-tethered chain subjected to shear [11,26], the tethered system setup differs from a chain attached to a filament [27] in that it is the shearing flow due to rotation that leads to deformation. At low rotation rates (Ω < Ω * Wi ), the tethered DNA is expected to remain relaxed and simply rotate with the anchoring point. Simulation snapshots at low rotation rates ( Fig. 1d; Wi 10 −0.3 ) show that the tethered DNA molecules follow the rotation in a quasi-static manner. They remain in a thermalized state throughout the complete simulation runs because the shear is too low to stretch the DNA into an outof-equilibrium conformation. The critical rotation rate for deformation is predicted to scale as Ω * Wi τ b N −3ν , which is the same as in the untethered state and is confirmed by DPD simulations of tethered DNA (Fig. 5a).

Shofar Conformation
For greater rotation rates, the tethered polymer deforms into a conformational profile that decreases in size as a function of the angle φ from the free end (Fig. 1c). At a rotation rate Ω ≈ 130 rpm a configuration of varying transverse conformational breadth is readily achievable and an example of the growth of the conformational size away from the tethering point is qualitatively clear in Fig. 5b. If the planar case may be described as a "trumpet" conformation [28] then our experiments produce a cylindrical "shofar" conformation.
DPD snapshots also demonstrate this shofar conformation (Fig. 1d), showing that a moderate rotation rate is sufficient to deform a long polymer and wrap it around the microwire. There remain visible thermal fluctuations along the backbone, and the conformational size can be seen to increase with the distance to the tethering point, as in the experiments shown in Fig. 5b. The simulations can quantitatively measure the statistical shape of the shofar conformation ξ (n) as a function of segment index from the free-end (Fig. 5c). A free-end regime is clearly identifiable at small n in numerical measurements of the transverse conformational size ξ, while at moderate n the size decreases rapidly until the tethering point is approached. Although the range is limited by the contour length of the simulated chain, the mid-region conformation of the simulated chain appears to exhibit a rapid decaying scaling law.
In order to understand the observed shofar conformation, we extend the scaling theory for highly extended untethered DNA (Eq. 3) to a tethered chain, in which the tension in each Kuhn segment varies with the distance from the free end. Temporarily, consider a series of well-defined blobs for a sufficiently long genomic-length polymer: The i th blob is described by the radial distance r i = a + ξ i and arc length s i = r i φ (Fig. 1c). The drag force on each is f i ηv φ (r i ) ξ i but the tension on each segment is due to the total drag of all the preceding blobs. However, the tension does not vary stepwise through a procession of well-defined blobs but rather varies continuously from segment to segment as j f j = j ηv φ ξ j → ηv φ ξ (ds/ξ), which is still expected to be locally balanced by the thermal force k B T /ξ. Therefore, for sufficiently long polymers, the scaling theory predicts that the transverse conformational profile as a function of angle φ from the free end is This scaling result reproduces the classic planar-shear result when the shear rateγ replaces Ω and the linear distance replaces the arc length s = aφ [26,29]. Assuming the tension varies slowly along the chain, we can write the rate of change of the local extension s with segment index n counted from the free end to be ds/dn ≈ b (b/ξ) (1−ν)/ν [26], or ds/dn ≈ (a + ξ) (dφ/dn) by inverting Eq. 6. Assuming that a b and considering n sufficiently far from the free end, we equate the two forms and integrate to find (1−ν)/(3ν−1) n 2ν/(3ν−1) . Spooling tethered DNA strands. a. Simulations of the critical rotation rate for which deformation of tethered DNA occurs. The solid line shows the predicted scaling of Ω * Wi ∼ N −3ν . b. Conformation of a single λ-phage DNA strand in the period immediately after the microwire has stopped rotating at Ω = 130 rpm, showing the relaxation of the polymer (see Supplemental Video 2). The end of the polymer that appears towards the top of the image is tethered to the wire. The conformational size is observed to increase steadily away from the tethering point, as predicted for a tethered polymer in the shofar configuration. c. Frenchhorn conformations seen via DPD blob size ξ as a function of segment index from the free-end. Various rotation rates are collapsed according to Eq. 8 with blob size normalized by the free-end blob size ξ1. Dotted line shows scaling prediction ξ ∼ n −ν/(3ν−1) ∼ n −0.77 . Diamonds ( ) denote the predicted transition from shofar to stem at n * (τ b Ω) −1 . Inset. Wrapping angle from the free end, collapsed via Eq. 7. Dashed line shows the scaling φ ∼ n 2ν/(3ν−1) ∼ n 1.54 .
Within the shofar conformation and sufficiently far from the free end, simulations show that the angle from the free end grows with segment index as described by Eq. 7 (Fig. 5c inset; dashed line). As expected, the predicted scaling does not hold near the free end.
Further details about the conformation can be obtained from the blob size as a function of segment index from the free end. Substituted into Eq. 6, the angle ξ (n) produces which can be written ξ (n) /R g0 (Wi n/N ) −ν/(3ν−1) . As in the planar case, the blob framework (Eq. 8) predicts that the conformational profile diverges as n → 0. Therefore, it is truncated at the first blob φ 1 ≈ ξ 1 /a. This predicts ξ 1 /b (Ωτ b ) −1/3 , which is the same as the untethered conformational size (Eq. 3). The conformation crosses to the shofar conformation that was observed in the inset and is seen to be well described by Eq. 8 ( Fig. 5c; dotted line) for segments in the mid-region of the simulated chain. Within the blob framework, the conformation at the tethering point is predicted to be . However, as the tethered end is approached, the conformational breadth as a function of index begins to saturate to the segment size for many of rotation rates in the simulations (Fig. 5c). These curves represent a new conformation.

French-horn Conformation
If the rotation rate is increased to Ω * horn τ b N −1 then the conformational profile at the tethering point is predicted to approach the Kuhn size b and the blob-scaling regime ends. For higher rotation rates (Ω > Ω * horn ), the segments near the tethering point are strongly stretched and the chain adopts a "French-horn" conformation with a strongly stretched "stem" near the tethering point followed by a shofar-shaped horn conformation as the free end is approached (Fig. 1c). The onset of this Frenchhorn conformation occurs at substantially slower rotation rates than strongly stretched conformations of unthethered polymers. The segments within the stem and indeed those for which the characteristic blob size is not much larger than the double stranded DNA Kuhn length (ξ (n) b) may not be in the scaling regime. However, sufficiently far from from the stem the conformation of the "horn" portion can still be described by blob theory, and the segment index where the transition from shofar to stem occurs is n * Ω * stem /Ω (τ b Ω) −1 , and so the fraction of segments in the strongly stretched stem is Since the tension in the chain arises from the cumulated drag of the shofar conformation, obtaining λ → 1 requires unphysical rotation rates (Ω → Ω stem = τ −1 b ). Likewise, even strongly stretching 90% of the chain demands N Ωτ b = 10; however, strongly stretching half of the DNA requires a rotation rate that is decreased by nearly an order of magnitude to only N Ωτ b = 2. Due to limitations on the current apparatus (≈ 130 rpm maximum experimental rotation rate), the stem of the French-horn conformation has not yet been experimentally observed. Though the relaxation time of a single Kuhn segment (τ b ∼ 10 −4 s) is quite short, N can be extremely large for DNA, such that experimentally obtainable rotation rates are predicted to strongly stretch significant portions of genomic-length biopolymers.
A numerical rotation rate of Ω = 0.1 suffices to wrap the majority of the monomers in a single-file manner in simulations ( Fig. 1d; Wi 10 2 ). Here, the stemstate dominates with a small proportion of monomers forming a more relaxed conformation as the free end is approached. Sufficiently long chains encircle the microwire with a strongly stretched and ordered rim of genomic information. In simulations, the stem portion does not reach a constant cutoff due to the softness of the DPD beads and slight bond stretching, which results in a reduced but non-zero decrease in effective blob size (Fig. 5c). As in the untethered case, in the strongly stretched stem the DPD simulations provide a qualitative picture but not quantitative predictions for semiflexible double stranded DNA. However, the predicted crossover from the shofar to the stem of the French-horn conformation at n * is seen to approximate this point (Fig. 5c; ).
Our results reproduced the French-horn conformation for tethered chains, providing confidence that a significant portion of the DNA can be expected to be strongly stretched. Future research on rotationinduced macromolecular spooling should perform simulations with semiflexible polymer models that represent DNA dynamics at length scales below the Kuhn length in order to explore the stem region of the French-horn conformation in more detail. To ensure that the strongly stretched stem follows a deterministic screw-like, overlapfree organization, the rotation rate should be slowly increased. In this way, excluded volume effects act to prevent overlap such that, when the stem forms, it is stretched into a non-overlapping single-file conformation, an important requirement for sequential ordering. While genetic information in the shofar portion of the DNA remains disordered in space, it is arranged cylindrically in the stem, just as linearly advancing music notes are punched on the cylinder of a music box.

DISCUSSION
Organization of long, single strands of DNA sequentially in space must overcome the high entropy as-sociated with complex, folded polymer conformations. Rotation-induced macromolecular spooling presents a new single-molecule manipulation concept that explicitly exploits disordering thermal fluctuations through rotation-induced hydrodynamic drag, allowing DNA to be ordered into a curvilinear progression of base-pairs, arranged sequentially on the surface of a rotating microwire.
We observe untethered strands to be highly deformed, and cross streamlines towards the microwire, essentially generating a system capable of studying the singlemolecule origin of the Weissenberg effect in rod-climbing of polymeric liquids [30]. We have further found enhanced extension to be possible with tethering. It is compelling that DNA tethered to the rotating microwire takes various previously-unseen conformational states, including the experimentally observed shofar and computationally observed French-horn conformations, at moderate rotations rates. The portion of the French-horn conformation that is in the fully extended stem state increases with both rotation rate and contour length. Longer strands, such as genomic-length DNA, are seen to have a greater portion of their segments in the strongly stretched stem regime of the French-horn conformation. In fact, by anchoring one end of a DNA strand to the microwire after utilizing DNA ligase to conjoin a long "drag-tag" of additional DNA (that is not of interest) to the free end to act as the unordered shofar-shaped portion of the French-horn conformation, it is possible that future applications could tautly wrap entire DNA lengths into a deterministic single-file and unknotted stem conformation encircling and forming a genomic rim around the microwire.
Since the single-file stem conformation results from shear drag on the shofar end, rotation-induced macromolecular spooling of very long strands of fragile DNA may face limitations due to shear-induced fragmentation [31][32][33], overstretching [34][35][36] or tether failing. Within the stem, the tension is largest and the rate of tether breakage events may be increased or overstretching may occur at high rotation rates. Although we do not observe instances of fragmentation of tethered DNA in the current implementation, we do see tethered strands disappear after repeated spinning, indicating that the tether has failed. We speculate that gently ramping up the rotation rate will both reduce sudden fragmentation and allow excluded-volume interactions to reduce instances of entangled conformations or overlaps within the strongly stretched stem.
Furthermore, a subsequent experimental step could be implemented to immobilize DNA on the microwire, following its ordering on the surface, for storage and readout. In this way, the complete polymer could be organized into a linear progression of base-pairs that are not susceptible to thermal fluctuations as random coils. While the relatively short viral genomic DNA considered here (∼ 10 2 kbp; ∼ 10 1 µm) may be fully stretched in linear geometries at micron-scales, the distribution of genome sizes varies enormously [37], such that linearly stretching DNA is an impossibility at micronscales for many genomes -an impossibility rotationinduced spooling could elegantly overcome through compact curvilinear sequential organization. The genomic sequence could be stored on the wire for long periods and at high density, in analogy to punching sequential music notes on the surface of a music box cylinder.