Tuning of activation thresholds explains flexibility in the selection and development of T cells in the thymus

Immature CD4+CD8+ thymocytes expressing T-cell antigen receptors (TCR) are selected by TCR-mediated recognition of peptides associated with major histocompatibility complex molecules on thymic stromal cells. Selection ensures reactivity of the mature cells to foreign antigens and tolerance to self. Although much has been learned about the factors that determine whether a thymocyte with a given specificity will be positively or negatively selected, selection as an aspect of the developmental process as a whole is less well-understood. Here we invoke a model in which thymocytes tune their response characteristics individually and dynamically in the course of development. Cellular development and selection are driven by receptor-mediated metabolic perturbations. Perturbation is a measure of the net intracellular change induced by external stimulation. It results from the integration of several signals and countersignals over time and therefore depends on the environment and the maturation stage of the cell. Individual cell adaptation limits the range of perturbations. Such adaptation renders thymocytes less sensitive to the level of stimulation per se, but responsive to environmental changes in that level. This formulation begins to explain the mechanisms that link developmental and selection events to each other.

Evolutionary theory predicts late-life mortality plateaus

Most demographic data indicate a roughly exponential increase in adult mortality with age, a phenomenon that has been explained in terms of a decline in the force of natural selection acting on age-specific mortality. Scattered demographic findings suggest the existence of a late-life mortality plateau in both humans and dipteran insects, seemingly at odds with both prior data and evolutionary theory. Extensions to the evolutionary theory of aging are developed which indicate that such late-life mortality plateaus are to be expected when enough late-life data are collected. This expanded theory predicts late-life mortality plateaus, with both antagonistic pleiotropy and mutation accumulation as driving population genetic mechanisms.

Evolutionary conservation of apoptosis mechanisms: Lepidopteran and baculoviral inhibitor of apoptosis proteins are inhibitors of mammalian caspase-9

We cloned a new inhibitor of apoptosis protein (IAP) homolog, SfIAP, from Spodoptera frugiperda Sf-21 cells, a host of insect baculoviruses. SfIAP contains two baculovirus IAP repeat domains followed by a RING domain. SfIAP has striking amino acid sequence similarity with baculoviral IAPs, CpIAP and OpIAP, suggesting that baculoviral IAPs may be host-derived genes. SfIAP and baculoviral CpIAP inhibit Bax but not Fas-induced apoptosis in human cells. Their apoptosis-suppressing activity in mammalian cells requires both baculovirus IAP repeat and RING domains. Further biochemical data suggest that SfIAP and CpIAP are specific inhibitors of mammalian caspase-9, the pinnacle caspase in the mitochondria/cytochrome c pathway for apoptosis, but are not inhibitors of downstream caspase-3 and caspase-7. Thus the mechanisms by which insect and baculoviral IAPs suppress apoptosis may involve inhibition of an insect caspase-9 homologue. Peptides representing the IAP-binding domain of the Drosophila cell death protein Grim abrogated human caspase suppression by SfIAP and CpIAP, implying evolutionary conservation of the functions of IAPs and their inhibitors.

Cryptocyanin, a crustacean molting protein: Evolutionary link with arthropod hemocyanins and insect hexamerins

Cryptocyanin, a copper-free hexameric protein in crab (Cancer magister) hemolymph, has been characterized and the amino acid sequence has been deduced from its cDNA. It is markedly similar in sequence, size, and structure to hemocyanin, the copper-containing oxygen-transport protein found in many arthropods. Cryptocyanin does not bind oxygen, however, and lacks three of the six highly conserved copper-binding histidine residues of hemocyanin. Cryptocyanin has no phenoloxidase activity, although a phenoloxidase is present in the hemolymph. The concentration of cryptocyanin in the hemolymph is closely coordinated with the molt cycle and reaches levels higher than hemocyanin during premolt. Cryptocyanin resembles insect hexamerins in the lack of copper, molt cycle patterns of biosynthesis, and potential contributions to the new exoskeleton. Phylogenetic analysis of sequence similarities between cryptocyanin and other members of the hemocyanin gene family shows that cryptocyanin is closely associated with crustacean hemocyanins and suggests that cryptocyanin arose as a result of a hemocyanin gene duplication. The presence of both hemocyanin and cryptocyanin in one animal provides an example of how insect hexamerins might have evolved from hemocyanin. Our results suggest that multiple members of the hemocyanin gene family—hemocyanin, cryptocyanin, phenoloxidase, and hexamerins—may participate in two vital functions of molting animals, oxygen binding and molting. Cryptocyanin may provide important molecular data to further investigate evolutionary relationships among all molting animals.

Evolutionary relationships among diverse bacteriophages and prophages: All the world’s a phage

We report DNA and predicted protein sequence similarities, implying homology, among genes of double-stranded DNA (dsDNA) bacteriophages and prophages spanning a broad phylogenetic range of host bacteria. The sequence matches reported here establish genetic connections, not always direct, among the lambdoid phages of Escherichia coli, phage φC31 of Streptomyces, phages of Mycobacterium, a previously unrecognized cryptic prophage, φflu, in the Haemophilus influenzae genome, and two small prophage-like elements, φRv1 and φRv2, in the genome of Mycobacterium tuberculosis. The results imply that these phage genes, and very possibly all of the dsDNA tailed phages, share common ancestry. We propose a model for the genetic structure and dynamics of the global phage population in which all dsDNA phage genomes are mosaics with access, by horizontal exchange, to a large common genetic pool but in which access to the gene pool is not uniform for all phage.

Molecular phylogenetic analysis of evolutionary trends in stonefly wing structure and locomotor behavior

Insects in the order Plecoptera (stoneflies) use a form of two-dimensional aerodynamic locomotion called surface skimming to move across water surfaces. Because their weight is supported by water, skimmers can achieve effective aerodynamic locomotion even with small wings and weak flight muscles. These mechanical features stimulated the hypothesis that surface skimming may have been an intermediate stage in the evolution of insect flight, which has perhaps been retained in certain modern stoneflies. Here we present a phylogeny of Plecoptera based on nucleotide sequence data from the small subunit rRNA (18S) gene. By mapping locomotor behavior and wing structural data onto the phylogeny, we distinguish between the competing hypotheses that skimming is a retained ancestral trait or, alternatively, a relatively recent loss of flight. Our results show that basal stoneflies are surface skimmers, and that various forms of surface skimming are distributed widely across the plecopteran phylogeny. Stonefly wings show evolutionary trends in the number of cross veins and the thickness of the cuticle of the longitudinal veins that are consistent with elaboration and diversification of flight-related traits. These data support the hypothesis that the first stoneflies were surface skimmers, and that wing structures important for aerial flight have become elaborated and more diverse during the radiation of modern stoneflies.

Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships

Pairwise sequence comparison methods have been assessed using proteins whose relationships are known reliably from their structures and functions, as described in the scop database [Murzin, A. G., Brenner, S. E., Hubbard, T. & Chothia C. (1995) J. Mol. Biol. 247, 536–540]. The evaluation tested the programs blast [Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. 215, 403–410], wu-blast2 [Altschul, S. F. & Gish, W. (1996) Methods Enzymol. 266, 460–480], fasta [Pearson, W. R. & Lipman, D. J. (1988) Proc. Natl. Acad. Sci. USA 85, 2444–2448], and ssearch [Smith, T. F. & Waterman, M. S. (1981) J. Mol. Biol. 147, 195–197] and their scoring schemes. The error rate of all algorithms is greatly reduced by using statistical scores to evaluate matches rather than percentage identity or raw scores. The E-value statistical scores of ssearch and fasta are reliable: the number of false positives found in our tests agrees well with the scores reported. However, the P-values reported by blast and wu-blast2 exaggerate significance by orders of magnitude. ssearch, fasta ktup = 1, and wu-blast2 perform best, and they are capable of detecting almost all relationships between proteins whose sequence identities are >30%. For more distantly related proteins, they do much less well; only one-half of the relationships between proteins with 20–30% identity are found. Because many homologs have low sequence similarity, most distant relationships cannot be detected by any pairwise comparison method; however, those which are identified may be used with confidence.

Proposed mechanism for stability of proteins to evolutionary mutations

It is shown that the sequence-ordering tendencies induced by design into different fast-folding, thermally stable native structures interfere. This interference results in a type of quasiorthogonality between optimal native structures, which divides sequence space into fast-folding, thermally stable families surrounded by slow-folding, low stability shells. A concrete example of this effect is provided by using a simple α carbon type model in which a complete correspondence is established between sequence and structure. It is speculated that gaps can occur in the space of protein-like sequences separating the sequence families and resulting in a mechanism for stability and diversity of protein sequence information.

Inositol 1,4,5-tris-phosphate activation of inositol tris-phosphate receptor Ca2+ channel by ligand tuning of Ca2+ inhibition

Inositol 1,4,5-tris-phosphate (IP3) binding to its receptors (IP3R) in the endoplasmic reticulum (ER) activates Ca2+ release from the ER lumen to the cytoplasm, generating complex cytoplasmic Ca2+ concentration signals including temporal oscillations and propagating waves. IP3-mediated Ca2+ release is also controlled by cytoplasmic Ca2+ concentration with both positive and negative feedback. Single-channel properties of the IP3R in its native ER membrane were investigated by patch clamp electrophysiology of isolated Xenopus oocyte nuclei to determine the dependencies of IP3R on cytoplasmic Ca2+ and IP3 concentrations under rigorously defined conditions. Instead of the expected narrow bell-shaped cytoplasmic free Ca2+ concentration ([Ca2+]i) response centered at ≈300 nM–1 μM, the open probability remained elevated (≈0.8) in the presence of saturating levels (10 μM) of IP3, even as [Ca2+]i was raised to high concentrations, displaying two distinct types of functional Ca2+ binding sites: activating sites with half-maximal activating [Ca2+]i (Kact) of 210 nM and Hill coefficient (Hact) ≈2; and inhibitory sites with half-maximal inhibitory [Ca2+]i (Kinh) of 54 μM and Hill coefficient (Hinh) ≈4. Lowering IP3 concentration was without effect on Ca2+ activation parameters or Hinh, but decreased Kinh with a functional half-maximal activating IP3 concentration (KIP3) of 50 nM and Hill coefficient (HIP3) of 4 for IP3. These results demonstrate that Ca2+ is a true receptor agonist, whereas the sole function of IP3 is to relieve Ca2+ inhibition of IP3R. Allosteric tuning of Ca2+ inhibition by IP3 enables the individual IP3R Ca2+ channel to respond in a graded fashion, which has implications for localized and global cytoplasmic Ca2+ concentration signaling and quantal Ca2+ release.

Regulation of number and size of digits by posterior Hox genes: A dose-dependent mechanism with potential evolutionary implications

The proper development of digits, in tetrapods, requires the activity of several genes of the HoxA and HoxD homeobox gene complexes. By using a variety of loss-of-function alleles involving the five Hox genes that have been described to affect digit patterning, we report here that the group 11, 12, and 13 genes control both the size and number of murine digits in a dose-dependent fashion, rather than through a Hox code involving differential qualitative functions. A similar dose–response is observed in the morphogenesis of the penian bone, the baculum, which further suggests that digits and external genitalia share this genetic control mechanism. A progressive reduction in the dose of Hox gene products led first to ectrodactyly, then to olygodactyly and adactyly. Interestingly, this transition between the pentadactyl to the adactyl formula went through a step of polydactyly. We propose that in the distal appendage of polydactylous short-digited ancestral tetrapods, such as Acanthostega, the HoxA complex was predominantly active. Subsequent recruitment of the HoxD complex contributed to both reductions in digit number and increase in digit length. Thus, transition through a polydactylous limb before reaching and stabilizing the pentadactyl pattern may have relied, at least in part, on asynchronous and independent changes in the regulation of HoxA and HoxD gene complexes.

Evolutionary specialization of the nuclear targeting apparatus

The α- and β-karyopherins (Kaps), also called importins, mediate the nuclear transport of proteins. All α-Kaps contain a central domain composed of eight approximately 40 amino acid, tandemly arranged, armadillo-like (Arm) repeats. The number and order of these repeats have not changed since the common origin of fungi, plants, and mammals. Phylogenetic analysis suggests that the various α-Kaps fall into two groups, α1 and α2. Whereas animals encode both types, the yeast genome encodes only an α1-Kap. The β-Kaps are characterized by 14–15 tandemly arranged HEAT motifs. We show that the Arm repeats of α-Kaps and the HEAT motifs of β-Kaps are similar, suggesting that the α-Kaps and β-Kaps (and for that matter, all Arm and HEAT repeat-containing proteins) are members of the same protein superfamily. Phylogenetic analysis indicates that there are at least three major groups of β-Kaps, consistent with their proposed cargo specificities. We present a model in which an α-independent β-Kap progenitor gave rise to the α-dependent β-Kaps and the α-Kaps.

A fully automatic evolutionary classification of protein folds: Dali Domain Dictionary version 3

The Dali Domain Dictionary (http://www.ebi.ac.uk/dali/domain) is a numerical taxonomy of all known structures in the Protein Data Bank (PDB). The taxonomy is derived fully automatically from measurements of structural, functional and sequence similarities. Here, we report the extension of the classification to match the traditional four hierarchical levels corresponding to: (i) supersecondary structural motifs (attractors in fold space), (ii) the topology of globular domains (fold types), (iii) remote homologues (functional families) and (iv) homologues with sequence identity above 25% (sequence families). The computational definitions of attractors and functional families are new. In September 2000, the Dali classification contained 10 531 PDB entries comprising 17 101 chains, which were partitioned into five attractor regions, 1375 fold types, 2582 functional families and 3724 domain sequence families. Sequence families were further associated with 99 582 unique homologous sequences in the HSSP database, which increases the number of effectively known structures several-fold. The resulting database contains the description of protein domain architecture, the definition of structural neighbours around each known structure, the definition of structurally conserved cores and a comprehensive library of explicit multiple alignments of distantly related protein families.

Feeding specialization and host-derived chemical defense in Chrysomeline leaf beetles did not lead to an evolutionary dead end

Combination of molecular phylogenetic analyses of Chrysomelina beetles and chemical data of their defensive secretions indicate that two lineages independently developed, from an ancestral autogenous metabolism, an energetically efficient strategy that made the insect tightly dependent on the chemistry of the host plant. However, a lineage (the interrupta group) escaped this subordination through the development of a yet more derived mixed metabolism potentially compatible with a large number of new host-plant associations. Hence, these analyses on leaf beetles document a mechanism that can explain why high levels of specialization do not necessarily lead to “evolutionary dead ends.”

Lessons from the past: Evolutionary impacts of mass extinctions

Mass extinctions have played many evolutionary roles, involving differential survivorship or selectivity of taxa and traits, the disruption or preservation of evolutionary trends and ecosystem organization, and the promotion of taxonomic and morphological diversifications—often along unexpected trajectories—after the destruction or marginalization of once-dominant clades. The fossil record suggests that survivorship during mass extinctions is not strictly random, but it often fails to coincide with factors promoting survival during times of low extinction intensity. Although of very serious concern, present-day extinctions have not yet achieved the intensities seen in the Big Five mass extinctions of the geologic past, which each removed ≥50% of the subset of relatively abundant marine invertebrate genera. The best comparisons for predictive purposes therefore will involve factors such as differential extinction intensities among regions, clades, and functional groups, rules governing postextinction biotic interchanges and evolutionary dynamics, and analyses of the factors that cause taxa and evolutionary trends to continue unabated, to suffer setbacks but resume along the same trajectory, to survive only to fall into a marginal role or disappear (“dead clade walking”), or to undergo a burst of diversification. These issues need to be addressed in a spatially explicit framework, because the fossil record suggests regional differences in postextinction diversification dynamics and biotic interchanges. Postextinction diversifications lag far behind the initial taxonomic and morphological impoverishment and homogenization; they do not simply reoccupy vacated adaptive peaks, but explore opportunities as opened and constrained by intrinsic biotic factors and the ecological and evolutionary context of the radiation.

Disrupting evolutionary processes: The effect of habitat fragmentation on collared lizards in the Missouri Ozarks

Humans affect biodiversity at the genetic, species, community, and ecosystem levels. This impact on genetic diversity is critical, because genetic diversity is the raw material of evolutionary change, including adaptation and speciation. Two forces affecting genetic variation are genetic drift (which decreases genetic variation within but increases genetic differentiation among local populations) and gene flow (which increases variation within but decreases differentiation among local populations). Humans activities often augment drift and diminish gene flow for many species, which reduces genetic variation in local populations and prevents the spread of adaptive complexes outside their population of origin, thereby disrupting adaptive processes both locally and globally within a species. These impacts are illustrated with collared lizards (Crotaphytus collaris) in the Missouri Ozarks. Forest fire suppression has reduced habitat and disrupted gene flow in this lizard, thereby altering the balance toward drift and away from gene flow. This balance can be restored by managed landscape burns. Some have argued that, although human-induced fragmentation disrupts adaptation, it will also ultimately produce new species through founder effects. However, population genetic theory and experiments predict that most fragmentation events caused by human activities will facilitate not speciation, but local extinction. Founder events have played an important role in the macroevolution of certain groups, but only when ecological opportunities are expanding rather than contracting. The general impact of human activities on genetic diversity disrupts or diminishes the capacity for adaptation, speciation, and macroevolutionary change. This impact will ultimately diminish biodiversity at all levels.