Prediction and Analysis of Post-Translational Pyruvoyl Residue Modification Sites from Internal Serines in Proteins

Most of pyruvoyl-dependent proteins observed in prokaryotes and eukaryotes are critical regulatory enzymes, which are primary targets of inhibitors for anti-cancer and anti-parasitic therapy. These proteins undergo an autocatalytic, intramolecular self-cleavage reaction in which a covalently bound pyruvoyl group is generated on a conserved serine residue. Traditional detections of the modified serine sites are performed by experimental approaches, which are often labor-intensive and time-consuming. In this study, we initiated in an attempt for the computational predictions of such serine sites with Feature Selection based on a Random Forest. Since only a small number of experimentally verified pyruvoyl-modified proteins are collected in the protein database at its current version, we only used a small dataset in this study. After removing proteins with sequence identities >60%, a non-redundant dataset was generated and was used, which contained only 46 proteins, with one pyruvoyl serine site for each protein. Several types of features were considered in our method including PSSM conservation scores, disorders, secondary structures, solvent accessibilities, amino acid factors and amino acid occurrence frequencies. As a result, a pretty good performance was achieved in our dataset. The best 100.00% accuracy and 1.0000 MCC value were obtained from the training dataset, and 93.75% accuracy and 0.8441 MCC value from the testing dataset. The optimal feature set contained 9 features. Analysis of the optimal feature set indicated the important roles of some specific features in determining the pyruvoyl-group-serine sites, which were consistent with several results of earlier experimental studies. These selected features may shed some light on the in-depth understanding of the mechanism of the post-translational self-maturation process, providing guidelines for experimental validation. Future work should be made as more pyruvoyl-modified proteins are found and the method should be evaluated on larger datasets. At last, the predicting software can be downloaded from http://www.nkbiox.com/sub/pyrupred/index.html.     

Computational Analysis of the Predicted Evolutionary Conservation of Human Phosphorylation Sites

Protein kinase-mediated phosphorylation is among the most important post-translational modifications. However, few phosphorylation sites have been experimentally identified for most species, making it difficult to determine the degree to which phosphorylation sites are conserved. The goal of this study was to use computational methods to characterize the conservation of human phosphorylation sites in a wide variety of eukaryotes. Using experimentally-determined human sites as input, homologous phosphorylation sites were predicted in all 432 eukaryotes for which complete proteomes were available. For each pair of species, we calculated phosphorylation site conservation as the number of phosphorylation sites found in both species divided by the number found in at least one of the two species. A clustering of the species based on this conservation measure was concordant with phylogenies based on traditional genomic measures. For a subset of the 432 species, phosphorylation site conservation was compared to conservation of both protein kinases and proteins in general. Protein kinases exhibited the highest degree of conservation, while general proteins were less conserved and phosphorylation sites were least conserved. Although preliminary, these data tentatively suggest that variation in phosphorylation sites may play a larger role in explaining phenotypic differences among organisms than differences in the complements of protein kinases or general proteins.     

Evolutionary Dynamics of Variant Genomes of Human Papillomavirus Types 18, 45, and 97

Human papillomavirus type 18 (HPV18) and HPV45 account for approximately 20% of all cervix cancers. We show that HPV18, HPV45, and the recently discovered HPV97 comprise a clade sharing a most recent common ancestor within HPV α7 species. Variant lineages of these HPV types were classified by sequence analysis of the upstream regulatory region/E6 region among cervical samples from a population-based study in Costa Rica, and 27 representative genomes from each major variant lineage were sequenced. Nucleotide variation within HPV18 and HPV45 was 3.82% and 2.39%, respectively, and amino acid variation was 4.73% and 2.87%, respectively. Only 18 nucleotide variations, of which 10 were nonsynonymous, were identified among three HPV97 genomes. Full-genome comparisons revealed maximal diversity between HPV18 African and non-African variants (2.6% dissimilarity), whereas HPV18 Asian-American [E1 (AA)] and European (E2) variants were closely related (less than 0.5% dissimilarity); HPV45 genomes had a maximal difference of 1.6% nucleotides. Using a Bayesian Markov chain Monte Carlo (MCMC) method, the divergence times of HPV18, -45, and -97 from their most recent common ancestors indicated that HPV18 diverged approximately 7.7 million years (Myr) ago, whereas HPV45 and HPV97 split off around 5.7 Myr ago, in a period encompassing the divergence of the great ape species. Variants within the HPV18/45/97 lineages were estimated to have diverged from their common ancestors in the genus Homo within the last 1 Myr (<0.7 Myr). To investigate the molecular basis of HPV18, HPV45, and HPV97 evolution, regression models of codon substitution were used to identify lineages and amino acid sites under selective pressure. The E5 open reading frame (ORF) of HPV18 and the E4 ORFs of HPV18, HPV45, and HPV18/45/97 had nonsynonymous/synonymous substitution rate ratios (d_N/d_S) over 1 indicative of positive Darwinian selection. The L1 ORF of HPV18 genomes had an increased proportion of nonsynonymous substitutions (4.93%; average d_N/d_S ratio [M3] = 0.3356) compared to HPV45 (1.86%; M3 = 0.1268) and HPV16 (2.26%; M3 = 0.1330) L1 ORFs. In contrast, HPV18 and HPV16 genomes had similar amino acid substitution rates within the E1 ORF (2.89% and 3.24%, respectively), while HPV45 E1 was highly conserved (amino acid substitution rate was 0.77%). These data provide an evolutionary history of this medically important clade of HPVs and identify an unexpected divergence of the L1 gene of HPV18 that may have clinical implications for the long-term use of an L1-virus-like particle-based prophylactic vaccine.Papillomaviruses (PVs) are a large family of related viruses with circular double-stranded DNA genomes 8 kb in size. Some PV types cause epithelial hyperplasias ranging from benign exophytic warts to premalignant lesions that can progress to invasive cancer. Among the 61 currently designated alpha human PVs (HPVs), the majority have been isolated from the mucosal surface of the genital or oral region (8, 14). Of these, a select group have oncogenic potential and are associated with cervical cancer (11). Specifically, HPV type 16 (HPV16) and HPV18 have been identified in approximately two-thirds of cervical cancers, this tumor is the second most common cancer in women, and it is the principal cancer of women in developing countries (5, 24, 25, 30, 37).To date, studies of HPV18 variants have identified three lineages corresponding to the continental locations where the viral samples were obtained: European (E), Asian-American (AA), and African (Af) (29). The phylogeny of HPV18 variants is reflective of the migration patterns of Homo sapiens and suggests that HPV18 variant lineages might have diverged through genetic isolation at approximately the same time as Homo sapiens began establishing residence in different continental regions. Previous HPV18 intratypic phylogenetic analyses were limited to partial regions of the genome (3, 7, 29). Nevertheless, studies also suggest that HPV18 variants are associated with different levels of oncogenic potential and persistence and histological tumor types (1, 6, 35, 36, 46).HPV45 and HPV97 are the viral types most closely related to HPV18 and taken together form a clade and share a most recent common ancestor (MRCA). HPV97 is a recently described rare type (8, 17). HPV18 and HPV45 account for approximately 20% of all cervix cancers (25). Although HPV45 is a common type found in cervical cancer, its evolutionary history and sequence variability have not been extensively studied.In this report, 27 complete genomes representing the major variant lineages of HPV18, HPV45, and HPV97 were cloned and/or sequenced from clinical samples. Based on full genomes, the intratype/intertype evolutionary trees of HPV18, HPV45, and HPV97 were constructed. By examining the rate ratio of nonsynonymous (d_N) to synonymous (d_S) substitutions per site, diversifying selection acting on each of the eight protein-encoding regions of HPV18, HPV45, and HPV97 was evaluated. In addition, the times of divergence of HPV18/45/97 variants from their MRCA were investigated. These data provide an evolutionary history of this medically important clade of HPVs.     

Conservation and Periodicity of DNA Bend Sites in Eukaryotic Genomes

DNA bend sites appear every 680 bp on average in the human -and ß-globin gene regions. Although most of theirmolecular nature has not been unraveled, a potential bend coresequence A₂N₈A₂N₈A₂ (A/A/A) and its complementary T₂N₈T₂N₈T₂(T/T/T) appeared preferentially either in or very close to mostof the bend sites, whereas other combinations of A₂ and T₂ dinucleotides,A/T/T + A/A/T, T/T/A + T/A/A and A/T/A + T/A/T, did not. Thedistances between any two of the core sequences in the entireß-globin locus showed a strong bias to a length of701–800 bp and multiples thereof, suggesting that thereis periodicity throughout the locus. This bias was not foundfor other combinations of A₂ and T₂. Again, this periodicitywas identified in many eukaryotic genes, whereas the tendencywas absent in mRNAs and prokaryotic as well as viral genomes.     

Post-Translational Tubulin Modifications in Human Astrocyte Cultures

The cytoskeletal protein tubulin plays an integral role in the functional specialization of many cell types. In the central nervous system, post-translational modifications and the expression of specific tubulin isotypes in neurons have been analyzed in greater detail than in their astrocytic counterparts. In this study, we characterized post-translational specifications of tubulin in human astrocytes using the normal human astrocyte (NHA; Lonza) commercial cell line of fetal origin. Immunocytochemical techniques were implemented in conjunction with confocal microscopy to image class III β-tubulin (βIII-tubulin), acetylated tubulin, and polyglutamylated tubulin using fluorescent antibody probes. Fluorescent probe intensity differences and colocalization were quantitatively assessed with the ‘EBImage’ package for the statistical programming language R. Colocalization analysis revealed that, although both acetylated tubulin and polyglutamylated tubulin showed a high degree of correlation with βIII-tubulin, the correlation with acetylated tubulin was stronger. Quantification and statistical analysis of fluorescence intensity demonstrated that the fluorescence probe intensity ratio for acetylated tubulin/βIII-tubulin was greater than the ratio for polyglutamylated tubulin/βIII-tubulin. The open source GEODATA set GSE819950, comprising RNA sequencing data for the NHA cell line, was mined for the expression of enzymes responsible for tubulin modifications. Our analysis uncovered greater expression at the mRNA level for enzymes reported to function in acetylation and deacetylation as compared to enzymes implicated in glutamylation and deglutamylation. Taken together, the results represent a step toward unraveling the tubulin isotypic expression profile and post-translational modification patterns in astrocytes during human brain development.     

Comparisons of Two Large Phaeoviral Genomes and Evolutionary Implications

The evolution of viral genomes has recently attracted considerable attention. We compare the sequences of two large viral genomes, EsV-1 and FirrV-1, belonging to the family of phaeoviruses which infect different species of marine brown algae. Although their genomes differ substantially in size, these viruses share similar morphologies and similar latent infection cycles. In fact, sequence comparisons show that the viruses have more than 60% of their genes in common. However, the order of genes is completely different in the two genomes, suggesting that extensive recombinational events in addition to several large deletions had occurred during the separate evolutionary routes from a common ancestor. We investigated genes encoding components of signal transduction pathways and genes encoding replicative functions in more detail. We found that the two genomes possess different, although overlapping, sets of genes in both classes, suggesting that different genes from each class were lost, perhaps randomly, after the separate evolution from an ancestral genome. Random loss would also account for the fact that more than one-third of the genes in one viral genome has no counterparts in the other genome. We speculate that the ancestral genome belonged to a cellular organism that had once invaded a primordial brown algal host.     

Pathway Length and Evolutionary Constraint in Amino Acid Biosynthesis

The evolutionary properties of a metabolic network may be determined by the topology of the network. One attribute of pathways that make up the network is the number of enzymatic steps between initial substrates and final products. To determine the effect of pathway length on evolutionary lability of pathway structure, we examined amino acid biosynthetic pathways across 48 sequenced organisms. We demonstrate that longer pathways exhibit lower rates of change in pathway structure than shorter pathways. This finding suggests that increasing complexity may increase constraint on evolutionary change. (Matthew T. Rutter and Rebecca A. Zufall) Both authors contributed equally to this work.     

Post-Translational Modification of Xanthine Dehydrogenase in a Natural Population of DROSOPHILA MELANOGASTER

Second chromosomes of D. melanogaster were isolated from a single natural population, and 40 were analyzed by gel-sieving electrophoresis for the presence of polymorphic loci on chromosome 2 that act to modify xanthine dehydrogenase and/or aldehyde oxidase, whose structural genes map to chromosome 3. Clear evidence of polymorphism for one or more xanthine dehydrogenase modifier loci was obtained.     

基因组功能预测的进化印记方法

改善基因组功能预测方案是目前功能基因组学的迫切问题,生物进化历程会在分子序列上留下相应进化印记－直系同源簇的特异模体,在这一生物学事实的基础上,提出了一个新的基因缚功能预测方法,首先利用进化分析方法构建直系同源簇,再找到各直系同源簇的功能模体,这样可以形成特异的功能模体库,未知基因的功能预测可望通过搜索该功能模体库而得以高效,准确地完成,对５个家族的检验初步证实该方案是可行的。     

Genetic Variability among Complete Human Respiratory Syncytial Virus Subgroup A Genomes: Bridging Molecular Evolutionary Dynamics and Epidemiology

Human respiratory syncytial virus (RSV) is an important cause of severe lower respiratory tract infections in infants and the elderly. In the vast majority of cases, however, RSV infections run mild and symptoms resemble those of a common cold. The immunological, clinical, and epidemiological profile of severe RSV infections suggests a disease caused by a virus with typical seasonal transmission behavior, lacking clear-cut virulence factors, but instead causing disease by modifying the host’s immune response in a way that stimulates pathogenesis. Yet, the interplay between RSV-evoked immune responses and epidemic behavior, and how this affects the genomic evolutionary dynamics of the virus, remains poorly understood. Here, we present a comprehensive collection of 33 novel RSV subgroup A genomes from strains sampled over the last decade, and provide the first measurement of RSV-A genomic diversity through time in a phylodynamic framework. In addition, we map amino acid substitutions per protein to determine mutational hotspots in specific domains. Using Bayesian genealogical inference, we estimated the genomic evolutionary rate to be 6.47×10⁻⁴ (credible interval: 5.56×10⁻⁴, 7.38×10⁻⁴) substitutions/site/year, considerably slower than previous estimates based on G gene sequences only. The G gene is however marked by elevated substitution rates compared to other RSV genes, which can be attributed to relaxed selective constraints. In line with this, site-specific selection analyses identify the G gene as the major target of diversifying selection. Importantly, statistical analysis demonstrates that the immune driven positive selection does not leave a measurable imprint on the genome phylogeny, implying that RSV lineage replacement mainly follows nonselective epidemiological processes. The roughly 50 years of RSV-A genomic evolution are characterized by a constant population size through time and general co-circulation of lineages over many epidemic seasons – a conclusion that might be taken into account when developing future therapeutic and preventive strategies.     

Phylogenetic and Biological Significance of Evolutionary Elements from Metazoan Mitochondrial Genomes

The evolutionary history of living species is usually inferred through the phylogenetic analysis of molecular and morphological information using various mathematical models. New challenges in phylogenetic analysis are centered mostly on the search for accurate and efficient methods to handle the huge amounts of sequence data generated from newer genome sequencing. The next major challenge is the determination of relationships between the evolution of structural elements and their functional implementation, which is largely ignored in previous analyses. Here, we described the discovery of structural elements in metazoan mitochondrial genomes, termed key K-strings, that can serve as a basis for phylogenetic tree construction. Although comprising only a small fraction (0.73%) of all K-strings, these key K-strings are pivotal to the tree construction because they allow for a significant reduction in the computational time required to construct phylogenetic trees, and more importantly, they make significant improvement to the results of phylogenetic inference. The trees constructed from the key K-strings were consistent overall to our current view of metazoan phylogeny and exhibited a more rational topology than the trees constructed by using other conventional methods. Surprisingly, the key K-strings tended to accumulate in the conserved regions of the original sequences, which were most likely due to strong selection pressure. Furthermore, the special structural features of the key K-strings should have some potential applications in the study of the structures and functions relationship of proteins and in the determination of evolutionary trajectory of species. The novelty and potential importance of key K-strings lead us to believe that they are essential evolutionary elements. As such, they may play important roles in the process of species evolution and their physical existence. Further studies could lead to discoveries regarding the relationship between evolution and processes of speciation.     

Molecular Evolution of Human Coronavirus Genomes

Download video (16MB)Help with mp4 files