An in-silico study of alphaherpesviruses ICP0 genes: positive selection or strong mutational GC-pressure?

The purpose of our work was to analyze the case of the strong mutational GC-pressure influence on the ratio between nonsynonymous (DN) and synonymous (DS) distances (DN/DS ratio). We have used as the material the genes coding for ICP0 from five completely sequenced genomes of simplexviruses. DN/DS ratio, total GC-content (G + C), and GC-content in first, second, and third codon positions (1GC, 2GC, and 3GC, respectively) have been calculated separately for exon 2, nonconserved part of exon 3, and conserved part of exon 3 from ICP0 genes. Results showed that DN is more than DS only in the conserved part of exon 3 of ICP0 genes from cercopithecine herpesvirus 2 and cercopithecine herpesvirus 16. However, the cause of this result (DN/DS = 2.54) is the GC-pressure acting on the coding districts with 3GC = 99% rather than the biological process called positive selection. Only in these two viruses, because of the strong GC-pressure, 3GC has reached 99% in the conserved part of ICP0 exon 3, and so nucleotide substitutions that increase the GC-content practically cannot occur in third codon positions, where most substitutions are synonymous. In this case, GC-pressure has a substrate for nucleotide substitutions only in first and second codon positions, where most substitutions are nonsynonymous.     

Stabilization of secondary structure elements by specific combinations of hydrophilic and hydrophobic amino acid residues is more important for proteins encoded by GC-poor genes

Stabilization of secondary structure elements by specific combinations of hydrophobic and hydrophilic amino acids has been studied by the way of analysis of pentapeptide fragments from twelve partial bacterial proteomes. PDB files describing structures of proteins from species with extremely high and low genomic GC-content, as well as with average G + C were included in the study. Amino acid residues in 78,009 pentapeptides from alpha helices, beta strands and coil regions were classified into hydrophobic and hydrophilic ones. The common propensity scale for 32 possible combinations of hydrophobic and hydrophilic amino acid residues in pentapeptide has been created: specific pentapeptides for helix, sheet and coil were described. The usage of pentapeptides preferably forming alpha helices is decreasing in alpha helices of partial bacterial proteomes with the increase of the average genomic GC-content in first and second codon positions. The usage of pentapeptides preferably forming beta strands is increasing in coil regions and in helices of partial bacterial proteomes with the growth of the average genomic GC-content in first and second codon positions. Due to these circumstances the probability of coil-sheet and helix-sheet transitions should be increased in proteins encoded by GC-rich genes making them prone to form amyloid in certain conditions. Possible causes of the described fact that importance of alpha helix and coil stabilization by specific combinations of hydrophobic and hydrophilic amino acids is growing with the decrease of genomic GC-content have been discussed.     

Study of Completed Archaeal Genomes and Proteomes： Hypothesis of Strong Mutational AT Pressure Existed in Their Common Predecessor

The number of completely sequenced archaeal genomes has been sufficient for a large-scale bioinformatic study.We have conducted analyses for each coding region from 36 archaeal genomes using the original CGS algorithm by calculating the total GC content(G+C),GC content in first,second and third codon positions as well as in fourfold and twofold degenerated sites from third codon positions,levels of arginine codon usage(Arg2:AGA/G;Arg4:CGX),levels of amino acid usage and the entropy of amino acid content distribution.In archaeal genomes with strong GC pressure,arginine is coded preferably by GC-rich Arg4 codons,whereas in most of archaeal genomes with G+C0.6,arginine is coded preferably by AT-rich Arg2 codons.In the genome of Haloquadratum walsbyi,which is closely related to GC-rich archaea,GC content has decreased mostly in third codon positions,while Arg4Arg2 bias still persists.Proteomes of archaeal species carry characteristic amino acid biases:levels of isoleucine and lysine are elevated,while levels of alanine,histidine,glutamine and cytosine are relatively decreased.Numerous genomic and proteomic biases observed can be explained by the hypothesis of previously existed strong mutational AT pressure in the common predecessor of all archaea.     

Percent of highly immunogenic amino acid residues forming B-cell epitopes is higher in homologous proteins encoded by GC-rich genes

We analyzed the dependence of the percent of highly immunogenic amino acid residues included in B-cell epitopes of homologous proteins on the GC-content (G+C) of genes coding for them in twenty-seven lineages of proteins (and subsequent genes), which belong to seven Varicello and five Simplex viruses. We found out that proteins encoded by genes of a high GC-content usually contain more targets for humoral immune response than their homologs encoded by GC-poor genes. This tendency is characteristic not only to the lineages of glycoproteins, which are the main targets for humoral immune response against Simplex and Varicello viruses, but also to the lineages of capsid proteins and even "housekeeping" enzymes. The percent of amino acids included in linear B-cell epitopes has been predicted for 324 proteins by BepiPred algorithm (www.cbs.dtu.dk/services/BepiPred), the percent of highly immunogenic amino acids included in discontinuous B-cell epitopes and the percent of exposed amino acid residues have been predicted by Epitopia algorithm (http://epitopia.tau.ac.il/). Immunological consequences of the directional mutational GC-pressure are mostly due to the decrease in the total usage of highly hydrophobic amino acids and due to the increase in proline and glycine levels of usage in proteins. The weaker the negative selection on amino acid substitutions caused by symmetric mutational pressure, the higher the slope of direct dependence of the percent of highly immunogenic amino acids included in B-cell epitopes on G+C.     

Larvae of related Diptera species from thermally contrasting habitats exhibit continuous up-regulation of heat shock proteins and high thermotolerance

A population of Stratiomys japonica, a species belonging to the family Stratiomyidae (Diptera), common name ‘soldier flies’, occurs in a hot volcanic spring, which is apparently among the most inhospitable environments for animals because of chemical and thermal conditions. Larvae of this species, which naturally often experience temperatures more than 40 °C, have constitutively high concentrations of the normally inducible heat-shock protein Hsp70, but very low level of corresponding mRNA. Larvae of three other species of the same family, Stratiomys singularior, Nemotelus bipunctatus and Oxycera pardalina, are confined to different type semi-aquatic habitats with contrasting thermal regime. However, all of them shared the same pattern of Hsp70 expression. Interestingly, heat-shock treatment of S. japonica larvae activates heat-shock factor and significantly induces Hsp70 synthesis, whereas larvae of O. pardalina, a species from constant cold environment, produce significantly less Hsp70 in response to heat shock. Adults of the four species also exhibit lower, but detectable levels of Hsp70 without heat shock. Larvae of all species studied have very high tolerance to temperature stress in comparison with other Diptera species investigated, probably representing an inherent adaptive feature of all Stratiomyidae enabling successful colonization of highly variable and extreme habitats.     

Molecular phylogeny of fungus gnats (Diptera: Mycetophilidae) revisited: position of Manotinae,Metanepsiini, and other enigmatic taxa as inferred from multigene analysis

The phylogeny of selected genera from four subfamilies of fungus gnats (Diptera: Mycetophilidae) – Manotinae, Leiinae, Sciophilinae and Gnoristinae (including Metanepsiini) – is reconstructed based on the combined analysis of five mitochondrial (12S, 16S, COI, COII, cytB) and two nuclear (28S, ITS2) gene markers. Results of the different analyses all support Manotinae as a monophyletic group, with Leiinae as the sister group. Allactoneura DeMeijere is nested in the monophyletic and strongly supported clade of Leiinae. The tribe Metanepsiini is revealed as paraphyletic and the genera Metanepsia Edwards and Chalastonepsia Søli do not appear to be closely related. The genera Docosia Winnertz, Ectrepesthoneura Enderlein, Novakia Strobl and Syntemna Winnertz were placed with a group of genera included traditionally in the Gnoristinae. The monophyly of Dziedzickia Johannsen and Phthinia Winnertz is not supported. The genera of Sciophilinae (excluding Paratinia Mik but including Eudicrana Loew) form a monophyletic group in the Bayesian model.