Relative Substitution Rates Research Articles

Beyond Ramsey: Gender-Based Taxation with Non-Cooperative Couples

This paper explores the implications of gender-based income taxation in a non- cooperative model of a couple’s time allocation between market work and providing a household public good. We find that the optimal structure of differential taxation by gender is solely determined by spouses’ relative marginal rates of substitution between the public household good and private consumption. Breaking down this general rule into the primitives of the model, the spouse with a lower personal valuation of the public household good should be taxed at a higher rate. If these valuations are identical, a comparative advantage in home production relative to market work will imply a higher marginal tax rate. Using a realistic calibration, we show that these two results may combine to imply a higher optimal tax rate on female labor supply. This result stands in sharp contrast to previous models of gender-based taxation in which households select Pareto efficient allocations. Extending the model to include altruistic preferences, leisure, or human capital accumulation reduces optimal tax rates, while sequential labor supply decisions affect the optimal tax rate of the primary earner in an ambiguous direction.

Fast Side Chain Replacement in Proteins Using a Coarse-Grained Approach for Evaluating the Effects of Mutation During Evolution

For high-throughput structural genomic and evolutionary bioinformatics approaches, there is a clear need for fast methods to evaluate substitutions structurally. Coarse-grained methods are both powerful and fast, and a coarse-grained approach to position the substituted side chains is presented. Through the application of a coarse-grained method, a speed-up on the single- residue replacement, of at least sevenfold is achieved compared with modern all-atom approaches. At the same time, this approach maintains a small median RMSD from the leading all-atom approach (as measured in coarse-grained space), and predicts the conformation of point mutants with similar accuracy and generates biologically realistic side chain angles. This method is also substantially more predictable in its run time, making it useful for high-throughput studies of protein structural evolution. To demonstrate the utility of this method, it has been implemented in a forward simulation of sequences threaded through the SH2 domains, with selective pressures to fold and bind specifically. The relative substitution rates across the protein structure and at the binding interface are reflective of those observed in SH2 domain evolution. The algorithm has been implemented in C++, with the source code and binaries (currently supported for Linux systems) freely available as SARA at http://www.wyomingbioinformatics.org/LiberlesGroup/SARA .

Risk, ambiguity, and state-preference theory

The state-preference framework for modeling choice under uncertainty, in which objects of choice are allocations of wealth or commodities across states of the world, is a natural one for modeling “smooth” ambiguity-averse preferences. It does not require reference to objective probabilities, personalistic consequences, or counterfactual acts, and it allows for state dependence of utility and unobservable background risk. The decision maker’s local revealed beliefs are encoded in her risk-neutral probabilities (her relative marginal rates of substitution between states) and her local risk preferences are encoded in the matrix of derivatives of the risk-neutral probabilities. This matrix plays a central but generally unappreciated role in the modeling of risk attitudes in the state-preference framework. It can be computed by inverting a bordered Slutsky matrix and vice versa, it generalizes the Arrow–Pratt measure for approximating local risk premia, and its structure reveals whether the decision maker’s risk preferences are ambiguity averse as well as risk-averse. Two versions of the smooth ambiguity model are analyzed—the source-dependent risk aversion model and the second-order uncertainty (KMM) model—and it is shown that in both cases, the overall premium for local uncertainty can be decomposed as the sum of a risk premium and an ambiguity premium.

A Molecular Clock for Malaria Parasites

The evolutionary origins of new lineages of pathogens are fundamental to understanding emerging diseases. Phylogenetic reconstruction based on DNA sequences has revealed the sister taxa of human pathogens, but the timing of host-switching events, including the human malaria pathogen Plasmodium falciparum, remains controversial. Here, we establish a rate for cytochrome b evolution in avian malaria parasites relative to its rate in birds. We found that the parasite cytochrome b gene evolves about 60% as rapidly as that of host cytochrome b, corresponding to approximately 1.2% sequence divergence per million years. This calibration puts the origin of P. falciparum at 2.5 million years ago (Ma), the initial radiation of mammalian Plasmodium at 12.8 Ma, and the contemporary global diversity of the Haemosporida across terrestrial vertebrates at 16.2 Ma.

Difference in larval type explains patterns of nonsynonymous substitutions in two ancient paralogs of the histone H3 gene in sea stars

Paralogous genes frequently show differences in patterns and rates of substitution that are typically attributed to different selection regimes, mutation rates, or local recombination rates. Here, two anciently diverged paralogous copies of the histone H3 gene in sea stars, the tandem-repetitive early-stage gene and a newly isolated gene with lower copy number that was termed the "putative late-stage histone H3 gene" were analyzed in 69 species with varying mode of larval development. The two genes showed differences in relative copy number, overall substitution rates, nucleotide composition, and codon usage, but similar patterns of relative nonsynonymous substitution rates, when analyzed by the d(N)/d(S) ratio. Sea stars with a nonpelagic and nonfeeding larval type (i.e., brooding lineages) were observed to have d(N)/d(S) ratios that were larger than for nonbrooders but equal between the two paralogs. This finding suggested that demographic differences between brooding and nonbrooding lineages were responsible for the elevated d(N)/d(S) ratios observed for brooders and refuted a suggestion from a previous analysis of the early-stage gene that the excess nonsynonymous substitutions were due to either (1) gene expression differences at the larval stage between brooders and nonbrooders or (2) the highly repetitive structure of the early-stage histone H3 gene.

Uncorrected Nucleotide Bias in mtDNA Can Mimic the Effects of Positive Darwinian Selection

The relative rates of nucleotide substitution at synonymous and nonsynonymous sites within protein-coding regions have been widely used to infer the action of natural selection from comparative sequence data. It is known, however, that mutational and repair biases can affect rates of evolution at both synonymous and nonsynonymous sites. More importantly, it is also known that synonymous sites are particularly prone to the effects of nucleotide bias. This means that nucleotide biases may affect the calculated ratio of substitution rates at synonymous and nonsynonymous sites. Using a large data set of animal mitochondrial sequences, we demonstrate that this is, in fact, the case. Highly biased nucleotide sequences are characterized by significantly elevated dN/dS ratios, but only when the nucleotide frequencies are not taken into account. When the analysis is repeated taking the nucleotide frequencies at each codon position into account, such elevated ratios disappear. These results suggest that the recently reported differences in dN/dS ratios between vertebrate and invertebrate mitochondrial sequences could be explained by variations in mitochondrial nucleotide frequencies rather than the effects of positive Darwinian selection.

Acid-Catalyzed Enolization of β-Tetralone

This experiment allows students to use 1H NMR to directly compare the relative initial rates of substitution of the benzylic and non-benzylic α hydrogens of β-tetralone and correlate their findings with the predictions made by resonance theory. The experiment demonstrates that the benzylic hydrogens undergo α substitution much faster than the non-benzylic hydrogens. The described experiment can be used in the introductory two-semester organic chemistry course sequence to make the concepts of enolization and resonance concrete.

Testing candidate plant barcode regions in the Myristicaceae

The concept and practice of DNA barcoding have been designed as a system to facilitate species identification and recognition. The primary challenge for barcoding plants has been to identify a suitable region on which to focus the effort. The slow relative nucleotide substitution rates of plant mitochondria and the technical issues with the use of nuclear regions have focused attention on several proposed regions in the plastid genome. One of the challenges for barcoding is to discriminate closely related or recently evolved species. The Myristicaceae, or nutmeg family, is an older group within the angiosperms that contains some recently evolved species providing a challenging test for barcoding plants. The goal of this study is to determine the relative utility of six coding (Universal Plastid Amplicon - UPA, rpoB, rpoc1, accD, rbcL, matK) and one noncoding (trnH-psbA) chloroplast loci for barcoding in the genus Compsoneura using both single region and multiregion approaches. Five of the regions we tested were predominantly invariant across species (UPA, rpoB, rpoC1, accD, rbcL). Two of the regions (matK and trnH-psbA) had significant variation and show promise for barcoding in nutmegs. We demonstrate that a two-gene approach utilizing a moderately variable region (matK) and a more variable region (trnH-psbA) provides resolution among all the Compsonuera species we sampled including the recently evolved C. sprucei and C. mexicana. Our classification analyses based on nonmetric multidimensional scaling ordination, suggest that the use of two regions results in a decreased range of intraspecific variation relative to the distribution of interspecific divergence with 95% of the samples correctly identified in a sequence identification analysis.

Phylogenetic analysis of the corbiculate bee tribes based on 12 nuclear protein-coding genes (Hymenoptera: Apoidea: Apidae)

The corbiculate bees comprise four tribes, the advanced eusocial Apini and Meliponini, the primitively eusocial Bombini, and the solitary or communal Euglossini. Recovering a robust phylogeny for the four tribes is of considerable importance for understanding the evolution of eusociality, yet previous morphological and molecular studies reached strikingly different conclusions. We study an expanded data set consisting of 12 nuclear genes to explore lines of support for the molecular hypothesis. Results corroborate previous molecular studies; support increases as more genes are added. Across genes, support for the molecular hypothesis is positively correlated with the number of informative sites and the relative substitution rate. Phylogenetic signals supporting the molecular tree rest almost entirely upon synonymous changes at the first and third codon positions. We discuss possible future approaches for resolving the frustratingly persistent corbiculate bee controversy.

When is Vendor Managed Inventory Good for the Retailer? Impact of Relative Margins and Substitution Rates

When customers cannot find a particular item at a retailer because it is out of stock, they are likely, with some probability, to switch to a substitute product from another manufacturer at the same store. Analyzing a two-product full substitution case, this paper examines two beliefs argued in the literature: (1) that, under Vendor Managed Inventory (VMI), the retailers would benefit because manufacturers would increase stocking quantities to avoid losing sales to a competitor and (2) that substitution benefits retailers who make a sale regardless. We find that the first proposition, while appealing, is simply not true for many cases. We also find that the second proposition does not hold for a wide number of cases. We contribute to the understanding of the inherent tradeoffs involved in deciding to use RMI or VMI in the presence of competing, substitute products.

Relative Rates of Evolution in the Coding and Control Regions of African mtDNAs

Reduced median networks of African haplogroup L mitochondrial DNA (mtDNA) sequences were analyzed to determine the pattern of substitutions in both the noncoding control and coding regions. In particular, we attempted to determine the causes of the previously reported (Howell et al. 2004) violation of the molecular clock during the evolution of these sequences. In the coding region, there was a significantly higher rate of substitution at synonymous sites than at nonsynonymous sites as well as in the tRNA and rRNA genes. This is further evidence for the operation of purifying selection during human mtDNA evolution. For most sites in the control region, the relative rate of substitution was similar to the rate of neutral evolution (assumed to be most closely approximated by the substitution rate at 4-fold degenerate sites). However, there are a number of mutational hot spots in the control region, approximately 3% of the total sites, that have a rate of substitution greater than the neutral rate, at some sites by more than an order of magnitude. It is possible either that these sites are evolving under conditions of positive selection or that the substitution rate at some sites in the control region is strongly dependent upon sequence context. Finally, we obtained preliminary evidence for "nonideal" evolution in the control region, including haplogroup-specific substitution patterns and a decoupling between relative rates of substitution in the control and coding regions.

Different Mutations in the Genome-Linked Protein VPg of Potato virus Y Confer Virulence on the pvr23 Resistance in Pepper

Five different amino acid substitutions in the VPg of Potato virus Y were shown to be independently responsible for virulence toward pvr2(3) resistance gene of pepper. A consequence of these multiple mutations toward virulence involving single nucleotide substitutions is a particularly high frequency of resistance breaking (37% of inoculated plants from the first inoculation) and suggests a potentially low durability of pvr2(3) resistance. These five mutants were observed with significantly different frequencies, one of them being overrepresented. Genetic drift alone could not explain the observed distribution of virulent mutants. More plausible scenarios were obtained by taking into account either the relative substitution rates, the relative fitness of the mutants in pvr2(3) pepper plants, or both.

Analysis of the overdispersed clock in the short-term evolution of hepatitis C virus: Using the E1/E2 gene sequences to infer infection dates in a single source outbreak.

The assumption of a molecular clock for dating events from sequence information is often frustrated by the presence of heterogeneity among evolutionary rates due, among other factors, to positively selected sites. In this work, our goal is to explore methods to estimate infection dates from sequence analysis. One such method, based on site stripping for clock detection, was proposed to unravel the clocklike molecular evolution in sequences showing high variability of evolutionary rates and in the presence of positive selection. Other alternatives imply accommodating heterogeneity in evolutionary rates at various levels, without eliminating any information from the data. Here we present the analysis of a data set of hepatitis C virus (HCV) sequences from 24 patients infected by a single individual with known dates of infection. We first used a simple criterion of relative substitution rate for site removal prior to a regression analysis. Time was regressed on maximum likelihood pairwise evolutionary distances between the sequences sampled from the source individual and infected patients. We show that it is indeed the fastest evolving sites that disturb the molecular clock and that these sites correspond to positively selected codons. The high computational efficiency of the regression analysis allowed us to compare the site-stripping scheme with random removal of sites. We demonstrate that removing the fast-evolving sites significantly increases the accuracy of estimation of infection times based on a single substitution rate. However, the time-of-infection estimations improved substantially when a more sophisticated and computationally demanding Bayesian method was used. This method was used with the same data set but keeping all the sequence positions in the analysis. Consequently, despite the distortion introduced by positive selection on evolutionary rates, it is possible to obtain quite accurate estimates of infection dates, a result of especial relevance for molecular epidemiology studies.

Evolutionary Rates and Expression Level in Chlamydomonas

In many biological systems, especially bacteria and unicellular eukaryotes, rates of synonymous and nonsynonymous nucleotide divergence are negatively correlated with the level of gene expression, a phenomenon that has been attributed to natural selection. Surprisingly, this relationship has not been examined in many important groups, including the unicellular model organism Chlamydomonas reinhardtii. Prior to this study, comparative data on protein-coding sequences from C. reinhardtii and its close noninterfertile relative C. incerta were very limited. We compiled and analyzed protein-coding sequences for 67 nuclear genes from these taxa; the sequences were mostly obtained from the C. reinhardtii EST database and our C. incerta EST data. Compositional and synonymous codon usage biases varied among genes within each species but were highly correlated between the orthologous genes of the two species. Relative rates of synonymous and nonsynonymous substitution across genes varied widely and showed a strong negative correlation with the level of gene expression estimated by the codon adaptation index. Our comparative analysis of substitution rates in introns of lowly and highly expressed genes suggests that natural selection has a larger contribution than mutation to the observed correlation between evolutionary rates and gene expression level in Chlamydomonas.

The evolutionary radiation of Arvicolinae rodents (voles and lemmings): relative contribution of nuclear and mitochondrial DNA phylogenies

BackgroundMitochondrial and nuclear genes have generally been employed for different purposes in molecular systematics, the former to resolve relationships within recently evolved groups and the latter to investigate phylogenies at a deeper level. In the case of rapid and recent evolutionary radiations, mitochondrial genes like cytochrome b (CYB) are often inefficient for resolving phylogenetic relationships. One of the best examples is illustrated by Arvicolinae rodents (Rodentia; Muridae), the most impressive mammalian radiation of the Northern Hemisphere which produced voles, lemmings and muskrats. Here, we compare the relative contribution of a nuclear marker – the exon 10 of the growth hormone receptor (GHR) gene – to the one of the mitochondrial CYB for inferring phylogenetic relationships among the major lineages of arvicoline rodents.ResultsThe analysis of GHR sequences improves the overall resolution of the Arvicolinae phylogeny. Our results show that the Caucasian long-clawed vole (Prometheomys schaposnikowi) is one of the basalmost arvicolines, and confirm that true lemmings (Lemmus) and collared lemmings (Dicrostonyx) are not closely related as suggested by morphology. Red-backed voles (Myodini) are found as the sister-group of a clade encompassing water vole (Arvicola), snow vole (Chionomys), and meadow voles (Microtus and allies). Within the latter, no support is recovered for the generic recognition of Blanfordimys, Lasiopodomys, Neodon, and Phaiomys as suggested by morphology. Comparisons of parameter estimates for branch lengths, base composition, among sites rate heterogeneity, and GTR relative substitution rates indicate that CYB sequences consistently exhibit more heterogeneity among codon positions than GHR. By analyzing the contribution of each codon position to node resolution, we show that the apparent higher efficiency of GHR is due to their third positions. Although we focus on speciation events spanning the last 10 million years (Myr), CYB sequences display highly saturated codon positions contrary to the nuclear exon. Lastly, variable length bootstrap predicts a significant increase in resolution of arvicoline phylogeny through the sequencing of nuclear data in an order of magnitude three to five times greater than the size of GHR exon 10.ConclusionOur survey provides a first resolved gene tree for Arvicolinae. The comparison of CYB and GHR phylogenetic efficiency supports recent assertions that nuclear genes are useful for resolving relationships of recently evolved animals. The superiority of nuclear exons may reside both in (i) less heterogeneity among sites, and (ii) the presence of highly informative sites in third codon positions, that evolve rapidly enough to accumulate synapomorphies, but slow enough to avoid substitutional saturation.

Adaptive Molecular Evolution in the Opsin Genes of Rapidly Speciating Cichlid Species

Cichlid fish inhabit a diverse range of environments that vary in the spectral content of light available for vision. These differences should result in adaptive selective pressure on the genes involved in visual sensitivity, the opsin genes. This study examines the evidence for differential adaptive molecular evolution in East African cichlid opsin genes due to gross differences in environmental light conditions. First, we characterize the selective regime experienced by cichlid opsin genes using a likelihood ratio test format, comparing likelihood models with different constraints on the relative rates of amino acid substitution, across sites. Second, we compare turbid and clear lineages to determine if there is evidence of differences in relative rates of substitution. Third, we present evidence of functional diversification and its relationship to the photic environment among cichlid opsin genes. We report statistical evidence of positive selection in all cichlid opsin genes, except short wavelength-sensitive 1 and short wavelength-sensitive 2b. In all genes predicted to be under positive selection, except short wavelength-sensitive 2a, we find differences in selective pressure between turbid and clear lineages. Potential spectral tuning sites are variable among all cichlid opsin genes; however, patterns of substitution consistent with photic environment-driven evolution of opsin genes are observed only for short wavelength-sensitive 1 opsin genes. This study identifies a number of promising candidate-tuning sites for future study by site-directed mutagenesis. This work also begins to demonstrate the molecular evolutionary dynamics of cichlid visual sensitivity and its relationship to the photic environment.

Relationships among the three major lineages of the Acari (Arthropoda:Arachnida) inferred from small subunit rRNA: paraphyly of the Parasitiformes with respect to the Opilioacariformes and relative rates of nucleotide substitution

We inferred phylogeny among the three major lineages of the Acari (mites) from the small subunit rRNA gene. Our phylogeny indicates that the Opilioacariformes is the sister-group to the Ixodida+Holothyrida, not the Ixodida+Mesostigmata+Holothyrida, as previously thought. Support for this relationship increased when sites with the highest rates of nucleotide substitution, and thus the greatest potential for saturation with nucleotide substitutions, were removed. Indeed, the increase in support (and resolution) was despite a 70% reduction in the number of parsimony-informative sites from 408 to 115. This shows that rather than ‘noisy’ sites having no impact on resolution of deep branches, ‘noisy’ sites have the potential to obscure phylogenetic relationships. The arrangement, Ixodida+Holothyrida+Opilioacariformes, however, may be an artefact of long-branch attraction since relative-rate tests showed that the Mesostigmata have significantly faster rates of nucleotide substitution than other parasitiform mites. Thus, the fast rates of nucleotide substitution of the Mesostigmata might have caused the Mesostigmata to be attracted to the outgroup in our trees. We tested the hypothesis that the high rate of nucleotide substitution in some mites was related to their short generation times. The Acari species that have high nucleotide substitution rates usually have short generation times; these mites also tend to be more active and thus have higher metabolic rates than other mites. Therefore, more than one factor may affect the rate of nucleotide substitution in these mites.

Estimating Tradeoffs Among Health Care System’s Objectives

In this paper, we assess the applicability of the World Health Organization’s (WHO) weighting scheme based on expert’s opinions used in ranking 191 countries’ health care system. Whereas the WHO used a composite index aggregated over the five outputs, we consider the disaggregated measures since the preferences for each health system output may vary across countries. Analyzing the same data, using an expanded data envelopment analysis (DEA) approach, we cluster all 191 countries into groups regarding their relative substitution rates on the five disaggregated health objectives. The different tradeoffs among groups and objectives are determined and compared to the WHO weighting scheme. Our findings demonstrate that three groups emerged that are closely related to an indicator of the countries’ wealth with each group having very specific marginal rates of substitution among outputs and inputs. We conclude that utilizing a uniform weighting system does not capture the priorities different groups of countries pursue.

Relative rates of nucleotide substitution in frogs.

Accurate estimation of relative mutation rates of mitochondrial DNA (mtDNA) and single-copy nuclear DNA (scnDNA) within lineages contributes to a general understanding of molecular evolutionary processes and facilitates making demographic inferences from population genetic data. The rate of divergence at synonymous sites ( K(s)) may be used as a surrogate for mutation rate. Such data are available for few organisms and no amphibians. Relative to mammals and birds, amphibian mtDNA is thought to evolve slowly, and the K(s) ratio of mtDNA to scnDNA would be expected to be low as well. Relative K(s) was estimated from a mitochondrial gene, ND2, and a nuclear gene, c-myc, using both "approximate" and likelihood methods. Three lineages of congeneric frogs were studied and this ratio was found to be approximately 16, the highest of previously reported ratios. No evidence of a low K(s) in the nuclear gene was found: c-myc codon usage was not biased, the K(s) was double the intron divergence rate, and the absolute K(s) was similar to estimates obtained here for other genes from other frog species. A high K(s) in mitochondrial vs. nuclear genes was unexpected in light of previous reports of a slow rate of mtDNA evolution in amphibians. These results highlight the need for further investigation of the effects of life history on mutation rates.

Likelihood analysis of asymmetrical mutation bias gradients in vertebrate mitochondrial genomes.

Protein-coding genes in mitochondrial genomes have varying degrees of asymmetric skew in base frequencies at the third codon position. The variation in skew among genes appears to be caused by varying durations of time that the heavy strand spends in the mutagenic single-strand state during replication (D(ssH)). The primary data used to study skew have been the gene-by-gene base frequencies in individual taxa, which provide little information on exactly what kinds of mutations are responsible for the base frequency skew. To assess the contribution of individual mutation components to the ancestral vertebrate substitution pattern, here we analyze a large data set of complete vertebrate mitochondrial genomes in a phylogeny-based likelihood context. This also allows us to evaluate the change in skew continuously along the mitochondrial genome and to directly estimate relative substitution rates. Our results indicate that different types of mutation respond differently to the D(ssH) gradient. A primary role for hydrolytic deamination of cytosines in creating variance in skew among genes was not supported, but rather linearly increasing rates of mutation from adenine to hypoxanthine with D(ssH) appear to drive regional differences in skew. Substitutions due to hydrolytic deamination of cytosines, although common, appear to quickly saturate, possibly due to stabilization by the mitochondrial DNA single-strand-binding protein. These results should form the basis of more realistic models of DNA and protein evolution in mitochondria.