Equilibrium distributions of topological states in circular DNA: Interplay of supercoiling and knotting

Two variables define the topological state of closed double-stranded DNA: the knot type, K, and ΔLk, the linking number difference from relaxed DNA. The equilibrium distribution of probabilities of these states, P(ΔLk, K), is related to two conditional distributions: P(ΔLk|K), the distribution of ΔLk for a particular K, and P(K|ΔLk) and also to two simple distributions: P(ΔLk), the distribution of ΔLk irrespective of K, and P(K). We explored the relationships between these distributions. P(ΔLk, K), P(ΔLk), and P(K|ΔLk) were calculated from the simulated distributions of P(ΔLk|K) and of P(K). The calculated distributions agreed with previous experimental and theoretical results and greatly advanced on them. Our major focus was on P(K|ΔLk), the distribution of knot types for a particular value of ΔLk, which had not been evaluated previously. We found that unknotted circular DNA is not the most probable state beyond small values of ΔLk. Highly chiral knotted DNA has a lower free energy because it has less torsional deformation. Surprisingly, even at |ΔLk| > 12, only one or two knot types dominate the P(K|ΔLk) distribution despite the huge number of knots of comparable complexity. A large fraction of the knots found belong to the small family of torus knots. The relationship between supercoiling and knotting in vivo is discussed.

The Malthusian parameter of ascents: What prevents the exponential increase of one’s ancestors?

The reason that the indefinite exponential increase in the number of one’s ancestors does not take place is found in the law of sibling interference, which can be expressed by the following simple equation:\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\begin{matrix}{\mathit{N}}_{{\mathit{n}}} \enskip & \\ {\mathit{{\blacksquare}}} \enskip & \\ {\mathit{ASZ}} \enskip & \end{matrix} {\mathrm{\hspace{.167em}{\times}\hspace{.167em}2\hspace{.167em}=\hspace{.167em}}}{\mathit{N_{n+1},}}\end{equation*}\end{document} where Nn is the number of ancestors in the nth generation, ASZ is the average sibling size of these ancestors, and Nn+1 is the number of ancestors in the next older generation (n + 1). Accordingly, the exponential increase in the number of one’s ancestors is an initial anomaly that occurs while ASZ remains at 1. Once ASZ begins to exceed 1, the rate of increase in the number of ancestors is progressively curtailed, falling further and further behind the exponential increase rate. Eventually, ASZ reaches 2, and at that point, the number of ancestors stops increasing for two generations. These two generations, named AN SA and AN SA + 1, are the most critical in the ancestry, for one’s ancestors at that point come to represent all the progeny-produced adults of the entire ancestral population. Thereafter, the fate of one’s ancestors becomes the fate of the entire population. If the population to which one belongs is a successful, slowly expanding one, the number of ancestors would slowly decline as you move toward the remote past. This is because ABZ would exceed 2. Only when ABZ is less than 2 would the number of ancestors increase beyond the AN SA and AN SA + 1 generations. Since the above is an indication of a failing population on the way to extinction, there had to be the previous AN SA involving a far greater number of individuals for such a population. Simulations indicated that for a member of a continuously successful population, the AN SA ancestors might have numbered as many as 5.2 million, the AN SA generation being the 28th generation in the past. However, because of the law of increasingly irrelevant remote ancestors, only a very small fraction of the AN SA ancestors would have left genetic traces in the genome of each descendant of today.

Invariant size–frequency distributions along a latitudinal gradient in marine bivalves

In the most extensive analysis of body size in marine invertebrates to date, we show that the size–frequency distributions of northeastern Pacific bivalves at the provincial level are surprisingly invariant in modal and median size as well as size range, despite a 4-fold change in species richness from the tropics to the Arctic. The modal sizes and shapes of these size–frequency distributions are consistent with the predictions of an energetic model previously applied to terrestrial mammals and birds. However, analyses of the Miocene–Recent history of body sizes within 82 molluscan genera show little support for the expectation that the modal size is an evolutionary attractor over geological time.

Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations

We describe and test a Markov chain model of microsatellite evolution that can explain the different distributions of microsatellite lengths across different organisms and repeat motifs. Two key features of this model are the dependence of mutation rates on microsatellite length and a mutation process that includes both strand slippage and point mutation events. We compute the stationary distribution of allele lengths under this model and use it to fit DNA data for di-, tri-, and tetranucleotide repeats in humans, mice, fruit flies, and yeast. The best fit results lead to slippage rate estimates that are highest in mice, followed by humans, then yeast, and then fruit flies. Within each organism, the estimates are highest in di-, then tri-, and then tetranucleotide repeats. Our estimates are consistent with experimentally determined mutation rates from other studies. The results suggest that the different length distributions among organisms and repeat motifs can be explained by a simple difference in slippage rates and that selective constraints on length need not be imposed.

The estimation of statistical parameters for local alignment score distributions

The distribution of optimal local alignment scores of random sequences plays a vital role in evaluating the statistical significance of sequence alignments. These scores can be well described by an extreme-value distribution. The distribution’s parameters depend upon the scoring system employed and the random letter frequencies; in general they cannot be derived analytically, but must be estimated by curve fitting. For obtaining accurate parameter estimates, a form of the recently described ‘island’ method has several advantages. We describe this method in detail, and use it to investigate the functional dependence of these parameters on finite-length edge effects.

Characteristics of the deep ocean carbon system during the past 150,000 years: ΣCO2 distributions, deep water flow patterns, and abrupt climate change

Studies of carbon isotopes and cadmium in bottom-dwelling foraminifera from ocean sediment cores have advanced our knowledge of ocean chemical distributions during the late Pleistocene. Last Glacial Maximum data are consistent with a persistent high-ΣCO2 state for eastern Pacific deep water. Both tracers indicate that the mid-depth North and tropical Atlantic Ocean almost always has lower ΣCO2 levels than those in the Pacific. Upper waters of the Last Glacial Maximum Atlantic are more ΣCO2-depleted and deep waters are ΣCO2-enriched compared with the waters of the present. In the northern Indian Ocean, δ13C and Cd data are consistent with upper water ΣCO2 depletion relative to the present. There is no evident proximate source of this ΣCO2-depleted water, so I suggest that ΣCO2-depleted North Atlantic intermediate/deep water turns northward around the southern tip of Africa and moves toward the equator as a western boundary current. At long periods (>15,000 years), Milankovitch cycle variability is evident in paleochemical time series. But rapid millennial-scale variability can be seen in cores from high accumulation rate series. Atlantic deep water chemical properties are seen to change in as little as a few hundred years or less. An extraordinary new 52.7-m-long core from the Bermuda Rise contains a faithful record of climate variability with century-scale resolution. Sediment composition can be linked in detail with the isotope stage 3 interstadials recorded in Greenland ice cores. This new record shows at least 12 major climate fluctuations within marine isotope stage 5 (about 70,000–130,000 years before the present).