Evolutionary stability of MHC class II haplotypes in diverse rhesus macaque populations

A thoroughly characterized breeding colony of 172 pedigreed rhesus macaques was used to analyze exon 2 of the polymorphic Mamu-DPB1, -DQA1, -DQB1, and -DRB loci. Most of the monkeys or their ancestors originated in India, though the panel also included animals from Burma and China, as well as some of unknown origin and mixed breeds. In these animals, mtDNA appears to correlate with the aforementioned geographic origin, and a large number of Mamu class II alleles were observed. The different Mamu-DPB1 alleles were largely shared between monkeys of different origin, whereas in humans particular alleles appear to be unique for ethnic populations. In contrast to Mamu-DPB1, the highly polymorphic -DQA1/DQB1 alleles form tightly linked pairs that appear to be about two-thirds population specific. For most of the DQA1/DQB1 pairs, Mamu-DRB region configurations present on the same chromosome have been ascertained, resulting in 41 different -DQ/DRB haplotypes. These distinct DQ/DRB haplotypes seem to be specific for monkeys of a determined origin. Thus, in evolutionary terms, the Mamu-DP, -DQ, and -DR regions show increasing instability with regard to allelic polymorphism, such as for -DP/DQ, or gene content and allelic polymorphism, such as for -DR, resulting in population-specific class II haplotypes. Furthermore, novel haplotypes are generated by recombination-like events. The results imply that mtDNA analysis in combination with Mhc typing is a helpful tool for selecting animals for biomedical experiments.The sequences reported in this paper have been deposited in the EMBL database (accession nos. AJ534296–AJ534304, AJ 564564, and AJ557455–AJ557511)     

Evolutionary determinants of population differences in population growth rate × habitat temperature interactions in Chironomus riparius

Little is known about intraspecific variation in fitness performance in response to thermal stress among natural populations and how this relates to evolutionary aspects of species ecology. In this study, population growth rate (PGR; a composite fitness measure) varied among five natural Chironomus riparius populations sampled across a climatic gradient when subjected to three temperature treatments reflecting the typical range of summer habitat temperatures (20, 24 and 28 °C). The variation could be explained by a complex model including effects of genetic drift, genetic diversity and adaptation to average temperature during the warmest month, in addition to experimental temperature. All populations suffered a decrease in PGR from 20 to 28 °C and ΔPGR was significantly correlated with the respective average habitat temperature in the warmest month—populations from warmer areas showing lower ΔPGR. This implies that long-term exposure to higher temperatures in the warmest month (the key reproductive period for C. riparius) is likely to be a key selective force influencing fitness at higher temperatures. A comparison of phenotypic divergence and neutral genetic differentiation revealed that one phenotypic trait—the number of fertile egg masses per female—appeared to be under positive selection in some populations. Our findings support a role for response to temperature selection along a climatic gradient and suggest population history is a key determinant of intraspecific fitness variation. We stress the importance of integrating different types of data (climatic, experimental, genetic) in order to understand the effects of global climate change on biodiversity.     

Evolutionary Behaviour,Trade-Offs and Cyclic and Chaotic Population Dynamics

Many studies of the evolution of life-history traits assume that the underlying population dynamical attractor is stable point equilibrium. However, evolutionary outcomes can change significantly in different circumstances. We present an analysis based on adaptive dynamics of a discrete-time demographic model involving a trade-off whose shape is also an important determinant of evolutionary behaviour. We derive an explicit expression for the fitness in the cyclic region and consequently present an adaptive dynamic analysis which is algebraic. We do this fully in the region of 2-cycles and (using a symbolic package) almost fully for 4-cycles. Simulations illustrate and verify our results. With equilibrium population dynamics, trade-offs with accelerating costs produce a continuously stable strategy (CSS) whereas trade-offs with decelerating costs produce a non-ES repellor. The transition to 2-cycles produces a discontinuous change: the appearance of an intermediate region in which branching points occur. The size of this region decreases as we move through the region of 2-cycles. There is a further discontinuous fall in the size of the branching region during the transition to 4-cycles. We extend our results numerically and with simulations to higher-period cycles and chaos. Simulations show that chaotic population dynamics can evolve from equilibrium and vice-versa.     

Evolutionary modifications of molecular structure of ATP-synthase γ-subunit

The ATP-synthase γ-subunit (FoF1) belongs to the rotor part of this oligomeric complex. Catalytic hydrolysis of adenosine triphosphate (ATP) is accompanied by rotation of γ-polypeptide inside the sphere formed by six subunits (αβ)₃ of the enzyme. The γ-subunit regulates ATPase and ATP-synthase activities of the FoF1. In the present work, evolutionary and reverse changes of this regulatory polypeptide and their effect on properties of the enzyme are studied. It is suggested that elongation of the γ-subunit globular part had resulted from the atpC intragene duplication in the process of adaptive evolution. The evolved fragment participates in light regulation of the chloroplast ATP-synthase.     

Evolutionary Psychology��s Place in Evolutionary Studies: a Tale of Promise and Challenge

The Evolutionary Studies (EvoS) Consortium and the academic programs born from its creation have been wildly successful in their initial ventures. These achievements are marked by feedback from across the EvoS campuses, the resultant scholarly work produced by participating students, and faculty collaborations spurred by exposure to the organization. The success of EvoS is probably best marked by the recent National Science Foundation grant (CCLI Award #0817337), awarded jointly to SUNY New Paltz and Binghamton University, with the purpose of expanding EvoS beyond the bounds of these two institutions. A particularly noteworthy element of many EvoS programs is the role of Evolutionary Psychology (EP), a perspective in the behavioral sciences that addresses questions of human behavior from the perspective of evolution. In light of several forms of data, including analyses of a variety of disciplines drawn on from evolutionary psychologists in their work, we argue that evolutionary psychologists may well be the most naturally interdisciplinary scholars within the behavioral sciences, making them highly appropriate for inclusion in EvoS. But our research shows not only promise regarding the relationship between EP and EvoS—challenges are raised as well. We present additional data showing that EP is currently represented disproportionately within the EvoS world—a fact that clearly shows that there are currently limitations to the potential impact of EvoS in modern academia. Scholars from other disciplines, particularly within the humanities and social sciences, seem to be missing the evolution revolution. Implications regarding how EvoS can broaden its scope to be even more powerful in its integrative scope are discussed.     

Evolutionary changes ofα-crystallin and the phylogeny of mammalian orders

Summary The sequences of the A chains of the eye lens protein-crystallin from seventeen mammalian species were compared. They showed a generally slow rate of evolution, but with marked variations in different lineages. Most substitutions have occurred in the C-terminal part of the chain, which probably forms part of the surface of the-crystallin aggregate. The ancestral sequence method of Dayhoff revealed interesting indications about the phylogenetic relationships between the eleven mammalian orders that were represented by the investigated species. Most evident was the divergence of marsupial and placental orders. A notable resemblance between the hyrax and elephant sequences was observed, setting them apart from the ungulates, including whale. Primates, rodents, lagomorphs, insectivores and tupaiids seem to derive from a common stem group. These phylogenetic inferences are discussed in relation to current palaeontological and taxonomical opinions, and compared to evidence from other protein sequence data.     

Conditions for Global Dynamic Stability of a Class of Resource-Bounded Model Ecosystems

This paper studies a class of dynamical systems that model multi-species ecosystems. These systems are ‘resource bounded’ in the sense that species compete to utilize an underlying limiting resource or substrate. This boundedness means that the relevant state space can be reduced to a simplex, with coordinates representing the proportions of substrate utilized by the various species. If the vector field is inward pointing on the boundary of the simplex, the state space is forward invariant under the system flow, a requirement that can be interpreted as the presence of non-zero exogenous recruitment. We consider conditions under which these model systems have a unique interior equilibrium that is globally asymptotically stable. The systems we consider generalize classical multi-species Lotka–Volterra systems, the behaviour of which is characterized by properties of the community (or interaction) matrix. However, the more general systems considered here are not characterized by a single matrix, but rather a family of matrices. We develop a set of ‘explicit conditions’ on the basis of a notion of ‘uniform diagonal dominance’ for such a family of matrices, that allows us to extract a set of sufficient conditions for global asymptotic stability based on properties of a single, derived matrix. Examples of these explicit conditions are discussed.     

Environmental Fluctuations and Level of Density-Compensation Strongly Affects the Probability of Fixation and Fixation Times

The probability of, and time to, fixation of a mutation in a population has traditionally been studied by the classic Wright–Fisher model where population size is constant. Recent theoretical expansions have covered fluctuating populations in various ways but have not incorporated models of how the environment fluctuates in combination with different levels of density-compensation affecting fecundity. We tested the hypothesis that the probability of, and time to, fixation of neutral, advantageous and deleterious mutations is dependent on how the environment fluctuates over time, and on the level of density-compensation. We found that fixation probabilities and times were dependent on the pattern of autocorrelation of carrying capacity over time and interacted with density-compensation. The pattern found was most pronounced at small population sizes. The patterns differed greatly depending on whether the mutation was neutral, advantageous, or disadvantageous. The results indicate that the degree of mismatch between carrying capacity and population size is a key factor, rather than population size per se, and that effective population sizes can be very low also when the census population size is far above the carrying capacity. This study highlights the need for explicit population dynamic models and models for environmental fluctuations for the understanding of the dynamics of genes in populations.     

Evolutionary diversification of Mesobuthus α-scorpion toxins affecting sodium channels

α-Scorpion toxins constitute a family of peptide modulators that induce a prolongation of the action potential of excitable cells by inhibiting voltage-gated sodium channel inactivation. Although they all adopt a conserved structural scaffold, the potency and phylogentic preference of these toxins largely vary, which render them an intriguing model for studying evolutionary diversification among family members. Here, we report molecular characterization of a new multigene family of α-toxins comprising 13 members (named MeuNaTxα-1 to MeuNaTxα-13) from the scorpion Mesobuthus eupeus. Of them, five native toxins (MeuNaTxα-1 to -5) were purified to homogeneity from the venom and the solution structure of MeuNaTxα-5 was solved by nuclear magnetic resonance. A systematic functional evaluation of MeuNaTxα-1, -2, -4, and -5 was conducted by two-electrode voltage-clamp recordings on seven cloned mammalian voltage-gated sodium channels (Na(v)1.2 to Na(v)1.8) and the insect counterpart DmNa(v)1 expressed in Xenopus oocytes. Results show that all these four peptides slow inactivation of DmNa(v)1 and are inactive on Na(v)1.8 at micromolar concentrations. However, they exhibit differential specificity for the other six channel isoforms (Na(v)1.2 to Na(v)1.7), in which MeuNaTxα-4 shows no activity on these isoforms and thus represents the first Mesobuthus-derived insect-selective α-toxin identified so far with a half maximal effective concentration of 130 ± 2 nm on DmNa(v)1 and a half maximal lethal dose of about 200 pmol g(-1) on the insect Musca domestica; MeuNaTxα-2 only affects Na(v)1.4; MeuNaTxα-1 and MeuNaTxα-5 have a wider range of channel spectrum, the former active on Na(v)1.2, Na(v)1.3, Na(v)1.6, and Na(v)1.7, whereas the latter acting on Na(v)1.3-Na(v)1.7. Remarkably, MeuNaTxα-4 and MeuNaTxα-5 are two nearly identical peptides differing by only one point mutation at site 50 (A50V) but exhibit rather different channel subtype selectivity, highlighting a switch role of this site in altering the target specificity. By the maximum likelihood models of codon substitution, we detected nine positively selected sites (PSSs) that could be involved in functional diversification of Mesobuthus α-toxins. The PSSs include site 50 and other seven sites located in functional surfaces of α-toxins. This work represents the first thorough investigation of evolutionary diversification of α-toxins derived from a specific scorpion lineage from the perspectives of sequence, structure, function, and evolution.     

Close Evolutionary Relatedness of α-Amylases from Archaea and Plants

The amino acid sequences of 22 α-amylases from family 13 of glycosyl hydrolases were analyzed with the aim of revealing the evolutionary relationships between the archaeal α-amylases and their eubacterial and eukaryotic counterparts. Two evolutionary distance trees were constructed: (i) the first one based on the alignment of extracted best-conserved sequence regions (58 residues) comprising β2, β3, β4, β5, β7, and β8 strand segments of the catalytic (α/β)₈-barrel and a short conserved stretch in domain B protruding out of the barrel in the β3 →α3 loop, and (ii) the second one based on the alignment of the substantial continuous part of the (α/β)₈-barrel involving the entire domain B (consensus length: 386 residues). With regard to archaeal α-amylases, both trees compared brought, in fact, the same results; i.e., all family 13 α-amylases from domain Archaea were clustered with barley pI isozymes, which represent all plant α-amylases. The enzymes from Bacillus licheniformis and Escherichia coli, representing liquefying and cytoplasmic α-amylases, respectively, seem to be the further closest relatives to archaeal α-amylases. This evolutionary relatedness clearly reflects the discussed similarities in the amino acid sequences of these α-amylases, especially in the best-conserved sequence regions. Since the results for α-amylases belonging to all three domains (Eucarya, Eubacteria, Archaea) offered by both evolutionary trees are very similar, it is proposed that the investigated conserved sequence regions may indeed constitute the ``sequence fingerprints' of a given α-amylase. Received: 3 June 1998 / Accepted: 20 August 1998     

Modeling and Estimation of Kinetic Parameters and Replicative Fitness of HIV-1 from Flow-Cytometry-Based Growth Competition Experiments

Growth competition assays have been developed to quantify the relative fitness of HIV-1 mutants. In this article, we develop mathematical models to describe viral/cellular dynamic interactions in the assay system from which the competitive fitness indices or parameters are defined. In our previous HIV-viral fitness experiments, the concentration of uninfected target cells was assumed to be constant (Wu et al. 2006). But this may not be true in some experiments. In addition, dual infection may frequently occur in viral fitness experiments and may not be ignorable. Here, we relax these two assumptions and extend our earlier viral fitness model (Wu et al. 2006). The resulting models then become nonlinear ODE systems for which closed-form solutions are not achievable. In the new model, the viral relative fitness is a function of time since it depends on the target cell concentration. First, we studied the structure identifiability of the nonlinear ODE models. The identifiability analysis showed that all parameters in the proposed models are identifiable from the flow-cytometry-based experimental data that we collected. We then employed a global optimization approach (the differential evolution algorithm) to directly estimate the kinetic parameters as well as the relative fitness index in the nonlinear ODE models using nonlinear least square regression based on the experimental data. Practical identifiability was investigated via Monte Carlo simulations.     

Serological identification and molecular characterization of B (A) 02 subtype in patients and blood donors from Eastern China

This study was designed to identify the rare?type?ABO?blood?groups, B(A) 02, from Eastern China. Three samples with discordant serological results during routine blood type identification and four samples from one sample’ family were selected. All of them were detected by serological method. The exon 6 and 7 of the ABO genes were amplified by PCR and sequenced. They were typed as AsubB by serology and as BO by genotype. In AsubB samples, nt 700C>G mutation was detected in B gene, which was previously defined as B(A)02 alleles. In these seven samples, six showed B(A)02/O01 and one showed B(A)02/O02.B(A)02 allele was found to be more common in this study than B(A)04 which is considered to be more frequent than B(A)02. The careful identification of rare blood types is important for the safety of clinical blood transfusion.     

Evolutionary loss of parasitism by nematodes? Discovery of a free-living filaroid nematode

A cattle-drinking pool in nature reserve "Zwin" on the Belgian coast contained free-living third-stage infective filaroid juveniles. These juveniles clearly differ morphologically from all known nematodes. Morphological and molecular analyses indicate a position within the Filaroidea. The aberrant biology of this nematode, namely, a free-living stage in an aquatic environment, is unknown within this superfamily, and the evolution of the parasitic phenotype to a free-living state is generally thought to be unlikely. However, the obtained placement in the small subunit molecular phylogenetic tree suggests that this free-living stage is most likely a secondary adaptation. It is reasonable to assert that nematodes with complex life cycles still have the genetic potential for a reversion from parasitism to a (partial) free-living stage.     

Evolutionary analysis of prokaryotic heat-shock transcription regulatory protein σ³²

Heat-stress to any living cell is known to trigger a universal defense response, called heat-shock response, with rapid induction of tens of different heat-shock proteins. Bacterial heat-shock genes are transcribed by the σ³²-bound RNA polymerase instead of the normal σ⁷⁰-bound RNA polymerase. In this study, the diversity in sequence, variation in secondary structure and function amongst the different functional regions of the proteobacterial σ³² family of proteins, and their phylogenetic relationships have been analyzed. Bacterial σ³² proteins can be subdivided into different functional regions which are referred to as regions 2, 3, and 4. There is a great deal of sequence conservation among the functional regions of proteobacterial σ³² family of proteins though some mutations are also present in these regions. Region 2 is the most conserved one, while region 4 has comparatively more variable sequences. In the present work, we tried to explore the effects of mutations in these regions. Our study suggests that the sequence diversities due to natural mutations in the different regions of proteobacterial σ³² family lead to different functions. So far, this study is the first bioinformatic approach towards the understanding of the mechanistic details of σ³² family of proteins using the protein sequence information only. This study therefore may help in elucidating the hitherto unknown molecular mechanism of the functionalities of σ³²family of proteins.     

The Role of Cohesiveness in the Permeability of the Spatial Assemblies of FG Nucleoporins

Nuclear pore complexes (NPCs) conduct selective, bidirectional transport across the nuclear envelope. The NPC passageway is lined by intrinsically disordered proteins that contain hydrophobic phenylalanine-glycine (FG) motifs, known as FG nucleoporins (FG nups), that play the key role in the NPC transport mechanism. Cohesive interactions among the FG nups, which arise from the combination of hydrophobic, electrostatic, and other forces, have been hypothesized to control the morphology of the assemblies of FG nups in the NPC, as well as their permeability with respect to the transport proteins. However, the role of FG nup cohesiveness is still vigorously debated. Using coarse-grained polymer theory and numerical simulations, we study the effects of cohesiveness on the selective permeability of in vitro FG nup assemblies in different geometries that have served as proxies for the morphological and transport properties of the NPC. We show that in high-density FG nup assemblies, increase in cohesiveness leads to the decrease in their permeability, in accordance with the accepted view. On the other hand, the permeability of low-density assemblies is a nonmonotonic function of the cohesiveness, and a moderate increase in cohesiveness can enhance permeability. The density- and cohesiveness-dependent effects on permeability are explained by considering the free-energy cost associated with penetrating the FG nup assemblies. We discuss the implications of these findings for the organization and function of the NPC.     

Evolutionary responses of foraging-related traits in unstable predator–prey systems

The evolutionary responses of predators to prey and of prey to predators are analysed using models for the dynamics of a quantitative trait that determines the capture rate of prey by an average searching predator. Unlike previous investigations, the analysis centres on models and/or parameter values for which the two-species equilibrium is locally unstable. The instability in some models is driven by the predators non-linear functional response to prey; in other models, the cycles are a direct consequence of evolutionary response to selection acting on the trait. When the values of predator and prey traits combine multiplicatively to determine the capture rate, the predators trait shows only a transient response to changes in the preys trait in stable systems. However, when the population densities exhibit sustained oscillations, predators often evolve an increased long-term mean capture rate in response to an increased prey escape ability. Under the multiplicative model, prey in stable systems always evolve increased escape ability in response to an increased predator capture a     

Evolutionary implications of Avian Infectious Bronchitis Virus （AIBV） analysis

For developing efficient vaccines, it is essential to identify which amino acid changes are most important to the survival of the virus. We investigate the amino acid substitution features in the Avian Infectious Bronchitis Virus （AIBV） antigenic domain of a vaccine serotype （DE072） and a virulent viral strain （GA98） to better understand adaptive evolution of AIBV. In addition, the SARS Coronavirus （SARS-CoV） was also analyzed in the same way. It is interesting to find that extreme comparability exists between AIBV and SARS in amino acid substitution pattern. It suggests that amino acid changes that result in overall shift of residue charge and polarity should be paid special attention to during the development of vaccines.     

Bistability Analysis of an Apoptosis Model in the Presence of Nitric Oxide

Bistability in apoptosis, or programmed cell death, is crucial for the healthy functioning of multicellular organisms. The aim in this study is to show the presence of bistability in a mitochondria-dependent apoptosis model under nitric oxide effects using chemical reaction network theory. The model equations are a set of coupled ordinary differential equations arising from the assumed mass action kinetics. Whether these equations have a capacity for bistability (cell survival and apoptosis) is determined using a modular approach in which the model is decomposed into modules. Each module contains only a subset of the whole model and is analyzed separately. It is seen that bistability in a module is preserved throughout the whole model after adding the remaining reactions in the pathway on these modules. It is also found that inhibitor effect of some proteins and the appearance of a reacting protein in a later stage as a product is a desired feature but not sufficient for bistability (in the absence of cooperativity effects). On the whole model, two apoptotic and two cell survival states are obtained depending on the initial cell conditions. The results suggest that the antiapoptotic effects of nitric oxide species are responsible for the bistable character of the apoptotic pathway when cooperativity is not assumed in the apoptosome formation.     

Evolutionary History of Eukaryotic α-Glucosidases from the α-Amylase Family

Although some α-glucosidases from the α-amylase family (glycoside hydrolase family GH13) have been studied extensively, their exact number, organization on the chromosome, and orthology/paralogy relationship were unknown. This was true even for important disease vectors where gut α-glucosidase is known to be receptor for the Bin toxin used to control the population of some mosquito species. In some cases orthologs from related species were studied intensively, while potentially important paralogs were omitted. We have, therefore, used a bioinformatics approach to identify all family GH13 α-glucosidases from the selected species from Metazoa (including three mosquito species: Aedes aegypti, Anopheles gambiae, and Culex quinquefasciatus) as well as from Fungi in an effort to characterize their arrangement on the chromosome and evolutionary relationships among orthologs and among paralogs. We also searched for pseudogenes and genes coding for enzymatically inactive proteins with a possible new function. We have found GH13 α-glucosidases mostly in Arthropoda and Fungi where they form gene families, as a result of multiple lineage-specific gene duplications. In mosquito species we have identified 14 α-glucosidase (Aglu) genes of which only five have been biochemically characterized so far, two are putative pseudogenes and the rest remains uncharacterized. We also revealed quite a complex evolutionary history of the eukaryotic α-glucosidases probably involving multiple losses of genes or horizontal gene transfer from bacteria.