Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level

Detecting positive Darwinian selection at the DNA sequence level has been a subject of considerable interest. However, positive selection is difficult to detect because it often operates episodically on a few amino acid sites, and the signal may be masked by negative selection. Several methods have been developed to test positive selection that acts on given branches (branch methods) or on a subset of sites (site methods). Recently, Yang, Z., and R. Nielsen (2002. Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. Mol. Biol. Evol. 19:908-917) developed likelihood ratio tests (LRTs) based on branch-site models to detect positive selection that affects a small number of sites along prespecified lineages. However, computer simulations suggested that the tests were sensitive to the model assumptions and were unable to distinguish between relaxation of selective constraint and positive selection (Zhang, J. 2004. Frequent false detection of positive selection by the likelihood method with branch-site models. Mol. Biol. Evol. 21:1332-1339). Here, we describe a modified branch-site model and use it to construct two LRTs, called branch-site tests 1 and 2. We applied the new tests to reanalyze several real data sets and used computer simulation to examine the performance of the two tests by examining their false-positive rate, power, and robustness. We found that test 1 was unable to distinguish relaxed constraint from positive selection affecting the lineages of interest, while test 2 had acceptable false-positive rates and appeared robust against violations of model assumptions. As test 2 is a direct test of positive selection on the lineages of interest, it is referred to as the branch-site test of positive selection and is recommended for use in real data analysis. The test appeared conservative overall, but exhibited better power in detecting positive selection than the branch-based test. Bayes empirical Bayes identification of amino acid sites under positive selection along the foreground branches was found to be reliable, but lacked power.     

Positive selection, relaxation, and acceleration in the evolution of the human and chimp genome

For years evolutionary biologists have been interested in searching for the genetic bases underlying humanness. Recent efforts at a large or a complete genomic scale have been conducted to search for positively selected genes in human and in chimp. However, recently developed methods allowing for a more sensitive and controlled approach in the detection of positive selection can be employed. Here, using 13,198 genes, we have deduced the sets of genes involved in rate acceleration, positive selection, and relaxation of selective constraints in human, in chimp, and in their ancestral lineage since the divergence from murids. Significant deviations from the strict molecular clock were observed in 469 human and in 651 chimp genes. The more stringent branch-site test of positive selection detected 108 human and 577 chimp positively selected genes. An important proportion of the positively selected genes did not show a significant acceleration in rates, and similarly, many of the accelerated genes did not show significant signals of positive selection. Functional differentiation of genes under rate acceleration, positive selection, and relaxation was not statistically significant between human and chimp with the exception of terms related to G-protein coupled receptors and sensory perception. Both of these were over-represented under relaxation in human in relation to chimp. Comparing differences between derived and ancestral lineages, a more conspicuous change in trends seems to have favored positive selection in the human lineage. Since most of the positively selected genes are different under the same functional categories between these species, we suggest that the individual roles of the alternative positively selected genes may be an important factor underlying biological differences between these species.     

Evidence for positive selection on the leptin gene in Cetacea and Pinnipedia

The leptin gene has received intensive attention and scientific investigation for its importance in energy homeostasis and reproductive regulation in mammals. Furthermore, study of the leptin gene is of crucial importance for public health, particularly for its role in obesity, as well as for other numerous physiological roles that it plays in mammals. In the present work, we report the identification of novel leptin genes in 4 species of Cetacea, and a comparison with 55 publicly available leptin sequences from mammalian genome assemblies and previous studies. Our study provides evidence for positive selection in the suborder Odontoceti (toothed whales) of the Cetacea and the family Phocidae (earless seals) of the Pinnipedia. We also detected positive selection in several leptin gene residues in these two lineages. To test whether leptin and its receptor evolved in a coordinated manner, we analyzed 24 leptin receptor gene (LPR) sequences from available mammalian genome assemblies and other published data. Unlike the case of leptin, our analyses did not find evidence of positive selection for LPR across the Cetacea and Pinnipedia lineages. In line with this, positively selected sites identified in the leptin genes of these two lineages were located outside of leptin receptor binding sites, which at least partially explains why co-evolution of leptin and its receptor was not observed in the present study. Our study provides interesting insights into current understanding of the evolution of mammalian leptin genes in response to selective pressures from life in an aquatic environment, and leads to a hypothesis that new tissue specificity or novel physiologic functions of leptin genes may have arisen in both odontocetes and phocids. Additional data from other species encompassing varying life histories and functional tests of the adaptive role of the amino acid changes identified in this study will help determine the factors that promote the adaptive evolution of the leptin genes in marine mammals.     

Statistical properties of the branch-site test of positive selection

The branch-site test is a likelihood ratio test to detect positive selection along prespecified lineages on a phylogeny that affects only a subset of codons in a protein-coding gene, with positive selection indicated by accelerated nonsynonymous substitutions (with ω = d(N)/d(S) > 1). This test may have more power than earlier methods, which average nucleotide substitution rates over sites in the protein and/or over branches on the tree. However, a few recent studies questioned the statistical basis of the test and claimed that the test generated too many false positives. In this paper, we examine the null distribution of the test and conduct a computer simulation to examine the false-positive rate and the power of the test. The results suggest that the asymptotic theory is reliable for typical data sets, and indeed in our simulations, the large-sample null distribution was reliable with as few as 20-50 codons in the alignment. We examined the impact of sequence length, the strength of positive selection, and the proportion of sites under positive selection on the power of the branch-site test. We found that the test was far more powerful in detecting episodic positive selection than branch-based tests, which average substitution rates over all codons in the gene and thus miss the signal when most codons are under strong selective constraint. Recent claims of statistical problems with the branch-site test are due to misinterpretations of simulation results. Our results, as well as previous simulation studies that have demonstrated the robustness of the test, suggest that the branch-site test may be a useful tool for detecting episodic positive selection and for generating biological hypotheses for mutation studies and functional analyses. The test is sensitive to sequence and alignment errors and caution should be exercised concerning its use when data quality is in doubt.     

Positive Darwinian selection at the pantophysin (Pan I) locus in marine gadid fishes

Maximum-likelihood models of codon substitution were used to test for positive Darwinian selection at the vesicle protein pantophysin in two allelic lineages segregating in the Atlantic cod Gadus morhua and in 18 related species of marine gadid fishes. Positive selection was detected in the two intravesicular loops of the integral membrane protein but not in four membrane-spanning regions or the 3' cytoplasmic tail. The proportion of positively selected sites (24.9%) and the mean nonsynonymous/synonymous rate ratio (omega = d(N)/d(S) = 5.35) were both greater in the first intravesicular (IV1) domain compared with the second intravesicular (IV2) domain (11.0% positively selected sites with mean omega = 3.76). Likelihood ratio tests comparing models that assume identical omega ratios along all branches of the phylogeny to those that allow omega ratios to vary among lineages were not significant for either the IV1 or IV2 domains, indicating that the selective pressures favoring amino acid replacements have operated consistently in both regions during the diversification of the group. Positive selection was observed in the IV1 domain in both G. morhua allelic lineages, and, although three of the four codons that differ between alleles were targets of positive selection in the broader group, no similar polymorphisms were detected in other taxa. The two G. morhua Pan I alleles appeared to have evolved before the speciation event separating it from its sister taxon, Theragra chalcogramma, and on the basis of a standard mtDNA clock are estimated to be at least 2 Myr old. Although the function of pantophysin remains unknown, the strong signal of positive selection at specific sites in the IV1 and IV2 domains may help clarify its role in cellular trafficking pathways.     

False-positive results obtained from the branch-site test of positive selection

Natural selection operating at the amino acid sequence level can be detected by comparing the rates of synonymous (r(S)) and nonsynonymous (r(N)) nucleotide substitutions, where r(N)/r(S) (omega) > 1 and omega < 1 suggest positive and negative selection, respectively. The branch-site test has been developed for detecting positive selection operating at a group of amino acid sites for a pre-specified (foreground) branch of a phylogenetic tree by taking into account the heterogeneity of omega among sites and branches. Here the performance of the branch-site test was examined by computer simulation, with special reference to the false-positive rate when the divergence of the sequences analyzed was small. The false-positive rate was found to inflate when the assumptions made on the omega values for the foreground and other (background) branches in the branch-site test were violated. In addition, under a similar condition, false-positive results were often obtained even when Bonferroni correction was conducted and the false-discovery rate was controlled in a large-scale analysis. False-positive results were also obtained even when the number of nonsynonymous substitutions for the foreground branch was smaller than the minimum value required for detecting positive selection. The existence of a codon site with a possibility of occurrence of multiple nonsynonymous substitutions for the foreground branch often caused the branch-site test to falsely identify positive selection. In the re-analysis of orthologous trios of protein-coding genes from humans, chimpanzees, and macaques, most of the genes previously identified to be positively selected for the human or chimpanzee branch by the branch-site test contained such a codon site, suggesting a possibility that a significant fraction of these genes are false-positives.     

The Avian XPR1 Gammaretrovirus Receptor Is under Positive Selection and Is Disabled in Bird Species in Contact with Virus-Infected Wild Mice

Xenotropic mouse leukemia viruses (X-MLVs) are broadly infectious for mammals except most of the classical strains of laboratory mice. These gammaretroviruses rely on the XPR1 receptor for entry, and the unique resistance of laboratory mice is due to two mutations in different putative XPR1 extracellular loops. Cells from avian species differ in susceptibility to X-MLVs, and 2 replacement mutations in the virus-resistant chicken XPR1 (K496Q and Q579E) distinguish it from the more permissive duck and quail receptors. These substitutions align with the two mutations that disable the laboratory mouse XPR1. Mutagenesis of the chicken and duck genes confirms that residues at both sites are critical for virus entry. Among 32 avian species, the 2 disabling XPR1 mutations are found together only in the chicken, an omnivorous, ground-dwelling fowl that was domesticated in India and/or Southeast Asia, which is also where X-MLV-infected house mice evolved. The receptor-disabling mutations are also present separately in 5 additional fowl and raptor species, all of which are native to areas of Asia populated by the virus-infected subspecies Mus musculus castaneus. Phylogenetic analysis showed that the avian XPR1 gene is under positive selection at sites implicated in receptor function, suggesting a defensive role for XPR1 in the avian lineage. Contact between bird species and virus-infected mice may thus have favored selection of mouse virus-resistant receptor orthologs in the birds, and our data suggest that similar receptor-disabling mutations were fixed in mammalian and avian species exposed to similar virus challenges.     

Positive diversifying selection in avian Mx genes

Mx proteins are interferon-induced GTPases that confer antiviral activities against RNA viruses. We analysed the molecular evolution of the Mx gene in birds using data on interspecific divergence in anseriform and galliform birds, and on intraspecific diversity in commercial chicken lines, local Chinese chicken breeds as well as in the mallard. The overall ratio of non-synonymous to synonymous substitution was unusually high, 0.80, indicating relaxed constraint or positive selection. Evidence for the latter was provided by that a total of 11–18 codons were found to have evolved under positive selection. The great majority of these codons are located in a region unique to birds at the N-terminal end of the Mx protein. We found an excess of non-synonymous polymorphisms relative to synonymous variants in all comparisons. This, together with positive Tajima’s D values in the local Chinese chicken breeds and in the mallard suggests that balancing selection is acting in avian Mx genes. As such, Mx mimics the major histocompatibility complex system, indicating that heterozygous individuals are better off withstanding pathogen attack. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. S.B. and L.Q. contributed equally to this work.     

Episodic positive selection during the evolution of naphthalene dioxygenase to nitroarene dioxygenase

Using different maximum-likelihood models of adaptive evolution, signatures of natural selective pressure, operating across the naphthalene family of dioxygenases, were examined. A lineage- and branch-site specific combined analysis revealed that purifying selection pressure dominated the evolutionary history of the enzyme family. Specifically, episodic positive Darwinian selection pressure, affecting only a few sites in a subset of lineages, was found to be responsible for the evolution of nitroarene dioxygenases (NArDO) from naphthalene dioxygenase (NDO). Site-specific analysis confirmed the absence of diversifying selection pressure at any particular site. Different sets of positively selected residues, obtained from branch-site specific analysis, were detected for the evolution of each NArDO. They were mainly located around the active site, the catalytic pocket and their adjacent regions, when mapped onto the 3D structure of the α-subunit of NDO. The present analysis enriches the current understanding of adaptive evolution and also broadens the scope for rational alteration of substrate specificity of enzyme by directed evolution.     

Multiple hypothesis testing to detect lineages under positive selection that affects only a few sites

Detection of positive Darwinian selection has become ever more important with the rapid growth of genomic data sets. Recent branch-site models of codon substitution account for variation of selective pressure over branches on the tree and across sites in the sequence and provide a means to detect short episodes of molecular adaptation affecting just a few sites. In likelihood ratio tests based on such models, the branches to be tested for positive selection have to be specified a priori. In the absence of a biological hypothesis to designate so-called foreground branches, one may test many branches, but a correction for multiple testing becomes necessary. In this paper, we employ computer simulation to evaluate the performance of 6 multiple test correction procedures when the branch-site models are used to test every branch on the phylogeny for positive selection. Four of the methods control the familywise error rates (FWERs), whereas the other 2 control the false discovery rate (FDR). We found that all correction procedures achieved acceptable FWER except for extremely divergent sequences and serious model violations, when the test may become unreliable. The power of the test to detect positive selection is influenced by the strength of selection and the sequence divergence, with the highest power observed at intermediate divergences. The 4 correction procedures that control the FWER had similar power. We recommend Rom's procedure for its slightly higher power, but the simple Bonferroni correction is useable as well. The 2 correction procedures that control the FDR had slightly more power and also higher FWER. We demonstrate the multiple test procedures by analyzing gene sequences from the extracellular domain of the cluster of differentiation 2 (CD2) gene from 10 mammalian species. Both our simulation and real data analysis suggest that the multiple test procedures are useful when multiple branches have to be tested on the same data set.     

裸子植物中光敏色素PHY-PAS1结构域的适应性进化

光敏色素是一类红光／远红光受体,在植物种子萌发到成熟的整个生长发育过程中均起重要的调节作用。光敏色素PHY-PAS1结构域存在于光敏色素基因家族的所有成员中,对调节发色团的光谱特性和光信号转导非常关键。光敏色素基因家族通过基因重复产生,而基因重复可能与物种形成有关。PHYP基因是裸子植物光敏色素基因家族发生第1次重复后产生的,并且以单拷贝形式存在。为了研究不同裸子植物PHYP基因编码蛋白的PHY-PAS1结构域在进化过程中是否受到相同的选择压力以及是否发生了适应性进化,该研究利用分支模型、位点模型以及分支．位点模型对裸子植物31条PHYP基因序列编码蛋白的PHY-PAS1结构域所受到的选择压力进行了分析。结果表明,在由PHY-PAS1结构域序列构建的系统树中,多数分支处于强烈的负选择压力下（ω〈1）：有14个分支处于正选择压力下（ω〉1）,其中13个分支发生在属内种间;与之相比,在较为古老的谱系中相对缺少这种正选择压力。     

Evidence for positive selection on the floral scent gene isoeugenol-O-methyltransferase

Isoeugenol-O-methyltransferase (IEMT) is an enzyme involved in the production of the floral volatile compounds methyl eugenol and methyl isoeugenol in Clarkia breweri (Onagraceae). IEMT likely evolved by gene duplication from caffeic acid-O-methyltransferase followed by amino acid divergence, leading to the acquisition of its novel function. To investigate the selective context under which IEMT evolved, maximum likelihood methods that estimate variable d(N)/d(S) ratios among lineages, among sites, and among a combination of both lineages and sites were utilized. Statistically significant support was obtained for a hypothesis of positive selection driving the evolution of IEMT since its origin. Subsequent Bayesian analyses identified several sites in IEMT that have experienced positive selection. Most of these positions are in the active site of IEMT and have been shown by site-directed mutagenesis to have large effects on substrate specificity. Although the selective agent is unknown, the adaptive evolution of this gene may have resulted in increased effectiveness of pollinator attraction or herbivore repellence.     

裸子植物中光敏色素PHY-PAS1 结构域的适应性进化

光敏色素是一类红光/远红光受体, 在植物种子萌发到成熟的整个生长发育过程中均起重要的调节作用。光敏色素PHYPAS1 结构域存在于光敏色素基因家族的所有成员中, 对调节发色团的光谱特性和光信号转导非常关键。光敏色素基因家族通过基因重复产生, 而基因重复可能与物种形成有关。PHYP基因是裸子植物光敏色素基因家族发生第1次重复后产生的, 并且以单拷贝形式存在。为了研究不同裸子植物PHYP基因编码蛋白的PHY-PAS1结构域在进化过程中是否受到相同的选择压力以及是否发生了适应性进化, 该研究利用分支模型、位点模型以及分支-位点模型对裸子植物31条PHYP基因序列编码蛋白的PHY-PAS1结构域所受到的选择压力进行了分析。结果表明, 在由PHY-PAS1结构域序列构建的系统树中, 多数分支处于强烈的负选择压力下 (w<1); 有14个分支处于正选择压力下 (w>1), 其中13个分支发生在属内种间; 与之相比, 在较为古老的谱系中相对缺少这种正选择压力。     

Frequent false detection of positive selection by the likelihood method with branch-site models

Positive Darwinian selection promotes fixations of advantageous mutations during gene evolution and is probably responsible for most adaptations. Detecting positive selection at the DNA sequence level is of substantial interest because such information provides significant insights into possible functional alterations during gene evolution as well as important nucleotide substitutions involved in adaptation. Efficient detection of positive selection, however, has been difficult because selection often operates on only a few sites in a short period of evolutionary time. A likelihood-based method with branch-site models was recently introduced to overcome such difficulties. Here I examine the accuracy of the method using computer simulation. I find that the method detects positive selection in 20%-70% of cases when the DNA sequences are generated by computer simulation under no positive selection. Although the frequency of such false detection varies depending on, among other things, the tree topology, branch length, and selection scheme, the branch-site likelihood method generally gives misleading results. Thus, detection of positive selection by this method alone is unreliable. This unreliability may have resulted from its over-sensitivity to violations of assumptions made in the method, such as certain distributions of selective strength among sites and equal transition/transversion ratios for synonymous and nonsynonymous substitutions.     

Heterogeneous Selective Pressure Acting on Influenza B Victoria- and Yamagata-Like Hemagglutinins

As a consequence of immune pressure, influenza virus hemagglutinin presents some of its amino acids under positive selection. Several authors have reported the existence of influenza A hemagglutinin codons under positive selective pressure (PSP). In this framework, the present work objectives were to demonstrate the presence of PSP and evaluate its effects on Victoria- and Yamagata-like influenza B viruses. Methodology adopted consisted in estimating the acceptance rate of nonsynonymous substitutions (ω = d_N/d_S) that describe the strength of selective pressure and identifying codons that may be positively selected, applying a set of continuous-time Markov chain codon-substitution models. Two groups of HA1 sequences (140 from Yamagata and 60 from Victoria lineage) were used. All the model maximum-likelihood estimates were obtained using codeml software application (PAML 3.15). The hypothesis of no existence of sites under PSP was rejected for both lineages (p < 0.001), using likelihood ratio tests. These results demonstrate the presence of positive selection acting on hemagglutinin of both Yamagata- and Victoria-like influenza B viruses. Several different sites were identified to be under PSP on Yamagata and Victoria hemagglutinins. Sites found with a posterior probability > 0.95 were codons 197 and 199 in both lineages, codon 75 in the Yamagata lineage, and codon 129 in the Victoria lineage. The detected amino acids are located at or near antigenic sites in influenza A virus H3 hemagglutinin. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A random effects branch-site model for detecting episodic diversifying selection

Adaptive evolution frequently occurs in episodic bursts, localized to a few sites in a gene, and to a small number of lineages in a phylogenetic tree. A popular class of "branch-site" evolutionary models provides a statistical framework to search for evidence of such episodic selection. For computational tractability, current branch-site models unrealistically assume that all branches in the tree can be partitioned a priori into two rigid classes--"foreground" branches that are allowed to undergo diversifying selective bursts and "background" branches that are negatively selected or neutral. We demonstrate that this assumption leads to unacceptably high rates of false positives or false negatives when the evolutionary process along background branches strongly deviates from modeling assumptions. To address this problem, we extend Felsenstein's pruning algorithm to allow efficient likelihood computations for models in which variation over branches (and not just sites) is described in the random effects likelihood framework. This enables us to model the process at every branch-site combination as a mixture of three Markov substitution models--our model treats the selective class of every branch at a particular site as an unobserved state that is chosen independently of that at any other branch. When benchmarked on a previously published set of simulated sequences, our method consistently matched or outperformed existing branch-site tests in terms of power and error rates. Using three empirical data sets, previously analyzed for episodic selection, we discuss how modeling assumptions can influence inference in practical situations.     

Varying signals of the effects of natural selection during teleost growth hormone gene evolution.

The growth hormone (GH) gene of teleost fish exhibits a higher degree of variability compared with other vertebrate groups. However, the different selective constraints at the sequence level are not well understood. In this study, maximum-likelihood (ML) models of codon substitutions were used to investigate Darwinian adaptive evolution of the GH gene in teleost fishes. Complete GH gene sequences of 54 fish species were classified into 4 orders, and the variable nature of GH was examined by determining the dN and dS rate variation and the rates of molecular evolution for each teleost order. The results indicate that although the overall evolution rate for teleost GH is high ((1.15 +/- 0.01) x 10(-9) substitutions/(aa site x y)) compared with the "slow phases" in mammals ((0.21 to 0.28 +/- 0.05) x 10(-9)), the vital structure of this gene has been retained. While the majority of the amino acid changes appear to be due to relaxation of purifying selection, some positively selected sites were detected in regions with no specifically identified role in protein function. The positively selected regions observed in salmoniformes lineage suggests a possible role for positive selection driving functional divergence in paralogous forms of the GH gene after whole-genome duplication in this lineage.     

Molecular evolution of genes in avian genomes

Background

Obtaining a draft genome sequence of the zebra finch (Taeniopygia guttata), the second bird genome to be sequenced, provides the necessary resource for whole-genome comparative analysis of gene sequence evolution in a non-mammalian vertebrate lineage. To analyze basic molecular evolutionary processes during avian evolution, and to contrast these with the situation in mammals, we aligned the protein-coding sequences of 8,384 1:1 orthologs of chicken, zebra finch, a lizard and three mammalian species.

Results

We found clear differences in the substitution rate at fourfold degenerate sites, being lowest in the ancestral bird lineage, intermediate in the chicken lineage and highest in the zebra finch lineage, possibly reflecting differences in generation time. We identified positively selected and/or rapidly evolving genes in avian lineages and found an over-representation of several functional classes, including anion transporter activity, calcium ion binding, cell adhesion and microtubule cytoskeleton.

Conclusions

Focusing specifically on genes of neurological interest and genes differentially expressed in the unique vocal control nuclei of the songbird brain, we find a number of positively selected genes, including synaptic receptors. We found no evidence that selection for beneficial alleles is more efficient in regions of high recombination; in fact, there was a weak yet significant negative correlation between ω and recombination rate, which is in the direction predicted by the Hill-Robertson effect if slightly deleterious mutations contribute to protein evolution. These findings set the stage for studies of functional genetics of avian genes.     

Maximum‐likelihood approaches reveal signatures of positive selection in BMP15 and GDF9 genes modulating ovarian function in mammalian female fertility

Bone morphogenetic proteins (BMPs) and the growth factors (GDFs) play an important role in ovarian folliculogenesis and essential regulator of processes of numerous granulosa cells. BMP15 gene variations linked to various ovarian phenotypic consequences subject to the species, from infertility to improved prolificacy in sheep, primary ovarian insufficiency in women or associated with minor subfertility in mouse. To study the evolving role of BMP15 and GDF9, a phylogenetic analysis was performed. To find out the candidate gene associated with prolificacy in mammals, the nucleotide sequence of BMP15 and GDF9 genes was recognized under positive selection in various mammalian species. Maximum‐likelihood approaches used on BMP15 and GDF9 genes exhibited a robust divergence and a prompted evolution as compared to other TGFβ family members. Furthermore, among 32 mammalian species, we identified positive selection signals in the hominidae clade resulting to 132D, 147E, 163Y, 191W, and 236P codon sites of BMP15 and 162F, 188K, 206R, 240A, 244L, 246H, 248S, 251D, 253L, 254F and other codon sites of GDF9. The positively selected amino acid sites such as Alanine, Lucien, Arginine, and lysine are important for signaling. In conclusion, this study evidences that GDF9 and BMP15 genes have rapid evolution than other TGFß family members and was subjected to positive selection in the mammalian clade. Selected sites under the positive selection are of remarkable significance for the particular functioning of the protein and consequently for female fertility.     

The Cytoplasmic Location of Chicken Mx Is Not the Determining Factor for Its Lack of Antiviral Activity

Background

Chicken Mx belongs to the Mx family of interferon-induced dynamin-like GTPases, which in some species possess potent antiviral properties. Conflicting data exist for the antiviral capability of chicken Mx. Reports of anti-influenza activity of alleles encoding an Asn631 polymorphism have not been supported by subsequent studies. The normal cytoplasmic localisation of chicken Mx may influence its antiviral capacity. Here we report further studies to determine the antiviral potential of chicken Mx against Newcastle disease virus (NDV), an economically important cytoplasmic RNA virus of chickens, and Thogoto virus, an orthomyxovirus known to be exquisitely sensitive to the cytoplasmic MxA protein from humans. We also report the consequences of re-locating chicken Mx to the nucleus.

Methodology/Principal Findings

Chicken Mx was tested in virus infection assays using NDV. Neither the Asn631 nor Ser631 Mx alleles (when transfected into 293T cells) showed inhibition of virus-directed gene expression when the cells were subsequently infected with NDV. Human MxA however did show significant inhibition of NDV-directed gene expression. Chicken Mx failed to inhibit a Thogoto virus (THOV) minireplicon system in which the cytoplasmic human MxA protein showed potent and specific inhibition. Relocalisation of chicken Mx to the nucleus was achieved by inserting the Simian Virus 40 large T antigen nuclear localisation sequence (SV40 NLS) at the N-terminus of chicken Mx. Nuclear re-localised chicken Mx did not inhibit influenza (A/PR/8/34) gene expression during virus infection in cell culture or influenza polymerase activity in A/PR/8/34 or A/Turkey/50-92/91 minireplicon systems.

Conclusions/Significance

The chicken Mx protein (Asn631) lacks inhibitory effects against THOV and NDV, and is unable to suppress influenza replication when artificially re-localised to the cell nucleus. Thus, the natural cytoplasmic localisation of the chicken Mx protein does not account for its lack of antiviral activity.