Global Shifts in Genome and Proteome Composition Are Very Tightly Coupled

High-Throughput Sequencing Enables Rapid Analyses of Nematode Mitochondrial Genomes from an Environmental Sample

Mitochondrial genomes serve as essential tools in evolutionary biology, phylogenetics, and population genetics due to their maternal inheritance, lack of recombination, and conserved structure. Traditional morphological methods for identifying nematodes are often insufficient for distinguishing cryptic species complexes. This study highlights recent advancements in nematode mitochondrial genome research, particularly the impact of long-read sequencing technologies such as Oxford Nanopore. These technologies have facilitated the assembly of mitochondrial genomes from mixed soil samples, overcoming challenges associated with designing specific primers for long PCR amplification across different groups of parasitic nematodes. In this study, we successfully recovered and assembled eleven nematode mitochondrial genomes using long-read sequencing, including those of two plant-parasitic nematode species. Notably, we detected Heterodera cruciferae in Victoria, expanding its known geographic range within Australia. Additionally, short-read sequencing data from a previous draft genome study revealed the presence of the mitochondrial genome of Heterodera filipjevi. Comparative analyses of Heterodera mitogenomes revealed conserved protein-coding genes essential for oxidative phosphorylation, as well as gene rearrangements and variations in transfer RNA placement, which may reflect adaptations to parasitic lifestyles. The consistently high A+T content and strand asymmetry observed across species align with trends reported in related genera. This study demonstrates the utility of long-read sequencing for identifying coexisting nematode species in agricultural fields, providing a rapid, accurate, and comprehensive alternative to traditional diagnostic methods. By incorporating non-target endemic species into public databases, this approach enhances biodiversity records and informs biosecurity strategies. These findings reinforce the potential of mitochondrial genomics to strengthen Australia's as well as the global biosecurity framework against plant-parasitic nematode threats.

The protein domains of vertebrate species in which selection is more effective have greater intrinsic structural disorder

The nearly neutral theory of molecular evolution posits variation among species in the effectiveness of selection. In an idealized model, the census population size determines both this minimum magnitude of the selection coefficient required for deleterious variants to be reliably purged, and the amount of neutral diversity. Empirically, an 'effective population size' is often estimated from the amount of putatively neutral genetic diversity and is assumed to also capture a species' effectiveness of selection. A potentially more direct measure of the effectiveness of selection is the degree to which selection maintains preferred codons. However, past metrics that compare codon bias across species are confounded by among-species variation in %GC content and/or amino acid composition. Here, we propose a new Codon Adaptation Index of Species (CAIS), based on Kullback-Leibler divergence, that corrects for both confounders. We demonstrate the use of CAIS correlations, as well as the Effective Number of Codons, to show that the protein domains of more highly adapted vertebrate species evolve higher intrinsic structural disorder.

Omics data analysis reveals the system-level constraint on cellular amino acid composition

Proteins play a pivotal role in coordinating the functions of organisms, essentially governing their traits, as the dynamic arrangement of diverse amino acids leads to a multitude of folded configurations within peptide chains. Despite dynamic changes in amino acid composition of an individual protein (referred to as AAP) and great variance in protein expression levels under different conditions, our study, utilizing transcriptomics data from four model organisms uncovers surprising stability in the overall amino acid composition of the total cellular proteins (referred to as AACell). Although this value may vary between different species, we observed no significant differences among distinct strains of the same species. This indicates that organisms enforce system-level constraints to maintain a consistent AACell, even amid fluctuations in AAP and protein expression. Further exploration of this phenomenon promises insights into the intricate mechanisms orchestrating cellular protein expression and adaptation to varying environmental challenges.

The protein domains of vertebrate species in which selection is more effective have greater intrinsic structural disorder

The nearly neutral theory of molecular evolution posits variation among species in the effectiveness of selection. In an idealized model, the census population size determines both this minimum magnitude of the selection coefficient required for deleterious variants to be reliably purged, and the amount of neutral diversity. Empirically, an "effective population size" is often estimated from the amount of putatively neutral genetic diversity and is assumed to also capture a species' effectiveness of selection. A potentially more direct measure of the effectiveness of selection is the degree to which selection maintains preferred codons. However, past metrics that compare codon bias across species are confounded by among-species variation in %GC content and/or amino acid composition. Here we propose a new Codon Adaptation Index of Species (CAIS), based on Kullback-Leibler divergence, that corrects for both confounders. We demonstrate the use of CAIS correlations, as well as the Effective Number of Codons, to show that the protein domains of more highly adapted vertebrate species evolve higher intrinsic structural disorder.

Comparison of cysteine content in whole proteomes across the three domains of life

An empirical observation suggests that Giardia lamblia proteins have larger cysteine content than their counterparts in other organisms. As this parasite lacks conventional antioxidant stress systems, it is generally accepted that high cysteine content helps G. lamblia cope with oxygen toxicity, a strategy apparently shared by other organisms. Here, we question whether the high cysteine content in some organisms is genuine or just a simple assumption based on singular observations. To this end, we analyzed the cysteine content in 78 proteomes of organisms spanning the three domains of life. The results indicate that the cysteine content in eukaryota is approximately double that in archaea and bacteria, with G. lamblia among the highest. Atypical cysteine contents were found in a few organisms correlating with specific environmental conditions, supporting the evolutionary amino acid-level selection of amino acid composition.

PaxDb 5.0: Curated Protein Quantification Data Suggests Adaptive Proteome Changes in Yeasts

The "Protein Abundances Across Organisms" database (PaxDb) is an integrative metaresource dedicated to protein abundance levels, in tissue-specific or whole-organism proteomes. PaxDb focuses on computing best-estimate abundances for proteins in normal/healthy contexts and expresses abundance values for each protein in "parts per million" in relation to all other protein molecules in the cell. The uniform data reprocessing, quality scoring, and integrated orthology relations have made PaxDb one of the preferred tools for comparisons between individual datasets, tissues, or organisms. In describing the latest version 5.0 of PaxDb, we particularly emphasize the data integration from various types of raw data and how we expanded the number of organisms and tissue groups as well as the proteome coverage. The current collection of PaxDb includes 831 original datasets from 170 species, including 22 Archaea, 81 Bacteria, and 67 Eukaryota. Apart from detailing the data update, we also present a comparative analysis of the human proteome subset of PaxDb against the two most widely used human proteome data resources: Human Protein Atlas and Genotype-Tissue Expression. Lastly, through our protein abundance data, we reveal an evolutionary trend in the usage of sulfur-containing amino acids in the proteomes of Fungi.

Determinants of associations between codon and amino acid usage patterns of microbial communities and the environment inferred based on a cross-biome metagenomic analysis

Codon and amino acid usage were associated with almost every aspect of microbial life. However, how the environment may impact the codon and amino acid choice of microbial communities at the habitat level is not clearly understood. Therefore, in this study, we analyzed codon and amino acid usage patterns of a large number of environmental samples collected from diverse ecological niches. Our results suggested that samples derived from similar environmental niches, in general, show overall similar codon and amino acid distribution as compared to samples from other habitats. To substantiate the relative impact of the environment, we considered several factors, such as their similarity in GC content, or in functional or taxonomic abundance. Our analysis demonstrated that none of these factors can fully explain the trends that we observed at the codon or amino acid level implying a direct environmental influence on them. Further, our analysis demonstrated different levels of selection on codon bias in different microbial communities with the highest bias in host-associated environments such as the digestive system or oral samples and the lowest level of selection in soil and water samples. Considering a large number of metagenomic samples here we showed that microorganisms collected from similar environmental backgrounds exhibit similar patterns of codon and amino acid usage irrespective of the location or time from where the samples were collected. Thus our study suggested a direct impact of the environment on codon and amino usage of microorganisms that cannot be explained considering the influence of other factors.

Protein structure, amino acid composition and sequence determine proteome vulnerability to oxidation-induced damage

Oxidative stress alters cell viability, from microorganism irradiation sensitivity to human aging and neurodegeneration. Deleterious effects of protein carbonylation by reactive oxygen species (ROS) make understanding molecular properties determining ROS susceptibility essential. The radiation‐resistant bacterium Deinococcus radiodurans accumulates less carbonylation than sensitive organisms, making it a key model for deciphering properties governing oxidative stress resistance. We integrated shotgun redox proteomics, structural systems biology, and machine learning to resolve properties determining protein damage by γ‐irradiation in Escherichia coli and D. radiodurans at multiple scales. Local accessibility, charge, and lysine enrichment accurately predict ROS susceptibility. Lysine, methionine, and cysteine usage also contribute to ROS resistance of the D. radiodurans proteome. Our model predicts proteome maintenance machinery, and proteins protecting against ROS are more resistant in D. radiodurans. Our findings substantiate that protein‐intrinsic protection impacts oxidative stress resistance, identifying causal molecular properties.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

When acting as a reproductive barrier for sympatric speciation, hybrid sterility can only be primary

Animal gametes unite to form a zygote that develops into an adult with gonads that, in turn, produce gametes. Interruption of this germinal cycle by prezygotic or postzygotic reproductive barriers can result in two cycles, each with the potential to evolve into a new species. When the speciation process is complete, members of each species are fully reproductively isolated from those of the other. During speciation a primary barrier may be supported and eventually superceded by a later-appearing secondary barrier. For those holding certain cases of prezygotic isolation to be primary (e.g. elephant cannot copulate with mouse), the onus is to show that they had not been preceded over evolutionary time by periods of postzygotic hybrid inviability (genically determined) or sterility (genically or chromosomally determined). Likewise, the onus is upon those holding cases of hybrid inviability to be primary (e.g. Dobzhansky–Muller epistatic incompatibilities) to show that they had not been preceded by periods, however brief, of hybrid sterility. The latter, when acting as a sympatric barrier causing reproductive isolation, can only be primary. In many cases, hybrid sterility may result from incompatibilities between parental chromosomes that attempt to pair during meiosis in the gonad of their offspring (Winge-Crowther-Bateson incompatibilities). While such incompatibilities have long been observed on a microscopic scale, there is growing evidence for a role of dispersed finer DNA sequence differences (i.e. in base k-mers).

Success of alignment-free oligonucleotide (k-mer) analysis confirms relative importance of genomes not genes in speciation and phylogeny

The utility of DNA sequence substrings (k-mers) in alignment-free phylogenetic classification, including that of bacteria and viruses, is increasingly recognized. However, its biological basis eludes many 21st century practitioners. A path from the 19th century recognition of the informational basis of heredity to the modern era can be discerned. Crick’s DNA ‘unpairing postulate’ predicted that recombinational pairing of homologous DNAs during meiosis would be mediated by short k-mers in the loops of stem-loop structures extruded from classical duplex helices. The complementary ‘kissing’ duplex loops – like tRNA anticodon–codon k-mer duplexes – would seed a more extensive pairing that would then extend until limited by lack of homology or other factors. Indeed, this became the principle behind alignment-based methods that assessed similarity by degree of DNA–DNA reassociation in vitro. These are now seen as less sensitive than alignment-free methods that are closely consistent, both theoretically and mechanistically, with chromosomal anti-recombination models for the initiation of divergence into new species. The analytical power of k-mer differences supports the theses that evolutionary advance sometimes serves the needs of nucleic acids (genomes) rather than proteins (genes), and that such differences can play a role in early speciation events.

The landscape of microbial phenotypic traits and associated genes

Bacteria and Archaea display a variety of phenotypic traits and can adapt to diverse ecological niches. However, systematic annotation of prokaryotic phenotypes is lacking. We have therefore developed ProTraits, a resource containing ∼545 000 novel phenotype inferences, spanning 424 traits assigned to 3046 bacterial and archaeal species. These annotations were assigned by a computational pipeline that associates microbes with phenotypes by text-mining the scientific literature and the broader World Wide Web, while also being able to define novel concepts from unstructured text. Moreover, the ProTraits pipeline assigns phenotypes by drawing extensively on comparative genomics, capturing patterns in gene repertoires, codon usage biases, proteome composition and co-occurrence in metagenomes. Notably, we find that gene synteny is highly predictive of many phenotypes, and highlight examples of gene neighborhoods associated with spore-forming ability. A global analysis of trait interrelatedness outlined clusters in the microbial phenotype network, suggesting common genetic underpinnings. Our extended set of phenotype annotations allows detection of 57 088 high confidence gene-trait links, which recover many known associations involving sporulation, flagella, catalase activity, aerobicity, photosynthesis and other traits. Over 99% of the commonly occurring gene families are involved in genetic interactions conditional on at least one phenotype, suggesting that epistasis has a major role in shaping microbial gene content.

Proteomic indicators of oxidation and hydration state in colorectal cancer

New integrative approaches are needed to harness the potential of rapidly growing datasets of protein expression and microbial community composition in colorectal cancer. Chemical and thermodynamic models offer theoretical tools to describe populations of biomacromolecules and their relative potential for formation in different microenvironmental conditions. The average oxidation state of carbon (Z_C) can be calculated as an elemental ratio from the chemical formulas of proteins, and water demand per residue (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\overline{n}}_{{\mathrm{H}}_{2}\mathrm{O}}$\end{document}n¯H2O) is computed by writing the overall formation reactions of proteins from basis species. Using results reported in proteomic studies of clinical samples, many datasets exhibit higher mean Z_C or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\overline{n}}_{{\mathrm{H}}_{2}\mathrm{O}}$\end{document}n¯H2O of proteins in carcinoma or adenoma compared to normal tissue. In contrast, average protein compositions in bacterial genomes often have lower Z_C for bacteria enriched in fecal samples from cancer patients compared to healthy donors. In thermodynamic calculations, the potential for formation of the cancer-related proteins is energetically favored by changes in the chemical activity of H₂O and fugacity of O₂ that reflect the compositional differences. The compositional analysis suggests that a systematic change in chemical composition is an essential feature of cancer proteomes, and the thermodynamic descriptions show that the observed proteomic transformations in host tissue could be promoted by relatively high microenvironmental oxidation and hydration states.

The Code of Silence: Widespread Associations Between Synonymous Codon Biases and Gene Function

Some mutations in gene coding regions exchange one synonymous codon for another, and thus do not alter the amino acid sequence of the encoded protein. Even though they are often called 'silent,' these mutations may exhibit a plethora of effects on the living cell. Therefore, they are often selected during evolution, causing synonymous codon usage biases in genomes. Comparative analyses of bacterial, archaeal, fungal, and human cancer genomes have found many links between a gene's biological role and the accrual of synonymous mutations during evolution. In particular, highly expressed genes in certain functional categories are enriched with optimal codons, which are decoded by the abundant tRNAs, thus enhancing the speed and accuracy of the translating ribosome. The set of genes exhibiting codon adaptation differs between genomes, and these differences show robust associations to organismal phenotypes. In addition to selection for translation efficiency, other distinct codon bias patterns have been found in: amino acid starvation genes, cyclically expressed genes, tissue-specific genes in animals and plants, oxidative stress response genes, cellular differentiation genes, and oncogenes. In addition, genomes of organisms harboring tRNA modifications exhibit particular codon preferences. The evolutionary trace of codon bias patterns across orthologous genes may be examined to learn about a gene's relevance to various phenotypes, or, more generally, its function in the cell.