Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach

A maximum likelihood estimator based on the coalescent for unequal migration rates and different subpopulation sizes is developed. The method uses a Markov chain Monte Carlo approach to investigate possible genealogies with branch lengths and with migration events. Properties of the new method are shown by using simulated data from a four-population n-island model and a source–sink population model. Our estimation method as coded in migrate is tested against genetree; both programs deliver a very similar likelihood surface. The algorithm converges to the estimates fairly quickly, even when the Markov chain is started from unfavorable parameters. The method was used to estimate gene flow in the Nile valley by using mtDNA data from three human populations.

Ligand occupancy of the alpha-V-beta3 integrin is necessary for smooth muscle cells to migrate in response to insulin-like growth factor.

Smooth muscle cells (SMCs) have been shown to migrate in response to insulin-like growth factor I (IGF-I). However, the mechanism mediating this response has not been determined. The migration rates of porcine and human vascular SMCs were assessed in a monolayer wounding assay. IGF-I and IGF-II induced increases of 141% and 97%, respectively, in the number of cells that migrated in 4 days. The presence of 0.2% fetal bovine serum in the culture medium was necessary for the IGFs to stimulate migration over uncoated plastic surfaces. However, if vitronectin was used as the substratum, IGF-I stimulated migration by 162% even in the absence of serum. To determine the role of integrins in mediating this migration, SMC surface proteins were labeled with 125I and immunoprecipitated with specific anti-integrin antibodies. Integrins containing alpha-V (vitronectin receptor), alpha5 (fibronectin receptor), and alpha3 (collagen/laminin receptor) subunits were the most abundant. IGF-I treatment caused a 73% reduction in alpha5-integrin subunit protein and a 25% increase in alpha-V subunit. More importantly, ligand binding of alpha-V-beta3 was increased by 2.4-fold. We therefore examined whether the function of the alpha-V-beta3 integrin was important for IGF-I-mediated migration. The disintegrin kistrin was shown by affinity crosslinking to specifically bind with high affinity to alpha-V-beta3 and not to alpha5-beta1 or other abundant integrins. The related disintegrin echistatin specifically inhibited 125I-labeled kistrin binding to alpha-V-beta3, while a structurally distinct disintegrin, decorsin, had 1000-fold lower affinity. The addition of increasing concentrations of either kistrin or echistatin inhibited IGF-I-induced migration, whereas decorsin had a minimal effect. The potency of these disintegrins in inhibiting IGF-I-induced migration paralleled their apparent affinity for the alpha-V integrin. Furthermore, an alpha-V-beta3 blocking antibody inhibited SMC migration by 80%. In summary, vitronectin receptor activation is a necessary component of IGF-I-mediated stimulation of smooth muscle migration, and alpha-V-beta3 integrin antagonists appear to be important reagents for modulating this process.

Biogeographic range expansion into South America by Coccidioides immitis mirrors New World patterns of human migration

Long-distance population dispersal leaves its characteristic signature in genomes, namely, reduced diversity and increased linkage between genetic markers. This signature enables historical patterns of range expansion to be traced. Herein, we use microsatellite loci from the human pathogen Coccidioides immitis to show that genetic diversity in this fungus is geographically partitioned throughout North America. In contrast, analyses of South American C. immitis show that this population is genetically depauperate and was founded from a single North American population centered in Texas. Variances of allele distributions show that South American C. immitis have undergone rapid population growth, consistent with an epidemic increase in postcolonization population size. Herein, we estimate the introduction into South America to have occurred within the last 9,000–140,000 years. This range increase parallels that of Homo sapiens. Because of known associations between Amerindians and this fungus, we suggest that the colonization of South America by C. immitis represents a relatively recent and rapid codispersal of a host and its pathogen.

Neural codes: Firing rates and beyond

Computational neuroscience has contributed significantly to our understanding of higher brain function by combining experimental neurobiology, psychophysics, modeling, and mathematical analysis. This article reviews recent advances in a key area: neural coding and information processing. It is shown that synapses are capable of supporting computations based on highly structured temporal codes. Such codes could provide a substrate for unambiguous representations of complex stimuli and be used to solve difficult cognitive tasks, such as the binding problem. Unsupervised learning rules could generate the circuitry required for precise temporal codes. Together, these results indicate that neural systems perform a rich repertoire of computations based on action potential timing.

Dynamic evolution of plant mitochondrial genomes: Mobile genes and introns and highly variable mutation rates

We summarize our recent studies showing that angiosperm mitochondrial (mt) genomes have experienced remarkably high rates of gene loss and concomitant transfer to the nucleus and of intron acquisition by horizontal transfer. Moreover, we find substantial lineage-specific variation in rates of these structural mutations and also point mutations. These findings mostly arise from a Southern blot survey of gene and intron distribution in 281 diverse angiosperms. These blots reveal numerous losses of mt ribosomal protein genes but, with one exception, only rare loss of respiratory genes. Some lineages of angiosperms have kept all of their mt ribosomal protein genes whereas others have lost most of them. These many losses appear to reflect remarkably high (and variable) rates of functional transfer of mt ribosomal protein genes to the nucleus in angiosperms. The recent transfer of cox2 to the nucleus in legumes provides both an example of interorganellar gene transfer in action and a starting point for discussion of the roles of mechanistic and selective forces in determining the distribution of genetic labor between organellar and nuclear genomes. Plant mt genomes also acquire sequences by horizontal transfer. A striking example of this is a homing group I intron in the mt cox1 gene. This extraordinarily invasive mobile element has probably been acquired over 1,000 times separately during angiosperm evolution via a recent wave of cross-species horizontal transfers. Finally, whereas all previously examined angiosperm mtDNAs have low rates of synonymous substitutions, mtDNAs of two distantly related angiosperms have highly accelerated substitution rates.

Crystal structure at 2.6-A resolution of human macrophage migration inhibitory factor.

Macrophage migration inhibitory factor (MIF) was the first cytokine to be described, but for 30 years its role in the immune response remained enigmatic. In recent studies, MIF has been found to be a novel pituitary hormone and the first protein identified to be released from immune cells on glucocorticoid stimulation. Once secreted, MIF counterregulates the immunosuppressive effects of steroids and thus acts as a critical component of the immune system to control both local and systemic immune responses. We report herein the x-ray crystal structure of human MIF to 2.6 angstrom resolution. The protein is a trimer of identical subunits. Each monomer contains two antiparallel alpha-helices that pack against a four-stranded beta-sheet. The monomer has an additional two beta-strands that interact with the beta-sheets of adjacent subunits to form the interface between monomers. The three beta-sheets are arranged to form a barrel containing a solvent-accessible channel that runs through the center of the protein along a molecular 3-fold axis. Electrostatic potential maps reveal that the channel has a positive potential, suggesting that it binds negatively charged molecules. The elucidated structure for MIF is unique among cytokines or hormonal mediators, and suggests that this counterregulator of glucocorticoid action participates in novel ligand-receptor interactions.

Positive selection and rates of evolution in immunodeficiency viruses from humans and chimpanzees.

Evolutionary theory predicts the recent spread of primate immunodeficiency viruses (PIVs) to new human populations to be accompanied by positive selection in response to new host environments and/or by random genetic drift. I assess evidence for positive selection in human and chimpanzee PIVs type I (PIV1s), using ratios of synonymous to nonsynonymous nucleotide change based on branch lengths and outgroup rooting. Ratios are smaller for PIV1s from humans than for PIV1 from a chimpanzee for the pol, gag, and env glycoprotein 120 (gp120) regions, indicating greater effects of positive selection in PIV1s from humans. Parsimony-based relative rate tests for amino acid changes showed significant differences between PIV1s from humans and chimpanzees in 18 of 48 pairwise comparisons, with all 18 showing faster rates of change in PIV1s from humans. This study indicates that in some instances, the recent evolution of human PIV1s follows a speciational pattern, in which increased diversification of taxa is correlated with greater amounts of character change appearing and being maintained through time. This extends the generality of the speciational pattern to a group of organisms (viruses) having the fastest known rates of anagenetic change for nucleotide characters and indicates that comprehensive understanding of PIV1 evolution requires consideration of both anagenetic change within viral lineages and the relative historical success of different viral clades. Phylogenetic analyses show that neither PIV1s infecting humans nor those infecting chimpanzees represent monophyletic groups and suggest multiple host-species shifts for PIV1s.

Optimization of rates of protein folding: the nucleation-condensation mechanism and its implications.

Small, single-module proteins that fold in a single cooperative step may be paradigms for understanding early events in protein-folding pathways generally. Recent experimental studies of the 64-residue chymotrypsin inhibitor 2 (CI2) support a nucleation mechanism for folding, as do some computer stimulations. CI2 has a nucleation site that develops only in the transition state for folding. The nucleus is composed of a set of adjacent residues (an alpha-helix), stabilized by long-range interactions that are formed as the rest of the protein collapses around it. A simple analysis of the optimization of the rate of protein folding predicts that rates are highest when the denatured state has little residual structure under physiological conditions and no intermediates accumulate. This implies that any potential nucleation site that is composed mainly of adjacent residues should be just weakly populated in the denatured state and become structured only in a high-energy intermediate or transition state when it is stabilized by interactions elsewhere in the protein. Hierarchical mechanisms of folding in which stable elements of structure accrete are unfavorable. The nucleation-condensation mechanism of CI2 fulfills the criteria for fast folding. On the other hand, stable intermediates do form in the folding of more complex proteins, and this may be an unavoidable consequence of increasing size and nucleation at more than one site.