Amino acid compositions contribute to the proteins' evolution under the influence of their abundances and genomic GC content

Using yeasts for the studies of nonfunctional factors in protein evolution

The evolution of protein sequence is driven not only by factors directly related to protein function and shape but also by nonfunctional factors. Such factors in protein evolution might be categorized as those connected to energetic costs, synthesis efficiency, and avoidance of misfolding and toxicity. A common approach to studying them is correlational analysis contrasting them with some characteristics of the protein, like amino acid composition, but these features are interdependent. To avoid possible bias, empirical studies are needed, and not enough work has been done to date. In this review, we describe the role of nonfunctional factors in protein evolution and present an experimental approach using yeast as a suitable model organism. The focus of the proposed approach is on the potential negative impact on the fitness of mutations that change protein properties not related to function and the frequency of mutations that change these properties. Experimental results of testing the misfolding avoidance hypothesis as an explanation for why highly expressed proteins evolve slowly are inconsistent with correlational research results. Therefore, more efforts should be made to empirically test the effects of nonfunctional factors in protein evolution and to contrast these results with the results of the correlational analysis approach.

The genome of Candidatus phytoplasma ziziphi provides insights into their biological characteristics

Phytoplasmas are obligate cell wall-less prokaryotic bacteria that primarily multiply in plant phloem tissue. Jujube witches' broom (JWB) associated with phytoplasma is a destructive disease of jujube (Ziziphus jujuba Mill.). Here we report the complete 'Candidatus Phytoplasma ziziphi' chromosome of strain Hebei-2018, which is a circular genome of 764,108-base pairs with 735 predicted CDS. Notably, extra 19,825 bp (from 621,995 to 641,819 bp) compared to the previously reported one complements the genes involved in glycolysis, such as pdhA, pdhB, pdhC, pdhD, ackA, pduL and LDH. The synonymous codon usage bias (CUB) patterns by using comparative genomics analysis among the 9 phytoplasmas were similar for most codons. The ENc-GC3s analysis among the 9 phytoplasmas showed a greater effect under the selection on the CUBs of phytoplasmas genes than mutation and other factors. The genome exhibited a strongly reduced ability in metabolic synthesis, while the genes encoding transporter systems were well developed. The genes involved in sec-dependent protein translocation system were also identified.The expressions of nine FtsHs encoding membrane associated ATP-dependent Zn proteases and Mn-SodA with redox capacity in the Ca. P. ziziphi was positively correlated with the phytoplasma concentration. Taken together, the genome will not only expand the number of phytoplasma species and provide some new information about Ca. P. ziziphi, but also contribute to exploring its pathogenic mechanism.

Metabolic Diversity and Aero-Tolerance in Anammox Bacteria from Geochemically Distinct Aquifers

Anaerobic ammonium oxidation (anammox) is important for converting bioavailable nitrogen into dinitrogen gas, particularly in carbon-poor environments. However, the diversity and prevalence of anammox bacteria in the terrestrial subsurface-a typically oligotrophic environment-are little understood. To determine the distribution and activity of anammox bacteria across a range of aquifer lithologies and physicochemistries, we analyzed 16S rRNA genes and quantified hydrazine synthase genes and transcripts sampled from 59 groundwater wells and metagenomes and metatranscriptomes from an oxic-to-dysoxic subset. Data indicate that anammox and anammox-associated bacteria (class "Candidatus Brocadiae") are prevalent in the aquifers studied, and that anammox community composition is strongly differentiated by dissolved oxygen (DO), but not ammonia/nitrite. While "Candidatus Brocadiae" diversity decreased with increasing DO, "Candidatus Brocadiae" 16S rRNA genes and hydrazine synthase (hzsB) genes and transcripts were detected across a wide range of bulk groundwater DO concentrations (0 to 10 mg/L). Anammox genes and transcripts correlated significantly with those involved in aerobic ammonia oxidation (amoA), potentially representing a major source of nitrite for anammox. Eight "Candidatus Brocadiae" genomes (63 to 95% complete), representing 2 uncharacterized families and 6 novel species, were reconstructed. Six genomes have genes characteristic of anammox, all for chemolithoautotrophy. Anammox and aerotolerance genes of up to four "Candidatus Brocadiae" genomes were transcriptionally active under oxic and dysoxic conditions, although activity was highest in dysoxic groundwater. The coexpression of nrfAH nitrite reductase genes by "Candidatus Brocadiae" suggests active regeneration of ammonia for anammox. Our findings indicate that anammox bacteria contribute to loss of fixed N across diverse anoxic-to-oxic aquifer conditions, which is likely supported by nitrite from aerobic ammonia oxidation. IMPORTANCE Anammox is increasingly shown to play a major role in the aquatic nitrogen cycle and can outcompete heterotrophic denitrification in environments low in organic carbon. Given that aquifers are characteristically oligotrophic, anammox may represent a major route for the removal of fixed nitrogen in these environments, including agricultural nitrogen, a common groundwater contaminant. Our research confirms that anammox bacteria and the anammox process are prevalent in aquifers and occur across diverse lithologies (e.g., sandy gravel, sand-silt, and volcanic) and groundwater physicochemistries (e.g., various oxygen, carbon, nitrate, and ammonium concentrations). Results reveal niche differentiation among anammox bacteria largely driven by groundwater oxygen contents and provide evidence that anammox is supported by proximity to oxic niches and handoffs from aerobic ammonia oxidizers. We further show that this process, while anaerobic, is active in groundwater characterized as oxic, likely due to the availability of anoxic niches.

Ubiquitousness of Haloferax and Carotenoid Producing Genes in Arabian Sea Coastal Biosystems of India

This study presents a comparative analysis of halophiles from the global open sea and coastal biosystems through shotgun metagenomes (n = 209) retrieved from public repositories. The open sea was significantly enriched with Prochlorococcus and Candidatus pelagibacter. Meanwhile, coastal biosystems were dominated by Marinobacter and Alcanivorax. Halophilic archaea Haloarcula and Haloquandratum, predominant in the coastal biosystem, were significantly (p < 0.05) enriched in coastal biosystems compared to the open sea. Analysis of whole genomes (n = 23,540), retrieved from EzBioCloud, detected crtI in 64.66% of genomes, while cruF was observed in 1.69% Bacteria and 40.75% Archaea. We further confirmed the viability and carotenoid pigment production by pure culture isolation (n = 1351) of extreme halophiles from sediments (n = 410 × 3) sampling at the Arabian coastline of India. All red-pigmented isolates were represented exclusively by Haloferax, resistant to saturated NaCl (6 M), and had >60% G + C content. Multidrug resistance to tetracycline, gentamicin, ampicillin, and chloramphenicol were also observed. Our study showed that coastal biosystems could be more suited for bioprospection of halophiles rather than the open sea.

Molecular footprints of selection effects and whole genome duplication (WGD) events in three blueberry species: detected by transcriptome dataset

BACKGROUND

Both selection effects and whole genome duplication played very important roles in plant speciation and evolution, and to decipher the corresponding molecular footprint has always been a central task of geneticists. Vaccinium is species rich genus that comprised of about 450 species, and blueberry is one of the most important species of Vaccinium genus, which is gaining popularity because of high healthful value. In this article, we aimed to decipher the molecular footprints of natural selection on the single copy genes and WGD events occur in the evolutionary history of blueberry species.

RESULTS

We identified 30,143, 29,922 and 28,891 putative protein coding sequences from 45,535, 42,914 and 43,630 unigenes assembled from the leaves' transcriptome assembly of 19 rabbiteye (T1), 13 southern highbush (T2) and 22 northern highbush (T3) blueberry cultivars. A total of 17, 21 and 27 single copy orthologs were found to undergone positive selection in T1 versus T2, T1 versus T3, and T2 versus T3, respectively, and these orthologs were enriched in metabolic pathways including "Terpenoid backbone biosynthesis", "Valine, leucine and isoleucine biosynthesis", "Butanoate metabolism", "C5-Branched dibasic acid metabolism" "Pantothenate and CoA biosynthesis". We also detected significant molecular footprints of a recent (about 9.04 MYA), medium (about 43.44 MYA) and an ancient (about 116.39 MYA) WGD events that occurred in the evolutionary history of three blueberry species.

CONCLUSION

Some important functional genes revealed positive selection effect in blueberry. At least three rounds of WGD events were detected in the evolutionary history of blueberry species. Our work provides insights about the genetic mechanism of adaptive evolution in blueberry and species radiation of Vaccinium in short geological scale time.

Predictions of Apoptosis Proteins by Integrating Different Features Based on Improving Pseudo-Position-Specific Scoring Matrix

Apoptosis proteins are strongly related to many diseases and play an indispensable role in maintaining the dynamic balance between cell death and division in vivo. Obtaining localization information on apoptosis proteins is necessary in understanding their function. To date, few researchers have focused on the problem of apoptosis data imbalance before classification, while this data imbalance is prone to misclassification. Therefore, in this work, we introduce a method to resolve this problem and to enhance prediction accuracy. Firstly, the features of the protein sequence are captured by combining Improving Pseudo-Position-Specific Scoring Matrix (IM-Psepssm) with the Bidirectional Correlation Coefficient (Bid-CC) algorithm from position-specific scoring matrix. Secondly, different features of fusion and resampling strategies are used to reduce the impact of imbalance on apoptosis protein datasets. Finally, the eigenvector adopts the Support Vector Machine (SVM) to the training classification model, and the prediction accuracy is evaluated by jackknife cross-validation tests. The experimental results indicate that, under the same feature vector, adopting resampling methods remarkably boosts many significant indicators in the unsampling method for predicting the localization of apoptosis proteins in the ZD98, ZW225, and CL317 databases. Additionally, we also present new user-friendly local software for readers to apply; the codes and software can be freely accessed at https://github.com/ruanxiaoli/Im-Psepssm.

Detection, characterization and evolutionary aspects of S54LP of SP (SAP54 Like Protein of Sesame Phyllody): a phytoplasma effector molecule associated with phyllody development in sesame (Sesamum indicum L.)

SAP54, an effector protein secreted by phytoplasmas has been reported to induce phyllody. S54LP of SP (SAP54 Like Protein of Sesame Phyllody), a SAP54 ortholog from phyllody and witches' broom affected sesame (Sesamum indicum L.) was amplified, cloned and sequenced. Comparative sequence and phylogenetic analysis of diverse phytoplasma strains was carried out to delineate the evolution of S54LP of SP. The degree of polymorphism across SAP54 orthologs and the evolutionary forces acting on this effector protein were ascertained. Site-specific selection across SAP54 orthologs was estimated using Fixed Effects Likelihood (FEL) approach. Nonsynonymous substitutions were detected in the SAP54 orthologs' sequences from phytoplasmas belonging to same (sub) group. Phylogenetic analysis based on S54LP of SP grouped phytoplasmas belonging to same 16SrDNA (sub) groups into different clusters. Analysis of selection forces acting on SAP54 orthologs from nine different phytoplasma (sub)groups, affecting plant species belonging to twelve different families across ten countries showed the orthologs to be under purifying (negative) selection. One amino acid residue was found to be under pervasive diversifying (positive) selection and a total of three amino acid sites were found to be under pervasive purifying (negative) selection. The location of these amino acids in the signal peptide and mature protein was studied with an aim to understand their role in protein-protein interaction. Asparagine residues (at positions 68 and 84) were found to be under pervasive purifying selection suggesting their functional importance in the effector protein. Our study suggests lack of coevolution between SAP54 and 16SrDNA. Signal peptide appears to evolve at a rate slightly higher than the mature protein. Overall, SAP54 and its orthologs are evolving under purifying selection confirming their functional importance in phytoplasma virulence.

Comparative Mitogenomic Analysis of the Eurydema Genus in the Context of Representative Pentatomidae (Hemiptera: Heteroptera) Taxa

The family Pentatomidae, the largest within the superfamily Pentatomoidae, comprises about 5,000 species; many of which are economically important pests. Although the phylogeny of Pentatomidae species has been studied using various molecular markers, their phylogenetic relationships remain controversial. Recently, mitochondrial genomes (mitogenomes) have been extensively employed to examine the phylogenetics and evolution of different insects, and in this study, we sequenced complete/near-complete mitochondrial genomes from five shield bug species of Eurydema to gain a better understanding of phylogenetic relationships in the Pentatomidae. The five mitogenomes ranged in length from 15,500 to 16,752 bp and comprised 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region. We compared mitogenomic characteristics of the Pentatomidae and constructed phylogenetic trees using Bayesian inference and maximum likelihood methods. Our results showed that gene arrangements, base composition, start/stop codons, gene overlaps, and RNA structures were conserved within the Pentatomidae and that congeneric species shared more characteristics. Saturation and heterogeneity analyses revealed that our PCGs and PCGRNA datasets were valid for phylogenetic analysis. Phylogenetic analyses showed consistent topologies based on BI and ML methods. These analyses strongly supported that Eurydema species belong to the tribe Strachiini, and formed a sister group with Pentatomini. The relationships among Eurydema species were shown to be consistent with their morphological features. (Strachiini + Pentatomini) was found to be a stable sibling of the clade comprising Cappaeini, Graphosomini, and Carpocorini. Furthermore, our results indicated that Graphosoma rubrolineatum (Heteroptera: Pentatomidae) belongs to the Pentatominae and not the Podopinae.