Less is more: an adaptive branch-site random effects model for efficient detection of episodic diversifying selection

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.

The genomic and cellular basis of biosynthetic innovation in rove beetles

How evolution at the cellular level potentiates change at the macroevolutionary level is a major question in evolutionary biology. With >66,000 described species, rove beetles (Staphylinidae) comprise the largest metazoan family. Their exceptional radiation has been coupled to pervasive biosynthetic innovation whereby numerous lineages bear defensive glands with diverse chemistries. Here, we combine comparative genomic and single-cell transcriptomic data from across the largest rove beetle clade, Aleocharinae. We retrace the functional evolution of two novel secretory cell types that together comprise the tergal gland-a putative catalyst behind Aleocharinae's megadiversity. We identify key genomic contingencies that were critical to the assembly of each cell type and their organ-level partnership in manufacturing the beetle's defensive secretion. This process hinged on evolving a mechanism for regulated production of noxious benzoquinones that appears convergent with plant toxin release systems, and synthesis of an effective benzoquinone solvent that weaponized the total secretion. We show that this cooperative biosynthetic system arose at the Jurassic-Cretaceous boundary, and that following its establishment, both cell types underwent ∼150 million years of stasis, their chemistry and core molecular architecture maintained almost clade-wide as Aleocharinae radiated globally into tens of thousands of lineages. Despite this deep conservation, we show that the two cell types have acted as substrates for the emergence of adaptive, biochemical novelties-most dramatically in symbiotic lineages that have infiltrated social insect colonies and produce host behavior-manipulating secretions. Our findings uncover genomic and cell type evolutionary processes underlying the origin, functional conservation and evolvability of a chemical innovation in beetles.

Mapping the Dynamics of Contemporary PRRSV-2 Evolution and Its Emergence and Spreading Hotspots in the U.S. Using Phylogeography

The repeated emergence of new genetic variants of PRRSV-2, the virus that causes porcine reproductive and respiratory syndrome (PRRS), reflects its rapid evolution and the failure of previous control efforts. Understanding spatiotemporal heterogeneity in variant emergence and spread is critical for future outbreak prevention. Here, we investigate how the pace of evolution varies across time and space, identify the origins of sub-lineage emergence, and map the patterns of the inter-regional spread of PRRSV-2 Lineage 1 (L1)-the current dominant lineage in the U.S. We performed comparative phylogeographic analyses on subsets of 19,395 viral ORF5 sequences collected across the U.S. and Canada between 1991 and 2021. The discrete trait analysis of multiple spatiotemporally stratified sampled sets (n = 500 each) was used to infer the ancestral geographic region and dispersion of each sub-lineage. The robustness of the results was compared to that of other modeling methods and subsampling strategies. Generally, the spatial spread and population dynamics varied across sub-lineages, time, and space. The Upper Midwest was a main spreading hotspot for multiple sub-lineages, e.g., L1C and L1F, though one of the most recent emergence events (L1A(2)) spread outwards from the east. An understanding of historical patterns of emergence and spread can be used to strategize disease control and the containment of emerging variants.

Utilising natural diversity of kinases to rationally engineer interactions with the angiosperm immune receptor ZAR1

The highly conserved angiosperm immune receptor HOPZ-ACTIVATED RESISTANCE1 (ZAR1) recognises the activity of diverse pathogen effector proteins by monitoring the ZED1-related kinase (ZRK) family. Understanding how ZAR1 achieves interaction specificity for ZRKs may allow for the expansion of the ZAR1-kinase recognition repertoire to achieve novel pathogen recognition outside of model species. We took advantage of the natural diversity of Arabidopsis thaliana kinases to probe the ZAR1-kinase interaction interface and found that A. thaliana ZAR1 (AtZAR1) can interact with most ZRKs, except ZRK7. We found evidence of alternative splicing of ZRK7, resulting in a protein that can interact with AtZAR1. Despite high sequence conservation of ZAR1, interspecific ZAR1-ZRK pairings resulted in the autoactivation of cell death. We showed that ZAR1 interacts with a greater diversity of kinases than previously thought, while still possessing the capacity for specificity in kinase interactions. Finally, using AtZAR1-ZRK interaction data, we rationally increased ZRK10 interaction strength with AtZAR1, demonstrating the feasibility of the rational design of a ZAR1-interacting kinase. Overall, our findings advance our understanding of the rules governing ZAR1 interaction specificity, with promising future directions for expanding ZAR1 immunodiversity.

Evolution of insect innate immunity through domestication of bacterial toxins

Toxin cargo genes are often horizontally transferred by phages between bacterial species and are known to play an important role in the evolution of bacterial pathogenesis. Here, we show how these same genes have been horizontally transferred from phage or bacteria to animals and have resulted in novel adaptations. We discovered that two widespread bacterial genes encoding toxins of animal cells, cytolethal distending toxin subunit B (cdtB) and apoptosis-inducing protein of 56 kDa (aip56), were captured by insect genomes through horizontal gene transfer from bacteria or phages. To study the function of these genes in insects, we focused on Drosophila ananassae as a model. In the D. ananassae subgroup species, cdtB and aip56 are present as singular (cdtB) or fused copies (cdtB::aip56) on the second chromosome. We found that cdtB and aip56 genes and encoded proteins were expressed by immune cells, some proteins were localized to the wasp embryo's serosa, and their expression increased following parasitoid wasp infection. Species of the ananassae subgroup are highly resistant to parasitoid wasps, and we observed that D. ananassae lines carrying null mutations in cdtB and aip56 toxin genes were more susceptible to parasitoids than the wild type. We conclude that toxin cargo genes were captured by these insects millions of years ago and integrated as novel modules into their innate immune system. These modules now represent components of a heretofore undescribed defense response and are important for resistance to parasitoid wasps. Phage or bacterially derived eukaryotic toxin genes serve as macromutations that can spur the instantaneous evolution of novelty in animals.

Positive and Relaxed Selective Pressures Have Both Strongly Influenced the Evolution of Cryonotothenioid Fishes during Their Radiation in the Freezing Southern Ocean

Evolution in the chronic cold of the Southern Ocean has had a profound influence on the physiology of cryonotothenioid fishes. However, the suite of genetic changes underlying the physiological gains and losses in these fishes is still poorly surveyed. By identifying the genomic signatures of selection, this study aims to identify the functional classes of genes that have been changed following two major physiological transitions: the onset of freezing temperatures and the loss of hemoproteins. Looking at the changes that followed the onset of freezing temperatures, positive selective pressure was found among a set of broadly acting gene regulatory factors, suggesting a route through which cryonotothenioid gene expression has been retooled for life in the cold. Further, genes related to the cell cycle and cellular adhesion were found under positive selection suggesting that both present key challenges to life in freezing waters. By contrast, genes showing signatures of the relaxation of selective pressure showed a narrower biological impact, acting on genes related to mitochondrial function. Finally, although chronic cold-water temperatures appear correlated with substantial genetic change, the loss of hemoproteins resulted in little observable change in protein-coding genes relative to their red-blooded relatives. Combined, the influence of positive and relaxed selection shows that long-term exposure to cold has led to profound changes in cryonotothenioid genomes that may make it challenging for them to adapt to a rapidly changing climate.

Signatures of Adaptation, Constraints, and Potential Redundancy in the Canonical Immune Genes of a Key Pollinator

All organisms require an immune system to recognize, differentiate, and defend against pathogens. From an evolutionary perspective, immune systems evolve under strong selective pressures exerted by fast-evolving pathogens. However, the functional diversity of the immune system means that different immune components and their associated genes may evolve under varying forms of selection. Insect pollinators, which provide essential ecosystem services, are an important system in which to understand how selection has shaped immune gene evolution as their populations are experiencing declines with pathogens highlighted as a potential contributing factor. To improve our understanding of the genetic variation found in the immune genes of an essential pollinator, we performed whole-genome resequencing of wild-caught Bombus terrestris males. We first assessed nucleotide diversity and extended haplotype homozygosity for canonical immune genes finding the strongest signatures of positive selection acting on genes involved in pathogen recognition and antiviral defense, possibly driven by growing pathogen spread in wild populations. We also identified immune genes evolving under strong purifying selection, highlighting potential constraints on the bumblebee immune system. Lastly, we highlight the potential loss of function alleles present in the immune genes of wild-caught haploid males, suggesting that such genes are potentially less essential for development and survival and represent redundancy in the gene repertoire of the bumblebee immune system. Collectively, our analysis provides novel insights into the recent evolutionary history of the immune system of a key pollinator, highlighting targets of selection, constraints to adaptation, and potential redundancy.

Aminoacyl-tRNA Synthetase Evolution within the Dynamic Tripartite Translation System of Plant Cells

Eukaryotes maintain separate protein translation systems for nuclear and organellar genes, including distinct sets of tRNAs and aminoacyl-tRNA synthetases (aaRSs). In animals, mitochondrial-targeted aaRSs are expressed at lower levels and are less conserved in sequence than cytosolic aaRSs involved in translation of nuclear mRNAs, likely reflecting lower translational demands in mitochondria. In plants, translation is further complicated by the presence of plastids, which share most aaRSs with mitochondria. In addition, plant mitochondrial tRNA pools have a dynamic history of gene loss and functional replacement by tRNAs from other compartments. To investigate the consequences of these distinctive features of translation in plants, we analyzed sequence evolution in angiosperm aaRSs. In contrast to previously studied eukaryotic systems, we found that plant organellar and cytosolic aaRSs exhibit only a small difference in expression levels, and organellar aaRSs are slightly more conserved than cytosolic aaRSs. We hypothesize that these patterns result from high translational demands associated with photosynthesis in mature chloroplasts. We also investigated aaRS evolution in Sileneae, an angiosperm lineage with extensive mitochondrial tRNA replacement and aaRS retargeting. We predicted positive selection for changes in aaRS sequence resulting from these recent changes in subcellular localization and tRNA substrates but found little evidence for accelerated sequence divergence. Overall, the complex tripartite translation system in plant cells appears to have imposed more constraints on the long-term evolutionary rates of organellar aaRSs compared with other eukaryotic lineages, and plant aaRS protein sequences appear largely robust to more recent perturbations in subcellular localization and tRNA interactions.

Beginner's Guide on the Use of PAML to Detect Positive Selection

The CODEML program in the PAML package has been widely used to analyze protein-coding gene sequences to estimate the synonymous and nonsynonymous rates (dS and dN) and to detect positive Darwinian selection driving protein evolution. For users not familiar with molecular evolutionary analysis, the program is known to have a steep learning curve. Here, we provide a step-by-step protocol to illustrate the commonly used tests available in the program, including the branch models, the site models, and the branch-site models, which can be used to detect positive selection driving adaptive protein evolution affecting particular lineages of the species phylogeny, affecting a subset of amino acid residues in the protein, and affecting a subset of sites along prespecified lineages, respectively. A data set of the myxovirus (Mx) genes from ten mammal and two bird species is used as an example. We discuss a new feature in CODEML that allows users to perform positive selection tests for multiple genes for the same set of taxa, as is common in modern genome-sequencing projects. The PAML package is distributed at https://github.com/abacus-gene/paml under the GNU license, with support provided at its discussion site (https://groups.google.com/g/pamlsoftware). Data files used in this protocol are available at https://github.com/abacus-gene/paml-tutorial.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

The phylogenomics and evolutionary dynamics of the organellar genomes in carnivorous Utricularia and Genlisea species (Lentibulariaceae)

Utricularia and Genlisea are highly specialized carnivorous plants whose phylogenetic history has been poorly explored using phylogenomic methods. Additional sampling and genomic data are needed to advance our phylogenetic and taxonomic knowledge of this group of plants. Within a comparative framework, we present a characterization of plastome (PT) and mitochondrial (MT) genes of 26 Utricularia and six Genlisea species, with representatives of all subgenera and growth habits. All PT genomes maintain similar gene content, showing minor variation across the genes located between the PT junctions. One exception is a major variation related to different patterns in the presence and absence of ndh genes in the small single copy region, which appears to follow the phylogenetic history of the species rather than their lifestyle. All MT genomes exhibit similar gene content, with most differences related to a lineage-specific pseudogenes. We find evidence for episodic positive diversifying selection in PT and for most of the Utricularia MT genes that may be related to the current hypothesis that bladderworts' nuclear DNA is under constant ROS oxidative DNA damage and unusual DNA repair mechanisms, or even low fidelity polymerase that bypass lesions which could also be affecting the organellar genomes. Finally, both PT and MT phylogenetic trees were well resolved and highly supported, providing a congruent phylogenomic hypothesis for Utricularia and Genlisea clade given the study sampling.

Barbel regeneration and function divergence in red-tail catfish (Hemibagrus wyckioides) based on the chromosome-level genomes and comparative transcriptomes

Catfish (Siluriformes) are one of the most diverse vertebrate orders and are characterized by whisker-like barbels, which are important sensory organs in most of teleosts. However, their specific biological functions are still unclear. Red-tail catfish (Hemibagrus wyckioides) is well-known catfish species with four pairs of barbels, of which the maxillary barbels reach two-thirds of the body length. In this study, a 776.58 Mb high-quality chromosome-level genome was assembled into 29 chromosomes. Comparative genome data indicated that the barbeled regeneration gene ccl33 has expanded into 11 tandemly duplicated copies. Transcriptome data revealed the functional differentiation of different barbels and suggested that the maxillary barbel might be necessary for water temperature perception. Taste receptor genes were also characterized in teleosts with different food habits. Selection pressures were revealed to affect the sugar-based solute transport domain of the sweet taste receptor gene t1r2 in carnivorous fishes. In addition, the bitter taste receptor gene t2r200 was found to be lost from the genomes of four catfish species. Therefore, our study provides a genomic foundation for understanding the regeneration and functional differentiation of barbels in red-tail catfish and also reveals novel insights into the feeding evolution of fish species with different feeding habits.

Backbone phylogeny and adaptive evolution of Pleurospermum s. l.: New insights from phylogenomic analyses of complete plastome data

Pleurospermum is a taxonomically challenging taxon of Apiaceae, as its circumscription and composition remain controversial for morphological similarities with several related genera, leading to a dispute between Pleurospermum in the broad sense and strict sense. While evidence from previous molecular studies recognized plural branching lineages within the Pleurospermum s. l., it did not support the latest delimitation of Pleurospermum s. str. by only two closely related northern species. So far, no proper delimitation for Pleurospermum has come up, and many of the plural taxa in Pleurospermum s. l. remain unresolved, which may be due to poor phylogenetic resolution yielded barely from ITS sequences. Herein, we newly assembled 40 complete plastomes from 36 species of Pleurospermum s. l. and related genera, 34 of which were first reported and generated a well-resolved backbone phylogeny in a framework of the subfamily Apioideae. From the phylogeny with greatly improved resolution, a total of six well-supported monophyletic lineages within Pleurospermum s. l. were recognized falling in different major clades of Apioideae. Combining morphological characteristics with phylogenetic inference, we suggested to re-delimit the Pleurospermum s. str. by introducing nine species mainly from the Himalayan regions and proposed its boundary features; the remaining species were suggested to be excluded from Pleurospermum to incorporate into their more related taxa being revealed. On this basis, the plastome comparison revealed not only the high conservatism but also the mild differences among lineages in plastome structure and gene evolution. Overall, our study provided a backbone phylogeny essential for further studies of the taxonomically difficult taxa within Pleurospermum s. l.

Comparative genomics reveals the diversification of triterpenoid biosynthesis and origin of ocotillol-type triterpenes in Panax

Gene duplication is assumed to be the major force driving the evolution of metabolite biosynthesis in plants. Freed from functional burdens, duplicated genes can mutate toward novelties until fixed due to selective fitness. However, the extent to which this mechanism has driven the diversification of metabolite biosynthesis remains to be tested. Here we performed comparative genomics analysis and functional characterization to evaluate the impact of gene duplication on the evolution of triterpenoid biosynthesis using Panax species as models. We found that whole-genome duplications (WGDs) occurred independently in Araliaceae and Apiaceae lineages. Comparative genomics revealed the evolutionary trajectories of triterpenoid biosynthesis in plants, which was mainly promoted by WGDs and tandem duplication. Lanosterol synthase (LAS) was likely derived from a tandem duplicate of cycloartenol synthase that predated the emergence of Nymphaeales. Under episodic diversifying selection, the LAS gene duplicates produced by γ whole-genome triplication have given rise to triterpene biosynthesis in core eudicots through neofunctionalization. Moreover, functional characterization revealed that oxidosqualene cyclases (OSCs) responsible for synthesizing dammarane-type triterpenes in Panax species were also capable of producing ocotillol-type triterpenes. Genomic and biochemical evidence suggested that Panax genes encoding the above OSCs originated from the specialization of one OSC gene duplicate produced from a recent WGD shared by Araliaceae (Pg-β). Our results reveal the crucial role of gene duplication in diversification of triterpenoid biosynthesis in plants and provide insight into the origin of ocotillol-type triterpenes in Panax species.

Molecular evolution of the ependymin-related gene epdl2 in African weakly electric fish

Gene duplication and subsequent molecular evolution can give rise to taxon-specific gene specializations. In previous work, we found evidence that African weakly electric fish (Mormyridae) may have as many as three copies of the epdl2 gene, and the expression of two epdl2 genes is correlated with electric signal divergence. Epdl2 belongs to the ependymin-related family (EPDR), a functionally diverse family of secretory glycoproteins. In this study, we first describe vertebrate EPDR evolution and then present a detailed evolutionary history of epdl2 in Mormyridae with emphasis on the speciose genus Paramormyrops. Using Sanger sequencing, we confirm three apparently functional epdl2 genes in Paramormyrops kingsleyae. Next, we developed a nanopore-based amplicon sequencing strategy and bioinformatics pipeline to obtain and classify full-length epdl2 gene sequences (N = 34) across Mormyridae. Our phylogenetic analysis proposes three or four epdl2 paralogs dating from early Paramormyrops evolution. Finally, we conducted selection tests which detected positive selection around the duplication events and identified ten sites likely targeted by selection in the resulting paralogs. These sites' locations in our modeled 3D protein structure involve four sites in ligand binding and six sites in homodimer formation. Together, these findings strongly imply an evolutionary mechanism whereby epdl2 genes underwent selection-driven functional specialization after tandem duplications in the rapidly speciating Paramormyrops. Considering previous evidence, we propose that epdl2 may contribute to electric signal diversification in mormyrids, an important aspect of species recognition during mating.

Evolution of the Growth Hormone Gene Duplication in Passerine Birds

Birds of the order Passeriformes represent the most speciose order of land vertebrates. Despite strong scientific interest in this super-radiation, genetic traits unique to passerines are not well characterized. A duplicate copy of growth hormone (GH) is the only gene known to be present in all major lineages of passerines, but not in other birds. GH genes plausibly influence extreme life history traits that passerines exhibit, including the shortest embryo-to-fledging developmental period of any avian order. To unravel the implications of this GH duplication, we investigated the molecular evolution of the ancestral avian GH gene (GH or GH1) and the novel passerine GH paralog (GH2), using 497 gene sequences extracted from 342 genomes. Passerine GH1 and GH2 are reciprocally monophyletic, consistent with a single duplication event from a microchromosome onto a macrochromosome in a common ancestor of extant passerines. Additional chromosomal rearrangements have changed the syntenic and potential regulatory context of these genes. Both passerine GH1 and GH2 display substantially higher rates of nonsynonymous codon change than non-passerine avian GH, suggesting positive selection following duplication. A site involved in signal peptide cleavage is under selection in both paralogs. Other sites under positive selection differ between the two paralogs, but many are clustered in one region of a 3D model of the protein. Both paralogs retain key functional features and are actively but differentially expressed in two major passerine suborders. These phenomena suggest that GH genes may be evolving novel adaptive roles in passerine birds.

A de novo assembled genome of the Tibetan Partridge (Perdix hodgsoniae) and its high-altitude adaptation

The Tibetan Partridge (Perdix hodgsoniae) is an endemic species distributed in high-altitude areas of 3600-5600 m on the Qinghai-Tibet Plateau. To explore how the species is adapted to the high elevation environment, we assembled a draft genome based on both the Illumina and PacBio sequencing platforms with its population genetics and genomics analysis. In total, 134.74 Gb short reads and 30.81 Gb long reads raw data were generated. The 1.05-Gb assembled genome had a contig N50 of 4.56 Mb, with 91.94% complete BUSCOs. The 17 457 genes were annotated, and 11.35% of the genome was composed of repeat sequences. The phylogenetic tree showed that P. hodgsoniae was located at the basal position of the clade, including Golden Pheasant (Chrysolophus pictus), Common Pheasant (Phasianus colchicus), and Mikado Pheasant (Syrmaticus mikado). We found that 1014, 2595, and 2732 of the 6641 one-to-one orthologous genes were under positive selection in P. hodgsoniae, detected using PAML, BUSTED, and aBSREL programs, respectively, of which 965 genes were common under positive selection with 3 different programs. Several positively selected genes and immunity pathways relevant to high-altitude adaptation were detected. Gene family evolution showed that 99 gene families experienced significant expansion events, while 6 gene families were under contraction. The total number of olfactory receptor genes was relatively low in P. hodgsoniae. Genomic data provide an important resource for a further study on the evolutionary history of P. hodgsoniae, which provides a new insight into its high-altitude adaptation mechanisms.

Evaluation of the Function of ASFV Gene E66L in the Process of Virus Replication and Virulence in Swine

African swine fever virus (ASFV) is the etiological agent of an economically important disease of swine currently affecting large areas of Africa, Eurasia and the Caribbean. ASFV has a complex structure harboring a large dsDNA genome which encodes for more than 160 proteins. One of the proteins, E66L, has recently been involved in arresting gene transcription in the infected host cell. Here, we investigate the role of E66L in the processes of virus replication in swine macrophages and disease production in domestic swine. A recombinant ASFV was developed (ASFV-G-∆E66L), from the virulent parental Georgia 2010 isolate (ASFV-G), harboring the deletion of the E66L gene as a tool to assess the role of the gene. ASFV-G-∆E66L showed that the E66L gene is non-essential for ASFV replication in primary swine macrophages when compared with the parental highly virulent field isolate ASFV-G. Additionally, domestic pigs infected with ASFV-G-∆E66L developed a clinical disease undistinguishable from that produced by ASFV-G. Therefore, E66L is not involved in virus replication or virulence in domestic pigs.

That H9N2 avian influenza viruses circulating in different regions gather in the same live-poultry market poses a potential threat to public health

H9N2 avian influenza viruses are endemic and persistent in China, but those that are prevalent in different provinces are also causes of wide epidemics, related to the spread of wild birds and the cross-regional trade in live poultry. For the past 4 years, beginning in 2018, we have sampled a live-poultry market in Foshan, Guangdong, in this ongoing study. In addition to the prevalence of H9N2 avian influenza viruses in China during this period, we identified isolates from the same market belonging to clade A and clade B, which diverged in 2012-2013, and clade C, which diverged in 2014-2016, respectively. An analysis of population dynamics revealed that, after a critical divergence period from 2014 to 2016, the genetic diversity of H9N2 viruses peaked in 2017. Our spatiotemporal dynamics analysis found that clade A, B, and C, which maintain high rates of evolution, have different prevalence ranges and transmission paths. Clades A and B were mainly prevalent in East China in the early stage, and then spread to Southern China, becoming epidemic with clade C. Strains from different regions converge at the same live-poultry market to communicate, which may be one reasons the H9N2 viruses are difficult to eradicate and increasingly dominant throughout China. Selection pressure and molecular analysis have demonstrated that single amino acid polymorphisms at key receptor binding sites 156, 160, and 190 under positive selection pressure, suggesting that H9N2 viruses are undergoing mutations to adapt to new hosts. Live-poultry markets are important because people who visit them have frequent contact with poultry, H9N2 viruses from different regions converge at these markets and spread through contact between live birds and humans, generating increased risks of human exposure to these viruses and threatening public health safety. Thus, it is important to reducing the cross-regional trade of live poultry and strengthening the monitoring of avian influenza viruses in live-poultry markets to reduce the spread of avian influenza viruses.

Accelerated evolution of SARS-CoV-2 in free-ranging white-tailed deer

While SARS-CoV-2 has sporadically infected a wide range of animal species worldwide1, the virus has been repeatedly and frequently detected in white-tailed deer in North America2â€"7. The zoonotic origins of this pandemic virus highlight the need to fill the vast gaps in our knowledge of SARS-CoV-2 ecology and evolution in non-human hosts. Here, we detected SARS-CoV-2 was introduced from humans into white-tailed deer more than 30 times in Ohio, USA during November 2021-March 2022. Subsequently, deer-to-deer transmission persisted for 2-8 months, which disseminated across hundreds of kilometers. We discovered that alpha and delta variants evolved in white-tailed deer at three-times the rate observed in humans. Newly developed Bayesian phylogenetic methods quantified how SARS-CoV-2 evolution is not only faster in white-tailed deer but driven by different mutational biases and selection pressures. White-tailed deer are not just short-term recipients of human viral diversity but serve as reservoirs for alpha and other variants to evolve in new directions after going extinct in humans. The long-term effect of this accelerated evolutionary rate remains to be seen as no critical phenotypic changes were observed in our animal model experiments using viruses isolated from white-tailed deer. Still, SARS-CoV-2 viruses have transmitted in white-tailed deer populations for a relatively short duration, and the risk of future changes may have serious consequences for humans and livestock.

Conserved chloroplast genome sequences of the genus Clerodendrum Linn. (Lamiaceae) as a super-barcode

BACKGROUND

The plants of the genus Clerodendrum L. have great potential for development as an ornamental and important herbal resource. There is no significant morphological difference among many species of the genus Clerodendrum, which will lead to confusion among the herbs of this genus and ultimately affect the quality of the herbs. The chloroplast genome will contribute to the development of new markers used for the identification and classification of species.

METHODS AND RESULTS

Here, we obtained the complete chloroplast genome sequences of Clerodendrum chinense (Osbeck) Mabberley and Clerodendrum thomsoniae Balf.f. using the next generation DNA sequencing technology. The chloroplast genomes of the two species all encode a total of 112 unique genes, including 80 protein-coding, 28 tRNA, and four rRNA genes. A total of 44-42 simple sequence repeats, 19-16 tandem repeats and 44-44 scattered repetitive sequences were identified. Phylogenetic analyses showed that the nine Clerodendrum species were classified into two clades and together formed a monophyletic group. Selective pressure analyses of 77 protein-coding genes showed that there was no gene under positive selection in the Clerodendrum branch. Analyses of sequence divergence found two intergenic regions: trnH-GUG-psbA, nhdD-psaC, exhibiting a high degree of variations. Meanwhile, there was no hypervariable region identified in protein coding genes. However, the sequence identities of these two intergenic spacers (IGSs) are greater than 99% among some species, which will result in the two IGSs not being used to distinguish Clerodendrum species. Analysis of the structure at the LSC (Large single copy) /IR (Inverted repeat) and SSC (Small single copy)/IR boundary regions showed dynamic changes. The above results showed that the complete chloroplast genomes can be used as a super-barcode to identify these Clerodendrum species. The study lay the foundation for the understanding of the evolutionary process of the genus Clerodendrum.

Early Genomic Surveillance and Phylogeographic Analysis of Getah Virus, a Reemerging Arbovirus, in Livestock in China

Getah virus (GETV) mainly causes disease in livestock and may pose an epidemic risk due to its expanding host range and the potential of long-distance dispersal through animal trade. Here, we used metagenomic next-generation sequencing (mNGS) to identify GETV as the pathogen responsible for reemerging swine disease in China and subsequently estimated key epidemiological parameters using phylodynamic and spatially-explicit phylogeographic approaches. The GETV isolates were able to replicate in a variety of cell lines, including human cells, and showed high pathogenicity in a mouse model, suggesting the potential for more mammal hosts. We obtained 16 complete genomes and 79 E2 gene sequences from viral strains collected in China from 2016 to 2021 through large-scale surveillance among livestock, pets, and mosquitoes. Our phylogenetic analysis revealed that three major GETV lineages are responsible for the current epidemic in livestock in China. We identified three potential positively selected sites and mutations of interest in E2, which may impact the transmissibility and pathogenicity of the virus. Phylodynamic inference of the GETV demographic dynamics identified an association between livestock meat consumption and the evolution of viral genetic diversity. Finally, phylogeographic reconstruction of GETV dispersal indicated that the sampled lineages have preferentially circulated within areas associated with relatively higher mean annual temperature and pig population density. Our results highlight the importance of continuous surveillance of GETV among livestock in southern Chinese regions associated with relatively high temperatures. IMPORTANCE Although livestock is known to be the primary reservoir of Getah virus (GETV) in Asian countries, where identification is largely based on serology, the evolutionary history and spatial epidemiology of GETV in these regions remain largely unknown. Through our sequencing efforts, we provided robust support for lineage delineation of GETV and identified three major lineages that are responsible for the current epidemic in livestock in China. We further analyzed genomic and epidemiological data to reconstruct the recent demographic and dispersal history of GETV in domestic animals in China and to explore the impact of environmental factors on its genetic diversity and its diffusion. Notably, except for livestock meat consumption, other pig-related factors such as the evolution of live pig transport and pork production do not show a significant association with the evolution of viral genetic diversity, pointing out that further studies should investigate the potential contribution of other host species to the GETV outbreak. Our analysis of GETV demonstrates the need for wider animal species surveillance and provides a baseline for future studies of the molecular epidemiology and early warning of emerging arboviruses in China.

Killer Knots: Molecular Evolution of Inhibitor Cystine Knot Toxins in Wandering Spiders (Araneae: Ctenidae)

Venom expressed by the nearly 50,000 species of spiders on Earth largely remains an untapped reservoir of a diverse array of biomolecules with potential for pharmacological and agricultural applications. A large fraction of the noxious components of spider venoms are a functionally diverse family of structurally related polypeptides with an inhibitor cystine knot (ICK) motif. The cysteine-rich nature of these toxins makes structural elucidation difficult, and most studies have focused on venom components from the small handful of medically relevant spider species such as the highly aggressive Brazilian wandering spider Phoneutria nigriventer. To alleviate difficulties associated with the study of ICK toxins in spiders, we devised a comprehensive approach to explore the evolutionary patterns that have shaped ICK functional diversification using venom gland transcriptomes and proteomes from phylogenetically distinct lineages of wandering spiders and their close relatives. We identified 626 unique ICK toxins belonging to seven topological elaborations. Phylogenetic tests of episodic diversification revealed distinct regions between cysteine residues that demonstrated differential evidence of positive or negative selection, which may have structural implications towards the specificity and efficacy of these toxins. Increased taxon sampling and whole genome sequencing will provide invaluable insights to further understand the evolutionary processes that have given rise to this diverse class of toxins.

Comparative genomics of Hox and ParaHox genes among major lineages of Branchiopoda with emphasis on tadpole shrimps

Hox and ParaHox genes (HPHGs) are key developmental genes that pattern regional identity along the anterior–posterior body axis of most animals. Here, we identified HPHGs in tadpole shrimps (Pancrustacea, Branchiopoda, Notostraca), an iconic example of the so-called “living fossils” and performed a comparative genomics analysis of HPHGs and the Hox cluster among major branchiopod lineages. Notostraca possess the entire Hox complement, and the Hox cluster seems to be split into two different subclusters, although we were not able to support this finding with chromosome-level assemblies. However, the genomic structure of Hox genes in Notostraca appears more derived than that of Daphnia spp., which instead retains the plesiomorphic condition of a single compact cluster. Spinicaudata and Artemia franciscana show instead a Hox cluster subdivided across two or more genomic scaffolds with some orthologs either duplicated or missing. Yet, branchiopod HPHGs are similar among the various clades in terms of both intron length and number, as well as in their pattern of molecular evolution. Sequence substitution rates are in fact generally similar for most of the branchiopod Hox genes and the few differences we found cannot be traced back to natural selection, as they are not associated with any signals of diversifying selection or substantial switches in selective modes. Altogether, these findings do not support a significant stasis in the Notostraca Hox cluster and further confirm how morphological evolution is not tightly associated with genome dynamics.

The molecular evolution of genes previously associated with large sizes reveals possible pathways to cetacean gigantism

Cetaceans are a group of aquatic mammals with the largest body sizes among living animals, including giant representatives such as blue and fin whales. To understand the genetic bases of gigantism in cetaceans, we performed molecular evolutionary analyses on five genes (GHSR, IGF2, IGFBP2, IGFBP7, and EGF) from the growth hormone/insulin-like growth factor axis, and four genes (ZFAT, EGF, LCORL, and PLAG1) previously described as related to the size of species evolutionarily close to cetaceans, such as pigs, cows, and sheep. Our dataset comprised 19 species of cetaceans, seven of which are classified as giants because they exceed 10 m in length. Our results revealed signs of positive selection in genes from the growth hormone/insulin-like growth factor axis and also in those related to body increase in cetacean-related species. In addition, pseudogenization of the EGF gene was detected in the lineage of toothless cetaceans, Mysticeti. Our results suggest the action of positive selection on gigantism in genes that act both in body augmentation and in mitigating its consequences, such as cancer suppression when involved in processes such as division, migration, and cell development control.

Stabilization of membrane topologies by proteinaceous remorin scaffolds

In plants, the topological organization of membranes has mainly been attributed to the cell wall and the cytoskeleton. Additionally, few proteins, such as plant-specific remorins have been shown to function as protein and lipid organizers. Root nodule symbiosis requires continuous membrane re-arrangements, with bacteria being finally released from infection threads into membrane-confined symbiosomes. We found that mutations in the symbiosis-specific SYMREM1 gene result in highly disorganized perimicrobial membranes. AlphaFold modelling and biochemical analyses reveal that SYMREM1 oligomerizes into antiparallel dimers and may form a higher-order membrane scaffolding structure. This was experimentally confirmed when expressing this and other remorins in wall-less protoplasts is sufficient where they significantly alter and stabilize de novo membrane topologies ranging from membrane blebs to long membrane tubes with a central actin filament. Reciprocally, mechanically induced membrane indentations were equally stabilized by SYMREM1. Taken together we describe a plant-specific mechanism that allows the stabilization of large-scale membrane conformations independent of the cell wall.

Convergent genomics of longevity in rockfishes highlights the genetics of human life span variation

Longevity is a defining, heritable trait that varies dramatically between species. To resolve the genetic regulation of this trait, we have mined genomic variation in rockfishes, which range in longevity from 11 to over 205 years. Multiple shifts in rockfish longevity have occurred independently and in a short evolutionary time frame, thus empowering convergence analyses. Our analyses reveal a common network of genes under convergent evolution, encompassing established aging regulators such as insulin signaling, yet also identify flavonoid (aryl-hydrocarbon) metabolism as a pathway modulating longevity. The selective pressures on these pathways indicate the ancestral state of rockfishes was long lived and that the changes in short-lived lineages are adaptive. These pathways were also used to explore genome-wide association studies of human longevity, identifying the aryl-hydrocarbon metabolism pathway to be significantly associated with human survival to the 99th percentile. This evolutionary intersection defines and cross-validates a previously unappreciated genetic architecture that associates with the evolution of longevity across vertebrates.

Novel genome sequences and evolutionary dynamics of the North American anopheline species Anopheles freeborni, Anopheles crucians, Anopheles quadrimaculatus, and Anopheles albimanus

Anopheles mosquitoes are the principal vectors for malaria and lymphatic filariasis, and evidence for arboviral transmission under laboratory and natural contexts has been demonstrated. Vector management approaches require an understanding of the ecological, epidemiological, and biological contexts of the species in question, and increased interest in gene drive systems for vector control applications has resulted in an increased need for genome assemblies from understudied mosquito vector species. In this study, we present novel genome assemblies for Anopheles crucians, Anopheles freeborni, Anopheles albimanus, and Anopheles quadrimaculatus and examine the evolutionary relationship between these species. We identified 790 shared single-copy orthologs between the newly sequenced genomes and created a phylogeny using 673 of the orthologs, identifying 289 orthologs with evidence for positive selection on at least 1 branch of the phylogeny. Gene ontology terms such as calcium ion signaling, histone binding, and protein acetylation identified as being biased in the set of selected genes. These novel genome sequences will be useful in developing our understanding of the diverse biological traits that drive vectorial capacity in anophelines.

Wide distribution of D-xylose dehydrogenase in yeasts reveals a new element in the D-xylose metabolism for bioethanol production

D-xylose utilization by yeasts is an essential feature for improving second-generation ethanol production. However, industrial yeast strains are incapable of consuming D-xylose. Previous analyzes of D-xylose-consuming or fermenting yeast species reveal that the genomic features associated with this phenotype are complex and still not fully understood. Here we present a previously neglected yeast enzyme related to D-xylose metabolism, D-xylose dehydrogenase (XylDH), which is found in at least 105 yeast genomes. By analyzing the XylDH gene family, we brought evidence of gene evolution marked by purifying selection on codons and positive selection evidence in D-xylose-consuming and fermenting species, suggesting the importance of XylDH for D-xylose-related phenotypes in yeasts. Furthermore, although we found no putative metabolic pathway for XylDH in yeast genomes, namely the absence of three bacterial known pathways for this enzyme, we also provide its expression profile on D-xylose media following D-xylose reductase for two yeasts with publicly available transcriptomes. Based on these results, we suggest that XylDH plays an important role in D-xylose usage by yeasts, likely being involved in a cofactor regeneration system by reducing cofactor imbalance in the D-xylose reductase pathway.

Evolution of the orthopoxvirus core genome

Orthopoxviruses comprise several relevant pathogens, including the causative agent of smallpox and monkeypox virus. Analysis of orthopoxvirus genome evolution mainly focused on gene gains/losses. We instead analyzed core genes, which are conserved in all orthopoxviruses. We show that, despite their strong constraint, some genes involved in viral morphogenesis and transcription/replication were targets of pervasive positive selection, which was relatively uncommon in immunomodulatory genes. However at least three of the positively selected genes, E3L, A24R, and H3L, might have evolved in response to immune selection. Episodic positive selection was particularly common on the internal branches of the orthopox phylogeny and on the monkeypox virus lineage. The latter showed evidence of episodic positive selection at the D14L gene, which encodes a modulator of complement activation (MOPICE). Notably, two genes (B1R and A33R) targeted by episodic selection on more than one branch are involved in forms of intra-genomic conflict. Finally, we found that, in orthopoxvirus proteomes, intrinsically disordered regions (IDRs) tend to be less constrained and are common targets of positive selection. Extension of our analysis to all poxviruses showed no evidence that the IDR fraction differs with host range. Conversely, we found a strong effect of base composition, which was however not sufficient to explain IDR fraction. We thus suggest that, in poxviruses, the IDR fraction is maintained by modulating GC content to accommodate disorder-promoting codons. Overall, our data provide novel insight in orthopoxvirus evolution and provide a list of genes and sites that are expected to modulate viral phenotypes.

Prediction of effector protein structures from fungal phytopathogens enables evolutionary analyses

Elucidating the similarity and diversity of pathogen effectors is critical to understand their evolution across fungal phytopathogens. However, rapid divergence that diminishes sequence similarities between putatively homologous effectors has largely concealed the roots of effector evolution. Here we modelled the structures of 26,653 secreted proteins from 14 agriculturally important fungal phytopathogens, six non-pathogenic fungi and one oomycete with AlphaFold 2. With 18,000 successfully predicted folds, we performed structure-guided comparative analyses on two aspects of effector evolution: uniquely expanded sequence-unrelated structurally similar (SUSS) effector families and common folds present across the fungal species. Extreme expansion of lineage-specific SUSS effector families was found only in several obligate biotrophs, Blumeria graminis and Puccinia graminis. The highly expanded effector families were the source of conserved sequence motifs, such as the Y/F/WxC motif. We identified new classes of SUSS effector families that include known virulence factors, such as AvrSr35, AvrSr50 and Tin2. Structural comparisons revealed that the expanded structural folds further diversify through domain duplications and fusion with disordered stretches. Putatively sub- and neo-functionalized SUSS effectors could reconverge on regulation, expanding the functional pools of effectors in the pathogen infection cycle. We also found evidence that many effector families could have originated from ancestral folds conserved across fungi. Collectively, our study highlights diverse effector evolution mechanisms and supports divergent evolution as a major force in driving SUSS effector evolution from ancestral proteins.

Historic samples reveal loss of wild genotype through domestic chicken introgression during the Anthropocene

Human activities have precipitated a rise in the levels of introgressive gene flow among animals. The investigation of conspecific populations at different time points may shed light on the magnitude of human-mediated introgression. We used the red junglefowl Gallus gallus, the wild ancestral form of the chicken, as our study system. As wild junglefowl and domestic chickens readily admix, conservationists fear that domestic introgression into junglefowl may compromise their wild genotype. By contrasting the whole genomes of 51 chickens with 63 junglefowl from across their natural range, we found evidence of a loss of the wild genotype across the Anthropocene. When comparing against the genomes of junglefowl from approximately a century ago using rigorous ancient-DNA protocols, we discovered that levels of domestic introgression are not equal among and within modern wild populations, with the percentage of domestic ancestry around 20-50%. We identified a number of domestication markers in which chickens are deeply differentiated from historic junglefowl regardless of breed and/or geographic provenance, with eight genes under selection. The latter are involved in pathways dealing with development, reproduction and vision. The wild genotype is an allelic reservoir that holds most of the genetic diversity of G. gallus, a species which is immensely important to human society. Our study provides fundamental genomic infrastructure to assist in efforts to prevent a further loss of the wild genotype through introgression of domestic alleles.

Positive selection on ADAM10 builds species recognition in the synchronous spawning coral Acropora

The reef-building coral Acropora is a broadcast spawning hermaphrodite including more than 110 species in the Indo-Pacific. In addition, many sympatric species show synchronous spawning. The released gametes need to mate with conspecifics in the mixture of the gametes of many species for their species boundaries. However, the mechanism underlying the species recognition of conspecifics at fertilization remains unknown. We hypothesized that rapid molecular evolution (positive selection) in genes encoding gamete-composing proteins generates polymorphic regions that recognize conspecifics in the mixture of gametes from many species. We identified gamete proteins of Acropora digitifera using mass spectrometry and screened the genes that support branch site models that set the "foreground" branches showing strict fertilization specificity. ADAM10, ADAM17, Integrin α9, and Tetraspanin4 supported branch-site model and had positively selected site(s) that produced polymorphic regions. Therefore, we prepared antibodies against the proteins of A. digitifera that contained positively selected site(s) to analyze their functions in fertilization. The ADAM10 antibody reacted only with egg proteins of A. digitifera, and immunohistochemistry showed ADAM10 localized around the egg surface. Moreover, the ADAM10 antibody inhibited only A. digitifera fertilization but not the relative synchronous spawning species A. papillare. This study indicates that ADAM10 has evolved to gain fertilization specificity during speciation and contributes to species boundaries in this multi-species, synchronous-spawning, and species-rich genus.

Detecting signals of adaptive evolution in grape plastomes with a focus on the Cretaceous-Palaeogene (K/Pg) transition

BACKGROUND AND AIMS

Although plastid genes are widely used in phylogenetic studies, signals of positive selection have been scarcely investigated in the grape family. The plastomes from 91 accessions of Vitaceae were examined to understand the extent to which positive selection is present and to identify which genes are involved. Moreover, the changes through time of genes under episodic positive selection were investigated and the hypothesis of an adaptive process following the Cretaceous-Palaeogene (K/Pg) transition about 66 million years ago was tested.

METHODS

Different codon-substitution models were used to assess pervasive and episodic positive selection events on 70 candidate plastid genes. Divergence times between lineages were estimated and stochastic character mapping analysis was used to simulate variation over time of the genes found to be under episodic positive selection.

KEY RESULTS

A total of 20 plastid genes (29 %) showed positive selection. Among them, 14 genes showed pervasive signatures of positive selection and nine genes showed episodic signatures of positive selection. In particular, four of the nine genes (psbK, rpl20, rpoB, rps11) exhibited a similar pattern showing an increase in the rate of variation close to the K/Pg transition.

CONCLUSION

Multiple analyses have shown that the grape family has experienced ancient and recent positive selection events and that the targeted genes are involved in essential functions such as photosynthesis, self-replication and metabolism. Our results are consistent with the idea that the K/Pg transition has favoured an increased rate of change in some genes. Intense environmental perturbations have influenced the rapid diversification of certain lineages, and new mutations arising on some plastid genes may have been fixed by natural selection over the course of many generations.

Genomic analysis of hypoxia inducible factor alpha in ray-finned fishes reveals missing Ohnologs and evidence of widespread positive selection

As aquatic hypoxia worsens on a global scale, fishes will become increasingly challenged by low oxygen, and understanding the molecular basis of their response to hypoxia may help to better define the capacity of fishes to cope with this challenge. The hypoxia inducible factor (HIF) plays a critical role in the molecular response to hypoxia by activating the transcription of genes that serve to improve oxygen delivery to the tissues or enhance the capacity of tissues to function at low oxygen. The current study examines the molecular evolution of genes encoding the oxygen-dependent HIFα subunit (HIFA) in the ray-finned fishes (Actinopterygii). Genomic analyses demonstrate that several lineages retain four paralogs of HIFA predicted from two rounds of genome duplication at the base of vertebrate evolution, broaden the known distribution of teleost-specific HIFA paralogs, and provide evidence for salmonid-specific HIFA duplicates. Evolution of the HIFA gene family is characterized by widespread episodic positive selection at amino acid sites that potentially mediate protein stability, protein-protein interactions, and transcriptional regulation. HIFA transcript abundance depends upon paralog, tissue, and fish lineage. A phylogenetically-informed gene nomenclature is proposed along with avenues for future research on this critical family of transcription factors.

NSD1 gene evolves under episodic selection within primates and mutations of specific exons in humans cause Sotos syndrome

BACKGROUND

Modern human brains and skull shapes differ from other hominids. Brain growth disorders as micro- (ASPM, MCPH1) and macrocephaly (NFIX, GLI3) have been highlighted as relevant for the evolution in humans due to the impact in early brain development. Genes associated with macrocephaly have been reported to cause this change, for example NSD1 which causes Sotos syndrome.

RESULTS

In this study we performed a systematic literature review, located the reported variants associated to Sotos syndrome along the gene domains, compared the sequences with close primates, calculated their similarity, Ka/Ks ratios, nucleotide diversity and selection, and analyzed the sequence and structural conservation with distant primates. We aimed to understand if NSD1 in humans differs from other primates since the evolution of NSD1 has not been analyzed in primates, nor if the localization of the mutations is limited to humans. Our study found that most variations causing Sotos syndrome are in exon 19, 22 and 10. In the primate comparison we did not detect Ka/Ks ratios > 1, but a high nucleotide diversity with non-synonymous variations in exons 10, 5, 9, 11 and 23, and sites under episodic selection in exon 5 and 23, and human, macaque/colobus/tarsier/galago and tarsier/lemur/colobus. Most of the domains are conserved in distant primates with a particular progressive development from a simple PWWP1 in O. garnetti to a complex structure in Human.

CONCLUSION

NSD1 is a chromatin modifier that suggests that the selection could influence brain development during modern human evolution and is not present in other primates; however, nowadays the nucleotide diversity is associated with Sotos syndrome.

Elevated antioxidant defence in the brain of deep-diving pinnipeds

While foraging, marine mammals undertake repetitive diving bouts. When the animal surfaces, reperfusion makes oxygen readily available for the electron transport chain, which leads to increased production of reactive oxygen species and risk of oxidative damage. In blood and several tissues, such as heart, lung, muscle and kidney, marine mammals generally exhibit an elevated antioxidant defence. However, the brain, whose functional integrity is critical to survival, has received little attention. We previously observed an enhanced expression of several antioxidant genes in cortical neurons of hooded seals (Cystophora cristata). Here, we studied antioxidant gene expression and enzymatic activity in the visual cortex, cerebellum and hippocampus of harp seals (Pagophilus groenlandicus) and hooded seals. Moreover, we tested several genes for positive selection. We found that antioxidants in the first line of defence, such as superoxide dismutase (SOD), glutathione peroxidase (GPX) and glutathione (GSH) were constitutively enhanced in the seal brain compared to mice (Mus musculus), whereas the glutaredoxin and thioredoxin systems were not. Possibly, the activity of the latter systems is stress-induced rather than constitutively elevated. Further, some, but not all members, of the glutathione-s-transferase (GST) family appear more highly expressed. We found no signatures of positive selection, indicating that sequence and function of the studied antioxidants are conserved in pinnipeds.

Cannabinoid receptor 2 evolutionary gene loss makes parrots more susceptible to neuroinflammation

In vertebrates, cannabinoids modulate neuroimmune interactions through two cannabinoid receptors (CNRs) conservatively expressed in the brain (CNR1, syn. CB1) and in the periphery (CNR2, syn. CB2). Our comparative genomic analysis indicates several evolutionary losses in the CNR2 gene that is involved in immune regulation. Notably, we show that the CNR2 gene pseudogenized in all parrots (Psittaciformes). This CNR2 gene loss occurred because of chromosomal rearrangements. Our positive selection analysis suggests the absence of any specific molecular adaptations in parrot CNR1 that would compensate for the CNR2 loss in the modulation of the neuroimmune interactions. Using transcriptomic data from the brains of birds with experimentally induced sterile inflammation we highlight possible functional effects of such a CNR2 gene loss. We compare the expression patterns of CNR and neuroinflammatory markers in CNR2-deficient parrots (represented by the budgerigar, Melopsittacus undulatus and five other parrot species) with CNR2-intact passerines (represented by the zebra finch, Taeniopygia guttata). Unlike in passerines, stimulation with lipopolysaccharide resulted in neuroinflammation in the parrots linked with a significant upregulation of expression in proinflammatory cytokines (including interleukin 1 beta (IL1B) and 6 (IL6)) in the brain. Our results indicate the functional importance of the CNR2 gene loss for increased sensitivity to brain inflammation.

Bayesian Molecular Dating Analyses Combined with Mutational Profiling Suggest an Independent Origin and Evolution of SARS-CoV-2 Omicron BA.1 and BA.2 Sub-Lineages

The ongoing evolution of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) has resulted in the recent emergence of a highly divergent variant of concern (VOC) defined as Omicron or B.1.1.529. This VOC is of particular concern because it has the potential to evade most therapeutic antibodies and has undergone a sustained genetic evolution, resulting in the emergence of five distinct sub-lineages. However, the evolutionary dynamics of the initially identified Omicron BA.1 and BA.2 sub-lineages remain poorly understood. Herein, we combined Bayesian phylogenetic analysis, mutational profiling, and selection pressure analysis to track the virus's genetic changes that drive the early evolutionary dynamics of the Omicron. Based on the Omicron dataset chosen for the improved temporal signals and sampled globally between November 2021 and January 2022, the most recent common ancestor (tMRCA) and substitution rates for BA.1 were estimated to be that of 18 September 2021 (95% highest posterior density (HPD), 4 August-22 October 2021) and 1.435 × 10-3 (95% HPD  =  1.021 × 10-3 - 1.869 × 10-3) substitution/site/year, respectively, whereas 3 November 2021 (95% highest posterior density (HPD) 26 September-28 November 2021) and 1.074 × 10-3 (95% HPD  =  6.444 × 10-4 - 1.586 × 10-3) substitution/site/year were estimated for the BA.2 sub-lineage. The findings of this study suggest that the Omicron BA.1 and BA.2 sub-lineages originated independently and evolved over time. Furthermore, we identified multiple sites in the spike protein undergoing continued diversifying selection that may alter the neutralization profile of BA.1. This study sheds light on the ongoing global genomic surveillance and Bayesian molecular dating analyses to better understand the evolutionary dynamics of the virus and, as a result, mitigate the impact of emerging variants on public health.

Jack of all trades: Genome assembly of Wild Jack and comparative genomics of Artocarpus

Artocarpus (Moraceae), known as breadfruits for their diverse nutritious fruits, is prized for its high-quality timber, medicinal value, and economic importance. Breadfruits are native to Southeast Asia but have been introduced to other continents. The most commonly cultivated species are Artocarpus heterophyllus (Jackfruit) and Artocarpus altilis (Breadfruit). With numerous smaller but nutritionally comparable fruits on a larger tree, Artocarpus hirsutus, also called "Wild Jack" or "Ayani", is an elusive forest species endemic to Indian Western Ghats. In this study, we sequenced and assembled the whole genome of Artocarpus hirsutus sampled from the sacred groves of Coorg, India. To decipher demographic and evolutionary history, we compared our Wild Jack genome with previously published Jackfruit and Breadfruit genomes. Demographic history reconstruction indicates a stronger effect of habitat rather than phylogeny on the population histories of these plants. Repetitive genomic regions, especially LTR Copia, strongly affected the demographic trajectory of A. heterophyllus. Upon further investigation, we found a recent lineage-specific accumulation of LTR Copia in A. heterophyllus, which had a major contribution to its larger genome size. Several genes from starch, sucrose metabolism, and plant hormone signal transduction pathways, in Artocarpus species had signatures of selection and gene family evolution. Our comparative genomic framework provides important insights by incorporating endemic species such as the Wild Jack.

Divergent SARS-CoV-2 variant emerges in white-tailed deer with deer-to-human transmission

Wildlife reservoirs of broad-host-range viruses have the potential to enable evolution of viral variants that can emerge to infect humans. In North America, there is phylogenomic evidence of continual transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from humans to white-tailed deer (Odocoileus virginianus) through unknown means, but no evidence of transmission from deer to humans. We carried out an observational surveillance study in Ontario, Canada during November and December 2021 (n = 300 deer) and identified a highly divergent lineage of SARS-CoV-2 in white-tailed deer (B.1.641). This lineage is one of the most divergent SARS-CoV-2 lineages identified so far, with 76 mutations (including 37 previously associated with non-human mammalian hosts). From a set of five complete and two partial deer-derived viral genomes we applied phylogenomic, recombination, selection and mutation spectrum analyses, which provided evidence for evolution and transmission in deer and a shared ancestry with mink-derived virus. Our analysis also revealed an epidemiologically linked human infection. Taken together, our findings provide evidence for sustained evolution of SARS-CoV-2 in white-tailed deer and of deer-to-human transmission.

Selection-driven adaptation to the extreme Antarctic environment in the Emperor penguin

The eco-evolutionary history of penguins is characterised by shifting from temperate to cold environments. Breeding in Antarctica, the Emperor penguin appears as an extreme outcome of this process, with unique features related to insulation, heat production and energy management. However, whether this species actually diverged from a less cold-adapted ancestor, more ecologically similar to its sister species, the King penguin, is still an open question. As the Antarctic colonisation likely resulted in vast changes in selective pressure experienced by the Emperor penguin, the relative quantification of the genomic signatures of selection, unique to each sister species, could answer this question. Applying phylogeny-based selection tests on 7651 orthologous genes, we identified a more pervasive selection shift in the Emperor penguin than in the King penguin, supporting the hypothesis that its extreme cold adaptation is a derived state. Furthermore, among candidate genes under selection, four (TRPM8, LEPR, CRB1, and SFI1) were identified before in other cold-adapted homeotherms, like the woolly Mammoth, while other 161 genes can be assigned to biological functions relevant to cold adaptation identified in previous studies. Location and structural effects of TRPM8 substitutions in Emperor and King penguin lineages support their functional role with putative divergent effects on thermal adaptation. We conclude that extreme cold adaptation in the Emperor penguin largely involved unique genetic options which, however, affect metabolic and physiological traits common to other cold-adapted homeotherms.

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.

Adaptive duplication and genetic diversification of protein kinase R contribute to the specificity of bat-virus interactions

Several bat species act as asymptomatic reservoirs for many viruses that are highly pathogenic in other mammals. Here, we have characterized the functional diversification of the protein kinase R (PKR), a major antiviral innate defense system. Our data indicate that PKR has evolved under positive selection and has undergone repeated genomic duplications in bats in contrast to all studied mammals that have a single copy of the gene. Functional testing of the relationship between PKR and poxvirus antagonists revealed how an evolutionary conflict with ancient pathogenic poxviruses has shaped a specific bat host-virus interface. We determined that duplicated PKRs of the Myotis species have undergone genetic diversification, allowing them to collectively escape from and enhance the control of DNA and RNA viruses. These findings suggest that viral-driven adaptations in PKR contribute to modern virus-bat interactions and may account for bat-specific immunity.

Positive selection-driven fixation of a hominin-specific amino acid mutation related to dephosphorylation in IRF9

The arms race between humans and pathogens drives the evolution of the human genome. It is thus expected that genes from the interferon-regulatory factors family (IRFs), a critical family for anti-viral immune response, should be undergoing episodes of positive selection. Herein, we tested this hypothesis and found multiple lines of evidence for positive selection on the amino acid site Val129 (NP_006075.3:p.Ser129Val) of human IRF9. Interestingly, the ancestral reconstruction and population distribution analyses revealed that the ancestral state (Ser129) is conserved among mammals, while the derived positively selected state (Val129) was fixed before the “out-of-Africa” event ~ 500,000 years ago. The motif analysis revealed that this young amino acid (Val129) may serve as a dephosphorylation site of IRF9. Structural parallelism between homologous genes further suggested the functional effects underlying the dephosphorylation that may affect the immune activity of IRF9. This study provides a model in which a strong positive Darwinian selection drives a recent fixation of a hominin-specific amino acid leading to molecular adaptation involving dephosphorylation in an immune-responsive gene.

Widespread convergent evolution of alpha-neurotoxin resistance in African mammals

Convergent evolution is central to the study of adaptation and has been used to understand both the limits of evolution and the diverse patterns and processes which result in adaptive change. Resistance to snake venom alpha-neurotoxins (αNTXs) is a case of widespread convergence having evolved several times in snakes, lizards and mammals. Despite extreme toxicity of αNTXs, substitutions in its target, the nicotinic acetylcholine receptor (nAChR), prevent αNTX binding and render species resistant. Recently, the published meerkat (Herpestidae) genome revealed that meerkats have the same substitutions in nAChR as the venom-resistant Egyptian mongoose (Herpestidae), suggesting that venom-resistant nAChRs may be ancestral to Herpestids. Like the mongoose, many other species of feliform carnivores prey on venomous snakes, though their venom resistance has never been explored. To evaluate the prevalence and ancestry of αNTX resistance in mammals, we generate a dataset of mammalian nAChR using museum specimens and public datasets. We find five instances of convergent evolution within feliform carnivores, and an additional eight instances across all mammals sampled. Tests of selection show that these substitutions are evolving under positive selection. Repeated convergence suggests that this adaptation played an important role in the evolution of mammalian physiology and potentially venom evolution.

Lineage-specific protein repeat expansions and contractions reveal malleable regions of immune genes

Functional diversification, a higher evolutionary rate, and intense positive selection help a limited number of immune genes interact with many pathogens. Repeats in protein-coding regions are a well-known source of functional diversification, adaptive variation, and evolutionary novelty in a short time. Repeats play a crucial role in biochemical functions like functional diversification of transcription regulation, protein kinases, cell adhesion, signaling pathways, morphogenesis, DNA repair, recombination, and RNA processing. Repeat length variation can change the associated protein's interaction, efficacy, and overall protein network. Repeats have an intrinsic unstable nature and can potentially evolve rapidly and expedite the acquisition of complex phenotypic traits and functions. Because of their ability to generate rapid, adaptive variations over short evolutionary distances, repeats are considered "tuning knobs." Repeat length variation in specific genes, like RUNX2 and ALX4, is associated with morphological and physiological changes across vertebrates. Here we study repeat length variation as a potent source of species-specific immune diversification across several clades of tetrapods. Moreover, we provide a clade-wise comprehensive list of immune genes with repeat types for future studies of morphological/evolutionary changes within species groups. We observe significant repeat length variation of FASLG and C1QC in Rodentia and Primates' contrasting species groups, respectively.

Evolutionary analyses reveal immune cell receptor GPR84 as a conserved receptor for bacteria-derived molecules

The G protein-coupled receptor 84 (GPR84) is found in immune cells and its expression is increased under inflammatory conditions. Activation of GPR84 by medium-chain fatty acids results in pro-inflammatory responses. Here, we screened available vertebrate genome data and found that GPR84 is present in vertebrates for more than 500 million years but absent in birds and a pseudogene in bats. Cloning and functional characterization of several mammalian GPR84 orthologs in combination with evolutionary and model-based structural analyses revealed evidence for positive selection of bear GPR84 orthologs. Naturally occurring human GPR84 variants are most frequent in Asian populations causing a loss of function. Further, we identified cis- and trans-2-decenoic acid, both known to mediate bacterial communication, as evolutionary highly conserved ligands. Our integrated set of approaches contributes to a comprehensive understanding of GPR84 in terms of evolutionary and structural aspects, highlighting GPR84 as a conserved immune cell receptor for bacteria-derived molecules.

Characterization and evolutionary diversification of the phospholipase D gene family in mosses

Plant phospholipase D (PLD) exerts important roles in various biological processes, such as intracellular signaling and morphological development. Our knowledge about early land plant PLDs is still underdeveloped. In this study, we identified 84 PLD genes in six mosses, i.e., Physcomitrella patens, Ceratodon purpureus, Fontinalis antipyretica, Pleurozium schreberi, Sphagnum magellanicum, and Sphagnum fallax. These PLDs were classified into four clades (I-IV). We showed that PLD underwent rapid expansion in mosses. A total of six conserved domains and two core HKD motifs were detected. Structure analysis uncovered that the moss PLDs from within a clade generally exhibited similar exon-intron organization. Cis-elements prediction and expression analyses indicated that P. patens PLDs had key roles in stress responsiveness and plant development. Particularly, about half of the P. patens PLDs (e.g., PpPLD1, PpPLD2, and PpPLD5) were differentially expressed under biotic and abiotic stresses. We also determined the expression pattern of P. patens PLD genes in various tissues and at different stages of development. Although the moss, clubmoss, liverwort, and fern PLDs evolved largely under functional constraints, we found episodic positive selection in the moss PLDs, e.g., C. purpureus PLD2 and P. patens PLD11. We infer that the evolutionary force acting on the PLDs may have facilitated moss colonization of land. Our work provides valuable insights into the diversification of moss PLD genes, and can be used for future studies of their functions.

Enterovirus 3C Protease Cleaves TRIM7 To Dampen Its Antiviral Activity

Mammalian TRIM7 is an antiviral protein that inhibits multiple human enteroviruses by degrading the viral 2BC protein. Whether TRIM7 is reciprocally targeted by enteroviruses is not known. Here, we report that the 3C protease (3Cpro) from two enteroviruses, coxsackievirus B3 (CVB3) and poliovirus, targets TRIM7 for cleavage. CVB3 3Cpro cleaves TRIM7 at glutamine 24 (Q24), resulting in a truncated TRIM7 that fails to inhibit CVB3 due to dampened E3 ubiquitin ligase activity. TRIM7 Q24 is highly conserved across mammals, except in marsupials, which instead have a naturally occurring histidine (H24) that is not subject to 3Cpro cleavage. Marsupials also express two isoforms of TRIM7, and the two proteins from koalas have distinct antiviral activities. The longer isoform contains an additional exon due to alternate splice site usage. This additional exon contains a unique 3Cpro cleavage site, suggesting that certain enteroviruses may have evolved to target marsupial TRIM7 even if the canonical Q24 is missing. Combined with computational analyses indicating that TRIM7 is rapidly evolving, our data raise the possibility that TRIM7 may be targeted by enterovirus evasion strategies and that evolution of TRIM7 across mammals may have conferred unique antiviral properties. IMPORTANCE Enteroviruses are significant human pathogens that cause viral myocarditis, pancreatitis, and meningitis. Knowing how the host controls these viruses and how the viruses may evade host restriction is important for understanding fundamental concepts in antiviral immunity and for informing potential therapeutic interventions. In this study, we demonstrate that coxsackievirus B3 uses its virally encoded protease to target the host antiviral protein TRIM7 for cleavage, suggesting a potential mechanism of viral immune evasion. We additionally show that TRIM7 has evolved in certain mammalian lineages to express protein variants with distinct antiviral activities and susceptibilities to viral protease-mediated cleavage.