The potentiality of bacteria to drive SARS-CoV-2 mutation

A recent study published in mBio by Cao et al. revealed the crucial roles of bacteria in benefitting SARS-CoV-2 mutations (B. Cao, X. Wang, W. Yin, Z. Gao, and B. Xia, mBio e3187-23, 2024, https://doi.org/10.1128/mbio.03187-23). Understanding the underlying mechanisms driving the evolution of SARS-CoV-2 is crucial for predicting the future trajectory of the COVID-19 pandemic and developing preventive and treatment strategies. This study provides important insights into the rapid and complex evolution of viruses facilitated by bacterial-virus interactions.

Persistence and evolution of Pseudomonas aeruginosa following initiation of highly effective modulator therapy in cystic fibrosis

Today, more than 90% of people with cystic fibrosis (pwCF) are eligible for the highly effective cystic fibrosis transmembrane conductance regulator (CFTR) modulator therapy called elexacaftor/tezacaftor/ivacaftor (ETI) and its use is widespread. Given the drastic respiratory symptom improvement experienced by many post-ETI, clinical studies are already underway to reduce the number of respiratory therapies, including antibiotic regimens, that pwCF historically relied on to combat lung disease progression. Early studies suggest that bacterial burden in the lungs is reduced post-ETI, yet it is unknown how chronic Pseudomonas aeruginosa populations are impacted by ETI. We found that pwCF remain infected throughout their upper and lower respiratory tract with their same strain of P. aeruginosa post-ETI, and these strains continue to evolve in response to the newly CFTR-corrected airway. Our work underscores the continued importance of CF airway microbiology in the new era of highly effective CFTR modulator therapy.

IMPORTANCE

The highly effective cystic fibrosis transmembrane conductance regulator modulator therapy Elexakaftor/Tezacaftor/Ivacaftor (ETI) has changed cystic fibrosis (CF) disease for many people with cystic fibrosis. While respiratory symptoms are improved by ETI, we found that people with CF remain infected with Pseudomonas aeruginosa. How these persistent and evolving bacterial populations will impact the clinical manifestations of CF in the coming years remains to be seen, but the role and potentially changing face of infection in CF should not be discounted in the era of highly effective modulator therapy.

Examination of the role of mutualism in immune evasion

Though the earliest stages of oncogenesis, post initiation, are not well understood, it is generally appreciated that a successful transition from a collection of dysregulated cells to an aggressive tumour requires complex ecological interactions between cancer cells and their environment. One key component of tumorigenesis is immune evasion. To investigate the interplay amongst the ecological behaviour of mutualism and immune evasion, we used a computational simulation framework. Sensitivity analyses of the growth of a virtual tumour implemented as a 2D-hexagonal lattice model suggests tumour survival depends on the interplay between growth rates, mutualism and immune evasion. In 60% of simulations, cancer clones with low growth rates, but exhibiting mutualism were able to evade the immune system and continue progressing suggesting that tumours with equivalent growth rates and no mutualism are more likely to be eliminated than tumours with mutualism. Tumours with faster growth rates showed a lower dependence upon mutualism for progression. Geostatistical analysis showed decreased spatial heterogeneity over time for polyclonal tumours with a high division rate. Overall, these results suggest that in slow growing tumours, mutualism is critical for early tumorigenesis.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Buoyancy control and air breathing in royal knifefish (Chitala blanci) and a new hypothesis for the early evolution of vertebrate air-breathing behaviors

We present the first description of inspiration-first air breaths in royal knifefish, Chitala blanci, a ray-finned fish known to use four-stroke air breaths. Four-stroke breaths are used by nearly all ray-finned fish species that use their gas bladder to breathe air and are the ancestral breath type of ray-finned fishes. Interestingly, one such species, Amia calva, is known to perform two distinct breath types. Amia use four-stroke breaths when they need more oxygen and performs inspiration-first breaths to restore buoyancy. We observed that C. blanci also performs inspiration-first breaths and tested whether the two breath types are performed for the same functions in C. blanci as they are in Amia. We recorded the frequency of each breath type when exposed to aquatic hypoxia and two conditions of oxygen availability. We found that C. blanci performed more four-stroke breaths (81% ± 15% of total breaths) than inspiration-first breaths when exposed to aerial normoxia but performed more inspiration-first breaths (72% ± 40%) than four-stroke breaths when exposed to aerial hyperoxia. These patterns match those described for Amia and indicate that C. blanci performs four-stroke breaths in response to oxygen depletion and performs inspiration-first breaths to maintain buoyancy. Few studies have examined the role of air-breathing in buoyancy regulation. Decreasing buoyancy, rather than oxygen availability, to stimulate air breaths may reveal that inspiration-first breaths are more common among fishes than we are aware. We consider this possibility and present a new hypothesis for the origin and early evolution of air breathing in vertebrates.

Macroalgal deep genomics illuminate multiple paths to aquatic, photosynthetic multicellularity

Macroalgae are multicellular, aquatic autotrophs that play vital roles in global climate maintenance and have diverse applications in biotechnology and eco-engineering, which are directly linked to their multicellularity phenotypes. However, their genomic diversity and the evolutionary mechanisms underlying multicellularity in these organisms remain uncharacterized. In this study, we sequenced 110 macroalgal genomes from diverse climates and phyla, and identified key genomic features that distinguish them from their microalgal relatives. Genes for cell adhesion, extracellular matrix formation, cell polarity, transport, and cell differentiation distinguish macroalgae from microalgae across all three major phyla, constituting conserved and unique gene sets supporting multicellular processes. Adhesome genes show phylum- and climate-specific expansions that may facilitate niche adaptation. Collectively, our study reveals genetic determinants of convergent and divergent evolutionary trajectories that have shaped morphological diversity in macroalgae and provides genome-wide frameworks to understand photosynthetic multicellular evolution in aquatic environments.

Autocatalytic Selection as a Driver for the Origin of Life

Darwin's theory of evolution by natural selection was revolutionary because it provided a mechanism by which variation could be selected. This mechanism can only operate on living systems and thus cannot be applied to the origin of life. Here, we propose a viable alternative mechanism for prebiotic systems: autocatalytic selection, in which molecules catalyze reactions and processes that lead to increases in their concentration. Crucially, this provides a driver for increases in concentrations of molecules to a level that permits prebiotic metabolism. We show how this can produce high levels of amino acids, sugar phosphates, nucleotides and lipids and then lead on to polymers. Our outline is supported by a set of guidelines to support the identification of the most likely prebiotic routes. Most of the steps in this pathway are already supported by experimental results. These proposals generate a coherent and viable set of pathways that run from established Hadean geochemistry to the beginning of life.

Early Animal Origin of BACE1 APP/Aβ Proteolytic Function

Alzheimer's disease is characterized, in part, by the accumulation of β-amyloid (Aβ) in the brain. Aβ is produced via the proteolysis of APP by BACE1 and γ-secretase. Since BACE1 is the rate-limiting enzyme in the production of Aβ, and a target for therapeutics, it is of interest to know when its proteolytic function evolved and for what purpose. Here, we take a functional evolutionary approach to show that BACE1 likely evolved from a gene duplication event near the base of the animal clade and that BACE1 APP/Aβ proteolytic function evolved during early animal diversification, hundreds of millions of years before the evolution of the APP/Aβ substrate. Our examination of BACE1 APP/Aβ proteolytic function includes cnidarians, ctenophores, and choanoflagellates. The most basal BACE1 ortholog is found in cnidarians, while ctenophores, placozoa, and choanoflagellates have genes equally orthologous to BACE1 and BACE2. BACE1 from a cnidarian (Hydra) can cleave APP to release Aβ, pushing back the date of the origin of its function to near the origin of animals. We tested more divergent BACE1/2 genes from a ctenophore (Mnemiopsis) and a choanoflagellate (Monosiga), and neither has this activity. These findings indicate that the specific proteolytic function of BACE1 evolved during the very earliest diversification of animals, most likely after a gene-duplication event.

Long-term evolution of human seasonal influenza virus A(H3N2) is associated with an increase in polymerase complex activity

Since the influenza pandemic in 1968, influenza A(H3N2) viruses have become endemic. In this state, H3N2 viruses continuously evolve to overcome immune pressure as a result of prior infection or vaccination, as is evident from the accumulation of mutations in the surface glycoproteins hemagglutinin (HA) and neuraminidase (NA). However, phylogenetic studies have also demonstrated ongoing evolution in the influenza A(H3N2) virus RNA polymerase complex genes. The RNA polymerase complex of seasonal influenza A(H3N2) viruses produces mRNA for viral protein synthesis and replicates the negative sense viral RNA genome (vRNA) through a positive sense complementary RNA intermediate (cRNA). Presently, the consequences and selection pressures driving the evolution of the polymerase complex remain largely unknown. Here, we characterize the RNA polymerase complex of seasonal influenza A(H3N2) viruses representative of nearly 50 years of influenza A(H3N2) virus evolution. The H3N2 polymerase complex is a reassortment of human and avian influenza virus genes. We show that since 1968, influenza A(H3N2) viruses have increased the transcriptional activity of the polymerase complex while retaining a close balance between mRNA, vRNA, and cRNA levels. Interestingly, the increased polymerase complex activity did not result in increased replicative ability on differentiated human airway epithelial (HAE) cells. We hypothesize that the evolutionary increase in polymerase complex activity of influenza A(H3N2) viruses may compensate for the reduced HA receptor binding and avidity that is the result of the antigenic evolution of influenza A(H3N2) viruses.

The Patterns of Codon Usage between Chordates and Arthropods are Different but Co-evolving with Mutational Biases

Different frequencies amongst codons that encode the same amino acid (i.e. synonymous codons) have been observed in multiple species. Studies focused on uncovering the forces that drive such codon usage showed that a combined effect of mutational biases and translational selection works to produce different frequencies of synonymous codons. However, only few have been able to measure and distinguish between these forces that may leave similar traces on the coding regions. Here, we have developed a codon model that allows the disentangling of mutation, selection on amino acids and synonymous codons, and GC-biased gene conversion (gBGC) which we employed on an extensive dataset of 415 chordates and 191 arthropods. We found that chordates need 15 more synonymous codon categories than arthropods to explain the empirical codon frequencies, which suggests that the extent of codon usage can vary greatly between animal phyla. Moreover, methylation at CpG sites seems to partially explain these patterns of codon usage in chordates but not in arthropods. Despite the differences between the two phyla, our findings demonstrate that in both, GC-rich codons are disfavored when mutations are GC-biased, and the opposite is true when mutations are AT-biased. This indicates that selection on the genomic coding regions might act primarily to stabilize its GC/AT content on a genome-wide level. Our study shows that the degree of synonymous codon usage varies considerably among animals, but is likely governed by a common underlying dynamic.

Crystal structure of human Cep57 C-terminal domain reveals the presence of leucine zipper and the potential microtubule binding region

Cep57, a vital centrosome-associated protein, recruits essential regulatory enzymes for centriole duplication. Its dysfunction leads to anomalies, including reduced centrioles and mosaic-variegated aneuploidy syndrome. Despite functional investigations, understanding structural aspects and their correlation with functions is partial till date. We present the structure of human Cep57 C-terminal microtubule binding (MT-BD) domain, revealing conserved motifs ensuring functional preservation across evolution. A leucine zipper, with an adjacent possible microtubule-binding region, potentially forms a stabilizing scaffold for microtubule nucleation-accommodating pulling and tension from growing microtubules. This study highlights conserved structural features of Cep57 protein, compares them with other analogous proteins, and explores how protein function is maintained across diverse organisms.

Senolytics and cell senescence: historical and evolutionary perspectives

Senolytics are a new class of anti-aging drugs developed to selectively kill 'senescent' cells that are considered harmful in normal aging. More than 20 drug trials are ongoing with diverse 'senolytic cocktails'. This commentary on recent reviews of senolytics gives a historical context of mammalian cell senescence that enabled these new drugs. While cell senescence is considered harmful to aging tissues, many studies show its essential role in some regenerative and developmental processes for which senolytic drugs may interfere. Longer-term studies of side effects are needed before senolytics are considered for general clinical practice. The wide occurrence of cell senescence in eukaryotes, yeast to fish to humans, and suggests an ancient eukaryotic process that evolved multiple phenotypes.

Threat, safety, safeness and social safeness 30 years on: Fundamental dimensions and distinctions for mental health and well-being

In 1993, the British Journal of Clinical Psychology published my paper titled 'Defence and safety: Their function in social behaviour and psychopathology'. The paper highlights that to understand people's sensitivity to threat, we also need to understand their ability to identify what is safe. This paper offers an update on these concepts, highlighting distinctions that were implicit but not clearly defined at the time. Hence, the paper seeks to clarify distinctions between: (i) threat detection and response, (ii) safety and safety seeking, (iii) safeness and (iv) their social and non-social functions and forms. Threat detection and response are to prevent or minimize harm (e.g., run from a predator or fire). Safety checking relates to monitoring for the absence and avoidance of threat, while safety seeking links to the destination of the defensive behaviour (e.g., running home). Safety seeking also relates to maintaining vigilance to the appearance of potential harms and doing things believed to avoid harm. Threat-defending and safety checking and seeking are regulated primarily through evolved threat processing systems that monitor the nature, presence, controllability and/or absence of threat (e.g., amygdala and sympathetic nervous system). Safeness uses different monitoring systems via different psychophysiological systems (e.g., prefrontal cortex, parasympathetic system) for the presence of internal and external resources that support threat-coping, risk-taking, resource exploration. Creating brain states that recruit safeness processing can impact how standard evidence-based therapies (e.g., exposure, distress tolerance and reappraisal) are experienced and produce long-term change.

Paedomorphosis and retention of juvenile diet lead speciation in a group of Neotropical snakes (Colubroides-Philodryadini)

Dipsadidae is one of the largest clades of extant reptiles, showing an impressive morphological and ecological diversity. Despite this fact, the developmental processes behind its diversity are still largely unknown. In this study, we used 3D reconstructions based on micro-CT data and geometric morphometrics to evaluate the skull morphology of Philodryas agassizii, a small, surface-dwelling dipsadid that consume spiders. Adult individuals of P. agassizii exhibit a cranial morphology frequently observed in juveniles of other surface-dwelling colubroideans, represented in our analysis by its close relative Philodryas patagoniensis. Large orbits, gibbous neurocranial roof and a relatively short jaw complex are features present in juveniles of the latter species. Furthermore, we performed an extensive survey about diet of P. patagoniensis in which we detected an ontogenetic dietary shift, indicating that arthropods are more frequently consumed by juveniles of this dietary generalist. Thus, we infer that P. agassizzii retained not only the ancestral juvenile skull morphology but also dietary preferences. This study reveals that morphological changes driven by heterochronic changes, specifically paedomorphosis, influenced the retention of ancestral life history traits in P. agassizii, and therefore promoted cladogenesis. In this way, we obtained first evidence that heterochronic processes lead speciation in the snake megadiverse clade Dipsadidae.

SAFB restricts contact domain boundaries associated with L1 chimeric transcription

Long interspersed element-1 (LINE-1 or L1) comprises 17% of the human genome, continuously generates genetic variations, and causes disease in certain cases. However, the regulation and function of L1 remain poorly understood. Here, we uncover that L1 can enrich RNA polymerase IIs (RNA Pol IIs), express L1 chimeric transcripts, and create contact domain boundaries in human cells. This impact of L1 is restricted by a nuclear matrix protein scaffold attachment factor B (SAFB) that recognizes transcriptionally active L1s by binding L1 transcripts to inhibit RNA Pol II enrichment. Acute inhibition of RNA Pol II transcription abolishes the domain boundaries associated with L1 chimeric transcripts, indicating a transcription-dependent mechanism. Deleting L1 impairs domain boundary formation, and L1 insertions during evolution have introduced species-specific domain boundaries. Our data show that L1 can create RNA Pol II-enriched regions that alter genome organization and that SAFB regulates L1 and RNA Pol II activity to preserve gene regulation.

Gene Function is a Driver of Activin Signaling Pathway Evolution Following Whole-Genome Duplication in Rainbow Trout (Oncorhynchus mykiss)

The genomes of plant and animal species are influenced by ancestral whole-genome duplication (WGD) events, which have profound impacts on the regulation and function of gene networks. To gain insight into the consequences of WGD events, we characterized the sequence conservation and expression patterns of ohnologs in the highly duplicated activin receptor signaling pathway in rainbow trout (RBT). The RBT activin receptor signaling pathway is defined by tissue-specific expression of inhibitors and ligands and broad expression of receptors and Co-Smad signaling molecules. Signaling pathway ligands exhibited shared expression, while inhibitors and Smad signaling molecules primarily express a single dominant ohnolog. Our findings suggest that gene function influences ohnolog evolution following duplication of the activin signaling pathway in RBT.

Refining taxonomic identification of microalgae through molecular and genetic evolution: a case study of Prorocentrum lima and Prorocentrum arenarium

Species delimitation based on lineage definition has become increasingly popular. However, these methods have been limited, especially for species that lack genomic data and are morphologically similar. The trickiest part for the species identification is that the interspecific and intraspecific boundaries are vague. Taking Prorocentrum (Dinophyta) as an example, analysis of cell morphology, growth, and toxin synthesis in both species of P. lima and P. arenarium does not provide a reliable basis for species delineation. However, through phylogenetic and genetic distance analyses of their ITS and LSU sequences, establishment of evolutionary tree based on orthologous gene sequences, and combining the results of automatic barcode gap discovery and Poisson tree processes models, it was sustained that P. arenarium does not belong to the P. lima complex and should be considered as an independent species. Interspecies genetic evolution analysis revealed that P. lima and P. arenarium may contribute to evolutionary direction that favors combating reverse environmental factors. In P. lima, viral invasion may be one of the reasons for its large genome size. In the study, P. lima complex has been selected as an example to enhance the taxonomic identification of microalgae through molecular and genetic evolution, offering valuable insights into refining taxonomic identification and promoting microbial biodiversity research in other species.IMPORTANCEMicroalgae, especially the species known as Prorocentrum, have received significant attention due to their ability to trigger harmful algal blooms and produce toxins. However, the boundaries between species and within species are ambiguous. Clear and comprehensive species delineation indicates that Prorocentrum arenarium should be considered as an independent species, separate from the Prorocentrum lima complex. Improving the classification and identification of microalgae through molecular and genetic evolution will provide reference points for other cryptic species. Prorocentrum occupy multiple ecological niches in marine environments, and studying their evolutionary direction contributes to understanding their ecological adaptations and community succession.

Insertion-Deletion Events Are Depleted in Protein Regions with Predicted Secondary Structure

A fundamental goal in evolutionary biology and population genetics is to understand how selection shapes the fate of new mutations. Here, we test the null hypothesis that insertion-deletion (indel) events in protein-coding regions occur randomly with respect to secondary structures. We identified indels across 11,444 sequence alignments in mouse, rat, human, chimp, and dog genomes and then quantified their overlap with four different types of secondary structure-alpha helices, beta strands, protein bends, and protein turns-predicted by deep-learning methods of AlphaFold2. Indels overlapped secondary structures 54% as much as expected and were especially underrepresented over beta strands, which tend to form internal, stable regions of proteins. In contrast, indels were enriched by 155% over regions without any predicted secondary structures. These skews were stronger in the rodent lineages compared to the primate lineages, consistent with population genetic theory predicting that natural selection will be more efficient in species with larger effective population sizes. Nonsynonymous substitutions were also less common in regions of protein secondary structure, although not as strongly reduced as in indels. In a complementary analysis of thousands of human genomes, we showed that indels overlapping secondary structure segregated at significantly lower frequency than indels outside of secondary structure. Taken together, our study shows that indels are selected against if they overlap secondary structure, presumably because they disrupt the tertiary structure and function of a protein.

Copper-based grape pest management has impacted wine aroma

Despite the high energetic cost of the reduction of sulfate to H2S, required for the synthesis of sulfur-containing amino acids, some wine Saccharomyces cerevisiae strains have been reported to produce excessive amounts of H2S during alcoholic fermentation, which is detrimental to wine quality. Surprisingly, in the presence of sulfite, used as a preservative, wine strains produce more H2S than wild (oak) or wine velum (flor) isolates during fermentation. Since copper resistance caused by the amplification of the sulfur rich protein Cup1p is a specific adaptation trait of wine strains, we analyzed the link between copper resistance mechanism, sulfur metabolism and H2S production. We show that a higher content of copper in the must increases the production of H2S, and that SO2 increases the resistance to copper. Using a set of 51 strains we observed a positive and then negative relation between the number of copies of CUP1 and H2S production during fermentation. This complex pattern could be mimicked using a multicopy plasmid carrying CUP1, confirming the relation between copper resistance and H2S production. The massive use of copper for vine sanitary management has led to the selection of resistant strains at the cost of a metabolic tradeoff: the overproduction of H2S, resulting in a decrease in wine quality.

Genome-wide patterns of non-coding and protein-coding sequence variation in the major fungal pathogen Aspergillus fumigatus

A. fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related non-pathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the non-coding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. In general, A. fumigatus non-coding regions showed higher levels of sequence variation compared to their corresponding protein-coding regions. Focusing on 2,482 genes whose protein-coding sequence identity scores ranged between 75% and 99%, we identified 478 total genes with signatures of positive selection only in their non-coding regions and 65 total genes with signatures only in their protein-coding regions. 28 of the 478 non-coding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Non-coding region variation between A. fumigatus strains included single nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that non-coding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

Update on feline calicivirus: viral evolution, pathogenesis, epidemiology, prevention and control

Feline calicivirus (FCV) is a prevalent and impactful viral pathogen affecting domestic cats. As an RNA virus, FCV exhibits high mutability and genetic plasticity, enabling its persistence within cat populations. Viral genetic diversity is associated with a broad spectrum of clinical manifestations, ranging from asymptomatic infections and mild oral and upper respiratory tract diseases to the potential development of virulent systemic, and even fatal conditions. This diversity poses distinctive challenges in diagnosis, treatment, and prevention of diseases caused by FCV. Over the past four decades, research has significantly deepened understanding of this pathogen, with an emphasis on molecular biology, evolutionary dynamics, vaccine development, and disease management strategies. This review discusses various facets of FCV, including its genomic structure, evolution, innate immunity, pathogenesis, epidemiology, and approaches to disease management. FCV remains a complex and evolving concern in feline health, requiring continuous research to enhance understanding of its genetic diversity, to improve vaccine efficacy, and to explore novel treatment options.

Multiscale selection in spatially structured populations

The spatial structure of populations is key to many (eco-)evolutionary processes. In such cases, the strength and sign of selection on a trait may depend on the spatial scale considered. An example is the evolution of altruism: selection in local environments often favours cheaters over altruists, but this can be outweighed by selection at larger scales, favouring clusters of altruists over clusters of cheaters. For populations subdivided into distinct groups, this effect is described formally by multilevel selection theory. However, many populations do not consist of non-overlapping groups but rather (self-)organize into other ecological patterns. We therefore present a mathematical framework for multiscale selection. This framework decomposes natural selection into two parts: local selection, acting within environments of a certain size, and interlocal selection, acting among them. Varying the size of the local environments subsequently allows one to measure the contribution to selection of each spatial scale. To illustrate the use of this framework, we apply it to models of the evolution of altruism and pathogen transmissibility. The analysis identifies how and to what extent ecological processes at different spatial scales contribute to selection and compete, thus providing a rigorous underpinning to eco-evolutionary intuitions.

Phylogenetic Analysis of Wall-Associated Kinase Genes in Triticum Species and Characterization of TaWAK7 Involved in Wheat Powdery Mildew Resistance

Wall-associated kinases (WAKs), a group of receptor-like kinases, have been found to play important roles in defending against pathogens and in various developmental processes. However, the importance of this family in wheat remains largely unknown. Wheat powdery mildew is caused by Blumeria graminis f. sp. tritici (Bgt), which initiates infection on the cell surface and forms haustoria inside the cell; therefore, the defense to Bgt involves extracellular and subsequently intracellular signals. In this study, WAKs were identified genome-wide and analyzed phylogenetically, and then a transmembrane WAK gene that putatively participated in pathogen-associated molecular pattern-triggered immunity and effector-triggered immunity to Bgt was functionally and evolutionarily investigated. In total, 1,193 WAKs were identified from wheat and its Gramineae relatives. Phylogenetic analysis indicated that WAKs expanded through tandem duplication or segment duplication. TaWAK7, from chromosome 2A, was identified as a Bgt-inducible gene both in susceptible and resistant materials, but it showed distinct responsive patterns. Functional analysis showed that TaWAK7 was involved in both the basal and resistance gene-mediated resistances. The specific gene structures and protein characteristics of TaWAK7, along with its orthologs, were characterized both in subgenomes of Triticum spp. and in the A genome of multiple wheat accessions, which revealed that TaWAK7 orthologs underwent complex evolution with frequent gene fusion and domain deletion. In addition, three cytoplasmic proteins interacting with TaWAK7 were indicated by yeast two-hybrid and bimolecular fluorescence complementation assays. Binding of TaWAK7 with these proteins could change its subcellular localization from the plasma membrane to the cytoplasm. This study provides a better understanding of the evolution of WAKs at the genomic level and TaWAK7 at the gene level and provides useful clues for further investigation of how WAKs transmit the extracellular signals to the cytoplasm to activate defense responses.

Immune-tumor interaction dictates spatially directed evolution of esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a poor-prognostic cancer type with extensive intra- and inter-patient heterogeneity in both genomic variations and tumor microenvironment (TME). However, the patterns and drivers of spatial genomic and microenvironmental heterogeneity of ESCC remain largely unknown. Here, we generated a spatial multi-omic atlas by whole-exome, transcriptome, and methylome sequencing of 507 tumor samples from 103 patients. We identified a novel tumor suppressor PREX2, accounting for 22% of ESCCs with frequent somatic mutations or hyper-methylation, which promoted migration and invasion of ESCC cells in vitro. Analysis of the TME and quantification of subclonal expansion indicated that ESCCs undergo spatially directed evolution, where subclones mostly originated from the tumor center but had a biased clonal expansion to the upper direction of the esophagus. Interestingly, we found upper regions of ESCCs often underwent stronger immunoediting with increased selective fitness, suggesting more stringent immune selection. In addition, distinct TMEs were associated with variable genomic and clinical outcomes. Among them, hot TME was associated with high immune evasion and subclonal heterogeneity. We also found that immunoediting, instead of CD8+ T cell abundance, acts as an independent prognostic factor of ESCCs. Importantly, we found significant heterogeneity in previously considered potential therapeutic targets, as well as BRCAness characteristics in a subset of patients, emphasizing the importance of focusing on heterogeneity in ESCC targeted therapy. Collectively, these findings provide novel insights into the mechanisms of the spatial evolution of ESCC and inform precision therapeutic strategies.