DeepFlu: a deep learning approach for forecasting symptomatic influenza A infection based on pre-exposure gene expression Background and ObjectiveNot everyone gets sick after an exposure to influenza A viruses (IAV). Although KLRD1 has been identified as a potential biomarker for influenza susceptibility, it remains unclear whether forecasting symptomatic flu infection based on pre-exposure host gene expression might be possible.

MethodTo examine this hypothesis, we developed DeepFlu using the state-of-the-art deep learning approach on the human gene expression data infected with IAV subtype H1N1 or H3N2 viruses to forecast who would catch the flu prior to an exposure to IAV.

ResultsThe results indicated that such forecast is possible and, in other words, gene expression could reflect the strength of host immunity. In the leave-one-person-out cross-validation, DeepFlu based on deep neural network outperformed the models using convolutional neural network, random forest, or support vector machine, achieving 70.0% accuracy, 0.787 AUROC, and 0.758 AUPR for H1N1 and 73.8% accuracy, 0.847 AUROC, and 0.901 AUPR for H3N2. In the external validation, DeepFlu also reached 71.4% accuracy, 0.700 AUROC, and 0.723 AUPR for H1N1 and 73.5% accuracy, 0.732 AUROC, and 0.749 AUPR for H3N2, surpassing the KLRD1 biomarker. In addition, DeepFlu which was trained only by pre-exposure data worked the best than by other time spans and mixed training data of H1N1 and H3N2 did not necessarily enhance prediction. DeepFlu is available at https://github.com/ntou-compbio/DeepFlu.

ConclusionsDeepFlu is a prognostic tool that can moderately recognize individuals susceptible to the flu and may help prevent the spread of IAV.
Assessment and Optimization of the Interpretability of Machine Learning Models Applied to Transcriptomic Data Explainable artificial intelligence aims to interpret how the machine learning models make decisions, and many model explainers have been developed in the computer vision field. However, the understandings of the applicability of these model explainers to biological data are still lacking. In this study, we comprehensively evaluated multiple explainers by interpreting pretrained models of predicting tissue types from transcriptomic data, and by identifying top contributing genes from each sample with the greatest impacts on model prediction. To improve the reproducibility and interpretability of results generated by model explainers, we proposed a series of optimization strategies for each explainer on two different model architectures of Multilayer Perceptron (MLP) and Convolutional Neural Network (CNN). We observed three groups of explainer and model architecture combinations with high reproducibility. Group II, which contains three model explainers on aggregated MLP models, identified top contributing genes in different tissues that exhibited tissue-specific manifestation and were potential cancer biomarkers. In summary, our work provides novel insights and guidance for exploring biological mechanisms using explainable machine learning models.
Electrophysiological correlates of mood and reward dynamics in human adolescents Despite its omnipresence in everyday interactions and its importance for mental health, mood and its neuronal underpinnings are poorly understood. Computational models can help identify parameters affecting self-reported mood during mood induction tasks. Here we test if computationally modelled dynamics of self-reported mood during monetary gambling can be used to identify trial-by-trial variations in neuronal activity. To this end, we shifted mood in healthy (N=24) and depressed (N=30) adolescents by delivering individually tailored reward prediction errors whilst recording magnetoencephalography (MEG) data. Following a pre-registered analysis, we hypothesize that expectation (defined by previous reward outcomes) would be predictive of beta-gamma oscillatory power (25-40Hz), a frequency shown to modulate to reward feedback. We also hypothesize that trial variations in the evoked response to the presentation of gambling options and in source localized responses to reward feedback. Through our multilevel statistical analysis, we found confirmatory evidence that beta-gamma power is positively related to reward expectation during mood shifts, with possible localized sources in the posterior cingulate cortex. We also confirmed reward prediction error to be predictive of trial-level variations in the response of the paracentral lobule and expectation to have an effect on the cerebellum after presentation of gambling options. To our knowledge, this is the first study to relate fluctuations in mood on a minute timescale to variations in neural oscillations with noninvasive electrophysiology.

Significance StatementBrain mechanisms underlying mood and its relationship with changes in reward contingencies in the environment are still elusive but could have a strong impact on our understanding and treatment of debilitating mood disorders. Building on a previously proposed computational mood model we use multilevel statistical models to find relationship between trial-by-trial variations in model components of mood and neural responses to rewards measured with non-invasive electrophysiology (MEG). Through confirmatory analysis we show that it is possible to observe relationships between trial variations in neural responses and computational parameters describing mood dynamics. Identifying the dynamics of mood and the neural processes it affects could pave the way for more effective neuromodulation treatments. Superinfection exclusion: a viral strategy with short-term benefits and long-term drawbacks Viral superinfection occurs when multiple viral particles subsequently infect the same host. In nature, several viral species are found to have evolved diverse mechanisms to prevent superinfection (superinfection exclusion) but how this strategic choice impacts the fate of mutations in the viral population remains unclear. Using stochastic simulations, we find that genetic drift is suppressed when superinfection occurs, thus facilitating the fixation of beneficial mutations and the removal of deleterious ones. Interestingly, we also find that the competitive (dis)advantage associated with variations in life history parameters is not necessarily captured by the viral growth rate for either infection strategy. Putting these together, we then show that a mutant with superinfection exclusion will easily overtake a superinfecting population even if the latter has a much higher growth rate. Our findings suggest that while superinfection exclusion can negatively impact the long-term adaptation of a viral population, in the short-term it is ultimately a winning strategy.

Author summaryViral social behaviour has recently been receiving increasing attention in the context of ecological and evolutionary dynamics of viral populations. One fascinating and still relatively poorly understood example is superinfection or co-infection, which occur when multiple viruses infect the same host. Among bacteriophages, a wide range of mechanisms have been discovered that enable phage to prevent superinfection (superinfection exclusion) even at the cost of using precious resources for this purpose. What is the evolutionary impact of this strategic choice and why do so many phages exhibit this behaviour? Here, we conduct an extensive simulation study of a phage population to address this question. In particular, we investigate the fate of viral mutations arising in an environment with a constant supply of bacterial hosts designed to mimic a "turbidostat," as these are increasingly being used in laboratory evolution experiments. Our results show that allowing superinfection in the long-term yields a population which is more capable of adapting to changes in the environment. However, when in direct competition, mutants capable of preventing superinfection experience a very large advantage over their superinfecting counterparts, even if this ability comes at a significant cost to their growth rate. This indicates that while preventing superinfection can negatively impact the long-term prospects of a viral population, in the short-term it is ultimately a winning strategy. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. Mechanisms of splicing substrate definition by modulation of pre-mRNA 3D structural scaffolds Coordination of different serine-arginine-rich (SR) proteins - a class of critical splicing activators - facilitates recognition of the highly degenerate mammalian splice signal sequences against the background sequences. Yet, the mechanistic details of their actions remain unclear. Here we show that cooperative binding of SR proteins to exonic and intronic motifs remodels the pre-mRNA 3D structural scaffold. The scaffold generated by the pre-mRNA-specific combinations of different SR proteins in an appropriate stoichiometry is recognized by U1 snRNP - the 5' splice site recognizing factor. A large excess of U1 snRNP particles displaces majority of the bound SR protein molecules from the remodeled pre-mRNA. A higher than optimal stoichiometry of SR proteins occludes the binding sites on the pre-mRNA, raising the U1 snRNP levels required for SR protein displacement and potentially impeding spliceosome assembly. This novel splicing-related function of U1 snRNP is important for distinguishing the substrate and the non-substrate by U2AF65 - the primary 3' splice site-recognizing factor. Overall, this work elucidates early regulatory steps of mammalian splicing substrate definition by SR proteins. Modelling and assessing additional transmission routes for porcine reproductive and respiratory syndrome virus: vehicle movements and feed ingredients Accounting for multiple modes of livestock disease dissemination in epidemiological models remains a challenge. We developed and calibrated a mathematical model for transmission of porcine reproductive and respiratory syndrome virus (PRRSV), tailored to fit nine modes of between-farm transmission pathways including: farm-to-farm proximity (local transmission), contact network of batches of pigs transferred between farms (pig movements), re-break probabilities for farms with previous PRRSV outbreaks, with the addition of four different contact networks of transportation vehicles (vehicles to transport pigs to farms, pigs to markets, feed and crew) and the amount of animal by-products within feed ingredients (e.g. animal fat or meat and bone meal). The model was calibrated on weekly PRRSV outbreaks data. We assessed the role of each transmission pathway considering the dynamics of specific types of production (i.e., sow, nursery). Although our results estimated that the networks formed by transportation vehicles were more densely connected than the network of pigs transported between-farms, pig movements and farm proximity were the main PRRSV transmission routes regardless of farm types. Among the four vehicle networks, vehicles transporting pigs to farms explained a large proportion of infections, sow = 20.9%; nursery = 15%; and finisher = 20.6%. The animal by-products showed a limited association with PRRSV outbreaks through descriptive analysis, and our model results showed that the contribution of animal fat contributed only 2.5% and meat and bone meal only 0.03% of the infected sow farms. Our work demonstrated the contribution of multiple routes of PRRSV dissemination, which has not been deeply explored before. It also provides strong evidence to support the need for cautious, measured PRRSV control strategies for transportation vehicles and further research for feed by-products modeling. Finally, this study provides valuable information and opportunities for the swine industry to focus effort on the most relevant modes of PRRSV between-farm transmission. Vesalius: high-resolution in silico anatomization of Spatial Transcriptomic data using Image Analysis Characterization of tissue architecture promises to deliver insights into development, cell communication and disease. In silico spatial domain retrieval methods have been developed for spatial transcriptomics (ST) data assuming transcriptional similarity of neighboring barcodes. However, domain retrieval approaches with this assumption cannot work in complex tissues composed of multiple cell types. This task becomes especially challenging in cellular resolution ST methods. We developed Vesalius to decipher tissue anatomy from ST data by applying image processing technology. Vesalius uniquely detected territories composed of multiple cell types and successfully recovered tissue structures in high-resolution ST data including in mouse brain, embryo, liver, and colon. Utilizing this tissue architecture, Vesalius identified tissue morphology specific gene expression and regional specific gene expression changes for astrocytes, interneuron, oligodendrocytes, and entorhinal cells in the mouse brain. Effect of a low concentration in chloride, sodium and potassium on oocyte maturation, oocyte hydration, ovulation and egg quality in rainbow trout. Water salinity is an important environmental factor known to have detrimental effects on salmonid reproduction, mostly when migrating female broodfish are held in sea water. In contrast, data obtained in freshwater are scarce and the impact of low water salinity, and more specifically of low water concentrations in sodium, chloride and potassium, during reproduction in freshwater is currently unknown. For this reason, and because ion and water fluxes are critical for the final steps of the female gamete formation, including oocyte hydration and ovulation, the aim of the present study was to investigate the impact of low salinity water on final oocyte maturation, ovulation and, ultimately, on egg quality, using rainbow trout as a physiological model and relevant aquaculture species.

Fish from the same commercial strain were raised either in a site characterized by low concentrations of Na+, K+, and Cl- ions in the water or in a closely located control site exhibiting higher concentration in these elements. Egg quality and duration of final oocyte maturation were investigated using innovative phenotyping tools such as automatic assessment of egg viability using the VisEgg system and non-invasive echograph-based monitoring of final oocyte maturation duration, respectively. Oocyte hydration during final oocyte maturation and after ovulation was also investigated. Finally, molecular phenotyping was performed using real-time PCR-based monitoring of several key players of final oocyte maturation and ovulation associated with ion and water transport, inflammation, proteolytic activity, and coagulation. Oocyte hydration and gene expression data were analyzed in the light of the duration of final oocyte maturation.

Here we show that low water salinity (i.e., low water concentration in chloride, sodium and potassium) negatively influences final oocyte maturation, ovulation and, ultimately, egg quality. Low water salinity triggered delayed ovulation and lower oocyte viability. When investigating the impact of low water salinity on final oocyte maturation duration, individuals presenting the most severe phenotypes exhibited impaired oocyte hydration and abnormally reduced gene expression levels of several key players of the ovulatory process. While the under expression of water (i.e., aquaporins) and ion (i.e., solute carriers) transporters is consistent with impaired oocyte hydration, our observations also indicate that the overall ovulatory gene expression program is disrupted. Our results raise the question of the mechanisms underlying the negative influence of low salinity water on the dynamics of the preovulatory phase, on the control of the oocyte homeostasis, including hydration, and on the overall success of the maturation-ovulation process.

HighlightsO_LILow water salinity impairs final oocyte maturation and egg quality in rainbow trout
C_LIO_LILow water salinity induces delayed ovulation and impaired oocyte hydration
C_LIO_LILow water salinity induces a dysregulation of several key ovulatory genes
C_LIO_LIMonitoring of final oocyte maturation can be performed using ultrasound staging
C_LI Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI Thermal modulation of Zebrafish exploratory statistics reveals constraints on individual behavioral variability BackgroundVariability is a hallmark of animal behavior. It contributes to survival by endowing individuals and populations with the capacity to adapt to ever-changing environmental conditions. Intra-individual variability is thought to reflect both endogenous and exogenous modulations of the neural dynamics of the central nervous system. However, how variability is internally regulated and modulated by external cues remains elusive. Here we address this question by analyzing the statistics of spontaneous exploration of freely swimming zebrafish larvae, and by probing how these locomotor patterns are impacted when changing the water temperatures within an ethologically relevant range.

ResultsWe show that, for this simple animal model, five short-term kinematic parameters - interbout interval, turn amplitude, travelled distance, turn probability and orientational flipping rate - together control the long-term exploratory dynamics. We establish that the bath temperature consistently impacts the means of these parameters, but leave their pairwise covariance unchanged. These results indicate that the temperature merely controls the sampling statistics within a well-defined kinematic space delineated by this robust statistical structure. At a given temperature, individual animals explore the behavioral space over a timescale of tens of minutes, suggestive of a slow internal state modulation that could be externally biased through the bath temperature. By combining these various observations into a minimal stochastic model of navigation, we show that this thermal modulation of locomotor kinematics results in a thermophobic behavior, complementing direct gradient-sensing mechanisms.

ConclusionsThis study establishes the existence of a well-defined locomotor space accessible to zebrafish larvae during spontaneous exploration, and quantifies self-generated modulation of locomotor patterns. Intra-individual variability reflects a slow diffusive-like probing of this space by the animal. The bath temperature in turn restricts the sampling statistics to sub-regions, endowing the animal with basic thermophobicity. This study suggests that in Zebrafish, as well as in other ectothermic animals, ambient temperature could be used to efficiently manipulate internal states in a simple and ethological way. Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. The "Alu-ome" shapes the epigenetic environment of regulatory elements controlling cellular defense Promoters and enhancers are sites of transcription initiation (TSSs) and carry active histone modifications, including H3K4me1, H3K4me3, and H3K27ac. Yet, the principles governing the boundaries of such regulatory elements are still poorly characterized. Alu elements are good candidates for a boundary function, being highly abundant in gene-rich regions, while essentially excluded from regulatory elements. Here, we show that the interval from the TSS to the first upstream Alu accommodates all H3K4me3 and most H3K27ac marks, while excluding DNA methylation. Remarkably, the average length of this intervals greatly varies in-between tissues, being longer in stem-and shorter in immune-cells. The very shortest TSS-to-Alu intervals were observed at promoters active in T cells, particularly at immune genes, correlating with RNA polymerase II transcription through the first Alu and accumulation of H3K4me1 signal on this first Alu. Finally, DNA methylation at first-Alus was found to evolved with age, regressing from young to middle-aged, then recovering later in life. Thus, the first Alus upstream of TSSs appear as dynamic boundaries marking the transition from DNA methylation to active histone modifications at regulatory elements, while also participating in the recording of immune gene transcriptional events by positioning H3K4me1-modified nucleosomes. 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. A bacterial effector counteracts host autophagy by promoting degradation of an autophagy component Beyond its role in cellular homeostasis, autophagy plays anti- and pro-microbial roles in host-microbe interactions, both in animals and plants. One prominent role of anti-microbial autophagy is to degrade intracellular pathogens or microbial molecules, in a process termed xenophagy. Consequently, microbes evolved mechanisms to hijack or modulate autophagy to escape elimination. Although well-described in animals, the extent to which xenophagy contributes to plant-bacteria interactions remains unknown. Here, we provide evidence that Xanthomonas campestris pv. vesicatoria (Xcv) suppresses host autophagy by utilizing type-III effector XopL. XopL interacts with and degrades the autophagy component SH3P2 via its E3 ligase activity to promote infection. Intriguingly, XopL is targeted for degradation by defense-related selective autophagy mediated by NBR1/Joka2, revealing a complex antagonistic interplay between XopL and the host autophagy machinery. Our results implicate plant antimicrobial autophagy in depletion of a bacterial virulence factor and unravels an unprecedented pathogen strategy to counteract defense-related autophagy. Transcriptomic analysis of temporal shifts in berry development between two grapevine cultivars of the Pinot family reveals potential genes controlling ripening time BackgroundGrapevine cultivars of the Pinot family represent clonally propagated mutants with major phenotypic and physiological differences, such as different colour or shifted ripening time, as well as changes in important viticultural traits. Specifically, the cultivars  Pinot Noir (PN) and  Pinot Noir Precoce (PNP, early ripening) flower at the same time, but vary in the beginning of berry ripening (veraison) and, consequently, harvest time. In addition to genotype, seasonal climatic conditions (i.e. high temperatures) also affect ripening times. To reveal possible regulatory genes that affect the timing of veraison onset, we investigated differences in gene expression profiles between PN and PNP throughout berry development with a closely meshed time series and over two separate years.

ResultsThe difference in the duration of berry formation between PN and PNP was quantified to be approximately two weeks under the growth conditions applied, using plant material with a proven PN and PNP clonal relationship. Clusters of co-expressed genes and differentially expressed genes (DEGs) were detected which reflect the shift in the timing of veraison onset. Functional annotation of these DEGs fit to observed phenotypic and physiological changes during berry development. In total, we observed 3,342 DEGs in 2014 and 2,745 DEGs in 2017 between PN and PNP, with 1,923 DEGs across both years. Among these, 388 DEGs were identified as veraison-specific and 12 were considered as berry ripening time regulatory candidates. The expression profiles revealed two candidate genes for ripening time control which we designated VviRTIC1 and VviRTIC2 (VIT_210s0071g01145 and VIT_200s0366g00020, respectively). These genes likely contribute the phenotypic differences observed between PN and PNP.

ConclusionsMany of the 1,923 DEGs show highly similar expression profiles in both cultivars if the patterns are aligned according to developmental stage. In our work, putative genes differentially expressed between PNP and PN which could control ripening time as well as veraison-specific genes were identified. We point out connections of these genes to molecular events during berry development and discuss potential candidate genes which may control ripening time. Two of these candidates were observed to be differentially expressed in the early berry development phase. Several down-regulated genes during berry ripening are annotated as auxin response factors / ARFs. Conceivably, general changes in auxin signaling may cause the earlier ripening phenotype of PNP. Cellular feedback to organotelluranes displays modulation of antioxidant proteins gene expression Organotelluranes RT3 and RT4 are thiol reagents that induce mitochondrial transition pore (MTP) opening in a sensitive and insensitive manner to cyclosporin A. Although RT3 and RT4 promote glutathione depletion, paradoxically, they are also an efficient antioxidant for membrane lipids. These compounds antagonistic effects elicited the challenging question of how the gene expression of antioxidant enzymes would respond to treatment with these compounds. The influence of RT3 and RT4 on antioxidant enzyme expression was investigated in cultured aortic smooth muscle cells (ASMC). RT3 and RT4 promoted disruption of ionic calcium homeostasis, mitochondrial transmembrane potential ({Delta}{Psi}), and cell death in a dose-dependent manner. The cell death mechanisms responded qualitatively to the increase of the organotellurane concentration and changed from apoptosis to necrosis. RT3 and RT4 increased the expression of thioredoxin significantly. RT3 also increased the expression of glutaredoxin and glutathione peroxidase, slightly the catalase expression without significant effects on SOD expression. The results are consistent with GSH and protein thiol depletion and discussed based on the cell toxicity mechanism exhibited by these compounds.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=150 SRC="FIGDIR/small/425411v4_ufig1.gif" ALT="Figure 1">
View larger version (73K):
org.highwire.dtl.DTLVardef@52bb82org.highwire.dtl.DTLVardef@19476d9org.highwire.dtl.DTLVardef@6aa99org.highwire.dtl.DTLVardef@1aa605e_HPS_FORMAT_FIGEXP  M_FIG C_FIG Metabolic reconstitution by a gnotobiotic microbiota varies over the circadian cycle The capacity of the intestinal microbiota to degrade otherwise indigestible diet components is known to greatly improve the recovery of energy from food. This has led to the hypothesis that increased digestive efficiency may underlie the contribution of the microbiome to obesity. OligoMM12-colonized gnotobiotic mice have a consistently higher fat-mass than germ-free or fully colonized counterparts. We therefore investigated their food intake, digestion efficiency, energy expenditure and respiratory quotient using a novel isolator-housed metabolic cage system which allows long-term measurements without contamination risk. This demonstrated that microbiota- released calories are perfectly balanced by decreased food intake in fully colonized versus gnotobiotic OligoMM12 and germ-free mice fed a standard chow diet, i.e., microbiota-released calories can in fact be well-integrated into appetite control. We also observed no significant difference in energy expenditure per gram lean mass between the different microbiota groups, suggesting that cumulative very small differences in energy balance, or altered energy storage must underlie fat accumulation in OligoMM12 mice. Consistent with altered energy storage, major differences were observed in the type of respiratory substrates used in metabolism over the circadian cycle: in germ-free mice the respiratory exchange ratio was consistently lower than that of fully colonized mice at all times of day, indicative of more reliance on fat and less on glucose metabolism. Intriguingly the RER of OligoMM12-colonized gnotobiotic mice phenocopied fully colonized mice during the dark (active/eating) phase but phenocopied germ-free mice during the light (fasting/resting) phase. Further, OligoMM12-colonized mice showed a germ-free-like drop in liver glycogen storage during the light cycle and both liver and plasma metabolomes of OligoMM12 mice clustered closely with germ-free mice. This implies the existence of microbiota functions that are required to maintain normal host metabolism during the resting/fasting phase of circadian cycle, and which are absent in the OligoMM12 consortium. Omicron booster in ancestral strain vaccinated mice augments protective immunities against both the Delta and Omicron variants A booster vaccination is called for constraining the evolving epidemic of SARS-CoV-2. However, the necessity of a new COVID-19 vaccine is currently unclear. To compare the effect of an Omicron-matched S DNA vaccine and an ancestral S DNA vaccine in boosting cross-reactive immunities, we firstly immunized mice with two-dose of a DNA vaccine encoding the spike protein of the ancestral Wuhan strain. Then the mice were boosted with DNA vaccines encoding spike proteins of either the Wuhan strain or the Omicron variant. Specific antibody and T cell responses were measured at 4 weeks post boost. Our data showed that the Omicron-matched vaccine efficiently boosted RBD binding antibody and neutralizing antibody responses against both the Delta and the Omicron variants. Of note, antibody responses against the Omicron variant elicited by the Omicron-matched vaccine were much stronger than those induced by the ancestral S DNA vaccine. Meanwhile, CD8+ T cell responses against both the ancestral Wuhan strain and the Omicron strain also tended to be higher in mice boosted by the Omicron-matched vaccine than those in mice boosted with the ancestral S DNA vaccine, albeit no significant difference was observed. Our findings suggest that an Omicron-matched vaccine is preferred for boosting cross-reactive immunities. Force correction analysis method for derivation of multidimensional free energy landscapes from adaptively-biased replica simulations A methodology is proposed for the calculation of multidimensional free-energy landscapes of molecular systems, based on analysis of multiple Molecular Dynamics trajectories wherein adaptive biases have been applied to enhance the sampling of different collective variables. In this approach, which we refer to as Force Correction Analysis Method (FCAM), local averages of the total and biasing forces are evaluated post-hoc, and the latter are subtracted from the former to obtain unbiased estimates of the mean force across collective-variable space. Multidimensional free-energy surfaces and minimum free-energy pathways are then derived from integration of the mean force landscape through kinetic Monte Carlo algorithm. To evaluate the proposed method, a series of numerical tests and comparisons with existing approaches were carried out for small molecules, peptides, and proteins, based on all-atom trajectories generated with standard, concurrent and replica-exchange Metadynamics in collective-variable spaces ranging from one- to six-dimensional. The tests confirm the correctness of the FCAM formulation and demonstrate that calculated mean forces and free energies converge rapidly and accurately, outperforming other methods used to unbias this kind of simulation data.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=102 SRC="FIGDIR/small/431654v5_ufig1.gif" ALT="Figure 1">
View larger version (45K):
org.highwire.dtl.DTLVardef@1412289org.highwire.dtl.DTLVardef@1425d53org.highwire.dtl.DTLVardef@1fcf5d0org.highwire.dtl.DTLVardef@6e6941_HPS_FORMAT_FIGEXP  M_FIG TOC/Abstract Graphic

C_FIG 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 two cancer stem cell-like populations in triple-negative breast cancer xenografts Gene expression analysis at the single-cell level by next generation sequencing has revealed the existence of clonal dissemination and microheterogeneity in cancer metastasis. The current spatial analysis technologies can elucidate the heterogeneity of cell-cell interactions in situ. To reveal the regional and expressional heterogeneity in primary tumors and metastases, we performed transcriptomic analysis of microtissues dissected from a triple-negative breast cancer (TNBC) cell line MDA-MB-231 xenograft model with our automated tissue microdissection punching technology. This multiple-microtissue transcriptome analysis revealed three cancer cell-type clusters in the primary tumor and axillary lymph node metastasis, two of which were cancer stem cell (CSC)-like clusters (CD44/MYC-high, HMGA1-high). Reanalysis of public single-cell RNA-seq (scRNA-seq) datasets confirmed that the two CSC-like populations existed both in TNBC xenograft models and TNBC patients. In addition, the gene signature of the HMGA1-high CSC-like cluster has the potential to serve as a novel biomarker for diagnosis. The diversity of these multiple CSC-like populations may cause differential anticancer drug resistance, increasing the difficulty of curing this cancer. 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. Location and function of Plasmodium kinesins: key roles in parasite proliferation, polarity, and transmission Kinesins are microtubule-based motors important in cell division, motility, polarity, and intracellular transport in many eukaryotes. However, they are poorly studied in the divergent eukaryotic pathogens-Plasmodium spp., the causative agents of malaria, which manifest atypical aspects of cell division and plasticity of morphology throughout the lifecycle in both mammalian and mosquito hosts. Here we describe a genome-wide screen of Plasmodium kinesins, revealing diverse subcellular locations and functions in spindle assembly, axoneme formation and cell morphology. Surprisingly, only kinesin-13 is essential for growth in the mammalian host while the other eight kinesins are required during the proliferative and invasive stages of parasite transmission through the mosquito vector. In-depth analyses of kinesin-13 and kinesin-20 revealed functions in microtubule dynamics during apical cell polarity formation, spindle assembly, and axoneme biogenesis. These findings help us to understand the importance of microtubule motors and may be exploited to discover new therapeutic interventions against malaria.