Orthogonal combinatorial mutagenesis: a codon-level combinatorial mutagenesis method useful for low multiplicity and amino acid-scanning protocols

We describe here a method to generate combinatorial libraries of oligonucleotides mutated at the codon-level, with control of the mutagenesis rate so as to create predictable binomial distributions of mutants. The method allows enrichment of the libraries with single, double or larger multiplicity of amino acid replacements by appropriate choice of the mutagenesis rate, depending on the concentration of synthetic precursors. The method makes use of two sets of deoxynucleoside-phosphoramidites bearing orthogonal protecting groups [4,4′-dimethoxytrityl (DMT) and 9-fluorenylmethoxycarbonyl (Fmoc)] in the 5′ hydroxyl. These phosphoramidites are divergently combined during automated synthesis in such a way that wild-type codons are assembled with commercial DMT-deoxynucleoside-methyl-phosphoramidites while mutant codons are assembled with Fmoc-deoxynucleoside-methyl-phosphoramidites in an NNG/C fashion in a single synthesis column. This method is easily automated and suitable for low mutagenesis rates and large windows, such as those required for directed evolution and alanine scanning. Through the assembly of three oligonucleotide libraries at different mutagenesis rates, followed by cloning at the polylinker region of plasmid pUC18 and sequencing of 129 clones, we concluded that the method performs essentially as intended.

Intron “sliding” and the diversity of intron positions

Alignments of homologous genes typically reveal a great diversity of intron locations, far more than could fit comfortably in a single gene. Thus, a minority of these intron positions could be inherited from a single ancestral gene, but the larger share must be attributed to subsequent events of intron gain or intron “sliding” (movement from one position to another within a gene). Intron sliding has been argued from cases of discordant introns and from putative spatial clustering of intron positions. A list of 32 cases of discordant introns is presented here. Most of these cases are found to be artefactual. The spatial and phylogenetic distributions of intron positions from five published compilations of gene data, comprising 205 intron positions, have been examined systematically for evidence of intron sliding. The results suggest that sliding, if it occurs at all, has contributed little to the diversity of intron positions.

Population dynamics, demographic stochasticity, and the evolution of cooperation

A basic evolutionary problem posed by the Iterated Prisoner’s Dilemma game is to understand when the paradigmatic cooperative strategy Tit-for-Tat can invade a population of pure defectors. Deterministically, this is impossible. We consider the role of demographic stochasticity by embedding the Iterated Prisoner’s Dilemma into a population dynamic framework. Tit-for-Tat can invade a population of defectors when their dynamics exhibit short episodes of high population densities with subsequent crashes and long low density periods with strong genetic drift. Such dynamics tend to have reddened power spectra and temporal distributions of population size that are asymmetric and skewed toward low densities. The results indicate that ecological dynamics are important for evolutionary shifts between adaptive peaks.

Polyploid formation created unique avenues for response to selection in Gossypium (cotton)

A detailed restriction fragment length polymorphism map was used to determine the chromosomal locations and subgenomic distributions of quantitative trait loci (QTLs) segregating in a cross between cultivars of allotetraploid (AADD) Gossypium hirsutum (“Upland” cotton) and Gossypium barbadense (“Sea Island,” “Pima,” or “Egyptian” cotton) that differ markedly in the quality and quantity of seed epidermal fibers. Most QTLs influencing fiber quality and yield are located on the “D” subgenome, derived from an ancestor that does not produce spinnable fibers. D subgenome QTLs may partly account for the fact that domestication and breeding of tetraploid cottons has resulted in fiber yield and quality levels superior to those achieved by parallel improvement of “A” genome diploid cottons. The merger of two genomes with different evolutionary histories in a common nucleus appears to offer unique avenues for phenotypic response to selection. This may partly compensate for reduction in quantitative variation associated with polyploid formation and be one basis for the prominence of polyploids among extant angiosperms. These findings impel molecular dissection of the roles of divergent subgenomes in quantitative inheritance in many other polyploids and further exploration of both “synthetic” polyploids and exotic diploid genotypes for agriculturally useful variation.

Theoretical investigation of the [1,2]-sigmatropic hydrogen migration in the mechanism of oxidation of 2-aminobenzoyl-CoA by 2-aminobenzoyl-CoA monooxygenase/reductase

The flavin hydroperoxide at the active site of the mixed-function oxidase 2-aminobenzoyl-CoA monooxygenase/reductase (Azoarcus evansii) transfers an oxygen to the 5-position of the 2-aminobenzoyl-CoA substrate to provide the alkoxide intermediate II−. Hydrogen migration from C5 to C6 follows this monooxygenation. The nature of the monooxygenation intermediate and plausible competing reactions leading to hydrogen migration have been considered. Ab initio molecular orbital theory has been used to calculate structures and electron distributions in intermediate and transition state structures. Electrostatic potential surface calculations establish that the transition state and product, associated with the C5 to C6 hydrogen transfer, are stabilized by electron distribution to the benzoyl-CoA thioester carbonyl oxygen. This is not so for the transition state and product associated with hydrogen transfer from C5 to C4. The activation energy for the 5,6-shift is 2.5 kcal/mol lower than that for the 5,4-shift. In addition, the product of the hydrogen 5,6-shift is more stable than is the product of the hydrogen 5,4-shift, by ≈6 kcal/mol. These results explain why only the shift of hydrogen from C5 to C6 is observed experimentally. Oxygen transfer and hydrogen migration almost coincide in the gas phase (activation energy of ≈0.6 kcal/mol, equivalent to a single bond vibration). Enzymatic formation of alkoxide II− requires its stabilization; thus, the rate constant for its breakdown would be slower than in the gas phase.

Characterization of Cell–Matrix Adhesion Requirements for the Formation of Fascin Microspikes

Cell adhesion to thrombospondin-1 (TSP-1) correlates with assembly of cell–substratum contact structures that contain fascin microspikes. In this analysis, cell-matrix requirements for assembly of fascin microspikes were examined in detail. In six cell lines, cell spreading on a TSP-1 substratum correlated with expression of fascin protein and formation of fascin microspikes. Microspikes were not formed by H9c2 cells adherent on fibronectin, vitronectin, collagen IV, or platelet factor 4. However, both fascin microspikes and focal contacts were assembled by cells adherent on laminin-1. Using mixed substrata containing different proportions of TSP-1, and fibronectin, fascin microspike formation by H9c2 and C2C12 cells was found to be reduced on substrata containing 25% fibronectin and abolished on substrata containing 75% fibronectin. Adhesion to intermediate mixtures of TSP-1 and fibronectin resulted in coassembly of fascin microspikes and focal contacts, colocalization of fascin with actin stress fiber bundles and altered distributions of β1 integrins, cortical α-actinin, and tropomyosin. In cells adherent on 50% TSP-1:50% fibronectin, GRGDSP peptide treatment decreased focal contact assembly and altered cytoskeletal organization but did not inhibit microspike assembly. Treatment with chondroitin sulfate A or p-nitrophenol β-d-xylopyranoside decreased microspike formation and modified cytoskeletal organization but did not inhibit focal contact formation. In polarized migratory and postmitotic C2C12 cells, fascin microspikes and ruffles were localized at leading edges and TSP matrix deposition was also concentrated in this region. Depletion of matrix TSP by heparin treatment correlated with decreased microspike formation and cell motility. Thus, the balance of adhesive receptors ligated at the cell surface during initial cell–matrix attachment serves to regulate the type of substratum adhesion contact assembled and subsequent cytoskeletal organization. A role for fascin microspikes in cell motile behavior is indicated.

Particles move along actin filament bundles in nerve growth cones.

Organelle movement along actin filaments has been demonstrated in dissociated squid axoplasm [Kurznetsov, S. A., Langford, G.M. & Weiss, D. G. (1992) Nature (London) 356, 722-725 and Bearer, E.L., DeGiorgis, J.A., Bodner, R.A., Kao, A.W. & Reese, T.S. (1993) Proc. Natl. Acad. Sci. USA 90, 11252-11256] but has not been shown to occur in intact neurons. Here we demonstrate that intracellular transport occurs along actin filament bundles in intact neuronal growth cones. We used video-enhanced differential interference contrast microscopy to observe intracellular transport in superior cervical ganglion neurons cultured under conditions that enhance the visibility of actin bundles within growth cone lamellipodia. Intracellular particles, ranging in size from < 0.5-1.5 microns, moved along linear structures (termed transport bundles) at an average maximum rate of 0.48 micron/sec. After particle movement had been viewed, cultures were preserved by rapid perfusion with chemical fixative. To determine whether particle transport occurred along actin, we then used fluorescence microscopy to correlate this movement with actin and microtubule distributions in the same growth cones. The observed transport bundles colocalized with actin but not with microtubules. The rates of particle movement and the association of moving particles with actin filament bundles suggest that myosins may participate in the transport of organelles (or other materials) in intact neurons.

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.

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.

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.