Robust effects of corticothalamic feedback and behavioral state on movie responses in mouse dLGN Neurons in the dorsolateral geniculate nucleus (dLGN) of the thalamus receive a substantial proportion of modulatory inputs from corticothalamic (CT) feedback and brain stem nuclei. Hypothesizing that these modulatory influences might be differentially engaged depending on the visual stimulus and behavioral state, we performed in vivo extracellular recordings from mouse dLGN while optogenetically suppressing CT feedback and monitoring behavioral state by locomotion and pupil dilation. For naturalistic movie clips, we found CT feedback to consistently increase dLGN response gain and promote tonic firing. In contrast, for gratings, CT feedback effects on firing rates were mixed. For both stimulus types, the neural signatures of CT feedback closely resembled those of behavioral state, yet effects of behavioral state on responses to movies persisted even when CT feedback was suppressed. We conclude that CT feedback modulates visual information on its way to cortex in a stimulus-dependent manner, but largely independently of behavioral state.
Sex differences in behavioral responding and dopamine release during Pavlovian learning Learning associations between cues and rewards requires the mesolimbic dopamine system. The dopamine response to cues signals differences in reward value in well-trained animals. However, these value-related dopamine responses are absent during early training sessions when cues signal differences in the reward rate. These findings suggest cue-evoked dopamine release conveys differences between outcomes only after extensive training, though it is unclear if this is unique to when cues signal differences in reward rate, or if this is also evident when cues signal differences in other value-related parameters such as reward size. To address this, we utilized a Pavlovian conditioning task in which one audio cue was associated with a small reward (one pellet) and another audio cue was associated with a large reward (three pellets). We performed fast-scan cyclic voltammetry to record changes in dopamine release in the nucleus accumbens of male and female rats throughout learning. While female rats exhibited higher levels of conditioned responding, a faster latency to respond, and elevated post-reward head entries relative to male rats, there were no sex differences in the dopamine response to cues. Multiple training sessions were required before cue-evoked dopamine release signaled differences in reward size. Reward-evoked dopamine release scaled with reward size, though females displayed lower reward-evoked dopamine responses relative to males. Conditioned responding related to the decrease in the peak reward-evoked dopamine response and not to cue-evoked dopamine release. Collectively these data illustrate sex differences in behavioral responding as well as in reward-evoked dopamine release during Pavlovian learning.
Is Vaccination a Viable Method to Control Johne's Disease Caused by Mycobacterium avium subsp. paratuberculosis? Data from 12 Million Ovine Vaccinations and 7.6 Million Carcass Examinations in New South Wales, Australia from 1999-2009 BackgroundMycobacterium avium subsp. paratuberculosis (MAP) causes Johnes disease (or paratuberculosis), a chronic wasting disease of ruminants and other animals resulting from granulomatous enteritis. There are increasing concerns that MAP is zoonotic. The prevalence of Johnes disease is increasing worldwide. In an attempt to control an epidemic of ovine Johnes disease (OJD) in New South Wales (NSW), a government/industry sponsored voluntary vaccination/on-farm management program commenced in 2000. We report herein an observational study of changes in disease prevalence as vaccination progressed, based on abattoir surveillance data for OJD from 1999 to 2009. We also discuss the epidemiological, policy, regulatory, research, economic and sociological elements that contributed to the development of a mature control program, whose aim was to halt the epidemic spread of OJD in a naive sheep population.

MethodsNSW was divided into areas of "High" (HPA), "Medium" (MPA) and "Low" (LPA) OJD prevalence. A killed whole cell vaccine (Gudair(R)) was administered to sheep from 2000 to 2009. Trained examiners evaluated the viscera of adult sheep carcasses at slaughter for gross evidence of OJD. MAP infection was confirmed by histopathology.

Principal findingsFrom 2000-2009, 12 million vaccine doses were administered in NSW (91%; 10.9 million in the HPA). Many of the vaccinated flocks were suffering > 5% annual mortality in adult sheep, with some individual flocks with 10-15% losses attributable to OJD. A total of 7.6 million carcasses were examined (38%; 2.9 million from the HPA). Overall, 16% of slaughter consignments (sheep consigned to the abattoir from a single vendor) were positive for OJD, of which 94% were from the HPA. In the HPA, the percentage of animals with lesions attributable to OJD at slaughter fell progressively from 2.4% (10,406/432,860) at commencement of vaccination in 2000 to 0.8% (1,573/189,564) by 2009. Herd immunity from vaccination in the HPA was estimated at 70% by 2009. This coincided with a progressive decrease in reports of clinical disease and mortalities in vaccinated flocks.

SignificanceWe show a decrease in the prevalence of lesions attributable to OJD in NSW concomitant with initiation of voluntary vaccination, on-farm management plans, abattoir monitoring and feedback of animal prevalence data to sheep producers. We conclude that a target of [&le;] 1% regional prevalence of OJD affected sheep at slaughter is achievable using these interventions. Cascade screening following a polygenic risk score test: what is the risk of a relative conditional on a high score of a proband? Polygenic risk scores (PRSs) for predicting disease risk have become increasingly accurate, leading to rising popularity of PRS tests. Consider an individual whose PRS test has placed him/her at the top q-quantile of genetic risk. Recently, Reid et al. (Circ Genom Precis Med. 2021;14:e003262) have investigated whether such a finding should motivate cascade screening in the probands siblings. Specifically, using data from the UK biobank, Reid et al. computed the empirical probability of a sibling of the proband to also have a PRS at the top q-quantile. In this short note, I use the liability threshold model to compute this probability analytically (for either a sibling of the proband or for a more distant relative), showing excellent agreement with the empirical results of Reid et al., including that this probability is disease-independent. Further, I compute the probability of the relative of the proband to be affected, as a function of the quantile threshold q, the proportion of variance explained by the score, and the disease prevalence. Honeybee cognition as a tool for scientific engagement In addition to the crowdsourcing of experimental data, citizen science, and scientific engagement more broadly, serve as a bridge between researchers and the wider community. This serves to foster a greater understanding of the scientific method and science-based solutions generally. Apis mellifera (honeybees) are a well-established model for the study of learning and cognition and can serve as an engaging outreach system for this wider community. Here, we developed and implemented a protocol using well established honeybee conditioning protocols to safely study the effects of caffeine and dopamine on learning performance. Using this protocol, a group of high-school aged students as part of the Ryukyu Girls program demonstrated that caffeine, but not dopamine, significantly reduced the number of trials required for a successful conditioning response. This allowed these students to explore the scientific method in a relatable and engaging way.

Simple SummaryGlobal scientific literacy can be improved through widespread and effective community engagement by researchers. We propose Apis mellifera (honeybee) as an public engagement tool due to widespread awareness of colony collapse and the bees importance in food production. Moreover, their cognitive abilities make for engaging experiments. Their relative ease of cultivation means that studies can be performed cost-effectively, especially when partnering with local aperists. Using a proxy for honeybee learning, a group of non-specialist high-school-aged participants obtained data suggesting that caffeine, but not dopamine, improved learning. This hands-on experience facilitated student understanding of the scientific method, factors that shape learning and the importance of learning for hive health. Mechanistic basis for understanding the dual activities of the bifunctional Azotobacter vinelandii mannuronan C-5 epimerase and alginate lyase AlgE7 The functional properties of alginates are dictated by the monomer composition and molecular weight distribution. Mannuronan C-5 epimerases determine the former by epimerizing {beta}-D-mannuronic acid residues (M) into -L-guluronic acid residues (G). The molecular weight is affected by alginate lyases, which cleave alginate chains through {beta}-elimination. The reaction mechanisms for the epimerization and cleavage are similar and some enzymes can perform both. These dualistic enzymes share high sequence identity with mannuronan C-5 epimerases without lyase activity, and the mechanism behind their activity as well as the amino acids responsible for it are still unknown. In this study, we investigate mechanistic determinants of the bifunctional epimerase and lyase activity of AlgE7 from Azotobacter vinelandii. Based on sequence analyses, a range of AlgE7 variants were constructed and subjected to activity assays and product characterization by NMR. Our results show that the lyase activity of AlgE7 is regulated by the type of ion present: Calcium promotes it, whereas NaCl reduces it. By using defined poly-M and poly-MG substrates, the preferred cleavage sites of AlgE7 were found to be M{downarrow}XM and G{downarrow}XM, where X can be either M or G. By studying AlgE7 mutants, R148 was identified as an important residue for the lyase activity, and the point mutant R148G resulted in an enzyme with only epimerase activity. Based on the results obtained in the present study we suggest a unified catalytic reaction mechanism for both epimerase and lyase activities where H154 functions as the catalytic base and Y149 as the catalytic acid.

ImportancePost-harvest valorisation and upgrading of algal constituents is a promising strategy in the development of a sustainable bioeconomy based on algal biomass. In this respect, alginate epimerases and lyases are valuable enzymes for tailoring of the functional properties of alginate, a polysaccharide extracted from brown seaweed with numerous applications in food, medicine and material industries. By providing a better understanding of the reaction mechanism and of how the two enzyme reactions can be altered by changes in reaction conditions, this study opens for further applications of bacterial epimerases and lyases in enzymatic tailoring of alginate polymers. Engineered disulfide reveals structural dynamics of locked SARS-CoV-2 spike The spike (S) protein of SARS-CoV-2 has been observed in three distinct pre-fusion conformations: locked, closed and open. Of these, the function of the locked conformation remains poorly understood. Here we engineered a SARS-CoV-2 S protein construct "S-R/x3" to arrest SARS-CoV-2 spikes in the locked conformation by a disulfide bond. Using this construct we determined high-resolution structures confirming that the x3 disulfide bond has the ability to stabilize the otherwise transient locked conformations. Structural analyses reveal that wild-type SARS-CoV-2 spike can adopt two distinct locked-1 and locked-2 conformations. For the D614G spike, based on which all variants of concern were evolved, only the locked-2 conformation was observed. Analysis of the structures suggests that rigidified domain D in the locked conformations interacts with the hinge to domain C and thereby restrains RBD movement. Structural change in domain D correlates with spike conformational change. We propose that the locked-1 and locked-2 conformations of S are present in the acidic high-lipid cellular compartments during virus assembly and egress. In this model, release of the virion into the neutral pH extracellular space would favour transition to the closed or open conformations. The dynamics of this transition can be altered by mutations that modulate domain D structure, as is the case for the D614G mutation, leading to changes in viral fitness. The S-R/x3 construct provides a tool for the further structural and functional characterization of the locked conformations of S, as well as how sequence changes might alter S assembly and regulation of receptor binding domain dynamics. Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent global expansion A novel coronavirus, SARS-CoV-2, has caused over 190 million cases and over 4 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e. M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic one year later. By analyzing 4,013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of initial phase of global ignition over a timespan of 14-week. 11 major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted a half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implying that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261,323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e. including all the publicly available viral genomes) in order to recapitulate our findings over one-year timespan. By 25 December 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the four variants of concern (VOC) are the descendants of WE1 type. This study demonstrates the viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time. SARS-CoV-2 infection induces soluble platelet activation markers and PAI-1 in the early moderate stage of COVID-19 IntroductionCoagulation dysfunction and thromboembolism emerge as strong comorbidity factors in severe COVID-19. However, it is unclear when particularly platelet activation markers and coagulation factors dysregulated during the pathogenesis of COVID-19. Here, we sought to assess the levels of coagulation and platelet activation markers at moderate and severe stages of COVID-19 to understand the pathogenesis.

MethodsTo understand this, hospitalized COVID-19 patients with (severe cases that required intensive care) or without pneumonia (moderate cases) were recruited. Phenotypic and molecular characterizations were performed employing basic coagulation tests including PT, APTT, D-Dimer and TFPI. The flow cytometry-based multiplex assays were performed to assess FXI, anti-thrombin, prothrombin, fibrinogen, FXIII, P-selectin, sCD40L, plasminogen, tissue-plasminogen activator (tPA), plasminogen activator inhibitor-1 (PAI-1) and D-Dimer.

ResultsThe investigations revealed induction of plasma P-selectin and CD40 ligand (sCD40L) in moderate COVID-19 cases which were significantly abolished with the progression of COVID-19 severity. Moreover, a profound reduction in plasma tissue factor pathway inhibitor (TFPI) and FXIII were identified particularly in the severe COVID-19. Further analysis revealed fibrinogen induction in both moderate and severe patients. Interestingly, an elevated PAI-1 more prominently in moderate, and tPA particularly in severe COVID-19 cases were observed. Particularly, the levels of fibrinogen and tPA directly correlated with the severity of the disease.

ConclusionsIn summary, induction of soluble P-selectin, sCD40L, fibrinogen and PAI-1 in moderate COVID-19 cases suggests the activation of platelets and coagulation system before patients require intensive care. These findings would help in designing better thromboprophylaxis to limit the COVID-19 severity. A nontoxic tetrodotoxin analog, 5,6,11-trideoxytetrodotoxin, as an olfactory chemoattractant for the grass puffer Toxic puffers contain the potent neurotoxin, tetrodotoxin (TTX). Although TTX is considered to serve as a defense substance, previous behavioral studies have demonstrated that TTX acts as an attractive pheromone for some toxic puffers. To elucidate the physiological mechanism of putative pheromonal action of TTX, we examined whether grass puffers Takifugu alboplumbeus can detect TTX. Electroolfactogram (EOG) results suggest that the olfactory epithelium (OE) of grass puffers responded to a type of TTX analog (5,6,11-trideoxyTTX), although it did not respond to TTX. We also examined the attractive action of 5,6,11-trideoxyTTX on grass puffers by recording their swimming behavior under dark conditions. Grass puffers preferred to stay on the side of the aquarium where 5,6,11-trideoxyTTX was administered, and their swimming speed decreased. Additionally, odorant-induced labeling of olfactory sensory neurons by immunohistochemistry against neural activity marker (phosphorylated extracellular signal regulated kinase; pERK) revealed that labeled olfactory sensory neurons were localized in the region surrounding "islets" where there was considered as nonsensory epithelium. 5,6,11-trideoxyTTX has been known to accumulate in grass puffers, but its toxicity is much lower (almost nontoxic) than TTX. Our results suggest that toxic puffers may positively use this unknown function TTX analog as an olfactory chemoattractant. The BAF A12T mutation associated with premature aging impedes lamin A/C recruitment to sites of nuclear rupture, contributing to nuclear envelope fragility. The premature aging disorder Nestor Guillermo Progeria Syndrome (NGPS) is caused by a homozygous Alanine to Threonine mutation at position 12 (A12T) in Barrier-to- Autointegration Factor (BAF). BAF is a small essential protein that binds to DNA and nuclear envelope proteins. It contributes to important cellular processes including transcription regulation and nuclear envelope reformation after mitosis. More recently, BAF was identified as an important factor for nuclear envelope repair upon rupture in interphase. However, the mechanism by which the BAF A12T mutation causes NGPS has remained unclear. To investigate the effects of this mutation on nuclear envelope integrity, we used NGPS-derived patient cells and engineered an isogenic cell line by reversing the BAF A12T homozygous mutation using CRISPR/Cas9. Using a combination of cellular models, structural data and in vitro assays, we identified that the A12T mutation reduces the affinity of BAF for lamin A/C by tenfold. As a result, BAF A12T is unable to recruit lamin A/C to sites of nuclear envelope rupture. This leads to persistent lamin A/C gaps at sites of ruptures, and contributes to nuclear fragility in NGPS patient cells, which show increased frequency of nuclear envelope re- rupturing. Overexpression of wild-type BAF in a NGPS context rescues lamin A/C recruitment to sites of nuclear rupture, which could explain why the heterozygous A12T mutation does not cause premature aging. Identification and characterization of the WYL BrxR protein and its gene asseparable regulatory elements of a BREX phage restriction system Bacteriophage exclusion ( BREX) phage restriction systems are found in a wide range of bacteria. Various BREX systems encode unique combinations of proteins that usually include a site-specific methyltransferase; none appear to contain a nuclease. Here we describe the identification and characterization of a Type I BREX system from Acinetobacter and the effect of deleting each BREX ORF on growth, methylation, and restriction. We identified a previously uncharacterized gene in the BREX operon that is dispensable for methylation but involved in restriction. Biochemical and crystallographic analyses of this factor, which we term BrxR ( BREX Regulator), demonstrate that it forms a homodimer and specifically binds a DNA target site upstream of its transcription start site. Deletion of the BrxR gene causes cell toxicity, reduces restriction, and significantly increases the expression of BrxC. In contrast, the introduction of a premature stop codon into the BrxR gene, or a point mutation blocking its DNA binding ability, has little effect on restriction, implying that the BrxR coding sequence and BrxR protein play independent functional roles. We speculate that the BrxR coding sequence is involved in cis regulation of anti-phage activity, while the BrxR protein plays an additional regulatory role, perhaps during horizontal transfer. Coupled myovascular expansion directs cardiac growth and regeneration BackgroundHeart regeneration requires multiple cell types to enable cardiomyocyte (CM) proliferation. How these cells interact to create growth niches is unclear. Here we profile proliferation kinetics of cardiac endothelial cells (CECs) and CMs in the neonatal mouse heart and find that they are spatiotemporally coupled.

Methods and ResultsWe show that coupled myovascular expansion during cardiac growth or regeneration is dependent upon VEGF-VEGFR2 signaling, as genetic deletion of Vegfr2 from CECs or inhibition of VEGFA abrogates both CEC and CM proliferation. Repair of cryoinjury displays poor spatial coupling of CEC and CM proliferation. Boosting CEC density after cryoinjury with virus encoding Vegfa enhances regeneration. Using Mendelian randomization, we demonstrate that circulating VEGFA levels are positively linked with human myocardial mass, suggesting that Vegfa can stimulate human cardiac growth.

ConclusionsOur work demonstrates the importance of coupled CEC and CM expansion and reveals a myovascular niche that may be therapeutically targeted for heart regeneration. Anchorage by seed mucilage prevents seed dislodgement in high surface flow: a mechanistic investigation O_LIBackground and Aims: Seed mucilage is a common and highly diverse trait shared among thousands of angiosperm species. While long recognized that mucilage allows seeds to anchor to substrates (antitelechory), resisting abiotic and biotic dislodgement, we still lack a mechanistic understanding of this process.
C_LIO_LIMethods: We propose a mechanistic model of how mucilage affects substrate anchorage and fluid resistance, ultimately contributing to dislodgement resistance. To test this model, we subjected mucilaginous seeds of 52 species, varying in eight measured seed traits, to seven days of continuous water flow at a range of dislodgement potentials.
C_LIO_LIKey Results: Supporting our model, mucilage mass increased force necessary to dislodge both dry and wet seeds; our measurement of the dislodgement force of dry mucilage explained time to dislodgement well. The effect size was remarkably large; increasing the standardized mucilage mass by one standard deviation resulted in a 280-fold increase in the time to dislodgement. Fluid resistance was largely dependent on speed of water flow and the seeds modified drag coefficient, but not seed traits. Neither mucilage expansion speed nor mucilage decay rate explained dislodgement potential well.
C_LIO_LIConclusions: Our results suggest that the degree of anchorage to substrate, measured with a simple dislodgement force assay, is highly predictive of mucilaginous seed retention in highly erosive environments. In contrast, we found that other seed and mucilage traits are of lesser importance to anchorage.
C_LI spatialGE: Quantification and visualization of the tumor microenvironment heterogeneity using spatial transcriptomics SummarySpatially-resolved transcriptomics promises to increase our understanding of the tumor microenvironment and improve cancer prognosis and therapies. Nonetheless, analytical methods to explore associations between the spatial heterogeneity of the tumor and clinical data are not available. Hence, we have developed spatialGE, a software that provides visualizations and quantification of the tumor microenvironment heterogeneity through gene expression surfaces, spatial heterogeneity statistics (SThet) that can be compared against clinical information, spot-level cell deconvolution, and spatially-informed clustering (STclust), all using a new data object to store data and resulting analyses simultaneously.

Availability and implementationThe R package and tutorial/vignette are available at https://github.com/FridleyLab/spatialGE. A script to reproduce the analyses in this manuscript is available in Supplementary information.

ContactFridley.Lab@Moffitt.org or Brooke.Fridley@Moffitt.org

Supplementary informationAvailable at Bioinformatics online.



O_FIG O_LINKSMALLFIG WIDTH=198 HEIGHT=200 SRC="FIGDIR/small/454023v2_figG1.gif" ALT="Figure 1">
View larger version (73K):
org.highwire.dtl.DTLVardef@191cbe7org.highwire.dtl.DTLVardef@e1b6c7org.highwire.dtl.DTLVardef@a6f68forg.highwire.dtl.DTLVardef@1853601_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOGraphical abstractC_FLOATNO Overview of spatialGE features. A. The STList data object from spatialGE can be creared from several sources, including comma- or tab-separated files containing gene counts and spatial coordinates. The object can also be created directly from Visium outputs, Seurat objects, or GeoMx outputs. B. Users can optionally provide a metadata file, containing information associated with each sample (one row per sample, or per ROI if GeoMx data). C. Methods for quality control of data are provided by spatialGE, including visualizations of counts and genes per spot, as well as filtering of spots or genes within user-determined thresholds. D. A novel method (STclust) performs spatially informed clustering of spots and tissue domain identification. E. spatialGE provides different types of data visualization, including gene expression at each spot ("quilt plots"), as well as adaptation of spatial interpolation ("kriging") to spatial transcriptomics data (transcriptomic surface). F. spatialGE also leverages spatial statistics (Morans I, Gearys C, Getis-Ord Gi) to quantitatively describe heterogeneity within the tumor microenvironment and to explore associations between spatial heterogeneity and clinical oucomes. G. Gene expression deconvolution can also be applied to each spot to detect immune cell types (xCell) and classification of spots as tumor or stroma (ESTIMATE).

C_FIG High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Membrane remodeling due to mixture of multiple curvature proteins We present an extension of the Monte Carlo based mesoscopic membrane model, where the membrane is represented as a dynamically triangulated surface and the proteins are modeled as anisotropic inclusions formulated as in-plane nematic field variables adhering to the deformable elastic sheet. In the extended model, we have augmented the Hamiltonian to study membrane deformation due to a mixture of multiple types of curvature generating proteins. This feature opens the door for understanding how multiple kinds of curvature-generating proteins may be working in a coordinated manner to induce desired membrane morphologies. For example, among other things, we study membrane deformations and tubulation due to a mixture of positive and negative curvature proteins as mimics of various proteins from BAR domain family working together for curvature formation and stabilization. We also study the effect of membrane anisotropy, which manifests as membrane localization and differential binding affinity of a given curvature protein, leading to insights into the tightly regulated cargo sorting and transport processes. Our simulation results show different morphology of deformed vesicles that depend on the curvatures and densities of the participating proteins as well as on the protein-protein and membrane-proteins interactions.

SignificanceMembrane remodeling requires highly orchestrated interactions between several types of lipids and curvature proteins. Experimentally probing the membrane deformation intermediates is non-trivial. For example, it is not known exactly how two or more different kinds of BAR domain proteins work in concert to induce and stabilize certain membrane curvatures. In this work, we use mesoscopic continuum modeling and explore the factors that induce and stabilize a range of membrane deformation due to a mixture of positive and negative curvature proteins. Small extracellular vesicles from young mice prevent frailty, improve healthspan and decrease epigenetic age in old mice. Aging is associated with an increased risk of frailty, disability, comorbidities, institutionalization, falls, fractures, hospitalization, and mortality. Searching for strategies to delay the degenerative changes associated with aging and frailty is interesting. We treated old animals intravenously with small extracellular vesicles (sEVs) derived from adipose mesenchymal stem cells (ADSCs) of young animals, and we found an improvement of several functional parameters usually altered with aging, such as motor coordination, grip strength, fatigue resistance, fur regeneration, and renal function. Frailty index analysis showed that 40% of old control mice were frail, whereas none of the old ADSCs-sEVs treated mice were. Molecular and structural benefits in muscle and kidney accompanied this functional improvement. ADSCs-sEVs induced pro-regenerative effects and a decrease in oxidative stress, inflammation, and senescence markers. Moreover, predicted epigenetic age was lower in tissues of old mice treated with ADSCs-sEVs and their metabolome changed to a youth-like pattern. Finally, we gained some insight into the miRNAs contained in sEVs that might be, at least in part, responsible for the effects observed. We propose that young sEVs treatment can be beneficial against frailty and therefore can promote healthy aging. High performance single-cell gene regulatory network inference atscale: The Inferelator 3.0 MotivationGene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above.

ResultsIn this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data.

AvailabilityThe inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/). The limits of evolutionary convergence in sympatry: reproductive interference and historical constraints leading to local diversity in warning traits Mutualistic interactions between defended species represent a striking case of evolutionary convergence in sympatry, driven by the increased protection against predators brought by mimicry in warning traits. However, such convergence is often limited: sympatric defended species frequently display different or imperfectly similar warning traits. The phylogenetic distance between sympatric species may indeed prevent evolution towards the exact same signal. Moreover, warning traits are also involved in mate recognition, so that trait convergence might result in heterospecific courtship and mating. Here, we develop a mathematical model to investigate the strength and direction of evolution of warning trait in defended species with different ancestral traits. Specifically, we determine the effect of phenotypic distances between ancestral trait states of sympatric defended species and of costs of heterospecific sexual interactions on imperfect mimicry and trait divergence. Our analytical results confirm that reproductive interference and historical constraints limit the convergence of warning traits, leading to either complete divergence, or imperfect mimicry. Our model reveals that imperfect mimicry evolves only when ancestral trait values differ between species, due to historical constraints and highlights the importance of female and predator discriminations in the evolution of such imperfect mimicry. Our study thus provides new predictions on how reproductive interference interacts with historical constraints and may promote the emergence of novel warning traits, enhancing mimetic diversity. Comparative toxicogenomics of glyphosate and Roundup herbicides by mammalian stem cell-based genotoxicity assays and molecular profiling in Sprague-Dawley rats Whether glyphosate-based herbicides (GBHs) are more potent than glyphosate alone at activating cellular mechanisms, which drive carcinogenesis remains controversial. As GBHs are more cytotoxic that glyphosate, we reasoned they may also be more capable of activating carcinogenic pathways. We tested this hypothesis by comparing the effects of glyphosate with Roundup GBHs both in vitro and in vivo. First, glyphosate was compared with representative GBHs namely MON 52276 (EU), MON 76473 (UK) and MON 76207 (USA) using the mammalian stem cell-based ToxTracker system. Here, MON 52276 and MON 76473, but not glyphosate and MON 76207, activated oxidative stress and unfolded protein responses. Second, molecular profiling of liver was performed in female Sprague-Dawley rats exposed to glyphosate or MON 52276 (both at 0.5, 50, 175 mg/kg bw/day glyphosate) for 90 days. MON 52276 but not glyphosate increased hepatic steatosis and necrosis. MON 52276 and glyphosate altered the expression of genes in liver reflecting TP53 activation by DNA damage and circadian rhythm regulation. Genes most affected in liver were similarly altered in kidneys. Small RNA profiling in liver showed decreased amounts of miR-22 and miR-17 from MON 52276 ingestion. Glyphosate decreased mir-30 while miR-10 levels were increased. DNA methylation profiling of liver revealed 5,727 and 4,496 differentially methylated CpG sites between the control and glyphosate and MON 52276 exposed animals respectively. Apurinic/apyrimidinic DNA damage formation in liver was increased with glyphosate exposure. Altogether, our results show that Roundup formulations cause more biological changes linked with carcinogenesis than glyphosate. 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. Limited selection towards lighter seeds in Danish grasslands over an eight-year period Time-series analysis of cover data from 236 Danish grassland sites demonstrated community selection on seed mass. Across four grassland habitat types and during an eight-year period there was significant community selection towards lighter seeds. In the analysis, the continuous seed size variable was used to group plant species into functional types. This method was chosen in order to account for the sampling process of the abundance data, and is in contrast to most other analyses of trait selection, where the community weighted mean of the traits is used as the dependent variable. Generally, little is known on the underlying causes of possible plant trait community selection in semi-natural grasslands that are not subjected to successional processes. Programmed knockout mutation of liver fluke granulin, Ov-grn-1, impedes malignant transformation during chronic opisthorchiasis Infection with the food-borne liver fluke Opisthorchis viverrini is the principal risk factor for cholangiocarcinoma in the Mekong Basin countries of Thailand, Lao PDR, Vietnam, Myanmar and Cambodia. Using a novel model of CCA, involving infection with gene-edited liver flukes in the hamster during concurrent exposure to a nitrosamine, we explored the role of the fluke granulin-like growth factor Ov-GRN-1 in malignancy. We produced programmed gene knockout flukes ({Delta}Ov-grn-1) by delivery of a CRISPR/Cas9/gRNA system by electroporation. Genome sequencing confirmed Cas9-catalyzed mutations in the targeted genes, which was accompanied by rapid depletion of transcripts and the cognate proteins. Whereas Ov-grn-1 gene-edited parasites colonized the biliary tract and developed into adult flukes, less hepatobiliary tract disease manifested during chronic infection with{Delta} Ov-grn-1 worms in comparison to hamsters infected with control parasites. Specifically, immunohistochemical analysis of thin sections of livers revealed markedly less periductal fibrosis surrounding the flukes and less liver fibrosis globally during infection with{Delta} Ov-grn-1 genotype worms, minimal biliary epithelial cell proliferation, and markedly fewer mutations of TP53 in biliary epithelial cells. Moreover, fewer hamsters developed high-grade cholangiocarcinoma when infected with the{Delta} Ov-grn-1 flukes compared to controls. The clinically-relevant, pathophysiological phenotype of the hepatobiliary tract confirmed a role for this secreted growth factor in malignancy and morbidity during opisthorchiasis. Interfacial tension driven open droplet microfluidics We developed an open channel droplet microfluidic system that autonomously generates droplets at low Ca (~10-4-10-3) by leveraging competing hydrostatic and capillary pressure. With only our open channel polytetrafluoroethylene (PTFE) device, pipettes, and commercially available carrier fluid, we produce hundreds of microliter droplets; tubing, electronics, or pumps are not required, making droplet technology feasible for research labs without external flow generators. Furthermore, we demonstrated conceptual applications that showcase the process of droplet generation, splitting, transport, incubation, mixing, and sorting in our system. Unlike conventional droplet microfluidics, the open nature of the device enables the use of physical tools such as tweezers and styli to directly access the system; with this, we developed a new method of droplet sorting and transfer that capitalizes on the Cheerios effect, the aggregation of buoyant objects along a liquid interface. Our platform offers enhanced usability, direct access to the droplet contents, easy manufacturability, compact footprint, and high customizability. This design is a first step in exploring the space of power-free open droplet microfluidic systems and provide design rules for similar channel designs. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector.