Enhancement of CRISPR-mediated homologous recombination-based gene knock-in efficiency by simultaneous HDR overexpression and NHEJ suppression in poplar 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.
High expression of VRT2 increases the number of rudimentary basal spikelets in wheat Spikelets are the fundamental building blocks of Poaceae inflorescences and their development and branching patterns determine the various inflorescence architectures and grain yield of grasses. In wheat, the central spikelets produce the most and largest grains, while spikelet size gradually decreases acro- and basipetally, giving rise to the characteristic lanceolate shape of wheat spikes. The acropetal gradient correlates with the developmental age of spikelets, however the basal spikelets are developed first and the cause of their small size and rudimentary development is unclear. Here, we adapted G&T-seq, a low-input transcriptomics approach, to characterise gene expression profiles within spatial sections of individual spikes before and after the establishment of the lanceolate shape. We observed larger differences in gene expression profiles between the apical, central and basal sections of a single spike than between any section belonging to consecutive developmental timepoints. We found that SVP MADS-box transcription factors, including VRT-A2, are expressed highest in the basal section of the wheat spike and display the opposite expression gradient to flowering E-class SEP1 genes. Based on multi-year field trials and transgenic lines, we show that higher expression of VRT-A2 in the basal sections of the spike is associated with increased numbers of rudimentary basal spikelets. Our results, supported by computational modelling, suggest that the delayed transition of basal spikelets from vegetative to floral developmental programmes results in the lanceolate shape of wheat spikes. This study highlights the value of spatially resolved transcriptomics to gain new insights into developmental genetics pathways of grass inflorescences.

One sentence summaryLarge transcriptional gradients exist within a wheat spike and are associated with rudimentary basal spikelet development, resulting in the characteristic lanceolate shape of wheat spikes.
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. Genome editing in the unicellular holozoan Capsaspora owczarzaki suggests a premetazoan role for the Hippo pathway in multicellular morphogenesis Animal development is mediated by a surprisingly small set of canonical signaling pathways such as Wnt, Hedgehog, TGF-beta, Notch and Hippo pathways. Although once thought to be present only in animals, recent genome sequencing has revealed components of these pathways in the closest unicellular relatives of animals. These findings raise questions about the ancestral functions of these developmental pathways and their potential role in the emergence of animal multicellularity. Here, we provide the first functional characterization of any of these developmental pathways in unicellular organisms by developing techniques for genetic manipulation in Capsaspora owczarzaki, a close unicellular relative of animals that displays aggregative multicellularity. We then use these tools to characterize the Capsaspora ortholog of the Hippo signaling nuclear effector YAP/TAZ/Yorkie (coYki), a key regulator of tissue size in animals. In contrast to what might be expected based on studies in animals, we show that coYki is dispensable for cell proliferation but regulates cytoskeletal dynamics and the three-dimensional (3D) shape of multicellular structures. We further demonstrate that the cytoskeletal abnormalities of individual coYki mutant cells underlie the abnormal 3D shape of coYki mutant aggregates. Taken together, these findings implicate an ancestral role for the Hippo pathway in cytoskeletal dynamics and multicellular morphogenesis predating the origin of animal multicellularity, which was co-opted during evolution to regulate cell proliferation. Variation in female songbird state determines signal strength needed to evoke copulation It is the female response to male signals that determines courtship success. In most songbirds, females control reproduction via the copulation solicitation display (CSD), an innate, stereotyped posture produced in direct response to male displays. Because CSD can be elicited in the absence of males by the presentation of recorded song, CSD production enables investigations into the effects of underlying signal features and behavioral state on female mating preferences. Using computer vision to quantify CSD trajectory in female brown-headed cowbirds (Molothrus ater), we show that both song quality and a females internal state predict CSD production, as well as the onset latency and duration of the display. We also show that CSD can be produced in a graded fashion based on both signal strength and internal state. These results emphasize the importance of underlying receiver state in determining behavioral responses and suggest that female responsiveness acts in conjunction with male signal strength to determine the efficacy of male courtship. An Evolutionary Process Without Variation and Selection Natural selection successfully explains how organisms accumulate adaptive change despite that traits acquired over a lifetime are eliminated at the end of each generation. However, in some domains that exhibit cumulative, adaptive change--e.g., cultural evolution, and earliest life--acquired traits are retained; these domains do not face the problem that Darwins theory was designed to solve. Lack of transmission of acquired traits occurs when germ cells are protected from environmental change, due to a self-assembly code used in two distinct ways: (i) actively interpreted during development to generate a soma, and (ii) passively copied without interpretation during reproduction to generate germ cells. Early life and cultural evolution appear not to involve a self-assembly code used in these two ways. We suggest that cumulative, adaptive change in these domains is due to a lower-fidelity evolutionary process, and model it using Reflexively Autocatalytic and Foodset-generated networks. We refer to this more primitive evolutionary process as Self-Other Reorganisation (SOR) because it involves internal self-organising and self-maintaining processes within entities, as well as interaction between entities. SOR encompasses learning but in general operates across groups. We discuss the relationship between SOR and Lamarckism, and illustrate a special case of SOR without variation. JACKIE: Fast enumeration of genomic single- and multi-copy target sites and their off-targets for CRISPR and other engineered nuclease systems ZFP-, TALE-, and CRISPR-based methods for genome, epigenome editing and imaging have provided powerful tools to interrogate functions of genomes. Targeting sequence design is vital to the success of these experiments. While existing design software mainly focus on designing target sequence for specific elements, we report here the implementation of JACKIE (Jackie and Alberts Comprehensive K-mer Instances Enumerator), a suite of software for enumerating all single- and multi-copy sites in the genome that can be incorporated for genome-scale designs as well as loaded onto genome browsers alongside other tracks for convenient web-based graphic-user-interface (GUI)-enabled design. We also implement fast algorithms to identify sequence neighborhoods or off-target counts of targeting sequences so that designs with low probability of off-target can be identified among millions of design sequences in reasonable time. We demonstrate the application of JACKIE-designed CRISPR site clusters for genome imaging. Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes IntroductionIncreasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero.

MethodsWe propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F -penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation.

ResultsWe evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F -penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls.

Discussion/ConclusionWe have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Human foraging for primary rewards is guided by fronto-hippocampal dynamics of future anticipation Foraging is a fundamental food-seeking behavior in a wide range of species that enables survival in an uncertain world. During foraging, behavioral agents constantly face a trade-off between staying in their current location or exploring another. Despite ethological generality and importance of foraging, it remains unclear how the human brain guides continuous decision in such situations. Here we show that anticipatory activity dynamics in the anterior prefrontal cortex (aPFC) and hippocampus underpin foraging for primary rewards. While functional MRI was performed, humans foraged for real liquid rewards available after tens of seconds, and continuous decision during foraging was tracked by a dynamic pattern of brain activity that reflected anticipation of a future reward. When the dynamic anticipatory activity in the aPFC was enhanced, humans remained in their current environment, but when this activity diminished, they explored a new environment. Moreover, the anticipatory activity in the aPFC and hippocampus was associated with distinct decision strategies: aPFC activity was enhanced in humans adopting an exploratory strategy, whereas those remaining stationary showed enhanced activity in the hippocampus. Our results suggest that anticipatory dynamics in the fronto-hippocampal mechanisms underlie continuous decision-making during human foraging. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Illuminating the plant transposon insertion landscape in real time using Cas9-targeted Nanopore sequencing and a novel pipeline Transposable elements (TEs) are a major source of genomic novelty, contributing to plant adaptation, speciation and new cultivar production. The often large, complex genomes of plants make identifying TE insertions from short reads challenging, while whole-genome sequencing remains expensive. In addition, rare somatic insertions reflecting mobilome dynamic are difficult to track using short reads. To overcome these challenges, we applied the Cas9-targeted Nanopore sequencing (CANS) approach and a novel pipeline, NanoCasTE, to trace both genetically inherited and somatic TE insertions in plants. We performed CANS of EVD retrotransposon in wild-type Arabidopsis thaliana and obtained up to 40x coverage of the targets after only a few hours of sequencing on a MinION sequencer. To estimate the ability to detect new TE and T-DNA insertions, we exploited the A. thaliana ddm1 mutant, which has elevated TE activity. Using CANS, we detect the hemizygous T-DNA insertion site and genetically inherited EVD insertions. CANS also detected >800 somatic TE insertions allowing monitoring of TE insertions in real time. By this, we built a high-density physical map of EVD insertion preferences in the A. thaliana ddm1 and found that EVD insertions are biased toward pericentromeric regions in ddm1 that differ from wild-type A. thaliana ecotypes suggesting a pivotal role of DNA methylation in the shaping of EVD insertion preferences. Using CANS, we were able to detect transposition events of ONSEN transposon family that occurred in response to heat stress and found hotspots of ONSEN insertions. Our results indicate that CANS and NanoCasTE are effective tools for the detection of TEIs in real time, opening an avenue for obtaining a deeper understanding of mobilome organization in response to stress and in different genetic backgrounds. The ability of CANS to detect genetically inherited TEIs and T-DNA integration sites make it useful for screening insertion mutants and transgenic plants. Fatal neuroinvasion and SARS-CoV-2 tropism in K18-hACE2 mice is partially independent on hACE2 expression Animal models recapitulating distinctive features of severe COVID-19 are critical to enhance our understanding of SARS-CoV-2 pathogenesis. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. The precise mechanisms of lethality in this mouse model remain unclear. Here, we evaluated the spatiotemporal dynamics of SARS-CoV-2 infection for up to 14 days post-infection. Despite infection and moderate pneumonia, rapid clinical decline or death of mice was invariably associated with viral neuroinvasion and direct neuronal injury (including brain and spinal neurons). Neuroinvasion was observed as early as 4 dpi, with virus initially restricted to the olfactory bulb supporting axonal transport via the olfactory neuroepithelium as the earliest portal of entry. No evidence of viremia was detected suggesting neuroinvasion occurs independently of entry across the blood brain barrier. SARS-CoV-2 tropism was not restricted to ACE2-expressing cells (e.g., AT1 pneumocytes), and some ACE2-positive lineages were not associated with the presence of viral antigen (e.g., bronchiolar epithelium and brain capillaries). Detectable ACE2 expression was not observed in neurons, supporting overexpression of ACE2 in the nasal passages and neuroepithelium as more likely determinants of neuroinvasion in the K18-hACE2 model. Although our work incites caution in the utility of the K18-hACE2 model to study global aspects of SARS-CoV-2 pathogenesis, it underscores this model as a unique platform for exploring the mechanisms of SARS-CoV-2 neuropathogenesis that may have clinical relevance acknowledging the growing body of evidence that suggests COVID-19 may result in long-standing neurologic consequences.

IMPORTANCECOVID-19 is predominantly a respiratory disease caused by SARS-CoV-2 that has infected more than 191 million people with over 4 million fatalities (2021-07-20). The development of animal models recapitulating distinctive features of severe COVID-19 is critical to enhancing our understanding of SARS-CoV-2 pathogenesis and in the evaluation of vaccine and therapeutic efficacy. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. Here, we show lethality of this model is invariably associated with viral neuroinvasion linked with viral replication and assembly. Importantly, pneumonia albeit invariably present was generally moderate with the absence of culturable infectious virus at peak neuroinvasion. The dynamics of viral neuroinvasion and pneumonia were only partially dependent on hACE2. Overall, this study provides an in-depth sequential characterization of the K18-hACE2 model following SARS-CoV-2 infection, highlighting its significance to further study the mechanisms of SARS-CoV-2 neuropathogenesis. Cell Rearrangement Generates Pattern Emergence as a Function of Temporal Morphogen Exposure. As tissues elongate, cell rearrangement alters positional information in manner that must be accounted for to generate gene expression pattern. How this is achieved during paraxial mesoderm elongation is unknown. By reverse-engineering gene regulatory networks that predict single cell expression trajectories across the tissue, we find a network capable of recapitulating the full range of dynamic differentiation profiles observed both in vivo and in vitro. Simulating gene expression profiles on in toto cell tracking data sets reveal that temporal exposure to Wnt and FGF is generated by cell movement. The absence of reversal in gene expression towards a more premature gene expression state predicts the generation of aberrant tbx6 expression in the posterior progenitor zone that we then confirm by quantitative single cell imaging. Taken together, these results demonstrate cell rearrangement tunes the dynamics of mesoderm progenitor differentiation to generate pattern emergence as a function of temporal Wnt and FGF exposure. Staphylococcus aureus-serine protease-like protein B (SplB) activates PAR2 and induces endothelial barrier dysfunction Staphylococcus aureus (S. aureus) is a major cause of life-threatening systemic infection in humans. To cause blood stream infections such as sepsis and endocarditis, the bacteria must overcome the hosts endothelial barrier. The serine protease-like proteins (Spls) of S. aureus are known to contribute to pneumonia and allergic airway inflammation in animal models, but their role in endothelial damage is unknown. Here we demonstrate that SplB induces proinflammatory cytokine release in primary human vascular endothelial cells (HUVECs) in vitro. Mechanistically, we show that SplB selectively cleaves and activates human proteinase-activated receptor-2 (PAR2), and induces biased signaling via G13, Gi1/oB (NF-k{beta}), and {beta}-arrestin-1/-2. This activation did not trigger Gq/11-mediated calcium release nor ERK phosphorylation. Inhibition of PAR2 in HUVECs reduced the SplB-mediated cytokine release. Intravital microscopy of cremaster muscles in mice demonstrated that administration of SplB causes microvascular leakage. Genetic deletion of PAR2 in mice or neutralization of SplB with a monoclonal antibody preserved the endothelial barrier. This study identifies PAR2 as a receptor and substrate for SplB and highlights its role in mediating endothelial damage.



O_FIG O_LINKSMALLFIG WIDTH=199 HEIGHT=200 SRC="FIGDIR/small/424670v2_ufig1.gif" ALT="Figure 1">
View larger version (32K):
org.highwire.dtl.DTLVardef@16a60a2org.highwire.dtl.DTLVardef@ee4e41org.highwire.dtl.DTLVardef@153cd7borg.highwire.dtl.DTLVardef@d746eb_HPS_FORMAT_FIGEXP  M_FIG Graphical abstract C_FIG Robust retrieval of data stored in DNA by de Bruijn graph-based de novo strand assembly DNA data storage is a rapidly developing technology with great potential due to its high density, long-term durability, and low maintenance cost. The major technical challenges include various errors, such as the strand breaks, rearrangements, and indels that frequently arise during DNA synthesis, amplification, sequencing, and preservation. Through a de novo assembly strategy, we developed an algorithm based on the de Bruijn graph and greedy path search (DBGPS) to address these issues. DBGPS shows distinct advantages in handling DNA breaks, rearrangements, and indels. The robustness of DBGPS is demonstrated by accelerated aging, multiple independent data retrievals, deep error-prone PCR, and large data scale simulations. Remarkably, 6.8 MB of data can be retrieved accurately from a seriously corrupted sample that has been treated at 70 {degrees}C for 70 days. With DBGPS, we were able to achieve a logical density of 1.30 bits/cycle and a physical density of 295 PB/g.

One-Sentence SummaryA de Bruijn graph-based de novo assembly algorithm for DNA data storage enables fast and robust data readouts even with DNA samples that have been severely corrupted. First detection of herpesvirus and frequency of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. A CTP-dependent gating mechanism enables ParB spreading on DNA Proper chromosome segregation is essential in all living organisms. The ParA-ParB-parS system is widely employed for chromosome segregation in bacteria. Previously, we showed that Caulobacter crescentus ParB requires cytidine triphosphate to escape the nucleation site parS and spread by sliding to the neighboring DNA 1. Here, we provide the structural basis for this transition from nucleation to spreading by solving co-crystal structures of a C-terminal domain truncated C. crescentus ParB with parS and with a CTP analog. Nucleating ParB is an open clamp, in which parS is captured at the DNA-binding domain (the DNA-gate). Upon binding CTP, the N-terminal domain (NTD) self-dimerizes to close the NTD-gate of the clamp. The DNA-gate also closes, thus driving parS into a compartment between the DNA-gate and the C-terminal domain. CTP hydrolysis and/or the release of hydrolytic products are likely associated with re-opening of the gates to release DNA and to recycle ParB. Overall, we suggest a CTP-operated gating mechanism that regulates ParB nucleation, spreading, and recycling. Relationship between shear wave velocity and muscle activation is inconsistent across different muscle types There is an increasing use of shear wave ultrasound elastography to quantify mechanical properties of muscles under various conditions such as changes muscle length and levels of activation in healthy and pathological muscle. However, little is known about the variability in shear wave velocity among muscles as most studies investigate one specific muscle. The purpose of this study was to determine if the relationship between SWV and muscle activation is consistent across muscles with different architectures: biceps brachii, tibialis anterior, and medial gastrocnemius, All measures were made at matching levels of activation and approximately at the optimal length for each muscle to control for length-dependent changes in the relationship between activation and force or stiffness. We also conducted a control experiment to determine how the passive force within a muscle alters the relationship between muscle activation and shear wave velocity. The relationship between shear wave velocity-squared and activation above 10% MVC differed across muscles, with biceps brachii and medial gastrocnemius showing a lower slope than tibialis anterior. Shear wave velocity-squared also differed between muscles at the shortest length (p<0.001) and the increase in shear wave velocity-squared with muscle lengthening also differed between muscle types (p = 0.005) Muscle-specific differences could not be explained by the architectural features such as pennation angle, fiber length, and physiological cross-sectional area. Our results demonstrate that there is not a unique relationship between muscle activation and shear wave velocity, highlighting the importance of understanding the many factors contributing to shear wave propagation in muscle before elastography can be used to make quantitative comparisons across muscle types. Heterogeneity of circulating epithelial cells in breast cancer at single-cell resolution: identifying tumor and hybrid cells Circulating tumor cells and hybrid cells formed by the fusion of tumor cells with normal cells are leading players in metastasis and have prognostic relevance. Circulating tumor cells and hybrid cells are identified as CD45-negative and CD45-positive epithelial cells. However, such an approach is challenging because epithelial cells are observed in the blood of healthy individuals. In this study, we applied single-cell RNA sequencing to profile CD45-negative and CD45-positive circulating epithelial cells (CECs) in 20 breast cancer patients and one healthy donor. DNA ploidy analysis was used to identify the tumor and hybrid cells among CD45[boxh] and CD45+ CECs in patients, respectively. Functional enrichment analysis was applied to characterize aneuploid and diploid cells. Diploid cells were also annotated to generate cell-type candidates and analyzed for copy-number aberrations (CNAs) to confirm or refute their tumor origin. CD45[boxh] and CD45+ CECs were found in cancer patients (25.5 (range 0-404) and median 6.5 (0-147)) and the healthy donor (8 and 11 cells) and divided into three clusters. Two CD45[boxh] CEC clusters were predominantly aneuploid (97% and 98%), but one cluster contained more diploid (59%) than aneuploid cells. CD45+ CECs were mostly diploid: only clusters 1 and 2 had aneuploid cells (16% and 2%). Diploid CD45[boxh] and CD45+ CECs were annotated as different immune cells and surprisingly harbored many CNAs. Cancer-associated signaling pathways were found only in aneuploid cells of CD45[boxh] CEC cluster 1 and diploid cells of CD45+ CEC cluster 1. Thus, our findings suggest that CECs in breast cancer patients are a highly heterogeneous population comprising aneuploid (tumor and hybrid) and diploid (normal) cells. DNA ploidy analysis is an effective instrument for identifying tumor and hybrid cells among CD45[boxh] and CD45+ CECs, respectively. Hematopoietic differentiation is characterized by a transient peak of entropy at a single cell level. Hematopoietic differentiation has been metaphorically represented as linear trajectories with discrete steps from hematopoietic stem cells to mature cells. While the transcriptional state of cells at the beginning or at the end of these trajectories are well described from bulk analysis, what happens in the intermediate states has remained elusive until the use of single cell approaches. Applying Shannon entropy to measure cell-to-cell variability among cells at the same stage of differentiation, we observed a transient peak of gene expression variability in all the hematopoietic differentiation pathways. Strikingly, genes with the highest entropy variation in a given differentiation pathway matched genes known as pathway-specific, whereas genes with the highest expression variation were common to all pathways. Finally, we showed that the level of cell-to-cell variation is increased in the most immature compartment of hematopoiesis in myelodysplastic syndromes. These data suggest that differentiation could be better conceptualized as a dynamical stochastic process with a transient stage of cellular indetermination. Inhibitory proteins block substrate access by occupying the active site cleft of Bacillus subtilis intramembrane metalloprotease SpoIVFB Intramembrane proteases function in numerous signaling pathways that impact health, but knowledge about regulation of intramembrane proteolysis is limited. We examined inhibition of intramembrane metalloprotease SpoIVFB by proteins BofA and SpoIVFA. We found that BofA residues in and near a predicted transmembrane segment are required for SpoIVFB inhibition, and cross-linking experiments indicated that this transmembrane segment occupies the SpoIVFB active site cleft. SpoIVFA is also required for SpoIVFB inhibition. The inhibitory proteins block access of the substrate N-terminal Proregion to the membrane-embedded SpoIVFB active site, based on additional cross-linking experiments; however, the inhibitory proteins did not prevent interaction between the substrate C-terminal region and the SpoIVFB soluble domain. A structural model was built of SpoIVFB in complex with BofA and parts of SpoIVFA and substrate, using partial homology and constraints from cross-linking and co-evolutionary analyses. The model predicts that conserved BofA residues interact to stabilize a transmembrane segment and a membrane-embedded C-terminal region. SpoIVFA is predicted to bridge the BofA C-terminal region and SpoIVFB, forming a membrane-embedded inhibition complex. Our results reveal a novel mechanism of intramembrane protease inhibition with clear implications for relief from inhibition in vivo and design of inhibitors as potential therapeutics. The renal inflammatory network of nephronophthisis STRUCTURED ABSTRACTO_ST_ABSBACKGROUNDC_ST_ABSThe majority of genetic kidney disease leading to kidney failure is caused by mutations in ciliary genes. How cilia malfunction leads to progressive kidney damage is poorly understood, but recent evidence links ciliopathy genes to CCL2 dependent macrophage recruitment in autosomal dominant polycystic kidney disease (ADPKD), the most studied renal ciliopathy. Whether or not renal inflammation is involved in other renal ciliopathies is unclear.

METHODSWe combined mice models with kidney biopsies and renal epithelial cells sampled from human urine to characterize the renal inflammatory network of nephronophthisis (NPH), the most frequent renal ciliopathy in children.

RESULTSIn human, mutations in cilia genes involved in NPH enhance urine excretion of the chemokine CCL2, causing abnormal macrophage recruitment in kidney tissues from NPH patients. Differing from ADPKD, inactivating Ccl2 specifically in mouse tubular cells does not rescue the NPH phenotype, suggesting that other inflammatory mediators are involved. Using transcriptional data from 2 NPH models, we identify a set of pro-inflammatory cytokines upregulated in this disease, independently of CCL2. The majority of detectable transcripts from this set are specifically upregulated in kidney cells from NPH patients. In line with the function of these cytokines, NPH kidneys show disproportionate neutrophils and T cells infiltrates compared to healthy subject or hypertensive and diabetic chronic kidney disease patients.

CONCLUSIONSThis study reveals that inflammation is a central aspect in human NPH and delineates a specific set of inflammatory mediators that regulates immune cell recruitment in human NPH.

SIGNIFICANCE STATEMENTMutations in genes encoding primary cilia proteins are the leading cause of genetic kidney failure. In autosomal dominant polycystic kidney disease (ADPKD), deregulated cilia signaling leads to kidney infiltration by macrophages through the chemokine CCL2. Little is known about renal inflammation in nephronophthisis (NPH), the most frequent pediatric renal ciliopathy. Using NPH mice models, tissues and cells from NPH patients, we unveil renal inflammation as preeminent feature of NPH. Remarkably, the renal inflammatory evoked by ciliary gene mutations in NPH does not overlap with ADPKD: it is CCL2 independent, involves a prominent recruitment of neutrophils and T cells and a specific cytokine signature. This unforeseen findings strengthen the link between primary cilia and renal inflammation. Functional binding dynamics relevant to the evolution of zoonotic spillovers in endemic and emergent Betacoronavirus strains Comparative functional analysis of the dynamic interactions between various Betacoronavirus mutant strains and broadly utilized target proteins such as ACE2 and CD26, is crucial for a more complete understanding of zoonotic spillovers of viruses that cause diseases such as COVID-19. Here, we employ machine learning to replicated sets of nanosecond scale GPU accelerated molecular dynamics simulations to statistically compare and classify atom motions of these target proteins in both the presence and absence of different endemic and emergent strains of the viral receptor binding domain (RBD) of the S spike glycoprotein. Machine learning was used to identify functional binding dynamics that are evolutionarily conserved from bat CoV-HKU4 to human endemic/emergent strains. Conserved dynamics regions of ACE2 involve both the N-terminal helices, as well as a region of more transient dynamics encompassing K353, Q325 and a novel motif AAQPFLL 386-92 that appears to coordinate their dynamic interactions with the viral RBD at N501. We also demonstrate that the functional evolution of Betacoronavirus zoonotic spillovers involving ACE2 interaction dynamics are likely pre-adapted from two precise and stable binding sites involving the viral bat progenitor strains interaction with CD26 at SAMLI 291-5 and SS 333-334. Our analyses further indicate that the human endemic strains hCoV-HKU1 and hCoV-OC43 have evolved more stable N-terminal helix interactions through enhancement of an interfacing loop region on the viral RBD, whereas the highly transmissible SARS-CoV-2 variants (B.1.1.7, B.1.351 and P.1) have evolved more stable viral binding via more focused interactions between the viral N501 and ACE2 K353 alone. Climate warming changes synchrony of plants and pollinators Climate warming changes the phenology of many species. When interacting organisms respond differently, climate change may disrupt their interactions and affect the stability of ecosystems. Here, we used GBIF occurrence records to examine phenology trends in plants and their associated insect pollinators in Germany since the 1980s. We found strong phenological advances in plants, but differences in the extent of shifts among pollinator groups. The temporal trends in plant and insect phenologies were generally associated with interannual temperature variation, and thus likely driven by climate change. When examining the synchrony of species-level plant-pollinator interactions, their temporal trends differed among pollinator groups. Overall, plant-pollinator interactions become more synchronized, mainly since the phenology of plants responded more strongly to climate change than that of the pollinators. However, if the observed trends continue, many interactions may become more asynchronous again in the future. Our study suggests that climate change affects the phenologies of both plants and insects, and that it also influences the synchrony of plant-pollinator interactions. Receptor-binding domain recombinant protein on alum-CpG induces broad protection against SARS-CoV-2 variants of concern We conducted preclinical studies in mice using a yeast-produced SARS-CoV-2 RBD subunit vaccine candidate formulated with aluminum hydroxide (alum) and CpG deoxynucleotides. This formulation is equivalent to the Corbevax vaccine that recently received emergency use authorization by the Drugs Controller General of India. We compared the immune response of mice vaccinated with RBD/alum to mice vaccinated with RBD/alum+CpG. We also evaluated mice immunized with RBD/alum+CpG and boosted with RBD/alum. Mice were immunized twice intramuscularly at a 21-day interval. Compared to two doses of the /alum formulation, the RBD/alum+CpG vaccine induced a stronger and more balanced Th1/Th2 cellular immune response, with high levels of neutralizing antibodies against the original Wuhan isolate of SARS-CoV-2 as well as the B.1.1.7 (Alpha), B. 1.351 (Beta), B. 1.617.2 and (Delta) variants. Neutralizing antibody titers against the B.1.1.529 (BA.1, Omicron) variant exceeded those in human convalescent plasma after Wuhan infection but were lower than against the other variants. Interestingly, the second dose did not benefit from the addition of CpG, possibly allowing dose-sparing of the adjuvant in the future. The data reported here reinforces that the RBD/alum+CpG vaccine formulation is suitable for inducing broadly neutralizing antibodies against SARS-CoV-2 including variants of concern.