CodonTest: Modeling Amino Acid Substitution Preferences in Coding Sequences

Codon models of evolution have facilitated the interpretation of selective forces operating on genomes. These models, however, assume a single rate of non-synonymous substitution irrespective of the nature of amino acids being exchanged. Recent developments have shown that models which allow for amino acid pairs to have independent rates of substitution offer improved fit over single rate models. However, these approaches have been limited by the necessity for large alignments in their estimation. An alternative approach is to assume that substitution rates between amino acid pairs can be subdivided into rate classes, dependent on the information content of the alignment. However, given the combinatorially large number of such models, an efficient model search strategy is needed. Here we develop a Genetic Algorithm (GA) method for the estimation of such models. A GA is used to assign amino acid substitution pairs to a series of rate classes, where is estimated from the alignment. Other parameters of the phylogenetic Markov model, including substitution rates, character frequencies and branch lengths are estimated using standard maximum likelihood optimization procedures. We apply the GA to empirical alignments and show improved model fit over existing models of codon evolution. Our results suggest that current models are poor approximations of protein evolution and thus gene and organism specific multi-rate models that incorporate amino acid substitution biases are preferred. We further anticipate that the clustering of amino acid substitution rates into classes will be biologically informative, such that genes with similar functions exhibit similar clustering, and hence this clustering will be useful for the evolutionary fingerprinting of genes

Modeling HIV-1 Drug Resistance as Episodic Directional Selection

The evolution of substitutions conferring drug resistance to HIV-1 is both episodic, occurring when patients are on antiretroviral therapy, and strongly directional, with site-specific resistant residues increasing in frequency over time. While methods exist to detect episodic diversifying selection and continuous directional selection, no evolutionary model combining these two properties has been proposed. We present two models of episodic directional selection (MEDS and EDEPS) which allow the a priori specification of lineages expected to have undergone directional selection. The models infer the sites and target residues that were likely subject to directional selection, using either codon or protein sequences. Compared to its null model of episodic diversifying selection, MEDS provides a superior fit to most sites known to be involved in drug resistance, and neither one test for episodic diversifying selection nor another for constant directional selection are able to detect as many true positives as MEDS and EDEPS while maintaining acceptable levels of false positives. This suggests that episodic directional selection is a better description of the process driving the evolution of drug resistance

HIV-Specific Probabilistic Models of Protein Evolution

Comparative sequence analyses, including such fundamental bioinformatics techniques as similarity searching, sequence alignment and phylogenetic inference, have become a mainstay for researchers studying type 1 Human Immunodeficiency Virus (HIV-1) genome structure and evolution. Implicit in comparative analyses is an underlying model of evolution, and the chosen model can significantly affect the results. In general, evolutionary models describe the probabilities of replacing one amino acid character with another over a period of time. Most widely used evolutionary models for protein sequences have been derived from curated alignments of hundreds of proteins, usually based on mammalian genomes. It is unclear to what extent these empirical models are generalizable to a very different organism, such as HIV-1–the most extensively sequenced organism in existence. We developed a maximum likelihood model fitting procedure to a collection of HIV-1 alignments sampled from different viral genes, and inferred two empirical substitution models, suitable for describing between-and within-host evolution. Our procedure pools the information from multiple sequence alignments, and provided software implementation can be run efficiently in parallel on a computer cluster. We describe how the inferred substitution models can be used to generate scoring matrices suitable for alignment and similarity searches. Our models had a consistently superior fit relative to the best existing models and to parameter-rich data-driven models when benchmarked on independent HIV-1 alignments, demonstrating evolutionary biases in amino-acid substitution that are unique to HIV, and that are not captured by the existing models. The scoring matrices derived from the models showed a marked difference from common amino-acid scoring matrices. The use of an appropriate evolutionary model recovered a known viral transmission history, whereas a poorly chosen model introduced phylogenetic error. We argue that our model derivation procedure is immediately applicable to other organisms with extensive sequence data available, such as Hepatitis C and Influenza A viruses

Adaptive Evolution of the Myo6 Gene in Old World Fruit Bats (Family: Pteropodidae)

PMCID: PMC3631194This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Understanding the molecular determinants driving the immunological specificity of the protective pilus 2a backbone protein of Group B Streptococcus

The pilus 2a backbone protein (BP-2a) is one of the most structurally and functionally characterized components of a potential vaccine formulation against Group B Streptococcus. It is characterized by six main immunologically distinct allelic variants, each inducing variant-specific protection. To investigate the molecular determinants driving the variant immunogenic specificity of BP-2a, in terms of single residue contributions, we generated six monoclonal antibodies against a specific protein variant based on their capability to recognize the polymerized pili structure on the bacterial surface. Three mAbs were also able to induce complement-dependent opsonophagocytosis killing of live GBS and target the same linear epitope present in the structurally defined and immunodominant domain D3 of the protein. Molecular docking between the modelled scFv antibody sequences and the BP-2a crystal structure revealed the potential role at the binding interface of some non-conserved antigen residues. Mutagenesis analysis confirmed the necessity of a perfect balance between charges, size and polarity at the binding interface to obtain specific binding of mAbs to the protein antigen for a neutralizing response

The Nonlinear Meissner Effect in Unconventional Superconductors

We examine the long-wavelength current response in anisotropic superconductors and show how the field-dependence of the Meissner penetration length can be used to detect the structure of the order parameter. Nodes in the excitation gap lead to a nonlinear current-velocity constitutive equation at low temperatures which is distinct for each symmetry class of the order parameter. The effective Meissner penetration length is linear in $H$ and exhibits a characteristic anisotropy for fields in the $ab$-plane that is determined by the positions of the nodes in momentum space. The nonlinear current-velocity relation also leads to an intrinsic magnetic torque for in-plane fields that are not parallel to a nodal or antinodal direction. The torque scales as $H^3$ for $T\rightarrow 0$ and has a characteristic angular dependence. We analyze the effects of thermal excitations, impurity scattering and geometry on the current response of a $d_{x^2-y^2}$ superconductor, and discuss our results in light of recent measurements of the low-temperature penetration length and in-plane magnetization of single-crystals of $YBa_2Cu_3O_{7-\delta}$ and $LuBa_2Cu_3O_{7-\delta}$.Comment: 30 pages, RevTeX file with 16 postscript figures. Submitted to Phys. Rev.

Probing a Complex of Cytochromecand Cardiolipin by Magnetic Circular Dichroism Spectroscopy: Implications for the Initial Events in Apoptosis

Oxidation of cardiolipin (CL) by its complex with cytochrome c (cyt c) plays a crucial role in triggering apoptosis. Through a combination of magnetic circular dichroism spectroscopy and potentiometric titrations, we show that both the ferric and ferrous forms of the heme group of a CL:cyt c complex exist as multiple conformers at a physiologically relevant pH of 7.4. For the ferric state, these conformers are His/Lys- and His/OH–-ligated. The ferrous state is predominantly high-spin and, most likely, His/–. Interconversion of the ferric and ferrous conformers is described by a single midpoint potential of -80 ± 9 mV vs SHE. These results suggest that CL oxidation in mitochondria could occur by the reaction of molecular oxygen with the ferrous CL:cyt c complex in addition to the well-described reaction of peroxides with the ferric form

Understanding the attitudes of the elderly towards enrolment into cancer clinical trials

BACKGROUND: The optimal cancer treatment for an older population is largely unknown because of the low numbers of elderly patients accrued into clinical trials. This project focuses on the attitudes of the elderly about participation in clinical trials to determine if this is one of the barriers to the involvement of this population in clinical trials. METHODS: The first phase of this study was a self-administered questionnaire mailed to 425 elderly persons with cancer, selected from Princess Margaret Hospital oncology clinics. The second phase consisted of individual semi-structured interviews with cancer patients to assess their attitudes towards cancer, its management and enrolment into cancer clinical trials. RESULTS: Ninety-four patients responded to the survey giving a response rate of 22.1%. Three quarters of respondents stated that they would be willing to participate in a clinical trial. The factors that most influenced older patients' willingness to participate in a cancer study were recommendations from a cancer doctor and the chance that the study treatment may help them feel better. Seventeen survey responders participated in interviews. Common themes from these interviews included patient-physician communication, the referral process, and the role of age in cancer care decision-making. CONCLUSION: Most elderly people, who responded to this survey, are willing to consider participation in cancer clinical trials however, elderly patients do not appear to actively seek clinical trials and few were informed of the availability of clinical trials. Physician barriers and availability of appropriate clinical trials may play a bigger role in preventing accrual of elderly cancer patients into trials

Evolutionary Dynamics and Emergence of Panzootic H5N1 Influenza Viruses

The highly pathogenic avian influenza (HPAI) H5N1 virus lineage has undergone extensive genetic reassortment with viruses from different sources to produce numerous H5N1 genotypes, and also developed into multiple genetically distinct sublineages in China. From there, the virus has spread to over 60 countries. The ecological success of this virus in diverse species of both poultry and wild birds with frequent introduction to humans suggests that it is a likely source of the next human pandemic. Therefore, the evolutionary and ecological characteristics of its emergence from wild birds into poultry are of considerable interest. Here, we apply the latest analytical techniques to infer the early evolutionary dynamics of H5N1 virus in the population from which it emerged (wild birds and domestic poultry). By estimating the time of most recent common ancestors of each gene segment, we show that the H5N1 prototype virus was likely introduced from wild birds into poultry as a non-reassortant low pathogenic avian influenza H5N1 virus and was not generated by reassortment in poultry. In contrast, more recent H5N1 genotypes were generated locally in aquatic poultry after the prototype virus (A/goose/Guangdong/1/96) introduction occurred, i.e., they were not a result of additional emergence from wild birds. We show that the H5N1 virus was introduced into Indonesia and Vietnam 3–6 months prior to detection of the first outbreaks in those countries. Population dynamics analyses revealed a rapid increase in the genetic diversity of A/goose/Guangdong/1/96 lineage viruses from mid-1999 to early 2000. Our results suggest that the transmission of reassortant viruses through the mixed poultry population in farms and markets in China has selected HPAI H5N1 viruses that are well adapted to multiple hosts and reduced the interspecies transmission barrier of those viruses

Phylodynamic Reconstruction Reveals Norovirus GII.4 Epidemic Expansions and their Molecular Determinants

Noroviruses are the most common cause of viral gastroenteritis. An increase in the number of globally reported norovirus outbreaks was seen the past decade, especially for outbreaks caused by successive genogroup II genotype 4 (GII.4) variants. Whether this observed increase was due to an upswing in the number of infections, or to a surveillance artifact caused by heightened awareness and concomitant improved reporting, remained unclear. Therefore, we set out to study the population structure and changes thereof of GII.4 strains detected through systematic outbreak surveillance since the early 1990s. We collected 1383 partial polymerase and 194 full capsid GII.4 sequences. A Bayesian MCMC coalescent analysis revealed an increase in the number of GII.4 infections during the last decade. The GII.4 strains included in our analyses evolved at a rate of 4.3–9.0×10−3 mutations per site per year, and share a most recent common ancestor in the early 1980s. Determinants of adaptation in the capsid protein were studied using different maximum likelihood approaches to identify sites subject to diversifying or directional selection and sites that co-evolved. While a number of the computationally determined adaptively evolving sites were on the surface of the capsid and possible subject to immune selection, we also detected sites that were subject to constrained or compensatory evolution due to secondary RNA structures, relevant in virus-replication. We highlight codons that may prove useful in identifying emerging novel variants, and, using these, indicate that the novel 2008 variant is more likely to cause a future epidemic than the 2007 variant. While norovirus infections are generally mild and self-limiting, more severe outcomes of infection frequently occur in elderly and immunocompromized people, and no treatment is available. The observed pattern of continually emerging novel variants of GII.4, causing elevated numbers of infections, is therefore a cause for concern