CBEA: Competitive balances for taxonomic enrichmentanalysis Research in human associated microbiomes often involves the analysis of taxonomic count tables generated via high-throughput sequencing. It is difficult to apply statistical tools as the data is high-dimensional, sparse, and compositional. An approachable way to alleviate high-dimensionality and sparsity is to aggregate variables into pre-defined sets. Set-based analysis is ubiquitous in the genomics literature, and has demonstrable impact in improving interpretability and power of downstream analysis. Unfortunately, there is a lack of sophisticated set-based analysis methods specific to microbiome taxonomic data, where current practice often employs abundance summation as a technique for aggregation. This approach prevents comparison across sets of different sizes, does not preserve inter-sample distances, and amplifies protocol bias. Here, we attempt to fill this gap with a new single sample taxon enrichment method that uses a novel log-ratio formulation based on the competitive null hypothesis commonly used in the enrichment analysis literature. Our approach, titled competitive balances for taxonomic enrichment analysis (CBEA), generates sample-specific enrichment scores as the scaled log ratio of the subcomposition defined by taxa within a set and the subcomposition defined by its complement. We provide sample-level significance testing by estimating an empirical null distribution of our test statistic with valid p-values. Herein we demonstrate using both real data applications and simulations that CBEA controls for type I error even under high sparsity and high inter-taxa correlation scenarios. Additionally, it provides informative scores that can be inputs to downstream analyses such as prediction tasks.

Author summaryThe study of human associated microbiomes relies on genomic surveys via high-throughput sequencing. However, microbiome taxonomic data is sparse and high dimensional which prevents the application of standard statistical techniques. One approach to address this problem is to perform analyses at the level of taxon sets. Set-based analysis has a long history in the genomics literature, with demonstrable impact in improving both power and interpretability. Unfortunately, there is limited interest in developing new set-based tools tailored for microbiome taxonomic data given its unique features compared to other  omics data types. We developed a new tool to generate taxon set enrichment scores at the sample level through a novel log-ratio formulation based on the competitive null hypothesis. Our scores can be used for statistical inference at both the sample and population levels, as well as inputs to other downstream analyses such as prediction models. We demonstrate the performance of our method against competing approaches across both real data analyses and simulation studies.
Deep Learning Achieves Neuroradiologist-Level Performance in Detecting Hydrocephalus Requiring Treatment PurposeIn large clinical centers a small subset of patients present with hydrocephalus that requires surgical treatment. We aimed to develop a screening tool to detect such cases from the head MRI with performance comparable to neuroradiologists.

MethodsWe leveraged 496 clinical MRI exams collected retrospectively at a single clinical site from patients referred for any reason. This diagnostic dataset was enriched to have 259 hydrocephalus cases. A 3D convolutional neural network was trained on 16 manually segmented exams (ten hydrocephalus) and subsequently used to automatically segment the remaining 480 exams and extract volumetric anatomical features. A linear classifier of these features was trained on 240 exams to detect cases of hydrocephalus that required treatment with surgical intervention. Performance was compared to four neuroradiologists on the remaining 240 exams. Performance was also evaluated on a separate screening dataset of 451 exams collected from a routine clinical population to predict the consensus reading from four neuroradiologists using images alone. The pipeline was also tested on an external dataset of 31 exams from a 2nd clinical site.

ResultsThe most discriminant features were the Magnetic Resonance Hydrocephalic Index (MRHI), ventricle volume, and the ratio between ventricle and brain volume. At matching sensitivity, the specificity of the machine and the neuroradiologists did not show significant differences for detection of hydrocephalus on either dataset (proportions test, p > 0.05). ROC performance compared favorably with the state-of-the-art (AUC 0.90-0.96), and replicated in the external validation.

ConclusionHydrocephalus cases requiring treatment can be detected automatically from MRI in a heterogeneous patient population based on quantitative characterization of brain anatomy with performance comparable to that of neuroradiologists.
RNF213 loss of function reshapes vascular transcriptome and spliceosome leading to disrupted angiogenesis and aggravated vascular inflammatory responses. RationaleMoyamoya disease (MMD) is a rare cerebrovascular occlusive disease that affects Asian population more often. The pathogenesis of MMD is related to mutation in RNF213 gene. However, why, and how RNF213 mutation leads to MMD is still not fully understood.

ObjectiveAnalyze the impact of RNF213 loss of function on vascular cells and the observed changes correlate with MMD.

Methods and resultsRNF213 KD was conducted in HUVEC and vascular smooth muscle cells (vSMCs). First, HUVEC cells showed alteration of angiogenesis, migration under LPS stimulation, Leukocyte endothelial transmigration, and endothelial-to-vSMCs communication.

Transcriptome analysis revealed downregulation of genes regulating cell division and mitosis. Interestingly, Alternative splicing (AS) analysis revealed hundreds of AS events to occur after RNF213 KD in various types of AS. The alternatively spliced genes showed minimal overlap with transcriptome profiling. Pathway analysis revealed many processes and pathways to be regulated by AS events observed. Transcriptome profiling was also performed after LPS treatment and revealed the basis for increased sensitivity to LPS observed in our analysis. AS changes were also observed after LPS treatment in RNF213 KD HUVEC. Different types of AS showed different patterns of convergence or divergence in terms of the regulated pathways after LPS treatment. Finally, transcriptome and AS analysis was performed in vSMCs and showed various processes that impact vSMCs function and phenotype in RNF213 loss of function.

ConclusionOur data provide a wealth of information on RNF213 gene function in vascular tissues and shed important light on different processes that contribute to MMD pathogenesis. Our results signify an immune-driven and hemodynamically linked process of MMD initiation and progression. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. The impact of obesity; key features of raised intracranial pressure and clinical sequalae in rodents Elevated intracranial pressure (ICP) is observed in many brain disorders. Obesity has been linked to ICP pathogenesis in disorders such as idiopathic intracranial pressure (IIH). We investigated the effect of diet induced obesity (DIO) on ICP and clinically relevant sequelae.

Rats were fed either a control or high fat diet. Following weight gain long term ICP, headache behavior, body composition and retinal outcome were examined. Post-hoc analysis of retinal histology and molecular analysis of choroid plexus and trigeminal ganglion (TG) were performed.

DIO rats demonstrated raised ICP by 55% which correlated with the abdominal fat percentage and increased non-respiratory slow waves, suggestive of altered cerebral compliance. Concurrently, DIO rats demonstrated a specific cephalic cutaneous allodynia which negatively correlated with the abdominal fat percentage. This sensitivity was associated with increased expression of headache markers in TG. Additionally, DIO rats had an in vivo increased retinal nerve fiber layer thickness associated with raised ICP with a subsequent post-hoc demonstration of neuroretinal degeneration.

This study demonstrates for the first time that DIO leads to raised ICP and subsequent clinically relevant symptom development. This novel model of non-traumatic raised ICP could expand the knowledge regarding disorders with elevated ICP such as IIH. Filopodia-mediated basement membrane assembly at pre-invasive tumor boundaries Ductal carcinoma in situ (DCIS) is a pre-invasive stage of breast cancer, where the tumor is encapsulated by a basement membrane (BM). At the invasive phase, the BM barrier is compromised enabling tumor cells to escape into the surrounding stroma. The molecular mechanisms that establish and maintain an epithelial BM barrier in vivo are poorly understood. Myosin-X (MYO10) is a filopodia-inducing motor protein implicated in metastasis and poor clinical outcome in patients with invasive breast cancer (IBC). We compared MYO10 expression in patient-matched normal breast tissue and DCIS lesions and found elevated MYO10 expression in DCIS samples, suggesting that MYO10 might facilitate the transition from DCIS to IBC. Indeed, MYO10 promoted the formation of filopodia and cell invasion in vitro and positively regulated the dissemination of individual cancer cells from IBC lesions in vivo. However, MYO10-depleted DCIS xenografts were, unexpectedly, more invasive. In these xenografts, MYO10 depletion compromised BM formation around the lesions resulting in poorly defined tumor borders and increased cancer cell dispersal into the surrounding stroma. Moreover, MYO10-depleted tumors showed increased EMT-marker-positive cells, specifically at the tumor periphery. We also observed cancer spheroids undergoing rotational motion and recruiting BM components in a filopodia-dependent manner to generate a near-continuous extracellular matrix boundary. Taken together, our data identify a protective role for MYO10 in early-stage breast cancer, where MYO10-dependent tumor cell protrusions support BM assembly at the tumor-stroma interface to limit cancer progression, and a pro-invasive role that facilitates cancer cell dissemination at later stages.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=123 SRC="FIGDIR/small/464987v2_ufig1.gif" ALT="Figure 1">
View larger version (28K):
org.highwire.dtl.DTLVardef@d90a4eorg.highwire.dtl.DTLVardef@133ba9borg.highwire.dtl.DTLVardef@1204a05org.highwire.dtl.DTLVardef@dba254_HPS_FORMAT_FIGEXP  M_FIG C_FIG Highlights- Filopodia sculpt the tumor-proximal stroma in pre-invasive ductal carcinoma in situ (DCIS).
- Filopodia-dependent basement membrane (BM) assembly limits invasive transition of DCIS-like tumors in vivo.
- Loss of MYO10-dependent filopodia impairs BM assembly and induces an EMT-like phenotype at the tumor-stroma interface in vivo.
- MYO10 filopodia are anti-invasive in DCIS but facilitate dissemination in invasive breast cancer. Antiviral activity of molnupiravir precursor NHC against SARS-CoV-2 Variants of Concern (VOCs) and its therapeutic window in a human lung cell model SynopsisO_ST_ABSBackgroundC_ST_ABSThe UK Medicines and Regulatory Healthcare Agency (MHRA) have recently licensed the anti-viral drug, molnupiravir, for use in patients with mild-moderate COVID-19 disease with one or more risk factors for serious illness. Treatment with anti-viral drugs is best initiated early to prevent progression to severe disease, although the therapeutic window for intervention has not yet been fully defined.

ObjectivesThis study aimed to determine the activity of the molnupiravir (NHC) to different SARS-CoV-2 Variants of Concern (VoCs), and to establish the therapeutic window in human lung cell model.

MethodsDose response assays were performed in parallel to determine the IC50 (the concentration of drug required to inhibit virus titre by 50%) of NHC against different variants. Human ACE-2 A549 cells were treated with NHC at different time points either before, during or after infection with SARS- CoV-2.

ResultsHere we demonstrate that {beta}-D-N4-hydroxycytidine (NHC), the active metabolite of molnupiravir, has equivalent activity against four variants of SARS-CoV-2 in a human lung cell line ranging 0.04-0.16{micro}M IC50. Furthermore, we demonstrate that in-vitro activity of the drug is reduced in cells exposed to drug 48 hours after infection.

ConclusionsOne of the main advantages of molnupiravir is that it can be administered orally, and thus given to patients in an out-patient setting. These results support giving the drug early on after diagnosis or even in prophylaxis for individuals with high risk of developing severe disease. Mathematical Modeling and Analysis of Mitochondrial Retrograde Signaling Dynamics Mitochondria, semi-autonomous eukaryotic organelles, participate in energy production and metabolism, making mitochondrial quality control crucial. As most mitochondrial proteins are encoded by nuclear genes, quality control depends on proper mitochondria-nucleus communication, designated mitochondrial retrograde signaling. Early studies focused on retrograde signaling participants and specific gene knockouts. However, mitochondrial signal modulation remains elusive. Using yeast, we simulated signal propagation following mitochondrial damage and proposed a mathematical model based ordinary differential equations. Mitochondrial retrograde signaling decisions were described using a Boolean model. Dynamics were analyzed through an ordinary differential equation-based model and extended to evaluate the model response to noisy damage signals. Simulation revealed localized protein concentration dynamics, including waveforms, frequency response, and robustness under noise. Retrograde signaling is bistable with three localized steady states, and increased damage compromises robustness. We elucidated mitochondrial retrograde signaling, thus providing a basis for drug design against yeast and fungi.

Author SummaryThe yeast RTG pathway regulates mitochondrial metabolism and mitochondrial quality through passing mitochondrial signal to the nucleus to modulate gene expressions. Using microscopic data of RTG proteins translocation from cytosol to nucleus, the parameters were found by fitting 16 knockout conditions simulated from the proposed differential equation-based model. Further, through dose response, ultrasensitivity, frequency response and noise interference, we demonstrated a switch-like property of RTG pathway activation, a capacity charging property of RTG protein translocation, a low pass filter property of signals in different frequencies, and noise amplification due to mitochondrial signal. The proposed mathematical model enables us to understand the dynamics and mechanisms in mitochondrial retrograde signaling in yeast and provides potential antifungal treatment strategies. Precartilage condensation during limb skeletogenesis occurs by tissue phase separation controlled by a bistable cell-state switch with suppressed oscillatory dynamics The tetrapod limb skeleton is initiated in unpatterned limb bud mesenchyme by the formation of precartilage condensations. Here, based on time-lapse videographic analysis of a forming condensation in a high-density culture of chicken limb bud mesenchyme, we observe a phase transition to a more fluidized state for cells within spatial compacted foci (protocondensations that will progress to condensations), as reflected in their spatial confinement, cell-substratum interaction and speed of motion. Previous work showed that galectin-8 and galectin-1A, two proteins of the galactoside-binding galectin family, are the earliest determinants of this process in the chicken limb bud, and that their interactions in forming skeletogenic patterns of condensations can be interpreted mathematically through a reaction-diffusion-adhesion framework. Based on this framework, we use an ordinary differential equation-based approach to analyze the core switching modality of the galectin reaction network and characterize the states of the network independent of the diffusive and adhesive arms of the patterning mechanism. We identify two steady states where the concentrations of both galectins are respectively, negligible, and very high. An explicit Lyapunov function shows that there are no periodic solutions. For sigmoidal galectin production terms, the model exhibits a bistable switch that arises from a monostable state via saddle-node bifurcation. Our model therefore predicts that the galectin network exists in low and high expression states separated in space or time without any intermediate states. This provides a causal basis for the observed outside vs. inside transition observed in the in vitro video data. We performed a quantitative analysis of the distribution of galectin-1A in cultures of condensing chick limb mesenchymal cells and found that the interior of the protocondensations had concentrations of this protein (compared to the immediate exterior) over and above that expected from its higher cell density, consistent with the models predictions. The galectin-based patterning network is thus suggested, on theoretical grounds, to incorporate a core switch independent of any spatial or temporal dynamics, that drives the chondrogenic cell state transition in limb skeletogenesis. Applications of the indole-alkaloid gramine modulate the assembly of individual members of the barley rhizosphere microbiota Microbial communities proliferating at the root-soil interface, collectively referred to as the rhizosphere microbiota, represent an untapped beneficial resource for plant growth, development and health. Integral to a rational manipulation of the microbiota for sustainable agriculture is the identification of the molecular determinants of these communities. In plants, biosynthesis of allelochemicals is centre stage in defining inter-organismal relationships in the environment. Intriguingly, this process has been moulded by domestication and breeding selection. The indole-alkaloid gramine, whose occurrence in barley (Hordeum vulgare L.) is widespread among wild genotypes but has been counter selected in several modern varieties, is a paradigmatic example of this phenomenon. This prompted us to investigate how exogenous applications of gramine impacted on the rhizosphere microbiota of two, gramine-free, elite barley varieties grown in a reference agricultural soil. High throughput 16S rRNA gene amplicon sequencing revealed that applications of gramine interfere with the proliferation of a subset of soil microbes with a relatively broad phylogenetic assignment. Strikingly, growth of these bacteria appeared to be rescued by barley plants in a genotype- and dosage-independent manner. In parallel, we discovered that host recruitment cues can interfere with the impact of gramine application in a host genotype-dependent manner. Interestingly, this latter effect displayed a bias for members of the phyla Proteobacteria. These initial observations indicate that gramine can act as a determinant of the prokaryotic communities inhabiting the root-soil interface. 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. Development and Characterization of PAR-Trackers: New Tools for Detecting Poly(ADP-ribose) In Vitro and In Vivo AO_SCPLOWBSTRACTC_SCPLOWADP-ribosylation (ADPRylation) is a reversible post-translation modification resulting in the covalent attachment of ADP-ribose (ADPR) moieties on substrate proteins. Naturally-occurring protein motifs and domains, including WWEs, PBZs, and macrodomains, act as "readers" for protein-linked ADPR. Although recombinant, antibody-like ADPR detection reagents containing these readers have facilitated the detection of ADPR, they are limited in their ability to capture the dynamic nature of ADPRylation. Herein, we describe and characterize a set of poly(ADP-ribose) (PAR) Trackers (PAR-Ts) - optimized dimerization-dependent or split-protein reassembly PAR sensors in which a naturally occurring PAR binding domain, WWE, was fused to both halves of dimerization-dependent GFP (ddGFP) or split Nano Luciferase (NanoLuc), respectively. We demonstrate that these new tools allow the detection and quantification of PAR levels in extracts, living cells, and living tissues with greater sensitivity, as well as temporal and spatial precision. Importantly, these sensors detect changes in cellular ADPR levels in response to physiological cues (e.g., hormone-dependent induction of adipogenesis without DNA damage), as well as xenograft tumor tissues in living mice. Our results indicate that PAR Trackers have broad utility for detecting ADPR in many different experimental and biological systems. Spontaneous differentiation leads to emergence of hybrid cell states relate to poor prognosis in oral cancer PurposeCellular dynamics between phenotypically heterogeneous subpopulations of cancer cells within individual tumor is shown to be responsible for drug tolerance and overall poor prognosis; however, evidences were largely missing in oral cancer. Therefore, this study was undertaken to describe the dynamic phenotypic states among oral cancer cells, its influence on transcriptomic heterogeneity as well as its clinical significance.

Experimental DesignWe multiplexed phenotypic markers of putative oral-stem-like cancer cells (SLCCs) and characterized diversity among CD44-positive oral cancer cell subpopulations with respect to distinct expression of CD24 and aldehyde dehydrogenase (ALDH)-activity in multiple cell lines. Population trajectories were characterized by Markov model and cell states were defined based on the population specific RNA sequencing (RNAseq). ssGSEA based gene expression signatures were explore for prognostic significance.

ResultsOral cancer cells followed two distinct patterns of spontaneous repopulation dynamics with stochastic inter-conversions on  ALDH-axis, however a strict non-interconvertible transition on  CD24-axis. Interestingly, plastic  ALDH-axis was harnessed to enrich ALDHHigh subpopulations in response to Cisplatin treatment, to adapt a drug tolerant state. Phenotype-specific RNAseq results suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate states of stemness within the differentiating oral cancer cells. Further, survival analysis with each subpopulation-specific gene signature strongly suggested that the cell-state dynamics may act as possible mechanism to drive ITH, resulting in poor prognosis in patient.

ConclusionsOur results emphasized the prognostic power of the population dynamics in oral cancer. Importantly, we have described the phenotypic-composition of heterogeneous subpopulations critical for global tumor behaviour in oral cancer; which is a prerequisite knowledge important for precision treatment, however largely lacking for most solid tumors.

Graphical AbstractWe have characterized diversity among CD44-positive oral cancer cells lines with respect to distinct expression of CD24 and ALDH-activity. Subpopulations showed stochastic inter-conversions on ALDH-axis but a strict non-interconvertible transition of CD24Low to CD24High phenotype, even in response to chemotherapy-induced stress. RNAseq study suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate alternate states of stemness within the differentiating oral cancer cells. The described population dynamics demonstrtaed influence tumor behaviour possibly by increasing intratumoral heterogeneity in aggressive oral tumors.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=112 SRC="FIGDIR/small/457509v3_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@1dc6a7org.highwire.dtl.DTLVardef@dd0948org.highwire.dtl.DTLVardef@18c9603org.highwire.dtl.DTLVardef@cca8fb_HPS_FORMAT_FIGEXP  M_FIG C_FIG Translational relevanceIntratumoral heterogeneity (ITH) has been the clinically important factor, impacting aggressive cancer behaviour, drug tolerance and overall poor prognosis. Recent high-throughput studies have provided better cellular and molecular resolution of ITH; however, the prerequisite knowledge which defines the composition of subpopulations critical for global tumor behaviour is majorly lacking for most of the solid tumors. By combining phenotypic markers, we have defined four subpopulations of oral cancer cells. These subpopulations showed stochastic inter-conversions as well as a strict non-interconvertible transition among them to acheive heterogeneity. Importantly, transcriptional states of each subpopulations indicated a clinically relevant signatures for patient prognosis. Also, we observed interconversions of these subpopulations in response to Cisplatin to accumulate drug-tolerant cell state, as rapid and reversible strategy to respond to chemotherapy induced stress. Thus, the characteristics of described phenotypic subgroups may be translated to the clinic for estimating the extent of intratumoral heterogeneity in oral cancer patients. Estimating animal density with camera traps: a practitioner's guide of the REST model Camera traps are increasingly popular in wildlife research and have the potential to be reliable and cost-effective for estimating animal density. Although density estimation using this automatic technique has long been restricted to species with individually recognizable markings, several analytical approaches have been proposed to target animals lacking such markings (i.e., unmarked populations). Among these approaches, the random encounter and staying time (REST) model may be an efficient and cost-effective approach, even though the procedures for the implementations have not yet been shared with researchers. This paper presents a working protocol for implementing the REST model our research group has been developing. We also present the R code to perform parameter estimation with a maximum likelihood and Bayesian approach. We suggest that this model has potential for further development. We strongly hope that this paper will encourage many researchers to use the REST for density estimation in a wide variety of species in various habitats and make a significant contribution to advancing wildlife conservation and management. Learning and motor inhibitory control in crows and domestic chickens Cognitive abilities allow animals to navigate through complex, fluctuating environments. In the present study, we tested the performance of a captive group of eight crows, Corvus corone, and ten domestic chickens, Gallus gallus domesticus, in the cylinder task, as a test of motor inhibitory control, and reversal learning as a measure of learning ability and behavioural flexibility. Four crows and nine chickens completed the cylinder task, eight crows and six chickens completed the reversal learning experiment. Crows performed better in the cylinder task compared to chickens. In the reversal learning experiment, species did not significantly differ in the number of trials until learning criterion was reached. Performance in the reversal learning experiment did not correlate with performance in the cylinder task in chickens. Our results suggest crows to possess better motor inhibitory control compared to chickens. In contrast, learning performance in a reversal learning task did not differ between the species, indicating similar levels of behavioural flexibility. Interestingly, we describe notable individual differences in performance. We stress the importance not only to compare cognitive performance between species but also between individuals of the same species when investigating the evolution of cognitive skills. 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. Identification of age and ethnicity specific gene expression biomarkers for immune aging Human immune system functions over an entire lifetime, yet how and why the immune system becomes less effective with age are not well understood. Here, we characterize peripheral blood mononuclear cells transcriptome from 172 healthy adults with 21~90 years of age using RNA-seq and the weighted gene correlation network analyses (WGCNA). These data reveal a set of insightful gene expression modules and representative gene biomarkers for human immune system aging from Asian and Caucasian ancestry, respectively. Among them, the aging-specific modules show an age-related gene expression variation spike around early-seventies. In addition, it is not known whether Asian and Caucasian immune systems go through similar gene expression changes throughout their lifespan, and to what extent these aging-associated changes are shared among ethnicities. We find the top hub genes including NUDT7, CLPB, OXNAD1 and MLLT3 are shared between Asian and Caucasian aging related modules and further validated in human PBMCs from different age groups. Overall, the impact of age and race on transcriptional variation elucidated from this study provide insights into the transcriptional driver of immune aging. Protein-coding potential of RNAs measured by potentially translated island scores Recent studies have identified numerous RNAs that are functionally both coding and noncoding. However, the sequence characteristics that determine bifunctionality remain largely unknown. In this study, we developed and tested a potentially translated island (PTI) score, defined as the occupancy of the longest open reading frame (ORF) among all putative ORFs. We found that this score correlated with translation, including noncoding RNAs. In bacteria and archaea, coding and noncoding transcripts had narrow distributions of high and low PTI scores, respectively, whereas those of eukaryotes showed relatively broader distributions, with considerable overlap between coding and noncoding transcripts. The extent of overlap positively and negatively correlated with the mutation rates of genomes and effective population sizes of species, respectively. These overlaps were significantly increased in threatened species. In macroevolution, the appearance of the nucleus and multicellularity seem to have influenced the overlap of PTI score distributions, so that the probability of the existence of bifunctional RNAs is increased in eukaryotes. In mammalian testes, we observed an enrichment of noncoding RNAs with high PTI scores, which are candidates for bifunctional RNAs. These results suggest that the decrease in population size and the emergence of testes in eukaryotic multicellular organisms allow for the stable existence of bifunctional RNAs, consequently increasing the probability of the birth of novel coding and non-coding RNAs. TAS-Seq: a robust and sensitive amplification method for bead-based scRNA-seq Single-cell RNA-sequencing (scRNA-seq) is valuable for analyzing cellular heterogeneity. Cell composition accuracy is critical for analyzing cell-cell interaction networks from scRNA-seq data. We developed terminator-assisted solid-phase cDNA amplification and sequencing (TAS-Seq) for scRNA-seq based on a terminator, terminal transferase, and nanowell/bead-based scRNA-seq platform; TAS-Seq showed high tolerance to variations in the terminal transferase reaction, which complicate the handling of existing terminal transferase-based scRNA-seq methods. In murine and human lung samples, TAS-Seq yielded scRNA-seq data that were highly correlated with flow-cytometric data, showing higher gene-detection sensitivity and more robust detection of important cell-cell interactions and expression of growth factors/interleukins in cell subsets than 10X Chromium v2 and Smart-seq2. Expanding TAS-Seq application will improve understanding and atlas construction of lung biology at the single-cell level. First detection of herpesvirus and prevalence 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. 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. Common variants in breast cancer risk loci predispose to distinct tumor subtypes BackgroundGenome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER), but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

MethodsAmong 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

ResultsEighty-five of 173 variants were associated with at least one tumor feature (false discovery rate <5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at P<0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

ConclusionThis report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction. DNA Strand Breaks and Gaps Target Retroviral Binding and Integration Integration into a host genome is essential for retrovirus infection and is catalyzed by a nucleoprotein complex (Intasome) containing the viral integrase (IN) and reverse transcribed (RT) copy DNA (cDNA). Previous studies demonstrated DNA site recognition limited intasome integration. Using single molecule Forster resonance energy transfer (smFRET), we show Prototype Foamy Virus (PFV) intasomes pause at DNA strand breaks and gaps. The break/gap discontinuities are similar to base excision repair (BER) lesion-processing intermediates, which affect retrovirus integration in vivo. Pausing targeted site-directed integration at the break/gap without inducing intasome conformational alterations. An 8-oxo-guanine lesion normally processes by BER and a G/T mismatch or a +T nucleotide insertion that induce flexibility or a bend in the DNA backbone did not promote intasome pausing or targeted integration. These results suggest that repair intermediates can modulate dynamic intasome-DNA interactions which target retroviral integration. 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. Production and characterization of virus-free, CRISPR-CAR T cells capable of inducing solid tumor regression BackgroundChimeric antigen receptor (CAR) T cells traditionally harbor viral vector-based sequences that encode the CAR transgene in the genome. These T cell products have yet to show consistent anti-tumor activity in patients with solid tumors. Further, viral vector manufacturing is resource intensive, suffers from batch-to-batch variability, and includes several animal components, adding regulatory and supply chain pressures.

MethodsAnti-GD2 CAR T cells were generated using CRISPR/Cas9 within nine days using recombinant Cas9 protein and nucleic acids, without any viral vectors or animal components. The CAR was specifically targeted to the T Cell Receptor Alpha Constant gene (TRAC). T cell products were characterized at the level of the genome, transcriptome, proteome, and secretome using CHANGE-seq, scRNA-seq, spectral cytometry, and ELISA assays. Functionality was evaluated in vivo in an NSG xenograft neuroblastoma model.

ResultsIn comparison to traditional retroviral CAR T cells, virus-free CRISPR CAR (VFC-CAR) T cells exhibit TRAC-targeted genomic integration of the CAR transgene, elevation of transcriptional and protein characteristics associated with a memory phenotype, and low tonic signaling prior to infusion arising in part from the the knockout of the TCR. Upon exposure to the GD2 target antigen, anti-GD2 VFC-CAR T cells exhibited specific cytotoxicity against GD2+ cells in vitro and induced solid tumor regression in vivo, with robust homing, persistence, and low exhaustion against a human neuroblastoma xenograft model.

ConclusionsThis proof-of-principle study leveraging virus-free genome editing technology could enable flexible manufacturing of clinically relevant, high-quality CAR T cells to treat cancers, including solid tumors. Fast ATP-dependent Subunit Rotation in Reconstituted FoF1-ATP Synthase Trapped in Solution FoF1-ATP synthases are ubiquitous membrane-bound, rotary motor enzymes that can catalyze ATP synthesis and hydrolysis. Their enzyme kinetics are controlled by internal subunit rotation, by substrate and product concentrations, by mechanical inhibitory mechanisms, but also by the electrochemical potential of protons across the membrane. Single-molecule Forster resonance energy transfer (smFRET) has been used to detect subunit rotation within FoF1-ATP synthases embedded in freely diffusing liposomes. We now report that kinetic monitoring of functional rotation can be prolonged from milliseconds to seconds by utilizing an Anti-Brownian electrokinetic trap (ABEL trap). These extended observation times allowed us to observe fluctuating rates of functional rotation for individual FoF1-liposomes in solution. Broad distributions of ATP-dependent catalytic rates were revealed. The buildup of an electrochemical potential of protons was confirmed to limit the maximum rate of ATP hydrolysis. In the presence of ionophores or uncouplers, the fastest subunit rotation speeds measured in single reconstituted FoF1-ATP synthases were 180 full rounds per second. This was much faster than measured by biochemical ensemble averaging, but not as fast as the maximum rotational speed reported previously for isolated single F1 complexes uncoupled from the membrane-embedded Fo complex. Further application of ABEL trap measurements should help resolve the mechanistic causes of such fluctuating rates of subunit rotation.

TOC GRAPHIC

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=108 SRC="FIGDIR/small/436299v5_ufig1.gif" ALT="Figure 1">
View larger version (35K):
org.highwire.dtl.DTLVardef@1bd800org.highwire.dtl.DTLVardef@228aeorg.highwire.dtl.DTLVardef@1bce39forg.highwire.dtl.DTLVardef@3f01a_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOFigure 1.C_FLOATNO C_FIG