A beginners guide to estimating the non-synonymous to synonymous rate ratio of all protein-coding genes in a genome

Unraveling the evolution of the ATB2 subgroup basic leucine zipper transcription factors in plants and decoding the positive effects of BdibZIP44 and BdibZIP53 on heat stress in Brachypodium distachyon

In plants, basic region/leucine zipper motif (bZIP) transcription factors (TFs) stand as pivotal regulators in a broad spectrum of developmental mechanisms and adaptive strategies against environmental pressures. However, the ancestral origins and the evolutionary progression of their functional diversity across plant species have yet to be thoroughly illuminated. This study delved into the ATB2 subgroup bZIP homologs, tracing them back to the ancestral charophyte lineage predating land plant emergence, and categorized them into four distinct phylogenetic clusters (Clades A to D). Of particular note, our findings highlighted bZIP44_GBF6 and bZIP53 orthologs as angiosperm-specific innovations, distinguished by the acquisition of novel protein motifs and an intensified regime of purifying selection, reflecting their specialized evolutionary trajectories. Through synteny analysis, we uncovered that whole-genome duplication (WGD) events, post-monocot/dicot split, have played independent yet pivotal roles in shaping the bZIP44_GBF6 and bZIP53 lineages. Furthermore, an assessment of codon usage patterns disclosed a conspicuous bias in monocots favoring G3s, C3s, Gc3s, and GC content, while demonstrating a relative avoidance of T3s, A3s, and Nc usage frequencies. Functionally, we demonstrated that BdibZIP44 and BdibZIP53, localized to the nucleus, possessed the capability to dimerize, both homotypically and heterotypically. These proteins exhibited inducible expression under heat stress conditions in Brachypodium distachyon, implicating them in thermotolerance mechanisms. Overexpression studies reinforced their positive regulatory influence on heat stress resilience by augmenting the enzymatic activity of antioxidants, including catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD), which collectively enhance the clearance of deleterious reactive oxygen species (ROS). Taken together, this research significantly advanced our understanding of the origins and the adaptive evolutionary journey of ATB2 subgroup bZIP homologs in the plant kingdom. Moreover, it elucidated the vital roles of BdibZIP44 and BdibZIP53 in orchestrating a robust defense against high-temperature stress, thereby contributing to the broader discourse on plant adaptation and survival strategies under changing climatic conditions.

The CpG Landscape of Protein Coding DNA in Vertebrates

DNA methylation has fundamental implications for vertebrate genome evolution by influencing the mutational landscape, particularly at CpG dinucleotides. Methylation-induced mutations drive a genome-wide depletion of CpG sites, creating a dinucleotide composition bias across the genome. Examination of the standard genetic code reveals CpG to be the only facultative dinucleotide; it is however unclear what specific implications CpG bias has on protein coding DNA. Here, we use theoretical considerations of the genetic code combined with empirical genome-wide analyses in six vertebrate species-human, mouse, chicken, great tit, frog, and stickleback-to investigate how CpG content is shaped and maintained in protein-coding genes. We show that protein-coding sequences consistently exhibit significantly higher CpG content than noncoding regions and demonstrate that CpG sites are enriched in genes involved in regulatory functions and stress responses, suggesting selective maintenance of CpG content in specific loci. These findings have important implications for evolutionary applications in both natural and managed populations: CpG content could serve as a genetic marker for assessing adaptive potential, while the identification of CpG-free codons provides a framework for genome optimization in breeding and synthetic biology. Our results underscore the intricate interplay between mutational biases, selection, and epigenetic regulation, offering new insights into how vertebrate genomes evolve under varying ecological and selective pressures.

Comparative analysis of rapid alkalinization factor peptide-triggered plant immunity in citrus and closely related species

Rapid alkalinization factors (RALFs) are plant-derived, cysteine-rich small peptides that play crucial roles in plant development and immunity. However, their function in citrus immunity remains unexplored. This study systematically identified RALF family members in sweet orange, a cultivated citrus species, and five closely related species known for their strong resistance to citrus diseases. Most species contained 13 RALF genes, and their corresponding RALF orthologs exhibited minimal variation in gene structure, sequence identity, conserved domains, and motifs. However, differences in cis-acting elements within promoter regions led to differences in their expression patterns under biotic stress. Compared with sweet orange, citrus-related species exhibited a stronger and more widespread induction of RALF genes. Additionally, RALF peptides, particularly RALF11, triggered more robust immune responses in these species, including ROS bursts, mitogen-activated protein kinase activation, and the upregulation of defense-related genes. Furthermore, the basal expression levels of CrRLK1Ls, the RALF receptor homologs, were higher in citrus-related species than in sweet orange. Amplifying the RALF-CrRLK1L pathway through CrRLK1L overexpression and RALF peptide application significantly enhanced sweet orange resistance to citrus canker. These findings suggest that citrus-related species have endured selective pressure, leading to greater promoter sequence diversity, which is conducive to fine-tune gene expression in response to environmental challenges.

Diversification of social complexity following a major evolutionary transition in bees

How social complexity evolved remains a long-standing enigma. In most animal groups, social complexity is typically classified into a few discrete classes. This approach is oversimplified and constrains our inference of social evolution to a narrow trajectory consisting of transitions between classes. Such categorical classifications also limit quantitative studies on the molecular and environmental drivers of social complexity. The recent accumulation of relevant quantitative data has set the stage to overcome these limitations. Here, we propose a data-driven, high-dimensional approach for studying the full diversity of social phenotypes. We curated and analyzed a comprehensive dataset encompassing 17 social traits across 80 species and studied the evolution of social complexity in bees. We found that honey bees, stingless bees, and bumble bees underwent a major evolutionary transition ∼80 mya, inconsistent with the stepwise progression of the social ladder conceptual framework. This major evolutionary transition was followed by a phase of substantial phenotypic diversification of social complexity. Other bee lineages display a continuum of social complexity, ranging from solitary to simple societies, but do not reach the levels of social complexity seen in honey bees, stingless bees, and bumble bees. Bee evolution, therefore, provides a remarkable demonstration of a macroevolutionary process in which a major transition removed biological constraints and opened novel evolutionary opportunities, driving the exploration of the landscape of social phenotypes. Our approach can be extended to incorporate additional data types and readily applied to illuminate the evolution of social complexity in other animal groups.

Bioluminescence and environmental light drive the visual evolution of deep-sea shrimp (Oplophoroidea)

Light functions as the universal language in the deep sea (>200 m). Both bioluminescent emissions and downwelling light sources dimly illuminate the water column and can drive sensory system evolution. In pelagic environments, vertically migrating animals can experience drastic changes to their lighting environment across depth, subjecting them to unique selective pressures, possibly to distinguish between changes in ambient light and bioluminescent sources. Here we show that visual opsin diversity across a group of variable vertical migrators -bioluminescent deep-sea shrimp belonging to the Superfamily Oplophoroidea- is higher among species who migrate to shallower waters with more variable light conditions. Further, we provide evidence for adaptive visual evolution among species who have evolved an additional mode of bioluminescence (photophores), including positive selection for a putative mid-wavelength sensitive opsin that may facilitate light source discrimination. Diversification of this opsin appears to play an important role in the visual ecologies of photophore-bearing shrimp with its diversification in Oplophoroidea likely playing a critical role in the fitness and evolutionary success of this group.

The lawc gene emerged de novo from conserved genomic elements and acquired a broad expression pattern in Drosophila

It has recently become evident that the de novo emergence of genes is widespread and documented for a variety of organisms. De novo genes frequently emerge in proximity to existing genes, forming gene overlaps. Here, we present an analysis of the evolutionary history of a putative de novo gene, lawc, which overlaps with the conserved Trf2 gene, which encodes a general transcription factor in Drosophila melanogaster. We demonstrate that lawc emerged approximately 68 million years ago in the 5'-untranslated region (UTR) of Trf2 and displays an extensive spatiotemporal expression pattern. One of the most remarkable features of the lawc evolutionary history is that its emergence was facilitated by the engagement of Drosophilidae-specific short, highly conserved regions located in Trf2 introns. This represents a unique example of putative de novo gene birth involving conserved DNA regions localized in introns of conserved genes. The observed lawc expression pattern may be due to the overlap of lawc with the 5'-UTR of Trf2. This study not only enriches our understanding of gene evolution but also highlights the complex interplay between genetic conservation and innovation.

Selection on synonymous codon usage in soybean (Glycine max) WRKY genes

The WRKY transcription factor gene family in soybean [Glycine max (L.) Merr.] (GmWRKY) is critical for the plant's development and stress responses. This study examines the evolutionary dynamics of the GmWRKY gene family, focusing on its synonymous codon usage bias (CUB) in a comprehensive set of 179 coding sequences. CUB was analyzed using various indices, revealing a preference for A/T-ending codons and relatively low codon bias. Codon adaptation index (CAI) analysis suggested that these genes are optimized for efficient translation despite relatively low bias, reflecting a balance between codon diversity and translation efficiency. Neutrality and NC plots indicated that selective forces dominate over mutational forces in shaping codon usage, while selection signature analysis showed purifying selection being prevalent across the gene family. However, episodic positive selection was also detected in certain clades, highlighting potential adaptive diversification in response to environmental stress. Additionally, promoter binding site analysis uncovered correlations between codon usage and transcriptional regulation, indicating a context-dependent relationship between CUB and gene expression. Phylogenetic analysis identified 11 well-supported clades in the modern GmWRKY gene family and ancestral sequence reconstruction revealed more relaxed codon preferences and reduced selection constraints in modern GmWRKY genes, potentially linked to neofunctionalization and adaptation to environmental changes. These findings provide a framework for optimizing gene expression in transgenic soybean crops with resilience. Further functional validation of positively selected genes is recommended to elucidate their role in stress responses.

Revisiting the functional annotation of TriTryp using sequence similarity tools

Trypanosomatids are the causative agents of deadly diseases in humans and livestock. Given the high phylogenetic distance of trypanosomatids from model organisms, these organisms have ample unannotated genes. Manual functional annotation is time-consuming, highlighting the importance of automated functional annotation tools. The development of automated functional tools is a hot research topic, and multiple tools have been developed for the task. PANNZER2 is an automated functional annotation tool that merely relies on the sequence similarity of the query to the annotated proteins. We tried PANNZER2 on Trypanosoma brucei, the most studied organism among trypanosomatids, to see if it could improve our knowledge of the functions of the genes. Even with the availability of automated annotation tools like InterPro2GO in databases such as TriTrypDB, PANNZER2 has made surprisingly confident predictions for some hypothetical proteins in T. brucei. In this study, we identify gaps in such annotations because of not employing pairwise sequence alignment tools in TriTrypDB's automated annotation process. Our findings demonstrate that even the use of stringent cutoffs can successfully annotate a significant number of proteins. Additionally, we discovered that adjusting the open reading frames in certain genes leads to sequences with increased sequence signature coverage-characterized by the length covered by at least one sequence signature-compared to the original sequences. This enhanced sequence signature coverage suggests these genomic fragments could be pseudogenes. To facilitate further exploration, we developed a script to help identify potential pseudogenes within an organism's genome, offering researchers a new tool for genomic analysis and understanding. We extended all our analysis to Trypanosoma cruzi and Leishmania major to assess the impact of this approach across different species. Our study demonstrates that by utilizing pairwise sequence similarity alignment, even with stringent cutoffs, we can attribute 2986, 3953, and 3798 new GO terms to the genomes of T. brucei, T. cruzi, and L. major. Additionally, we found that 210, 239, and 29 genes exhibit increased sequence signature coverage following frame correction, suggesting the presence of pseudogenes.

Number of human protein interactions correlates with structural, but not regulatory conservation of the respective genes

Introduction

The differential ratio of nonsynonymous to synonymous nucleotide substitutions (dN/dS) is a common measure of the rate of structural evolution in proteincoding genes. In addition, we recently suggested that the proportion of transposable elements in gene promoters that host functional genomic sites serves as a marker of the rate of regulatory evolution of genes. Such functional genomic regions may include transcription factor binding sites and modified histone binding loci.

Methods

Here, we constructed a model of the human interactome based on 600,136 documented molecular interactions and investigated the overall relationship between the number of interactions of each protein and the rate of structural and regulatory evolution of the corresponding genes.

Results

By evaluating a total of 4,505 human genes and 1,936 molecular pathways we found a general correlation between structural and regulatory evolution rate metrics (Spearman 0.08-0.16 and 0.25-0.37 for gene and pathway levels, respectively, p < 0.01). Further exploration revealed in the established human interactome model lack of correlation between the rate of gene regulatory evolution and the number of protein interactions on gene level, and weak negative correlation (∼0.15) on pathway level. We also found a statistically significant negative correlation between the rate of gene structural evolution and the number of protein interactions (Spearman -0.11 and -0.3 for gene and pathway levels, respectively, p < 0.01).

Discussion

Our result suggests stronger structural rather than regulatory conservation of genes whose protein products have multiple interaction partners.

Diverse signatures of convergent evolution in cactus-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently approximately 17 times. Using a machine learning-based approach, we further found that cactophily can be predicted with 76% accuracy from both functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which we found to be likely associated with altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall-degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved independently through disparate molecular mechanisms. Notably, we found that multiple cactophilic species and their close relatives have been reported as emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle-and perhaps more generally lifestyles favoring thermotolerance-might preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high-throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

Adherence and metal-ion acquisition gene expression increases during infection with Treponema phagedenis strains from bovine digital dermatitis

Digital dermatitis (DD) is an ulcerative foot lesion on the heel bulbs of dairy cattle. DD is a polymicrobial disease with no precise etiology, although Treponema spirochetes are found disproportionally abundant in diseased tissue. Within Treponema, several different species are found in DD; however, the species Treponema phagedenis is uniformly found in copious quantities and deep within the skin layers of the active, ulcerative stages of disease. The pathogenic mechanisms these bacteria use to persist in the skin and the precise role they play in the pathology of DD are widely unknown. To explore the pathogenesis and virulence of Treponema phagedenis, newly isolated strains of this species were investigated in a subcutaneous murine abscess model. In the first trial, a dosage study was conducted to compare the pathogenicity of different strains across three different treponemes per inoculum (TPI) doses based on abscess volumes. In the second trial, the expression levels of 11 putative virulence genes were obtained to gain insight into their involvement in pathogenesis. During the RT-qPCR analysis, it was determined that genes encoding for two metal-ion import lipoproteins and two adherence genes were found highly upregulated during infection. Conversely, two genes involved in motility and chemotaxis were found to not be significantly upregulated or utilized during infection. These results were supported by gene expression data from natural M2 lesions of dairy cattle. This gene expression analysis could highlight the preference in strategy for T. phagedenis to persist and adhere in the host rather than engage in motility and disseminate.

The Mitogenome of the Haecon-5 Strain of Haemonchus contortus and a Comparative Analysis of Its Nucleotide Variation with Other Laboratory Strains

Haemonchus contortus (the barber's pole worm)-a highly pathogenic gastric nematode of ruminants-causes significant economic losses in the livestock industry worldwide. H. contortus has become a valuable model organism for both fundamental and applied research (e.g., drug and vaccine discovery) because of the availability of well-defined laboratory strains (e.g., MHco3(ISE).N1 in the UK and Haecon-5 in Australia) and genomic, transcriptomic and proteomic data sets. Many recent investigations have relied heavily on the use of the chromosome-contiguous genome of MHco3(ISE).N1 in the absence of a genome for Haecon-5. However, there has been no genetic comparison of these and other strains to date. Here, we assembled and characterised the mitochondrial genome (14.1 kb) of Haecon-5 and compared it with that of MHco3(ISE).N1 and two other strains (i.e., McMaster and NZ_Hco_NP) from Australasia. We detected 276 synonymous and 25 non-synonymous single nucleotide polymorphisms (SNPs) within Haecon-5. Between the Haecon-5 and MHco3(ISE).N1 strains, we recorded 345 SNPs, 31 of which were non-synonymous and linked to fixed amino acid differences in seven protein-coding genes (nad5, nad6, nad1, atp6, nad2, cytb and nad4) between these strains. Pronounced variation (344 and 435 SNPs) was seen between Haecon-5 and each of the other two strains from Australasia. The question remains as to what impact these mitogenomic mutations might have on the biology and physiology of H. contortus, which warrants exploration. The high degree of mitogenomic variability recorded here among these strains suggests that further work should be undertaken to assess the nature and extent of the nuclear genomic variation within H. contortus.

The gradual establishment of complex coumarin biosynthetic pathway in Apiaceae

Complex coumarins (CCs) represent characteristic metabolites found in Apiaceae plants, possessing significant medical value. Their essential functional role is likely as protectants against pathogens and regulators responding to environmental stimuli. Utilizing genomes and transcriptomes from 34 Apiaceae plants, including our recently sequenced Peucedanum praeruptorum, we conduct comprehensive phylogenetic analyses to reconstruct the detailed evolutionary process of the CC biosynthetic pathway in Apiaceae. Our results show that three key enzymes - p-coumaroyl CoA 2'-hydroxylase (C2'H), C-prenyltransferase (C-PT), and cyclase - originated successively at different evolutionary nodes within Apiaceae through various means of gene duplications: ectopic and tandem duplications. Neofunctionalization endows these enzymes with novel functions necessary for CC biosynthesis, thus completing the pathway. Candidate genes are cloned for heterologous expression and subjected to in vitro enzymatic assays to test our hypothesis regarding the origins of the key enzymes, and the results precisely validate our evolutionary inferences. Among the three enzymes, C-PTs are likely the primary determinant of the structural diversity of CCs (linear/angular), due to divergent activities evolved to target different positions (C-6 or C-8) of umbelliferone. A key amino acid variation (Ala161/Thr161) is identified and proven to play a crucial role in the alteration of enzymatic activity, possibly resulting in distinct binding forms between enzymes and substrates, thereby leading to different products. In conclusion, this study provides a detailed trajectory for the establishment and evolution of the CC biosynthetic pathway in Apiaceae. It explains why only a portion, not all, of Apiaceae plants can produce CCs and reveals the mechanisms of CC structural diversity among different Apiaceae plants.

Genomic Gigantism is not Associated with Reduced Selection Efficiency in Neotropical Salamanders

Genome size variation in eukaryotes has myriad effects on organismal biology from the genomic to whole-organism level. Large genome size may be associated with lower selection efficiency because lower effective population sizes allow fixation of deleterious mutations via genetic drift, increasing genome size and decreasing selection efficiency. Because of a hypothesized negative relationship between genome size and recombination rate per base pair, increased genome size could also increase the effect of linked selection in the genome, decreasing the efficiency with which natural selection can fix or remove mutations. We used a transcriptomic dataset of 15 and a subset of six Neotropical salamander species ranging in genome size from 12 to 87 pg to study the relationship between genome size and efficiency of selection. We estimated dN/dS of salamanders with small and large genomes and tested for relaxation of selection in the larger genomes. Contrary to our expectations, we did not find a significant relationship between genome size and selection efficiency or strong evidence for higher dN/dS values in species with larger genomes for either species set. We also found little evidence for relaxation of selection in species with larger genomes. A positive correlation between genome size and range size (a proxy of population size) in this group disagrees with predictions of stronger drift in species with larger genomes. Our results highlight the complex interactions between the many forces shaping genomic variation in organisms with genomic gigantism.

Extant and extinct bilby genomes combined with Indigenous knowledge improve conservation of a unique Australian marsupial

Ninu (greater bilby, Macrotis lagotis) are desert-dwelling, culturally and ecologically important marsupials. In collaboration with Indigenous rangers and conservation managers, we generated the Ninu chromosome-level genome assembly (3.66 Gbp) and genome sequences for the extinct Yallara (lesser bilby, Macrotis leucura). We developed and tested a scat single-nucleotide polymorphism panel to inform current and future conservation actions, undertake ecological assessments and improve our understanding of Ninu genetic diversity in managed and wild populations. We also assessed the beneficial impact of translocations in the metapopulation (N = 363 Ninu). Resequenced genomes (temperate Ninu, 6; semi-arid Ninu, 6; and Yallara, 4) revealed two major population crashes during global cooling events for both species and differences in Ninu genes involved in anatomical and metabolic pathways. Despite their 45-year captive history, Ninu have fewer long runs of homozygosity than other larger mammals, which may be attributable to their boom-bust life history. Here we investigated the unique Ninu biology using 12 tissue transcriptomes revealing expression of all 115 conserved eutherian chorioallantoic placentation genes in the uterus, an XY1Y2 sex chromosome system and olfactory receptor gene expansions. Together, we demonstrate the holistic value of genomics in improving key conservation actions, understanding unique biological traits and developing tools for Indigenous rangers to monitor remote wild populations.

Strength of selection in lung tumors correlates with clinical features better than tumor mutation burden

Single nucleotide substitutions are the most common type of somatic mutations in cancer genome. The goal of this study was to use publicly available somatic mutation data to quantify negative and positive selection in individual lung tumors and test how strength of directional and absolute selection is associated with clinical features. The analysis found a significant variation in strength of selection (both negative and positive) among tumors, with median selection tending to be negative even though tumors with strong positive selection also exist. Strength of selection estimated as the density of missense mutations relative to the density of silent mutations showed only a weak correlation with tumor mutation burden. In the "all histology together" analysis we found that absolute strength of selection was strongly correlated with all clinically relevant features analyzed. In histology-stratified analysis selection was strongest in small cell lung cancer. Selection in adenocarcinoma was somewhat higher compared to squamous cell carcinoma. The study suggests that somatic mutation- based quantifying of directional and absolute selection in individual tumors can be a useful biomarker of tumor aggressiveness.

Genomic Analysis of Influenza A and B Viruses Carrying Baloxavir Resistance-Associated Substitutions Serially Passaged in Human Epithelial Cells

Baloxavir marboxil (baloxavir) is an FDA-approved inhibitor of the influenza virus polymerase acidic (PA) protein. Here, we used next-generation sequencing to compare the genomic mutational profiles of IAV H1N1 and H3N2, and IBV wild type (WT) and mutants (MUT) viruses carrying baloxavir resistance-associated substitutions (H1N1-PA I38L, I38T, and E199D; H3N2-PA I38T; and IBV-PA I38T) during passaging in normal human bronchial epithelial (NHBE) cells. We determined the ratio of nonsynonymous to synonymous nucleotide mutations (dN/dS) and identified the location and type of amino acid (AA) substitutions that occurred at a frequency of ≥30%. We observed that IAV H1N1 WT and MUT viruses remained relatively stable during passaging. While the mutational profiles for IAV H1N1 I38L, I38T, and E199D, and IBV I38T MUTs were relatively similar after each passage compared to the respective WTs, the mutational profile of the IAV H3N2 I38T MUT was significantly different for most genes compared to H3N2 WT. Our work provides insight into how baloxavir resistance-associated substitutions may impact influenza virus evolution in natural settings. Further characterization of the potentially adaptive mutations identified in this study is needed.

Evolution towards simplicity in bacterial small heat shock protein system

Evolution can tinker with multi-protein machines and replace them with simpler single-protein systems performing equivalent functions in an equally efficient manner. It is unclear how, on a molecular level, such simplification can arise. With ancestral reconstruction and biochemical analysis, we have traced the evolution of bacterial small heat shock proteins (sHsp), which help to refold proteins from aggregates using either two proteins with different functions (IbpA and IbpB) or a secondarily single sHsp that performs both functions in an equally efficient way. Secondarily single sHsp evolved from IbpA, an ancestor specialized in strong substrate binding. Evolution of an intermolecular binding site drove the alteration of substrate binding properties, as well as the formation of higher-order oligomers. Upon two mutations in the α-crystallin domain, secondarily single sHsp interacts with aggregated substrates less tightly. Paradoxically, less efficient binding positively influences the ability of sHsp to stimulate substrate refolding, since the dissociation of sHps from aggregates is required to initiate Hsp70-Hsp100-dependent substrate refolding. After the loss of a partner, IbpA took over its role in facilitating the sHsp dissociation from an aggregate by weakening the interaction with the substrate, which became beneficial for the refolding process. We show that the same two amino acids introduced in modern-day systems define whether the IbpA acts as a single sHsp or obligatorily cooperates with an IbpB partner. Our discoveries illuminate how one sequence has evolved to encode functions previously performed by two distinct proteins.

Functional and evolutionary analysis of host Synaptogyrin-2 in porcine circovirus type 2 susceptibility

Mammalian evolution has been influenced by viruses for millions of years, leaving signatures of adaptive evolution within genes encoding for viral interacting proteins. Synaptogyrin-2 (SYNGR2) is a transmembrane protein implicated in promoting bacterial and viral infections. A genome-wide association study of pigs experimentally infected with porcine circovirus type 2b (PCV2b) uncovered a missense mutation (SYNGR2 p.Arg63Cys) associated with viral load. In this study, CRISPR/Cas9-mediated gene editing of the porcine kidney 15 (PK15, wtSYNGR2+p.63Arg) cell line generated clones homozygous for the favorable SYNGR2 p.63Cys allele (emSYNGR2+p.63Cys). Infection of edited clones resulted in decreased PCV2 replication compared to wildtype PK15 (P<0.05), with consistent effects across genetically distinct PCV2b and PCV2d isolates. Sequence analyses of wild and domestic pigs (n>700) revealed the favorable SYNGR2 p.63Cys allele is unique to domestic pigs and more predominant in European than Asian breeds. A haplotype defined by the SYNGR2 p.63Cys allele was likely derived from an ancestral haplotype nearly fixed within European (0.977) but absent from Asian wild boar. We hypothesize that the SYNGR2 p.63Cys allele arose post-domestication in ancestral European swine. Decreased genetic diversity in homozygotes for the SYNGR2 p.63Cys allele compared to SYNGR2 p.63Arg, corroborates a rapid increase in frequency of SYGNR2 p.63Cys via positive selection. Signatures of adaptive evolution across mammalian species were also identified within SYNGR2 intraluminal loop domains, coinciding with the location of SYNGR2 p.Arg63Cys. Therefore, SYNGR2 may reflect a novel component of the host-virus evolutionary arms race across mammals with SYNGR2 p.Arg63Cys representing a species-specific example of putative adaptive evolution.

Mitochondrial heteroplasmic shifts reveal a positive selection of breast cancer

BACKGROUND

Breast cancer is, despite screening, not always detected early enough and is together with other tumor types known to shed genetic information in circulation. Unlike single-copy nuclear DNA, mitochondrial DNA (mtDNA) copies range from 100s to 10,000s per cell, thus providing a potentially alternative to identify potential missing cancer information in circulation at an early stage.

METHODS

To characterize mitochondrial mutation landscapes in breast cancer, whole mtDNA sequencing and bioinformatics analyses were performed on 86 breast cancer biopsies and 50 available matched baseline cancer-free whole blood samples from the same individuals, selected from a cohort of middle-aged women in Sweden. To determine whether the mutations can be detected in blood plasma prior to cancer diagnosis, we further designed a nested case-control study (n = 663) and validated the shortlisted mutations using droplet digital PCR.

RESULTS

We detected different mutation landscapes between biopsies and matched whole blood samples. Compared to whole blood samples, mtDNA from biopsies had higher heteroplasmic mutations in the D-loop region (P = 0.02), RNR2 (P = 0.005), COX1 (P = 0.037) and CYTB (P = 0.006). Furthermore, the germline mtDNA mutations had higher heteroplasmy level than the lost (P = 0.002) and de novo mutations (P = 0.04). The nonsynonymous to synonymous substitution ratio (dN/dS) was higher for the heteroplasmic mutations (P = 7.25 × 10-12) than that for the homoplasmic mutations, but the de novo (P = 0.06) and lost mutations (P = 0.03) had lower dN/dS than the germline mutations. Interestingly, we found that the critical regions for mitochondrial transcription: MT-HSP1 (odds ratio [OR]: 21.41), MT-TFH (OR: 7.70) and MT-TAS2 (OR: 3.62), had significantly higher heteroplasmic mutations than the rest of the D-loop sub-regions. Finally, we found that the presence of mt.16093T > C mutation increases 67% risk of developing breast cancer.

CONCLUSIONS

Our findings show that mitochondrial genetic landscape changes during cancer pathogenesis and positive selection of mtDNA heteroplasmic mutations in breast cancer. Most importantly, the mitochondrial mutations identified in biopsies can be traced back in matched plasma samples and could potentially be used as early breast cancer diagnostic biomarkers.

Evolution of the Tóxicos en Levadura 63 (TL63) gene family in plants and functional characterization of Arabidopsis thaliana TL63 under oxidative stress

MAIN CONCLUSION

TL63 orthologs were angiosperm specific and had undergone motifs loss and gain, and increased purifying selection. AtTL63 was involved in the response of yeast and Arabidopsis plants to oxidative stress. The Tóxicos en Levadura (TL) family, a class of E3 ubiquitin ligases with typical RING-H2 type zinc finger structure, plays a pivotal role in mediating physiological processes and responding to stress in plants. However, the evolution and function of TL63 remain unclear. In this study, TL63 homologs were dated roughly back to the origin of land plants and confirmed to have subjected to the gain and loss of motifs and increased purifying selection. Phylogenetic analysis displayed that 279 TL63s could be divided into four main clades (Clade A-D). Notably, the ancestral tandem TL40/41 cluster contributed to the expansion of modern Brassicaceae TL40/41. The substitution rate tests revealed that the TL63 lineage was evidently different from other lineages. The codon usage index exhibited that monocotyledons preferred to use not A3s and T3s, but C3s, G3s, CAI, CBI and Fop. Sequence analysis showed that the TL63 homologs had conserved TM and GLD motifs and RING-H2 domain whose key amino acid residues accounted for the high average abundance. Particularly, Arabidopsis thaliana TL63 (AtTL63) was located in the nuclei, cell membranes and peroxisomes and expressed universally and significantly throughout A. thaliana development. Under H2O2 treatment, low or moderate expression of the AtTL63 held beneficial effects on the growth and viability of yeast cells and the mutation or overexpression of the AtTL63 positively affected the growth of A. thaliana plants. In brief, this study could supply useful insight into the evolution of the plant TL63s and the AtTL63 functions under oxidative stress.

Diverse signatures of convergent evolution in cacti-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently ~17 times. Using machine-learning, we further found that cactophily can be predicted with 76% accuracy from functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which is likely associated with duplication and altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall-degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved through disparate molecular mechanisms. Remarkably, multiple cactophilic lineages and their close relatives are emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle-and perhaps more generally lifestyles favoring thermotolerance-may preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

FREEDA: An automated computational pipeline guides experimental testing of protein innovation

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that lead to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA that provides a simple graphical user interface requiring only a gene name; integrates widely used molecular evolution tools to detect positive selection in rodents, primates, carnivores, birds, and flies; and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 centromere proteins, we find statistical evidence of positive selection within loops and turns of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of mouse CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

Selection in coral mitogenomes, with insights into adaptations in the deep sea

Corals are a dominant benthic fauna that occur across a vast range of depths from just below the ocean's surface to the abyssopelagic zone. However, little is known about the evolutionary mechanisms that enable them to inhabit such a wide range of environments. The mitochondrial (mt) genome, which is involved in energetic pathways, may be subject to selection pressures at greater depths to meet the metabolic demands of that environment. Here, we use a phylogenomic framework combined with codon-based models to evaluate whether mt protein-coding genes (PCGs) associated with cellular energy functions are under positive selection across depth in three groups of corals: Octocorallia, Scleractinia, and Antipatharia. The results demonstrated that mt PCGs of deep- and shallow-water species of all three groups were primarily under strong purifying selection (0.0474 < ω < 0.3123), with the exception of positive selection in atp6 (ω = 1.3263) of deep-sea antipatharians. We also found evidence for positive selection at fifteen sites across cox1, mtMutS, and nad1 in deep-sea octocorals and nad3 of deep-sea antipatharians. These results contribute to our limited understanding of mt adaptations as a function of depth and provide insight into the molecular response of corals to the extreme deep-sea environment.

Genome-wide characterization and comparative analysis of the OSCA gene family and identification of its potential stress-responsive members in legumes

Cicer arietinum, Cajanus cajan, Vigna radiata, and Phaseolus vulgaris are economically important legume crops with high nutritional value. They are negatively impacted globally by different biotic and abiotic stresses. Hyperosmolality-gated calcium-permeable channels (OSCA) have been characterized as osmosensors in Arabidopsis thaliana but have not previously reported in legumes. This study provides a genome-wide identification, characterization, and comparative analysis of OSCA genes in legumes. Our study identified and characterized 13 OSCA genes in C. cajan, V. radiata, P. vulgaris, and 12 in C. arietinum, classified into four distinct clades. We found evidence to suggest that the OSCAs might be involved in the interaction between hormone signalling pathways and stress signalling pathways. Furthermore, they play a major role in plant growth and development. The expression levels of the OSCAs vary under different stress conditions in a tissue-specific manner. Our study can be used to develop a detailed understanding of stress regulatory mechanisms of the OSCA gene family in legumes.

FREEDA: an automated computational pipeline guides experimental testing of protein innovation by detecting positive selection

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that leads to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA (Finder of Rapidly Evolving Exons in De novo Assemblies) that provides a simple graphical user interface requiring only a gene name, integrates widely used molecular evolution tools to detect positive selection, and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 mouse centromere proteins, we find evidence of positive selection in intrinsically disordered regions of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

Evolution of the WRKY66 Gene Family and Its Mutations Generated by the CRISPR/Cas9 System Increase the Sensitivity to Salt Stress in Arabidopsis

Group Ⅲ WRKY transcription factors (TFs) play pivotal roles in responding to the diverse abiotic stress and secondary metabolism of plants. However, the evolution and function of WRKY66 remains unclear. Here, WRKY66 homologs were traced back to the origin of terrestrial plants and found to have been subjected to both motifs' gain and loss, and purifying selection. A phylogenetic analysis showed that 145 WRKY66 genes could be divided into three main clades (Clade A-C). The substitution rate tests indicated that the WRKY66 lineage was significantly different from others. A sequence analysis displayed that the WRKY66 homologs had conserved WRKY and C2HC motifs with higher proportions of crucial amino acid residues in the average abundance. The AtWRKY66 is a nuclear protein, salt- and ABA- inducible transcription activator. Simultaneously, under salt stress and ABA treatments, the superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) activities, as well as the seed germination rates of Atwrky66-knockdown plants generated by the clustered, regularly interspaced, short palindromic repeats/CRISPR-associated 9 (CRISPR/Cas9) system, were all lower than those of wild type (WT) plants, but the relative electrolyte leakage (REL) was higher, indicating the increased sensitivities of the knockdown plants to the salt stress and ABA treatments. Moreover, RNA-seq and qRT-PCR analyses revealed that several regulatory genes in the ABA-mediated signaling pathway involved in stress response of the knockdown plants were significantly regulated, being evidenced by the more moderate expressions of the genes. Therefore, the AtWRKY66 likely acts as a positive regulator in the salt stress response, which may be involved in an ABA-mediated signaling pathway.

Evolution of TOP1 and TOP1MT Topoisomerases in Chordata

Type IB topoisomerases relax the torsional stress associated with DNA metabolism in the nucleus and mitochondria and constitute important molecular targets of anticancer drugs. Vertebrates stand out among eukaryotes by having two Type IB topoisomerases acting specifically in the nucleus (TOP1) and mitochondria (TOP1MT). Despite their major importance, the origin and evolution of these paralogues remain unknown. Here, we examine the molecular evolutionary processes acting on both TOP1 and TOP1MT in Chordata, taking advantage of the increasing number of available genome sequences. We found that both TOP1 and TOP1MT evolved under strong purifying selection, as expected considering their essential biological functions. Critical active sites, including those associated with resistance to anticancer agents, were found particularly conserved. However, TOP1MT presented a higher rate of molecular evolution than TOP1, possibly related with its specialized activity on the mitochondrial genome and a less critical role in cells. We could place the duplication event that originated the TOP1 and TOP1MT paralogues early in the radiation of vertebrates, most likely associated with the first round of vertebrate tetraploidization (1R). Moreover, our data suggest that cyclostomes present a specialized mitochondrial Type IB topoisomerase. Interestingly, we identified two missense mutations replacing amino acids in the Linker region of TOP1MT in Neanderthals, which appears as a rare event when comparing the genome of both species. In conclusion, TOP1 and TOP1MT differ in their rates of evolution, and their evolutionary histories allowed us to better understand the evolution of chordates.

Does integument arise de novo or from pre-existing structures? ── Insights from the key regulatory genes controlling integument development

The origin of seeds is one of the key innovations in land plant evolution. Ovules are the developmental precursors of seeds. The integument is the envelope structure surrounding the nucellus within the ovule and developing into the seed coat when ovules mature upon fertilization. The question of whether the integument arise de novo or evolve from elaboration of pre-existing structures has caused much debate. By exploring the origin and evolution of the key regulatory genes controlling integument development and their functions during both individual and historical developmental processes, we showed the widespread presence of the homologs of ANT, CUC, BEL1, SPL, C3HDZ, INO, ATS, and ETT in seedless plant genomes. All of these genes have undergone duplication-divergence events in their evolutionary history, with most of the descendant paralogous suffering motif gain and/or loss in the coding regions. Expression and functional characterization have shown that these genes are key components of the genetic program that patterns leaf-like lateral organs. Serial homology can thus be postulated between integuments and other lateral organs in terms of the shared master regulatory genes. Given that the genetic program patterning leaf-like lateral organs formed in seedless plants, and was reused during seed origin, the integument is unlikely to arise de novo but evolved from the stem segment-specific modification of pre-existing serially homologous structures. The master 'switches' trigging the modification to specify the integument identity remain unclear. We propose a successive transformation model of integument origin.

Plastome sequencing of South American Podocarpus species reveals low rearrangement rates despite ancient gondwanan disjunctions

BACKGROUND

Historical reconstructions within Podocarpaceae have provided valuable information to disentangle biogeographic scenarios that begun 65 Mya. However, early molecular phylogenies of Podocarpaceae failed to agree on the intergeneric relationships within the family. The aims of this study were to test whether plastome organization is stable within the genus Podocarpus, to estimate the selective regimes affecting plastome protein-coding genes, and to strengthen our understanding of the phylogenetic relationships and biogeographic history.

METHODS AND RESULTS

We sequenced the plastomes of four South American species from Patagonia, southern Yungas, and Brazilian subtropical forests. We compared their plastomes to those published from Brazil, Africa, New Zealand, and Southeast Asia, along with representatives from other genera within Podocarpaceae as outgroups. The four newly sequenced plastomes ranged in size between 133,791 bp and 133,991 bp. Gene content and order among chloroplasts from South American, African and Asian Podocarpus were conserved and different from the plastome of P. totara, from New Zealand. Most genes showed substitution patterns consistent with a conservative selective regime. Phylogenies inferred from either complete sequences or protein coding regions were mostly congruent with previous studies, but showed earlier branching of P. salignus, P. totara and P. sellowii.

CONCLUSIONS

Highly similar and conserved plastomes of African, South American and Asian species suggest that P. totara plastome should be revised and compared to other species from Oceanic distribution. Furthermore, given such structural conservation, we suggest plastome sequencing is not useful to test whether genomic order can be climatically or geologically structured.

On the Evolutionary Trajectory of SARS-CoV-2: Host Immunity as a Driver of Adaptation in RNA Viruses

Host immunity can exert a complex array of selective pressures on a pathogen, which can drive highly mutable RNA viruses towards viral escape. The plasticity of a virus depends on its rate of mutation, as well as the balance of fitness cost and benefit of mutations, including viral adaptations to the host's immune response. Since its emergence, SARS-CoV-2 has diversified into genetically distinct variants, which are characterised often by clusters of mutations that bolster its capacity to escape human innate and adaptive immunity. Such viral escape is well documented in the context of other pandemic RNA viruses such as the human immunodeficiency virus (HIV) and influenza virus. This review describes the selection pressures the host's antiviral immunity exerts on SARS-CoV-2 and other RNA viruses, resulting in divergence of viral strains into more adapted forms. As RNA viruses obscure themselves from host immunity, they uncover weak points in their own armoury that can inform more comprehensive, long-lasting, and potentially cross-protective vaccine coverage.

Ten Plastomes of Crassula (Crassulaceae) and Phylogenetic Implications

The genus Crassula is the second-largest genus in the family Crassulaceae, with about 200 species. As an acknowledged super-barcode, plastomes have been extensively utilized for plant evolutionary studies. Here, we first report 10 new plastomes of Crassula. We further focused on the structural characterizations, codon usage, aversion patterns, and evolutionary rates of plastomes. The IR junction patterns-IRb had 110 bp expansion to rps19-were conservative among Crassula species. Interestingly, we found the codon usage patterns of matK gene in Crassula species are unique among Crassulaceae species with elevated ENC values. Furthermore, subgenus Crassula species have specific GC-biases in the matK gene. In addition, the codon aversion motifs from matK, pafI, and rpl22 contained phylogenetic implications within Crassula. The evolutionary rates analyses indicated all plastid genes of Crassulaceae were under the purifying selection. Among plastid genes, ycf1 and ycf2 were the most rapidly evolving genes, whereas psaC was the most conserved gene. Additionally, our phylogenetic analyses strongly supported that Crassula is sister to all other Crassulaceae species. Our findings will be useful for further evolutionary studies within the Crassula and Crassulaceae.

Demographic history and conservation genomics of caribou (Rangifer tarandus) in Québec

The loss of genetic diversity is a challenge many species are facing, with genomics being a potential tool to inform and prioritize decision-making. Most caribou (Rangifer tarandus) populations have experienced significant recent declines throughout Québec, Canada, and are considered of concern, threatened or endangered. Here, we calculated the ancestral and contemporary patterns of genomic diversity of five representative caribou populations and applied a comparative population genomics framework to assess the interplay between demographic events and genomic diversity. We first calculated a caribou specific mutation rate, μ, by extracting orthologous genes from related ungulates and estimating the rate of synonymous mutations. Whole genome re-sequencing was then completed on 67 caribou: from these data we calculated nucleotide diversity, θ π and estimated the coalescent or ancestral effective population size (N e), which ranged from 12,030 to 15,513. When compared to the census size, N C, the endangered Gaspésie Mountain caribou population had the highest ancestral N e:N C ratio which is consistent with recent work suggesting high ancestral N e:N C is of conservation concern. In contrast, values of contemporary N e, estimated from linkage-disequilibrium, ranged from 11 to 162, with Gaspésie having among the highest contemporary N e:N C ratio. Importantly, classic conservation genetics theory would predict this population to be of less concern based on this ratio. Interestingly, F varied only slightly between populations, and despite evidence of bottlenecks across the province, runs of homozygosity were not abundant in the genome. Tajima's D estimates mirrored the demographic models and current conservation status. Our study highlights how genomic patterns are nuanced and potentially misleading if viewed only through a contemporary lens; we argue a holistic conservation genomics view should integrate ancestral N e and Tajima's D into management decisions.

Comparative Analysis of Whole Chloroplast Genomes of Three Common Species of Echinochloa (Gramineae) in Paddy Fields

Echinochloa crus-galli var. crus-galli, E. crus-galli var. zelayensis, and E. glabrescens, morphologically similar at the seedling stage, are the most pernicious barnyard grass species in paddy fields worldwide. Chloroplast (cp) genomes could be conducive to their identification. In this study, we assembled the complete cp genome sequences of Echinochloa crus-galli var. crus-galli (139,856 bp), E. crus-galli var. zelayensis (139,874 bp), and E. glabrescens (139,874 bp), which exhibited a typical circular tetramerous structure, large and small single-copy regions, and a pair of inverted repeats. In Echinochloa crus-galli var. crus-galli, there were 136 simple sequence (SSRs) and 62 long (LRs) repeats, and in the other two species, 139 SSRs and 68 LRs. Each cp genome contains 92 protein-encoding genes. In Echinochloa crus-galli var. crus-galli and E. glabrescens, 321 and 1 single-nucleotide polymorphisms were detected compared to Echinochloa crus-galli var. zelayensis. IR expansion and contraction revealed small differences between the three species. The phylogenetic tree based on cp genomes demonstrated the phylogenetic relationship between ten barnyard grass species and other common Gramineae plants, showing new genetic relationships of the genus Echinochloa. This study provides valuable information on cp genomes, useful for identifying and classifying the genus Echinochloa and studying its phylogenetic relationships and evolution.

Molecular and expression analysis indicate the role of CBL interacting protein kinases (CIPKs) in abiotic stress signaling and development in chickpea

Calcineurin B-like proteins (CBL)-interacting protein kinases (CIPKs) regulate the developmental processes, hormone signal transduction and stress responses in plants. Although the genome sequence of chickpea is available, information related to the CIPK gene family is missing in this important crop plant. Here, a total of 22 CIPK genes were identified and characterized in chickpea. We found a high degree of structural and evolutionary conservation in the chickpea CIPK family. Our analysis showed that chickpea CIPKs have evolved with dicots such as Arabidopsis and soybean, and extensive gene duplication events have played an important role in the evolution and expansion of the CIPK gene family in chickpea. The three-dimensional structure of chickpea CIPKs was described by protein homology modelling. Most CIPK proteins are localized in the cytoplasm and nucleus, as predicted by subcellular localization analysis. Promoter analysis revealed various cis-regulatory elements related to plant development, hormone signaling, and abiotic stresses. RNA-seq expression analysis indicated that CIPKs are significantly expressed through a spectrum of developmental stages, tissue/organs that hinted at their important role in plant development. The qRT-PCR analysis revealed that several CaCIPK genes had specific and overlapping expressions in different abiotic stresses like drought, salt, and ABA, suggesting the important role of this gene family in abiotic stress signaling in chickpea. Thus, this study provides an avenue for detailed functional characterization of the CIPK gene family in chickpea and other legume crops.

Genome-wide identification and expression of TIFY family in cassava (Manihot esculenta Crantz)

Plant-specific TIFY [TIF(F/Y)XG] proteins serve important roles in the regulation of plant stress responses. This family encodes four subfamilies of proteins, JAZ (JASMONATE ZIM-domain), PPD (PEAPOD), ZML (Zinc-finger Inflorescence-like), and TIFY. In this work, a total of 16 JAZ, 3 PPD, 7 ZML, and 2 TIFY genes were found in cassava (Manihot esculenta Crantz) at the genome-wide level. The phylogenetics, exon-intron structure, motif organization, and conserved domains of these genes were analyzed to characterize the members of the JAZ, PPD, and ZML subfamilies. Chromosome location and synteny analyses revealed that 26 JAZ, PPD, and ZML genes were irregularly distributed across 14 of the 18 chromosomes, and 18 gene pairs were implicated in large-scale interchromosomal segmental duplication events. In addition, JAZ, PPD, and ZML gene synteny comparisons between cassava and three other plant species (Arabidopsis, Populus trichocarpa, and rice) uncovered vital information about their likely evolution. The prediction of protein interaction network and cis-acting elements reveal the function of JAZ, PPD, and ZML genes. Subsequently, expression patterns of JAZ, PPD, and ZML genes were validated by qRT-PCR as being expressed in response to osmotic, salt, and cadmium stress. Moreover, almost all JAZ subfamily genes were responsive to jasmonic acid (JA) treatment. In particular, MeJAZ1, MeJAZ13, and MeJAZ14, were highly up-regulated by three treatments, and these genes may deserve further study. This comprehensive study lays the groundwork for future research into TIFY family genes in cassava and may be valuable for genetic improvement of cassava and other related species.

The genetic mechanisms underlying the convergent evolution of pollination syndromes in the Neotropical radiation of Costus L

Selection together with variation in floral traits can act to mold floral form, often driven by a plant's predominant or most effective pollinators. To investigate the evolution of traits associated with pollination, we developed a phylogenetic framework for evaluating tempo and mode of pollination shifts across the genus Costus L., known for its evolutionary toggle between traits related to bee and bird pollination. Using a target enrichment approach, we obtained 957 loci for 171 accessions to expand the phylogenetic sampling of Neotropical Costus. In addition, we performed whole genome resequencing for a subset of 20 closely related species with contrasting pollination syndromes. For each of these 20 genomes, a high-quality assembled transcriptome was used as reference for consensus calling of candidate loci hypothesized to be associated with pollination-related traits of interest. To test for the role these candidate genes may play in evolutionary shifts in pollinators, signatures of selection were estimated as dN/dS across the identified candidate loci. We obtained a well-resolved phylogeny for Neotropical Costus despite conflict among gene trees that provide evidence of incomplete lineage sorting and/or reticulation. The overall topology and the network of genome-wide single nucleotide polymorphisms (SNPs) indicate that multiple shifts in pollination strategy have occurred across Costus, while also suggesting the presence of previously undetected signatures of hybridization between distantly related taxa. Traits related to pollination syndromes are strongly correlated and have been gained and lost in concert several times throughout the evolution of the genus. The presence of bract appendages is correlated with two traits associated with defenses against herbivory. Although labellum shape is strongly correlated with overall pollination syndrome, we found no significant impact of labellum shape on diversification rates. Evidence suggests an interplay of pollination success with other selective pressures shaping the evolution of the Costus inflorescence. Although most of the loci used for phylogenetic inference appear to be under purifying selection, many candidate genes associated with functional traits show evidence of being under positive selection. Together these results indicate an interplay of phylogenetic history with adaptive evolution leading to the diversification of pollination-associated traits in Neotropical Costus.

Complete mitochondrial genome of Pectocera sp. (Elateridae: Dendrometrinae: Oxynopterini) and its phylogenetic implications

The Elateridae family (click beetles) represents a highly diverse lineage that possesses a specialized clicking mechanism to startle predators. At present, however, phylogenetic relationships, especially among recognized subfamilies, remain contentious. Mitochondrial genomes (mitogenomes) can help resolve previously intractable phylogenetic relationships using morphological or limited molecular data. Here, we report the complete mitogenome of Pectocera sp. (Elateridae: Dendrometrinae: Oxynopterini), which was 15,962 bp in length and showed a typical gene number and order as most beetle mitogenomes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 noncoding control region (AT-rich region). Comparative genomic analyses showed a high degree of feature similarity among Pectocera sp. and other click beetles. Evolutionary analysis of all PCGs based on the nonsynonymous/synonymous substitution rate ratio (ω) indicated that cox1 and atp8 exhibited the lowest and highest evolutionary rates, respectively, and that the evolutionary rates of all PCGs, except for cox3, nad2, and nad3, were lower than the average ω of click beetles. Phylogenetic analyses based on concatenated and coalescent approaches indicated that Pectocera sp. was sister to Campsosternus auratus in the same tribe (Oxynopterini) with high support. This study offers insight into the mitogenomic basis of Pectocera sp. and provides an important data resource for exploring the taxonomy, phylogeny, and evolution of click beetles.

Metaeffector interactions modulate the type III effector-triggered immunity load of Pseudomonas syringae

The bacterial plant pathogen Pseudomonas syringae requires type III secreted effectors (T3SEs) for pathogenesis. However, a major facet of plant immunity entails the recognition of a subset of P. syringae’s T3SEs by intracellular host receptors in a process called Effector-Triggered Immunity (ETI). Prior work has shown that ETI-eliciting T3SEs are pervasive in the P. syringae species complex raising the question of how P. syringae mitigates its ETI load to become a successful pathogen. While pathogens can evade ETI by T3SE mutation, recombination, or loss, there is increasing evidence that effector-effector (a.k.a., metaeffector) interactions can suppress ETI. To study the ETI-suppression potential of P. syringae T3SE repertoires, we compared the ETI-elicitation profiles of two genetically divergent strains: P. syringae pv. tomato DC3000 (PtoDC3000) and P. syringae pv. maculicola ES4326 (PmaES4326), which are both virulent on Arabidopsis thaliana but harbour largely distinct effector repertoires. Of the 529 T3SE alleles screened on A. thaliana Col-0 from the P. syringae T3SE compendium (PsyTEC), 69 alleles from 21 T3SE families elicited ETI in at least one of the two strain backgrounds, while 50 elicited ETI in both backgrounds, resulting in 19 differential ETI responses including two novel ETI-eliciting families: AvrPto1 and HopT1. Although most of these differences were quantitative, three ETI responses were completely absent in one of the pathogenic backgrounds. We performed ETI suppression screens to test if metaeffector interactions contributed to these ETI differences, and found that HopQ1a suppressed AvrPto1m-mediated ETI, while HopG1c and HopF1g suppressed HopT1b-mediated ETI. Overall, these results show that P. syringae strains leverage metaeffector interactions and ETI suppression to overcome the ETI load associated with their native T3SE repertoires.

Prevalence, Genetics, and Evolutionary Properties of Eurasian Avian-Like H1N1 Swine Influenza Viruses in Liaoning

Swine influenza virus (SIV) is an important zoonosis pathogen. The 2009 pandemic of H1N1 influenza A virus (2009/H1N1) highlighted the importance of the role of pigs as intermediate hosts. Liaoning province, located in northeastern China, has become one of the largest pig-farming areas since 2016. However, the epidemiology and evolutionary properties of SIVs in Liaoning are largely unknown. We performed systematic epidemiological and genetic dynamics surveillance of SIVs in Liaoning province during 2020. In total, 33,195 pig nasal swabs were collected, with an SIV detection rate of 2%. Our analysis revealed that multiple subtypes of SIVs are co-circulating in the pig population in Liaoning, including H1N1, H1N2 and H3N2 SIVs. Furthermore, 24 H1N1 SIVs were confirmed to belong to the EA H1N1 lineage and divided into two genotypes. The two genotypes were both triple reassortant, and the predominant one with polymerase, nucleoprotein (NP), and matrix protein (M) genes originating from 2009/H1N1; hemagglutinin (HA) and neuraminidase (NA) genes originating from EA H1N1; and the nonstructural protein (NS) gene originating from triple reassortant H1N2 (TR H1N2) was detected in Liaoning for the first time. According to our evolutionary analysis, the EA H1N1 virus in Liaoning will undergo further genome variation.

Evolutionary History of TOPIIA Topoisomerases in Animals

TOPIIA topoisomerases are required for the regulation of DNA topology by DNA cleavage and re-ligation and are important targets of antibiotic and anticancer agents. Humans possess two TOPIIA paralogue genes (TOP2A and TOP2B) with high sequence and structural similarity but distinct cellular functions. Despite their functional and clinical relevance, the evolutionary history of TOPIIA is still poorly understood. Here we show that TOPIIA is highly conserved in Metazoa. We also found that TOPIIA paralogues from jawed and jawless vertebrates had different origins related with tetraploidization events. After duplication, TOP2B evolved under a stronger purifying selection than TOP2A, perhaps promoted by the more specialized role of TOP2B in postmitotic cells. We also detected genetic signatures of positive selection in the highly variable C-terminal domain (CTD), possibly associated with adaptation to cellular interactions. By comparing TOPIIA from modern and archaic humans, we found two amino acid substitutions in the TOP2A CTD, suggesting that TOP2A may have contributed to the evolution of present-day humans, as proposed for other cell cycle-related genes. Finally, we identified six residues conferring resistance to chemotherapy differing between TOP2A and TOP2B. These six residues could be targets for the development of TOP2A-specific inhibitors that would avoid the side effects caused by inhibiting TOP2B. Altogether, our findings clarify the origin, diversification and selection pressures governing the evolution of animal TOPIIA.

Validation of predicted anonymous proteins simply using Fisher's exact test

Motivation

Genomes sequencing has become the primary (and often the sole) experimental method to characterize newly discovered organisms, in particular from the microbial world (bacteria, archaea, viruses). This generates an ever increasing number of predicted proteins the existence of which is unwarranted, in particular among those without homolog in model organisms. As a last resort, the computation of the selection pressure from pairwise alignments of the corresponding 'Open Reading Frames' (ORFs) can be used to validate their existences. However, this approach is error-prone, as not usually associated with a significance test.

Results

We introduce the use of the straightforward Fisher's exact test as a postprocessing of the results provided by the popular CODEML sequence comparison software. The respective rates of nucleotide changes at the nonsynonymous versus synonymous position (as determined by CODEML) are turned into entries into a 2 × 2 contingency table, the probability of which is computed under the Null hypothesis that they should not behave differently if the ORFs do not encode actual proteins. Using the genome sequences of two recently isolated giant viruses, we show that strong negative selection pressures do not always provide a solid argument in favor of the existence of proteins.

The First Eight Mitogenomes of Leaf-Mining Dactylispa Beetles (Coleoptera: Chrysomelidae: Cassidinae) Shed New Light on Subgenus Relationships

The taxonomic classification of Dactylispa, a large genus of leaf-mining beetles, is problematic because it is currently based on morphology alone. Here, the first eight mitochondrial genomes of Dactylispa species, which were used to construct the first molecular phylogenies of this genus, are reported. The lengths of the eight mitogenomes range from 17,189 bp to 20,363 bp. All of the mitochondrial genomes include 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and 1 A + T-rich region. According to the nonsynonymous/synonymous mutation ratio (Ka/Ks) of all PCGs, the highest and the lowest evolutionary rates were found for atp8 and cox1, respectively, which is a common phenomenon among animals. According to relative synonymous codon usage, UUA(L) has the highest frequency. With two Gonophorini species as the outgroup, mitogenome-based phylogenetic trees of the eight Dactylispa species were constructed using maximum likelihood (ML) and Bayesian inference (BI) methods based on the PCGs, tRNAs, and rRNAs. Two DNA-based phylogenomic inferences and one protein-based phylogenomic inference support the delimitation of the subgenera Dactylispa s. str. and Platypriella as proposed in the system of Chen et al. (1986). However, the subgenus Triplispa is not recovered as monophyletic. The placement of Triplispa species requires further verification and testing with more species. We also found that both adult body shape and host plant relationship might explain the subgeneric relationships among Dactylispa beetles to a certain degree.

Molecular evolution of PCSK family: Analysis of natural selection rate and gene loss

Proprotein convertases subtilisin kexins are serine endoproteases, playing critical roles in the biological functions, including lipid, glucose, and bile acid metabolism, as well as cell proliferation, migration, and metastasis. Experimental studies have demonstrated the physiological functions of PCSKs and their association with diseases; however, studies on the evolutionary history and diversification of these proteins are missing. In the present research, a bioinformatics study was conducted on the molecular evolution of several PCSKs family members and gene loss events across placental mammalian. In order to detect evolutionary constraints and positive selection, the CodeML program of the PAML package was used. The results showed the positive selection to occur in PCSK1, PCSK3, PCSK5, and PCSK7. A decelerated rate of evolution was observed in PCSK7, PCSK3, and MBTPS1 in Carnivores compared to the rest of phylogeny, and an accelerated evolution of PCSK1, PCSK7, and MBTPS1 in Muridae family of rodents was found. Additionally, our results indicated pcsk9 gene loss in 12 species comprising Carnivores and bats (Chiroptera). Future studies are required to evaluate the functional relevance and selective evolutionary advantages associated with these modifications in PCSK proteins during evolution.

Unveiling of Evolution Pattern for HY12 Enterovirus Quasispecies and Pathogenicity Alteration

Enterovirus, like the majority of RNA viruses, evolves to survive the changeable environments by a variety of strategies. Here, we showed that HY12 virus evolved to alter its characteristics and pathogenicity by employing a non-synonymous mutation. Analyses of 5'UTR, VP1 and VP2 gene sequences revealed the existence of HY12 virus in an array of mutants defined as quasispecies. The determination of diversity and complexity showed that the mutation rate and complexity of HY12 virus quasispecies increased, while the proportion of HY12 VP1 and VP2 consensus (master) sequences decreased with increasing passages. Synonymous mutation and non-synonymous mutation analysis displayed a positive selection for HY12 quasispecies evolution. A comparison of HY12 virus in different passages demonstrated that HY12 virus altered its characteristic, phenotype, and pathogenicity via non-synonymous mutation. These findings revealed the evolution pattern for HY12 virus, and the alteration of HY12 virus characteristics and pathogenicity by mutation.

Functions of intrinsically disordered proteins through evolutionary lenses

Protein sequences are the result of an evolutionary process that involves the balancing act of experimenting with novel mutations and selecting out those that have an undesirable functional outcome. In the case of globular proteins, the function relies on a well-defined conformation, therefore, there is a strong evolutionary pressure to preserve the structure. However, different evolutionary rules might apply for the group of intrinsically disordered regions and proteins (IDR/IDPs) that exist as an ensemble of fluctuating conformations. The function of IDRs can directly originate from their disordered state or arise through different types of molecular recognition processes. There is an amazing variety of ways IDRs can carry out their functions, and this is also reflected in their evolutionary properties. In this chapter we give an overview of the different types of evolutionary behavior of disordered proteins and associated functions in normal and disease settings.

Divergent connectomic organization delineates genetic evolutionary traits in the human brain

The relationship between human brain connectomics and genetic evolutionary traits remains elusive due to the inherent challenges in combining complex associations within cerebral tissue. In this study, insights are provided about the relationship between connectomics, gene expression and divergent evolutionary pathways from non-human primates to humans. Using in vivo human brain resting-state data, we detected two co-existing idiosyncratic functional systems: the segregation network, in charge of module specialization, and the integration network, responsible for information flow. Their topology was approximated to whole-brain genetic expression (Allen Human Brain Atlas) and the co-localization patterns yielded that neuron communication functionalities-linked to Neuron Projection-were overrepresented cell traits. Homologue-orthologue comparisons using dN/dS-ratios bridged the gap between neurogenetic outcomes and biological data, summarizing the known evolutionary divergent pathways within the Homo Sapiens lineage. Evidence suggests that a crosstalk between functional specialization and information flow reflects putative biological qualities of brain architecture, such as neurite cellular functions like axonal or dendrite processes, hypothesized to have been selectively conserved in the species through positive selection. These findings expand our understanding of human brain function and unveil aspects of our cognitive trajectory in relation to our simian ancestors previously left unexplored.

Functional analysis of the odorant receptor coreceptor in odor detection in Grapholita molesta (lepidoptera: Tortricidae)

The olfactory system must detect and discriminate various semiochemicals in the environment. In response to such diversity, insects have evolved a family of odorant-gated ion channels composed of a common receptor (coreceptor, Orco) and a ligand-binding tuning odorant receptor (OR) that confers odour specificity. This study aims to examine the expression pattern of Orco gene of Grapholita molesta (GmolOrco) and to elucidate the role of GmolOrco in detecting G. molesta sex pheromone and green leaf volatiles by using gene silencing via RNA interference (RNAi) coupled antennal electrophysiological (EAG). Multiple sequence alignment showed that GmolOrco shared high sequence similarities with the Orco ortholog of lepidopterans. The results of real-time quantitative PCR detection demonstrated that GmolOrco was predominantly expressed in adult antennae and had the highest expression quantity in adult period among the different developmental stages. Compared with the noninjected controls, GmolOrco expression in GmolOrcodouble-stranded RNA (dsRNA)-injected males was reduced to 39.92% and that in females was reduced to 40.43%. EAG assays showed that the responses of GmolOrco-dsRNA injected males to sex pheromones (Z)-8-dodecenyl acetate (Z8-12:OAc) and (Z)-8-dodecenyl alcohol (Z8-12:OH) were significantly reduced, and the GmolOrco-dsRNA-injected female to green leaf volatile (Z)-3-hexenyl acetate also significantly declined. We inferred that Orco-mediated olfaction was different in male and female G. molesta adults and was mainly involved in recognizing the sex pheromones released by female moths.

Examining Signatures of Natural Selection in Antifungal Resistance Genes Across Aspergillus Fungi

Certain Aspergillus fungi cause aspergillosis, a set of diseases that typically affect immunocompromised individuals. Most cases of aspergillosis are caused by Aspergillus fumigatus, which infects millions of people annually. Some closely related so-called cryptic species, such as Aspergillus lentulus, can also cause aspergillosis, albeit at lower frequencies, and they are also clinically relevant. Few antifungal drugs are currently available for treating aspergillosis and there is increasing worldwide concern about the presence of antifungal drug resistance in Aspergillus species. Furthermore, isolates from both A. fumigatus and other Aspergillus pathogens exhibit substantial heterogeneity in their antifungal drug resistance profiles. To gain insights into the evolution of antifungal drug resistance genes in Aspergillus, we investigated signatures of positive selection in 41 genes known to be involved in drug resistance across 42 susceptible and resistant isolates from 12 Aspergillus section Fumigati species. Using codon-based site models of sequence evolution, we identified ten genes that contain 43 sites with signatures of ancient positive selection across our set of species. None of the sites that have experienced positive selection overlap with sites previously reported to be involved in drug resistance. These results identify sites that likely experienced ancient positive selection in Aspergillus genes involved in resistance to antifungal drugs and suggest that historical selective pressures on these genes likely differ from any current selective pressures imposed by antifungal drugs.

Comparative Genomic Analysis Reveals Genetic Variation and Adaptive Evolution in the Pathogenicity-Related Genes of Phytophthora capsici

Phytophthora capsici is an oomycete pathogen responsible for damping off, root rot, fruit rot, and foliar blight in popular vegetable and legume crops. The existence of distinct aggressiveness levels and physiological races among the P. capsici population is a major constraint to developing resistant varieties of host crops. In the present study, we compared the genomes of three P. capsici isolates with different aggressiveness levels to reveal their genomic differences. We obtained genome sequences using short-read and long-read technologies, which yielded an average genome size of 76 Mbp comprising 514 contigs and 15,076 predicted genes. A comparative genomic analysis uncovered the signatures of accelerated evolution, gene family expansions in the pathogenicity-related genes among the three isolates. Resequencing two additional P. capsici isolates enabled the identification of average 1,023,437 SNPs, revealing the frequent accumulation of non-synonymous substitutions in pathogenicity-related gene families. Furthermore, pathogenicity-related gene families, cytoplasmic effectors and ATP binding cassette (ABC) transporters, showed expansion signals in the more aggressive isolates, with a greater number of non-synonymous SNPs. This genomic information explains the plasticity, difference in aggressiveness levels, and genome structural variation among the P. capsici isolates, providing insight into the genomic features related to the evolution and pathogenicity of this oomycete pathogen.