Horizontal gene transfer among genomes: The complexity hypothesis

Increasingly, studies of genes and genomes are indicating that considerable horizontal transfer has occurred between prokaryotes. Extensive horizontal transfer has occurred for operational genes (those involved in housekeeping), whereas informational genes (those involved in transcription, translation, and related processes) are seldomly horizontally transferred. Through phylogenetic analysis of six complete prokaryotic genomes and the identification of 312 sets of orthologous genes present in all six genomes, we tested two theories describing the temporal flow of horizontal transfer. We show that operational genes have been horizontally transferred continuously since the divergence of the prokaryotes, rather than having been exchanged in one, or a few, massive events that occurred early in the evolution of prokaryotes. In agreement with earlier studies, we found that differences in rates of evolution between operational and informational genes are minimal, suggesting that factors other than rate of evolution are responsible for the observed differences in horizontal transfer. We propose that a major factor in the more frequent horizontal transfer of operational genes is that informational genes are typically members of large, complex systems, whereas operational genes are not, thereby making horizontal transfer of informational gene products less probable (the complexity hypothesis).

Toward an understanding of the biochemical and pharmacological complexity of the saliva of a hematophagous sand fly Lutzomyia longipalpis

The saliva of blood-sucking arthropods contains powerful pharmacologically active substances and may be a vaccine target against some vector-borne diseases. Subtractive cloning combined with biochemical approaches was used to discover activities in the salivary glands of the hematophagous fly Lutzomyia longipalpis. Sequences of nine full-length cDNA clones were obtained, five of which are possibly associated with blood-meal acquisition, each having cDNA similarity to: (i) the bed bug Cimex lectularius apyrase, (ii) a 5′-nucleotidase/phosphodiesterase, (iii) a hyaluronidase, (iv) a protein containing a carbohydrate-recognition domain (CRD), and (v) a RGD-containing peptide with no significant matches to known proteins in the blast databases. Following these findings, we observed that the salivary apyrase activity of L. longipalpis is indeed similar to that of Cimex apyrase in its metal requirements. The predicted isoelectric point of the putative apyrase matches the value found for Lutzomyia salivary apyrase. A 5′-nucleotidase, as well as hyaluronidase activity, was found in the salivary glands, and the CRD-containing cDNA matches the N-terminal sequence of the HPLC-purified salivary anticlotting protein. A cDNA similar to α-amylase was discovered and salivary enzymatic activity demonstrated for the first time in a blood-sucking arthropod. Full-length clones were also found coding for three proteins of unknown function matching, respectively, the N-terminal sequence of an abundant salivary protein, having similarity to the CAP superfamily of proteins and the Drosophila yellow protein. Finally, two partial sequences are reported that match possible housekeeping genes. Subtractive cloning will considerably enhance efforts to unravel the salivary pharmacopeia of blood-sucking arthropods.

Fundamental patterns underlying gene expression profiles: Simplicity from complexity

Analysis of previously published sets of DNA microarray gene expression data by singular value decomposition has uncovered underlying patterns or “characteristic modes” in their temporal profiles. These patterns contribute unequally to the structure of the expression profiles. Moreover, the essential features of a given set of expression profiles are captured using just a small number of characteristic modes. This leads to the striking conclusion that the transcriptional response of a genome is orchestrated in a few fundamental patterns of gene expression change. These patterns are both simple and robust, dominating the alterations in expression of genes throughout the genome. Moreover, the characteristic modes of gene expression change in response to environmental perturbations are similar in such distant organisms as yeast and human cells. This analysis reveals simple regularities in the seemingly complex transcriptional transitions of diverse cells to new states, and these provide insights into the operation of the underlying genetic networks.

Genetic complexity of pathogen perception by plants: The example of Rcr3, a tomato gene required specifically by Cf-2

Genetic analysis of plant–pathogen interactions has demonstrated that resistance to infection is often determined by the interaction of dominant plant resistance (R) genes and dominant pathogen-encoded avirulence (Avr) genes. It was postulated that R genes encode receptors for Avr determinants. A large number of R genes and their cognate Avr genes have now been analyzed at the molecular level. R gene loci are extremely polymorphic, particularly in sequences encoding amino acids of the leucine-rich repeat motif. A major challenge is to determine how Avr perception by R proteins triggers the plant defense response. Mutational analysis has identified several genes required for the function of specific R proteins. Here we report the identification of Rcr3, a tomato gene required specifically for Cf-2-mediated resistance. We propose that Avr products interact with host proteins to promote disease, and that R proteins “guard” these host components and initiate Avr-dependent plant defense responses.

Molecular phylogeny analysis of fiddler crabs: test of the hypothesis of increasing behavioral complexity in evolution.

The current phylogenetic hypothesis for the evolution and biogeography of fiddler crabs relies on the assumption that complex behavioral traits are assumed to also be evolutionary derived. Indo-west Pacific fiddler crabs have simpler reproductive social behavior and are more marine and were thought to be ancestral to the more behaviorally complex and more terrestrial American species. It was also hypothesized that the evolution of more complex social and reproductive behavior was associated with the colonization of the higher intertidal zones. Our phylogenetic analysis, based upon a set of independent molecular characters, however, demonstrates how widely entrenched ideas about evolution and biogeography led to a reasonable, but apparently incorrect, conclusion about the evolutionary trends within this pantropical group of crustaceans. Species bearing the set of "derived traits" are phylogenetically ancestral, suggesting an alternative evolutionary scenario: the evolution of reproductive behavioral complexity in fiddler crabs may have arisen multiple times during their evolution. The evolution of behavioral complexity may have arisen by coopting of a series of other adaptations for high intertidal living and antipredator escape. A calibration of rates of molecular evolution from populations on either side of the Isthmus of Panama suggest a sequence divergence rate for 16S rRNA of 0.9% per million years. The divergence between the ancestral clade and derived forms is estimated to be approximately 22 million years ago, whereas the divergence between the American and Indo-west Pacific is estimated to be approximately 17 million years ago.

Slip complexity in earthquake fault models.

We summarize studies of earthquake fault models that give rise to slip complexities like those in natural earthquakes. For models of smooth faults between elastically deformable continua, it is critical that the friction laws involve a characteristic distance for slip weakening or evolution of surface state. That results in a finite nucleation size, or coherent slip patch size, h*. Models of smooth faults, using numerical cell size properly small compared to h*, show periodic response or complex and apparently chaotic histories of large events but have not been found to show small event complexity like the self-similar (power law) Gutenberg-Richter frequency-size statistics. This conclusion is supported in the present paper by fully inertial elastodynamic modeling of earthquake sequences. In contrast, some models of locally heterogeneous faults with quasi-independent fault segments, represented approximately by simulations with cell size larger than h* so that the model becomes "inherently discrete," do show small event complexity of the Gutenberg-Richter type. Models based on classical friction laws without a weakening length scale or for which the numerical procedure imposes an abrupt strength drop at the onset of slip have h* = 0 and hence always fall into the inherently discrete class. We suggest that the small-event complexity that some such models show will not survive regularization of the constitutive description, by inclusion of an appropriate length scale leading to a finite h*, and a corresponding reduction of numerical grid size.

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.

Slip complexity in dynamic models of earthquake faults.

We summarize recent evidence that models of earthquake faults with dynamically unstable friction laws but no externally imposed heterogeneities can exhibit slip complexity. Two models are described here. The first is a one-dimensional model with velocity-weakening stick-slip friction; the second is a two-dimensional elastodynamic model with slip-weakening friction. Both exhibit small-event complexity and chaotic sequences of large characteristic events. The large events in both models are composed of Heaton pulses. We argue that the key ingredients of these models are reasonably accurate representations of the properties of real faults.

A complexity measure for selective matching of signals by the brain.

We have previously derived a theoretical measure of neural complexity (CN) in an attempt to characterize functional connectivity in the brain. CN measures the amount and heterogeneity of statistical correlations within a neural system in terms of the mutual information between subsets of its units. CN was initially used to characterize the functional connectivity of a neural system isolated from the environment. In the present paper, we introduce a related statistical measure, matching complexity (CM), which reflects the change in CN that occurs after a neural system receives signals from the environment. CM measures how well the ensemble of intrinsic correlations within a neural system fits the statistical structure of the sensory input. We show that CM is low when the intrinsic connectivity of a simulated cortical area is randomly organized. Conversely, CM is high when the intrinsic connectivity is modified so as to differentially amplify those intrinsic correlations that happen to be enhanced by sensory input. When the input is represented by an individual stimulus, a positive value of CM indicates that the limited mutual information between sensory sheets sampling the stimulus and the rest of the brain triggers a large increase in the mutual information between many functionally specialized subsets within the brain. In this way, a complex brain can deal with context and go "beyond the information given."

Complexity of the erythroid transcription factor NF-E2 as revealed by gene targeting of the mouse p18 NF-E2 locus.

High-level globin expression in erythroid precursor cells depends on the integrity of NF-E2 recognition sites, transcription factor AP-1-like protein-binding motifs, located in the upstream regulatory regions of the alpha- and beta-globin loci. The NF-E2 transcription factor, which recognizes these sites, is a heterodimer consisting of (i) p45 NF-E2 (the larger subunit), a hematopoietic-restricted basic leucine zipper protein, and (ii) a widely expressed basic leucine zipper factor, p18 NF-E2, the smaller subunit. p18 NF-E2 protein shares extensive homology with the maf protooncogene family. To determine an in vivo role for p18 NF-E2 protein we disrupted the p18 NF-E2-encoding gene by homologous recombination in murine embryonic stem cells and generated p18 NF-E2-/- mice. These mice are indistinguishable from littermates throughout all phases of development and remain healthy in adulthood. Despite the absence of expressed p18 NF-E2, DNA-binding activity with the properties of the NF-E2 heterodimer is present in fetal liver erythroid cells of p18 NF-E2-/- mice. We speculate that another member of the maf basic leucine zipper family substitutes for the p18 subunit in a complex with p45 NF-E2. Thus, p18 NF-E2 per se appears to be dispensable in vivo.

Complexity, contingency, and criticality.

Complexity originates from the tendency of large dynamical systems to organize themselves into a critical state, with avalanches or "punctuations" of all sizes. In the critical state, events which would otherwise be uncoupled become correlated. The apparent, historical contingency in many sciences, including geology, biology, and economics, finds a natural interpretation as a self-organized critical phenomenon. These ideas are discussed in the context of simple mathematical models of sandpiles and biological evolution. Insights are gained not only from numerical simulations but also from rigorous mathematical analysis.