Positive Selection Research Articles

Genes involved in the limited spread of SARS-CoV-2 in the lower respiratory airways of hamsters may be associated with adaptive evolution.

Novel respiratory viruses can cause a pandemic and then evolve to coexist with humans. The Omicron strain of severe acute respiratory syndrome coronavirus 2 has spread worldwide since its emergence in late 2021, and its sub-lineages are now established in human society. Compared to previous strains, Omicron is markedly less invasive in the lungs and causes less severe disease. One reason for this is that humans are acquiring immunity through previous infection and vaccination, but the nature of the virus itself is also changing. Using our newly established low-volume inoculation system, which reflects natural human infection, we show that the Omicron strain spreads less efficiently into the lungs of hamsters compared with an earlier Wuhan strain. Furthermore, by characterizing chimeric viruses with the Omicron gene in the Wuhan strain genetic background and vice versa, we found that viral genes downstream of ORF3a, but not the S gene, were responsible for the limited spread of the Omicron strain in the lower airways of the virus-infected hamsters. Moreover, molecular evolutionary analysis of SARS-CoV-2 revealed a positive selection of genes downstream of ORF3a (M and E genes). Our findings provide insight into the adaptive evolution of the virus in humans during the pandemic convergence phase.IMPORTANCEThe severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant has spread worldwide since its emergence in late 2021, and its sub-lineages are established in human society. Compared to previous strains, the Omicron strain is less invasive in the lower respiratory tract, including the lungs, and causes less severe disease; however, the mechanistic basis for its restricted replication in the lower airways is poorly understood. In this study, using a newly established low-volume inoculation system that reflects natural human infection, we demonstrated that the Omicron strain spreads less efficiently into the lungs of hamsters compared with an earlier Wuhan strain and found that viral genes downstream of ORF3a are responsible for replication restriction in the lower respiratory tract of Omicron-infected hamsters. Furthermore, we detected a positive selection of genes downstream of ORF3a (especially the M and E genes) in SARS-CoV-2, suggesting that these genes may undergo adaptive changes in humans.

Identification of positively selected genes in Mycobacterium tuberculosis from southern Xinjiang Uygur autonomous region of China.

Tuberculosis (TB), mainly caused by Mycobacterium tuberculosis (Mtb), remains a serious public health problem. Increasing evidence supports that selective evolution is an important force affecting genomic determinants of Mtb phenotypes. It is necessary to further understand the Mtb selective evolution and identify the positively selected genes that probably drive the phenotype of Mtb. This study mainly focused on the positive selection of 807 Mtb strains from Southern Xinjiang of China using whole genome sequencing (WGS). PAML software was used for identifying the genes and sites under positive selection in 807 Mtb strains. Lineage 2 (62.70%) strains were the dominant strains in this area, followed by lineage 3 (19.45%) and lineage 4 (17.84%) strains. There were 239 codons in 47 genes under positive selection, and the genes were majorly associated with the functions of transcription, defense mechanisms, and cell wall/membrane/envelope biogenesis. There were 28 codons (43 mutations) in eight genes (gyrA, rpoB, rpoC, katG, pncA, embB, gid, and cut1) under positive selection in multi-drug resistance (MDR) strains but not in drug-susceptible (DS) strains, in which 27 mutations were drug-resistant loci, 9 mutations were non-drug-resistant loci but were in drug-resistant genes, 2 mutations were compensatory mutations, and 5 mutations were in unknown drug-resistant gene of cut1. There was a codon in Rv0336 under positive selection in L3 strains but not in L2 and L4 strains. The epitopes of T and B cells were both hyper-conserved, particularly in the T-cell epitopes. This study revealed the ongoing selective evolution of Mtb. We found some special genes and sites under positive selection which may contribute to the advantage of MDR and L3 strains. It is necessary to further study these mutations to understand their impact on phenotypes for providing more useful information to develop new TB interventions.

Machine-learning analysis reveals an important role for negative selection in shaping cancer aneuploidy landscapes

BackgroundAneuploidy, an abnormal number of chromosomes within a cell, is a hallmark of cancer. Patterns of aneuploidy differ across cancers, yet are similar in cancers affecting closely related tissues. The selection pressures underlying aneuploidy patterns are not fully understood, hindering our understanding of cancer development and progression.ResultsHere, we apply interpretable machine learning methods to study tissue-selective aneuploidy patterns. We define 20 types of features corresponding to genomic attributes of chromosome-arms, normal tissues, primary tumors, and cancer cell lines (CCLs), and use them to model gains and losses of chromosome arms in 24 cancer types. To reveal the factors that shape the tissue-specific cancer aneuploidy landscapes, we interpret the machine learning models by estimating the relative contribution of each feature to the models. While confirming known drivers of positive selection, our quantitative analysis highlights the importance of negative selection for shaping aneuploidy landscapes. This is exemplified by tumor suppressor gene density being a better predictor of gain patterns than oncogene density, and vice versa for loss patterns. We also identify the importance of tissue-selective features and demonstrate them experimentally, revealing KLF5 as an important driver for chr13q gain in colon cancer. Further supporting an important role for negative selection in shaping the aneuploidy landscapes, we find compensation by paralogs to be among the top predictors of chromosome arm loss prevalence and demonstrate this relationship for one paralog interaction. Similar factors shape aneuploidy patterns in human CCLs, demonstrating their relevance for aneuploidy research.ConclusionsOur quantitative, interpretable machine learning models improve the understanding of the genomic properties that shape cancer aneuploidy landscapes.

Positive selection underlies repeated knockout of ORF8 in SARS-CoV-2 evolution

Knockout of the ORF8 protein has repeatedly spread through the global viral population during SARS-CoV-2 evolution. Here we use both regional and global pathogen sequencing to explore the selection pressures underlying its loss. In Washington State, we identified transmission clusters with ORF8 knockout throughout SARS-CoV-2 evolution, not just on novel, high fitness viral backbones. Indeed, ORF8 is truncated more frequently and knockouts circulate for longer than for any other gene. Using a global phylogeny, we find evidence of positive selection to explain this phenomenon: nonsense mutations resulting in shortened protein products occur more frequently and are associated with faster clade growth rates than synonymous mutations in ORF8. Loss of ORF8 is also associated with reduced clinical severity, highlighting the diverse clinical impacts of SARS-CoV-2 evolution.

Under pressure: evidence for selection on color-related genes in poison frogs of the genus Ranitomeya

AbstractAposematic organisms rely on their bright conspicuous coloration to communicate to potential predators that they are toxic and unpalatable. These aposematic phenotypes are strongly tied to survival and therefore make excellent opportunities to investigate the genetic underpinning of coloration. The genus Ranitomeya includes phenotypically diverse members of Neotropical aposematic poison frogs native to South America. Significant progress has been made in elucidating the molecular mechanisms responsible for aposematic coloration in poison frogs, which have paved the way for future studies to test hypotheses of the evolution of coloration across aposematic vertebrates. However, very little is known about whether these color related genes are under positive selection. We assembled transcriptomes from publicly available data reads sets for 9 different color morphs of poison frogs in the Ranitomeya genus that display bright conspicuous coloration (four morphs of R. imitator, two morphs of R. variabilis, two morphs of R. fantastica, one morph of R. summersi) to identify protein-coding genes responsible for color production that are under positive selection. Our results show that there are multiple genes under strong positive selection that are predicted to play roles in melanin synthesis (dct, tyrp1, irf4), iridophore development (fhl1), keratin metabolism (ovol1), pteridine synthesis (prps1, xdh), and carotenoid metabolism (adh1b, aldh2). The identification of positive selection affecting candidate color-pattern genes is consistent with the possibility that these genes mediate (in part) the molecular evolution of coloration. This may be attributed to aposematic phenotypes being directly tied to survival and reproduction in poison frogs.

Chromosome-level genome assembly of Hippophae tibetana provides insights into high-altitude adaptation and flavonoid biosynthesis

BackgroundAs an endemic shrub of the Qinghai-Tibetan Plateau (QTP), the distribution of Hippophae tibetana Schlecht. ranges between 2800 and 5200 m above sea level. As the most basal branch of the Hippophae genus, H. tibetana has an extensive evolutionary history. The H. tibetana is a valuable tree for studying the ecological evolution of species under extreme conditions.ResultsHere, we generated a high-quality chromosome-level genome of H. tibetana. The total size of the assembly genome is 917 Mb. The phylogenomic analysis of 1064 single-copy genes showed a divergence between 3.4 and 12.8 Mya for H. tibetana. Multiple gene families associated with DNA repair and disease resistance were significantly expanded in H. tibetana. We also identified many genes related to DNA repair with signs of positive selection. These results showed expansion and positive selection likely play important roles in H. tibetana’s adaptation to comprehensive extreme environments in the QTP. A comprehensive genomic and transcriptomic analysis identified 49 genes involved in the flavonoid biosynthesis pathway in H. tibetana. We generated transgenic sea buckthorn hairy root producing high levels of flavonoid.ConclusionsTaken together, this H. tibetana high-quality genome provides insights into the plant adaptation mechanisms of plant under extreme environments and lay foundation for the functional genomic research and molecular breeding of H. tibetana.

Small proteins, great promises: Geographic bioprospecting of Bowman–Birk protease inhibitors and domestication side‐effects in African cowpea (Vigna unguiculata L.)

Societal Impact StatementThe legume crop cowpea is grown worldwide, but 90% of the world's total share is produced in Africa. It is a promising species due to its resilience properties, balance of macro and micronutrients and presence of health‐promoting bioactive compounds. In African countries, cowpea has a crucial role in guaranteeing food security as a subsistence crop for families and commercial income for small farmers. The discovery of compounds with high nutraceutical value and bioactive properties supports socio‐economic policies to improve health and nutrition, especially in low‐ and middle‐income countries. In turn, this encourages biodiversity protection and crop enhancement programmes.Summary Bowman–Birk protease inhibitors (BBIs) are a restricted group of small proteins in plants mainly involved in defence mechanisms against pests. BBIs are demonstrated to be active components capable of reducing the viabilities of different cancer cell lines. BBI bioactivity is directly linked to the inhibition capacity, but the variability and the efficiency against the physiological targets of different BBI isoforms remain still unexplored. We analysed the natural genetic diversity of two main genes encoding BBI trypsin‐trypsin (BBI‐TT) and trypsin‐chymotrypsin (BBI‐TC) in wild and domesticated cowpea samples mainly spread in Sub‐Saharan Africa. We analysed DNA sequences and respective amino acidic isoforms/isoproteins to explore signs of natural selection and haplotype relationships. Moreover, we calculated the binding energy between BBIs and their biological targets to identify which are the most efficient inhibitors and their geographical locations. We found a high level of haplotype diversity for both genes, almost exclusively in wild accessions and detected positive and negative selection signals in the amino acid sequences. Furthermore, in the wild diversity pool, some BBI‐TT and BBI‐TC mature proteins were potentially better interactors with the physiological targets. The long interaction between plant‐pathogen has selected new and useful isoforms in wild lineages that have allowed the chances of survival of the species to improve. On the other hand, the domestication process has produced an intense bottleneck leaving only poorly efficient BBI variants. In addition to providing information on the natural diversity and evolution of BBIs, our work discusses the potential applications in agriculture and human health.

Single-Cell DNA Sequencing Reveals an Evolutionary Pattern of CHIP in Transplant Eligible Multiple Myeloma Patients

Clonal hematopoiesis of indeterminate potential (CHIP) refers to the phenomenon where a hematopoietic stem cell acquires fitness-increasing mutation(s), resulting in its clonal expansion. CHIP is frequently observed in multiple myeloma (MM) patients, and it is associated with a worse outcome. High-throughput amplicon-based single-cell DNA sequencing was performed on circulating CD34+ cells collected from twelve MM patients before autologous stem cell transplantation (ASCT). Moreover, in four MM patients, longitudinal samples either before or post-ASCT were collected. Single-cell sequencing and data analysis were assessed using the MissionBio Tapestri® platform, with a targeted panel of 20 leukemia-associated genes. We detected CHIP pathogenic mutations in 6/12 patients (50%) at the time of transplant. The most frequently mutated genes were TET2, EZH2, KIT, DNMT3A, and ASXL1. In two patients, we observed co-occurring mutations involving an epigenetic modifier (i.e., DNMT3A) and/or a gene involved in splicing machinery (i.e., SF3B1) and/or a tyrosine kinase receptor (i.e., KIT) in the same clone. Longitudinal analysis of paired samples revealed a positive selection of mutant high-fitness clones over time, regardless of their affinity with a major or minor sub-clone. Copy number analysis of the panel of all genes did not show any numerical alterations present in stem cell compartment. Moreover, we observed a tendency of CHIP-positive patients to achieve a suboptimal response to therapy compared to those without. A sub-clone dynamic of high-fitness mutations over time was confirmed.

Chloroplast genome of Ecbolium viride (Forssk.) Alston: plastome evolution and phylogenomics of Justiceae (Acanthaceae, Acanthoideae).

Justicieae is the most taxonomically complex tribe in Acanthaceae. Here, we sequenced the plastome of Ecbolium viride, a medicinally important species. The genome was analyzed with previously reported plastome of Justiceae. The plastome of E. viride has quadripartite structure with a length of 151 185bp. The comparative genomic analyses revealed no structural inversion in Justiceae and some regions (rpoC2, ycf2, ycf1, ndhH rps16-trnQ-UGG, and trnL-CAA-ycf15) exhibiting a significant level of nucleotide divergence. The positive selection analyses revealed that some species in the tribe have undergone adaptive evolution. The visualization of the boundaries between the single copy and inverted repeat regions revealed that Justiceae chloroplast genome experienced some levels of variation, which give an insight into the evolution of the species. The longest genome was in the earliest diverged taxa of the tribe Pseuderanthemum haikangense and from this genome, a series of contraction and expansion occurred contributing to the evolution of other lineages. The plastome-based phylogeny revealed and confirmed the monophyly of Justiceae, polyphyly of Justicia and supported the tribal classification Graptophyllinae, Tetrameriinae, and Isoglossinae. We proposed that Declipterinae should be treated as subtribe and the status of Justiciinae can only be confirmed after the resolution of the polyphyletic Justicia.

Advancing the Diagnosis of Periprosthetic Joint Infections: Integrated Microfluidic Platform for Alpha-Defensins-Specific Aptamer Selection and Its Analytical Applications.

Periprosthetic joint infections (PJIs) pose a significant challenge in orthopedic surgery, particularly total joint arthroplasty (TJA), due to the potential for implant failure and increased patient morbidity. Early and accurate detection of PJIs is crucial for timely intervention and better patient prognosis. Herein, we successfully screened a high-affinity aptamer targeting alpha-defensin complex human neutrophil protein 1-3 (HNP 1-3; potential PJI biomarkers in synovial fluid [SF]) for the first time using systematic evolution of ligands by exponential enrichment (SELEX) on an integrated microfluidic platform. The compact microfluidic device enabled efficient screening, with each round completed within <2 h, comprising five rounds of positive selection, two rounds of negative selection, and one round of competitive selection. A novel one-aptamer-one-antibody assay was further developed from the optimal aptamer screened, and it could accurately quantify HNP 1-3 in SF within 3 h with only ∼50 μL of SF. The assay demonstrated strong binding affinity and specificity for the target protein in SF. Thirteen PJI SF samples were accurately diagnosed and the assay was accurate over a wide dynamic range (0.32-100 mg/L). This study has showcased a rapid and accurate diagnostic tool for PJI detection, which should see widespread use in the clinic, holding promise for potential analytical applications in orthopedic surgery and improving patient care.

Genetic characterization, structural analysis, and detection of positive selection in small heat shock proteins of Cypriniformes and Clupeiformes.

In the current investigation, a total of 42 full-length, non-redundant small heat shock proteins (sHsp) were detected in Cyprinus carpio, Labeo rohita, Danio rerio, Salmo salar, Oncorhynchus mykiss, and Clupea harengus. The sHsp genes were classified into three groups based on phylogenetic analysis. All the sHsps were shown to have higher aliphatic index values, which is an indication that these proteins are more thermally stable. The hydrophilic nature of sHsps was deduced from the fact that all fish species had negative GRAVY scores. In all of the representative fish species, sHsp genes were assigned to distinct chromosomes in an inconsistent and unequal manner. Segmental duplications are the main events that have contributed to the expansion of the sHsp genes in all species. We were also able to determine the selective pressure that was placed on particular codons and discovered several significant coding sites within the coding region of sHsps. Eventually, diversifying positive selection was found to be connected with evolutionary changes in sHsp proteins, which showed that gene evolution controlled the fish adaption event in response to environmental conditions. Clarification of the links between sHsps and environmental stress in fish will be achieved through rigorous genomic comparison, which will also yield substantial new insights.

Genomics reveals heterogeneous Plasmodium falciparum transmission and selection signals in Zambia

BackgroundGenomic surveillance is crucial for monitoring malaria transmission and understanding parasite adaptation to interventions. Zambia lacks prior nationwide efforts in malaria genomic surveillance among African countries.MethodsWe conducted genomic surveillance of Plasmodium falciparum parasites from the 2018 Malaria Indicator Survey in Zambia, a nationally representative household survey of children under five years of age. We whole-genome sequenced and analyzed 241 P. falciparum genomes from regions with varying levels of malaria transmission across Zambia and estimated genetic metrics that are informative about transmission intensity, genetic relatedness between parasites, and selection.ResultsWe provide genomic evidence of widespread within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in Zambia. Our analysis reveals country-level clustering of parasites from Zambia and neighboring regions, with distinct separation in West Africa. Within Zambia, identity by descent (IBD) relatedness analysis uncovers local spatial clustering and rare cases of long-distance sharing of closely related parasite pairs. Genomic regions with large shared IBD segments and strong positive selection signatures implicate genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Furthermore, differences in selection signatures, including drug resistance loci, are observed between eastern and western Zambian parasite populations, suggesting variable transmission intensity and ongoing drug pressure.ConclusionsOur findings enhance our understanding of nationwide P. falciparum transmission in Zambia, establishing a baseline for analyzing parasite genetic metrics as they vary over time and space. These insights highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Regulation of BCR-mediated Ca2+ mobilization by MIZ1-TMBIM4 safeguards IgG1+ GC B cell-positive selection.

The transition from immunoglobulin M (IgM) to affinity-matured IgG antibodies is vital for effective humoral immunity. This is facilitated by germinal centers (GCs) through affinity maturation and preferential maintenance of IgG+ B cells over IgM+ B cells. However, it is not known whether the positive selection of the different Ig isotypes within GCs is dependent on specific transcriptional mechanisms. Here, we explored IgG1+ GC B cell transcription factor dependency using a CRISPR-Cas9 screen and conditional mouse genetics. We found that MIZ1 was specifically required for IgG1+ GC B cell survival during positive selection, whereas IgM+ GC B cells were largely independent. Mechanistically, MIZ1 induced TMBIM4, an ancestral anti-apoptotic protein that regulated inositol trisphosphate receptor (IP3R)-mediated calcium (Ca2+) mobilization downstream of B cell receptor (BCR) signaling in IgG1+ B cells. The MIZ1-TMBIM4 axis prevented mitochondrial dysfunction-induced IgG1+ GC cell death caused by excessive Ca2+ accumulation. This study uncovers a unique Ig isotype-specific dependency on a hitherto unidentified mechanism in GC-positive selection.

Abstract SY36-01: Somatic evolution, cancer, and our inevitable decline with age: Inextricably linked

Abstract Why do we get cancer? Why is cancer highly associated with old age? Aging is associated with the accumulation of more mutations; some of which can contribute to cancer phenotypes. However, we now understand that carcinogenesis is much more complex than originally appreciated. In particular, there are tissue environmental forces that both impede and promote cancer evolution. Just as organismal evolution is known to be driven by environmental changes, cellular (somatic) evolution in our bodies is similarly driven by changes in tissue environments, whether caused by the normal process of aging, by lifestyle choices, or by extrinsic exposures. Environmental change promotes selection for new phenotypes that are adaptive to the new context. In our tissues, aging or insult-driven alterations in tissues drives selection for adaptive mutations, and some of these mutations can confer malignant phenotypes. To better understand the evolutionary forces that control somatic cell evolution, we used mouse models of cancer initiation, mathematical models of cellular evolution, and analyses of human tissue samples. We have shown that aging- and inflammation-dependent changes in stem cells and their tissue environments dramatically dictate whether cancer-causing mutations are advantageous to stem cells in our tissues, starting the cells down the path to cancer and deepening our insight into cancer risk. Recent studies from many labs have shown how as we age our tissues become dominated by clones bearing mutations in known cancer-associated genes. I will present studies from our labs exploring how cigarette smoking and aging can alter mutational landscapes in the non-malignant lung. To observe somatic mutations in lung biopsies from individuals with different smoking histories, we used a rare-mutation detection technique called Duplex Sequencing to analyze somatic variants from over 200 lung samples, including from subjects from the BDRE (Colorado), PEACE, and TRACERx cohorts. Using this method, we have identified large numbers of mutations in each brushing or biopsy, many of which are cataloged in the Catalogue Of Somatic Mutations In Cancer (COSMIC), and most of these are predicted to disrupt protein function. We also observe pervasive positive selection acting on mutations in many but not all of the cancer-associated genes, with different patterns of selection in the lungs of ever-smokers and never-smokers. Ongoing studies are using primary human and mouse lung epithelial cultures and mouse models to study how these mutations interact with smoking and age to alter associated clonal expansions, tissue integrity, and malignant progression. We propose a model whereby aging and microenvironmental changes induced through lifestyle choices like smoking can promote selection for cells with adaptive mutations, which can contribute to not only cancer risk but also tissue aging. Thus, while young tissues impede selection for such adaptive mutations, old age is associated with a feed-forward loop of aging tissue-mediated selection for mutant clones that then increase tissue aging. Thus, understanding the forces controlling clonal selection as we age and due to lifetime exposures could be critical for controlling multiple diseases of old age. [EJE, EL and FM contributed equally.These studies were led by CS and JD.] Citation Format: Edward J. Evans, Emilia Lim, Fabio Marongiu, Shi Biao Chia, William Hill, Clare Weeden, Oriol Pich, Andrew Rowan, Faiz Jabbar, Cristina Naceur-Lombardelli, Raju Veeriah, Moumita Ghosh, Daniel Merrick, York E. Miller, Robert Keith, Mariam Jamal-Hanjani, Charles Swanton, James V. DeGregori. Somatic evolution, cancer, and our inevitable decline with age: Inextricably linked [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr SY36-01.

High-Voltage Toxin'Roll: Electrostatic Charge Repulsion as a Dynamic Venom Resistance Trait in Pythonid Snakes.

The evolutionary interplay between predator and prey has significantly shaped the development of snake venom, a critical adaptation for subduing prey. This arms race has spurred the diversification of the components of venom and the corresponding emergence of resistance mechanisms in the prey and predators of venomous snakes. Our study investigates the molecular basis of venom resistance in pythons, focusing on electrostatic charge repulsion as a defense against α-neurotoxins binding to the alpha-1 subunit of the postsynaptic nicotinic acetylcholine receptor. Through phylogenetic and bioactivity analyses of orthosteric site sequences from various python species, we explore the prevalence and evolution of amino acid substitutions that confer resistance by electrostatic repulsion, which initially evolved in response to predatory pressure by Naja (cobra) species (which occurs across Africa and Asia). The small African species Python regius retains the two resistance-conferring lysines (positions 189 and 191) of the ancestral Python genus, conferring resistance to sympatric Naja venoms. This differed from the giant African species Python sebae, which has secondarily lost one of these lysines, potentially due to its rapid growth out of the prey size range of sympatric Naja species. In contrast, the two Asian species Python brongersmai (small) and Python bivittatus (giant) share an identical orthosteric site, which exhibits the highest degree of resistance, attributed to three lysine residues in the orthosteric sites. One of these lysines (at orthosteric position 195) evolved in the last common ancestor of these two species, which may reflect an adaptive response to increased predation pressures from the sympatric α-neurotoxic snake-eating genus Ophiophagus (King Cobras) in Asia. All these terrestrial Python species, however, were less neurotoxin-susceptible than pythons in other genera which have evolved under different predatory pressure as: the Asian species Malayopython reticulatus which is arboreal as neonates and juveniles before rapidly reaching sizes as terrestrial adults too large for sympatric Ophiophagus species to consider as prey; and the terrestrial Australian species Aspidites melanocephalus which occupies a niche, devoid of selection pressure from α-neurotoxic predatory snakes. Our findings underline the importance of positive selection in the evolution of venom resistance and suggest a complex evolutionary history involving both conserved traits and secondary evolution. This study enhances our understanding of the molecular adaptations that enable pythons to survive in environments laden with venomous threats and offers insights into the ongoing co-evolution between venomous snakes and their prey.

On the estimation of genome-average recombination rates.

The rate at which recombination events occur in a population is an indicator of its effective population size and the organism's reproduction mode. It determines the extent of linkage disequilibrium along the genome and, thereby, the efficacy of both purifying and positive selection. The population recombination rate can be inferred using models of genome evolution in populations. Classic methods based on the patterns of linkage-disequilibrium provide the most accurate estimates, providing large sample sizes are used and the demography of the population is properly accounted for. Here, the capacity of approaches based on the sequentially Markov coalescent (SMC) to infer the genome-average recombination rate from as little as a single diploid genome is examined. SMC approaches provide highly accurate estimates even in the presence of changing population sizes, providing that (1) within genome heterogeneity is accounted for and (2) classic maximum-likelihood optimization algorithms are employed to fit the model. SMC-based estimates proved sensitive to gene conversion, leading to an overestimation of the recombination rate if conversion events are frequent. Conversely, methods based on the correlation of heterozygosity succeed in disentangling the rate of crossing over from that of gene conversion events, but only when the population size is constant and the recombination landscape homogeneous. These results call for a convergence of these two methods to obtain accurate and comparable estimates of recombination rates between populations.

Cross-Reactive Polyclonal Antibodies Raised Against GalNAc-Conjugated siRNA Recognize Mostly the GalNAc Moiety.

Small interfering RNA (siRNA) is gaining momentum as a therapeutic modality with six approved products. Since siRNA has the potential to elicit undesired immune responses in patients, immunogenicity assessment is required during clinical development by regulatory authorities. In this study, anti-siRNA polyclonal antibodies were generated through animal immunization. These cross-reactive polyclonal antibodies recognized mostly the N-acetylgalactosamine (GalNAc) moiety with a small fraction against sequence-independent epitopes. We demonstrate that the polyclonal antibodies can be utilized as immunogenicity assay positive controls for the same class of GalNAc-conjugated siRNAs. In addition, anti-GalNAc mAbs showed desired sensitivity and drug tolerance, supporting their use as alternative surrogate positive controls. These findings can guide positive control selection and immunogenicity assay development for GalNAc-conjugated siRNAs and other oligonucleotide therapeutics.

A positive selection at binding site 501 in the B.1 lineage might have triggered the highly infectious sub-lineages of SARS-CoV-2

The descendants of the B lineage are the most predominant variants among the SARS-CoV-2 virus due to the incorporation of new mutations augmenting the infectivity of the virus. There is a substantial increase in the transition transversion bias, nucleotide diversity and purifying selection on the spike protein in the descendants of the B lineage of the SARS-CoV-2 virus on a temporal scale. A strong bias for C-to-U substitutions is found in the genes encoding spike protein in this lineage. The positive selection has operated on the spike gene of B lineages and its sub-lineages. The B.1 lineage has undergone positive selection on site 501 of the receptor binding domain ultimately reflected in a key substitution N501Y in its three descendant lineages namely B.1.1.7, B.1.351 and P.1. The intensity of purifying selection on the multiple sites of the spike gene has increased substantially in the sub-lineages of B.1 in a timescale. The binding site 501 on the spike protein in B lineage is found to coevolve with other amino acid sites. This study sheds light on the evolutionary trajectory of the B lineage into highly infectious descendants in the recent past under the influence of positive and purifying selection exerted by natural immunity and vaccination of the host.

Conservation genomics provides insights into genetic resilience and adaptation of the endangered Chinese hazelnut, Corylus chinensis

Global climate change has increased concerns regarding biodiversity loss. However, many key conservation issues still required further research, including demographic history, deleterious mutation load, adaptive evolution, and putative introgression. Here we generated the first chromosome-level genome of the endangered Chinese hazelnut, Corylus chinensis, and compared the genomic signatures with its sympatric widespread C. kwechowensis–C. yunnanensis complex. We found large genome rearrangements across all Corylus species and identified species-specific expanded gene families that may be involved in adaptation. Population genomics revealed that both C. chinensis and the C. kwechowensis–C. yunnanensis complex had diverged into two genetic lineages, forming a consistent pattern of southwestern-northern differentiation. Population size of the narrow southwestern lineages of both species have decreased continuously since the late Miocene, whereas the widespread northern lineages have remained stable (C. chinensis) or have even recovered from population bottlenecks (C. kwechowensis–C. yunnanensis complex) during the Quaternary. Compared with C. kwechowensis–C. yunnanensis complex, C. chinensis showed significantly lower genomic diversity and higher inbreeding level. However, C. chinensis carried significantly fewer deleterious mutations than C. kwechowensis–C. yunnanensis complex, as more effective purging selection reduced the accumulation of homozygous variants. We also detected signals of positive selection and adaptive introgression in different lineages, which facilitated the accumulation of favorable variants and formation of local adaptation. Hence, both types of selection and exogenous introgression could have mitigated inbreeding and facilitated survival and persistence of C. chinensis. Overall, our study provides critical insights into lineage differentiation, local adaptation, and the potential for future recovery of endangered trees.

Whole-genome resequencing of Chinese indigenous sheep provides insight into the genetic basis underlying climate adaptation

BackgroundChinese indigenous sheep are valuable resources with unique features and characteristics. They are distributed across regions with different climates in mainland China; however, few reports have analyzed the environmental adaptability of sheep based on their genome. We examined the variants and signatures of selection involved in adaptation to extreme humidity, altitude, and temperature conditions in 173 sheep genomes from 41 phenotypically and geographically representative Chinese indigenous sheep breeds to characterize the genetic basis underlying environmental adaptation in these populations.ResultsBased on the analysis of population structure, we inferred that Chinese indigenous sheep are divided into four groups: Kazakh (KAZ), Mongolian (MON), Tibetan (TIB), and Yunnan (YUN). We also detected a set of candidate genes that are relevant to adaptation to extreme environmental conditions, such as drought-prone regions (TBXT, TG, and HOXA1), high-altitude regions (DYSF, EPAS1, JAZF1, PDGFD, and NF1) and warm-temperature regions (TSHR, ABCD4, and TEX11). Among all these candidate genes, eight ABCD4, CNTN4, DOCK10, LOC105608545, LOC121816479, SEM3A, SVIL, and TSHR overlap between extreme environmental conditions. The TSHR gene shows a strong signature for positive selection in the warm-temperature group and harbors a single nucleotide polymorphism (SNP) missense mutation located between positions 90,600,001 and 90,650,001 on chromosome 7, which leads to a change in the protein structure of TSHR and influences its stability.ConclusionsAnalysis of the signatures of selection uncovered genes that are likely related to environmental adaptation and a SNP missense mutation in the TSHR gene that affects the protein structure and stability. It also provides information on the evolution of the phylogeographic structure of Chinese indigenous sheep populations. These results provide important genetic resources for future breeding studies and new perspectives on how animals can adapt to climate change.