Class 3 Hox genes in insects and the origin of zen.

We have cloned, from a beetle and a locust, genes that are homologous to the class 3 Hox genes of vertebrates. Outside the homeobox they share sequence motifs with the Drosophila zerknüllt (zen) and z2 genes, and like zen, are expressed only in extraembryonic membranes. We conclude that the zen genes of Drosophila derive from a Hox class 3 sequence that formed part of the common ancestral Hox cluster, but that in insects this (Hox) gene has lost its role in patterning the anterio-posterior axis of the embryo, and acquired a new function. In the lineage leading to Drosophila, the zen genes have diverged particularly rapidly.

Exon shuffling and the origin of the mitochondrial targeting function in plant cytochrome c1 precursor.

Since most of the examples of "exon shuffling" are between vertebrate genes, the view is often expressed that exon shuffling is limited to the evolutionarily recent lineage of vertebrates. Although exon shuffling in plants has been inferred from the analysis of intron phases of plant genes [Long, M., Rosenberg, C. & Gilbert, W. (1995) Proc. Natl. Acad. Sci. USA 92, 12495-12499] and from the comparison of two functionally unknown sunflower genes [Domon, C. & Steinmetz, A. (1994) Mol. Gen. Genet. 244, 312-317], clear cases of exon shuffling in plant genes remain to be uncovered. Here, we report an example of exon shuffling in two important nucleus-encoded plant genes: cytosolic glyceraldehyde-3-phosphate dehydrogenase (cytosolic GAPDH or GapC) and cytochrome c1 precursor. The intron-exon structures of the shuffled region indicate that the shuffling event took place at the DNA sequence level. In this case, we can establish a donor-recipient relationship for the exon shuffling. Three amino terminal exons of GapC have been donated to cytochrome c1, where, in a new protein environment, they serve as a source of the mitochondrial targeting function. This finding throws light upon an old important but unsolved question in gene evolution: the origin of presequences or transit peptides that generally exist in nucleus-encoded organelle genes.

Dynamic friction and the origin of the complexity of earthquake sources.

We study a simple antiplane fault of finite length embedded in a homogeneous isotropic elastic solid to understand the origin of seismic source heterogeneity in the presence of nonlinear rate- and state-dependent friction. All the mechanical properties of the medium and friction are assumed homogeneous. Friction includes a characteristic length that is longer than the grid size so that our models have a well-defined continuum limit. Starting from a heterogeneous initial stress distribution, we apply a slowly increasing uniform stress load far from the fault and we simulate the seismicity for a few 1000 events. The style of seismicity produced by this model is determined by a control parameter associated with the degree of rate dependence of friction. For classical friction models with rate-independent friction, no complexity appears and seismicity is perfectly periodic. For weakly rate-dependent friction, large ruptures are still periodic, but small seismicity becomes increasingly nonstationary. When friction is highly rate-dependent, seismicity becomes nonperiodic and ruptures of all sizes occur inside the fault. Highly rate-dependent friction destabilizes the healing process producing premature healing of slip and partial stress drop. Partial stress drop produces large variations in the state of stress that in turn produce earthquakes of different sizes. Similar results have been found by other authors using the Burridge and Knopoff model. We conjecture that all models in which static stress drop is only a fraction of the dynamic stress drop produce stress heterogeneity.

Skeletal matrices, muci, and the origin of invertebrate calcification.

The sudden appearance of calcified skeletons among many different invertebrate taxa at the Precambrian-Cambrian transition may have required minor reorganization of preexisting secretory functions. In particular, features of the skeletal organic matrix responsible for regulating crystal growth by inhibition may be derived from mucous epithelial excretions. The latter would have prevented spontaneous calcium carbonate overcrusting of soft tissues exposed to the highly supersaturated Late Proterozoic ocean [Knoll, A. H., Fairchild, I. J. & Swett, K. (1993) Palaios 8, 512-525], a putative function for which we propose the term "anticalcification." We tested this hypothesis by comparing the serological properties of skeletal water-soluble matrices and mucous excretions of three invertebrates--the scleractinian coral Galaxea fascicularis and the bivalve molluscs Mytilus edulis and Mercenaria mercenaria. Crossreactivities recorded between muci and skeletal water-soluble matrices suggest that these different secretory products have a high degree of homology. Furthermore, freshly extracted muci of Mytilus were found to inhibit calcium carbonate precipitation in solution.

On the origin of the order Artiodactyla.

The first known members of the order Artiodactyla appeared suddenly throughout the Holarctic region at the beginning of the Eocene. They are characterized by distinctive cursorial skeletal specializations. Owing to their abrupt appearance and the lack of transitional forms, the origin of the order is problematic. Descent from a "condylarth," specifically the arctocyonid Chriacus, has been suggested based on dental resemblances, but until now postcranial anatomy seemed to preclude close relationship between Arctocyonidae and Artiodactyla. A middle Paleocene specimen of a small arctocyonid (?Chriacus) reported here is much more similar to the oldest artiodactyl, Diacodexis, in the derived condition of the hindlimb, reviving the possibility that Artiodactyla evolved from an arctocyonid.

Archaeal-eubacterial mergers in the origin of Eukarya: phylogenetic classification of life.

A symbiosis-based phylogeny leads to a consistent, useful classification system for all life. "Kingdoms" and "Domains" are replaced by biological names for the most inclusive taxa: Prokarya (bacteria) and Eukarya (symbiosis-derived nucleated organisms). The earliest Eukarya, anaerobic mastigotes, hypothetically originated from permanent whole-cell fusion between members of Archaea (e.g., Thermoplasma-like organisms) and of Eubacteria (e.g., Spirochaeta-like organisms). Molecular biology, life-history, and fossil record evidence support the reunification of bacteria as Prokarya while subdividing Eukarya into uniquely defined subtaxa: Protoctista, Animalia, Fungi, and Plantae.

A single-stranded DNA binding protein binds the origin of replication of the duplex kinetoplast DNA.

Replication of the kinetoplast DNA (kDNA) minicircle of trypanosomatids initiates at a conserved 12-nt sequence, 5'-GGGGTTGGTGTA-3', termed the universal minicircle sequence (UMS). A sequence-specific single-stranded DNA-binding protein from Crithidia fasciculata binds the heavy strand of the 12-mer UMS. Whereas this UMS-binding protein (UMSBP) does not bind a duplex UMS dodecamer, it binds the double-stranded kDNA minicircle as well as a duplex minicircle fragment containing the origin-associated UMS. Binding of the minicircle origin region by the single-stranded DNA binding protein suggested the local unwinding of the DNA double helix at this site. Modification of thymine residues at this site by KMnO4 revealed that the UMS resides within an unwound or otherwise sharply distorted DNA at the minicircle origin region. Computer analysis predicts the sequence-directed curving of the minicircle origin region. Electrophoresis of a minicircle fragment containing the origin region in polyacrylamide gels revealed a significantly lower electrophoretic mobility than expected from its length. The fragment anomalous electrophoretic mobility is displayed only in its native conformation and is dependent on temperature and gel porosity, indicating the local curving of the DNA double helix. We suggest that binding of UMSBP at the minicircle origin of replication is possible through local unwinding of the DNA double helix at the UMS site. It is hypothesized here that this local melting is initiated through the untwisting of unstacked dinucleotide sequences at the bent origin site.

Genetic absolute dating based on microsatellites and the origin of modern humans.

We introduce a new genetic distance for microsatellite loci, incorporating features of the stepwise mutation model, and test its performance on microsatellite polymorphisms in humans, chimpanzees, and gorillas. We find that it performs well in determining the relations among the primates, but less well than other distance measures (not based on the stepwise mutation model) in determining the relations among closely related human populations. However, the deepest split in the human phylogeny seems to be accurately reconstructed by the new distance and separates African and non-African populations. The new distance is independent of population size and therefore allows direct estimation of divergence times if the mutation rate is known. Based on 30 microsatellite polymorphisms and a recently reported average mutation rate of 5.6 x 10(-4) at 15 dinucleotide microsatellites, we estimate that the deepest split in the human phylogeny occurred about 156,000 years ago. Unlike most previous estimates, ours requires no external calibration of the rate of molecular evolution. We can use such calibrations, however, to test our estimate.

Genetics and the origin of European languages.

A new set of European genetic data has been analyzed to dissect independent patterns of geographic variation. The most important cause of European genetic variation has been confirmed to correspond to the migration of Neolithic farmers from the area of origin of agriculture in the Middle East. The next most important component of genetic variation is apparently associated with a north-south gradient possibly due to adaptation to cold climates but also to the differentiation of the Uralic and the Indo-European language-speaking people; however, the relevant correlations are not significantly different from zero after elimination of the spatial autocorrelation. The third component is highly correlated with the infiltration of the Yamna ("Kurgan") people, nomadic pastoralists who domesticated the horse and who have been claimed to have spread Indo-European languages to Europe; this association, which is statistically significant even when taking spatial autocorrelations into account, does not completely exclude the hypothesis of Indo-European as the language of Neolithic farmers. It is possible that both expansions were responsible for the spread of different subfamilies of Indo-European languages, but our genetic data cannot resolve their relative importance.

Molecular evidence for an African origin of the Hawaiian endemic Hesperomannia (Asteraceae)

Identification of the progenitors of plants endemic to oceanic islands often is complicated by extreme morphological divergence between island and continental taxa. This is especially true for the Hawaiian Islands, which are 3,900 km from any continental source. We examine the origin of Hesperomannia, a genus of three species endemic to Hawaii that always have been placed in the tribe Mutisieae of the sunflower family. Phylogenetic analyses of representatives from all tribes in this family using the chloroplast gene ndhF (where ndhF is the ND5 protein of chloroplast NADH dehydrogenase) indicate that Hesperomannia belongs to the tribe Vernonieae. Phylogenetic comparisons within the Vernonieae using sequences of both ndhF and the internal transcribed spacer regions of nuclear ribosomal DNA reveal that Hesperomannia is sister to African species of Vernonia. Long-distance dispersal northeastward from Africa to southeast Asia and across the many Pacific Ocean island chains is the most likely explanation for this unusual biogeographic connection. The 17- to 26-million-year divergence time between African Vernonia and Hesperomannia estimated by the DNA sequences predates the age of the eight existing Hawaiian Islands. These estimates are consistent with an hypothesis that the progenitor of Hesperomannia arrived at one of the low islands of the Hawaiian-Emperor chain between the late Oligocene and mid-Miocene when these islands were above sea level. Subsequent to its arrival the southeast Pacific island chains served as steppingstones for dispersal to the existing Hawaiian Islands.

Ancient origin of pathogen recognition specificity conferred by the tomato disease resistance gene Pto

We have investigated the origin of the Pto disease resistance (R) gene that was previously identified in the wild tomato species Lycopersicon pimpinellifolium and isolated by map-based cloning. Pto encodes a serine-threonine protein kinase that specifically recognizes strains of Pseudomonas syringae pv. tomato (Pst) that express the avirulence gene avrPto. We examined an accession of the distantly related wild species Lycopersicon hirsutum var. glabratum that exhibits avrPto-specific resistance to Pst. The Pst resistance of L. hirsutum was introgressed into a susceptible Lycopersicon esculentum background to create the near-isogenic line 96T133-3. Resistance to Pst(avrPto) in 96T133-3 was inherited as a single dominant locus and cosegregated with a restriction fragment length polymorphism detected by the Pto gene. This observation suggested that a member of the Pto gene family confers Pst(avrPto) resistance in this L. hirsutum line. Here we report the cloning and characterization of four members of the Pto family from 96T133-3. One gene (LhirPto) is 97% identical to Pto and encodes a catalytically active protein kinase that elicits a hypersensitive response when coexpressed with avrPto in leaves of Nicotiana benthamiana. In common with the Pto kinase, the LhirPto protein physically interacts with AvrPto and downstream members of the Pto signaling pathway. Our studies indicate that R genes of the protein kinase class may not evolve rapidly in response to pathogen pressure and rather that their ability to recognize specific Avr proteins can be highly conserved.

The chimeric eukaryote: Origin of the nucleus from the karyomastigont in amitochondriate protists

We present a testable model for the origin of the nucleus, the membrane-bounded organelle that defines eukaryotes. A chimeric cell evolved via symbiogenesis by syntrophic merger between an archaebacterium and a eubacterium. The archaebacterium, a thermoacidophil resembling extant Thermoplasma, generated hydrogen sulfide to protect the eubacterium, a heterotrophic swimmer comparable to Spirochaeta or Hollandina that oxidized sulfide to sulfur. Selection pressure for speed swimming and oxygen avoidance led to an ancient analogue of the extant cosmopolitan bacterial consortium “Thiodendron latens.” By eubacterial-archaebacterial genetic integration, the chimera, an amitochondriate heterotroph, evolved. This “earliest branching protist” that formed by permanent DNA recombination generated the nucleus as a component of the karyomastigont, an intracellular complex that assured genetic continuity of the former symbionts. The karyomastigont organellar system, common in extant amitochondriate protists as well as in presumed mitochondriate ancestors, minimally consists of a single nucleus, a single kinetosome and their protein connector. As predecessor of standard mitosis, the karyomastigont preceded free (unattached) nuclei. The nucleus evolved in karyomastigont ancestors by detachment at least five times (archamoebae, calonymphids, chlorophyte green algae, ciliates, foraminifera). This specific model of syntrophic chimeric fusion can be proved by sequence comparison of functional domains of motility proteins isolated from candidate taxa.

Measurement of the binding of tyrosyl phosphopeptides to SH2 domains: a reappraisal.

Src homology 2 (SH2) domain-mediated interactions with phosphotyrosine residues are critical in many intracellular signal transduction pathways. Attempts to understand the determinants of specificity and selectivity of these interactions have prompted many binding studies that have used several techniques. Some discrepancies, in both the absolute and relative values of the dissociation constants for particular interactions, are apparent. To establish the correct dissociation constants and to understand the origin of these differences, we have analyzed three previously determined interactions using the techniques of surface plasmon resonance and isothermal titration calorimetry. We find that the binding of SH2 domains to phosphopeptides is weaker than generally presumed. A phosphopeptide based on the hamster polyoma middle tumor antigen interacts with the SH2 domain from Src with an equilibrium dissociation constant (Kd) of 600 nM; a phosphopeptide based on one binding site from the platelet-derived growth factor receptor binds to the N-terminal SH2 domain of the 1-phosphatidylinositol 3-kinase p85 subunit with a Kd of 300 nM; and a phosphopeptide based on the C terminus of Lck binds to the SH2 domain of Lck with a Kd of 4 microM. In addition, we demonstrate that avidity effects that result from the dimerization of glutathione S-transferase fusion proteins with SH2 domains could be responsible for overestimates of affinities for these interactions previously studied by surface plasmon resonance.

Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

Desynchronization of cells on the developmental path triggers the formation of spiral waves of cAMP during Dictyostelium aggregation

Whereas it is relatively easy to account for the formation of concentric (target) waves of cAMP in the course of Dictyostelium discoideum aggregation after starvation, the origin of spiral waves remains obscure. We investigate a physiologically plausible mechanism for the spontaneous formation of spiral waves of cAMP in D. discoideum. The scenario relies on the developmental path associated with the continuous changes in the activity of enzymes such as adenylate cyclase and phosphodiesterase observed during the hours that follow starvation. These changes bring the cells successively from a nonexcitable state to an excitable state in which they relay suprathreshold cAMP pulses, and then to autonomous oscillations of cAMP, before the system returns to an excitable state. By analyzing a model for cAMP signaling based on receptor desensitization, we show that the desynchronization of cells on this developmental path triggers the formation of fully developed spirals of cAMP. Developmental paths that do not correspond to the sequence of dynamic transitions no relay-relay-oscillations-relay are less able or fail to give rise to the formation of spirals.