Nonsynonymous Substitution Rates Research Articles

Enhanced Detection and Molecular Modeling of Adaptive Mutations in SARS-CoV-2 Coding and Non-Coding Regions Using the c/µ Test

Abstract Accurately identifying mutations under beneficial selection in viral genomes is crucial for understanding their molecular evolution and pathogenicity. Traditional methods like the Ka/Ks test, which assesses non-synonymous (Ka) versus synonymous (Ks) substitution rates, assume that synonymous substitutions at synonymous sites are neutral and thus is equal to the mutation rate (µ). Yet, evidence suggests that synonymous sites in translated (TR) regions and untranslated regions (UTR) can be under strong beneficial selection (Ks>µ) and strongly conserved (Ks ≈ 0), leading to false predictions of adaptive mutations from codon-by-codon Ka/Ks analysis. Our previous work used a relative substitution rate test (c/µ, c: substitution rate in UTR/TR, and µ: mutation rate) to identify adaptive mutations in SARS-CoV-2 genome without the neutrality assumption of the synonymous sites. This study refines the c/µ test by optimizing µ value, leading to a smaller set of nucleotide and amino acid sites under beneficial selection in both UTR (11 sites with c/µ>3) and TR (69 nonsynonymous sites: c/µ>3 and Ka/Ks>2.5; 107 synonymous sites: Ks/µ>3). Encouragingly, the top 2 mutations in UTR and 70% of the top nonsynonymous mutations in TR had reported or predicted effects in the literature. Molecular modeling of top adaptive mutations for some critical proteins (S, NSP11 and NSP5) was carried out to elucidate the possible molecular mechanism of their adaptivity.

Comparative mitochondrial genomics of Thelebolaceae in Antarctica: insights into their extremophilic adaptations and evolutionary dynamics

Species of Antarctomyces and Thelebolus (Thelebolaceae), primarily found in Antarctic environments, exhibit psychrophilic adaptations, yet their mitochondrial genomes have not been extensively studied. Furthermore, few studies have compared the mitochondrial genomes of psychrophilic, psychrotrophic, and mesophilic fungi. After successful sequencing and assembly, this study annotated the mitochondrial genomes of Antarctomyces psychrotrophicus CPCC 401038 and Thelebolus microsporus CPCC 401041. We also performed a comparative analysis with the previously characterized mitochondrial genomes of psychrotrophic and mesophilic fungi. The analysis revealed that nad4L was the most conserved gene across the mitochondrial genomes, characterized by its synonymous and non-synonymous substitution rates (Ks and Ka), genetic distance, and GC content and skew within the protein-coding genes (PCGs). Additionally, the mitochondrial genomes of psychrophilic and psychrotrophic fungi showed a higher proportion of protein-coding regions and a lower GC content compared to those of mesophilic fungi, underscoring the genetic basis of cold adaptation. Phylogenetic analyses based on these mitochondrial genes also confirmed the phylogenetic relationships of Thelebolaceae in the class Leotiomycetes. These findings advance our understanding of the phylogenetic relationships and evolutionary dynamics within the family Thelebolaceae, highlighting how different environmental temperatures influence fungal mitochondrial genomic structure and adaptation.

Elucidating the Mitogenomic Blueprint of Pomadasys perotaei from the Eastern Atlantic: Characterization and Matrilineal Phylogenetic Insights into Haemulid Grunts (Teleostei: Lutjaniformes).

The parrot grunt fish, Pomadasys perotaei, has a limited distribution in the Eastern Atlantic Ocean and is an important species in marine capture fisheries across several West African countries. Despite its ecological and economic significance, the mitogenomic information for this species is lacking. This study utilized next-generation sequencing to generate the de novo mitogenome of P. perotaei from Eastern Atlantic. The resulting mitogenome is 16,691 base pairs and includes 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNAs, and an AT-rich control region (CR). Most of the PCGs exhibit nonsynonymous (Ka) and synonymous (Ks) substitution rates of less than '1', indicating strong negative selection across haemulid fishes. The control region of Pomadasys species contains four conserved domains, as seen in other teleost's, with polymorphic nucleotides that can be used to study population structures through the amplification of short mitochondrial gene fragments. Additionally, Bayesian phylogenetic analysis based on PCGs revealed a non-monophyletic clustering pattern of Pomadasys within the haemulid matrilineal tree. Overall, the structural characterization and phylogenetic analysis enhance our understanding of the genetic composition and evolutionary history of Pomadasys species from the Indo-West Pacific and Eastern Atlantic Oceans.

Haemoglobin Gene Repertoire in Teleost and Cichlid Fishes Shaped by Gene Duplications and Genome Rearrangements.

Haemoglobin is a key molecule for oxygen transport in vertebrates. It exhibits remarkable gene diversity in teleost fishes, reflecting adaptation to various aquatic environments. In this study, we present the dynamic evolution of haemoglobin subunit genes based on a comparison of high-quality genome assemblies of 24 vertebrate species, including 17 teleosts (of which six are cichlids). Our findings indicate that teleost genomes contain a range of haemoglobin genes, from as few as five in fugu to as many as 43 in salmon, with the latter being the largest repertoire found in vertebrates. We find evidence that the teleost ancestor had at least four Hbα and three or four Hbβ subunit genes, and that the current gene diversity emerged during teleost radiation, driven primarily by (tandem) gene duplications, genome compaction, and rearrangement dynamics. We provide insights into the genomic organisation of haemoglobin clusters in different teleost species. We further show that the evolution of paralogous rhbdf1 genes flanking both teleost clusters (LA and MN) supports the hypothesis for the origin of the LA cluster by rearrangement within teleosts, rather than by the teleost specific whole-genome duplication. We specifically focus on cichlid fishes, where adaptation to low oxygen environment plays role in species diversification. Our analysis of six cichlid genomes, including Pungu maclareni from the Barombi Mbo crater lake, for which we sequenced a representative genome, reveals 18-32 copies of the Hb genes, and elevated rates of non-synonymous substitutions compared to other teleosts. Overall, this work facilitates a deeper understanding of how haemoglobin genes contribute to the adaptive potential of teleosts.

Natural selection on apical membrane antigen 1 (AMA1) of an emerging zoonotic malaria parasite Plasmodium inui

Apical membrane antigen 1 (AMA1) of malaria parasites plays an important role in host cell invasion. Antibodies to AMA1 can inhibit malaria merozoite invasion of erythrocytes while vaccine-induced specific cytotoxic T cell responses to this protein are associated with clinical protection. Polymorphisms in AMA1 of Plasmodium falciparum (PfAMA1) and P. vivax (PvAMA1) are of concern for vaccine development. To date, little is known about sequence diversity in ama1 of P. inui (Piama1), an emerging zoonotic malaria parasite. In this study, 80 complete Piama1 coding sequences were obtained from 57 macaques in Thailand that defined 60 haplotypes clustering in two phylogenetic lineages. In total, 74 nucleotide substitutions were identified and distributed unevenly across the gene. Blockwise analysis of the rates of synonymous (dS) and nonsynonymous (dN) nucleotide substitutions did not show a significant deviation from neutrality among Thai isolates. However, significantly negative Tajima’s D values were detected in domain I and the loop region of domain II, implying purifying selection. Codon-based analysis of dN/dS has identified 12 and 14 codons under positive and negative selections, respectively. Meanwhile, 85 amino acid substitutions were identified among 80 Thai and 11 non-Thai PiAMA1 sequences. Of these, 48 substituted residues had a significant alteration in physicochemical properties, suggesting positive selection. More than half of these positively selected amino acids (32 of 48) corresponded to the predicted B-cell or T-cell epitopes, suggesting that selective pressure could be mediated by host immunity. Importantly, 14 amino acid substitutions were singletons and predicted to be deleterious that could be subject to ongoing purifying selection or elimination. Besides genetic drift and natural selection, intragenic recombination identified in domain II could generate sequence variation in Piama1. It is likely that malarial ama1 exhibits interspecies differences in evolutionary histories. Knowledge of the sequence diversity of the Piama1 locus further provides an evolutionary perspective of this important malaria vaccine candidate.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity.

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection

BackgroundNematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures.ResultsGenerally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection.ConclusionsNematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae)

BackgroundCorydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes.ResultsWe obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit β (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera.ConclusionsThe C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

GTDrift: a resource for exploring the interplay between genetic drift, genomic and transcriptomic characteristics in eukaryotes.

We present GTDrift, a comprehensive data resource that enables explorations of genomic and transcriptomic characteristics alongside proxies of the intensity of genetic drift in individual species. This resource encompasses data for 1506 eukaryotic species, including 1413 animals and 93 green plants, and is organized in three components. The first two components contain approximations of the effective population size, which serve as indicators of the extent of random genetic drift within each species. In the first component, we meticulously investigated public databases to assemble data on life history traits such as longevity, adult body length and body mass for a set of 979 species. The second component includes estimations of the ratio between the rate of non-synonymous substitutions and the rate of synonymous substitutions (dN/dS) in protein-coding sequences for 1324 species. This ratio provides an estimate of the efficiency of natural selection in purging deleterious substitutions. Additionally, we present polymorphism-derived N e estimates for 66 species. The third component encompasses various genomic and transcriptomic characteristics. With this component, we aim to facilitate comparative transcriptomics analyses across species, by providing easy-to-use processed data for more than 16 000 RNA-seq samples across 491 species. These data include intron-centered alternative splicing frequencies, gene expression levels and sequencing depth statistics for each species, obtained with a homogeneous analysis protocol. To enable cross-species comparisons, we provide orthology predictions for conserved single-copy genes based on BUSCO gene sets. To illustrate the possible uses of this database, we identify the most frequently used introns for each gene and we assess how the sequencing depth available for each species affects our power to identify major and minor splice variants.

Auto-phylo v2 and auto-phylo-pipeliner: building advanced, flexible, and reusable pipelines for phylogenetic inferences, estimation of variability levels and identification of positively selected amino acid sites.

The vast amount of genome sequence data that is available, and that is predicted to drastically increase in the near future, can only be efficiently dealt with by building automated pipelines. Indeed, the Earth Biogenome Project will produce high-quality reference genome sequences for all 1.8 million named living eukaryote species, providing unprecedented insight into the evolution of genes and gene families, and thus on biological issues. Here, new modules for gene annotation, further BLAST search algorithms, further multiple sequence alignment methods, the adding of reference sequences, further tree rooting methods, the estimation of rates of synonymous and nonsynonymous substitutions, and the identification of positively selected amino acid sites, have been added to auto-phylo (version 2), a recently developed software to address biological problems using phylogenetic inferences. Additionally, we present auto-phylo-pipeliner, a graphical user interface application that further facilitates the creation and running of auto-phylo pipelines. Inferences on S-RNase specificity, are critical for both cross-based breeding and for the establishment of pollination requirements. Therefore, as a test case, we develop an auto-phylo pipeline to identify amino acid sites under positive selection, that are, in principle, those determining S-RNase specificity, starting from both non-annotated Prunus genomes and sequences available in public databases.

Full-length transcriptome analysis of Ophioglossum vulgatum: effects of experimentally identified chloroplast gene clusters on expression and evolutionary patterns.

Genes with similar or related functions in chloroplasts are often arranged in close proximity, forming clusters on chromosomes. These clusters are transcribed coordinated to facilitate the expression of genes with specific function. Our previous study revealed a significant negative correlation between the chloroplast gene expression level of the rare medicinal fern Ophioglossum vulgatum and its evolutionary rates as well as selection pressure. Therefore, in this study, we employed a combination of SMRT and Illumina sequencing technology to analyze the full-length transcriptome sequencing of O. vulgatum for the first time. In particular, we experimentally identified gene clusters based on transcriptome data and investigated the effects of chloroplast gene clustering on expression and evolutionary patterns. The results revealed that the total sequenced data volume of the full-length transcriptome of O. vulgatum amounted to 71,950,652,163bp, and 110 chloroplast genes received transcript coverage. Nine different types of gene clusters were experimentally identified in their transcripts. The chloroplast cluster genes may cause a decrease in non-synonymous substitution rate and selection pressure, as well as a reduction in transversion rate, transition rate, and their ratio. While expression levels of chloroplast cluster genes in leaf, sporangium, and stem would be relatively elevated. The Mann-Whitney U test indicated statistically significant in the selection pressure, sporangia and leaves groups (P < 0.05). We have contributed novel full-length transcriptome data resources for ferns, presenting new evidence on the effects of chloroplast gene clustering on expression land evolutionary patterns, and offering new theoretical support for transgenic research through gene clustering.

Characterization of the complete mitochondrial genome of Paecilomyces variotii and comparative evolutionary mitochondriomics of 36 fungi.

Paecilomyces variotii has important economic value in stimulating crop growth, biodegradation, and other aspects. Up to now, there are no research reports on its mitochondrial genome. The mitochondrial genome of Paecilomyces variotii was determined with the next-generation sequencing method (Illumina, NovaSeq), and its characteristics were analyzed using various bioinformatics approaches. The length of complete mitochondrial genome sequence of P. variotii is 40,965bp and consists of 14 protein-coding genes, 2 ribosomal RNA genes, 1 ribosomal protein S3 gene, 26 transport RNA genes. The results of phylogenetics analysis using Bayesian inference and Maximum likelihood methods showed that P. variotii belongs to the Eurotiales order in the Thermoascaceae family, and 9 genera within the Eurotiomycetes class were effectively distinguished with high support rates (bootstrap value > 92% and posterior probabilities > 99%). The analysis of synonymous substitution rates and nonsynonymous substitution rates indicated that the Ka/Ks values of the 14 PCGs in the mitochondrial genomes of the two orders in the Eurotiomycetes class ranged from 0 to 0.4333. This study revealed the structural and sequence information characteristics of the mitochondrial genome of P. variotii, and the phylogenetic results strongly support its classification within the family Thermoascaceae, consistent with traditional morphological taxonomy studies. The 14 PCGs in the mitochondrial genomes of the two orders in the Eurotiomycetes class are subject to strong purifying (negative) selection. The results of this research provides an important molecular basis for the development of genomics, evolutionary genetics and molecular markers of P. variotii in the future.

Genome-wide analysis of the peanut CaM/CML gene family reveals that the AhCML69 gene is associated with resistance to Ralstonia solanacearum

BackgroundCalmodulins (CaMs)/CaM-like proteins (CMLs) are crucial Ca2+-binding sensors that can decode and transduce Ca2+ signals during plant development and in response to various stimuli. The CaM/CML gene family has been characterized in many plant species, but this family has not yet been characterized and analyzed in peanut, especially for its functions in response to Ralstonia solanacearum. In this study, we performed a genome-wide analysis to analyze the CaM/CML genes and their functions in resistance to R. solanacearum.ResultsHere, 67, 72, and 214 CaM/CML genes were identified from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively. The genes were divided into nine subgroups (Groups I-IX) with relatively conserved exon‒intron structures and motif compositions. Gene duplication, which included whole-genome duplication, tandem repeats, scattered repeats, and unconnected repeats, produced approximately 81 pairs of homologous genes in the AhCaM/CML gene family. Allopolyploidization was the main reason for the greater number of AhCaM/CML members. The nonsynonymous (Ka) versus synonymous (Ks) substitution rates (less than 1.0) suggested that all homologous pairs underwent intensive purifying selection pressure during evolution. AhCML69 was constitutively expressed in different tissues of peanut plants and was involved in the response to R. solanacearum infection. The AhCML69 protein was localized in the cytoplasm and nucleus. Transient overexpression of AhCML69 in tobacco leaves increased resistance to R. solanacearum infection and induced the expression of defense-related genes, suggesting that AhCML69 is a positive regulator of disease resistance.ConclusionsThis study provides the first comprehensive analysis of the AhCaM/CML gene family and potential genetic resources for the molecular design and breeding of peanut bacterial wilt resistance.

Comparative Mitogenome Analyses of Fifteen Ramshorn Snails and Insights into the Phylogeny of Planorbidae (Gastropoda: Hygrophila).

Ramshorn snails from the family Planorbidae are important freshwater snails due to their low trophic level, and some of them act as intermediate hosts for zoonotic trematodes. There are about 250 species from 40 genera of Planorbidae, but only 14 species from 5 genera (Anisus, Biomphalaria, Bulinus, Gyraulus, and Planorbella) have sequenced complete mitochondrial genomes (mitogenomes). In this study, we sequenced and assembled a high-quality mitogenome of a ramshorn snail, Polypylis sp. TS-2018, which represented the first mitogenome of the genus. The mitogenome of Polypylis sp. TS-2018 is 13,749 bp in length, which is shorter than that of most gastropods. It contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and 2 ribosomal RNA (rRNA). We compared mitogenome characteristics, selection pressure, and gene rearrangement among all of the available mitogenomes of ramshorn snails. We found that the nonsynonymous and synonymous substitution rates (Ka/Ks) of most PCGs indicated purifying and negative selection, except for atp8 of Anisus, Biomphalaria, and Gyraulus, which indicated positive selection. We observed that transpositions and reverse transpositions occurred on 10 tRNAs and rrnS, which resulted in six gene arrangement types. We reconstructed the phylogenetic trees using the sequences of PCGs and rRNAs and strongly supported the monophyly of each genus, as well as three tribes in Planorbidae. Both the gene rearrangement and phylogenetic results suggested that Polypylis had a close relationship with Anisus and Gyraulus, while Bulinus was the sister group to all of the other genera. Our results provide useful data for further investigation of species identification, population genetics, and phylogenetics among ramshorn snails.

Genome reduction occurred in early Prochlorococcus with an unusually low effective population size.

In the oligotrophic sunlit ocean, the most abundant free-living planktonic bacterial lineages evolve convergently through genome reduction. The cyanobacterium Prochlorococcus responsible for 10% global oxygen production is a prominent example. The dominant theory known as "genome streamlining" posits that they have extremely large effective population sizes (Ne) such that selection for metabolic efficiency acts to drive genome reduction. Because genome reduction largely took place anciently, this theory builds on the assumption that their ancestors' Ne was similarly large. Constraining Ne for ancient ancestors is challenging because experimental measurements of extinct organisms are impossible and alternatively reconstructing ancestral Ne with phylogenetic models gives large uncertainties. Here, we develop a new strategy that leverages agent-based modeling to simulate the changes in the genome-wide ratio of radical to conservative nonsynonymous nucleotide substitution rate (dR/dC) in a possible range of Ne in ancestral populations. This proxy shows expected increases with decreases of Ne only when Ne falls to about 10k - 100k or lower, magnitudes characteristic of Ne of obligate endosymbiont species where drift drives genome reduction. Our simulations therefore strongly support a scenario where the primary force of Prochlorococcus genome reduction is drift rather than selection.

Large-scale Pan Genomic Analysis of Mycobacterium tuberculosis Reveals Key Insights Into Molecular Evolutionary Rate of Specific Processes and Functions.

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), an infectious disease that is a major killer worldwide. Due to selection pressure caused by the use of antibacterial drugs, Mtb is characterised by mutational events that have given rise to multi drug resistant (MDR) and extensively drug resistant (XDR) phenotypes. The rate at which mutations occur is an important factor in the study of molecular evolution, and it helps understand gene evolution. Within the same species, different protein-coding genes evolve at different rates. To estimate the rates of molecular evolution of protein-coding genes, a commonly used parameter is the ratio dN/dS, where dN is the rate of non-synonymous substitutions and dS is the rate of synonymous substitutions. Here, we determined the estimated rates of molecular evolution of select biological processes and molecular functions across 264 strains of Mtb. We also investigated the molecular evolutionary rates of core genes of Mtb by computing the dN/dS values, and estimated the pan genome of the 264 strains of Mtb. Our results show that the cellular amino acid metabolic process and the kinase activity function evolve at a significantly higher rate, while the carbohydrate metabolic process evolves at a significantly lower rate for M. tuberculosis. These high rates of evolution correlate well with Mtb physiology and pathogenicity. We further propose that the core genome of M. tuberculosis likely experiences varying rates of molecular evolution which may drive an interplay between core genome and accessory genome during M. tuberculosis evolution.

Genetic Diversity and Sequence Conservation of Peptide-Binding Regions of MHC Class I Genes in Pig, Cattle, Chimpanzee, and Human.

Comparative analyses of MHC gene diversity and evolution across different species could offer valuable insights into the evolution of MHC genes. Intra- and inter-species sequence diversity and conservation of 12 classical major histocompatibility complex (MHC) class I genes from cattle, chimpanzees, pigs, and humans was analyzed using 20 representative allelic groups for each gene. The combined analysis of paralogous loci for each species revealed that intra-locus amino-acid sequence variations in the peptide-binding region (PBR) of MHC I genes did not differ significantly between species, ranging from 8.44% for SLA to 10.75% for BoLA class I genes. In contrast, intraspecies differences in the non-PBRs of these paralogous genes were more pronounced, varying from 4.59% for SLA to 16.89% for HLA. Interestingly, the Shannon diversity index and rate of nonsynonymous substitutions for PBR were significantly higher in SLA and BoLA than those in Patr and HLA. Analysis of peptide-binding pockets across all analyzed MHC class I genes of the four species indicated that pockets A and E showed the lowest and highest diversity, respectively. The estimated divergence times suggest that primate and artiodactyl MHC class I genes diverged 60.41 Mya, and BoLA and SLA genes diverged 35.34 Mya. These results offer new insights into the conservation and diversity of MHC class I genes in various mammalian species.

Genomic estimates of mutation and substitution rates contradict the evolutionary speed hypothesis of the latitudinal diversity gradient.

The latitudinal diversity gradient (LDG) refers to a decrease in biodiversity from the equator to the poles. The evolutionary speed hypothesis, backed by the metabolic theory of ecology, asserts that nucleotide mutation and substitution rates per site per year are higher and thereby speciation rates are higher at higher temperatures, generating the LDG. However, prior empirical investigations of the relationship between the temperature and mutation or substitution rate were based on a few genes and the results were mixed. We here revisit this relationship using genomic data. No significant correlation between the temperature and mutation rate is found in 13 prokaryotes or in 107 eukaryotes. An analysis of 234 diverse trios of bacterial taxa indicates that the synonymous substitution rate is not significantly associated with the growth temperature. The same data, however, reveal a significant negative association between the nonsynonymous substitution rate and temperature, which is explainable by a larger fraction of detrimental nonsynonymous mutations at higher temperatures due to a stronger demand for protein stability. We conclude that the evolutionary speed hypothesis of the LDG is unsupported by genomic data and advise that future mechanistic studies of the LDG should focus on other hypotheses.

Adaptive evolution in virulence effectors of the rice blast fungus Pyricularia oryzae.

Plant pathogens secrete proteins called effectors that target host cellular processes to promote disease. Recently, structural genomics has identified several families of fungal effectors that share a similar three-dimensional structure despite remarkably variable amino-acid sequences and surface properties. To explore the selective forces that underlie the sequence variability of structurally-analogous effectors, we focused on MAX effectors, a structural family of effectors that are major determinants of virulence in the rice blast fungus Pyricularia oryzae. Using structure-informed gene annotation, we identified 58 to 78 MAX effector genes per genome in a set of 120 isolates representing seven host-associated lineages. The expression of MAX effector genes was primarily restricted to the early biotrophic phase of infection and strongly influenced by the host plant. Pangenome analyses of MAX effectors demonstrated extensive presence/absence polymorphism and identified gene loss events possibly involved in host range adaptation. However, gene knock-in experiments did not reveal a strong effect on virulence phenotypes suggesting that other evolutionary mechanisms are the main drivers of MAX effector losses. MAX effectors displayed high levels of standing variation and high rates of non-synonymous substitutions, pointing to widespread positive selection shaping the molecular diversity of MAX effectors. The combination of these analyses with structural data revealed that positive selection acts mostly on residues located in particular structural elements and at specific positions. By providing a comprehensive catalog of amino acid polymorphism, and by identifying the structural determinants of the sequence diversity, our work will inform future studies aimed at elucidating the function and mode of action of MAX effectors.