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.

Genomic Analyses Uncover Evolutionary Features of Influenza A/H3N2 Viruses in Yunnan Province, China, from 2017 to 2022

Influenza A viruses evolve at a high rate of nucleotide substitution, thereby requiring continuous monitoring to determine the efficacy of vaccines and antiviral drugs. In the current study, we performed whole-genome sequencing analyses of 253 influenza A/H3N2 strains from Yunnan Province, China, during 2017-2022. The hemagglutinin (HA) segments of Yunnan A/H3N2 strains isolated during 2017-2018 harbored a high genetic diversity due to heterogeneous distribution across branches. The mutation regularity of the predominant antigenic epitopes of HA segments in Yunnan was inconsistent in different years. Some important functional mutations in gene segments associated with viral adaptation and drug tolerance were revealed. The rapid genomic evolution of Yunnan A/H3N2 strains from 2017 to 2022 mainly concentrated on segments, i.e., matrix protein 2 (M2), non-structural protein 1 (NS1), neuraminidase (NA), NS2, and HA, with a high overall non-synonymous/synonymous substitution ratio (dN/dS). Our results highlighted a decline in vaccine efficacy against the A/H3N2 circulating strains, particularly against the Yunnan 2021-2022 A/H3N2 strains. These findings aid our understanding of evolutionary characteristics and epidemiological monitoring of the A/H3N2 viruses and provide in-depth insights into the protective efficacy of influenza vaccines.

Venomous Noodles: The Evolution of Toxins in Nemertea through Positive Selection and Gene Duplication

Some, probably most and perhaps all, members of the phylum Nemertea are poisonous, documented so far from marine and benthic specimens. Although the toxicity of these animals has been long known, systematic studies on the characterization of toxins, mechanisms of toxicity, and toxin evolution for this group are scarce. Here, we present the first investigation of the molecular evolution of toxins in Nemertea. Using a proteo-transcriptomic approach, we described toxins in the body and poisonous mucus of the pilidiophoran Lineus sanguineus and the hoplonemertean Nemertopsis pamelaroeae. Using these new and publicly available transcriptomes, we investigated the molecular evolution of six selected toxin gene families. In addition, we also characterized in silico the toxin genes found in the interstitial hoplonemertean, Ototyphlonemertes erneba, a meiofaunal taxa. We successfully identified over 200 toxin transcripts in each of these species. Evidence of positive selection and gene duplication was observed in all investigated toxin genes. We hypothesized that the increased rates of gene duplications observed for Pilidiophora could be involved with the expansion of toxin genes. Studies concerning the natural history of Nemertea are still needed to understand the evolution of their toxins. Nevertheless, our results show evolutionary mechanisms similar to other venomous groups.

Fitness effects of mutations to SARS-CoV-2 proteins

Knowledge of the fitness effects of mutations to SARS-CoV-2 can inform assessment of new variants, design of therapeutics resistant to escape, and understanding of the functions of viral proteins. However, experimentally measuring effects of mutations is challenging: we lack tractable lab assays for many SARS-CoV-2 proteins, and comprehensive deep mutational scanning has been applied to only two SARS-CoV-2 proteins. Here, we develop an approach that leverages millions of publicly available SARS-CoV-2 sequences to estimate effects of mutations. We first calculate how many independent occurrences of each mutation are expected to be observed along the SARS-CoV-2 phylogeny in the absence of selection. We then compare these expected observations to the actual observations to estimate the effect of each mutation. These estimates correlate well with deep mutational scanning measurements. For most genes, synonymous mutations are nearly neutral, stop-codon mutations are deleterious, and amino acid mutations have a range of effects. However, some viral accessory proteins are under little to no selection. We provide interactive visualizations of effects of mutations to all SARS-CoV-2 proteins (https://jbloomlab.github.io/SARS2-mut-fitness/). The framework we describe is applicable to any virus for which the number of available sequences is sufficiently large that many independent occurrences of each neutral mutation are observed.

Insights into early evolutionary adaptations of the Akkermansia genus to the vertebrate gut

Akkermansia, a relevant mucin degrader from the vertebrate gut microbiota, is a member of the deeply branched Verrucomicrobiota, as well as the only known member of this phylum to be described as inhabitants of the gut. Only a few Akkermansia species have been officially described so far, although there is genomic evidence addressing the existence of more species-level variants for this genus. This niche specialization makes Akkermansia an interesting model for studying the evolution of microorganisms to their adaptation to the gastrointestinal tract environment, including which kind of functions were gained when the Akkermansia genus originated or how the evolutionary pressure functions over those genes. In order to gain more insight into Akkermansia adaptations to the gastrointestinal tract niche, we performed a phylogenomic analysis of 367 high-quality Akkermansia isolates and metagenome-assembled genomes, in addition to other members of Verrucomicrobiota. This work was focused on three aspects: the definition of Akkermansia genomic species clusters and the calculation and functional characterization of the pangenome for the most represented species; the evolutionary relationship between Akkermansia and their closest relatives from Verrucomicrobiota, defining the gene families which were gained or lost during the emergence of the last Akkermansia common ancestor (LAkkCA) and; the evaluation of the evolutionary pressure metrics for each relevant gene family of main Akkermansia species. This analysis found 25 Akkermansia genomic species clusters distributed in two main clades, divergent from their non-Akkermansia relatives. Pangenome analyses suggest that Akkermansia species have open pangenomes, and the gene gain/loss model indicates that genes associated with mucin degradation (both glycoside hydrolases and peptidases), (micro)aerobic metabolism, surface interaction, and adhesion were part of LAkkCA. Specifically, mucin degradation is a very ancestral innovation involved in the origin of Akkermansia. Horizontal gene transfer detection suggests that Akkermansia could receive genes mostly from unknown sources or from other Gram-negative gut bacteria. Evolutionary metrics suggest that Akkemansia species evolved differently, and even some conserved genes suffered different evolutionary pressures among clades. These results suggest a complex evolutionary landscape of the genus and indicate that mucin degradation could be an essential feature in Akkermansia evolution as a symbiotic species.

Phylogenomic analysis of the Porphyromonas gingivalis - Porphyromonas gulae duo: approaches to the origin of periodontitis

Porphyromonas gingivalis is an oral human pathogen associated with the onset and progression of periodontitis, a chronic immune-inflammatory disease characterized by the destruction of the teeth-supporting tissue. P. gingivalis belongs to the genus Porphyromonas, which is characterized by being composed of Gram-negative, asaccharolytic, non-spore-forming, non-motile, obligatory anaerobic species, inhabiting niches such as the oral cavity, urogenital tract, gastrointestinal tract and infected wound from different mammals including humans. Among the Porphyromonas genus, P. gingivalis stands out for its specificity in colonizing the human oral cavity and its keystone pathogen role in periodontitis pathogenesis. To understand the evolutionary process behind P. gingivalis in the context of the Pophyoromonas genus, in this study, we performed a comparative genomics study with publicly available Porphyromonas genomes, focused on four main objectives: (A) to confirm the phylogenetic position of P. gingivalis in the Porphyromonas genus by phylogenomic analysis; (B) the definition and comparison of the pangenomes of P. gingivalis and its relative P. gulae; and (C) the evaluation of the gene family gain/loss events during the divergence of P. gingivalis and P. gulae; (D) the evaluation of the evolutionary pressure (represented by the calculation of Tajima-D values and dN/dS ratios) comparing gene families of P. gingivalis and P. gulae. Our analysis found 84 high-quality assemblies representing P. gingivalis and 14 P. gulae strains (from a total of 233 Porphyromonas genomes). Phylogenomic analysis confirmed that P. gingivalis and P. gulae are highly related lineages, close to P. loveana. Both organisms harbored open pangenomes, with a strong core-to-accessory ratio for housekeeping genes and a negative ratio for unknown function genes. Our analyses also characterized the gene set differentiating P. gulae from P. gingivalis, mainly associated with unknown functions. Relevant virulence factors, such as the FimA, Mfa1, and the hemagglutinins, are conserved in P. gulae, P. gingivalis, and P. loveana, suggesting that the origin of those factors occurred previous to the P. gulae - P. gingivalis divergence. These results suggest an unexpected evolutionary relationship between the P. gulae - P. gingivalis duo and P. loveana, showing more clues about the origin of the role of those organisms in periodontitis.

The Multifaceted Role of Homologous Recombination in a Fastidious Bacterial Plant Pathogen

Homologous recombination plays a key function in the evolution of bacterial genomes. Within Xylella fastidiosa, an emerging plant pathogen with increasing host and geographic ranges, it has been suggested that homologous recombination facilitates host switching, speciation, and the development of virulence. We used 340 whole-genome sequences to study the relationship between inter- and intrasubspecific homologous recombination, random mutation, and natural selection across individual X. fastidiosa genes. Individual gene orthologs were identified and aligned, and a maximum likelihood (ML) gene tree was generated. Each gene alignment and tree pair were then used to calculate gene-wide and branch-specific r/m values (relative effect of recombination to mutation), gene-wide and branch-site nonsynonymous over synonymous substitution rates (dN/dS values; episodic selection), and branch length (as a proxy for mutation rate). The relationships between these variables were evaluated at the global level (i.e., for all genes among and within a subspecies), among specific functional classes (i.e., COGs), and between pangenome components (i.e., accessory versus core genes). Our analysis showed that r/m varied widely among genes as well as across X. fastidiosa subspecies. While r/m and dN/dS values were positively correlated in some instances (e.g., core genes in X. fastidiosa subsp. fastidiosa and both core and accessory genes in X. fastidiosa subsp. multiplex), low correlation coefficients suggested no clear biological significance. Overall, our results indicate that, in addition to its adaptive role in certain genes, homologous recombination acts as a homogenizing and a neutral force across phylogenetic clades, gene functional groups, and pangenome components. IMPORTANCE There is ample evidence that homologous recombination occurs frequently in the economically important plant pathogen Xylella fastidiosa. Homologous recombination has been known to occur among sympatric subspecies and is associated with host-switching events and virulence-linked genes. As a consequence, is it generally assumed that recombinant events in X. fastidiosa are adaptive. This mindset influences expectations of how homologous recombination acts as an evolutionary force as well as how management strategies for X. fastidiosa diseases are determined. Yet, homologous recombination plays roles beyond that of a source for diversification and adaptation. Homologous recombination can act as a DNA repair mechanism, as a means to facilitate nucleotide compositional change, as a homogenization mechanism within populations, or even as a neutral force. Here, we provide a first assessment of long-held beliefs regarding the general role of recombination in adaptation for X. fastidiosa. We evaluate gene-specific variations in homologous recombination rate across three X. fastidiosa subspecies and its relationship to other evolutionary forces (e.g., natural selection, mutation, etc.). These data were used to assess the role of homologous recombination in X. fastidiosa evolution.

Detection of positive selection acting on protein surfaces at the whole-genome scale in the human malaria parasite Plasmodium falciparum

The host-parasite evolutionary arms race is a fundamental process with medical implications. During this process, the host develops parasite resistance, and the parasite develops host immune evasion strategies. Thus, this process accelerates relevant protein evolution. This study test hypothesizes that proteins subject to sequence evolution structural constraints play a crucial role and that these constraints hinder the modification of such proteins in this process. These hypotheses were tested using Plasmodium falciparum model and evaluated protein structures predicted for the entire proteome by the AlphaFold method. Based on dN/dS test results and P. falciparum and P. reichenowi comparisons, the presented approach identified proteins subject to purifying selection acting on the whole sequence and buried residues (dN < dS) and positive selection on nonburied residues. Of the 26 proteins, some known antigens (ring-exported protein 3, RAP protein, erythrocyte binding antigen-140, and protein P47) targeted by the host immune system are promising vaccine candidates. The set also contained 11 enzymes, including FIKK kinase, which modifies host proteins. This set was compared with genes for which the dN/dS test suggested that positive selection acts on the whole gene (i.e., dN > dS). The present study found that such genes encode enzymes and antigenic vaccine candidates less frequently than genes for which evolution is not subject to selection constraints and positive selection acts on only exposed residues. The analysis was repeated comparing P. falciparum with P. alderi, which is more distantly related. The study discusses the potential implications of the presented methodology for rational vaccine design and the parasitology and evolutionary biology fields.

Cyanobacterial Blooms Are Not a Result of Positive Selection by Freshwater Eutrophication

Long-standing cyanobacterial harmful algal blooms (CyanoHABs) are known to result from synergistic interaction between elevated nutrients and superior ecophysiology of cyanobacteria. However, it remains to be determined whether CyanoHABs are a result of positive selection by eutrophic waters. To address this, we conducted molecular evolutionary analyses on the genomes of 9 bloom-forming cyanobacteria, combined with pangenomics and metatranscriptomics. The results showed no positive selection by water eutrophication. Instead, all homologous genes in the species are under strong purifying selection based on the ratio of divergence at nonsynonymous and synonymous sites (dN/dS) and phylogeny. The dN/dS < 0.85 (median = 0.3) for all homologous genes are similar between the genes in the pathways driving CyanoHABs and housekeeping functions. Phylogenetic support for non-positive selection comes from the mixed clustering of strains: strains of the same species from diverse geographic origins form the same clusters, while strains from the same origins form different clusters. Further support lies in the codon adaptation index (CAI) and single nucleotide polymorphism (SNP). The CAI ranged from 0.42 to 0.9 (mean = 0.75), which indicates high-level codon usage bias; the pathways for CyanoHABs and housekeeping functions showed a similar CAI. Interestingly, CAI was negatively correlated with gene expression in 3 metatranscriptomes. The numbers of SNPs were concentrated around 5 to 50. As the SNP number increases, the gene expression level decreases. These negative correlations agree with the population-level dN/dS and phylogeny in supporting purifying selection in bloom-forming cyanobacteria. In summary, superior ecophysiology appears to be acquired prior to water eutrophication. IMPORTANCE CyanoHABs are global environmental hazards, and their mechanisms of action are being intensively investigated. On an ecological scale, CyanoHABs are consequences of synergistic interactions between biological functions and elevated nutrients in eutrophic waters. On an evolutionary scale, one important question is how bloom-forming cyanobacteria acquire these superior biological functions. There are several possibilities, including adaptive evolution and horizontal gene transfer. Here, we explored the possibility of positive selection. We reasoned that there are two possible periods for cyanobacteria to acquire these functions: before the onset of water eutrophication or during water eutrophication. Either way, there should be molecular signatures in protein sequences for positive selection. Interestingly, we found no positive selection by water eutrophication, but strong purifying selection instead on nearly all the genes, suggesting these superior functions aiding CyanoHABs are acquired prior to water eutrophication.

The Fast and the Furriest: Investigating the Rate of Selection on Mammalian Toxins

The evolution of venom and the selection pressures that act on toxins have been increasingly researched within toxinology in the last two decades, in part due to the exceptionally high rates of diversifying selection observed in animal toxins. In 2015, Sungar and Moran proposed the 'two-speed' model of toxin evolution linking evolutionary age of a group to the rates of selection acting on toxins but due to a lack of data, mammals were not included as less than 30 species of venomous mammal have been recorded, represented by elusive species which produce small amounts of venom. Due to advances in genomics and transcriptomics, the availability of toxin sequences from venomous mammals has been increasing. Using branch- and site-specific selection models, we present the rates of both episodic and pervasive selection acting upon venomous mammal toxins as a group for the first time. We identified seven toxin groups present within venomous mammals, representing Chiroptera, Eulipotyphla and Monotremata: KLK1, Plasminogen Activator, Desmallipins, PACAP, CRiSP, Kunitz Domain One and Kunitz Domain Two. All but one group (KLK1) was identified by our results to be evolving under both episodic and pervasive diversifying selection with four toxin groups having sites that were implicated in the fitness of the animal by TreeSAAP (Selection on Amino Acid Properties). Our results suggest that venomous mammal ecology, behaviour or genomic evolution are the main drivers of selection, although evolutionary age may still be a factor. Our conclusion from these results indicates that mammalian toxins are following the two-speed model of selection, evolving predominately under diversifying selection, fitting in with other younger venomous taxa like snakes and cone snails-with high amounts of accumulating mutations, leading to more novel adaptions in their toxins.

Gene transcriptional profiles in gonads of Bacillus taxa (Phasmida) with different cytological mechanisms of automictic parthenogenesis

The evolution of automixis - i.e., meiotic parthenogenesis - requires several features, including ploidy restoration after meiosis and maintenance of fertility. Characterizing the relative contribution of novel versus pre-existing genes and the similarities in their expression and sequence evolution is fundamental to understand the evolution of reproductive novelties. Here we identify gonads-biased genes in two Bacillus automictic stick-insects and compare their expression profile and sequence evolution with a bisexual congeneric species. The two parthenogens restore ploidy through different cytological mechanisms: in Bacillus atticus, nuclei derived from the first meiotic division fuse to restore a diploid egg nucleus, while in Bacillus rossius, diploidization occurs in some cells of the haploid blastula through anaphase restitution. Parthenogens' gonads transcriptional program is found to be largely assembled from genes that were already present before the establishment of automixis. The three species transcriptional profiles largely reflect their phyletic relationships, yet we identify a shared core of genes with gonad-biased patterns of expression in parthenogens which are either male gonads-biased in the sexual species or are not differentially expressed there. At the sequence level, just a handful of gonads-biased genes were inferred to have undergone instances of positive selection exclusively in the parthenogen species. This work is the first to explore the molecular underpinnings of automixis in a comparative framework: it delineates how reproductive novelties can be sustained by genes whose origin precedes the establishment of the novelty itself and shows that different meiotic mechanisms of reproduction can be associated with a shared molecular ground plan.

Nonadaptive molecular evolution of plastome during the speciation of Actaea purpurea and its relatives

We have seen an explosive increase of plant plastid genome (plastome) sequences in the last decade, and the view that sequence variation in plastomes is maintained by the mutation-drift balance has been challenged by new evidence. Although comparative genomic and population-level studies provided us with evidence for positive evolution of plastid genes at both the macro- and micro-evolution levels, less studies have systematically investigated how plastomes have evolved during the speciation process. We here sequenced 13 plastomes of Actaea purpurea (P.K. Hsiao) J. Compton, and its closest relatives, and conducted a systematic survey of positive selection in their plastid genes using the McDonald-Kreitman test and codon-based methods using maximum likelihood to estimate the ratio of nonsynonymous to synonymous substitutions (ω) across a phylogeny. We found that during the speciation of A. purpurea and its relatives, all plastid genes evolved neutrally or were under purifying selection. Genome size, gene order, and number were highly conserved. Comparing to A. purpurea, plastomes of Actaea japonica and Actaea biternata had low genetic diversity, consistent with previous studies. Our work not only sheds important light on the evolutionary history of A. purpurea and its kin, but also on the evolution of plastomes during plant speciation.

Non-synonymous to synonymous substitutions suggest that orthologs tend to keep their functions, while paralogs are a source of functional novelty

Orthologs separate after lineages split from each other and paralogs after gene duplications. Thus, orthologs are expected to remain more functionally coherent across lineages, while paralogs have been proposed as a source of new functions. Because protein functional divergence follows from non-synonymous substitutions, we performed an analysis based on the ratio of non-synonymous to synonymous substitutions (dN/dS), as proxy for functional divergence. We used five working definitions of orthology, including reciprocal best hits (RBH), among other definitions based on network analyses and clustering. The results showed that orthologs, by all definitions tested, had values of dN/dS noticeably lower than those of paralogs, suggesting that orthologs generally tend to be more functionally stable than paralogs. The differences in dN/dS ratios remained suggesting the functional stability of orthologs after eliminating gene comparisons with potential problems, such as genes with high codon usage biases, low coverage of either of the aligned sequences, or sequences with very high similarities. Separation by percent identity of the encoded proteins showed that the differences between the dN/dS ratios of orthologs and paralogs were more evident at high sequence identity, less so as identity dropped. The last results suggest that the differences between dN/dS ratios were partially related to differences in protein identity. However, they also suggested that paralogs undergo functional divergence relatively early after duplication. Our analyses indicate that choosing orthologs as probably functionally coherent remains the right approach in comparative genomics.

Comparative Analysis of SARS-CoV-2 Variants of Concern, Including Omicron, Highlights Their Common and Distinctive Amino Acid Substitution Patterns, Especially at the Spike ORF

In order to gain a deeper understanding of the recently emerged and highly divergent Omicron variant of concern (VoC), a study of amino acid substitution (AAS) patterns was performed and compared with those of the other four successful variants of concern (Alpha, Beta, Gamma, Delta) and one closely related variant of interest (VoI-Lambda). The Spike ORF consistently emerges as an AAS hotspot in all six lineages, but in Omicron this enrichment is significantly higher. The progenitors of each of these VoC/VoI lineages underwent positive selection in the Spike ORF. However, once they were established, their Spike ORFs have been undergoing purifying selection, despite the application of global vaccination schemes from 2021 onwards. Our analyses reject the hypothesis that the heavily mutated receptor binding domain (RBD) of the Omicron Spike was introduced via recombination from another closely related Sarbecovirus. Thus, successive point mutations appear as the most parsimonious scenario. Intriguingly, in each of the six lineages, we observed a significant number of AAS wherein the new residue is not present at any homologous site among the other known Sarbecoviruses. Such AAS should be further investigated as potential adaptations to the human host. By studying the phylogenetic distribution of AAS shared between the six lineages, we observed that the Omicron (BA.1) lineage had the highest number (8/10) of recurrent mutations.

Faster Rates of Molecular Sequence Evolution in Reproduction-Related Genes and in Species with Hypodermic Sperm Morphologies

Sexual selection drives the evolution of many striking behaviors and morphologies and should leave signatures of selection at loci underlying these phenotypes. However, although loci thought to be under sexual selection often evolve rapidly, few studies have contrasted rates of molecular sequence evolution at such loci across lineages with different sexual selection contexts. Furthermore, work has focused on separate sexed animals, neglecting alternative sexual systems. We investigate rates of molecular sequence evolution in hermaphroditic flatworms of the genus Macrostomum. Specifically, we compare species that exhibit contrasting sperm morphologies, strongly associated with multiple convergent shifts in the mating strategy, reflecting different sexual selection contexts. Species donating and receiving sperm in every mating have sperm with bristles, likely to prevent sperm removal. Meanwhile, species that hypodermically inject sperm lack bristles, potentially as an adaptation to the environment experienced by hypodermic sperm. Combining functional annotations from the model, Macrostomum lignano, with transcriptomes from 93 congeners, we find genus-wide faster sequence evolution in reproduction-related versus ubiquitously expressed genes, consistent with stronger sexual selection on the former. Additionally, species with hypodermic sperm morphologies had elevated molecular sequence evolution, regardless of a gene's functional annotation. These genome-wide patterns suggest reduced selection efficiency following shifts to hypodermic mating, possibly due to higher selfing rates in these species. Moreover, we find little evidence for convergent amino acid changes across species. Our work not only shows that reproduction-related genes evolve rapidly also in hermaphroditic animals, but also that well-replicated contrasts of different sexual selection contexts can reveal underappreciated genome-wide effects.

BASE: A novel workflow to integrate nonubiquitous genes in comparative genomics analyses for selection

Inferring the selective forces that orthologous genes underwent across different lineages can help us understand the evolutionary processes that have shaped their extant diversity and the phenotypes they underlie. The most widespread metric to estimate the selection regimes of coding genes-across sites and phylogenies-is the ratio of nonsynonymous to synonymous substitutions (dN/dS, also known as ω). Nowadays, modern sequencing technologies and the large amount of already available sequence data allow the retrieval of thousands of orthologous genes across large numbers of species. Nonetheless, the tools available to explore selection regimes are not designed to automatically process all genes, and their practical usage is often restricted to the single-copy ones which are found across all species considered (i.e., ubiquitous genes). This approach limits the scale of the analysis to a fraction of single-copy genes, which can be as low as an order of magnitude in respect to those which are not consistently found in all species considered (i.e., nonubiquitous genes). Here, we present a workflow named BASE that-leveraging the CodeML framework-eases the inference and interpretation of gene selection regimes in the context of comparative genomics. Although a number of bioinformatics tools have already been developed to facilitate this kind of analyses, BASE is the first to be specifically designed to allow the integration of nonubiquitous genes in a straightforward and reproducible manner. The workflow-along with all relevant documentation-is available at github.com/for-giobbe/BASE.

Rates and patterns of molecular evolution in bryophyte genomes, with focus on complex thalloid liverworts, Marchantiopsida

Plants commonly referred to as "bryophytes" belong to three major lineages of non-vascular plants: the liverworts, the hornworts and the mosses. They are unique among land plants in having a dominant haploid generation and a short-lived diploid sporophytic generation. The dynamics of selection acting on a haploid genome differs from those acting on a diploid genome: new mutations are directly exposed to selection. The general aim of this paper is to investigate the diversification rateof bryophytes - measured as silent site substitution rate representing neutral evolution (mutation rate) and the nonsynonymous to synonymous substitution rate ratio (dN/dS) representing selective evolution - and compare it with earlier studies on vascular plants. Results show that the silent site substitution rate is lower for liverworts as compared to angiosperms, but not as low as for gymnosperms. The selection pressure, measured as dN/dS, isnot remarkably lower for bryophytes as compared to other diploid dominant plants as would be expected by the masking hypothesis, indicating that other factors are more important than ploidy.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Evolutionary dynamics of sex-biased genes expressed in cricket brains and gonads

Sex-biased gene expression, particularly sex-biased expression in the gonad, has been linked to rates of protein sequence evolution (nonsynonymous to synonymous substitutions, dN/dS) in animals. However, in insects, sex-biased expression studies remain centred on a few holometabolous species. Moreover, other major tissue types such as the brain remain underexplored. Here, we studied sex-biased gene expression and protein evolution in a hemimetabolous insect, the cricket Gryllus bimaculatus. We generated novel male and female RNA-seq data for two sexual tissue types, the gonad and somatic reproductive system, and for two core components of the nervous system, the brain and ventral nerve cord. From a genome-wide analysis, we report several core findings. Firstly, testis-biased genes had accelerated evolution, as compared to ovary-biased and unbiased genes, which was associated with positive selection events. Secondly, although sex-biased brain genes were much less common than for the gonad, they exhibited a striking tendency for rapid protein sequence evolution, an effect that was stronger for the female than male brain. Further, some sex-biased brain genes were linked to sexual functions and mating behaviours, which we suggest may have accelerated their evolution via sexual selection. Thirdly, a tendency for narrow cross-tissue expression breadth, suggesting low pleiotropy, was observed for sex-biased brain genes, suggesting relaxed purifying selection, which we speculate may allow enhanced freedom to evolve adaptive protein functional changes. The findings of rapid evolution of testis-biased genes and male and female-biased brain genes are discussed with respect to pleiotropy, positive selection and the mating biology of this cricket.

Weak selection on synonymous codons substantially inflates dN/dS estimates in bacteria

Synonymous codon substitutions are not always selectively neutral as revealed by several types of analyses, including studies of codon usage patterns among genes. We analyzed codon usage in 13 bacterial genomes sampled from across a large order of bacteria, Enterobacterales, and identified presumptively neutral and selected classes of synonymous substitutions. To estimate substitution rates, given a neutral/selected classification of synonymous substitutions, we developed a flexible [Formula: see text] substitution model that allows multiple classes of synonymous substitutions. Under this multiclass synonymous substitution (MSS) model, the denominator of [Formula: see text] includes only the strictly neutral class of synonymous substitutions. On average, the value of [Formula: see text] under the MSS model was 80% of that under the standard codon model in which all synonymous substitutions are assumed to be neutral. The indication is that conventional [Formula: see text] analyses overestimate these values and thus overestimate the frequency of positive diversifying selection and underestimate the strength of purifying selection. To quantify the strength of selection necessary to explain this reduction, we developed a model of selected compensatory codon substitutions. The reduction in synonymous substitution rate, and thus the contribution that selection makes to codon bias variation among genes, can be adequately explained by very weak selection, with a mean product of population size and selection coefficient, [Formula: see text].

GenomegaMap: Within-Species Genome-Wide dN/dS Estimation from over 10,000 Genomes

The dN/dS ratio provides evidence of adaptation or functional constraint in protein-coding genes by quantifying the relative excess or deficit of amino acid-replacing versus silent nucleotide variation. Inexpensive sequencing promises a better understanding of parameters, such as dN/dS, but analyzing very large data sets poses a major statistical challenge. Here, I introduce genomegaMap for estimating within-species genome-wide variation in dN/dS, and I apply it to 3,979 genes across 10,209 tuberculosis genomes to characterize the selection pressures shaping this global pathogen. GenomegaMap is a phylogeny-free method that addresses two major problems with existing approaches: 1) It is fast no matter how large the sample size and 2) it is robust to recombination, which causes phylogenetic methods to report artefactual signals of adaptation. GenomegaMap uses population genetics theory to approximate the distribution of allele frequencies under general, parent-dependent mutation models. Coalescent simulations show that substitution parameters are well estimated even when genomegaMap’s simplifying assumption of independence among sites is violated. I demonstrate the ability of genomegaMap to detect genuine signatures of selection at antimicrobial resistance-conferring substitutions in Mycobacterium tuberculosis and describe a novel signature of selection in the cold-shock DEAD-box protein A gene deaD/csdA. The genomegaMap approach helps accelerate the exploitation of big data for gaining new insights into evolution within species.

Comparative Gene Expression Analysis Reveals Mechanism of Pinus contorta Response to the Fungal Pathogen Dothistroma septosporum

Many conifers have distributions that span wide ranges in both biotic and abiotic conditions, but the basis of response to biotic stress has received much less attention than response to abiotic stress. In this study, we investigated the gene expression response of lodgepole pine (Pinus contorta) to attack by the fungal pathogen Dothistroma septosporum, which causes Dothistroma needle blight, a disease that has caused severe climate-related outbreaks in northwestern British Columbia. We inoculated tolerant and susceptible pines with two D. septosporum isolates and analyzed the differentially expressed genes (DEGs), differential exon usage, and coexpressed gene modules using RNA-sequencing data. We found a rapid and strong transcriptomic response in tolerant lodgepole pine samples inoculated with one D. septosporum isolate, and a late and weak response in susceptible samples inoculated with another isolate. We mapped 43 of the DEG- or gene module-identified genes to the reference plant-pathogen interaction pathway deposited in the Kyoto Encyclopedia of Genes and Genomes database. These genes are present in PAMP-triggered and effector-triggered immunity pathways. Genes comprising pathways and gene modules had signatures of strong selective constraint, while the highly expressed genes in tolerant samples appear to have been favored by selection to counterattack the pathogen. We identified candidate resistance genes that may respond to D. septosporum effectors. Taken together, our results show that gene expression response to D. septosporum infection in lodgepole pine varies both among tree genotypes and pathogen strains and involves both known candidate genes and a number of genes with previously unknown functions.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Mitogenomics of the endemic Ethiopian rats: looking for footprints of adaptive evolution in sky islands

Organisms living in high altitude must adapt to environmental conditions with hypoxia and low temperature, e.g. by changes in the structure and function of proteins associated with oxidative phosphorylation in mitochondria. Here we analysed the signs of adaptive evolution in 27 mitogenomes of endemic Ethiopian rats (Stenocephalemys), where individual species adapted to different elevation. Significant signals of positive selection were detected in 10 of the 13 mitochondrial protein-coding genes, with a majority of functional substitutions in the NADH dehydrogenase complex. Higher frequency of positively selected sites was found in phylogenetic lineages corresponding to Afroalpine specialists.

Nonsynonymous Polymorphism Counts in Bacterial Genomes: a Comparative Examination

Genomic data reveal single-nucleotide polymorphisms (SNPs) that may carry information about the evolutionary history of bacteria. However, it remains unclear what inferences about selection can be made from genomic SNP data. Bacterial species are often sampled during epidemic outbreaks or within hosts during the course of chronic infections. SNPs obtained from genomic analysis of these data are not necessarily fixed. Treating them as fixed during analysis by using measures such as the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS) may lead to incorrect inferences about the strength and direction of selection. In this study, we consider data from a range of whole-genome sequencing studies of bacterial pathogens and explore patterns of nonsynonymous variation to assess whether evidence of selection can be identified by investigating SNP counts alone across multiple WGS studies. We visualize these SNP data in ways that highlight their relationship to neutral baseline expectations. These neutral expectations are based on a simple model of mutation, from which we simulate SNP accumulation to investigate how SNP counts are distributed under alternative assumptions about positive and negative selection. We compare these patterns with empirical SNP data and illustrate the general difficulty of detecting positive selection from SNP data. Finally, we consider whether SNP counts observed at the between-host population level differ from those observed at the within-host level and find some evidence that suggests that dynamics across these two scales are driven by different underlying processes.IMPORTANCE Identifying selection from SNP data obtained from whole-genome sequencing studies is challenging. Some current measures used to identify and quantify selection acting on genomes rely on fixed differences; thus, these are inappropriate for SNP data where variants are not fixed. With the increase in whole-genome sequencing studies, it is important to consider SNP data in the context of evolutionary processes. How SNPs are counted and analyzed can help in understanding mutation accumulation and trajectories of strains. We developed a tool for identifying possible evidence of selection and for comparative analysis with other SNP data. We propose a model that provides a rule-of-thumb guideline and two new visualization techniques that can be used to interpret and compare SNP data. We quantify the expected proportion of nonsynonymous SNPs in coding regions under neutrality and demonstrate its use in identifying evidence of positive and negative selection from simulations and empirical data.

Comparative genome analyses of Mycobacteroides immunogenum reveals two potential novel subspecies

Mycobacteroides immunogenum is an emerging opportunistic pathogen implicated in nosocomial infections. Comparative genome analyses may provide better insights into its genomic structure, functions and evolution. The present analysis showed that M. immunogenum has an open pan-genome. Approximately 36.8% of putative virulence genes were identified in the accessory regions of M. immunogenum. Phylogenetic analyses revealed two potential novel subspecies of M. immunogenum, supported by evidence from ANIb (average nucleotide identity using blast) and GGDC (Genome to Genome Distance Calculator) analyses. We identified 74 genomic islands (GIs) in Subspecies 1 and 23 GIs in Subspecies 2. All Subspecies 2-harboured GIs were not found in Subspecies 1, indicating that they might have been acquired by Subspecies 2 after their divergence. Subspecies 2 has more defence genes than Subspecies 1, suggesting that it might be more resistant to the insertion of foreign DNA and probably explaining why Subspecies 2 has fewer GIs. Positive selection analysis suggest that M. immunogenum has a lower selection pressure compared to non-pathogenic mycobacteria. Thirteen genes were positively selected and many were involved in virulence.

The Organelle Genomes in the Photosynthetic Red Algal Parasite Pterocladiophila hemisphaerica (Florideophyceae, Rhodophyta) Have Elevated Substitution Rates and Extreme Gene Loss in the Plastid Genome

Comparative organelle genome studies of parasites can highlight genetic changes that occur during the transition from a free-living to a parasitic state. Our study focuses on a poorly studied group of red algal parasites, which are often closely related to their red algal hosts and from which they presumably evolved. Most of these parasites are pigmented and some show photosynthetic capacity. Here, we assembled and annotated the complete organelle genomes of the photosynthetic red algal parasite, Pterocladiophila hemisphaerica. The plastid genome is the smallest known red algal plastid genome at 68,701 bp. The plastid genome has many genes missing, including all photosynthesis-related genes. In contrast, the mitochondrial genome is similar in architecture to that of other free-living red algae. Both organelle genomes show elevated mutation rates and significant changes in patterns of selection, measured as dN/dS ratios. This caused phylogenetic analyses, even of multiple aligned proteins, to be unresolved or give contradictory relationships. Full plastid datasets interfered by selected best gene evolution models showed the supported relationship of P. hemisphaerica within the Ceramiales, but the parasite was grouped with support as sister to the Gracilariales when interfered under the GHOST model. Nuclear rDNA showed a supported grouping of the parasite within a clade containing several red algal orders including the Gelidiales. This photosynthetic parasite, which is unable to photosynthesize with its own plastid due to the total loss of all photosynthesis genes, raises intriguing questions on parasite-host organelle genome capabilities and interactions.

Comparative Transcriptomics Reveals Gene Families Associated with Predatory Behavior in Photuris femme fatale Fireflies

Identifying the basis of phenotypic variation is a key objective of genetics. This work has been mostly limited to model systems with a plethora of genetic manipulation and functional characterization tools. With the development of high-throughput sequencing and new computational tools, it is possible to identify candidate genes related to phenotypic variation in non-model organisms. Fireflies are excellent for studying phenotypic variation because of their diverse and well-characterized behaviors. Most adult fireflies emit a single mating flash pattern and do not eat. In contrast, adult females of many species in the genus Photuris employ multiple flash patterns and prey upon mate-seeking males of other firefly species. To investigate the genetic basis for this variation, we used comparative transcriptomics to identify positively selected genes between a predatory firefly, Photuris sp., and a non-predatory relative, Photuris frontalis, controlling for genes generally under selection in fireflies by comparing to a Photinus firefly. Nine gene families were identified under positive selection in the predatory versus non-predatory Photuris comparison, including genes involved in digestion, detoxification, vision, reproduction, and neural processes. These results generate intriguing hypotheses about the genetic basis for insect behavior and highlight the utility of comparative transcriptomic tools to investigate complex behaviors in non-model systems.

Bioinformatic Detection of Positive Selection Pressure in Plant Pathogens: The Neutral Theory of Molecular Sequence Evolution in Action

The genomes of plant pathogenic fungi and oomycetes are often exposed to strong positive selection pressure. During speciation, shifts in host range and preference can lead to major adaptive changes. Furthermore, evolution of total host resistance to most isolates can force rapid evolutionary changes in host-specific pathogens. Crop pathogens are subjected to particularly intense selective pressures from monocultures and fungicides. Detection of the footprints of positive selection in plant pathogen genomes is a worthwhile endeavor as it aids understanding of the fundamental biology of these important organisms. There are two main classes of test for detection of positively selected alleles. Tests based on the ratio of non-synonymous to synonymous substitutions per site detect the footprints of multiple fixation events between divergent lineages. Thus, they are well-suited to the study of ancient adaptation events spanning speciations. On the other hand, tests that scan genomes for local fluctuations in allelic diversity within populations are suitable for detection of recent positive selection in populations. In this review, I briefly describe some of the more widely used tests of positive selection and the theory underlying them. I then discuss various examples of their application to plant pathogen genomes, emphasizing the types of genes that are associated with signatures of positive selection. I conclude with a discussion of the practicality of such tests for identification of pathogen genes of interest and the important features of pathogen ecology that must be taken into account for accurate interpretation.

Patterns of gene evolution following duplications and speciations in vertebrates

BACKGROUND

Eukaryotic genes typically form independent evolutionary lineages through either speciation or gene duplication events. Generally, gene copies resulting from speciation events (orthologs) are expected to maintain similarity over time with regard to sequence, structure and function. After a duplication event, however, resulting gene copies (paralogs) may experience a broader set of possible fates, including partial (subfunctionalization) or complete loss of function, as well as gain of new function (neofunctionalization). This assumption, known as the Ortholog Conjecture, is prevalent throughout molecular biology and notably plays an important role in many functional annotation methods. Unfortunately, studies that explicitly compare evolutionary processes between speciation and duplication events are rare and conflicting.

METHODS

To provide an empirical assessment of ortholog/paralog evolution, we estimated ratios of nonsynonymous to synonymous substitutions (ω = dN/dS) for 251,044 lineages in 6,244 gene trees across 77 vertebrate taxa.

RESULTS

Overall, we found ω to be more similar between lineages descended from speciation events (p < 0.001) than lineages descended from duplication events, providing strong support for the Ortholog Conjecture. The asymmetry in ω following duplication events appears to be largely driven by an increase along one of the paralogous lineages, while the other remains similar to the parent. This trend is commonly associated with neofunctionalization, suggesting that gene duplication is a significant mechanism for generating novel gene functions.

Molecular convergence and positive selection associated with the evolution of symbiont transmission mode in stony corals

Heritable symbioses have been critical for the evolution of life. The genetic consequences of evolving a heritable symbiosis from the perspective of the symbiont are well established, but concomitant changes in the host remain unresolved. In stony corals, heritable, vertical transmission has evolved repeatedly, providing a unique opportunity to investigate the genomic basis of this complex trait. We conducted a comparative analysis of 25 coral transcriptomes to identify orthologous genes exhibiting signatures of positive selection and convergent amino acid substitutions in vertically transmitting lineages. The frequency of convergence events tends to be higher among vertically transmitting lineages, consistent with the proposed role of selection in driving the evolution of convergent transmission mode phenotypes. Of 10 774 orthologous genes, 403 exhibited at least one molecular convergence event and evidence of positive selection in at least one vertically transmitting lineage. Functional enrichments among these top candidate genes include processes previously implicated in symbiosis including endocytosis, immune response, cytoskeletal protein binding and cytoplasmic membrane-bounded vesicles. Finally, several novel candidates were identified among 100 genes showing evidence of positive selection at the particular convergence event, highlighting the value of our approach for generating new insight into host mechanisms associated with the evolution of heritable symbioses.

Absence of a Faster-X Effect in Beetles (Tribolium, Coleoptera)

The faster-X effect, namely the rapid evolution of protein-coding genes on the X chromosome, has been widely reported in metazoans. However, the prevalence of this phenomenon across diverse systems and its potential causes remain largely unresolved. Analysis of sex-biased genes may elucidate its possible mechanisms: for example, in systems with X/Y males a more pronounced faster-X effect in male-biased genes than in female-biased or unbiased genes may suggest fixation of recessive beneficial mutations rather than genetic drift. Further, theory predicts that the faster-X effect should be promoted by X chromosome dosage compensation. Here, we asked whether we could detect a faster-X effect in genes of the beetle Tribolium castaneum (and T. freemani orthologs), which has X/Y sex-determination and heterogametic males. Our comparison of protein sequence divergence (dN/dS) on the X chromosome vs. autosomes indicated a rarely observed absence of a faster-X effect in this organism. Further, analyses of sex-biased gene expression revealed that the X chromosome was particularly highly enriched for ovary-biased genes, which evolved slowly. In addition, an evaluation of male X chromosome dosage compensation in the gonads and in non-gonadal somatic tissues indicated a striking lack of compensation in the testis. This under-expression in testis may limit fixation of recessive beneficial X-linked mutations in genes transcribed in these male sex organs. Taken together, these beetles provide an example of the absence of a faster-X effect on protein evolution in a metazoan, that may result from two plausible factors, strong constraint on abundant X-linked ovary-biased genes and a lack of gonadal dosage compensation.

Polymorphism Data Assist Estimation of the Nonsynonymous over Synonymous Fixation Rate Ratio ω for Closely Related Species

The ratio of nonsynonymous over synonymous sequence divergence, dN/dS, is a widely used estimate of the nonsynonymous over synonymous fixation rate ratio ω, which measures the extent to which natural selection modulates protein sequence evolution. Its computation is based on a phylogenetic approach and computes sequence divergence of protein-coding DNA between species, traditionally using a single representative DNA sequence per species. This approach ignores the presence of polymorphisms and relies on the indirect assumption that new mutations fix instantaneously, an assumption which is generally violated and reasonable only for distantly related species. The violation of the underlying assumption leads to a time-dependence of sequence divergence, and biased estimates of ω in particular for closely related species, where the contribution of ancestral and lineage-specific polymorphisms to sequence divergence is substantial. We here use a time-dependent Poisson random field model to derive an analytical expression of dN/dS as a function of divergence time and sample size. We then extend our framework to the estimation of the proportion of adaptive protein evolution α. This mathematical treatment enables us to show that the joint usage of polymorphism and divergence data can assist the inference of selection for closely related species. Moreover, our analytical results provide the basis for a protocol for the estimation of ω and α for closely related species. We illustrate the performance of this protocol by studying a population data set of four corvid species, which involves the estimation of ω and α at different time-scales and for several choices of sample sizes.

Relaxation of selective constraints shapes variation of toll-like receptors in a colonial waterbird, the black-headed gull

Nonspecific innate immune response is activated by toll-like receptors (TLRs), which recognize conserved molecular motifs characteristic for a broad spectrum of pathogens. In this study, we examined nucleotide substitution patterns and allelic diversity at five TLR genes in a wild nonpasserine bird, the black-headed gull Chroicocephalus ridibundus. We hypothesized that balancing selection can maintain high allelic diversity of TLR genes in the black-headed gull because of its ecological characteristics, coloniality, and migratoriness, which are associated with increased exposure and transmission of pathogens. Although we found moderately high levels of sequence polymorphism (8-49 haplotypes retrieved per locus within a sample of 60 individuals), most of these haplotypes were recorded at low frequencies within our study population. At the same time, we found no convincing evidence for the role of balancing selection in the maintenance of this variation (Tajima's D < 0.5), and sites with a significant excess of nonsynonymous mutations (dN/dS > 1) were recorded only at two loci (TLR5 and TLR7). This pattern is consistent with relaxation of selective constraints, where most mutations are slightly deleterious and usually removed by purifying selection. No differences in the diversity and nucleotide substitution rates were found between endosomal loci responsible for viral RNA sensing and loci responsible for the recognition of extracellular pathogens. Our study provides the first information on evolutionary mechanisms shaping polymorphism of TLRs in a species from Lari suborder (gulls and allies) and suggests that TLR genes may be poorly responsive to ecological and life-history characteristics of hosts.

Analysis of tomato spotted wilt virus RNA-dependent RNA polymerase adaptative evolution and constrained domains using homology protein structure modelling

Tomato spotted wilt virus (TSWV; genus Orthotospovirus, family Tospoviridae) has a huge impact on a large range of plants worldwide. In this study, we determined the sequence of the large (L) RNA segment that encodes the RNA-dependent RNA polymerase (RdRp) from a TSWV isolate (LYE51) collected in the south of France. Analysis of the phylogenetic relationships of TSWV-LYE51 with other TSWV isolates shows that it is closely related to other European isolates. A 3D model of TSWV-LYE51 RdRp was built by homology with the RdRp structure of the La Crosse virus (genus Orthobunyavirus, family Peribunyaviridae). Finally, an analysis of positive and negative selection was carried out on 30 TSWV full-length RNA L sequences and compared with the phylogeny and the protein structure data. We showed that the seven codons that are under positive selection are distributed all along the RdRp gene. By contrast, the codons associated with negative selection are especially concentrated in three highly constrained domains: the endonuclease in charge of the cap-snatching mechanism, the thumb domain and the mid domain. Those three domains could constitute good candidates to look for host targets on which genetic resistance by loss of susceptibility could be developed.

Genetic assessment of the internal transcribed spacer region (ITS1.2) in Mangifera indica L. landraces

Mango (Mangifera indica) is one of the most important tropical fruits in the world. Twenty-two genotypes of native mangoes from different regions of southern Iran (Hormozgan and Kerman) were collected and analyzed for the ribosomal genes. GC content was found to be 55.5%. Fu and Li's D* test statistic (0.437), Fu and Li's F* test statistic (0.500) and Tajima's D (1.801) were positive and nonsignificant. A total of 769 positions were identified (319 with insertion or deletion including 250 polymorphic and 69 monomorphic loci; 450 loci without any insertion or deletion including 35 Singletons and 22 haplotypes). Nucleotide diversity of 0.309 and a high genetic differentiation including Chi square of 79.8; P value of 0.3605 and df value of 76 was observed among mango genotypes studied. The numerical value of the ratio dN/dS (0.45) indicated a pure selection in the examined gene and the absence of any key changes. Cluster analysis differentiated the mango used in this research (M. indica L.) into two genotypes but could not differentiate their geographical locations. The results of this study indicated that a high genetic distance exists between HajiGholam (Manojan) and Arbabi (Rodan) genotypes and showed higher genetic diversity in mango of Rodan region. Results of present study suggested that for successful breeding, the genotypes of Rodan region mango especially Arbabi mango can be used as a gene donor and ITS can be a suitable tool for genetic evaluations of inter and intra species.

Data for positive selection test and co-evolutionary analysis on mammalian cereblon

Cereblon (CRBN) is a substrate recognition subunit of the CRL4 E3 ubiquitin ligase complex, directly binding to specific substrates for poly-ubiquitination followed by proteasome-dependent degradation of proteins. Cellular CRBN is responsible for energy metabolism, ion-channel activation, and cellular stress response through binding to proteins related to the respective pathways. As CRBN binds to various proteins, the selective pressure at the interacting surface is expected to result in functional divergence. Here, we present two mammalian CRBN datasets of molecular evolutionary analyses. (1) The multiple sequence alignment data shows that positive selection occurred, determined with a dN/dS calculation. (2) Data on co-evolutionary analysis between vertebrate CRBN and related proteins are represented by calculating the correlation coefficient based on the comparison of phylogenetic trees. Co-evolutionary analysis shows the similarity of evolutionary traits of two proteins. Further molecular, functional interpretation of these analyses is explained in ‘Positive selection of Cereblon modified function including its E3 Ubiquitin Ligase activity and binding efficiency with AMPK’ (W. Onodera, T. Asahi, N. Sawamura, Positive selection of cereblon modified function including its E3 ubiquitin ligase activity and binding efficiency with AMPK. Mol Phylogenet Evol. (2019) 135:78-85. [1]).

Characterization of the melanopsin gene (Opn4x) of diurnal and nocturnal snakes

BACKGROUND

A number of non-visual responses to light in vertebrates, such as circadian rhythm control and pupillary light reflex, are mediated by melanopsins, G-protein coupled membrane receptors, conjugated to a retinal chromophore. In non-mammalian vertebrates, melanopsin expression is variable within the retina and extra-ocular tissues. Two paralog melanopsin genes were classified in vertebrates, Opn4x and Opn4m. Snakes are highly diversified vertebrates with a wide range of daily activity patterns, which raises questions about differences in structure, function and expression pattern of their melanopsin genes. In this study, we analyzed the melanopsin genes expressed in the retinas of 18 snake species from three families (Viperidae, Elapidae, and Colubridae), and also investigated extra-retinal tissue expression.

RESULTS

Phylogenetic analysis revealed that the amplified gene belongs to the Opn4x group, and no expression of the Opn4m was found. The same paralog is expressed in the iris, but no extra-ocular expression was detected. Molecular evolutionary analysis indicated that melanopsins are evolving primarily under strong purifying selection, although lower evolutionary constraint was detected in snake lineages (ω = 0.2), compared to non-snake Opn4x and Opn4m (ω = 0.1). Statistical analysis of selective constraint suggests that snake phylogenetic relationships have driven stronger effects on melanopsin evolution, than the species activity pattern. In situ hybridization revealed the presence of melanopsin within cells in the outer and inner nuclear layers, in the ganglion cell layer, and intense labeling in the optic nerve.

CONCLUSIONS

The loss of the Opn4m gene and extra-ocular photosensitive tissues in snakes may be associated with a prolonged nocturnal/mesopic bottleneck in the early history of snake evolution. The presence of melanopsin-containing cells in all retinal nuclear layers indicates a globally photosensitive retina, and the expression in classic photoreceptor cells suggest a regionalized co-expression of melanopsin and visual opsins.

Genetic evolution and codon usage analysis of NKX-2.5 gene governing heart development in some mammals

NKX-2.5 gene is responsible for cardiac development and its targeted disruption apprehends cardiac development at the linear heart tube stage. Bioinformatic analysis was employed to investigate the codon usage pattern and dN/dS of mammalian NKX-2.5 gene. The relative synonymous codon usage analysis revealed variation in codon usage and two synonymous codons namely ATA (Ile) and GTA (Val) were absent in NKX-2.5 gene across selected mammalian species suggesting that these two codons were possibly selected against during evolution. Parity rule 2 analysis of two and four fold amino acids showed CT bias whereas six-fold amino acids revealed GA bias. Neutrality analysis suggests that selection played a prominent role while mutation had a minor role. The dN/dS analysis suggests synonymous substitution played a significant role and it negatively correlated with p-distance of the gene. Purifying natural selection played a dominant role in the genetic evolution of NKX-2.5 gene in mammals.

N-glycoproteins exhibit a positive expression level-evolutionary rate correlation

The different proteins of any proteome evolve at enormously different rates. One of the primary factors influencing rates of protein evolution is expression level, with highly expressed proteins tending to evolve at slow rates. This phenomenon, known as the expression level-evolutionary rate (E-R) anticorrelation, has been attributed to the abundance-dependent deleterious effects of misfolding or misinteraction. We have recently shown that secreted proteins either lack an E-R anticorrelation or exhibit a significantly reduced E-R anticorrelation. This effect may be due to the strict quality control to which secreted proteins are subject in the endoplasmic reticulum (which is expected to reduce the rate of misfolding and its deleterious effects) or to their extracellular location (expected to reduce the rate of misinteraction and its deleterious effects). Among secreted proteins, N-glycosylated ones are under particularly strong quality control. Here, we investigate how N-linked glycosylation affects the E-R anticorrelation. Strikingly, we observe a positive E-R correlation among N-glycosylated proteins. That is, N-glycoproteins that are highly expressed evolve at faster rates than lowly expressed N-glycoproteins, in contrast to what is observed among intracellular proteins.

Positive selection of cereblon modified function including its E3 ubiquitin ligase activity and binding efficiency with AMPK

Cereblon (CRBN) is a substrate receptor for an E3 ubiquitin ligase that directly binds to target proteins resulting in cellular activities, such as energy metabolism, membrane potential regulation, and transcription factor degradation. Genetic mutations in human CRBN lead to intellectual disabilities. In addition, it draws pathological attention because direct binding with immunomodulatory drugs can cure multiple myeloma (MM) and lymphocytic leukemia. To further explore the function of CRBN, we focused on its molecular evolution. Since CRBN interacts directly with its substrates and is widely conserved in vertebrates, evolutionary study to identify the selective pressure on CRBN that occur during CRBN-substrate interaction is an effective approach to search for a novel active site. Using mammalian CRBN sequences, dN/dS analysis was conducted to detect positive selection. By multiple sequence alignment we found that the residue at position 366 was under positive selection. This residue is present in the substrate-binding domain of CRBN. Most mammals harbor cysteine at position 366, whereas rodents and chiroptera have serine at this site. Subsequently, we constructed a C366S human CRBN to confirm the potential of positive selection. Auto-ubiquitination activity occurs in E3 ubiquitin ligases, including CRBN, and increased in C366S CRBN, which lead to the conclusion that E3 ubiquitin ligase activity may have changed over the course of mammalian evolution. Furthermore, binding with AMP-activated protein kinase was augmented when the substitution was present, which is supported by coevolution analysis. These results suggest that the molecular evolution of CRBN occurred through codon-based positive selection, providing a new approach to investigate CRBN function.

A Single Adaptive Mutation in Sodium Taurocholate Cotransporting Polypeptide Induced by Hepadnaviruses Determines Virus Species Specificity

Hepatitis B virus (HBV) and its hepadnavirus relatives infect a wide range of vertebrates, from fish to human. Hepadnaviruses and their hosts have a long history of acquiring adaptive mutations. However, there are no reports providing direct molecular evidence for such a coevolutionary "arms race" between hepadnaviruses and their hosts. Here, we present evidence suggesting that the adaptive evolution of the sodium taurocholate cotransporting polypeptide (NTCP), an HBV receptor, has been influenced by virus infection. Evolutionary analysis of the NTCP-encoding genes from 20 mammals showed that most NTCP residues are highly conserved among species, exhibiting evolution under negative selection (dN/dS ratio [ratio of nonsynonymous to synonymous evolutionary changes] of <1); this observation implies that the evolution of NTCP is restricted by maintaining its original protein function. However, 0.7% of NTCP amino acid residues exhibit rapid evolution under positive selection (dN/dS ratio of >1). Notably, a substitution at amino acid (aa) 158, a positively selected residue, converting the human NTCP to a monkey-type sequence abrogated the capacity to support HBV infection; conversely, a substitution at this residue converting the monkey Ntcp to the human sequence was sufficient to confer HBV susceptibility. Together, these observations suggested a close association of the aa 158 positive selection with the pressure by virus infection. Moreover, the aa 158 sequence determined attachment of the HBV envelope protein to the host cell, demonstrating the mechanism whereby HBV infection would create positive selection at this NTCP residue. In summary, we provide the first evidence in agreement with the function of hepadnavirus as a driver for inducing adaptive mutation in host receptor.IMPORTANCE HBV and its hepadnavirus relatives infect a wide range of vertebrates, with a long infectious history (hundreds of millions of years). Such a long history generally allows adaptive mutations in hosts to escape from infection while simultaneously allowing adaptive mutations in viruses to overcome host barriers. However, there is no published molecular evidence for such a coevolutionary arms race between hepadnaviruses and hosts. In the present study, we performed coevolutionary phylogenetic analysis between hepadnaviruses and the sodium taurocholate cotransporting polypeptide (NTCP), an HBV receptor, combined with virological experimental assays for investigating the biological significance of NTCP sequence variation. Our data provide the first molecular evidence supporting that HBV-related hepadnaviruses drive adaptive evolution in the NTCP sequence, including a mechanistic explanation of how NTCP mutations determine host viral susceptibility. Our novel insights enhance our understanding of how hepadnaviruses evolved with their hosts, permitting the acquisition of strong species specificity.

Looking for Darwin in Genomic Sequences: Validity and Success Depends on the Relationship Between Model and Data

Codon substitution models (CSMs) are commonly used to infer the history of natural section for a set of protein-coding sequences, often with the explicit goal of detecting the signature of positive Darwinian selection. However, the validity and success of CSMs used in conjunction with the maximum likelihood (ML) framework is sometimes challenged with claims that the approach might too often support false conclusions. In this chapter, we use a case study approach to identify four legitimate statistical difficulties associated with inference of evolutionary events using CSMs. These include: (1) model misspecification, (2) low information content, (3) the confounding of processes, and (4) phenomenological load, or PL. While past criticisms of CSMs can be connected to these issues, the historical critiques were often misdirected, or overstated, because they failed to recognize that the success of any model-based approach depends on the relationship between model and data. Here, we explore this relationship and provide a candid assessment of the limitations of CSMs to extract historical information from extant sequences. To aid in this assessment, we provide a brief overview of: (1) a more realistic way of thinking about the process of codon evolution framed in terms of population genetic parameters, and (2) a novel presentation of the ML statistical framework. We then divide the development of CSMs into two broad phases of scientific activity and show that the latter phase is characterized by increases in model complexity that can sometimes negatively impact inference of evolutionary mechanisms. Such problems are not yet widely appreciated by the users of CSMs. These problems can be avoided by using a model that is appropriate for the data; but, understanding the relationship between the data and a fitted model is a difficult task. We argue that the only way to properly understand that relationship is to perform in silico experiments using a generating process that can mimic the data as closely as possible. The mutation-selection modeling framework (MutSel) is presented as the basis of such a generating process. We contend that if complex CSMs continue to be developed for testing explicit mechanistic hypotheses, then additional analyses such as those described in here (e.g., penalized LRTs and estimation of PL) will need to be applied alongside the more traditional inferential methods.

Protocols for the Molecular Evolutionary Analysis of Membrane Protein Gene Duplicates

Gene duplication is an important process in the evolution of gene content in eukaryotic genomes. Understanding when gene duplicates contribute new molecular functions to genomes through molecular adaptation is one important goal in comparative genomics. In large gene families, however, characterizing adaptation and neofunctionalization across species is challenging, as models have traditionally quantified the timing of duplications without considering underlying gene trees. This protocol combines multiple approaches to detect adaptation in protein duplicates at a phylogenetic scale. We include a description of models for gene tree-species tree reconciliation that enable different types of inference, as well as a practical guide to their use. Although simulation-based approaches successfully detect shifts in the rate of duplication/retention, the conflation between the duplication and retention processes, the distinct trajectories of duplicates under non-, sub-, and neofunctionalization, as well as dosage effects offer hitherto unexplored analytical avenues. We introduce mathematical descriptions of these probabilities and offer a road map to computational implementation whose starting point is parsimony reconciliation. Sequence evolution information based on the ratio of nonsynonymous to synonymous nucleotide substitution rates (dN/dS) can be combined with duplicate survival probabilities to better predict the emergence of new molecular functions in retained duplicates. Together, these methods enable characterization of potentially adaptive candidate duplicates whose neofunctionalization may contribute to phenotypic divergence across species.

Evolutionary divergence of mitochondrial genomes in two Tetranychus species distributed across different climates

There is increasing evidence that mitochondrial genomes (mitogenomes) can be under selection, whereas the selective regimes shaping mitogenome evolution remain largely unclear. To test for mitogenome evolution in relation to the climate adaptation, we explored mtDNA variation in two spider mite (Tetranychus) species that distribute across different climates. We sequenced 26 complete mitogenomes of Tetranychus truncates, which occurs in both warm and cold regions, and nine complete mitogenomes of Tetranychus pueraricola, which is restricted to warm regions. Patterns of evolution in the two species' mitogenomes were compared through a series of d_N /d_S methods and physicochemical profiles of amino acid replacements. We found that: (1) the mitogenomes of both species were under widespread purifying selection; (2) elevated directional adaptive selection was observed in the T. truncatus mitogenome, perhaps linked to the cold climates adaptation of T. truncatus; and (3) the strength of selection varied across genes, and diversifying positive selection detected on ND4 and ATP6 pointed to their crucial roles during adaptation to different climatic conditions. This study gained insight into the mitogenome evolution in relation to the climate adaptation.

Selective Pressures on Human Cancer Genes along the Evolution of Mammals

Cancer is a disease driven by both somatic mutations that increase survival and proliferation of cell lineages and the evolution of genes associated with cancer risk in populations. Several genes associated with cancer in humans, hereafter cancer genes, show evidence of germline positive selection among species. Taking advantage of a large collection of mammalian genomes, we systematically looked for signatures of germline positive selection in 430 cancer genes available in COSMIC. We identified 40 cancer genes with a robust signal of positive selection in mammals. We found evidence for fewer selective constraints—higher number of non-synonymous substitutions per non-synonymous site to the number of synonymous substitutions per synonymous site (dN/dS)—and higher incidence of positive selection—more positively selected sites—in cancer genes bearing germline and recessive mutations that predispose to cancer. This finding suggests a potential association between relaxed selection, positive selection, and risk of hereditary cancer. On the other hand, we did not find significant differences in terms of tissue or gene type. Human cancer genes under germline positive selection in mammals are significantly enriched in the processes of DNA repair, with high presence of Fanconi anaemia/Breast Cancer A (FA/BRCA) pathway components and T cell proliferation genes. We also show that the inferred positively selected sites in the two genes with the strongest signal of positive selection, i.e., BRCA2 and PTPRC, are in regions of functional relevance, which could be relevant to cancer susceptibility.

Genome Sequence of Peacock Reveals the Peculiar Case of a Glittering Bird

The unique ornamental features and extreme sexual traits of Peacock have always intrigued scientists and naturalists for centuries. However, the genomic basis of these phenotypes are yet unknown. Here, we report the first genome sequence and comparative analysis of peacock with the high quality genomes of chicken, turkey, duck, flycatcher and zebra finch. Genes involved in early developmental pathways including TGF-β, BMP, and Wnt signaling, which have been shown to be involved in feather patterning, bone morphogenesis, and skeletal muscle development, revealed signs of adaptive evolution and provided useful clues on the phenotypes of peacock. Innate and adaptive immune genes involved in complement system and T-cell response also showed signs of adaptive evolution in peacock suggesting their possible role in building a robust immune system which is consistent with the predictions of the Hamilton–Zuk hypothesis. This study provides novel genomic and evolutionary insights into the molecular understanding toward the phenotypic evolution of Indian peacock.

Correlated Selection on Amino Acid Deletion and Replacement in Mammalian Protein Sequences

A low ratio of nonsynonymous and synonymous substitution rates (dN/dS) at a codon is an indicator of functional constraint caused by purifying selection. Intuitively, the functional constraint would also be expected to prevent such a codon from being deleted. However, to the best of our knowledge, the correlation between the rates of deletion and substitution has never actually been estimated. Here, we use 8595 protein-coding region sequences from nine mammalian species to examine the relationship between deletion rate and dN/dS. We find significant positive correlations at the levels of both sites and genes. We compared our data against controls consisting of simulated coding sequences evolving along identical phylogenetic trees, where deletions occur independently of substitutions. A much weaker correlation was found in the corresponding simulated sequences, probably caused by alignment errors. In the real data, the correlations cannot be explained by alignment errors. Separate investigations on nonsynonymous (dN) and synonymous (dS) substitution rates indicate that the correlation is most likely due to a similarity in patterns of selection rather than in mutation rates.

Position Matters: Network Centrality Considerably Impacts Rates of Protein Evolution in the Human Protein-Protein Interaction Network

The proteins of any organism evolve at disparate rates. A long list of factors affecting rates of protein evolution have been identified. However, the relative importance of each factor in determining rates of protein evolution remains unresolved. The prevailing view is that evolutionary rates are dominantly determined by gene expression, and that other factors such as network centrality have only a marginal effect, if any. However, this view is largely based on analyses in yeasts, and accurately measuring the importance of the determinants of rates of protein evolution is complicated by the fact that the different factors are often correlated with each other, and by the relatively poor quality of available functional genomics data sets. Here, we use correlation, partial correlation and principal component regression analyses to measure the contributions of several factors to the variability of the rates of evolution of human proteins. For this purpose, we analyzed the entire human protein–protein interaction data set and the human signal transduction network—a network data set of exceptionally high quality, obtained by manual curation, which is expected to be virtually free from false positives. In contrast with the prevailing view, we observe that network centrality (measured as the number of physical and nonphysical interactions, betweenness, and closeness) has a considerable impact on rates of protein evolution. Surprisingly, the impact of centrality on rates of protein evolution seems to be comparable, or even superior according to some analyses, to that of gene expression. Our observations seem to be independent of potentially confounding factors and from the limitations (biases and errors) of interactomic data sets.

Correlates of evolutionary rates in the murine sperm proteome

Background

Protein-coding genes expressed in sperm evolve at different rates. To gain deeper insight into the factors underlying this heterogeneity we examined the relative importance of a diverse set of previously described rate correlates in determining the evolution of murine sperm proteins.

Results

Using partial rank correlations we detected several major rate indicators: Phyletic gene age, numbers of protein-protein interactions, and survival essentiality emerged as particularly important rate correlates in murine sperm proteins. Tissue specificity, numbers of paralogs, and untranslated region lengths also correlate significantly with sperm genes’ evolutionary rates, albeit to a lesser extent. Multifunctionality, coding sequence or average intron lengths, and mean expression level have insignificant or virtually no independent effects on evolutionary rates in murine sperm genes. Gene ontology enrichment analyses of three equally sized murine sperm protein groups classified based on their evolutionary rates indicate strongest sperm-specific functional specialization in the most quickly evolving gene class.

Conclusions

We propose a model according to which slowly evolving murine sperm proteins tend to be constrained by factors such as survival essentiality, network connectivity, and/or broad expression. In contrast, evolutionary change may arise especially in less constrained sperm proteins, which might, moreover, be prone to specialize to reproduction-related functions. Our results should be taken into account in future studies on rate variations of reproductive genes.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1157-6) contains supplementary material, which is available to authorized users.

Disulfide Bonds Enable Accelerated Protein Evolution

The different proteins of any proteome evolve at enormously different rates. What factors contribute to this variability, and to what extent, is still a largely open question. We hypothesized that disulfide bonds, by increasing protein stability, should make proteins' structures relatively independent of their amino acid sequences, thus acting as buffers of deleterious mutations and enabling accelerated sequence evolution. In agreement with this hypothesis, we observed that membrane proteins with disulfide bonds evolved 88% faster than those without disulfide bonds, and that extracellular proteins with disulfide bonds evolved 49% faster than those without disulfide bonds. In addition, genes encoding proteins with disulfide bonds exhibit an increased likelihood of showing signatures of positive selection. Multivariate analyses indicate that the trend is independent of a number of potentially confounding factors. The effect, however, is not observed among the longest proteins, which can become stabilized by mechanisms other than disulfide bonds.