We use Monte-Carlo simulations to study the effect of variable rigidity on plectoneme formation and localization in supercoiled dsDNA. We show that the presence of soft sequences increases the number of plectoneme branches and that the edges of the branches tend to be localized at these sequences. We propose an experimental approach to test our results in vitro, and discuss the possible role played by plectoneme localization in the search process of transcription factors for their targets (promoter regions) on the bacterial genome.The conformational properties of double-stranded (ds) DNA molecules are usually modeled by treating these biopolymers as semi-flexible chains with uniform rigidity that can be represented by a single persistence length (rigidity is proportional to persistence length), l p 50nm in physiological conditions (1, 2). This is justified when one is interested in large-scale properties of ds-DNA for which one can replace the sequence-dependent distribution of the elastic rigidity by its average over the chain. When one is interested in small and intermediatescale phenomena one has to consider the full sequencedependence of the rigidity. While some studies suggested that the rigidity of bare dsDNA varies across a limited range 40nm < l p < 75nm (3, 4), experiments on cyclization of short dsDNA fragments (≈ 100bp) reported much higher cyclization ratios than expected (5, 6). This led to the proposal of DNA kinks -pointlike highly flexible domains (7) -perhaps due to formation of "DNA bubbles" (8, 9). Such sequence-dependent rigidity is a property of bare dsDNA and it has been suggested that the effect can be utilized for the design of promoter sequences in order to control the DNA binding affinity of transcription factors that are sensitive to DNA bendability (10). It may also arise as the consequence of the binding of proteins to specific DNA sequences; indeed, in vivo DNA is partially covered by proteins that affect its flexibility (3,(11)(12)(13)(14) and/or its local curvature (15,16). For example, RecA bacterial proteins polymerize along DNA to give an effective persistence length of hundreds of nm to the 14). Other positivelycharged proteins (e.g., HMGB) and polyamines, increase DNA's flexibility significantly (11,17,18). Thus, the variability of DNA rigidity may be even higher in vivo than that of bare DNA in vitro.In this work we study the interplay between local rigidity and plectoneme localization in supercoiled dsDNA. When circular DNA is subjected to sufficiently large torsional stress, the minimization of the free energy yields strongly-writhed conformations known as plectonemes (see e.g. Fig. 1). We show that the number of plectonemic branches and the locations of the edges (end loops of the branches in Fig. 1) are affected by nonuniform DNA rigidity. This applies not only to circular chains (e.g., bacterial and mitochondrial DNA, and DNA p occupies 1/3 of the chain and is divided into 4 equally-spaced domains of identical size, l (h) p occupies the remaining 2/3 of the chain, and t...