Omicron escapes the majority of existing SARS-CoV-2 neutralizing antibodies The SARS-CoV-2 B.1.1.529 variant (Omicron) contains 15 mutations on the receptor-binding domain (RBD). How Omicron would evade RBD neutralizing antibodies (NAbs) requires immediate investigation. Here, we used high-throughput yeast display screening1,2 to determine the RBD escaping mutation profiles for 247 human anti-RBD NAbs and showed that the NAbs could be unsupervised clustered into six epitope groups (A-F), which is highly concordant with knowledge-based structural classifications3-5. Strikingly, various single mutations of Omicron could impair NAbs of different epitope groups. Specifically, NAbs in Group A-D, whose epitope overlap with ACE2-binding motif, are largely escaped by K417N, G446S, E484A, and Q493R. Group E (S309 site)6 and F (CR3022 site)7 NAbs, which often exhibit broad sarbecovirus neutralizing activity, are less affected by Omicron, but still, a subset of NAbs are escaped by G339D, N440K, and S371L. Furthermore, Omicron pseudovirus neutralization showed that single mutation tolerating NAbs could also be escaped due to multiple synergetic mutations on their epitopes. In total, over 85% of the tested NAbs are escaped by Omicron. Regarding NAb drugs, the neutralization potency of LY-CoV016/LY-CoV555, REGN10933/REGN10987, AZD1061/AZD8895, and BRII-196 were greatly reduced by Omicron, while VIR-7831 and DXP-604 still function at reduced efficacy. Together, data suggest Omicron would cause significant humoral immune evasion, while NAbs targeting the sarbecovirus conserved region remain most effective. Our results offer instructions for developing NAb drugs and vaccines against Omicron and future variants.
Adaptive anti-tumor immunity is orchestrated by a population of CCL5-producing tissue-resident NK cells T cell-directed cancer immunotherapy often fails to generate lasting tumor control. Harnessing additional effectors of the immune response against tumors may strengthen the clinical benefit of immunotherapies. Here, we demonstrate that therapeutic targeting of the IFN{gamma}-IL-12 pathway relies on the ability of a population of tissue-resident NK (trNK) cells to orchestrate an anti-tumor microenvironment. Particularly, utilizing an engineered adenoviral platform, we show that paracrine IL-12 enhances functional DC-CD8 T cell interactions to generate adaptive anti-tumor immunity. This effect depends on the abundance of trNK cells and specifically their capacity to produce the cDC1-chemoattractant CCL5. Failure to respond to IL-12 and other IFN{gamma}-inducing therapies such as immune checkpoint blockade in tumors with low trNK cell infiltration could be overcome by intra-tumoral delivery of CCL5. Our findings reveal a novel barrier for T cell-focused therapies and offer mechanistic insights into how T cell-NK cell-DC crosstalk can be enhanced to promote anti-tumor immunity and overcome resistance.

SignificanceWe identified the lack of CCL5-producing, tissue-resident NK (trNK) cells as a barrier to T cell-focused therapies. While IL-12 induces anti-tumoral DC-T cell crosstalk in trNK cellrich tumors, resistance to IL-12 or anti-PD-1 in trNK cellpoor tumors can be overcome by the additional delivery of CCL5.
Deep learning: step forward to high-resolution in vivo shortwave infrared imaging Shortwave infrared window (SWIR: 1000-1700 nm) represents a major improvement compared to the NIR-I region (700-900 nm) in terms of temporal and spatial resolutions in depths down to 4 mm. SWIR is a fast and cheap alternative to more precise methods such as X-ray and opto-acoustic imaging. Main obstacles in SWIR imaging are the noise and scattering from tissues and skin that reduce the precision of the method. We demonstrate that the combination of SWIR in vivo imaging in the NIRIIb region (1500-1700 nm) with advanced deep learning image analysis allows to overcome these obstacles and making a large step forward to high resolution imaging: it allows to precisely segment vessels from tissues and noise, provides morphological structure of the vessels network, with learned pseudo-3D shape, their relative position, dynamic information of blood vascularization in depth in small animals and distinguish the vessels types: artieries and veins. For demonstration we use neural network IterNet that exploits structural redundancy of the blood vessels, which provides a useful analysis tool for raw SWIR images. F1F0 ATP Hydrolysis is a Determinant of Metabolic Rate, a Correlate of Lifespan, and a Weakness of Cancer Dramatic rectal temperature decrease, in mice administered with a drug that specifically inhibits F1F0 ATP hydrolysis, suggests that F1F0 ATP hydrolysis is a major determinant of metabolic rate in vivo. Across twelve investigated species, less F1F0 ATP hydrolysis correlates with greater maximal lifespan. Specific drug inhibition of F1F0 ATP hydrolysis exerts potent anti-cancer activity in vitro. Mouse lemur, a new primate model to investigate cardiac pacemaker activity. BackgroundGray mouse lemurs (Microcebus murinus) are one of the smallest and most ancestral primates known. Their size falls in between that of mice and rats, and their genetic proximity to humans makes them an emerging model to study age-related neurodegeneration. Since Microcebus murinus replicates similar senescence processes of humans, it constitutes a useful model also for studying cardiovascular dysfunctions. However, its cardiac physiology is unknown. Thus, we investigated sinoatrial node (SAN) pacemaker activity of mouse lemurs, presenting the first characterization of heart automaticity in nonhuman primates.

Methods and ResultsWe recorded the heart automaticity in freely moving mouse lemurs and in their SAN tissues and pacemaker myocytes. Mouse lemurs have a heart rate (HR) in between those of mice and rats and a similar generation of electrical impulses in the SAN tissue. Conversely, their ventricular depolarization is more similar to large mammals. Mouse lemur SAN myocytes express the main pacemaker currents at densities similar to mice: the hyperpolarization-activated current (If) and the L-type (ICa,L) and T-type (ICa,T) calcium currents. Moreover, mouse lemur SAN myocytes display spontaneous ryanodine receptor (RyR)-dependent Ca2+ transients. Using muscle-derived stem cells (MDSCs) from the Microcebus murinus skeletal muscle, we also differentiated pacemaker-like (PML) cells showing spontaneous automaticity and expressing markers of native SAN myocytes (HCN4 and connexin-45).

ConclusionsIn conclusion, our characterization of heart automaticity in Microcebus murinus provides new opportunities for comparative cardio-physiology studies in primates and with humans, and for pharmacologic testing of molecules that could assist in modulating age-related dysfunctions of automaticity. 'rbioacc': an R-package to analyse toxicokinetic data O_LI rbioacc is an  R-package dedicated to the analysis of experimental data collected from bioaccumulation tests during which organisms are exposed to a chemical (exposure phase) and then put into a clean media (depuration phase). Internal concentrations are regularly measured over time all along the experiment.
C_LIO_LI rbioacc provides ready-to-use functions to visualize and fully analyze such data. Under a Bayesian framework, this package fits a generic one-compartment toxicokinetic (TK) model automatically built from the data. It provides TK parameter estimates (appropriate uptake and elimination rates) and bioaccumulation metrics (e.g., BCF, BSAF, BMF). All parameter estimates, bioaccumulation metrics as well as predictions of internal concentrations into organisms are delivered with their uncertainty.
C_LIO_LIBioaccumulation metrics are primarily provided in support of environmental risk assessment, in full compliance with regulatory requirements required to approve marketing applications of chemical substances.
C_LIO_LIThis paper gives brief worked examples of the use of  rbioacc from data collected through standard bioaccumulation tests, and publicly available within the scientific literature. These examples constitute step-by-step user-guides to analyze any new data set, uploaded in the right format.
C_LI Precise exogenous insertion and sequence replacements in poplar by simultaneous HDR overexpression and NHEJ suppression using CRISPR-Cas9 CRISPR-mediated genome editing has become a powerful tool for genetic modification of biological traits. However, developing an efficient, site-specific, gene knock-in system based on homology-directed DNA repair (HDR) remains a significant challenge in plants, especially in woody species like poplar. Here, we show that simultaneous inhibition of non-homologous end joining (NHEJ) recombination cofactor XRCC4 and overexpression of HDR enhancer factors CtIP and MRE11 can improve the HDR efficiency for gene knock-in. Using this approach, the BleoR gene was integrated onto the 3' end of the MKK2 MAP Kinase gene to generate a BleoR-MKK2 fusion protein. Based on exogenous BleoR expression, the HDR-mediated knock-in efficiency was up to [~]40-fold greater when using a XRCC4 silencing incorporated with a combination of CtIP and MRE11 overexpression compared to no HDR enhancement or NHEJ silencing. Furthermore, this corporation of HDR enhancer overexpression and NHEJ repression also resulted in 7-fold fewer CRISPR-induced Insertions and Deletions (InDels), resulting in no functional effects on MKK2-based salt stress responses in poplar. Therefore, this approach may be useful not only in poplar and plants or crops but also in mammalians for improving CRISPR-mediated gene knock-in efficiency. Phylogenomic discordance suggests polytomies along the backbone of the large genus Solanum Premise of the studyEvolutionary studies require solid phylogenetic frameworks, but increased volumes of phylogenomic data have revealed incongruent topologies among gene trees in many organisms both between and within genomes. Some of these incongruences indicate polytomies that may remain impossible to resolve. Here we investigate the degree of gene-tree discordance in Solanum, one of the largest flowering plant genera that includes the cultivated potato, tomato, and eggplant, as well as 24 minor crop plants.

MethodsA densely sampled species-level phylogeny of Solanum is built using unpublished and publicly available Sanger sequences comprising 60% of all accepted species (742 spp.) and nine regions (ITS, waxy, and seven plastid markers). The robustness of this topology is tested by examining a full plastome dataset with 140 species and a nuclear target-capture dataset with 39 species of Solanum (Angiosperms353 probe set).

Key resultsWhile the taxonomic framework of Solanum remained stable, gene tree conflicts and discordance between phylogenetic trees generated from the target-capture and plastome datasets were observed. The latter correspond to regions with short internodal branches, and network analysis and polytomy tests suggest the backbone is composed of three polytomies found at different evolutionary depths. The strongest area of discordance, near the crown node of Solanum, could potentially represent a hard polytomy.

ConclusionsWe argue that incomplete lineage sorting due to rapid diversification is the most likely cause for these polytomies, and that embracing the uncertainty that underlies them is crucial to understand the evolution of large and rapidly radiating lineages. A noble TGF beta biogenesis inhibitor exhibits both potent anti-fibrotic and anti-inflammatory capabilities AbstractsIdiopathy pulmonary fibrosis (IPF) is an intractable and fatal human disorder. Our previous study showed that eupatilin exerted a potent anti-fibrotic effect on both in vitro fibrogenesis and bleomycin-induced lung fibrosis model (BLM). Subsequently, an analog called ONG41008 had been identified as a more potent anti-fibrotic than eupatilin and also showed a potent anti-inflammatory capability. Orally administered ONG41008 significantly improved onset of BLM in both prophylactic and therapeutic model and its therapeutic efficacy was similarly compared to or better than pirfenidone by measuring production of collagen and hydroxyproline. Staining collagen or SMA corroborated these results.

As in vitro fibrogenesis models, DHLF (Diseased Human Lung Fibroblasts from IPF patients) and HSC (hepatic stellate cells) were used for direct effects of ONG41008 on pivotal cellular and molecular functions associated with pathogenic myofibroblasts; ONG41008 dismantled latent TGFb complex (LTC), generating inactive forms of TGF{beta}, likely limiting TGF{beta} to TGF{beta} receptor via depolymerization of F-actin and this blunted SMAD2/SMAD3 phosphorylation, thereby reprogramming EMT. A set of cell imaging studies and transcriptomic analysis were conducted to explore how ONG41008 elicited both anti-fibrotic and anti-inflammatory capabilities. Elastin (ELN) seemed to be a pioneering pharmacodynamic marker. It was also found that NOX4 played an important role in anti- fibrosis because it was functionally connected to major central nod proteins such as lysyl- oxidase (LOX) and numerous collagen family members in an ONG41008-specific fibrogenic interactome. Human NOX4 was significantly induced by TGF{beta} and completely knocked down by ONG41008. It has been shown that production of reactive oxygen species (ROS) led to activation of inflammasome. ONG41008 may be likely related to anti-inflammation, leading to a key protective effect on fibrogenesis. Concomitant with downregulation of NOX4, expression of macrophages homing chemokines, CCL2 and CCL7 were significantly attenuated by ONG41008. In vitro anti-inflammatory activities of ONG41008 were investigated in RAW264.7 cells, a mouse monocytic cell line stimulated with LPS. ONG41008 substantially attenuated TNF, CXCL10, CCL2 and CCL7, which are proinflammatory cytokine and important chemokines influencing T cells or macrophages. TNF was situated at the central nod in LPS-treated macrophages via an ONG41008-specific interactome analysis.

Taken together, ONG41008 is a TGF{beta} biogenesis inhibitor, being a potent drug for a broad range of fibrotic diseases and could antagonize inflammatory diseases as well. Flexible reprogramming of Pristionchus pacificus motivation for attacking Caenorhabditis elegans in predator-prey competition Animals with diverse diets must adapt their food priorities to a wide variety of environmental conditions. This diet optimization problem is especially complex for predators that compete with prey for food. Although predator-prey competition is widespread and ecologically critical, it remains difficult to disentangle predatory and competitive motivations for attacking competing prey. Here, we dissect the foraging decisions of the omnivorous nematode Pristionchus pacificus to reveal that its seemingly failed predatory attempts against Caenorhabditis elegans are actually motivated acts of efficacious territorial aggression. While P. pacificus easily kills and eats larval C. elegans with a single bite, adult C. elegans typically survives and escapes from bites. However, nonfatal biting can provide competitive benefits by reducing access of adult C. elegans and its progeny to bacterial food that P. pacificus also eats. We show that P. pacificus considers costs and benefits of both predatory and territorial outcomes to decide which food goal, prey or bacteria, should guide its motivation for biting. These predatory and territorial motivations impose different sets of rules for adjusting willingness to bite in response to changes in bacterial abundance. In addition to biting, predatory and territorial motivations also influence which search tactic P. pacificus uses to increase encounters with C. elegans. When treated with an octopamine receptor antagonist, P. pacificus switches from territorial to predatory motivation for both biting and search. Overall, we demonstrate that P. pacificus assesses alternate outcomes of attacking C. elegans and flexibly reprograms its foraging strategy to prioritize either prey or bacterial food. An Arginine-Rich Motif in the ORF2 Capsid Protein Regulates the Hepatitis E Virus Lifecycle and Interactions with the Host Cell Hepatitis E virus (HEV) infection is the most common cause of acute viral hepatitis worldwide. Hepatitis E is usually asymptomatic and self-limiting but it can become chronic in immunocompromised patients and is associated with increased fulminant hepatic failure and mortality rates in pregnant women. HEV genome encodes three proteins including the ORF2 protein that is the viral capsid protein. Interestingly, HEV produces 3 isoforms of the ORF2 capsid protein which are partitioned in different subcellular compartments and perform distinct functions in the HEV lifecycle. Notably, the infectious ORF2 (ORF2i) protein is the structural component of virions, whereas the genome-free secreted and glycosylated ORF2 proteins likely act as a humoral immune decoy. Here, by using a series of ORF2 capsid protein mutants expressed in the infectious genotype 3 p6 HEV strain as well as chimeras between ORF2 and the CD4 glycoprotein, we demonstrated how an Arginine-Rich Motif (ARM) located in the ORF2 N-terminal region controls the fate and functions of ORF2 isoforms. We showed that the ARM controls ORF2 nuclear translocation, promoting regulation of host antiviral responses. This motif also regulates the dual topology and functionality of ORF2 signal peptide, leading to the production of either cytosolic infectious ORF2i or reticular non-infectious glycosylated ORF2 forms. It serves as maturation site of glycosylated ORF2 by cellular proteases, and promotes ORF2-host cell membrane interactions. The identification of ORF2 ARM as a unique central regulator of the HEV lifecycle uncovers how viruses settle strategies to condense their genetic information and hijack cellular processes.

Author summaryHepatitis E virus (HEV) is the major cause of acute viral hepatitis worldwide. Although infection with HEV is usually self-resolving, it can cause up to 30% mortality in pregnant women in the third trimester. There is no specific treatment nor universal vaccine to fight against HEV. In our study, we focused on the HEV ORF2 capsid protein which is produced as different forms that perform distinct functions in the HEV lifecycle. The infectious ORF2i form is the structural component of virions, while the other forms likely act as an immune decoy. Herein, we deciphered the molecular determinants of ORF2 multifunctionality. We identified an Arginine-Rich Motif (ARM) located in the ORF2 N-terminal region that controls the subcellular localization, the fate and functions of ORF2 forms. It also promotes ORF2-host cell interactions. Our observations highlight the ORF2 ARM as a unique central regulator of ORF2 addressing that finely controls the HEV lifecycle. Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed. Designing viral diagnostics with model-based optimization Diagnostics, particularly for rapidly evolving viruses, stand to benefit from a principled, measurement-driven design that harnesses machine learning and vast genomic data--yet the capability for such design has not been previously built. Here, we develop and extensively validate an approach to designing viral diagnostics that applies a learned model within a combinatorial optimization framework. Concentrating on CRISPR-based diagnostics, we screen a library of 19,209 diagnostic-target pairs and train a deep neural network that predicts, from RNA sequence alone, diagnostic signal better than contemporary techniques. Our model then makes it possible to design assays that are maximally sensitive over the spectrum of a viruss genomic variation. We introduce ADAPT (https://adapt.guide), a system for fully-automated design, and use ADAPT to design optimal diagnostics for the 1,933 vertebrate-infecting viral species within 2 hours for most species and 24 hours for all but 3. We experimentally show ADAPTs designs are sensitive and specific down to the lineage level, including against viruses that pose challenges involving genomic variation and specificity. ADAPTs designs exhibit significantly higher fluorescence and permit lower limits of detection, across a viruss entire variation, than the outputs of standard design techniques. Our model-based optimization strategy has applications broadly to viral nucleic acid diagnostics and other sequence-based technologies, and, paired with clinical validation, could enable a critically-needed, proactive resource of assays for surveilling and responding to pathogens. Royal knifefish generate powerful suction feeding through large neurocranial elevation and high epaxial muscle power Suction feeding in ray-finned fishes involves powerful buccal cavity expansion to accelerate water and food into the mouth. Previous XROMM studies in largemouth bass (Micropterus salmoides), bluegill sunfish (Lepomis macrochirus), and channel catfish (Ictalurus punctatus) have shown that more than 90% of suction power in high performance strikes comes from the axial musculature. Thus, the shape of the axial muscles and skeleton may impact suction feeding mechanics. Royal knifefish (Chitala blanci) have an unusual postcranial morphology, with a ventrally flexed vertebral column and relatively large mass of epaxial muscle. Based on their body shape, we hypothesized that royal knifefish would generate high power strikes by utilizing large neurocranial elevation, vertebral column extension, and epaxial shortening. As predicted, C. blanci generated high suction expansion power compared to the other three species studied to date (up to 160 W), which was achieved by increasing both the rate of volume change and the intraoral subambient pressure. The large epaxial muscle (25% of body mass) shortened at high velocities to produce large neurocranial elevation and vertebral extension (up to 41 deg, combined), as well as high muscle mass-specific power (up to 800 W kg-1). For the highest power strikes, axial muscles generated 95% of the power, and 64% of the axial muscle mass consisted of the epaxial muscles. The epaxial-dominated suction expansion of royal knifefish supports our hypothesis that postcranial morphology may be a strong predictor of suction feeding biomechanics.

SUMMARY STATEMENTRoyal knifefish rely on their distinct postcranial morphology--with a curved vertebral column and large dorsal body muscles--to produce large neurocranial elevation and powerful suction feeding. Topaz1, an essential gene for murine spermatogenesis, down-regulates the expression of numerous testis-specific long non-coding RNAs. Spermatogenesis involves coordinated processes, including meiosis, to produce functional gametes. We previously reported Topaz1 as a germ cell-specific gene highly conserved in vertebrates. Topaz1 knockout males are sterile with testes that lack haploid germ cells because of meiotic arrest after prophase I. To better characterize Topaz1-/- testes, we used RNA-sequencing analyses at two different developmental stages (P16 and P18). The absence of TOPAZ1 disturbed the expression of genes involved in microtubule and/or cilium mobility, which was consistent with testicular histology showing the disruption of microtubules and centrosomes. Moreover, a quarter of P18 dysregulated genes are long non-coding RNAs (lncRNAs), and three of them are testis-specific and located in spermatocytes, their expression starting between P11 and P15. The suppression of one of them, 4939463O16Rik, did not alter fertility although sperm parameters were disturbed and sperm concentration fell. The transcriptome of P18-4939463O16Rik-/- testes was altered and the molecular pathways affected included microtubule-based processes, the regulation of cilium movement and spermatogenesis. The absence of TOPAZ1 protein or 4930463O16Rik produced the same enrichment clusters in mutant testes despite a contrasted phenotype on male fertility. In conclusion, TOPAZ1 appeared to stabilize the expression of numerous lncRNAs. Its suppression is not essential for fertility but required during the terminal differentiation of male gametes.

Author SummaryThe Topaz1 gene was initially characterized during the initiation of meiosis in the sheep fetal ovary. In order to determine its function, a KO of the murine gene was performed. In this species, only males were sterile and spermatogenesis was blocked before the first meiotic division. Here, we show that cytoskeletal elements are markedly disturbed in mutant testes, indicating that these elements play an important function in spermatogenesis. While the mitotic spindle of spermatogonia was normal, the meiotic spindle of spermatocytes was hemi-spindle-shaped and the homologous chromosome pairs could position themselves on the equatorial plate. In addition, lncRNAs account for 25% of genes whose expression in testes varies significantly in the absence of Topaz1. This suggests a key role for these factors in spermatogenesis. Largely testis-specific, they may be involved in spermatogenesis and play a more or less critical role in mouse fertility, which probably also depends on their redundancies. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. P2X2 receptor subunit interfaces are missense variant hotspots where mutations tend to increase apparent ATP affinity Background and PurposeP2X receptors (P2XRs) are trimeric ligand-gated ion channels (LGICs) that open a cation-selective pore in response to ATP binding to their large extracellular domain (ECD). The seven known P2XR subtypes can assemble as homo- or heterotrimeric complexes and contribute to numerous physiological functions, including nociception, inflammation and hearing. The overall structure of P2XRs is well established, but little is known about the spectrum and prevalence of human genetic variations and the functional implications in specific domains.

Experimental ApproachHere we examine the impact of P2X2 receptor (P2X2R) inter-subunit interface missense variants identified in the human population or through structural predictions. We test both single and double mutants through electrophysiological and biochemical approaches.

Key resultsWe demonstrate that predicted ECD inter-subunit interfaces display a higher-than-expected density of missense variations and that the majority of mutations that disrupt putative inter-subunit interactions result in channels with higher apparent ATP affinity. Lastly, we show that double mutants at the subunit interface show significant energetic coupling, especially if located in close proximity.

Conclusions and ImplicationsWe provide the first structural mapping of the mutational burden across the human population in a LGIC and show that the density of missense mutations is constrained between protein domains, indicating evolutionary selection at the domain level. Our data may indicate that, unlike other LGICs, P2X2Rs have evolved an intrinsically high threshold for activation, possibly to allow for additional modulation or as a cellular protection mechanism against overstimulation.

Bullet point summary What is already known: O_LIP2X2 receptors are ATP-activated ion channels implicated in hearing and nociceptice pathways
C_LI
 What this study adds: O_LIA structural mapping of missense variants observed in the human population
C_LIO_LIWe identify the intersubunit-interface as a variant hotspot and decipher functional impact of mutations
C_LI
 Clinical significance: O_LIThe development of both inhibitors and activators of P2X2 receptor function may be required
C_LI Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Advanced paternal age diversifies individual trajectories of vocalization patterns in neonatal mice Infant crying is an innate communicative behavior that is frequently impaired in certain neurodevelopmental disorders (NDDs). Since advanced paternal age is a reported risk factor for NDDs in offspring, we evaluated the impact of a fathers age on early vocal development in C57BL/6J mice. We recorded and applied a unique combination of computational analyses to ultrasonic vocalizations (USVs) emitted by mouse pups sired by young and aged fathers. Our data showed that advanced paternal age reduced the number and duration of USVs, and altered the syllable composition in pups. Moreover, pups born to young fathers showed convergent vocal characteristics with a rich repertoire during postnatal development, while those born to aged fathers exhibited more divergent vocal patterns with limited repertoire. Principal component analysis in conjunction with clustering analysis demonstrated that pups from aged fathers deviated from typical trajectories of vocal development, which were considered as atypical individuals. Thus, our study indicates that advanced paternal age has a significant effect on offsprings early vocal development. It is suggested that the trajectories of vocal development could be a useful marker of the NDD-like phenotype associated with the advanced paternal age. In addition, our comprehensive computational analysis described here is an effective approach to characterize the altered individual diversity relevant to neurodevelopmental disorders.

One Sentence SummaryAdvanced paternal age affects vocal development in early postnatal mice, with more pups showing atypical developmental trajectories. Creatine transporter deficient rat model show motor dysfunction, cerebellar alterations and muscle creatine deficiency without muscle atrophy Creatine (Cr) is a nitrogenous organic acid and plays roles as fast phosphate energy buffer to replenish ATP, osmolyte, antioxidant, neuromodulator, and as a compound with anabolic and ergogenic properties in muscle. Cr is taken from the diet or endogenously synthetized by the enzymes AGAT and GAMT, and specifically taken up by the transporter SLC6A8. Loss-of-function mutations in the genes encoding for the enzymes or the transporter cause Cerebral Creatine Deficiency Syndromes (CCDS). CCDS are characterized by brain Cr deficiency, intellectual disability with severe speech delay, behavioral troubles, epilepsy and motor dysfunction. Among CCDS, the X-linked Cr transporter deficiency (CTD) is the most prevalent with no efficient treatment so far. Different animal models of CTD show reduced brain Cr levels, cognitive deficiencies and together they cover other traits similar to those of patients. However, motor function was poorly explored in CTD models and some controversies in the phenotype exist in comparison with CTD patients. Our recently described Slc6a8Y389C knock-in (KI) rat model of CTD showed mild impaired motor function, morphological alterations in cerebellum, reduced muscular mass, Cr deficiency and increased guanidinoacetate content in muscle, although no consistent signs of muscle atrophy. Our results indicate that such motor dysfunction co-occurred with both nervous and muscle dysfunction, suggesting that muscle strength and performance as well as neuronal connectivity might be affected by this Cr deficiency in muscle and brain. A combined human in silico and CRISPR/Cas9-mediated in vivo zebrafish based approach for supporting gene target validation in early drug discovery The clinical heterogeneity of heart failure has challenged our understanding of the underlying genetic mechanisms of this disease. In this respect, large-scale patient DNA sequencing studies have become an invaluable strategy for identifying potential genetic contributing factors. The complex aetiology of heart failure, however, also means that in vivo models are vital to understand the links between genetic perturbations and functional impacts. Traditional approaches (e.g. genetically-modified mice) are optimal for assessing small numbers of proposed target genes, but less practical when multiple targets are identified. The zebrafish, in contrast, offers great potential for higher throughput in vivo gene functional assessment to aid target prioritisation and support definitive studies undertaken in mice. Here we used whole-exome sequencing and bioinformatics on human patient data to identify 3 genes (API5, HSPB7, and LMO2) suggestively associated with heart failure that were also predicted to play a broader role in disease aetiology. The role of these genes in cardiovascular system development and function was then further investigated using in vivo CRISPR/Cas9-mediated gene mutation analysis in zebrafish. We observed multiple impacts in F0 knockout zebrafish embryos (crispants) following effective somatic mutation, including reductions in ventricle size, pericardial oedema, and chamber malformation. In the case of lmo2, there was also a significant impact on cardiovascular function as well as an expected reduction in erythropoiesis. The data generated from both the human in silico and zebrafish in vivo assessments undertaken supports roles for API5, HSPB7, and LMO2 in human cardiovascular disease and identifies them as potential drug targets for further investigation. The data presented also supports the use of human in silico genetic variant analysis, in combination with zebrafish crispant phenotyping, as a powerful approach for assessing gene function as part of an integrated multi-level drug target validation strategy. Parent-offspring inference in inbred populations Genealogical relationships are fundamental components of genetic studies. However, it is often challenging to infer correct and complete pedigrees even when genome-wide information is available. For example, inbreeding can obfuscate genetic differences between individuals, making it difficult to even distinguish first-degree relatives such as parent-offspring from full siblings. Similarly, genotyping errors can interfere with the detection of genetic similarity between parents and their offspring. Inbreeding is common in natural, domesticated, and experimental populations and genotyping of these populations often has more errors than in human datasets, so efficient methods for building pedigrees under these conditions are necessary. Here, we present a new method for parent-offspring inference in inbred pedigrees called SPORE (Specific Parent-Offspring Relationship Estimation). SPORE is vastly superior to existing pedigree-inference methods at detecting parent-offspring relationships, in particular when inbreeding is high or in the presence of genotyping errors, or both. SPORE therefore fills an important void in the arsenal of pedigree inference tools.

Author SummaryKnowing the genealogical relationships among individuals is critical for genetic analyses, such as for identifying the mutations that cause diseases or that contribute to valuable agricultural traits such as milk production. Although many tools exist for establishing pedigrees using genetic information, these tools fail when individuals are highly inbred, such as in domesticated animals, or in groups of people in which consanguineous matings are common. Furthermore, existing tools do not work well when genetic information has errors at levels observed in modern datasets. Here, we present a novel approach to solve these problems. Our method is significantly more accurate than existing tools and more tolerant of errors in the genetic data. We expect that our method, which is simple to use and computationally efficient, will be a useful tool in a diversity of settings, from the studies of human and natural populations, to agricultural and experimental settings. Human 14-3-3 proteins site-selectively bind the mutational hotspot region of SARS-CoV-2 nucleoprotein modulating its phosphoregulation Phosphorylated within its Ser/Arg-rich region, the SARS-CoV-2 nucleoprotein (N) recruits the phosphopeptide-binding human 14-3-3 proteins that play a well-recognized role in replication of many viruses. Here we use genetic code expansion to demonstrate that phosphorylation of SARS-CoV-2 N at either of two pseudo-repeats centered at Ser197 and Thr205 is sufficient for 14-3-3 binding. According to fluorescence anisotropy, the pT205-motif, present in SARS-CoV-2 but not in SARS-CoV, is preferred over the pS197-motif by all seven human 14-3-3 isoforms, which display unforeseen pT205/pS197 binding selectivity hierarchy. Crystal structures demonstrate that pS197 and pT205 are mutually exclusive 14-3-3-binding sites, whereas SAXS and biochemical data indicate 14-3-3 binding occludes the Ser/Arg-rich region, inhibiting its dephosphorylation. This Ser/Arg-rich region of N is highly prone to mutations, as exemplified by the Omicron and Delta variants, with our data suggesting how the strength of its 14-3-3 binding can be linked with the replicative fitness of the virus. Engineered bacteria detect tumor DNA Advances in bacterial engineering have catalysed the development of living cell diagnostics and therapeutics1-3, including microbes that respond to diseases such as gut inflammation4, intestinal bleeding5, pathogens6 and hypoxic tumors7. Bacteria can easily access the entire gastrointestinal tract via oral administration8, and they can produce outputs that can be noninvasively measured in stool4 or urine7. Cellular memory, such as bistable switches4,9,10 or genomic rearrangement11, has been used to allow bacteria to store information over time. However, living biosensors have not yet been engineered to detect specific DNA sequences or mutations from outside the cell. Here, we engineer naturally competent Acinetobacter baylyi to detect donor DNA from the genomes of colorectal cancer (CRC) cells and organoids. We characterize the functionality of the biosensors in vitro with co-culture assays and then validate in vivo with sensor bacteria delivered orally or rectally into mice injected with orthotopic donor CRC organoids. We observe horizontal gene transfer from the tumor to the sensor bacteria in vivo, allowing their detection in stool. The sensor bacteria achieved 100% discrimination between mice with and without CRC using both delivery methods. Our findings establish a framework for biosensing applications that require the detection of mutations or organisms within environments that are difficult to sample. In addition, the platform can be readily expanded to include in situ production and delivery of therapeutic payloads at the detection site. ksrates: positioning whole-genome duplications relative to speciation events in KS distributions SummaryWe present ksrates, a user-friendly command-line tool to position ancient whole-genome duplication (WGD) events with respect to speciation events in a phylogeny by comparing paralog and ortholog KS distributions derived from genomic or transcriptomic sequences, while adjusting for substitution rate differences among the lineages involved.

Availability and implementationksrates is implemented in Python 3 and as a Nextflow pipeline. The source code, Singularity and Docker containers, documentation and tutorial are available via https://github.com/VIB-PSB/ksrates.

Contactsteven.maere@ugent.vib.be, rolf.lohaus@ugent.vib.be