The bioinformatic analysis of cannabinoid receptors (CBRs) CB1 and CB2 reveals a detailed picture of their structure, evolution, and physiological significance within the endocannabinoid system (ECS). The study highlights the evolutionary conservation of these receptors evidenced by sequence alignments across diverse species including humans, amphibians, and fish. Both CBRs share a structural hallmark of seven transmembrane (TM) helices, characteristic of class A G-protein-coupled receptors (GPCRs), which are critical for their signalling functions. The study reports a similarity of 44.58 % between both CBR sequences, which suggests that while their evolutionary paths and physiological roles may differ, there is considerable conservation in their structures. Pathway databases like KEGG, Reactome, and WikiPathways were employed to determine the involvement of the receptors in various signalling pathways. The pathway analyses integrated within this study offer a detailed view of the CBRs interactions within a complex network of cannabinoid-related signalling pathways. High-resolution crystal structures (PDB ID: 5U09 for CB1 and 5ZTY for CB2) provided accurate structural information, showing the binding pocket volume and surface area of the receptors, essential for ligand interaction. The comparison between these receptors\' natural sequences and their engineered pseudo-CBRs (p-CBRs) showed a high degree of sequence identity, confirming the validity of using p-CBRs in receptor-ligand interaction studies. This comprehensive analysis enhances the understanding of the structural and functional dynamics of cannabinoid receptors, highlighting their physiological roles and their potential as therapeutic targets within the ECS.

The increasing resistance to polymyxins in Acinetobacter baumannii has made it even more urgent to develop new treatments. Anti-virulence compounds have been researched as a new solution. Here, we evaluated the modification of virulence features of A. baumannii after acquiring resistance to polymyxin B. The results showed lineages attaining unstable resistance to polymyxin B, except for Ab7 (A. baumannii polymyxin B resistant lineage), which showed stable resistance without an associated fitness cost. Analysis of virulence by a murine sepsis model indicated diminished virulence in Ab7 (A. baumannii polymyxin B resistant lineage) compared with Ab0 (A. baumannii polymyxin B susceptible lineage). Similarly, downregulation of virulence genes was observed by qPCR at 1 and 3 h of growth. However, an increase in bauE, abaI, and pgAB expression was observed after 6 h of growth. Comparison analysis of Ab0, Ab7, and Pseudomonas aeruginosa suggested no biofilm formation by Ab7. In general, although a decrease in virulence was observed in Ab7 when compared with Ab0, some virulence feature that enables infection could be maintained. In light of this, virulence genes bauE, abaI, and pgAB showed a potential relevance in the maintenance of virulence in polymyxin B-resistant strains, making them promising anti-virulence targets.

Phenological escape, whereby species alter the timing of life-history events to avoid seasonal antagonists, is usually analyzed either as a potential evolutionary outcome given current selection coefficients, or as a realized outcome in response to known enemies. We here gain mechanistic insights into the evolutionary trajectory of phenological escape in the brassicaceous herb Cardamine pratensis, by comparing the flowering schedules of two sympatric ecotypes in different stages of a disruptive response to egg-laying pressure imposed by the pierid butterfly Anthocharis cardamines, whose larvae are pre-dispersal seed predators (reducing realized fecundity by ~70%). When the focal point of highest intensity selection (peak egg-laying) occurs early in the flowering schedule, selection for late flowering dependent on reduced egg-laying combined with selection for early flowering dependent on reduced predator survival results in a symmetrical bimodal flowering curve; when the focal point occurs late, an asymmetrical flowering curve results with a large early flowering mode due to selection for reduced egg-laying augmented by selection for infested plants to outrun larval development and dehisce prior to seed-pod consumption. Unequal selection pressures on high and low fecundity ramets, due to asynchronous flowering and morphologically targeted (size-dependent) egg-laying, constrain phenological escape, with bimodal flowering evolving primarily in response to disruptive selection on high fecundity phenotypes. These results emphasize the importance of analyzing variation in selection coefficients among morphological phenotypes over the entire flowering schedule to predict how populations will evolve in response to altered phenologies resulting from climate change.

HIV-1 matrix protein p17 variants (vp17s), characterized by amino acid insertions at the COOH-terminal region of the viral protein, have been recently identified and studied for their biological activity. Different from their wild-type counterpart (refp17), vp17s display a potent B cell growth and clonogenic activity. Recent data have highlighted the higher prevalence of vp17s in people living with HIV-1 (PLWH) with lymphoma compared with those without lymphoma, suggesting that vp17s may play a key role in lymphomagenesis. Molecular mechanisms involved in vp17 development are still unknown. Here we assessed the efficiency of HIV-1 Reverse Transcriptase (RT) in processing this genomic region and highlighted the existence of hot spots of mutation in Gag, at the end of the matrix protein and close to the matrix-capsid junction. This is possibly due to the presence of inverted repeats and palindromic sequences together with a high content of Adenine in the 322-342 nucleotide portion, which constrain HIV-1 RT to pause on the template. To define the recombinogenic properties of hot spots of mutation in the matrix gene, we developed plasmid vectors expressing Gag and a minimally modified Gag variant, and measured homologous recombination following cell co-nucleofection by next-generation sequencing. Data obtained allowed us to show that a wide range of recombination events occur in concomitance with the identified hot spots of mutation and that imperfect events may account for vp17s generation.

Despite extensive scientific efforts directed toward the evolutionary trajectory of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans at the beginning of the COVID-19 epidemic, it remains unclear how the virus jumped into and evolved in humans so far. Herein, we recruited almost all adult coronavirus disease 2019 (COVID-19) cases appeared locally or imported from abroad during the first 8 months of the outbreak in Shanghai. From these patients, SARS-CoV-2 genomes occupying the important phylogenetic positions in the virus phylogeny were recovered. Phylogenetic and mutational landscape analyses of viral genomes recovered here and those collected in and outside of China revealed that all known SARS-CoV-2 variants exhibited the evolutionary continuity despite the co-circulation of multiple lineages during the early period of the epidemic. Various mutations have driven the rapid SARS-CoV-2 diversification, and some of them favor its better adaptation and circulation in humans, which may have determined the waxing and waning of various lineages.

Evolutionary tempo and mode summarize ancient and controversial subjects of theoretical biology such as gradualism, convergence, contingence, trends, and entrenchment. We employed an integrative methodological approach to explore the evolutionary tempo and mode of Lepidosaurian phalangeal formulae (PFs). This approach involves quantifying the frequencies of morphological changes along an evolutionary trajectory. The five meristic characters encoded by PFs are particularly valuable in revealing evolutionary patterns, owing to their discrete nature and extensive documentation in the literature. Based on a pre-existing dataset of PFs from 649 taxa (35 Lepidosauria families, including fossils), from which there exists a unique repertoire of 53 formulations, our approach simultaneously considers phenetic and phylogenetic data. This culminates in a diagram accounting for the phylogenetic dynamic of evolution traversing across different regions of morphospace. The method involves enumerating phenotypical options, reconstructing phenotypes across the phylogeny, projecting phenotypes onto a morphospace, and constructing a flow network from the frequency of evolutionary transitions between unique phenotypic conditions. This approach links Markovian chains and evolutionary trajectories to formally define parameters that describe the underlying transitions of morphological change. Among other results, we found that (a) PF evolution exhibits a clear trend towards reduction in the phalangeal count and that (b) evolutionary change tends to occur significantly between morphologically similar PFs. Notwithstanding, although minor but not trivial, transitions between distant formulas -jumps- occur. Our results support a pluralistic view including stasis, gradualism, and saltationism discriminating their prevalence in a target character evolution.

Adenosine-to-inosine (A-to-I) RNA editing, mediated by metazoan ADAR enzymes, is a prevalent post-transcriptional modification that diversifies the proteome and promotes adaptive evolution of organisms. The Drosophila Adar gene has an auto-recoding site (termed S>G site) that forms a negative-feedback loop and stabilizes the global editing activity. However, the evolutionary trajectory of Adar S>G site in many other insects remains largely unknown, preventing us from a deeper understanding on the significance of this auto-editing mechanism. In this study, we retrieved the well-annotated genomes of 375 arthropod species including the five major insect orders (Lepidoptera, Diptera, Coleoptera, Hymenoptera and Hemiptera) and several outgroup species. We performed comparative genomic analysis on the Adar auto-recoding S>G site. We found that the ancestral state of insect S>G site was an uneditable serine codon (unSer) and that this state was largely maintained in Hymenoptera. The editable serine codon (edSer) appeared in the common ancestor of Lepidoptera, Diptera and Coleoptera and was almost fixed in the three orders. Interestingly, Hemiptera species possessed comparable numbers of unSer and edSer codons, and a few \'intermediate codons\', demonstrating a multi-step evolutionary trace from unSer-to-edSer with non-synchronized mutations at three codon positions. We argue that the evolution of Adar S>G site is the best genomic evidence supporting the \'proteomic diversifying hypothesis\' of RNA editing. Our work deepens our understanding on the evolutionary significance of Adar auto-recoding site which stabilizes the global editing activity and controls transcriptomic diversity.

JUJUNCAO (Cenchrus fungigraminus; 2n = 4x = 28) is a Cenchrus grass with the highest biomass production among cultivated plants, and it can be used for mushroom cultivation, animal feed, and biofuel production. Here, we report a nearly complete genome assembly of JUJUNCAO and reveal that JUJUNCAO is an allopolyploid that originated ∼2.7 million years ago (mya). Its genome consists of two subgenomes, and subgenome A shares high collinear synteny with pearl millet. We also investigated the genome evolution of JUJUNCAO and suggest that the ancestral karyotype of Cenchrus split into the A and B ancestral karyotypes of JUJUNCAO. Comparative transcriptome and DNA methylome analyses revealed functional divergence of homeologous gene pairs between the two subgenomes, which was a further indication of asymmetric DNA methylation. The three types of centromeric repeat in the JUJUNCAO genome (CEN137, CEN148, and CEN156) may have evolved independently within each subgenome, with some introgressions of CEN156 from the B to the A subgenome. We investigated the photosynthetic characteristics of JUJUNCAO, revealing its typical C4 Kranz anatomy and high photosynthetic efficiency. NADP-ME and PEPCK appear to cooperate in the major C4 decarboxylation reaction of JUJUNCAO, which is different from other C4 photosynthetic subtypes and may contribute to its high photosynthetic efficiency and biomass yield. Taken together, our results provide insights into the highly efficient photosynthetic mechanism of JUJUNCAO and provide a valuable reference genome for future genetic and evolutionary studies, as well as genetic improvement of Cenchrus grasses.

Nonkeratinizing nasopharyngeal carcinoma (NK-NPC) has a strong association with Epstein-Barr virus (EBV) infection. The role of NK cells and the tumor cell evolutionary trajectory in NK-NPC remain unclear. In this study, we aim to investigate the function of NK cell and the evolutionary trajectory of tumor cells in NK-NPC by single-cell transcriptomic analysis, proteomics and immunohistochemistry.
NK-NPC (n = 3) and normal nasopharyngeal mucosa cases (n = 3) were collected for proteomic analysis. Single-cell transcriptomic data of NK-NPC (n = 10) and nasopharyngeal lymphatic hyperplasia (NLH, n = 3) were obtained from Gene Expression Omnibus (GSE162025 and GSE150825). Quality control, dimension reduction and clustering were based on Seurat software (v4.0.2) process and batch effects were removed by harmony (v0.1.1) software. Normal cells of nasopharyngeal mucosa and tumor cells of NK-NPC were identified using copykat software (v1.0.8). Cell-cell interactions were explored using CellChat software (v1.4.0). Tumor cell evolutionary trajectory analysis was performed using SCORPIUS software (v1.0.8). Protein and gene function enrichment analyses were performed using clusterProfiler software (v4.2.2).
A total of 161 differentially expressed proteins were obtained between NK-NPC (n = 3) and normal nasopharyngeal mucosa (n = 3) by proteomics (log2 fold change > 0.5 and P value < 0.05). Most of proteins associated with the nature killer cell mediated cytotoxicity pathway were downregulated in the NK-NPC group. In single cell transcriptomics, we identified three NK cell subsets (NK1-3), among which NK cell exhaustion was identified in the NK3 subset with high ZNF683 expression (a signature of tissue-resident NK cell) in NK-NPC. We demonstrated the presence of this ZNF683 + NK cell subset in NK-NPC but not in NLH. We also performed immunohistochemical experiments with TIGIT and LAG3 to confirm NK cell exhaustion in NK-NPC. Moreover, the trajectory analysis revealed that the evolutionary trajectory of NK-NPC tumor cells was associated with the status of EBV infection (active or latent). The analysis of cell-cell interactions uncovered a complex network of cellular interactions in NK-NPC.
This study revealed that the NK cell exhaustion might be induced by upregulation of inhibitory receptors on the surface of NK cells in NK-NPC. Treatments for the reversal of NK cell exhaustion may be a promising strategy for NK-NPC. Meanwhile, we identified a unique evolutionary trajectory of tumor cells with active status of EBV-infection in NK-NPC for the first time. Our study may provide new immunotherapeutic targets and new sight of evolutionary trajectory involving tumor genesis, development and metastasis in NK-NPC.

Diabetic nephropathy (DN) is the leading cause of death and the greatest risk to the lives of people with advanced diabetes. Yet, the molecular mechanisms underlying its development and progression remain unknown. In this research, we studied the primary pathways driving DN using transcriptome sequencing and immune repertoire analysis. Firstly, we found that the diversity and abundance of the immune repertoire in late DN were significantly increased, while there was no significant change in early DN. Furthermore, B cell-mediated antibody responses may be the leading cause of DN progression. By analyzing master regulators, we found the key DN-driving transcription factors. In the late stage of DN, immune cells, fibroblasts, and epithelial cells were abundant, but other stromal cells were few. Early DN kidneys had a higher tissue stemness score than normal and advanced DN kidneys. We showed that DN progression involves proximal tubular metabolic reprogramming and stemness restoration using Monocle3. Through WGCNA, we found that co-expression modules that regulate DN progression and immune repertoire diversity mainly regulate immune-related signaling pathways. In addition, we also found that early DN had apparent activation of immune-related signaling pathways mainly enriched in immune cells. Finally, we found that activation of fibroblasts is typical of early DN. These results provide a research basis for further exploring the molecular biology and cellular mechanisms of the occurrence and development of DN and provide a theoretical basis for the prevention and treatment of DN.