The evolution of sexual size dimorphism in the house finch. IV. Population divergence in ontogeny

Differences among taxa in sexual size dimorphism of adults can be produced by changes in distinct developmental processes and thus may reflect different evolutionary histories. Here we examine whether divergence in sexual dimorphism of adults between recently established Montana and Alabama populations of the house finch (Carpodacus mexicanus) can be attributed to population differences in growth of males and females. In both populations, males and females were similar at hatching, but as a result of sex-specific growth attained sexual size dimorphism by the time of independence. Timing and extent of growth varied between the sexes: Females maintained maximum rates of growth for a longer time than males, whereas males had higher initial growth rates and achieved maximum growth earlier and at smaller sizes than females. Ontogeny of sexual dimorphism differed between populations, but in each population, sexual dimorphism in growth parameters and sexual dimorphism at the time of nest leaving were similar to sexual dimorphism of adults. Variation in growth of females contributed more to population divergence than did growth of males. In each population, we found close correspondence between patterns of sexual dimorphism in growth and population divergence in morphology of adults: Traits that were the most sexually dimorphic in growth in each population contributed the most to population divergence in both sexes. We suggest that sex-specific expression of phenotypic and genetic variation throughout the ontogeny of house finches can result in different responses to selection between males and females of the same age, and thus produce fast population divergence in the sexual size dimorphism.     

Rapid evolution in the Nebria gregaria group (Coleoptera: Carabidae) and the paleogeography of the Queen Charlotte Islands

Morphological differentiation in the ground beetles of the Nebria gregaria group, found on the Queen Charlotte Islands, has been used as support for the glacial refugium proposed for the northwest coast of North America. Two members of this species group, N. charlottae and N. louiseae, are restricted to cobble beaches in this archipelago. A third, N. haida, is found only in alpine regions of the archipelago and the adjacent mainland. The remaining two species of the gregaria group, N. lituyae and N. gregaria, show highly restricted distributions in the mountains of the Alaska panhandle and on the beaches of the Aleutian Islands, respectively. To determine the relationships of the five species, we conducted phylogenetic analyses on nucleotide sequence data obtained from five regions of the mitochondrial DNA. In total, 1835 bp were analyzed. The results suggest that one species, N. lituyae, does not belong in the gregaria group, and that only seven mutations separated the two most divergent of the four remaining species. We also conducted random amplified polymorphic DNA fingerprinting analyses on genomic DNA extracted from the five species. Analyses of genetic diversity revealed a lack of molecular differentiation among the Queen Charlotte species, suggesting that these populations may be postglacial in origin and that together N. gregaria, N. charlottae, N. louiseae, and N. haida might represent local variations of a single species. These results are consistent with conclusions derived for the morphological and genetical differentiation among Gasterosteus populations in the archipelago.     

Protein structure determination using a database of interatomic distance probabilities

The accelerated pace of genomic sequencing has increased the demand for structural models of gene products. Improved quantitative methods are needed to study the many systems (e.g., macromolecular assemblies) for which data are scarce. Here, we describe a new molecular dynamics method for protein structure determination and molecular modeling. An energy function, or database potential, is derived from distributions of interatomic distances obtained from a database of known structures. X-ray crystal structures are refined by molecular dynamics with the new energy function replacing the Van der Waals potential. Compared to standard methods, this method improved the atomic positions, interatomic distances, and side-chain dihedral angles of structures randomized to mimic the early stages of refinement. The greatest enhancement in side-chain placement was observed for groups that are characteristically buried. More accurate calculated model phases will follow from improved interatomic distances. Details usually seen only in high-resolution refinements were improved, as is shown by an R-factor analysis. The improvements were greatest when refinements were carried out using X-ray data truncated at 3.5 A. The database potential should therefore be a valuable tool for determining X-ray structures, especially when only low-resolution data are available.     

Acetylcholinesterase of Schistosoma mansoni--functional correlates. Contributed in honor of Professor Hans Neurath's 90th birthday

Acetylcholinesterase (AChE) is an enzyme broadly distributed in many species, including parasites. It occurs in multiple molecular forms that differ in their quaternary structure and mode of anchoring to the cell surface. This review summarizes biochemical and immunological investigations carried out in our laboratories on AChE of the helmint, Schistosoma mansoni. AChE appears in S. mansoni in two principal molecular forms, both globular, with sedimentation coefficients of approximately 6.5 and 8 S. On the basis of their substrate specificity and sensitivity to inhibitors, both are "true" acetylcholinesterases. Approximately half of the AChE activity of S. mansoni is located on the outer surface of the parasite, attached to the tegumental membrane via a covalently attached glycosylphosphatidylinositol anchor. The remainder is located within the parasite, mainly associated with muscle tissue. Whereas the internal enzyme is most likely involved in termination of neurotransmission at cholinergic synapses, the role of the surface enzyme remains to be established; there are, however, indications that it is involved in signal transduction. The two forms of AChE differ in their heparin-binding properties, only the internal 8 S form of the AChE being retained on a heparin column. The two forms differ also in their immunological specificity, since they are selectively recognized by different monoclonal antibodies. Polyclonal antibodies raised against S. mansoni AChE purified by affinity chromatography are specific for the parasite AChE, reacting with both molecular forms, but do not recognize AChE from other species. They interact with the surface-localized enzyme on the intact organism, and produce almost total complement-dependent killing of the parasite. S. mansoni AChE is thus demonstrated to be a functional protein, involved in multifaceted activities, which can serve as a suitable candidate for diagnostic purposes, vaccine development, and drug design.     

Protein-protein recognition: an experimental and computational study of the R89K mutation in Raf and its effect on Ras binding

Binding of the protein Raf to the active form of Ras promotes activation of the MAP kinase signaling pathway, triggering cell growth and differentiation. Raf/Arg89 in the center of the binding interface plays an important role determining Ras-Raf binding affinity. We have investigated experimentally and computationally the Raf-R89K mutation, which abolishes signaling in vivo. The binding to [gamma-35S]GTP-Ras of a fusion protein between the Raf-binding domain (RBD) of Raf and GST was reduced at least 175-fold by the mutation, corresponding to a standard binding free energy decrease of at least 3.0 kcal/mol. To compute this free energy and obtain insights into the microscopic interactions favoring binding, we performed alchemical simulations of the RBD, both complexed to Ras and free in solution, in which residue 89 is gradually mutated from Arg into Lys. The simulations give a standard binding free energy decrease of 2.9+/-1.9 kcal/mol, in agreement with experiment. The use of numerous runs with three different force fields allows insights into the sources of uncertainty in the free energy and its components. The binding decreases partly because of a 7 kcal/mol higher cost to desolvate Lys upon binding, compared to Arg, due to better solvent interactions with the more concentrated Lys charge in the unbound state. This effect is expected to be general, contributing to the lower propensity of Lys to participate in protein-protein interfaces. Large contributions to the free energy change also arise from electrostatic interactions with groups up to 8 A away, namely residues 37-41 in the conserved effector domain of Ras (including 4 kcal/mol from Ser39 which loses a bifurcated hydrogen bond to Arg89), the conserved Lys84 and Lys87 of Raf, and 2-3 specific water molecules. This analysis will provide insights into the large experimental database of Ras-Raf mutations.     

Computational studies on mutant protein stability: The correlation between surface thermal expansion and protein stability

Thermal stability of mutant proteins has been investigated using temperature dependent molecular dynamics (MD) simulations in vacuo. The numerical modeling was aimed at mimicking protein expansion upon heating. After the conditions for an expanding protein accessible surface area were established for T4 lysozyme and barnase wild-type proteins, MD simulations were carried out under the same conditions using the crystal structures of several mutant proteins. The computed thermal expansion of the accessible surface area of mutant proteins was found to be strongly correlated with their experimentally measured stabilities. A similar, albeit weaker, correlation was observed for model mutant proteins. This opens the possibility of obtaining stability information directly from protein structure.     

A 4D TROSY-based pulse scheme for correlating 1HNi,15Ni,13Cαi,13C′i−1 chemical shifts in high molecular weight, 15N,13C, 2H labeled proteins

A 4D TROSY-based triple resonance experiment, 4D-HNCO_i–1CA_i, is presented which correlates intra-residue ¹HN, ¹⁵N, ¹³ C chemical shifts with the carbonyl (¹³C) shift of the preceding residue. The experiment is best used in concert with recently described 4D TROSY-HNCOCA and -HNCACO experiments [Yang, D. and Kay, L.E. (1999) J. Am. Chem. Soc., 121, 2571–2575]. In cases where degeneracy of (¹HN,¹⁵N) spin pairs precludes assignment using the HNCOCA and HNCACO, the HNCO_i–1CA_i often allows resolution of the ambiguity by linking the ¹³C and ¹³C spins surrounding the (¹HN,¹⁵N) pair. The experiment is demonstrated on a sample of ¹⁵N, ¹³C, ² H labeled maltose binding protein in complex with -cyclodextrin that tumbles with a correlation time of 46 ns.     

A rapid mini-prep method for isolation of ectomycorrhizal DNA from Tuber species

A rapid procedure has been developed to isolate DNA from the ectomycorrhizae of Tuber spp. for use in PCR experiments. The method described is fast and sensitive and can overcome the amplification problems that can arise in the presence of inhibitors. For this reason it can be used to type ectomycorrhizae even starting from a single root tip and make mycorrhizae identification much more rapid.     

RAPD divergence caused by microsite edaphic selection in wild barley

Random amplified polymorphic DNA polymerase chain reaction (RAPDPCR) was used to assess genetic diversity in four subpopulations (86 individuals) of wild barley, Hordeum spontaneum, sampled from Tabigha microsite near the Sea of Galilee, Israel. The microsite consists of two 100 m transects that are topographically separated by 100 m, each equally subdivided into 50 m of basalt and terra rossa soil types. Despite the same macroclimate characterizing the area around the Sea of Galilee, the microsite offers two edaphically different microhabitats, with basalt being a more ecologically heterogeneous and broader-niche than the relatively drier but more homogeneous and narrow-niche terra rossa. Analysis of 118 putative loci revealed significant (P<0.05) genetic differentiation in polymorphism (P0.05) between the two soils across the transects with P being higher in the more heterogeneous basalt (mean P0.05 = 0.902), than in terra rossa (mean P0.05 = 0.820). Gene diversity (He) was higher in basalt (mean He=0.371), than in terra rossa (mean He=0.259). Furthermore, unique alleles were confined to one soil type, either in one or both transects. Rare alleles were observed more frequently in terra rossa than basalt, and in transect II only. Gametic phase disequilibria showed a larger multilocus association of alleles in basalt than terra rossa, and in transect I than II. Spearman rank correlation (r_s) revealed a strong association between specific loci and soil types, and transects. Also, analysis of multilocus organization revealed soil-specific multilocus-genotypes. Therefore, our results suggest an edaphically differentiated genetic structure, which corroborates the niche width-variation hypothesis, and can be explained, in part, by natural selection. This pattern of RAPD diversity is in agreement with allozyme and hordein protein diversities in the same subpopulations studied previously. This revised version was published online in July 2006 with corrections to the Cover Date.     

A plasmid-based vector system for the cloning and expression of Helicobacter pylori genes encoding outer membrane proteins

Helicobacter pylori produces a number of proteins associated with the outer membrane, including adhesins and the vacuolating cytotoxin. We observed that the functional expression of such proteins is deleterious to Escherichia coli, the host bacterium used for gene cloning. Therefore, a general method was developed for the functional expression of such genes on a shuttle vector in H. pylori, which has been termed SOMPES (Shuttle vector-based Outer Membrane Protein Expression System). The intact, active gene is reconstituted by recombination in H. pylori from partial gene sequences cloned on an E. coli-H. pylori shuttle vector. This system was established in an H. pylori strain carrying a precise, unmarked chromosomal deletion of the vacA gene, which was constructed by adapting the streptomycin sensitivity system to H. pylori. It is based on the expression of the H. pylori rpsL gene as a counterselectable marker in the genetic background of an rpsL mutant. The utility of this approach is demonstrated by the expression of a recombinant gene encoding vacuolating cytotoxin (vacA) and a recombinant gene encoding an adherence-associated outer membrane protein (alpA) in H. pylori. Received: 10 May 1999 / Accepted: 7 July 1999