Variation in per capita interaction strength: Thresholds due to nonlinear dynamics and nonequilibrium conditions

I measured the strength of interaction between a marine herbivore and its growing resource over a realistic range of absolute and relative abundances. The herbivores (hermit crabs: Pagurus spp.) have slow and/or weak functional and numerical responses to epiphytic diatoms (Isthmia nervosa), which show logistic growth in the absence of consumers. By isolating this interaction in containers in the field, I mimicked many of the physical and biological variables characteristic of the intertidal while controlling the densities of focal species. The per capita effects of consumers on the population dynamics of their resource (i.e., interaction strength) were defined by using the relationship between hermit crab density and proportional change in the resource. When this relationship is fit by a Weibull function, a single parameter distinguishes constant interaction strength from one that varies as a function of density. Constant interaction strength causes the proportion of diatoms to fall linearly or proportionally as hermit crab density increases whereas per capita effects that increase with density cause an accelerating decline. Although many mathematical models of species interactions assume linear dynamics and invariant parameters, at least near equilibrium, the per capita effects of hermit crabs on diatoms varied substantially, apparently crossing a threshold from weak to strong when consumption exceeded resource production. This threshold separates a domain of coexistence from one of local extinction of the resource. Such thresholds may help explain trophic cascades, resource compensation, and context-dependent interaction strengths, while indicating a way to predict trophic effects, despite nonlinearities, as a function of vital rates.

Temporal variation, dormancy, and coexistence: A field test of the storage effect

Theoretical models suggest that overlapping generations, in combination with a temporally fluctuating environment, may allow the persistence of competitors that otherwise would not coexist. Despite extensive theoretical development, this “storage effect” hypothesis has received little empirical attention. Herein I present the first explicit mathematical analysis of the contribution of the storage effect to the dynamics of competing natural populations. In Oneida Lake, NY, data collected over the past 30 years show a striking negative correlation between the water-column densities of two species of suspension-feeding zooplankton, Daphnia galeata mendotae and Daphnia pulicaria. I have demonstrated competition between these two species and have shown that both possess long-lived eggs that establish overlapping generations. Moreover, recruitment to this long-lived stage varies annually, so that both daphnids have years in which they are favored (for recruitment) relative to their competitor. When the long-term population growth rates are calculated both with and without the effects of a variable environment, I show that D. galeata mendotae clearly cannot persist without the environmental variation and prolonged dormancy (i.e., storage effect) whereas D. pulicaria persists through consistently high per capita recruitment to the long-lived stage.

Phylogenetic assessment of length variation at a microsatellite locus

Sixty-six haplotypes at a locus containing a simple dinucleotide (CA)n microsatellite repeat were isolated by PCR–single-strand conformational polymorphism from populations of the horseshoe crab Limulus polyphemus. These haplotypes were sequenced to assess nucleotide variation directly. Thirty-four distinct sequences (alleles) were identified in a region 570 bp long that included the microsatellite motif. In the repeat region itself, CA-number varied in integer values from 5 to 11 across alleles, except that a (CA)8 class was not observed. Differences among alleles were due also to polymorphisms at 22 sites in regions immediately flanking the microsatellite repeats. Nucleotide substitutions in these regions were used to estimate phylogenetic relationships among alleles, and the gene phylogeny was used to trace the evolution of length variation and CA repeat numbers. A low correlation between size variation and genealogical relationships among alleles suggests that absolute fragment size (as normally scored in microsatellite assays) is an unreliable indicator of historical affinities among alleles. This finding on the molecular fine structure of microsatellite variation suggests the need for caution in the use of repeat counts at microsatellite loci as secure indicators of allelic relationships.

Evolution of plastid gene rps2 in a lineage of hemiparasitic and holoparasitic plants: Many losses of photosynthesis and complex patterns of rate variation

The plastid genomes of some nonphotosynthetic parasitic plants have experienced an extreme reduction in gene content and an increase in evolutionary rate of remaining genes. Nothing is known of the dynamics of these events or whether either is a direct outcome of the loss of photosynthesis. The parasitic Scrophulariaceae and Orobanchaceae, representing a continuum of heterotrophic ability ranging from photosynthetic hemiparasites to nonphotosynthetic holoparasites, are used to investigate these issues. We present a phylogenetic hypothesis for parasitic Scrophulariaceae and Orobanchaceae based on sequences of the plastid gene rps2, encoding the S2 subunit of the plastid ribosome. Parasitic Scrophulariaceae and Orobanchaceae form a monophyletic group in which parasitism can be inferred to have evolved once. Holoparasitism has evolved independently at least five times, with certain holoparasitic lineages representing single species, genera, and collections of nonphotosynthetic genera. Evolutionary loss of the photosynthetic gene rbcL is limited to a subset of holoparasitic lineages, with several holoparasites retaining a full length rbcL sequence. In contrast, the translational gene rps2 is retained in all plants investigated but has experienced rate accelerations in several hemi- as well as holoparasitic lineages, suggesting that there may be substantial molecular evolutionary changes to the plastid genome of parasites before the loss of photosynthesis. Independent patterns of synonymous and nonsynonymous rate acceleration in rps2 point to distinct mechanisms underlying rate variation in different lineages. Parasitic Scrophulariaceae (including the traditional Orobanchaceae) provide a rich platform for the investigation of molecular evolutionary process, gene function, and the evolution of parasitism.

Phase variation of the lpf operon is a mechanism to evade cross-immunity between Salmonella serotypes

Conventional wisdom holds that phase variation is a mechanism for immune evasion. However, despite fimbrial phase variation, mice previously exposed to Salmonella typhimurium are protected against a subsequent challenge. We evaluated whether lpf phase variation instead may be a mechanism to evade cross-immunity between Salmonella serotypes. Mice were immunized orally with S. typhimurium aroA mutants either that expressed the lpf operon (phase-on variant) or in which the entire lpf operon had been removed by deletion. During a subsequent challenge with virulent Salmonella enteritidis a selection against lpf phase-on variants was observed in mice previously exposed to S. typhimurium long polar fimbriae. Vaccination with S. typhimurium did not confer protection against challenge with S. enteritidis, presumably because lpf phase-off variants were able to evade cross-immunity. We propose that lpf phase variation is a mechanism to evade cross-immunity between Salmonella serotypes, thereby allowing their coexistence in a host population.

Allelic variation of a dehydrin gene cosegregates with chilling tolerance during seedling emergence

Dehydrins (DHNs, LEA D-11) are plant proteins present during environmental stresses associated with dehydration or low temperatures and during seed maturation. Functions of DHNs have not yet been defined. Earlier, we hypothesized that a ≈35-kDa DHN and membrane properties that reduce electrolyte leakage from seeds confer chilling tolerance during seedling emergence of cowpea (Vigna unguiculata L. Walp.) in an additive and independent manner. Evidence for this hypothesis was not rigorous because it was based on correlations of presence/absence of the DHN and slow electrolyte leakage with chilling tolerance in closely related cowpea lines that have some other genetic differences. Here, we provide more compelling genetic evidence for involvement of the DHN in chilling tolerance of cowpea. We developed near-isogenic lines by backcrossing. We isolated and determined the sequence of a cDNA corresponding to the ≈35-kDa DHN and used gene-specific oligonucleotides derived from it to test the genetic linkage between the DHN presence/absence trait and the DHN structural gene. We tested for association between the DHN presence/absence trait and both low-temperature seed emergence and electrolyte leakage. We show that allelic differences in the Dhn structural gene map to the same position as the DHN protein presence/absence trait and that the presence of the ≈35-kDa DHN is indeed associated with chilling tolerance during seedling emergence, independent of electrolyte leakage effects. Two types of allelic variation in the Dhn gene were identified in the protein-coding region, deletion of one Φ-segment from the DHN-negative lines and two single amino acid substitutions.

Global variation in the genetic and biochemical basis of diamondback moth resistance to Bacillus thuringiensis

Insecticidal proteins from the soil bacterium Bacillus thuringiensis (Bt) are becoming a cornerstone of ecologically sound pest management. However, if pests quickly adapt, the benefits of environmentally benign Bt toxins in sprays and genetically engineered crops will be short-lived. The diamondback moth (Plutella xylostella) is the first insect to evolve resistance to Bt in open-field populations. Here we report that populations from Hawaii and Pennsylvania share a genetic locus at which a recessive mutation associated with reduced toxin binding confers extremely high resistance to four Bt toxins. In contrast, resistance in a population from the Philippines shows multilocus control, a narrower spectrum, and for some Bt toxins, inheritance that is not recessive and not associated with reduced binding. The observed variation in the genetic and biochemical basis of resistance to Bt, which is unlike patterns documented for some synthetic insecticides, profoundly affects the choice of strategies for combating resistance.

Sequence variation in the Fanconi anemia gene FAA

Fanconi anemia (FA) is a genetically heterogeneous autosomal recessive syndrome associated with chromosomal instability, hypersensitivity to DNA crosslinking agents, and predisposition to malignancy. The gene for FA complementation group A (FAA) recently has been cloned. The cDNA is predicted to encode a polypeptide of 1,455 amino acids, with no homologies to any known protein that might suggest a function for FAA. We have used single-strand conformational polymorphism analysis to screen genomic DNA from a panel of 97 racially and ethnically diverse FA patients from the International Fanconi Anemia Registry for mutations in the FAA gene. A total of 85 variant bands were detected. Forty-five of the variants are probably benign polymorphisms, of which nine are common and can be used for various applications, including mapping studies for other genes in this region of chromosome 16q. Amplification refractory mutation system assays were developed to simplify their detection. Forty variants are likely to be pathogenic mutations. Seventeen of these are microdeletions/microinsertions associated with short direct repeats or homonucleotide tracts, a type of mutation thought to be generated by a mechanism of slipped-strand mispairing during DNA replication. A screening of 350 FA probands from the International Fanconi Anemia Registry for two of these deletions (1115–1118del and 3788–3790del) revealed that they are carried on about 2% and 5% of the FA alleles, respectively. 3788–3790del appears in a variety of ethnic groups and is found on at least two different haplotypes. We suggest that FAA is hypermutable, and that slipped-strand mispairing, a mutational mechanism recognized as important for the generation of germ-line and somatic mutations in a variety of cancer-related genes, including p53, APC, RB1, WT1, and BRCA1, may be a major mechanism for FAA mutagenesis.

Drought-induced shift of a forest–woodland ecotone: Rapid landscape response to climate variation

In coming decades, global climate changes are expected to produce large shifts in vegetation distributions at unprecedented rates. These shifts are expected to be most rapid and extreme at ecotones, the boundaries between ecosystems, particularly those in semiarid landscapes. However, current models do not adequately provide for such rapid effects—particularly those caused by mortality—largely because of the lack of data from field studies. Here we report the most rapid landscape-scale shift of a woody ecotone ever documented: in northern New Mexico in the 1950s, the ecotone between semiarid ponderosa pine forest and piñon–juniper woodland shifted extensively (2 km or more) and rapidly (<5 years) through mortality of ponderosa pines in response to a severe drought. This shift has persisted for 40 years. Forest patches within the shift zone became much more fragmented, and soil erosion greatly accelerated. The rapidity and the complex dynamics of the persistent shift point to the need to represent more accurately these dynamics, especially the mortality factor, in assessments of the effects of climate change.

A single response mechanism is responsible for evolutionary adaptive variation in a bird’s laying date

Many proximate factors determine a bird’s laying date, including environmental and social stimuli as well as individual responses to internal and external factors. However, the relative importance of these factors has not been experimentally demonstrated. Here we show that (i) large differences in the onset of first clutches between different populations result from variation in different responses to photoperiod and not from variation in responses to any other proximate factors and (ii) the same response mechanism causes maladaptive laying dates in habitats modified by humans. We present, to our knowledge, the first experimental demonstration that a single response mechanism is responsible for evolutionary adaptive intraspecific variation in a vertebrate life history trait.

Permanence or change? The meaning of genetic variation

Selected aspects of the evolutionary process and more specifically of the genetic variation are considered, with an emphasis in studies performed by my group. One key aspect of evolution seems to be the concomitant occurrence of dichotomic, contradictory (dialect) processes. Genetic variation is structured, and the dynamics of change at one level is not necessarily paralleled by that in another. The pathogenesis-related protein superfamily can be cited as an example in which permanence (the maintenance of certain key genetic features) coexists with change (modifications that led to different functions in different classes of organisms). Relationships between structure and function are exemplified by studies with hemoglobin Porto Alegre. The genetic structure of tribal populations may differ in important aspects from that of industrialized societies. Evolutionary histories also may differ when considered through the investigation of patrilineal or matrilineal lineages. Global evaluations taking into consideration all of these aspects are needed if we really want to understand the meaning of genetic variation.

Functional and physiological consequences of genetic variation at phosphoglucose isomerase: Heat shock protein expression is related to enzyme genotype in a montane beetle

Allele frequency variation at the phosphoglucose isomerase (PGI) locus in Californian populations of the beetle Chrysomela aeneicollis suggests that PGI may be undergoing natural selection. We quantified (i) apparent Michaelis-Menten constant (Km) of fructose 6-phosphate at different temperatures and (ii) thermal stability for three common PGI genotypes (1–1, 1–4, and 4–4). We also measured air temperature (Ta) and beetle body temperature (Tb) in three montane drainages in the Sierra Nevada, California. Finally, we measured 70-kDa heat shock protein (Hsp70) expression in field-collected and laboratory-acclimated beetles. We found that PGI allele 1 predominated in the northernmost drainage, Rock Creek (RC), which was also significantly cooler than the southernmost drainage, Big Pine Creek (BPC), where PGI allele 4 predominated. Allele frequencies and air temperatures were intermediate in the middle drainage, Bishop Creek (BC). Differences among genotypes in Km (1–1 > 1–4 > 4–4) and thermal stability (4–4 > 1–4 > 1–1) followed a pattern consistent with temperature adaptation. In nature, Tb was closely related to Ta. Hsp70 expression in adult beetles decreased with elevation and differed among drainages (BPC > BC > RC). After laboratory acclimation (8 days, 20°C day, 4°C night) and heat shock (4 h, 28–36°C), Hsp70 expression was greater for RC than BPC beetles. In RC, field-collected beetles homozygous for PGI 1–1 had higher Hsp70 levels than heterozygotes or a 4–4 homozygote. These results reveal functional and physiological differences among PGI genotypes, which suggest that montane populations of this beetle are locally adapted to temperature.

The nucleotide changes governing cuticular hydrocarbon variation and their evolution in Drosophila melanogaster

The cuticular hydrocarbon (CH) pheromones in Drosophila melanogaster exhibit strong geographic variation. African and Caribbean populations have a high ratio of 5,9 heptacosadiene/7,11 heptacosadiene (the “High” CH type), whereas populations from all other areas have a low ratio (“Low” CH type). Based on previous genetic mapping, DNA markers were developed that localized the genetic basis of this CH polymorphism to within a 13-kb region. We then carried out a hierarchical search for diagnostic nucleotide sites starting with four lines, and increasing to 24 and 43 lines from a worldwide collection. Within the 13-kb region, only one variable site shows a complete concordance with the CH phenotype. This is a 16-bp deletion in the 5′ region of a desaturase gene (desat2) that was recently suggested to be responsible for the CH polymorphism on the basis of its expression [Dallerac, R., Labeur, C., Jallon, J.-M., Knipple, D. C., Roelofs, W. L. & Wicker-Thomas, C. (2000) Proc. Natl. Acad. Sci. 97, 9449–9454]. The cosmopolitan Low type is derived from the ancestral High type, and DNA sequence variations suggest that the former spread worldwide with the aid of positive selection. Whether this CH variation could be a component of the sexual isolation between Zimbabwe and other cosmopolitan populations remains an interesting and unresolved question.

Global survey of genetic variation in CCR5, RANTES, and MIP-1α: Impact on the epidemiology of the HIV-1 pandemic

Expression of CC chemokine receptor 5 (CCR5), the major coreceptor for HIV-1 cell entry, and its ligands (e.g., RANTES and MIP-1α) is widely regarded as central to the pathogenesis of HIV-1 infection. By surveying nearly 3,000 HIV+ and HIV− individuals from worldwide populations for polymorphisms in the genes encoding RANTES, MIP-1α, and CCR5, we show that the evolutionary histories of human populations have had a significant impact on the distribution of variation in these genes, and that this may be responsible, in part, for the heterogeneous nature of the epidemiology of the HIV-1 pandemic. The varied distribution of RANTES haplotypes (AC, GC, and AG) associated with population-specific HIV-1 transmission- and disease-modifying effects is a striking example. Homozygosity for the AC haplotype was associated with an increased risk of acquiring HIV-1 as well as accelerated disease progression in European Americans, but not in African Americans. Yet, the prevalence of the ancestral AC haplotype is high in individuals of African origin, but substantially lower in non-Africans. In a Japanese cohort, AG-containing RANTES haplotype pairs were associated with a delay in disease progression; however, we now show that their contribution to HIV-1 pathogenesis and epidemiology in other parts of the world is negligible because the AG haplotype is infrequent in non-Far East Asians. Thus, the varied distribution of RANTES, MIP-1α, and CCR5 haplotype pairs and their population-specific phenotypic effects on HIV-1 susceptibility and disease progression results in a complex pattern of biological determinants of HIV-1 epidemiology. These findings have important implications for the design, assessment, and implementation of effective HIV-1 intervention and prevention strategies.

Flower color variation: A model for the experimental study of evolution

We review the study of flower color polymorphisms in the morning glory as a model for the analysis of adaptation. The pathway involved in the determination of flower color phenotype is traced from the molecular and genetic levels to the phenotypic level. Many of the genes that determine the enzymatic components of flavonoid biosynthesis are redundant, but, despite this complexity, it is possible to associate discrete floral phenotypes with individual genes. An important finding is that almost all of the mutations that determine phenotypic differences are the result of transposon insertions. Thus, the flower color diversity seized on by early human domesticators of this plant is a consequence of the rich variety of mobile elements that reside in the morning glory genome. We then consider a long history of research aimed at uncovering the ecological fate of these various flower phenotypes in the southeastern U.S. A large body of work has shown that insect pollinators discriminate against white phenotypes when white flowers are rare in populations. Because the plant is self-compatible, pollinator bias causes an increase in self-fertilization in white maternal plants, which should lead to an increase in the frequency of white genes, according to modifier gene theory. Studies of geographical distributions indicate other, as yet undiscovered, disadvantages associated with the white phenotype. The ultimate goal of connecting ecology to molecular genetics through the medium of phenotype is yet to be attained, but this approach may represent a model for analyzing the translation between these two levels of biological organization.