Characterising ontogeny of quantity discrimination in zebrafish A sense of non-symbolic numerical magnitudes is widespread in the animal kingdom and has been documented in adult zebrafish. Here we investigated the ontogeny of this ability using a group size preference task in juvenile zebrafish. Fish showed group size preference from 21 days post fertilization (dpf) and reliably chose the larger group when presented with discriminations of between 1 vs. 3, 2 vs. 5 and 2 vs. 3 conspecifics but not 2 vs. 4 conspecifics. When the ratio between the number of conspecifics in each group was maintained at 1:2, fish could discriminate between 1 vs. 2 individuals and 3 vs. 6, but again, not when given a choice between 2 vs. 4 individuals. These findings are in agreement with studies in other species suggesting the systems involved in quantity representation do not operate separately from other cognitive mechanisms. Rather they suggest quantity processing in fish may be the result of an interplay between attentional, cognitive and memory-related mechanisms as in humans and other animals. Our results emphasise the potential of the use of zebrafish to explore the genetic and neural processes underlying the ontogeny and function of number cognition.
Accuracy and Power Analysis of Social Networks Built From Count Data O_LIPower analysis is used to estimate the probability of correctly rejecting a null hypothesis for a given statistical model and dataset. Conventional power analyses assume complete information, but the stochastic nature of behavioural sampling can mean that true and estimated networks are poorly correlated. Power analyses do not currently take the effect of sampling into account. This could lead to inaccurate estimates of statistical power, potentially yielding misleading results.
C_LIO_LIHere we develop a method for computing network correlation: the correlation between an estimated social network and its true network, using a Gamma-Poisson model of social event rates for networks constructed from count data. We use simulations to assess how the level of network correlation affects the power of nodal regression analyses. We also develop a generic method of power analysis applicable to any statistical test, based on the concept of diminishing returns.
C_LIO_LIWe demonstrate that our network correlation estimator is both accurate and moderately robust to its assumptions being broken. We show that social differentiation, mean social event rate, and the harmonic mean of sampling times positively impacts the strength of network correlation. We also show that the required level of network correlation to achieve a given power level depends on many factors, but that 0.80 network correlation usually corresponds to around 0.80 power for nodal regression in ideal circumstances.
C_LIO_LIWe provide guidelines for using our network correlation estimator to verify the accuracy of networks built from count data, and to conduct power analysis. This can be used prior to data collection, in post hoc analyses, or even for subsetting networks in dynamic network analysis. The network correlation estimator and custom power analysis methods have been made available as an R package.
C_LI
Fatal COVID-19 ARDS associated with incomplete AEC1 differentiation from the transitional state without senescence or fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. Assigning metabolic rate measurements to torpor and euthermy in heterothermic endotherms: 'torpor', a new package for R. Torpor is a state of controlled reduction of metabolic rate (M) in endotherms. Assigning measurements of M to torpor or euthermy can be challenging, especially when the difference between euthermic M and torpid M is small, in species defending a high minimal body temperature in torpor, in thermolabile species, and slightly below the thermoneutral zone (TNZ). Here, we propose a novel method for distinguishing torpor from euthermy. We use the variation in M measured during euthermic rest and torpor at varying ambient temperatures (Ta) to objectively estimate the lower critical temperature (Tlc) of the TNZ and to assign measurements to torpor, euthermic rest or rest within TNZ. In addition, this method allows the prediction of M during euthermic rest and torpor at varying Ta, including resting M within the TNZ. The present method has shown highly satisfactory results using 28 published sets of metabolic data obtained by respirometry on 26 species of mammals. Ultimately, this novel method aims to facilitate analysis of respirometry data in heterothermic endotherms. Finally, the development of the associated R-package (torpor) will enable widespread use of the method amongst biologists.

Summary statementThe presented method and its associated R-package (torpor) enable the assignment of metabolic rate measurements to torpor or euthermy, ultimately improving the standardization of respirometry analyses in heterotherms. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

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

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. A Machine Learning Approach to Identify Small Molecule Inhibitors of Secondary Nucleation inalpha-Synuclein Aggregation Drug development is an increasingly active area of application of machine learning methods, due to the need to overcome the high attrition rates of conventional drug discovery pipelines. This issue is especially pressing for neurodegenerative diseases where very few disease-modifying drugs have been approved. To address this problem, we describe a machine learning approach to identify specific inhibitors of the proliferation of -synuclein aggregates through secondary nucleation, a process that has been implicated in Parkinsons disease and related synucleinopathies. We use a combination of docking simulations followed by machine learning to first identify initial hit compounds and then explore the chemical space around these compounds. Our results demonstrate that this approach leads to the identification of novel chemical matter with an improved hit rate and potency over conventional similarity search approaches. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. Dynamic interaction of BcsD subunit in type I bacterial cellulose synthase Cellulose is a promising biological material for supporting sustainable human life. This natural polymer is synthesized by cellulose synthase, a protein complex in the cell membrane. Cellulose synthase in bacteria is a hetero-subunit complex, and its subunit organization varies widely depending on the species. In the type I bacterial cellulose synthase complex, the BcsD (bacterial cellulose synthase D) protein is believed to play an important role in producing cellulose with long slender fiber morphology and high crystallinity, given the phenotype of the bcsD-deficient mutant and the specific existence of the type I operon in bacterial species synthesizing crystalline cellulose microfibrils such as Acetobacter. In this study, we successfully established a heterogeneously expressed Bcs protein in Escherichia coli as an experimental model and conducted biochemical studies for the BcsD protein and the other three major subunits of bacterial cellulose synthase, BcsA, BcsB, and BcsC. It has been shown that the BcsD protein interacts with the functionally required minimal subunits of the BcsAB complex, as well as the BcsC protein. Furthermore, it was shown that BcsD interacts with the BcsAB complex in two modes: direct protein-protein interactions and indirect interactions through the product cellulose. The former and latter modes represent the basal and active states of the type I bacterial cellulose synthase, respectively. This dynamic behavior of the BcsD protein is important for the type I bacterial cellulose synthase complex to regulate the crystallization process of cellulose. Multifunctional RNA-binding proteins influence mRNA abundance and translational efficiency of distinct sets of target genes RNA-binding proteins (RBPs) can regulate more than a single aspect of RNA metabolism. We searched for such previously undiscovered multifunctionality within a set of 143 RBPs, by defining the predictive value of RBP abundance for the transcription and translation levels of known RBP target genes across 80 human hearts. This led us to newly associate 27 RBPs with cardiac translational regulation in vivo. Of these, 21 impacted both RNA expression and translation, albeit for virtually independent sets of target genes. We highlight a subset of these, including G3BP1, PUM1, UCHL5, and DDX3X, where dual regulation is achieved through differential affinity for target length, by which separate biological processes are controlled. Similar to the RNA helicase DDX3X, the known splicing factors EFTUD2 and PRPF8 - all identified as multifunctional RBPs by our analysis - selectively influence target translation rates depending on 5 UTR structure. Our analyses identify dozens of RBPs as being multifunctional and pinpoint potential novel regulators of translation, postulating unanticipated complexity of protein-RNA interactions at consecutive stages of gene expression. Is Pisum sativum a good model species for the study of root responses to neighbours and barriers in soil? A bayesian hierarchical meta-analysis Plant-plant competition is ubiquitous in nature. However, studying below ground behaviour of roots has always posed certain difficulties. Pea (Pisum sativum L.) has become a sort of model species for ecological studies about how plant roots respond to neighbouring plant roots and barriers in soil. However, published results point in several different directions. This has sometimes been interpreted as pea having sophisticated context dependent responses that can change in complex ways depending on its surroundings. To explore this further, here, we combine the result of five new experiments with published results to examine 18 unique experiments from 7 different studies for a total of 254 replicates. We used a Bayesian hierarchical meta-analysis approach to estimating the likely effect size from the available data, as well as quantify heterogeneity among different experiments, studies and cultivars. The posterior distributions show that, at the coarsest possible scale of total root production, it is unlikely that P. sativum root growth is influenced by either neighbours or barriers to root growth imposed by the walls of pots that vary in volume. We suggest that further work should consider repeating experiments that reported finer scale root plasticity in pea at the rhizosphere scale, and also consider alternative model species for the study of plant root responses to external cues. 3D Echo Planar Time-resolved Imaging (3D-EPTI) for ultrafast multi-parametric quantitative MRI Multi-parametric quantitative MRI has shown great potential to improve the sensitivity and specificity of clinical diagnosis and to enhance our understanding of complex brain processes, but suffers from long scan time especially at high spatial resolution. To address this long-standing challenge, we introduce a novel approach, termed 3D Echo Planar Time-resolved Imaging (3D-EPTI), which significantly increases the acceleration capacity of MRI sampling, and provides high acquisition efficiency for multi-parametric MRI. This is achieved by exploiting the spatiotemporal correlation of MRI data at multiple timescales through new encoding strategies within and between efficient continuous readouts. Specifically, an optimized spatiotemporal CAIPI encoding within the readouts combined with a radial-block sampling strategy across the readouts enables an acceleration rate of 800 folds in the k-t space. A subspace reconstruction was employed to resolve thousands of artifact-free high-quality multi-contrast images spaced at a time interval of ~1 ms. We have demonstrated the ability of 3D-EPTI to provide robust and repeatable whole-brain simultaneous T1, T2, T2*, PD and B1+ mapping at high isotropic resolution within minutes (e.g., 1-mm isotropic resolution in 3 minutes), and to enable submillimeter multi-parametric imaging to study detailed brain structures.

HighlightsO_LIUltra-fast acquisition for 3D multi-parametric quantitative MRI.
C_LIO_LISimultaneous T1 T2 T2* PD and B1+ mapping.
C_LIO_LI3-minute scan at 1-mm isotropic resolution with whole-brain coverage.
C_LIO_LIMulti-parametric mapping at 700-m isotropic resolution in 10 minutes.
C_LIO_LIRepeatable quantification and cortical-depth analysis.
C_LI The obstacles and potential clues of prime editing applications in tomato, a dicot plant Precision genome editing is highly desired for crop improvement. The recently emerged CRISPR/Cas technology offers great potential applications in precision plant genome engineering. A prime editing (PE) approach combining a reverse transcriptase (RT) with a Cas9 nickase and a "priming" extended guide RNA has shown a high frequency for precise genome modification in mammalian cells and several plant species. Nevertheless, the applications of the PE approach in dicot plants are still limited and inefficient. We designed and tested prime editors for precision editing of a synthetic sequence in a transient assay and for desirable alleles of 10 loci in tomato by stable transformation. Our data obtained by targeted deep sequencing also revealed only low PE efficiencies in both the tobacco and tomato systems. Further assessment of the activities of the PE components uncovered potential reasons for the inefficiency of the PE complexes. Modifying the pegRNA sequences by shortening or introducing mismatches to the primer binding sequences (PBS) in order to reduce their melting temperatures (Tm) did not enhance the PE efficiency at the SlBMP21, SlALC and SlALS1 loci. Our data show challenges of PE approach in tomato, indicating a further improvement of the PE system for the successful applications such as use of improved expression systems. Our work provides an important clue for the successful application of the PE approach in crop improvement. Loci for prediction of penicillin and tetracycline susceptibility in Neisseria gonorrhoeae: a genome wide association study BackgroundWhile Neisseria gonorrhoeae poses an urgent public health threat because of increasing antimicrobial resistance, much of the circulating population remains susceptible to historical treatment regimens. Point-of-care diagnostics that report susceptibility could allow for reintroduction of these regimens, but development of such diagnostics has been limited to ciprofloxacin, for which susceptibility can be predicted from a single locus.

MethodsWe assembled a dataset of 12,045 N. gonorrhoeae genomes with phenotypic resistance data for tetracycline (n = 3,611) and penicillin (n = 6,935). Using conditional genome wide association studies (GWAS), we sought to define genetic variants associated with susceptibility to penicillin and tetracycline. We evaluated the sensitivity and specificity of these variants for predicting susceptibility and non-resistance in our collection of gonococcal genomes.

FindingsIn our conditional penicillin GWAS, the presence of a genetic variant defined by a non-mosaic penA allele without an insertion at codon 345 was significantly associated with penicillin susceptibility and had the highest negative effect size of significant variants (p = 5.0 x 10-14, {beta} = -2.5). In combination with the absence of blaTEM, this variant predicted penicillin susceptibility with high specificity (99.8%) and modest sensitivity (36.7%). For tetracycline, the wild type allele at rpsJ codon 57, encoding valine, was significantly associated with tetracycline susceptibility (p = 5.6 x 10-16, {beta} = -1.6) after conditioning on the presence of tetM. The combination of rpsJ codon 57 allele and tetM absence predicted tetracycline susceptibility with high specificity (97.2%) and sensitivity (88.7%).

InterpretationAs few as two genetic loci can predict susceptibility to penicillin and tetracycline in N. gonorrhoeae with high specificity. Molecular point-of-care diagnostics targeting these loci have the potential to increase available treatments for gonorrhea.

FundingNational Institute of Allergy and Infectious Diseases, the National Science Foundation, and the Smith Family Foundation

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSWe searched PubMed with the terms "Neisseria gonorrhoeae" and "diagnostic" or "assay" plus "penicillin" or "tetracycline" for reports in any language published up to July 1, 2021. We additionally searched for "Neisseria gonorrhoeae" and "genome wide association study". We found that previously proposed molecular diagnostics for penicillin and tetracycline susceptibility either exclusively focused on plasmid-mediated resistance (i.e., targeting blaTEM or tetM) or did not include variants in genes encoding antibiotic targets (e.g., did not include penA or rpsJ). Targets for molecular surveillance have focused on resistance-associated alleles rather than susceptibility-associated alleles. We did not find any previous penicillin or tetracycline conditional genome wide association studies (GWAS) in N. gonorrhoeae.

Added value of this studyTo identify targets for molecular diagnostics that predict penicillin and tetracycline susceptibility, we conducted GWAS conditioning on the presence of plasmid-mediated resistance determinants to detect chromosomal loci with the highest association with susceptibility. We discovered a sequence (penA_01) that differentiates susceptible isolates from those with a resistance-associated insertion at codon 345 and from those with mosaic penA alleles, which is associated with penicillin susceptibility. We also found that rpsJ codon 57 was the chromosomal locus contributing the most to tetracycline susceptibility. The combination of these chromosomal loci and the absence of plasmid encoded determinants predicts penicillin and tetracycline susceptibility with high specificity in both a large global collection of N. gonorrhoeae and a validation dataset consisting of recently published genomes from CDCs Gonococcal Isolate Surveillance Program (GISP) surveillance collected in 2018.

Implications of all the available evidenceThe chromosomal loci penA_01 and rpsJ codon 57 in combination with plasmid loci blaTEM and tetM are candidates for the development of point-of-care molecular diagnostics for penicillin and tetracycline susceptibility. The loci may be combined with the currently available ciprofloxacin susceptibility diagnostics to predict susceptibility to multiple antibiotics. Additionally, our study suggests that conditional GWAS focused on variants associated with susceptibility may be a promising approach to identify minimal sets of loci for molecular diagnostics and surveillance. Differential metabolism of arsenicals regulates Fps1-mediated arsenite transport Arsenic is an environmental toxin that exists mainly as pentavalent arsenate and trivalent arsenite. Both forms activate the yeast SAPK Hog1, but with different consequences. We describe a mechanism by which cells distinguish between these arsenicals through one-step metabolism to differentially regulate the bidirectional glycerol channel Fps1, an adventitious port for arsenite. Cells exposed to arsenate reduce it to thiol-reactive arsenite, which modifies a set of cysteine residues in target proteins; whereas cells exposed to arsenite metabolize it to methylarsenite, which modifies an additional set of cysteine residues. Hog1 becomes arsenylated, which prevents it from closing Fps1. However, this block is overcome in cells exposed to arsenite through methylarsenylation of Acr3, an arsenite efflux pump that we found also regulates Fps1 directly. This adaptation allows cells to restrict arsenite entry through Fps1, but also allows its exit when produced from arsenate exposure. These results have broad implications for understanding how SAPKs activated by diverse stressors can drive stress-specific outputs.

Summary (for JCB table of contents)

Jongmin Lee, David E. Levin

Lee and Levin investigate the question of how different stressors can drive divergent outputs from an active SAPK. The work describes the mechanism by which two forms of arsenic that both activate the yeast SAPK Hog1 differentially regulate the glycerol channel Fps1. Identifying and correcting for misspecifications in GWAS summary statistics and polygenic scores Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, due to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched LD references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power. 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 of Offspring from Azoospermic Mice with Meiotic Failure: Precise Biparental Meiosis within Halved Oocytes While the large volume of mammalian oocytes is necessary for embryo development, it can lead to error-prone chromosomal segregation during meiosis. Conversely, we hypothesized that smaller oocytes would have a great unidentified potential to stabilize unstable meiosis and improve the development of the resultant embryos. Here, we show that reducing ooplasmic volume can rescue highly error-prone fertilization using primary spermatocytes by preventing segregation errors of chromosomes during biparental meiosis. High-resolution live-imaging analysis revealed that erroneous chromosome segregation occurred in most (90%) spermatocyte-injected oocytes of normal size, but could be ameliorated to 40% in halved oocytes. The birth rate improved remarkably from 1% to 19% (P < 0.0001). Importantly, this technique enabled the production of offspring from azoospermic mice with spermatocyte arrest caused by STX2 deficiency, an azoospermia factor also found in humans. Thus, contrary to popular opinion, oocytes inherently have a strong potential for precise meiotic divisions, which can be evoked by reduction of the ooplasmic volume. Their potential might help rescue cases of untreatable human azoospermia with spermatocyte arrest. Host phenology can drive the evolution of intermediate virulence strategies in some parasites AO_SCPLOWBSTRACTC_SCPLOWThe traditional mechanistic trade-offs resulting in a negative correlation between transmission and virulence are the foundation of nearly all current theory on the evolution of parasite virulence. Several ecological factors have been shown to modulate the optimal virulence strategies predicted from mechanistic trade-off models, but these ecological factors have not yet been shown to be sufficient to explain the intermediate virulence strategies observed in any natural system. The timing of seasonal activity, or phenology, is a common factor that influences the types and impact of many ecological interactions but is difficult to incorporate into virulence evolution studies. We develop a mathematical model of a disease system with seasonal host activity to study the evolutionary consequences of host phenology on the virulence of obligate-killer parasite. Results from this model demonstrated that seasonal host activity is sufficient to drive the evolution of intermediate parasite virulence in some types of natural disease systems, even when a traditional mechanistic trade-off between transmission and virulence is not assumed in the modeling framework. The optimal virulence strategy in these systems can be determined by both the duration of the host activity period as well as the variation in the host emergence timing. Parasites with low virulence strategies are favored in environments with long host activity periods and in environments in which hosts emerge synchronously. The results demonstrate that host phenology can be sufficient to select for intermediate optimal virulence strategies, providing an alternative mechanism to account for virulence evolution in some natural systems. A novel nongenetic murine model of hyperglycemia and hyperlipidemia-associated aggravated atherosclerosis ObjectiveAtherosclerosis, the main pathology underlying cardiovascular diseases is accelerated in diabetic patients. Genetic mouse models require breeding efforts which are time-consuming and costly. Our aim was to establish a new nongenetic model of inducible metabolic risk factors that mimics hyperlipidemia, hyperglycemia, or both and allows the detection of phenotypic differences dependent on the metabolic stressor(s).

Methods and ResultsWild-type mice were injected with gain-of-function PCSK9D377Y (proprotein convertase subtilisin/kexin type 9) mutant adeno-associated viral particles (AAV) and streptozotocin and fed either a high-fat diet (HFD) for 12 or 20 weeks or a high-cholesterol/high-fat diet (Paigen diet, PD) for 8 weeks. To evaluate atherosclerosis, two different vascular sites (aortic sinus and the truncus of the brachiocephalic artery) were examined in the mice. Combined hyperlipidemic and hyperglycemic (HGHCi) mice fed a HFD or PD displayed characteristic features of aggravated atherosclerosis when compared to hyperlipidemia (HCi HFD or PD) mice alone. Atherosclerotic plaques of HGHCi HFD animals were larger, showed a less stable phenotype (measured by the increased necrotic core area, reduced fibrous cap thickness, and less -SMA-positive area) and had more inflammation (increased plasma IL-1{beta} level, aortic pro-inflammatory gene expression, and MOMA-2-positive cells in the BCA) after 20 weeks of HFD. Differences between the HGHCi and HCi HFD models were confirmed using RNA-seq analysis of aortic tissue, revealing that significantly more genes were dysregulated in mice with combined hyperlipidemia and hyperglycemia than in the hyperlipidemia-only group. The HGHCi-associated genes were related to pathways regulating inflammation (increased Cd68, iNos, and Tnfa expression) and extracellular matrix degradation (Adamts4 and Mmp14). When comparing HFD with PD, the PD aggravated atherosclerosis to a greater extent in mice and showed plaque formation after 8 weeks. Hyperlipidemic and hyperglycemic mice fed a PD (HGHCi PD) showed less collagen (Sirius red) and increased inflammation (CD68-positive cells) within aortic plaques than hyperlipidemic mice (HCi PD). HGHCi-PD mice represent a directly inducible hyperglycemic atherosclerosis model compared with HFD-fed mice, in which atherosclerosis is severe by 8 weeks.

ConclusionWe established a nongenetically inducible mouse model allowing comparative analyses of atherosclerosis in HCi and HGHCi conditions and its modification by diet, allowing analyses of multiple metabolic hits in mice. Tobramycin-induced secretion of P. aeruginosa 5' tRNA-fMet halves suppresses lung inflammation via AGO2 gene silencing Although inhaled tobramycin increases lung function in people with cystic fibrosis (pwCF), the density of P. aeruginosa in the lungs is only modestly reduced by tobramycin; hence, the mechanism whereby tobramycin improves lung function is unclear. Here, we demonstrate that tobramycin increases the abundance of two 5' tRNA-fMet halves in outer membrane vesicles (OMVs) secreted by P. aeruginosa and that the 5' tRNA-fMet halves reduce IL-8 secretion by CF bronchial epithelial cells (CF-HBECs). In mouse lung, the 5' tRNA-fMet halves attenuate KC secretion and neutrophil recruitment. We also report that the 5' tRNA-fMet halves suppress pro-inflammatory network gene expression by an Argonaut 2 (AGO2)-mediated gene silencing mechanism, thereby reducing IL-8 secretion in CF-HBECs. Moreover, tobramycin reduces the IL-8 concentration and neutrophil content in bronchoalveolar lavage fluid of pwCF. Thus, we conclude that tobramycin improves lung function in part by reducing chronic inflammation and neutrophil-mediated lung damage in pwCF. Dynamics of SARS-CoV-2 mutations reveals regional-specificity and similar trends of N501 and High-Frequency mutation N501Y in different levels of control measures Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1 058 020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3 % of global COVID-19 cases and determined three types of mutation dynamics: High-Frequency, Medium-Frequency, and Low-Frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that Medium-Frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of High-Frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with High-Frequency or not-High-Frequency and both follows similar trends in different levels of stringency. Adversarial domain translation networks for fast and accurate integration of large-scale atlas-level single-cell datasets The rapid emergence of large-scale atlas-level single-cell RNA-seq datasets presents remarkable opportunities for broad and deep biological investigations through integrative analyses. However, harmonizing such datasets requires integration approaches to be not only computationally scalable, but also capable of preserving a wide range of fine-grained cell populations. We created Portal, a unified framework of adversarial domain translation to learn harmonized representations of datasets. With innovation in model and algorithm designs, Portal achieves superior performance in preserving biological variation during integration, while achieving integration of millions of cells in minutes with low memory consumption. We show that Portal is widely applicable to integrating datasets across samples, platforms and data types (including scRNA-seq, snRNA-seq and scATAC-seq). Finally, we demonstrate the power of Portal by applying it to the integration of cross-species datasets with limited shared information among them, elucidating biological insights into the similarities and divergences in the spermatogenesis process among mouse, macaque and human. Alpha Oscillatory Activity Causally Linked to Working Memory Retention: Insights from Online Phase-corrected Closed-loop Transcranial Alternating Current Stimulation (tACS) Although previous studies have reported correlations between alpha oscillations and the "retention" sub-process of working memory (WM), causal evidence has been limited in human neuroscience due to the lack of delicate modulation of human brain. Conventional tACS is not suitable for demonstrating the causal evidence for parietal alpha oscillations in WM retention because of its inability to modulate brain oscillations within a short period (i.e., the retention sub-process). Here, we developed an online phase-corrected closed-loop transcranial alternating current stimulation (tACS) system capable of precisely correcting for the phase differences between tACS and concurrent endogenous oscillations. This system permits both up- and down-regulation of brain oscillations at the target stimulation frequency within a short stimulation period, and is here applied to empirically demonstrate that parietal alpha oscillations causally relate to WM retention. Our experimental design included both in-phase and anti-phase alpha-tACS applied to 39 participants during the retention sub-processes of a modified Sternberg paradigm. Compared to in-phase alpha-tACS, anti-phase alpha-tACS decreased both WM performance and alpha activity. Moreover, the in-phase tACS-induced changes in WM performance were positively correlated with alpha oscillatory activity. These findings strongly support a causal link between alpha oscillations and WM retention, and illustrate the broad application prospects of phase-corrected tACS. Autophagosomes fuse to phagosomes and facilitate the degradation of apoptotic cells in Caenorhabditis elegans Autophagosomes are double-membrane intracellular vesicles that degrade protein aggregates, intracellular organelles, and other cellular components. During the development of the nematode Caenorhabditis elegans, many somatic and germ cells undergo apoptosis. These cells are engulfed and degraded by their neighboring cells. We discovered a novel role of autophagosomes in facilitating the degradation of apoptotic cells using a real-time imaging technique. Specifically, the double-membrane autophagosomes in engulfing cells are recruited to the surfaces of phagosomes containing apoptotic cells and subsequently fuse to phagosomes, allowing the inner vesicle to enter the phagosomal lumen. Mutants defective in the production of autophagosomes display significant defects in the degradation of apoptotic cells, demonstrating the importance of autophagosomes to this process. The signaling pathway led by the phagocytic receptor CED-1, the adaptor protein CED-6, and the large GTPase dynamin (DYN-1) promotes the recruitment of autophagosomes to phagosomes. Moreover, the subsequent fusion of autophagosomes with phagosomes requires the functions of the small GTPase RAB-7 and the HOPS complex components. Further observations suggest that autophagosomes provide apoptotic cell-degradation activities in addition to and in parallel of lysosomes. Our findings reveal that, unlike the single-membrane, LC3-associated phagocytosis (LAP) vesicles reported for mammalian phagocytes, the canonical double-membrane autophagosomes facilitate the clearance of C. elegans apoptotic cells. These findings add autophagosomes to the collection of intracellular organelles that contribute to phagosome maturation, identify novel crosstalk between the autophagy and phagosome maturation pathways, and discover the upstream signaling molecules that initiate this crosstalk. How do stomata evolve on the road to C4 photosynthesis? C4 photosynthesis optimizes plants carbon and water relations, allowing high photosynthetic rate with low stomatal conductance. Stomata have long been believed as a part of C4 syndrome. However, it remains unclear how stomata traits evolved along the path from C3 to C4. Stomatal patterning was examined in Flaveria genus, a model for studying C4 evolution. Comparative, transgenic and semi-in-vitro experiments were used to study molecular basis that underlies stomatal traits along C4 evolution. Novel results: the evolution from C3 to C4 species through intermediate species is accompanied by a stepwise rather than an abrupt change in the stomatal traits. The initial change occurs near Type II and dramatic change occurs at the C4-like species. On the road to C4, stomata become less in number but bigger in size and changes in stomatal density dominates changes in maximum stomatal conductance (gsmax). The reduction of FSTOMAGEN expression underlies altered gsmax between Flaveria species with different photosynthetic pathways and likely occurs in other C4 lineages. Our study provides insight into the pattern, mechanism and role of stomatal evolution along the road towards C4. This work highlights the stomatal traits in the current C4 evolutionary model and the co-evolution of photosynthetic pathway and stomata. CRISPR-based kinome-screening revealed MINK1 as a druggable player to rewire 5FU-resistance in OSCC through AKT/MDM2/p53 axis Cisplatin, 5FU and docetaxel (TPF) are the most common chemotherapy regimen used for advanced OSCC. However, many cancer patients experience relapse, continued tumor growth, and spread due to drug resistance, which leads to treatment failure and metastatic disease. Here, using a CRISPR/Cas9 based kinome knockout screening, Misshapen-like kinase 1 (MINK1) is identified as an important mediator of 5FU resistance in OSCC. Analysis of clinical samples demonstrated significantly higher MINK1 expression in the tumor tissues of chemotherapy non-responder as compared to chemotherapy responders. The nude mice and zebrafish xenograft experiments indicate that knocking out MINK1 restores 5FU mediated cell death in chemoresistant OSCC. An antibody based phosphorylation array screen revealed MINK1 as a negative regulator of p53. Mechanistically, MINK1 modulates AKT phosphorylation at Ser473, which enables p-MDM2 (Ser 166) mediated degradation of p53. We also identified lestaurtinib as a potent inhibitor of MINK1 kinase activity. The patient derived TPF resistant cell based xenograft data suggest that lestaurtinib restores 5FU sensitivity and facilitates a significant reduction of tumor burden. Overall, our study suggests that MINK1 is a major driver of 5FU resistance in OSCC. The novel combination of MINK1 inhibitor lestaurtinib and 5FU needs further clinical investigation in advanced OSCC.