Phylogenetic analysis of genes of the influenza virus. Relationship between adaptability and neutrality

Detailed phylogenetic analysis of the gene family of hemagglutinin H3 of influenza A-type virus was fulfilled, taking into account the domain structure of protein and positions of antigen determinants. The densities of distribution of fixed synonimic replacements between domains HA1 and HA2 were shown to be actually equal (rho (HA1) = rho (HA2], and those of nonsynonimic ones to be unequal: their ratios were rho (HA1): rho (HA2) = 2.8 for nonepidemic branches, and rho (HA1): rho (HA2) = 7.7 for epidemic ones. For the positions of antigen determinants (agd) these densities differ still stronger from HA2: rho (agd): rho (HA2) = 10 for nonepidemic branches, and rho (agd): rho (HA2) = 36 for epidemic ones. In total, the rate of fixation of nonsynonimic replacements per position of antigen determinants for epidemic branches is 32 times higher than for nonepidemic ones. The absolute value of this estimation is K(ns)d = (9.1 +/- 0.7).10(-3) of nonsynonimic replacements per nonsynonimic position per year and seems to be twice as much as the maximum rate of neutral fixations Ks = (4.28 +/- 0.68).10(-3). Therefore, the epidemic reproduction of influenza virus is highly adaptive, exactly being focused on positions of antigen determinants. The evolution of influenza virus is stochastic process, both with the neutral and adaptive fixations.     

SYMMETREE: whole-tree analysis of differential diversification rates

SymmeTREE implements several tests of differential diversification rates that exploit information on the topological distribution of species diversity throughout entire trees to address two general questions: (1) Has a given tree experienced significant variation in diversification rates among its branches? and (2) If so, along which branches have significant shifts in diversification rate occurred? These explicitly model-based methods are robust to uncertainty in estimates of branch length/duration and can accommodate incompletely resolved trees and other forms of phylogenetic uncertainty. AVAILABILITY: http://www.phylodiversity.net/bmoore/software.html CONTACT: brian.moore@yale.edu.     

An experimental study of colony-founding in pine saplings by queens of the arboreal ant, Crematogaster ashmeadi

Summary: Newly mated queens of the arboreal ant Crematogaster ashmeadi initiate colonies in old beetle galleries in the dead branches of longleaf pine trees. In a study by Hahn (1996), a number of tree characteristics were correlated with the number of newly-mated queens in those trees, with branch length the best indicator of queen presence. Three of these characteristics, tree height, dead branch length, and the number of dead branches were tested in an experiment to see which, if any, the queens were using to choose a tree. Both tree height and the number of dead branches significantly influenced queen choice: shorter trees (4-5 m) had more queens than tall ones (7-9 m), and trees with 8 branches had significantly more queens than trees with 2 branches. Branch length had no effect on the number of queens. These findings suggest that newly mated Crematogaster ashmeadi queens search for founding trees on the basis of the height of a sapling and its number of dead branches. Modes of searching are discussed.     

Positive selection on the H3 hemagglutinin gene of human influenza virus A.

The hemagglutinin (HA) gene of influenza viruses encodes the major surface antigen against which neutralizing antibodies are produced during infection or vaccination. We examined temporal variation in the HA1 domain of HA genes of human influenza A (H3N2) viruses in order to identify positively selected codons. Positive selection is defined for our purposes as a significant excess of nonsilent over silent nucleotide substitutions. If past mutations at positively selected codons conferred a selective advantage on the virus, then additional changes at these positions may predict which emerging strains will predominate and cause epidemics. We previously reported that a 38% excess of mutations occurred on the tip or terminal branches of the phylogenetic tree of 254 HA genes of influenza A (H3N2) viruses. Possible explanations for this excess include processes other than viral evolution during replication in human hosts. Of particular concern are mutations that occur during adaptation of viruses for growth in embryonated chicken eggs in the laboratory. Because the present study includes 357 HA sequences (a 40% increase), we were able to separately analyze those mutations assigned to internal branches. This allowed us to determine whether mutations on terminal and internal branches exhibit different patterns of selection at the level of individual codons. Additional improvements over our previous analysis include correction for a skew in the distribution of amino acid replacements across codons and analysis of a population of phylogenetic trees rather than a single tree. The latter improvement allowed us to ascertain whether minor variation in tree structure had a significant effect on our estimate of the codons under positive selection. This method also estimates that 75.6% of the nonsilent mutations are deleterious and have been removed by selection prior to sampling. Using the larger data set and the modified methods, we confirmed a large (40%) excess of changes on the terminal branches. We also found an excess of changes on branches leading to egg-grown isolates. Furthermore, 9 of the 18 amino acid codons, identified as being under positive selection to change when we used only mutations assigned to internal branches, were not under positive selection on the terminal branches. Thus, although there is overlap between the selected codons on terminal and internal branches, the codons under positive selection on the terminal branches differ from those on the internal branches. We also observed that there is an excess of positively selected codons associated with the receptor-binding site and with the antibody-combining sites. This association may explain why the positively selected codons are restricted in their distribution along the sequence. Our results suggest that future studies of positive selection should focus on changes assigned to the internal branches, as certain of these changes may have predictive value for identifying future successful epidemic variants.     

Modeling egg distribution of Tinocallis kahawaluokalani (Kirkaldy) (Hemiptera: Aphididae) on Lagerstroemia indica L. (Lythraceae)

This is the first record of winter eggs of the holocyclic monoeceous crapemyrtle aphid Tinocallis kahawaluokalani (Kirkaldy) on Lagerstroemia indica L., in Brazil. The shiny black eggs were observed since early autumn, laid on small folds and crevices of the branches. In order to evaluate and model the egg abundance and distribution, four branches from the cardinal points of 10 plants of two age groups, 5-10 and 20-30 years-old, were collected randomly and cut in eight segments of 10 cm and the number of eggs was registered, in the winter 2001. The eggs were laid mainly on the middle portion of the branch, from 40 cm to 60 cm from the apex on the older trees (54%) and on 30 cm to 60 cm on younger ones (58%). The data fit in a longitudinal regression model that expresses the tendency of the egg distribution on the branches. The number of eggs was greater on the 20-30 year-old plants (61%) than on younger ones (39%). The average number (+/- CI 95%) of eggs was 70.5 +/- 9.3 and 47.4 +/- 12.5, respectively, for the older and younger tree groups. There was no correlation between egg distribution and the cardinal positions of the branches.     

应用PV技术对7种针阔叶幼树抗旱性的研究

为进一步探讨PV(Pressure-Volume)技术在树木抗旱性研究中的应用前景,比较成、幼龄树木的抗旱性差别,在对樟子松等7种成龄树木抗旱性生理指标研究的基础上,对其幼树的抗旱性进行了研究。结果表明,对成、幼龄树木一年生小枝用PV曲线所测得的诸水分状况参数随年生长发育进程的的变化规律是一致的;用嫩枝生长初期和枝条完全木质化时期两个阶段的水分参数比较树木的抗旱性、具有较大的可靠性和实用性;幼树在嫩枝生长初期的抗旱性较成树弱,当新枝完全木质化之后,二者的抗旱性基本相近。     

Long-branch attraction and the rDNA model of early eukaryotic evolution

Phylogenetic analyses of ribosomal RNA genes have become widely accepted as a framework for understanding broad-scale eukaryotic evolution. Nevertheless, conflicts exist between the phylogenetic placement of certain taxa in rDNA trees and their expected position based on fossils, cytology, or protein-encoding gene sequences. For example, pelobiont amoebae appear to be an ancient group based on cytologic features, but they are not among the early eukaryotic brances in rDNA analyses. In this report, the derived position of pelobionts in rDNA trees is shown to be unreliable and likely due to long-branch attraction among more deeply branching sequences. All sequences that branch near the base of the tree suffer from relatively high apparent substitution rates and exhibit greater variation in ssu rDNA sequence length. Moreover, the order of the branches leading from the root of the eukaryotic tree to the base of the so-called "crown taxa" is consistent with a sequential attachment, due to "long-branch" effects, of sequences with increasing rates of evolution. These results suggest that the basal eurkaryotic topology drawn from rDNA analyses may be, in reality, an artifact of variation in the rate of molecular evolution among eukaryotic taxa.     

Influenza C virus hemagglutinin: comparison with influenza A and B virus hemagglutinins

The complete nucleotide sequence of the influenza C/California/78 virus RNA 4 was obtained by using cloned cDNA derived from the RNA segment. This gene is 2,071 nucleotides long and can code for a polypeptide of 654 amino acids. Although there are no convincing sequence homologies between RNA 4 and the hemagglutinin genes of influenza A and B viruses, we suggest, on the basis of structural features, that RNA 4 of the influenza C virus codes for the hemagglutinin. The structural features which are common to the hemagglutinins of influenza A, B, and C viruses include (i) a hydrophobic signal peptide, (ii) an arginine cleavage site between the hemagglutinin 1 and 2 subunits, (iii) hydrophobic regions at the amino and carboxyl termini of the hemagglutinin 2 subunit, and (iv) several conserved cysteine residues. Additional evidence that RNA 4 of influenza C virus codes for the hemagglutinin is that the tripeptide Ile-Phe-Gly, known to be present at the amino terminus of the hemagglutinin 2 subunit of influenza C virus, is encoded by RNA 4 at a point immediately adjacent to the presumptive arginine cleavage site. The lack of primary sequence homology between the influenza C virus hemagglutinin and the influenza A or B virus hemagglutinins, which all have similar functions, might be attributed to convergent rather than divergent evolution. However, the structural similarities among the influenza A, B, and C virus hemagglutinins strongly suggest that the three hemagglutinin genes have diverged from a common precursor.     

Linker and/or transmembrane regions of influenza A/Group-1, A/Group-2, and type B virus hemagglutinins are packed differently within trimers

Influenza virus hemagglutinin is a homotrimeric spike glycoprotein crucial for virions' attachment, membrane fusion, and assembly reactions. X-ray crystallography data are available for hemagglutinin ectodomains of various types/subtypes but not for anchoring segments. To get structural information for the linker and transmembrane regions of hemagglutinin, influenza A (H1-H16 subtypes except H8 and H15) and B viruses were digested with bromelain or subtilisin Carlsberg, either within virions or in non-ionic detergent micelles. Proteolytical fragments were analyzed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Within virions, hemagglutinins of most influenza A/Group-1 and type B virus strains were more susceptible to digestion with bromelain and/or subtilisin compared to A/Group-2 hemagglutinins. The cleavage sites were always located in the hemagglutinin linker sequence. In detergent, 1) bromelain cleaved hemagglutinin of every influenza A subtype in the linker region; 2) subtilisin cleaved Group-2 hemagglutinins in the linker region; 3) subtilisin cleaved Group-1 hemagglutinins in the transmembrane region; 4) both enzymes cleaved influenza B virus hemagglutinin in the transmembrane region. We propose that the A/Group-2 hemagglutinin linker and/or transmembrane regions are more tightly associated within trimers than type A/Group-1 and particularly type B ones. This hypothesis is underpinned by spatial trimeric structure modeling performed for transmembrane regions of both Group-1 and Group-2 hemagglutinin representatives. Differential S-acylation of the hemagglutinin C-terminal anchoring segment with palmitate/stearate residues possibly contributes to fine tuning of transmembrane trimer packing and stabilization since decreased stearate amount correlated with deeper digestion of influenza B and some A/Group-1 hemagglutinins.     

海南省2016-2018年A(H1N1)pdm09亚型流感病毒基因特性分析

2009年A(H1N1)pdm09亚型流感病毒在墨西哥暴发,之后在全世界流行。为了解海南省2016-2018年A(H1N1)pdm09亚型流感病毒流行态势,分析血凝素(HA)与神经氨酸酶(NA)基因遗传进化特征与变异情况,本研究从中国流感监测信息系统获取海南省2016-2018年流感病毒病原学监测数据,选取5家流感监测网络实验室分离鉴定的37株A(H1N1)pdm09亚型流感毒株进行HA与NA基因测序,利用MEGA 10.1.8构建HA与NA基因种系进化树,并分析其氨基酸变异情况。结果显示,2016-2018年共出现3次A(H1N1)pdm09亚型流感病毒活动高峰。2017年10月份以后的分离株(4/8)与2018年大部分分离株(21/22)独立于疫苗株A/Michigan/45/2015聚为一个小支,发生20余处HA与NA氨基酸位点变异。与疫苗株A/California/7/2009(2010-2016)相比,2016-2018年流感病毒分离株在HA基因抗原决定簇上发生7处氨基酸变异并有一个潜在糖基化位点,未发现HA基因受体结合位点变异与NA基因耐药性变异。本研究提示,2016-2018年,A(H1N1)pdm09亚型流感病毒逐步发生规律性进化,氨基酸变异频率有增加趋势,今后应持续加强流感病毒病原学监测,密切追踪A(H1N1)pdm09亚型流感病毒基因变异情况,为科学防控提供理论依据。     

The phylogeny of human globin genes investigated by the maximum parsimony method

Gene phylogenetic trees were constructed by the maximum parsimony method for various sets of ninety six globin chain amino acid sequences spanning plant and animal kingdoms. The method, executed by several computer programs, constructed ancestor and descendant globin messengers on tree topologies which required the least number of nucleotide replacements to account for the evolution of the globins. The human myoglobin-hemoglobin divergence was traced to a gene duplication which occurred either in the first vertebrates or earlier yet in the common ancestor of chordates and annelids, the alpha-beta divergence to a gene duplication in the common ancestor of teleosts and tetrapods, the gamma divergence from typical beta chains to a gene duplication in basal therian mammals, and the delta separation from beta to a duplication in the basal catarrhine primates. Evidence was provided by the globin phylogenies for the hominoid affinities of the gibbon and the close phyletic relationship of the African apes to man. Over the period of teleos-tetrapod divergence the globin messengers evolved at an average rate of 18.5 nucleotide replacements per 100 codons per 10⁸ years, a faster rate than most previous estimates. Very fast and very slow rates were encountered in different globin lineages and at different stages of descent, reducing the effectiveness of globins as molecular clocks. Rates increased with gene duplication and decreased after selection discovered useful specializations in the products of genes which had previously been freer to accept mutations. The early eutherian radiation was characterized by rapid rates of globin evolution, but the later hominoid radiation by extremely slow rates. This pattern was related to more complicated grades of internal organization evolving in human ancestors. The types of nucleotide replacements in the globin messengers over the long course of globin evolution did not seem indicative of any special mutational mechanisms.     

Determination of mutation trend in hemagglutinins by means of translation probability between RNA codons and mutated amino acids

In this study, we used the 183 translation probabilities between RNA codons and mutated amino acids to construct the theoretical distributions of mutated amino acids in hemagglutinins of influenza A virus. We then compared the actual distributions of mutated amino acids from 953 hemagglutinins with their theoretical ones. The results demonstrated that mutated amino acids generally follow the direction of the theoretical distributions governed by RNA codons. This, in turn, highlights the mutation trend of amino acids in hemagglutinins and provides a method for estimating possible mutations in a protein according to its theoretical distributions of mutated amino acids.     

A rate-independent technique for analysis of nucleic acid sequences: evolutionary parsimony

The method of evolutionary parsimony--or operator invariants--is a technique of nucleic acid sequence analysis related to parsimony analysis and explicitly designed for determining evolutionary relationships among four distantly related taxa. The method is independent of substitution rates because it is derived from consideration of the group properties of substitution operators rather than from an analysis of the probabilities of substitution in branches of a tree. In both parsimony and evolutionary parsimony, three patterns of nucleotide substitution are associated one-to-one with the three topologically linked trees for four taxa. In evolutionary parsimony, the three quantities are operator invariants. These invariants are the remnants of substitutions that have occurred in the interior branch of the tree and are analogous to the substitutions assigned to the central branch by parsimony. The two invariants associated with the incorrect trees must equal zero (statistically), whereas only the correct tree can have a nonzero invariant. The chi 2-test is used to ascertain the nonzero invariant and the statistically favored tree. Examples, obtained using data calculated with evolutionary rates and branchings designed to camouflage the true tree, show that the method accurately predicts the tree, even when substitution rates differ greatly in neighboring peripheral branches (conditions under which parsimony will consistently fail). As the number of substitutions in peripheral branches becomes fewer, the parsimony and the evolutionary-parsimony solutions converge. The method is robust and easy to use.      

Relationships between tree size, crown shape, gender segregation and sex allocation in Pinus halepensis, a Mediterranean pine tree

Background and Aims

Sex allocation has been studied mainly in small herbaceous plants but much less in monoecious wind-pollinated trees. The aim of this study was to explore changes in gender segregation and sex allocation by Pinus halepensis, a Mediterranean lowland pine tree, within tree crowns and between trees differing in their size or crown shape.

Methods

The production of new male and female cones and sex allocation of biomass, nitrogen and phosphorus were studied. The relationship between branch location, its reproductive status and proxies of branch vigour was also studied.

Key Results

Small trees produced only female cones, but, as trees grew, they produced both male and female cones. Female cones were produced mainly in the upper part of the crown, and male cones in its middle and lower parts. Lateral branch density was correlated with the number of male but not female cones; lateral branches were more dense in large than in small trees and even denser in hemispherical trees. Apical branches grew faster, were thicker and their phosphorus concentration was higher than in lateral shoots. Nitrogen concentration was higher in cone-bearing apical branches than in apical vegetative branches and in lateral branches with or without cones. Allocation to male relative to female function increased with tree size as predicted by sex allocation theory.

Conclusions

The adaptive values of sex allocation and gender segregation patterns in P. halepensis, in relation to its unique life history, are demonstrated and discussed. Small trees produce only female cones that have a higher probability of being pollinated than the probability of male cones pollinating; the female-first strategy enhances population spread. Hemispherical old trees are loaded with serotinous cones that supply enough seeds for post-fire germination; thus, allocation to males is more beneficial than to females.     

Evaluating a non-destructive method for calibrating tree biomass equations derived from tree branching architecture

Key message

Functional branch analysis (FBA) is a promising non-destructive method that can produce accurate tree biomass equations when applied to trees which exhibit fractal branching architecture.

Abstract

Functional branch analysis (FBA) is a promising non-destructive alternative to the standard destructive method of tree biomass equation development. In FBA, a theoretical model of tree branching architecture is calibrated with measurements of tree stems and branches to estimate the coefficients of the biomass equation. In this study, species-specific and mixed-species tree biomass equations were derived from destructive sampling of trees in Western Kenya and compared to tree biomass equations derived non-destructively from FBA. The results indicated that the non-destructive FBA method can produce biomass equations that are similar to, but less accurate than, those derived from standard methods. FBA biomass prediction bias was attributed to the fact that real trees diverged from fractal branching architecture due to highly variable length–diameter relationships of stems and branches and inaccurate scaling relationships for the lengths of tree crowns and trunks assumed under the FBA model.     

Phylogenetic extinction rates and comparative methodology

Species are not independent points for comparative analyses because closely related species share more evolutionary history and are therefore more similar to each other than distantly related species. The extent to which independent-contrast analysis reduces type I and type II statistical error in comparison with cross-species analysis depends on the relative branch lengths in the phylogenetic tree: as deeper branches get relatively long, cross-species analyses have more statistical type I and type II error. Phylogenetic trees reconstructed from extant species, under the assumptions of a branching process with speciation (branching) and extinction rates remaining constant through time, will have relatively longer deep branches as the extinction rate increases relative to the speciation rate. We compare the statistical performance of cross-species and independent-contrast analyses with varying relative extinction rates, and conclude that cross-species comparisons have unacceptable statistical performance, particularly when extinction rates are relatively high.     

Comparative analysis of evolutionary mechanisms of the hemagglutinin and three internal protein genes of influenza B virus: multiple cocirculating lineages and frequent reassortment of the NP, M, and NS genes

Phylogenetic profiles of the genes coding for the hemagglutinin (HA) protein, nucleoprotein (NP), matrix (M) protein, and nonstructural (NS) proteins of influenza B viruses isolated from 1940 to 1998 were analyzed in a parallel manner in order to understand the evolutionary mechanisms of these viruses. Unlike human influenza A (H3N2) viruses, the evolutionary pathways of all four genes of recent influenza B viruses revealed similar patterns of genetic divergence into two major lineages. Although evolutionary rates of the HA, NP, M, and NS genes of influenza B viruses were estimated to be generally lower than those of human influenza A viruses, genes of influenza B viruses demonstrated complex phylogenetic patterns, indicating alternative mechanisms for generation of virus variability. Topologies of the evolutionary trees of each gene were determined to be quite distinct from one another, showing that these genes were evolving in an independent manner. Furthermore, variable topologies were apparently the result of frequent genetic exchange among cocirculating epidemic viruses. Evolutionary analysis done in the present study provided further evidence for cocirculation of multiple lineages as well as sequestering and reemergence of phylogenetic lineages of the internal genes. In addition, comparison of deduced amino acid sequences revealed a novel amino acid deletion in the HA1 domain of the HA protein of recent isolates from 1998 belonging to the B/Yamagata/16/88-like lineage. It thus became apparent that, despite lower evolutionary rates, influenza B viruses were able to generate genetic diversity among circulating viruses through a combination of evolutionary mechanisms involving cocirculating lineages and genetic reassortment by which new variants with distinct gene constellations emerged.     

Long branch effects distort maximum likelihood phylogenies in simulations despite selection of the correct model

The aim of our study was to test the robustness and efficiency of maximum likelihood with respect to different long branch effects on multiple-taxon trees. We simulated data of different alignment lengths under two different 11-taxon trees and a broad range of different branch length conditions. The data were analyzed with the true model parameters as well as with estimated and incorrect assumptions about among-site rate variation. If length differences between connected branches strongly increase, tree inference with the correct likelihood model assumptions can fail. We found that incorporating invariant sites together with Γ distributed site rates in the tree reconstruction (Γ+I) increases the robustness of maximum likelihood in comparison with models using only Γ. The results show that for some topologies and branch lengths the reconstruction success of maximum likelihood under the correct model is still low for alignments with a length of 100,000 base positions. Altogether, the high confidence that is put in maximum likelihood trees is not always justified under certain tree shapes even if alignment lengths reach 100,000 base positions.     

Use of forest canopy by European red squirrelsSciurus vulgaris in Northern Greece: claws and the small branch niche

Moving and standing in trees impose multiple problems to arboreal mammals. Among them, the major ones are the negotiation of slender terminal branches and of large vertical supports. Both microhabitats are important as they have been linked alternatively to the evolutionary loss of claws in early primates. Therefore, rates of use of these different supports by claw-bearing arboreal mammals may offer insights to their actual significance in the adaptive evolution of early primates. In this context, canopy, tree crown, branch size, inclination, and texture use were recorded on four adult free ranging European red squirrelsSciurus vulgaris Linnaeus, 1758 in a mixed coniferous forest in northern Greece.S. vulgaris was mainly arboreal, exploiting the terminal branch zone, using frequently oblique and intermediately textured supports<5 cm and moderately large vertical branches. Furthermore, comparative data from other sciurid species and clawed primates showed positive correlations of small and horizontal support use, and negative ones of vertical support use to body mass. These findings show that keeled functional claws do not impede habitual use of slender branches and may not facilitate efficient climbing on large vertical trunks. These observations partly question the association between habitual use of the small branch niche and primate adaptations and lend support to alternative hypotheses, underscoring the importance of inquiring for more complex mechanisms that lead to the evolution of the unique set of primate morphological adaptations.     

The K tree score: quantification of differences in the relative branch length and topology of phylogenetic trees

SUMMARY: We introduce a new phylogenetic comparison method that measures overall differences in the relative branch length and topology of two phylogenetic trees. To do this, the algorithm first scales one of the trees to have a global divergence as similar as possible to the other tree. Then, the branch length distance, which takes differences in topology and branch lengths into account, is applied to the two trees. We thus obtain the minimum branch length distance or K tree score. Two trees with very different relative branch lengths get a high K score whereas two trees that follow a similar among-lineage rate variation get a low score, regardless of the overall rates in both trees. There are several applications of the K tree score, two of which are explained here in more detail. First, this score allows the evaluation of the performance of phylogenetic algorithms, not only with respect to their topological accuracy, but also with respect to the reproduction of a given branch length variation. In a second example, we show how the K score allows the selection of orthologous genes by choosing those that better follow the overall shape of a given reference tree. AVAILABILITY: http://molevol.ibmb.csic.es/Ktreedist.html