Model sharing in the human medial temporal lobe Effective planning involves knowing where different actions will take us. However natural environments are rich and complex, leading to an exponential increase in memory demand as a plan grows in depth. One potential solution to this problem is to share the neural state transition functions used for planning between similar contexts. Here, we asked human participants to perform a sequential decision making task designed so that knowledge could be shared between some contexts but not others. Computational modelling showed that participants shared a model of state transitions between contexts where appropriate. fMRI data identified the medial temporal lobe as a locus for learning of state transitions, and within the same region, correlated BOLD patterns were observed in contexts where state transition information were shared. Finally, we show that the transition model is updated more strongly following the receipt of positive compared to negative outcomes, a finding that challenges conventional theories of planning which assume knowledge about our environment is updated independently of outcomes received. Together, these findings propose a computational and neural account of how information relevant for planning can be shared between contexts.
Hierachical Resampling for Bagging in Multi-Study Prediction with Applications to Human Neurochemical Sensing We propose the "study strap ensemble," which combines advantages of two common approaches to fitting prediction models when multiple training datasets ("studies") are available: pooling studies and fitting one model versus averaging predictions from multiple models each fit to individual studies. The study strap ensemble fits models to bootstrapped datasets, or "pseudo-studies." These are generated by resampling from multiple studies with a hierarchical resampling scheme that generalizes the randomized cluster bootstrap. The study strap is controlled by a tuning parameter that determines the proportion of observations to draw from each study. When the parameter is set to its lowest value, each pseudo-study is resampled from only a single study. When it is high, the study strap ignores the multi-study structure and generates pseudo-studies by merging the datasets and drawing observations like a standard bootstrap. We empirically show the optimal tuning value often lies in between, and prove that special cases of the study strap draw the merged dataset and the set of original studies as pseudo-studies. We extend the study strap approach with an ensemble weighting scheme that utilizes information in the distribution of the covariates of the test dataset.

Our work is motivated by neuroscience experiments using real-time neurochemical sensing during awake behavior in humans. Current techniques to perform this kind of research require measurements from an electrode placed in the brain during awake neurosurgery and rely on prediction models to estimate neurotransmitter concentrations from the electrical measurements recorded by the electrode. These models are trained by combining multiple datasets that are collected in vitro under heterogeneous conditions in order to promote accuracy of the models when applied to data collected in the brain. A prevailing challenge is deciding how to combine studies or ensemble models trained on different studies to enhance model generalizability.

Our methods produce marked improvements in simulations and in this application. All methods are available in the studyStrap CRAN package.
Transcripts switched off at the stop of phloem unloading highlight the energy efficiency of sugar import in the ripening V. vinifera fruit. Transcriptomic changes at the cessation of sugar accumulation in the pericarp of Vitis vinifera were addressed on single berries re-synchronized according to their individual growth patterns. The net rates of water, sugars and K+ accumulation inferred from individual growth and solute concentration confirmed that these inflows stopped simultaneously in the ripe berry, while the small amount of malic acid remaining at this stage was still being oxidized at low rate. Re-synchronized individual berries displayed negligible variations in gene expression among triplicates. RNA-Seq studies revealed sharp reprogramming of cell wall enzymes and structural proteins at the stop of phloem unloading, associated with an 80% repression of multiple sugar transporters and aquaporins on the plasma or tonoplast membranes, with the noticeable exception of H+/sugar symporters, that were rather weakly and constitutively expressed. This was verified in three genotypes placed in contrasted thermo-hydric conditions. The prevalence of SWEET suggests that electrogenic transporters would play a minor role on the plasma membranes of SE/CC complex and the one of the flesh, while sucrose/H+ exchangers dominate on its tonoplast. Cis-regulatory elements present in their promoters allowed to sort these transporters in different groups, also including specific TIPs and PIPs paralogs, and cohorts of cell wall related genes. Together with simple thermodynamic considerations, these results lead to propose that H+/sugar exchangers at the tonoplast, associated with a considerably acidic vacuolar pH, may exhaust cytosolic sugars in the flesh and alleviate the need for supplementary energization of sugar transport at the plasma membrane. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. The compound of baicalein, wogonin, and oroxylin-A inhibits EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. On Periodic Gene Expression in Budding Yeast Metabolic oscillation percolate throughout cellular physiology. The transcriptome oscillates between gene cohorts encoding for biosynthesis and growth, and for catabolism and stress response. However, long protein half-lives dampen effects on protein abundance, and the function of periodic transcription remains elusive. To test three prevailing hypotheses, we performed RNA-seq analysis during a dynamic state of the system, just preceding a bifurcation of system dynamics to a longer period. A few enzymes with short protein half-live and high transcript abundance amplitude could play a direct role in metabolic dynamics (hypothesis H1), e.g., carbonic anhydrase, sulfate uptake and the glyoxylate cycle. Substantial changes in the relative duration of expression of the antagonistic cohorts precede the bifurcation, supporting the idea of a function in cellular resource allocation (H2). The pulse-width modulation model, a mathematical formulation of this idea, explains the dependence of the cycle period on the culture growth rate. This pulse-like model of cell growth provides a first theoretical framework, where the phenomenon is understood as a mechanism of cellular resource allocation and protein homeostasis. The universal temporal program of transcription encodes the spatial structure of the cell, and we further suggest a function in pattern formation of the growing cell (H3). Genome-wide association study of musical beat synchronization demonstrates high polygenicity Moving in synchrony to the beat is a fundamental component of musicality. Here, we conducted a genome-wide association study (GWAS) to identify common genetic variants associated with beat synchronization in 606,825 individuals. Beat synchronization exhibited a highly polygenic architecture, with sixty-nine loci reaching genome-wide significance (p<5x10-8) and SNP-based heritability (on the liability scale) of 13%-16%. Heritability was enriched for genes expressed in brain tissues, and for fetal and adult brain-specific gene regulatory elements, underscoring the role of central nervous system-expressed genes linked to the genetic basis of the trait. We performed validations of the self-report phenotype (through internet-based experiments) and of the GWAS (polygenic scores for beat synchronization were associated with patients algorithmically classified as musicians in medical records of a separate biobank). Genetic correlations with breathing function, motor function, processing speed, and chronotype suggest shared genetic architecture with beat synchronization and provide avenues for new phenotypic and genetic explorations. Dietary palmitic acid induces innate immune memory via ceramide production that enhances severity of acute septic shock and clearance of infection Trained immunity is an innate immune memory response that is induced by primary microbial or sterile stimuli that sensitizes monocytes and macrophages to a secondary pathogenic challenge, reprogramming the host response to infection and inflammatory disease. Nutritional components, such as dietary fatty acids, can act as inflammatory stimuli, but it is unknown if they can act as the primary stimuli in the context of innate immune memory. Here we find mice fed diets enriched in saturated fatty acids (SFAs) confer a hyper-inflammatory response to systemic lipopolysaccharide (LPS) and increased mortality, independent of diet-induced microbiome and glycemic modulation. Ex vivo, we show monocytes and splenocytes from mice fed enriched SFAs do not have altered baseline inflammation, but enhanced responses to a secondary inflammatory challenge. Lipidomics identified enhanced free palmitic acid (PA) and PA-associated lipids in SFA-fed mice serum. We found pre-treatment with physiologically relevant concentrations of PA alone reprograms macrophages to induce a hyper-inflammatory response to secondary challenge with LPS. This response was found to be dependent on the synthesis of ceramide, and reversible when treated with oleic acid, a mono-unsaturated FA that depletes intracellular ceramide. In vivo, we found systemic PA confers enhanced inflammation and mortality during an acute septic response to systemic LPS, which was not reversible for up to 7 days post-PA-exposure. While PA-treatment is harmful for acute septic shock outcome, we find PA exposure enhanced clearance of Candida albicans in RAG-/- mice. These are the first data to implicate enriched dietary SFAs, and specifically PA, in the induction of long-lived innate immune memory that is detrimental during an acute septic response, but beneficial for clearance of pathogens. Endogenous retroviruses co-opted as divergently transcribed regulatory elements shape the regulatory landscape of embryonic stem cells Transposable elements are an abundant source of transcription factor binding sites and favorable genomic integration may lead to their recruitment by the host genome for gene regulatory functions. However, it is unclear how frequent co-option of transposable elements as regulatory elements is, to which regulatory programs they contribute and how they compare to regulatory elements devoid of transposable elements. Here, we report a transcription initiation-centric, in-depth characterization of the transposon-derived regulatory landscape of mouse embryonic stem cells. We demonstrate that a substantial number of transposable elements, in particular endogenous retroviral elements, carry open chromatin regions that are divergently transcribed into unstable RNAs in a cell-type specific manner, and that these elements contribute to a sizable proportion of active enhancers and gene promoters. We further show that transposon subfamilies contribute differently and distinctly to the pluripotency regulatory program through their repertoires of transcription factor binding sites, shedding light on the formation of regulatory programs and the origins of regulatory elements. A dysmorphic mouse model reveals developmental interactions of chondrocranium and dermatocranium The cranial endo- and dermal skeletons, which comprise the vertebrate skull, evolved independently over 470 million years ago and form separately during embryogenesis. In mammals, much of the cartilaginous chondrocranium is transient, undergoing endochondral ossification or disappearing, so its role in skull morphogenesis is not well studied and it remains an enigmatic structure. We provide complete three-dimensional (3D) reconstructions of the laboratory mouse chondrocranium from embryonic day 13.5 through 17.5 using a novel methodology of uncertainty-guided segmentation of phosphotungstic enhanced 3D microcomputed tomography images with sparse annotation. We evaluate the embryonic mouse chondrocranium and dermatocranium in 3D and delineate the effects of a Fgfr2 variant on embryonic chondrocranial cartilages and on the association with forming dermal bones using the Fgfr2cC342Y/+ Crouzon syndrome mouse. We show that the dermatocranium develops outside of and in shapes that conform to the chondrocranium. Results reveal direct effects of the Fgfr2 variant on embryonic cartilage, on chondrocranium morphology, and on the association between chondrocranium and dermatocranium development. Histologically we observe a trend of relatively more chondrocytes, larger chondrocytes, and/or more matrix in the Fgfr2cC342Y/+ embryos at all timepoints before the chondrocranium begins to disintegrate at E16.5. The chondrocrania and forming dermatocrania of Fgfr2cC342Y/+ embryos are relatively large, but a contrasting trend begins at E16.5 and continues into early postnatal (P0 and P2) timepoints, with the skulls of older Fgfr2cC342Y/+ mice reduced in most dimensions compared to Fgfr2c+/+ littermates. Our findings have implications for the study and treatment of human craniofacial disease, for understanding the impact of chondrocranial morphology on skull growth, and potentially on the evolution of skull morphology. 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. 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. Competition and environmental gradients structure assemblages of Canidae across the planet The phylogenetic information of assemblages carries the signature of ecological and evolutionary processes that assembled these communities. Identifying the mechanisms that shape communities is not simple, as they can vary spatially. Here, we investigated how the phylogenetic structure of Canidae across the globe is affected by the environment and competition. We first identified phylogenetically clustered and overdispersed assemblages of canids over the planet taking into account regional pools of species in South and North America, Eurasia, and Africa. Then, we applied Structural Equation Models in these communities in order to identify the effect of temperature, vegetation cover, human impact, and body size dissimilarity on the spatial distribution of the phylogenetic information of canids. Only southern South America and the Middle East had strong phylogenetic clustered patterns, while the rest of the planet predominantly showed overdispersed assemblages. Body size dissimilarity and vegetation cover were the most important variables to explain the spatial patterns of phylogenetic structure in clustered and overdispersed communities, respectively. Canidae community composition across the world presents significant patterns of clustering and overdispersion, which vary following mainly the environmental gradient, suggesting habitat filtering as the main force acting on Canidae assemblages. Bayesian Additive Regression Trees for Genotype by Environment Interaction Models We propose a new class of models for the estimation of genotype by environment (GxE) interactions in plant-based genetics. Our approach, named AMBARTI, uses semi-parametric Bayesian additive regression trees to accurately capture marginal genotypic and environment effects along with their interaction in a fully Bayesian model. We demonstrate that our approach is competitive or superior to similar models widely used in the literature via both simulation and a real world data set. Furthermore, we introduce new types of visualisation to properly assess both the marginal and interactive predictions from the model. An R package that implements our approach is available at https://github.com/ebprado/ambarti. Daily torpor reduces the energetic consequences of habitat selection for a widespread bat Homeothermy requires increased metabolic rates as temperatures decline below the thermoneutral zone, so homeotherms typically select microhabitats within or near their thermoneutral zones during periods of inactivity. However, many mammals and birds are heterotherms that relax internal controls on body temperature when maintaining a high, stable body temperature is energetically costly. Such heterotherms should be less tied to microhabitats near their thermoneutral zones, and because heterotherms spend more time in torpor and expend less energy at colder temperatures, heterotherms may even select microhabitats in which temperatures are well below their thermoneutral zones. We studied how temperature and daily torpor influence selection of diurnal roosts by a heterothermic bat (Myotis thysanodes). We (1) quantified the relationship between ambient temperature and daily duration of torpor, (2) simulated daily energy expenditure over a range of microhabitat (roost) temperatures, and (3) quantified the influence of roost temperature on roost selection. While warm roosts substantially reduced energy expenditure of simulated homeothermic bats, heterothermic bats modulated their use of torpor to maintain a constant level of energy expenditure over the course of a day. Daily torpor expanded the range of energetically economical microhabitats, such that roost selection was independent of roost temperature. Our work adds to a growing literature documenting functions of torpor beyond its historical conceptualization as a last-resort measure to save energy during extended or acute energetic stress. Force profile of the two-handed hardstyle kettlebell swing performed by an RKC-certified Instructor BackgroundThe effects of hardstyle kettlebell training are increasingly cited in strength and conditioning research, yet reference data from a proficient swing is scarce. The aim of this exploratory study was to investigate the force profile of a two-handed hardstyle swing performed by a Russian Kettlebell Challenge (RKC) instructor.

MethodsThe subject is a 44-year-old male, body mass 75.6 kg, height 173.5 cm, with six years of regular hardstyle training experience. Two-handed hardstyle swings were performed with a series of incremental mass kettlebells (8-68 kg). Ground reaction force (GRFs) was obtained from a floor-mounted force platform. Force-time curves (FTCs), peak force, forward force, rate of force development (RFD) and swing cadence were investigated.

ResultsData revealed the FTC of a proficient swing is highly consistent and dominated by a single force peak (mean SD = 47 N), with a profile that remained largely unchanged to 24 kg. Pearson correlation analysis revealed a very strong positive correlation in peak force with kettlebell mass (r = 0.95), which increased disproportionately from the lightest to heaviest kettlebells; net peak force increased from 8.36 {+/-} 0.75 N.kg-1 (0.85 x BW) to 12.82 {+/-} 0.39 N.kg-1 (1.3x BW). There was a strong negative correlation between RFD and kettlebell mass (r = 0.82) that decreased from 39.2 N.s-1.kg-1 to 21.5 N.s-1.kg-1. There was a very strong positive correlation in forward ground reaction force with kettlebell mass (r = 0.99), expressed as a ratio of vertical ground reaction, that increased from 0.092 (9.2%) to 0.205 (20.5%). Swing cadence exceeded 40 swings per minute (SPM) with all kettlebells.

ConclusionOur findings challenge some of the popular beliefs of the hardstyle kettlebell swing. Consistent with hardstyle practice, and previous kinematic analysis of expert and novice, force-time curves show a characteristic single large force peak, differentiating passive from active shoulder flexion. Ground reaction force did not increase proportionate to kettlebell mass, with a magnitude of forward force smaller than described in practice. These results could be useful for coaches and trainers wanting to improve athletic performance, and healthcare providers using the kettlebell swing for therapeutic purposes. Findings from this study were used to inform the BELL Trial, a pragmatic controlled trial of kettlebell training with older adults. www.anzctr.org.au ACTRN12619001177145. A combined human in silico and CRISPR/Cas9-mediated in vivo zebrafish based approach for supporting gene target validation in early drug discovery The clinical heterogeneity of heart failure has challenged our understanding of the underlying genetic mechanisms of this disease. In this respect, large-scale patient DNA sequencing studies have become an invaluable strategy for identifying potential genetic contributing factors. The complex aetiology of heart failure, however, also means that in vivo models are vital to understand the links between genetic perturbations and functional impacts. Traditional approaches (e.g. genetically-modified mice) are optimal for assessing small numbers of proposed target genes, but less practical when multiple targets are identified. The zebrafish, in contrast, offers great potential for higher throughput in vivo gene functional assessment to aid target prioritisation and support definitive studies undertaken in mice. Here we used whole-exome sequencing and bioinformatics on human patient data to identify 3 genes (API5, HSPB7, and LMO2) suggestively associated with heart failure that were also predicted to play a broader role in disease aetiology. The role of these genes in cardiovascular system development and function was then further investigated using in vivo CRISPR/Cas9-mediated gene mutation analysis in zebrafish. We observed multiple impacts in F0 knockout zebrafish embryos (crispants) following effective somatic mutation, including reductions in ventricle size, pericardial oedema, and chamber malformation. In the case of lmo2, there was also a significant impact on cardiovascular function as well as an expected reduction in erythropoiesis. The data generated from both the human in silico and zebrafish in vivo assessments undertaken supports roles for API5, HSPB7, and LMO2 in human cardiovascular disease and identifies them as potential drug targets for further investigation. The data presented also supports the use of human in silico genetic variant analysis, in combination with zebrafish crispant phenotyping, as a powerful approach for assessing gene function as part of an integrated multi-level drug target validation strategy. Curation, inference, and assessment of a globally reconstructed gene regulatory network for Streptomyces coelicolor Streptomyces coelicolor A3(2) is a model microorganism for the study of Streptomycetes, antibiotic production, and secondary metabolism in general. Even though S. coelicolor has an outstanding variety of regulators among bacteria, little effort to globally study its transcription has been made. We manually curated 29 years of literature and databases to assemble a meta-curated experimentally-validated gene regulatory network (GRN) with 5386 genes and 9707 regulatory interactions ([~]41% of the total expected interactions). This provides the most extensive and up-to-date reconstruction available for the regulatory circuitry of this organism. Only [~]6% (534/9707) are supported by experiments confirming the binding of the transcription factor to the upstream region of the target gene, the so-called "strong" evidence. While for the remaining interactions there is no confirmation of direct binding. To tackle network incompleteness, we performed network inference using several methods (including two proposed here) for motif identification in DNA sequences and GRN inference from transcriptomics. Further, we contrasted the structural properties and functional architecture of the networks to assess the reliability of the predictions, finding the inference from DNA sequence data to be the most trustworthy approach. Finally, we show two applications of the inferred and the curated networks. The inference allowed us to propose novel transcription factors for the key Streptomyces antibiotic regulatory proteins (SARPs). The curated network allowed us to study the conservation of the system-level components between S. coelicolor and Corynebacterium glutamicum. There we identified the basal machinery as the common signature between the two organisms. The curated networks were deposited in Abasy Atlas (https://abasy.ccg.unam.mx/) while the inferences are available as Supplementary Material. A Novel Organoid Model of In Vitro Spermatogenesis Using Human Induced Pluripotent Stem Cells Infertility is thought to be caused by genetic mutations and dysfunction in the cellular niche where spermatogenesis takes place. An understanding of the specialized cellular processes which drive spermatogenesis is needed to develop treatments; however, the development of in vitro systems to study these cells has been hindered by our reliance on rarely available human testicular tissues for research. Human induced pluripotent stem cells (hiPSCs) can be used to derive human testicular-like cells, and thus provide an avenue for the development of in vitro testicular model systems. Therefore, this study set out to engineer a human testicular tissue model using hiPSCs for the first time. We demonstrate the ability of hiPSC-derived testicular cells to self-organize and mature into testicular-like tissues using organoid culture. Moreover, we show that hiPSC-derived testicular organoids promote testicular somatic cell maturation and spermatogenesis up to the post-meiotic spermatid stage. These hiPSC-derived testicular organoids have the potential to replace rarely available primary testicular tissues to further infertility research in an in vitro setting. A Substrate-triggered μ-Peroxodiiron(III) Intermediate in the 4-Chloro-L-Lysine-Fragmenting Heme-Oxygenase-like Diiron Oxidase (HDO) BesC: Substrate Dissociation from, and C4 Targeting by, the Intermediate The enzyme BesC from the {beta}-ethynyl-L-serine biosynthetic pathway in Streptomyces cattleya fragments 4-chloro-L-lysine (produced from L-Lysine by BesD) to ammonia, formaldehyde, and 4-chloro-L-allylglycine and can analogously fragment L-Lys itself. BesC belongs to the emerging family of O2-activating non-heme-diiron enzymes with the "heme-oxygenase-like" protein fold (HDOs). Here we show that binding of L-Lys or an analog triggers capture of O2 by the proteins diiron(II) cofactor to form a blue {micro}-peroxodiiron(III) intermediate analogous to those previously characterized in two other HDOs, the olefin-installing fatty acid decarboxylase, UndA, and the guanidino-N-oxygenase domain of SznF. The [~] 5- and [~] 30-fold faster decay of the intermediate in reactions with 4-thia-L-Lys and (4RS)-chloro-DL-lysine than in the reaction with L-Lys itself, and the primary deuterium kinetic isotope effects (D-KIEs) on decay of the intermediate and production of L-allylglycine in the reaction with 4,4,5,5-[2H]-L-Lys, imply that the peroxide intermediate or a successor complex with which it reversibly interconverts initiates the oxidative fragmentation by abstracting hydrogen from C4. Surprisingly, the sluggish substrate L-Lys can dissociate after triggering the intermediate to form, thereby allowing one of the better substrates to bind and react. Observed linkage between Fe(II) and substrate binding suggests that the triggering event involves a previously documented (in SznF) ordering of the dynamic HDO architecture that contributes one of the iron sites, a hypothesis consistent with the observation that the diiron(III) product cluster produced upon decay of the intermediate spontaneously degrades, as it has been shown to do in all other HDOs studied to date. Heterogeneity of circulating epithelial cells in breast cancer at single-cell resolution: identifying tumor and hybrid cells Circulating tumor cells and hybrid cells formed by the fusion of tumor cells with normal cells are leading players in metastasis and have prognostic relevance. Circulating tumor cells and hybrid cells are identified as CD45-negative and CD45-positive epithelial cells. However, such an approach is challenging because epithelial cells are observed in the blood of healthy individuals. In this study, we applied single-cell RNA sequencing to profile CD45-negative and CD45-positive circulating epithelial cells (CECs) in 20 breast cancer patients and one healthy donor. DNA ploidy analysis was used to identify the tumor and hybrid cells among CD45[boxh] and CD45+ CECs in patients, respectively. Functional enrichment analysis was applied to characterize aneuploid and diploid cells. Diploid cells were also annotated to generate cell-type candidates and analyzed for copy-number aberrations (CNAs) to confirm or refute their tumor origin. CD45[boxh] and CD45+ CECs were found in cancer patients (25.5 (range 0-404) and median 6.5 (0-147)) and the healthy donor (8 and 11 cells) and divided into three clusters. Two CD45[boxh] CEC clusters were predominantly aneuploid (97% and 98%), but one cluster contained more diploid (59%) than aneuploid cells. CD45+ CECs were mostly diploid: only clusters 1 and 2 had aneuploid cells (16% and 2%). Diploid CD45[boxh] and CD45+ CECs were annotated as different immune cells and surprisingly harbored many CNAs. Cancer-associated signaling pathways were found only in aneuploid cells of CD45[boxh] CEC cluster 1 and diploid cells of CD45+ CEC cluster 1. Thus, our findings suggest that CECs in breast cancer patients are a highly heterogeneous population comprising aneuploid (tumor and hybrid) and diploid (normal) cells. DNA ploidy analysis is an effective instrument for identifying tumor and hybrid cells among CD45[boxh] and CD45+ CECs, respectively. 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. Atypical integrative element with strand-biased circularization activity assists interspecies antimicrobial resistance gene transfer from Vibrio alfacsensis The exchange of antimicrobial resistance (AMR) genes between aquaculture and terrestrial microbial populations has emerged as a serious public health concern. However, the nature of the mobile genetic elements in marine bacteria is poorly documented. To gain insight into the genetic mechanisms underlying AMR gene transfer from marine bacteria, we mated a multidrug-resistant Vibrio alfacsensis strain with an Escherichia coli strain, and then determined the complete genome sequences of the donor and transconjugant strains. Sequence analysis revealed a multidrug resistance conjugative plasmid in the donor strain, which was integrated into the chromosome of the recipient. The plasmid backbone in the transconjugant chromosome was flanked by two copies of a 7.1 kb unclassifiable integrative element harboring a {beta}-lactamase gene. The 7.1 kb element (named SE-6945) and previously reported Tn6283 share four coding sequences, two of which encode the catalytic R-H-R-Y motif of tyrosine recombinases. Polymerase chain reaction and sequencing experiments revealed that SE-6945 generates a circular copy of its one specific strand without generating an empty donor site, in contrast to the movement of integron gene cassette or ICE/IMEs discovered to date. Copy numbers and location of SE-6945 in the chromosome affected antibiotic resistance levels of the transconjugants. SE-6945/Tn6283-related integrative elements are collectively termed SEs (strand-biased circularizing integrative elements). SEs were identified in the genomes of other Vibrio species. Overall, these results suggest that in addition to conjugative plasmids, SEs are involved in the spread of AMR genes among marine bacteria. 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. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

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

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Cryptic pathogen-sugar interactions revealed by universal saturation transfer analysis Many host pathogen interactions such as human viruses (including non-SARS-coronaviruses) rely on attachment to host cell-surface glycans. There are conflicting reports about whether the Spike protein of SARS-CoV-2 binds to sialic acid commonly found on host cell-surface N-linked glycans. In the absence of a biochemical assay, the ability to analyze the binding of glycans to heavily- modified proteins and resolve this issue is limited. Classical Saturation Transfer Difference (STD) NMR can be confounded by overlapping sugar resonances that compound with known experimental constraints. Here we present  universal saturation transfer analysis (uSTA), an NMR method that builds on existing approaches to provide a general and automated workflow for studying protein-ligand interactions. uSTA reveals that B-origin-lineage-SARS-CoV-2 spike trimer binds sialoside sugars in an  end on manner and modelling guided by uSTA localises binding to the spike N-terminal domain (NTD). The sialylated-polylactosamine motif is found on tetraantennary human N-linked-glycoproteins in deeper lung and may have played a role in zoonosis. Provocatively, sialic acid binding is abolished by mutations in some subsequent SARS- CoV-2 variants-of-concern. A very high resolution cryo-EM structure confirms the NTD location and  end on mode; it rationalises the effect of NTD mutations and the structure-activity relationship of sialic acid analogues. uSTA is demonstrated to be a robust, rapid and quantitative tool for analysis of binding, even in the most demanding systems.

Extended AbstractThe surface proteins found on both pathogens and host cells mediate entry (and exit) and influence disease progression and transmission. Both types can bear host-generated post- translational modifications such as glycosylation that are essential for function but can confound biophysical methods used for dissecting key interactions. Several human viruses (including non- SARS-coronaviruses) attach to host cell-surface N-linked glycans that include forms of sialic acid (sialosides). There remains, however, conflicting evidence as to if or how SARS-associated coronaviruses might use such a mechanism. Here, we demonstrate quantitative extension of  saturation transfer protein NMR methods to a complete mathematical model of the magnetization transfer caused by interactions between protein and ligand. The method couples objective resonance-identification via a deconvolution algorithm with Bloch-McConnell analysis to enable a structural, kinetic and thermodynamic analysis of ligand binding beyond previously-perceived limits of exchange rates, concentration or system. Using an automated and openly available workflow this  universal saturation transfer analysis (uSTA) can be readily-applied in a range of even heavily-modified systems in a general manner to now obtain quantitative binding interaction parameters (KD, kEx). uSTA proved critical in mapping direct interactions between natural sialoside sugar ligands and relevant virus-surface attachment glycoproteins - SARS-CoV-2-spike and influenza-H1N1-haemagglutinin variants - by quantitating ligand signal in spectral regions otherwise occluded by resonances from mobile protein glycans (that also include sialosides). In B- origin-lineage-SARS-CoV-2 spike trimer  end on-binding to sialoside sugars was revealed contrasting with  extended surface-binding for heparin sugar ligands; uSTA-derived constraints used in structural modelling suggested sialoside-glycan binding sites in a beta-sheet-rich region of spike N-terminal domain (NTD). Consistent with this, uSTA-glycan binding was minimally- perturbed by antibodies that neutralize the ACE2-binding domain (RBD) but strongly disrupted in spike from the B1.1.7/alpha and B1.351/beta variants-of-concern, which possess hotspot mutations in the NTD. Sialoside binding in B-origin-lineage-NTD was unequivocally pinpointed by cryo-EM to a site that is created from residues that are notably deleted in variants (e.g. H69,V70,Y145 in alpha). An analysis of beneficial genetic variances in cohorts of patients from early 2020 suggests a model in which this site in the NTD of B-origin-lineage-SARS-CoV-2 (but not in alpha/beta-variants) may have exploited a specific sialylated-polylactosamine motif found on tetraantennary human N-linked-glycoproteins in deeper lung. Together these confirm a novel binding mode mediated by the unusual NTD of SARS-CoV-2 and suggest how it may drive virulence and/or zoonosis via modulation of glycan attachment. Since cell-surface glycans are widely relevant to biology and pathology, uSTA can now provide ready, quantitative, widespread analysis of complex, host-derived and post-translationally modified proteins with putative ligands relevant to disease even in previously confounding complex systems.