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.

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).

Distinct desmocollin isoforms occur in the same desmosomes and show reciprocally graded distributions in bovine nasal epidermis.

The adhesive core of the desmosome is composed of cadherin-like glycoproteins of two families, desmocollins and desmogleins. Three isoforms of each are expressed in a tissue-specific and developmentally regulated pattern. In bovine nasal epidermis, the three desmocollin (Dsc) isoforms are expressed in overlapping domains; Dsc3 expression is strongest in the basal layer, while Dsc2 and Dsc1 are strongly expressed in the suprabasal layers. Herein we have investigated whether different isoforms are assembled into the same or distinct desmosomes by performing double immunogold labeling using isoform-specific antibodies directed against Dsc1 and Dsc3. The results show that individual desmosomes harbor both isoforms in regions where their expression territories overlap. Quantification showed that the ratio of the proteins in each desmosome altered gradually from basal to immediately suprabasal and upper suprabasal layers, labeling for Dsc1 increasing and Dsc3 decreasing. Thus desmosomes are constantly modified as cells move up the epidermis, with continuing turnover of the desmosomal glycoproteins. Statistical analysis of the quantitative data showed a possible relationship between the distributions of the two isoforms. This gradual change in desmosomal composition may constitute a vertical adhesive gradient within the epidermis, having important consequences for cell positioning and differentiation.