Geographic variation and positive selection on M7 lysin, an acrosomal sperm protein in mussels (Mytilus spp.)

Successful fertilization in free-spawning marine organisms depends on the interactions between genes expressed on the surfaces of eggs and sperm. Positive selection frequently characterizes the molecular evolution of such genes, raising the possibility that some common deterministic process drives the evolution of gamete recognition genes and may even be important for understanding the evolution of prezygotic isolation and speciation in the marine realm. One hypothesis is that gamete recognition genes are subject to selection for prezygotic isolation, namely reinforcement. In a previous study, positive selection on the gene coding for the acrosomal sperm protein M7 lysin was demonstrated among allopatric populations of mussels in the Mytilus edulis species group (M. edulis, M. galloprovincialis, and M. trossulus). Here, we expand sampling to include M7 lysin haplotypes from populations where mussel species are sympatric and hybridize to determine whether there is a pattern of reproductive character displacement, which would be consistent with reinforcement driving selection on this gene. We do not detect a strong pattern of reproductive character displacement; there are no unique haplotypes in sympatry nor is there consistently greater population structure in comparisons involving sympatric populations. One distinct group of haplotypes, however, is strongly affected by natural selection and this group of haplotypes is found within M. galloprovincialis populations throughout the Northern Hemisphere concurrent with haplotypes common to M. galloprovincialis and M. edulis. We suggest that balancing selection, perhaps resulting from sexual conflicts between sperm and eggs, maintains old allelic diversity within M. galloprovincialis.

Patterns of variation within self-incompatibility loci

Diverse self-incompatibility (SI) mechanisms permit flowering plants to inhibit fertilization by pollen that express specificities in common with the pistil. Characteristic of at least two model systems is greatly reduced recombination across large genomic tracts surrounding the S-locus, which regulates SI. In three angiosperm families, including the Solanaceae, the gene that controls the expression of gametophytic SI in the pistil encodes a ribonuclease (S-RNase). The gene that controls pollen SI expression is currently unknown, although several candidates have recently been proposed. Although each candidate shows a high level of polymorphism and complete allelic disequilibrium with the S-RNase gene, such properties may merely reflect tight linkage to the S-locus, irrespective of any functional role in SI. We analyzed the magnitude and nature of nucleotide variation, with the objective of distinguishing likely candidates for regulators of SI from other genes embedded in the S-locus region. We studied the S-RNase gene of the Solanaceae and 48A, a candidate for the pollen gene in this system, and we also conducted a parallel analysis of the regulators of sporophytic SI in Brassica, a system in which both the pistil and pollen genes are known. Although the pattern of variation shown by the pollen gene of the Brassica system is consistent with its role as a determinant of pollen specificity, that of 48A departs from expectation. Our analysis further suggests that recombination between 48A and S-RNase may have occurred during the interval spanned by the gene genealogy, another indication that 48A may not regulate SI expression in pollen.

Family Size Evolution in Drosophila Chemosensory Gene Families: A Comparative Analysis with a Critical Appraisal of Methods

Gene turnover rates and the evolution of gene family sizes are important aspects of genome evolution. Here, we use curated sequence data of the major chemosensory gene families from Drosophila-the gustatory receptor, odorant receptor, ionotropic receptor, and odorant-binding protein families-to conduct a comparative analysis among families, exploring different methods to estimate gene birth and death rates, including an ad hoc simulation study. Remarkably, we found that the state-of-the-art methods may produce very different rate estimates, which may lead to disparate conclusions regarding the evolution of chemosensory gene family sizes in Drosophila. Among biological factors, we found that a peculiarity of D. sechellia's gene turnover rates was a major source of bias in global estimates, whereas gene conversion had negligible effects for the families analyzed herein. Turnover rates vary considerably among families, subfamilies, and ortholog groups although all analyzed families were quite dynamic in terms of gene turnover. Computer simulations showed that the methods that use ortholog group information appear to be the most accurate for the Drosophila chemosensory families. Most importantly, these results reveal the potential of rate heterogeneity among lineages to severely bias some turnover rate estimation methods and the need of further evaluating the performance of these methods in a more diverse sampling of gene families and phylogenetic contexts. Using branch-specific codon substitution models, we find further evidence of positive selection in recently duplicated genes, which attests to a nonneutral aspect of the gene birth-and-death process.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a 3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.

Rapid evolution in plant chitinases: Molecular targets of selection in plant-pathogen coevolution

Many pathogen recognition genes, such as plant R-genes, undergo rapid adaptive evolution, providing evidence that these genes play a critical role in plant-pathogen coevolution. Surprisingly, whether rapid adaptive evolution also occurs in genes encoding other kinds of plant defense proteins is unknown. Unlike recognition proteins, plant chitinases attack pathogens directly, conferring disease resistance by degrading chitin, a component of fungal cell walls. Here, we show that nonsynonymous substitution rates in plant class I chitinase often exceed synonymous rates in the plant genus Arabis (Cruciferae) and in other dicots, indicating a succession of adaptively driven amino acid replacements. We identify individual residues that are likely subject to positive selection by using codon substitution models and determine the location of these residues on the three-dimensional structure of class I chitinase. In contrast to primate lysozymes and plant class III chitinases, structural and functional relatives of class I chitinase, the adaptive replacements of class I chitinase occur disproportionately in the active site cleft. This highly unusual pattern of replacements suggests that fungi directly defend against chitinolytic activity through enzymatic inhibition or other forms of chemical resistance and identifies target residues for manipulating chitinolytic activity. These data also provide empirical evidence that plant defense proteins not involved in pathogen recognition also evolve in a manner consistent with rapid coevolutionary interactions.

A generalized mechanistic codon model.

Models of codon evolution have attracted particular interest because of their unique capabilities to detect selection forces and their high fit when applied to sequence evolution. We described here a novel approach for modeling codon evolution, which is based on Kronecker product of matrices. The 61 × 61 codon substitution rate matrix is created using Kronecker product of three 4 × 4 nucleotide substitution matrices, the equilibrium frequency of codons, and the selection rate parameter. The entities of the nucleotide substitution matrices and selection rate are considered as parameters of the model, which are optimized by maximum likelihood. Our fully mechanistic model allows the instantaneous substitution matrix between codons to be fully estimated with only 19 parameters instead of 3,721, by using the biological interdependence existing between positions within codons. We illustrate the properties of our models using computer simulations and assessed its relevance by comparing the AICc measures of our model and other models of codon evolution on simulations and a large range of empirical data sets. We show that our model fits most biological data better compared with the current codon models. Furthermore, the parameters in our model can be interpreted in a similar way as the exchangeability rates found in empirical codon models.

Mapping the molecular characteristics of Brazilian human T-cell lymphotropic virus type 1 Env (gp46) and Pol amino acid sequences for vaccine design

This study was carried out to evaluate the molecular pattern of all available Brazilian human T-cell lymphotropic virus type 1 Env (n = 15) and Pol (n = 43) nucleotide sequences via epitope prediction, physico-chemical analysis, and protein potential sites identification, giving support to the Brazilian AIDS vaccine program. In 12 previously described peptides of the Env sequences we found 12 epitopes, while in 4 peptides of the Pol sequences we found 4 epitopes. The total variation on the amino acid composition was 9 and 17% for human leukocyte antigen (HLA) class I and class II Env epitopes, respectively. After analyzing the Pol sequences, results revealed a total amino acid variation of 0.75% for HLA-I and HLA-II epitopes. In 5 of the 12 Env epitopes the physico-chemical analysis demonstrated that the mutations magnified the antigenicity profile. The potential protein domain analysis of Env sequences showed the loss of a CK-2 phosphorylation site caused by D197N mutation in one epitope, and a N-glycosylation site caused by S246Y and V247I mutations in another epitope. Besides, the analysis of selection pressure have found 8 positive selected sites (w = 9.59) using the codon-based substitution models and maximum-likelihood methods. These studies underscore the importance of this Env region for the virus fitness, for the host immune response and, therefore, for the development of vaccine candidates.

Divergent evolution and purifying selection of the flaA gene sequences in Aeromonas

BACKGROUND: The bacterial flagellum is the most important organelle of motility in bacteria and plays a key role in many bacterial lifestyles, including virulence. The flagellum also provides a paradigm of how hierarchical gene regulation, intricate protein-protein interactions and controlled protein secretion can result in the assembly of a complex multi-protein structure tightly orchestrated in time and space. As if to stress its importance, plants and animals produce receptors specifically dedicated to the recognition of flagella. Aside from motility, the flagellum also moonlights as an adhesion and has been adapted by humans as a tool for peptide display. Flagellar sequence variation constitutes a marker with widespread potential uses for studies of population genetics and phylogeny of bacterial species. RESULTS: We sequenced the complete flagellin gene (flaA) in 18 different species and subspecies of Aeromonas. Sequences ranged in size from 870 (A. allosaccharophila) to 921 nucleotides (A. popoffii). The multiple alignment displayed 924 sites, 66 of which presented alignment gaps. The phylogenetic tree revealed the existence of two groups of species exhibiting different FlaA flagellins (FlaA1 and FlaA2). Maximum likelihood models of codon substitution were used to analyze flaA sequences. Likelihood ratio tests suggested a low variation in selective pressure among lineages, with an omega ratio of less than 1 indicating the presence of purifying selection in almost all cases. Only one site under potential diversifying selection was identified (isoleucine in position 179). However, 17 amino acid positions were inferred as sites that are likely to be under positive selection using the branch-site model. Ancestral reconstruction revealed that these 17 amino acids were among the amino acid changes detected in the ancestral sequence. CONCLUSION: The models applied to our set of sequences allowed us to determine the possible evolutionary pathway followed by the flaA gene in Aeromonas, suggesting that this gene have probably been evolving independently in the two groups of Aeromonas species since the divergence of a distant common ancestor after one or several episodes of positive selection. REVIEWERS: This article was reviewed by Alexey Kondrashov, John Logsdon and Olivier Tenaillon (nominated by Laurence D Hurst).

Mitochondrial sequence data and Dipsacales phylogeny: Mixed models, partitioned Bayesian analyses, and model selection

Phylogenetic analyses of chloroplast DNA sequences, morphology, and combined data have provided consistent support for many of the major branches within the angiosperm, clade Dipsacales. Here we use sequences from three mitochondrial loci to test the existing broad scale phylogeny and in an attempt to resolve several relationships that have remained uncertain. Parsimony, maximum likelihood, and Bayesian analyses of a combined mitochondrial data set recover trees broadly consistent with previous studies, although resolution and support are lower than in the largest chloroplast analyses. Combining chloroplast and mitochondrial data results in a generally well-resolved and very strongly supported topology but the previously recognized problem areas remain. To investigate why these relationships have been difficult to resolve we conducted a series of experiments using different data partitions and heterogeneous substitution models. Usually more complex modeling schemes are favored regardless of the partitions recognized but model choice had little effect on topology or support values. In contrast there are consistent but weakly supported differences in the topologies recovered from coding and non-coding matrices. These conflicts directly correspond to relationships that were poorly resolved in analyses of the full combined chloroplast-mitochondrial data set. We suggest incongruent signal has contributed to our inability to confidently resolve these problem areas. (c) 2007 Elsevier Inc. All rights reserved.

Inflation, currency substitution and dollarization: The case of Argentina

This dissertation discusses the relationship between inflation, currency substitution and dollarization that has taken place in Argentina for the past several decades.^ First, it is shown that when consumers are able to hold only domestic monetary balances (without capital mobility) an increase in the rate of inflation will produce a balance of payments deficit. We then look at the same issue but with heterogeneous consumers, this heterogeneity being generated by non-proportional lump-sum transfers.^ Second, we discussed some necessary assumptions related to currency substitution models and concluded that there was no a-priori conclusion on whether currencies should be assumed to be "cooperant" or "non-cooperant" in utility. That is to say, whether individuals held different currencies together or one instead of the other.^ Third, we went into discussing the issue of currency substitution as being a constraint on governments' inflationary objectives rather than a choice of those governments to avoid hyperinflations. We showed that imperfect substitutability between currencies does not "reduce the scope for rational (hyper)inflationary processes" as it had been previously argued. It will ultimately depend on the parametrization used and not on the intrinsic characteristics of imperfect substitutability between currencies.^ We further showed that in Argentina, individuals have been able to endogenize the money supply by holding foreign monetary balances. We argued that the decision to hold foreign monetary balances by individuals is always a second best due to the trade-off between holding foreign monetary balances and consumption. For some levels of income, consumption, and foreign inflation, individuals would prefer to hold domestic monetary balances rather than foreign ones.^ We then modeled the distinction between dollarization and currency substitution. We concluded that although dollarization is necessary for currency substitution to take place, the decision to use foreign monetary balances for transactions purposes is largely independent from the dollarization process.^ Finally, we concluded that Argentina should not fully dollarize its economy because dollarization is always a second best to using a domestic currency. Further, we argued that a fixed exchange system would be better than a flexible exchange rate or a "crawling-peg" system because of the characteristics of the political system and the possibilities of "mass praetorianism" to develop, which is intricately linked to "populist" solutions. ^

Inflation, currency substitution, and dollarization : the case of Argentina

This dissertation discusses the relationship between inflation, currency substitution and dollarization that has taken place in Argentina for the past several decades. First, it is shown that when consumers are able to hold only domestic monetary balances (without capital mobility) an increase in the rate of inflation will produce a balance of payments deficit. We then look at the same issue but with heterogeneous consumers, this heterogeneity being generated by non-proportional lump-sum transfers. Second, we discussed some necessary assumptions related to currency substitution models and concluded that there was no a-priori conclusion on whether currencies should be assumed to be "cooperant" or "non-cooperant" in utility. That is to say, whether individuals held different currencies together or one instead of the other. Third, we went into discussing the issue of currency substitution as being a constraint on governments inflationary objectives rather than a choice of those governments to avoid hyperinflations. We showed that imperfect substitutability between currencies does not "reduce the scope for rational (hyper)inflationary processes" as it had been previously argued. It will ultimately depend on the parametrization used and not on the intrinsic characteristics of imperfect substitutability between currencies. We further showed that in Argentina, individuals have been able to endogenize the money supply by holding foreign monetary balances. We argued that the decision to hold foreign monetary balances by individuals is always a second best due to the trade-off between holding foreign monetary balances and consumption. For some levels of income, consumption, and foreign inflation, individuals would prefer to hold domestic monetary balances rather than foreign ones. We then modeled the distinction between dollarization and currency substitution. We concluded that although dollarization is necessary for currency substitution to take place, the decision to use foreign monetary balances for transactions purposes is largely independent from the dollarization process. Finally, we concluded that Argentina should not fully dollarize its economy because dollarization is always a second best to using a domestic currency. Further, we argued that a fixed exchange system would be better than a flexible exchange rate or a "crawling-peg" system because of the characteristics of the political system and the possibilities of "mass praetorianism" to develop, which is intricately linked to "populist" solutions.

Towards an assessment of character interdependence in avian RNA phylogenetics: A general secondary structure model for the avian mitochondrial 16S rRNA

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Toward a resolution of Campanulid phylogeny, with special reference to the placement of Dipsacales

Broad-scale phylogenetic analyses of the angiosperms and of the Asteridae have failed to confidently resolve relationships among the major lineages of the campanulid Asteridae (i.e., the euasterid II of APG II, 2003). To address this problem we assembled presently available sequences for a core set of 50 taxa, representing the diversity of the four largest lineages (Apiales, Aquifoliales, Asterales, Dipsacales) as well as the smaller ""unplaced"" groups (e.g., Bruniaceae, Paracryphiaceae, Columelliaceae). We constructed four data matrices for phylogenetic analysis: a chloroplast coding matrix (atpB, matK, ndhF, rbcL), a chloroplast non-coding matrix (rps16 intron, trnT-F region, trnV-atpE IGS), a combined chloroplast dataset (all seven chloroplast regions), and a combined genome matrix (seven chloroplast regions plus 18S and 26S rDNA). Bayesian analyses of these datasets using mixed substitution models produced often well-resolved and supported trees. Consistent with more weakly supported results from previous studies, our analyses support the monophyly of the four major clades and the relationships among them. Most importantly, Asterales are inferred to be sister to a clade containing Apiales and Dipsacales. Paracryphiaceae is consistently placed sister to the Dipsacales. However, the exact relationships of Bruniaceae, Columelliaceae, and an Escallonia clade depended upon the dataset. Areas of poor resolution in combined analyses may be partly explained by conflict between the coding and non-coding data partitions. We discuss the implications of these results for our understanding of campanulid phylogeny and evolution, paying special attention to how our findings bear on character evolution and biogeography in Dipsacales.

Competitiveness and labour regulations

Includes bibliography