A dopamine-dependent decrease in dorsomedial striatum direct pathway neuronal activity is required for learned motor coordination It has been suggested that the dorsomedial striatum (DMS) facilitates the early stages of motor learning for goal-directed actions, whereas at later stages, control is transferred to the dorsolateral striatum (DLS), which enables the motor actions to become a skill or habit. It is unknown whether these striatal regions are simultaneously active while expertise is acquired during skill learning. To address this question, we developed the "treadmill training task" to track changes in mouse locomotor coordination during practice at running that simultaneously provides a means to measure local neuronal activity using photometry. We analyzed body position and paw movement to evaluate changes in motor coordination over practice sessions on the treadmill using DeepLabCut and custom-built code. By correlating improvements in motor coordination during training with simultaneous neuronal calcium activity in the striatum, we found that DMS direct pathway neurons exhibited decreased activity as the mouse gains proficiency at running. In contrast, direct pathway activity in the DLS was similar throughout training and did not correlate with learned skill proficiency. Pharmacological blockade of D1 dopamine receptors in these subregions during task performance confirmed that dopamine neurotransmission in the DMS direct pathway activity is necessary for efficient motor coordination learning, while dopamine signalling in the DLS is important for both coordination learning and maintenance of the acquired skill. These results provide new tools to measure changes in fine motor skills during simultaneous recordings of brain activity, revealing fundamental features of the neuronal substrates of motor learning.
DeCalciOn: A hardware system for real-time decoding of in-vivo calcium imaging data Epifluorescence miniature microscopes ("miniscopes") are widely used for in vivo calcium imaging of neural population activity. Imaging data is usually collected while subjects are engaged in a task and stored for later offline analysis, but emerging techniques for online imaging offer potential for novel real-time experiments in which closed-loop interventions (such as neurostimulation or sensory feedback) are triggered at short latencies in response to neural population activity. Here we introduce DeCalciOn, a plug-and-play hardware device for online population decoding of in vivo calcium signals that can trigger closed-loop feedback at millisecond latencies, and is compatible with miniscopes that use the UCLA Data Acquisition (DAQ) interface. In performance tests, the position of rats (n=13) on a linear track was decoded in real time from hippocampal CA1 population activity by 24 linear classifiers. DeCalciOn required <2.5 ms after each end-of-frame to decode up to 1,024 calcium traces and trigger TTL control outputs. Decoding was most efficient using a  contour-free method of extracting traces from ROIs that were unaligned with neurons in the image, but  contour-based extraction from neuronal ROIs is also supported. DeCalciOn is an easy-to-use system for real-time decoding of calcium fluorescence that enables closed-loop feedback experiments in behaving animals.
Conformist social learning inoculates against adverse risk aversion Given the ubiquity of potentially adverse behavioural bias owing to myopic trial-and-error learning, it seems paradoxical that improvements in decision-making performance through conformist social learning, a process widely considered to be bias amplification, still prevail in animal collective behaviour. Here we show, through model analyses and large-scale interactive behavioural experiments with 585 human subjects, that conformist influence can indeed promote favourable risk taking in repeated experience-based decision making, even though many individuals are systematically biased towards adverse risk aversion. Although strong positive feedback conferred by copying the majoritys behaviour could result in unfavourable informational cascades, our differential equation model of collective behavioural dynamics identified a key role for increasing exploration by negative feedback arising when a weak minority influence undermines the inherent behavioural bias. This  collective behavioural rescue, emerging through coordination of positive and negative feedback, highlights a benefit of collective learning in a broader range of environmental conditions than previously assumed and resolves the ostensible paradox of adaptive collective behavioural flexibility under conformist influences. Universality of evolutionary trajectories under arbitrary competition dynamics The assumption of constant population size is central in population genetics. It led to a large body of results, that are robust to modeling choices and that have proven successful to understand evolutionary dynamics. In reality, allele frequencies and population size are both determined by the interaction between a population and the environment. Relaxing the constant-population assumption have two big drawbacks. It increases the technical difficulty of the analysis, and it requires specifying a mechanism for the saturation of the population size, possibly making the results contingent on model details. Here, we develop a framework that encompasses a great variety of systems with an arbitrary mechanism for population growth limitation. By using techniques based on scale separation for stochastic processes, we are able to calculate analytically properties of evolutionary trajectories, such as the fixation probability. Remarkably, these properties assume a universal form with respect to our framework, which depends on only three parameters related to the inter-generation timescale, the invasion fitness, and the carrying capacity of the strains. In other words, different systems, such as Lotka-Volterra or a chemostat model (contained in our framework), share the same evolutionary outcomes after a proper re-mapping of their parameters. An important and surprising consequence of our results is that the direction of selection can be inverted, with a population evolving to reach lower values of invasion fitness. Comprehensive-immunopeptidomics analysis reveals presentation of potential neoantigens carrying cancer driver mutations Understanding the properties of human leukocyte antigen (HLA) peptides (HLAp) is essential for cancer precision medicine, while the direct identification of HLAp from a tiny piece of clinical tissues by mass spectrometry (MS) is still confronted with technical challenges. Here, to overcome that hindrances, high-field asymmetric waveform ion mobility spectrometry (FAIMS) is introduced to conduct differential ion mobility (DIM)-MS by seamless gas-phase fractionation optimal for scarce samples. The global-immunopeptidomics analysis enables in-depth identification of 44,785 HLAp, including 2 neoantigens with KRAS-G12V and CPPED1-R228Q from stage IV colorectal cancer (CRC) tissues obtained from 17 subjects. Comparison of tissue-based personal immunopeptidome disclosed cancer-specific processing of HLAp. Since the direct identification of neoantigens from CRC tissues suggests more potential neoantigens have yet to been identified, the targeted-immunopeptidomics by screening-oriented parallel reaction monitoring (PRM) has been established and identified two more neoantigens with oncogenic KRAS. The targeted-immunopeptidomics analysis clarifies the presentation of aimed neoantigens, including neoantigens that are unpredictable by computational algorithms. The comprehensive-immunopeptidomics approach combining the FAIMS-assisted DIM-MS and the screening-oriented PRM contributes the direct and effective identification of HLAp and emerging neoantigens from clinical tissue samples, as an advanced strategy to identify the authentic and recurrent neoantigens for cancer immunotherapy. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. Synthetic reconstruction of the hunchback promoter specifies the role of Bicoid, Zelda and Hunchback in the dynamics of its transcription For over 40 years, the Bicoid-hunchback (Bcd-hb) system in the fruit fly embryo has been used as the model to study how positional information in morphogen concentration gradients is robustly translated into step-like responses. A body of quantitative comparisons between theory and experiment have since questioned the initial paradigm that the sharp hb transcription pattern emerges solely from diffusive biochemical interactions between the Bicoid transcription factor and the gene promoter region. Several alternative mechanisms have been proposed, such as additional sources of positional information, positive feedback from Hb proteins or out-of-equilibrium transcription activation.

By using the MS2-MCP RNA-tagging system and analysing in real time, the transcription dynamics of synthetic reporters for Bicoid and/or its two partners Zelda and Hunchback, we show that all the hb expression pattern features and temporal dynamics are compatible with an equilibrium model with a short decay length Bicoid activity gradient as a sole source of positional information. Meanwhile, Bicoids partners speed-up the process by different means: Zelda lowers the Bicoid concentration threshold required for transcriptional activation while Hunchback reduces burstiness and increases the polymerase firing rate. Binding of single-mutant epidermal growth factor (EGF) ligands alter the stability of the EGF receptor dimer and promote growth signaling The Epidermal Growth Factor Receptor (EGFR) is a membrane-anchored tyrosine kinase that is able to selectively respond to multiple extra-cellular stimuli. Previous studies have indicated that the modularity of this system is affected by ligand-induced differences in the stability of the dimerized receptor in a process known as "Biased signaling". However, this hypothesis has not been explored using single-mutant ligands thus far. Herein, we developed a new approach to identify residues responsible for functional divergence combining the conservation and co-evolution information of ortholog and paralog genes encoding the epidermal growth factor (EGF) ligand. Then, we mutated these residues and assessed the mutants effects on the receptor by employing a combination of molecular dynamics (MD) and biochemical techniques. Although the EGF mutants had comparable binding affinities to the wild type ligand for EGFR, the EGF mutants induced a different phosphorylation and cell growth pattern in multiple cell lines. The MD simulations of the EGF mutants show a long-range effect on the receptor dimer interface. For the first time in this study, a single mutation in EGF is shown to be enough to alter the activation of the pathway at the cellular level. These results also support the theory of biased signaling in the tyrosine kinase receptor system and demonstrate a promising new way to study ligand-receptor interactions. Extent of Molecular Chaperone Association Might Determine Fates of Membraneless Organelles during Aging in C. elegans Proteome imbalance can lead to protein misfolding and aggregation which is associated with pathologies. Protein aggregation can also be an active, organized process and can be exploited by cells as a survival strategy. In adverse conditions, it is beneficial to deposit the proteins in a condensate rather degrading and resynthesizing. Membraneless organelles (MLOs) are biological condensates formed through liquid-liquid phase separation (LLPS), involving cellular components such as nucleic acids and proteins. LLPS is a regulated process, which when perturbed, can undergo a transition from a physiological liquid condensate to pathological solid-like protein aggregates.

To understand how the MLO-associated proteins (MLO-APs) behave during aging, we performed a comparative meta-analysis with age related proteome of C. elegans. We found that the MLO-APs are highly abundant throughout the lifespan. Interestingly, they are aggregating more in long-lived mutant worms compared to the age matched wildtype worms. GO term analysis revealed that the cell cycle and embryonic development are among the top enriched processes in addition to RNP components in insoluble proteome. Considering antagonistic pleotropic nature of these developmental genes and post mitotic status of C. elegans, we assume that these proteins phase transit during post development. As the organism ages, these MLO-APs either mature to become more insoluble or dissolve in uncontrolled manner. However, in the long-lived daf-2 mutant worms, the MLOs may attain protective states due to extended availability and association of molecular chaperones. N-VEGF, the forgotten isoform of VEGF-A, activate cellular stress survival circuits Vascular endothelial growth factor A (VEGF-A) is a secreted protein that stimulates angiogenesis in response to hypoxia. Under hypoxic conditions, a non-canonical long isoform called L-VEGF is concomitantly expressed with VEGF-A. Once translated, L-VEGF it is proteolytically cleaved to generate N-VEGF and VEGF-A. Interestingly, while VEGF-A is secreted and affects the surrounding cells, N-VEGF is mobilized to the nucleus. This suggests that N-VEGF participates in transcriptional response to hypoxia. In this study, we performed a series of complementary experiments to examine the functional role of N-VEGF. Strikingly, we found that the mere expression of N-VEGF followed by its hypoxia-independent mobilization to the nucleus was sufficient to induce key genes associated with angiogenesis, such as Hif1, VEGF-A isoforms, as well as genes associated with cell survival under hypoxia. Complementarily, when N-VEGF was genetically depleted, key hypoxia-induced genes were downregulated and cells were significantly susceptible to hypoxia-mediated apoptosis. This is the first reports of N-VEGF serving as an autoregulatory arm of VEGF-A. Further experiments will be needed to determine the role of N-VEGF in cancer and embryogenesis. The quantitative genetics of the prevalence of infectious diseases: hidden genetic variation due to Indirect Genetic Effects dominates heritable variation and response to selection Infectious diseases have profound effects on life, both in nature and agriculture. Despite the availability of well-established epidemiological theory, however, a quantitative genetic theory of the host population for the endemic prevalence of infectious diseases is almost entirely lacking. While several studies have demonstrated the relevance of the transmission dynamics of infectious diseases for heritable variation and response to selection of the host population, our current theoretical framework of quantitative genetics excludes these dynamics. As a consequence, we do not know which genetic effects of the host population determine the prevalence of an infection, and have no concepts of breeding value and heritable variation for endemic prevalence.

Here we integrate quantitative genetics and epidemiology, and propose a quantitative genetic theory for R0 and for the endemic prevalence of an infectious disease. We first identify the genetic factors that determine the prevalence of an infection, using an approach founded in epidemiological theory. Subsequently we investigate the population level consequences of individual genetic variation, both for R0 and for the endemic prevalence. Next, we present expressions for the breeding value and heritable variation, for both prevalence and individual binary disease status, and show that these parameters depend strongly on the level of the prevalence. Results show that heritable variation for endemic prevalence is substantially greater than currently believed, and increases when prevalence approaches zero, while heritability of individual disease status goes to zero. As a consequence, response of prevalence to selection accelerates considerably when prevalence goes down, in contrast to predictions from classical theory. Finally, we show that most of the heritable variation for the endemic prevalence of an infection is hidden due to indirect genetic effects, suggesting a key role for kin-group selection both in the evolutionary history of current populations and for genetic improvement strategies in animals and plants. A novel Boolean network inference strategy to model early hematopoiesis aging Hematopoietic stem cell (HSC) aging is a multifactorial event that leads to changes in HSC properties and function. These changes are intrinsically coordinated and affect the early hematopoiesis, involving hematopoietic stem and progenitor cells (HSPCs). The objective of this work is to better understand the mechanisms and factors controlling these changes. We have therefore developed an original strategy to construct a Boolean network of genes explaining the priming and homeostasis of HSCs (graphical abstract). Based on our previous scRNA-seq data, we performed an exhaustive analysis of the transcriptional network and identified active transcription modules or regulons along the differentiation trajectory of selected HSPC states. This global view of transcriptional regulation led us to focus on 15 components, 13 selected TFs (Tal1, Fli1, Gata2, Gata1, Zfpm1, Egr1, Junb, Ikzf1, Myc, Cebpa, Bclaf1, Klf1, Spi1) and 2 complexes regulating the ability of HSC to cycle (CDK4/6 - Cyclin D and CIP/KIP). We then defined the connections controlling the differentiation dynamics of HSC states and constructed an influence graph between the TFs involved in the dynamics by mixing observations from our scRNA-seq data and knowledge from the literature. Then, using answer set programming (ASP) and in silico perturbation analysis, we obtained a Boolean model which is the solution of a Boolean satisfiability problem. Finally, perturbation of the model based on age-related changes revealed important regulations, such as the overactivation of Egr1 and Junb or the loss of Cebpa activation by Gata2, which were found to be relevant for the myeloid bias of aged HSC. Our work shows the efficiency of the combination of manual and systematic methods to elaborate a Boolean model. The developed strategy led to the proposal of new regulatory mechanisms underlying the differentiation bias of aged HSCs, explaining the decreased transcriptional priming of HSCs to all mature cell types except megakaryocytes. High expression of VRT2 increases the number of rudimentary basal spikelets in wheat Spikelets are the fundamental building blocks of Poaceae inflorescences and their development and branching patterns determine the various inflorescence architectures and grain yield of grasses. In wheat, the central spikelets produce the most and largest grains, while spikelet size gradually decreases acro- and basipetally, giving rise to the characteristic lanceolate shape of wheat spikes. The acropetal gradient correlates with the developmental age of spikelets, however the basal spikelets are developed first and the cause of their small size and rudimentary development is unclear. Here, we adapted G&T-seq, a low-input transcriptomics approach, to characterise gene expression profiles within spatial sections of individual spikes before and after the establishment of the lanceolate shape. We observed larger differences in gene expression profiles between the apical, central and basal sections of a single spike than between any section belonging to consecutive developmental timepoints. We found that SVP MADS-box transcription factors, including VRT-A2, are expressed highest in the basal section of the wheat spike and display the opposite expression gradient to flowering E-class SEP1 genes. Based on multi-year field trials and transgenic lines, we show that higher expression of VRT-A2 in the basal sections of the spike is associated with increased numbers of rudimentary basal spikelets. Our results, supported by computational modelling, suggest that the delayed transition of basal spikelets from vegetative to floral developmental programmes results in the lanceolate shape of wheat spikes. This study highlights the value of spatially resolved transcriptomics to gain new insights into developmental genetics pathways of grass inflorescences.

One sentence summaryLarge transcriptional gradients exist within a wheat spike and are associated with rudimentary basal spikelet development, resulting in the characteristic lanceolate shape of wheat spikes. Fine-tuning the Spike: Role of the nature and topology of the glycan shield in the structure and dynamics of SARS-CoV-2 S The dense glycan shield is an essential feature of the SARS-CoV-2 spike (S) architecture, key to immune evasion and to the activation of the prefusion conformation. Recent studies indicate that the occupancy and structures of the SARS-CoV-2 S glycans depend not only on the nature of the host cell, but also on the structural stability of the trimer; a point that raises important questions about the relative competence of different glycoforms. Moreover, the functional role of the glycan shield in the SARS-CoV-2 pathogenesis suggests that the evolution of the sites of glycosylation is potentially intertwined with the evolution of the protein sequence to affect optimal activity. Our results from multi-microsecond molecular dynamics simulations indicate that the type of glycosylation at N234, N165 and N343 greatly affects the stability of the receptor binding domain (RBD) open conformation, and thus its exposure and accessibility. Furthermore, our results suggest that the loss of glycosylation at N370, a newly acquired modification in the SARS-CoV-2 S glycan shields topology, may have contributed to increase the SARS-CoV-2 infectivity as we find that N-glycosylation at N370 stabilizes the closed RBD conformation by binding a specific cleft on the RBD surface. We discuss how the absence of the N370 glycan in the SARS-CoV-2 S frees the RBD glycan binding cleft, which becomes available to bind cell-surface glycans, potentially increases host cell surface localization. Visualization of lake succession as indicated by macrophyte species composition and trophic levels in lakes subdivided by volcanic eruption at Akan Caldera, Hokkaido, Japan It is not feasible to continuously observe ecological succession in lakes because of the long time-scales generally involved. Thus, the process has been inductively deduced by comparing many lakes with different succession states, or indirectly simulated by tracking studies of smaller water bodies, experiments using microcosms or mesocosms, and reconstruction of lake history by sediment analysis. However, the reality of succession processes in large lakes with slow succession is not well understood, and new approaches are needed. Theoretically, in a group of large and small lakes of similar ages and with similar initial and watershed environments, the rate of nutrient accumulation in each lake depends on the ratio of watershed area to lake size, and the lakes are predicted to evolve to different trophic levels over time. Here, we tested this hypothesis on the 10 lakes of varying sizes in Akan Caldera, Japan, which were formed thousands of years ago by fragmentation due to volcanic eruptions within the caldera. Topographic and water quality studies showed that the ratio of accumulated watershed area to lake size (area and volume), expressed logarithmically, had a positive linear regression with the total phosphorus concentration, an indicator of trophic level. The trophic levels of the lakes were diverse, including oligotrophic, mesotrophic, and eutrophic types in the traditional "lake type" classification based on total phosphorus and chlorophyll-a concentrations. Furthermore, 21 species of aquatic macrophytes were observed by a diving survey, and the plant species composition was classified into five groups corresponding to the trophic status of the lakes, indicating a conventional "hydrarch succession". The diversity of water quality and aquatic vegetation in a group of lakes with similar origins paves the way for new comparative studies of lakes, including large lakes. 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. LPS-TLR4 pathway exaggerates alcohol-induced acute-on-chronic liver failure via provoking NETs formation Backgrounds: Intrahepatic infiltration of neutrophils is a character of alcoholic acute-on-chronic liver failure (AACLF) and neutrophil extracellular traps (NETs) are an important strategy for neutrophils to fix and kill invading microorganisms. Intestinal bacteria and the gut-liver axis have been thought to play a key role in many liver diseases also including AACLF. However, whether NETs appear in AACLF and play a role in AACLF is still unsure. Methods: WT, NE KO, and TLR4 KO mice were used to build the AACLF model, and the intestinal bacteria were eliminated at the same time and LPS was given. Then the formation of NETs and AACLF related markers were detected. Results: The serum MPO-DNA and LPS concentration was increased in AACLF patients and a correlation was revealed between these two indexes. More intrahepatic NETs formed in AACLF mice by testing MPO-DNA, Cit H3, and NE. These markers decreased with gut detergent and restored markers with gut detergent plus LPS supplement. While NETs formation failed to change with gut microbiome or combine LPS supplement in TLR4 KO mice. As we tested AACLF related characters, liver injury, intrahepatic fat deposition, inflammation, and fibrosis alleviated with depletion of NE. These related marks were also attenuated with gut sterilization by antibiotics and recovered with combined treatment with antibiotics plus LPS. But the liver injury, intrahepatic fat, fibro deposition, and liver inflammation-related markers did show a significant difference in TLR4 KO mice when they received the same treatment. Conclusion: Intestinal-derived LPS promotes NETs formation in AACLF through the TLR4 pathway and further accelerates the AACLF process by NETs. Cell-specific Bioorthogonal Tagging of Glycoproteins Altered glycosylation is an undisputed corollary of cancer development. Understanding these alterations is paramount but hampered by limitations underlying cellular model systems. For instance, the intricate interactions between tumour and host cannot be adequately recapitulated in monoculture of tumour-derived cell lines. More complex co-culture models usually rely on sorting procedures for proteome analyses and rarely capture the details of protein glycosylation. Here, we report a strategy termed Bio-Orthogonal Cell line-specific Tagging of Glycoproteins (BOCTAG). Cells are equipped by transfection with an artificial biosynthetic pathway that transforms bioorthogonally tagged sugars into the corresponding nucleotide-sugars. Only transfected cells incorporate bioorthogonal tags into glycoproteins in the presence of non-transfected cells. We employ BOCTAG as an imaging technique and to annotate cell-specific glycosylation sites in mass spectrometry-glycoproteomics. We demonstrate application in co-culture and mouse models, allowing for profiling of the glycoproteome as an important modulator of cellular function. Trypanosomal variant surface glycoprotein expression in human African trypanosomiasis patients Trypanosoma brucei gambiense, an extracellular protozoan parasite, is the primary causative agent of human African trypanosomiasis. T. b. gambiense is endemic to West and Central Africa, where it is transmitted by the bite of infected tsetse flies. In the bloodstream of an infected host, the parasite evades antibody recognition by altering the Variant Surface Glycoprotein (VSG) that forms a dense coat on its cell surface through a process known as antigenic variation. Each VSG has a variable N-terminal domain exposed to the host and a less variable C-terminal domain that is partially hidden from host antibodies. Our lab developed VSG-seq, a targeted RNA-seq method, to study VSG gene expression in T. brucei. Studies using VSG-seq to characterize antigenic variation in a mouse model have revealed marked diversity in VSG expression within parasite populations, but this finding has not yet been validated in a natural human infection. Here, we used VSG-seq to analyze VSGs expressed in the blood of 12 patients infected with T. b. gambiense. The number of VSGs identified per patient ranged from one to 14; notably, two VSGs were shared by more than one patient. Analysis of expressed VSG N-terminal domain types revealed that 82% of expressed VSGs encoded a type B N-terminus, a bias not seen in datasets from other T. brucei subspecies. C-terminal types in T. b. gambiense infection were also restricted. These results demonstrate a bias either in the underlying VSG repertoire of T. b. gambiense or in this subspecies expression of VSGs during infection. This work demonstrates the feasibility of using VSG-seq to study antigenic variation in human infections and highlights the importance of understanding VSG repertoires in the field.

Author SummaryHuman African Trypanosomiasis is a neglected tropical disease primarily caused by the extracellular parasite Trypanosoma brucei gambiense. To avoid elimination by the host, these parasites repeatedly replace their dense surface coat of Variant Surface Glycoprotein (VSG). Despite the important role of VSGs in prolonging infection, VSG expression during natural human infections is poorly understood. A better understanding of natural VSG gene expression dynamics can clarify the mechanisms that T. brucei uses to alter its VSG coat and improve how trypanosomiasis is diagnosed in humans. We analyzed the expressed VSGs detected in the blood of patients with trypanosomiasis. Our findings indicate that parasites in both natural and experimental infections express diverse sets of VSGs. Comparison study of sixteen differential abundance testing methods using two large Parkinson disease gut microbiome datasets BackgroundWhen studying the relationship between the microbiome and a disease, a common question asked is what individual microbes are differentially abundant between a disease and healthy state. Numerous differential abundance (DA) testing methods exist and range from standard statistical tests to methods specifically designed for microbiome data. Comparison studies of DA testing methods have been performed, but none were performed on microbiome datasets collected for the study of real, complex disease. Due to this, we performed DA testing of microbial genera using 16 DA methods in two large, uniformly collected gut microbiome datasets on Parkinson disease (PD), and compared their results.

ResultsPairwise concordances between methods ranged from 46%-99% similarity. Average pairwise concordance per dataset was 76%, and dropped to 62% when taking replication of signals across datasets into account. Certain methods consistently resulted in above average concordances (e.g. Kruskal-Wallis, ALDEx2, GLM with centered-log-ratio transform), while others consistently resulted in lower than average concordances (e.g. edgeR, fitZIG). Overall, [~]80% of genera tested were detected as differentially abundant by at least one method in each dataset. Requiring associations to replicate across datasets reduced significant signals by almost half. Further requirement of signals to be replicated by the majority of methods ([&ge;]8) yielded 19 associations. Only one genus (Agathobacter) was replicated by all methods. Use of hierarchical clustering revealed three groups of DA signatures that were (1) replicated by the majority of methods and included genera previously associated with PD, (2) replicated by few or no methods, and (3) replicated by a subset of methods and included rarer genera, all enriched in PD.

ConclusionsDifferential abundance tests yielded varied results. Using one method on one dataset may find true associations, but may also detect non-reproducible signals, adding to inconsistency in the literature. To help lower false positives, one might analyze data with two or more DA methods to gauge concordance, and use a built-in replication dataset to show reproducibility. This study corroborated previously reported microorganism associations in PD, and revealed a potential new group of microorganisms whose abundance is significantly elevated in PD, and might be worth pursuing in future investigations. 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. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. 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. Soluble ECM promotes organotypic formation in lung alveolar model Micropatterned suspension culture creates consistently sized and shaped cell aggregates but has not produced organotypic structures from stable cells, thus restricting its use in accurate disease modeling. Here, we show that organotypic structure is achieved in hybrid suspension culture via supplementation of soluble extracellular matrix (ECM). We created a viable lung organoid from epithelial, endothelial, and fibroblast human stable cell lines in suspension culture. We demonstrate the importance of soluble ECM in organotypic patterning with the emergence of lumen-like structures with airspace showing feasible gas exchange units, formation of branching, perfusable vasculature, and long-term 70-day maintenance of lumen structure. Our results show a dependent relationship between enhanced fibronectin fibril assembly and the incorporation of ECM in the organoid. We successfully applied this technology in modeling lung fibrosis via bleomycin induction and test a potential antifibrotic drug in vitro while maintaining fundamental cell-cell interactions in lung tissue. Our human fluorescent lung organoid (hFLO) model represents features of pulmonary fibrosis which were ameliorated by fasudil treatment. We also demonstrate a 3D culture method with potential of creating organoids from mature cells, thus opening avenues for disease modeling and regenerative medicine, enhancing understanding of lung cell biology in health and lung disease. TCAB1 prevents nucleolar accumulation of the telomerase RNA to promote telomerase assembly Localization of a wide variety of RNAs to non-membrane bound cellular compartments such as nucleoli, Cajal bodies, and stress-granules is critical for their function and stability. The molecular mechanisms that underly the recruitment and exclusion of specific RNAs from these phase-separated organelles is poorly understood. Telomerase is a ribonucleoprotein (RNP), that is composed of the reverse transcriptase protein TERT, the telomerase RNA (TR), and several auxiliary proteins that associate with TR, including TCAB1. Here we show that, that in the absence of TCAB1, TR is sequestered in the nucleolus, while TERT localizes to the nucleoplasm and is excluded from the nucleolus, which prevents telomerase assembly. Thus, nuclear compartmentalization by the non-membrane bound nucleolus counteracts telomerase assembly and TCAB1 is required to exclude the telomerase RNA from the nucleolus. Our work provides general insight into the mechanism and functional consequences of RNA recruitment to organelles formed by phase-separation and proposes a new model explaining the critical role of TCAB1 in telomerase function. Blocking antibodies against integrin-a3, integrin-aM, and integrin-aMb2 de-differentiate myofibroblasts and reverse lung and kidney fibroses in a mouse model Fibrosis is involved in 45% of deaths in the United States, and no treatment exists to reverse the progression of the disease. Myofibroblasts are key to the progression and maintenance of fibrosis. We investigated features of cell adhesion necessary for monocytes to differentiate into myofibroblasts, seeking to identify pathways key to myofibroblast differentiation. Blocking antibodies against integrins 3, M, and M{beta}2 de-differentiate myofibroblasts in vitro, lower the pro-fibrotic secretome of myofibroblasts, and reverse lung and kidney fibrosis in vivo. Decorins collagen-binding peptide directs blocking antibodies (against integrins-3, -M, -M{beta}2) to both fibrotic lungs and fibrotic kidneys, reducing the dose of antibody necessary to reverse fibrosis. This targeted immunotherapy blocking key integrins may be an effective therapeutic for the treatment and reversal of fibrosis.

SummaryBlocking antibodies against integrins-3, -M, and -M{beta}2 can be targeted to sites of fibrosis, reverse lung and kidney fibroses, and offer the potential to bring immunotherapy to fibrosis