Global structural stability and the role of cooperation in mutualistic system Dynamical systems on graphs allow to describe multiple phenomena from different areas of Science. In particular, many complex systems in Ecology are studied by this approach. In this paper we analize the mathematical framework for the study of the structural stability of each stationary point, feasible or not, introducing a generalization for this concept, defined as Global Structural Stability. This approach would fit with the proper mathematical concept of structural stability, in which we find a full description of the complex dynamics on the phase space due to nonlinear dynamics. This fact can be analyzed as an informational field grounded in a global attractor whose structure can be completely characterized. These attractors are stable under perturbation and suppose the minimal structurally stable sets. We also study in detail, mathematically and computationally, the zones characterizing different levels of biodiversity in bipartite graphs describing mutualistic antagonistic systems of population dynamics. In particular, we investigate the dependence of the region of maximal biodiversity of a system on its connectivity matrix. On the other hand, as the network topology does not completely determine the robustness of the dynamics of a complex network, we study the correlation between structural stability and several graph measures. A systematic study on synthetic and biological graphs is presented, including 10 mutualistic networks of plants and seed-dispersal and 1000 random synthetic networks. We compare the role of centrality measures and modularity, concluding the importance of just cooperation strength among nodes when describing areas of maximal biodiversity. Indeed, we show that cooperation parameters are the central role for biodiversity while other measures act as secondary supporting functions.

Author summaryWe introduce the concept of Global Structural Stability as a proper mathematical concept to fully understand biodiversity in some ecological systems. Our concept retakes the definitions in the classical works of R. Thom [1] and Andronov-Pontryagin [2]. Moreover, there exists a close relation between the structure of a complex network, described as a graph, and its associated dynamics. Mutualistic networks introduce cooperation links between two groups of species, as plant and pollinators or seed-dispersal. The understanding of organizational aspects leading to maximizing biodiversity is one of the more important research areas in Theoretical and Applied Ecology. In this work we introduce a systematic study on different graph measures in order to identify optimal organization for maximal biodiversity (defined as structural stability). Our results conclude that, for mutualistic systems, the strength in cooperation parameters are the core fact, i.e., cooperation is the real fact optimizing biodiversity among other possible structural configurations.
Destructive fishing - A ubiquitously used but vague term? Usage and impacts across academic research, media and policy The term "destructive fishing" appears in multiple international policy instruments intended to improve outcomes for marine biodiversity, coastal communities and sustainable fisheries. However, the meaning of "destructive fishing" is often vague, limiting effectiveness in policy. Therefore, in this study we systematically reviewed the use of "destructive fishing" in three record types: academic literature, media articles, and policy documents between 1976-2020. A more detailed analysis was performed on sub-sets of these records, considering the extent to which the term is characterised, geographic distribution of use, and specific impacts and practices associated with the term. We found that use of "destructive fishing" relative to the generic term  fisheries has increased since the 1990s. Records focused predominantly on fishing practices in South-eastern Asia, followed by Southern Asia and Europe. The term was characterised in detail in only 15% of records. Habitat damage and blast/poison fishing were the most associated ecological impacts and gear/practices, respectively. Bottom trawling and unspecified net fishing were regularly linked to destructive fishing. Importantly, the three record types use the term differently. Academic literature tends to specifically articulate the negative impacts, while media articles focus generally on associated gears/practices. Significant regional variation also exists in how the term is used and what phenomena it is applied to. This study provides evidence and recommendations to inform stakeholders in any future pursuit of a unified definition of "destructive fishing" to support more meaningful implementation of global sustainability goals.
Cis-Regulatory Hubs: a new 3D Model of Complex Disease Genetics with an Application to Schizophrenia The 3-dimensional (3D) conformation of the chromatin creates complex networks of noncoding regulatory regions (distal elements) and promoters impacting gene regulation. Despite the importance of the role of noncoding regions in complex diseases, little is known about their interplay within regulatory hubs and implication in multigenic diseases like schizophrenia. Here we show that cis-regulatory hubs (CRHs) in neurons highlight functional interactions between distal elements and promoters, providing a model to explain epigenetic mechanisms involved in complex diseases. CRHs represent a new 3D model, where distal elements interact to create a complex network of active genes. In a disease context, CRHs highlighted strong enrichments in schizophrenia-associated genes, schizophrenia-associated SNPs and schizophrenia heritability compared to equivalent structures. Finally, CRHs exhibit larger proportions of genes differentially expressed in schizophrenia compared to promoter-distal element pairs or TADs. CRHs thus capture causal regulatory processes improving the understanding of complex disease etiology such as schizophrenia. These multiple lines of genetic and statistical evidence support CRHs as 3D models to study dysregulation of gene expression in complex diseases more generally.

Summary BlurbGenes and their regulatory elements are organized within 3D networks in neurons which model functional structures and explain schizophrenia genetic etiology. Kidney failure alters parathyroid Pin1 phosphorylation and parathyroid hormone mRNA binding proteins leading to secondary hyperparathyroidism Parathyroid hormone (PTH) regulates calcium metabolism and bone strength. Chronic kidney disease leads to secondary hyperparathyroidism (SHP) which increases morbidity and mortality. In experimental SHP, the increased PTH gene expression is due to enhanced PTH mRNA stability, mediated by changes in its interaction with stabilizing AUF1 and destabilizing KSRP. Pin1 isomerizes target proteins and leads to KSRP dephosphorylation. In SHP, Pin1 isomerase activity is decreased and phosphorylated KSRP fails to bind PTH mRNA, resulting in high PTH mRNA stability and levels. However, the up- and down-stream mechanisms by which kidney failure stimulates the parathyroid remain elusive. We now reveal a pathway where kidney failure induces parathyroid Pin1 phosphorylation, linking Pin1, KSRP and PTH mRNA stability as key players for the high PTH expression in SHP. We identified by mass-spectrometry, changes in rat parathyroid proteome and phosphoproteome profiles induced by impaired renal function, including KSRP phosphorylation at Pin1 target sites. Furthermore, both acute and chronic kidney failure led to parathyroid-specific Pin1 Ser16 and Ser71 phosphorylation, which disrupts Pin1 activity. Accordingly, pharmacologic Pin1 inhibition, that mimics the decreased Pin1 activity in SHP, increased PTH expression ex-vivo in parathyroid organ cultures and in transfected cells, through the PTH mRNA protein-interacting element and KSRP phosphorylation at potential Pin1-binding motifs. Therefore, kidney failure leads to loss of parathyroid Pin1 activity by inducing Pin1 phosphorylation. This predisposes parathyroids to increase PTH production through impaired PTH mRNA decay that is dependent on KSRP phosphorylation. Pin1 and KSRP phosphorylation and the Pin1-KSRP-PTH mRNA axis thus drive SHP. The unfolded protein response triggers the immune deficiency pathway in ticks The insect immune deficiency (IMD) pathway is a defense mechanism that senses and responds to Gram negative bacteria. Ticks lack genes encoding upstream components that initiate the IMD pathway. Despite this deficiency, core signaling molecules are present and functionally restrict tick-borne pathogens. The molecular events preceding activation remain undefined. Here, we show that the Unfolded Protein Response (UPR) initiates the IMD network in Ixodes scapularis ticks. The endoplasmic reticulum (ER) stress receptor, IRE1, is phosphorylated in response to tick-borne bacteria, but does not splice the mRNA encoding XBP1. Instead, through protein modeling and reciprocal pulldowns, we show that Ixodes IRE1 complexes with TRAF2. Disrupting IRE1-TRAF2 signaling blocks IMD pathway activation and diminishes the production of reactive oxygen species. Through in vitro, in vivo, and ex vivo techniques we demonstrate that the UPR-IMD pathway circuitry limits the Lyme disease-causing spirochete Borrelia burgdorferi and the rickettsial agents Anaplasma phagocytophilum and A. marginale (anaplasmosis). Altogether, our study uncovers a novel linkage between the UPR and the IMD pathway in ticks. The adenoviral E1B-55k protein present in HEK293 cells mediates abnormal accumulation of key WNT signaling proteins in large cytoplasmic aggregates. HEK293 cells are one of the most widely used cell lines in research and HEK293 cells are frequently used as an in vitro model for studying the WNT signaling pathway. The HEK293 cell line was originally established by transfection of human embryonic kidney cells with sheared adenovirus 5 DNA and it is known that that HEK293 cells stably express the adenoviral E1A and E1B-55k proteins. Here we show that HEK293 cells display an unexpected distribution of key components of the WNT/{beta}-catenin signaling pathway where AXIN1, APC, DVL2 and tankyrase are all co-localized in large spherical cytoplasmic aggregates. The cytoplasmic aggregates are enclosed by a narrow layer of the adenoviral E1B-55k protein. Reduction of E1B-55k protein levels leads to disappearance of the cytoplasmic aggregates thus corroborating an essential role of the E1B-55k protein in mediating the formation of the aggregates. Furthermore, HEK293 cells with reduced E1B-55k protein levels display reduced levels of transcriptional activation of WNT/{beta}-catenin signaling upon stimulation by the Wnt3A agonist. The demonstrated influence of the E1B-55k protein on the cellular localization of WNT/{beta}-catenin signaling components and on transcriptional regulation of WNT/{beta}-catenin signaling asks for caution in the interpretation of data derived from the HEK293 cell line. Direct cochlear recordings in humans show a theta rhythmic modulation of the auditory nerve by selective attention The architecture of the efferent auditory system enables prioritization of strongly overlapping spatiotemporal cochlear activation patterns elicited by relevant and irrelevant inputs. So far, attempts at finding such attentional modulations of cochlear activity delivered indirect insights in humans or required direct recordings in animals. The extent to which spiral ganglion cells forming the human auditory nerve are sensitive to selective attention remains largely unknown. We investigated this question by testing the effects of attending to either the auditory or visual modality in human cochlear implant (CI) users (3 female, 13 male). Auditory nerve activity was directly recorded with standard CIs during a silent (anticipatory) cue-target interval. When attending the upcoming auditory input, ongoing auditory nerve activity within the theta range (5-8 Hz) was enhanced. Crucially, using the broadband signal (4-25 Hz), a classifier was even able to decode the attended modality from single-trial data. Follow-up analysis showed that the effect was not driven by a narrow frequency in particular. Using direct cochlear recordings from deaf individuals, our findings suggest that cochlear spiral ganglion cells are sensitive to top-down attentional modulations. Given the putatively broad hair-cell degeneration of these individuals, the effects are likely mediated by alternative efferent pathways as compared to previous studies using otoacoustic emissions. Successful classification of single-trial data could additionally have a significant impact on future closed-loop CI developments that incorporate real-time optimization of CI parameters based on the current mental state of the user.

Significance StatementThe efferent auditory system in principle allows top-down modulation of auditory nerve activity, however evidence for this is lacking in humans. Using cochlear recordings in participants performing an audiovisual attention task, we show that ongoing auditory nerve activity in the silent cue-target period is directly modulated by selective attention. Specifically, ongoing auditory nerve activity is enhanced within the theta range when attending upcoming auditory input. Furthermore, over a broader frequency range, the attended modality can be decoded from single-trial data. Demonstrating this direct top-down influence on auditory nerve activity substantially extends previous works that focus on outer hair cell activity. Generally, our work could promote the use of standard cochlear implant electrodes to study cognitive neuroscientific questions. Water levels primarily drive variation in photosynthesis and nutrient use of scrub Red Mangroves in the southeastern Florida Everglades We investigated how mangrove-island micro-elevation (i.e., habitat: center vs. edge) affects tree physiology in a scrub mangrove forest of the southeastern Everglades. We measured leaf gas exchange rates of scrub Rhizophora mangle L. trees monthly during 2019, hypothesizing that CO2 assimilation (Anet) and stomatal conductance (gsw) would decline with increasing water levels and salinity, expecting more-considerable differences at mangrove-island edges than centers, where physiological stress is greatest. Water levels varied between 0 and 60 cm from the soil surface, rising during the wet season (May-October) relative to the dry season (November-April). Porewater salinity ranged from 15 to 30 ppt, being higher at mangrove-island edges than centers. Anet maximized at 15.1 {micro}mol m-2 s-1, and gsw was typically <0.2 mol m-2 s-1, both of which were greater in the dry than the wet season and greater at island centers than edges, with seasonal variability being roughly equal to variation between habitats. After accounting for season and habitat, water level positively affected Anet in both seasons but did not affect gsw. Our findings suggest that inundation stress (i.e., water level) is the primary driver of variation in leaf gas exchange rates of scrub mangroves in the Florida Everglades, while also constraining Anet more than gsw. The interaction between inundation stress due to permanent flooding and habitat varies with season as physiological stress is alleviated at higher-elevation mangrove-island center habitats during the dry season. Freshwater inflows during the wet season, increase water levels and inundation stress at higher-elevation mangrove-island centers, but also potentially alleviate salt and sulfide stress in soils. Thus, habitat heterogeneity leads to differences in nutrient and water acquisition and use between trees growing in island centers versus edges, creating distinct physiological controls on photosynthesis, which likely affect carbon flux dynamics of scrub mangroves in the Everglades. 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/). Polycystins are upregulated in fibrotic kidneys and promote renal fibrosis Autosomal polycystic kidney disease (ADPKD) is caused by mutations in PKD1or PKD2 gene, which encodes polycystin-1 or polycystin-2 respectively. Paradoxically, polycystins are up-regulated in human ADPKD kidneys and overexpression of polycystin-1 even induced renal cyst formation in mice, suggesting that up-regulation of polycystin is pathogenic. In this study, we determined the role of polycystin in renal fibrosis and explored its underlying mechanisms. We found that the expression of polycystin-1 or polycystin-2 was positively correlated with the expression of collagen-I in the kidney of three different fibrotic mouse models. Pharmaceutical inhibition or genetic knockout of polycystin-2 reduced the expression of pro-fibrotic markers in fibrotic kidneys of two animal models. Moreover, conditional knockout of Pkd1 gene attenuated renal fibrosis in two mouse models. Enhancer of Zeste homolog 2 (EZH2) is a methyltransferase which plays important role in renal fibrosis. We found that EZH2 was positively correlated with the expression of polycystins in fibrotic kidneys. Conditional knockout of the Ezh2 attenuated the anti-fibrotic responses induced by Pkd1 deletion in fibrotic kidneys.

In conclusion, polycystins are up-regulated in injured kidneys and polycystins promote renal fibrosis through up-regulation of EZH2 in fibrotic kidneys, suggesting that renal fibrosis in polycystic kidneys is not the direct effect of PKD gene mutations.

SignificanceMutations of polycystins cause autosomal polycystic kidney disease (ADPKD) which is associated with increased renal fibrosis in cystic kidneys. However, the role of polycystin in renal fibrosis remains unclear. Here, we showed that the expression of polycystins is positively correlated with the development of renal fibrosis in injured kidneys. Pharmaceutical inhibition or genetic deletion of polycystin-2 attenuated renal fibrosis in fibrotic kidneys. Moreover, genetic deletion of polycystin-1 inhibited renal fibrosis in two different animal models. Finally, we showed that ablation of the methyltransferase EZH2 mitigated the anti-fibrotic effect of polycystin-1 deletion in fibrotic kidneys. We conclude that polycystin promotes through EZH2 and renal fibrosis in cystic kidney is not the direct effect of PKD gene mutations. Cellular Fitness Phenotypes of Cancer Target Genes from Oncobiology to Cancer Therapeutics To define the growing significance of cellular targets of cancer drugs, we examined the fitness dependency of cellular targets or effectors of cancer drug targets across human cancer cells from 19 cancer types. We observed that the deletion of 35 out of 47 cellular mediators or targets of oncology drugs did not result in the expected loss of cell fitness in appropriate cancer types for which drugs targeting or utilizing these molecules were approved. Additionally, our analysis recognized 43 cellular targets as fitness genes in several cancer types in which these drugs were not approved, and thus, providing clues repurposing approved oncology drugs in cancer types. For example, we found the widespread upregulation and fitness dependency of the components of the mevalonate and purine biosynthesis pathways (currently targeted by bisphosphonates, statins, and pemetrexed in certain cancers) and an association between the overexpression of these targets and reduction in the overall survival duration of patients with breast and other hard-to-treat cancers, for which such drugs are not approved. In brief, the present analysis raised cautions about off-target and undesirable effects of certain oncology drugs in a subset of cancers where the in-tended cellular effectors of drug might not be fitness genes and offers a potential rationale for repurposing certain approved oncology drugs for targeted therapeutics in additional cancer types. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. Conformist social learning leads to self-organised prevention against adverse bias in risky decision making 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. Spatial proteomics reveals subcellular reorganization in human keratinocytes exposed to UVA light The effects of UV light on the skin have been extensively investigated. However, systematic information about how exposure to UVA light, the least energetic but the most abundant UV radiation reaching the Earth, shapes the subcellular organization of proteins is lacking. Using subcellular fractionation, mass spectrometry-based proteomics, machine learning algorithms, immunofluorescence, and functional assays, we mapped the subcellular reorganization of the proteome of human keratinocytes in response to UVA light. Our workflow quantified and assigned subcellular localization for over 1600 proteins, of which about 200 were found to redistribute upon UVA exposure. Reorganization of the proteome affected modulators of signaling pathways, cellular metabolism, and DNA damage response. Strikingly, mitochondria were identified as one of the main targets of UVA-induced stress. Further investigation demonstrated that UVA induces mitochondrial fragmentation, up-regulates redox-responsive proteins and attenuates respiratory rates. These observations emphasize the role of this radiation as a potent metabolic stressor in the skin. Novel Integrative Modeling of Molecules and Morphology across Evolutionary Timescales Evolutionary models account for either population- or species-level processes, but usually not both. We introduce a new model, the FBD-MSC, which makes it possible for the first time to integrate both the genealogical and fossilization phenomena, by means of the multispecies coalescent (MSC) and the fossilized birth-death (FBD) processes. Using this model, we reconstruct the phylogeny representing all extant and many fossil Caninae, recovering both the relative and absolute time of speciation events. We quantify known inaccuracy issues with divergence time estimates using the popular strategy of concatenating molecular alignments, and show that the FBD-MSC solves them. Our new integrative method and empirical results advance the paradigm and practice of probabilistic total evidence analyses in evolutionary biology. Changes in Behavior and Salivary Serum Amyloid A in cows with Subclinical Mastitis The aim of this study was to identify detailed changes in behavior, and in salivary serum amyloid A (SAA), associated with subclinical mastitis. This included standard sickness behaviors (which form part of an adaptive response to conserve energy, minimize heat loss and facilitate recovery following infection and inflammation) and lesser-studied behaviors, that are biologically predicted to change at low-level infection, and therefore particularly relevant for detecting subclinical mastitis (here labelled  luxury behaviors). SAA is a physiological marker of systemic inflammation, with levels in milk and serum already known to increase during subclinical mastitis. We investigated whether the same was true for SAA in cow saliva. Data were collected for n = 34 commercial barn-housed dairy cows, tested in 17 matched-pairs. Each pair comprised a cow with subclinical mastitis (SCM) and a healthy control (CTRL), identified using somatic cell count (SCC) (SCM: SCC >200 x1000 cells/ml; CTRL: SCC <100 x1000 cells/ml). SCM cows were selected for study ad-hoc, at which point they were paired with a CTRL cow, based upon parity and calving date; consequently, the full data set was accrued over several months. Data were collected for each pair over 3 days: SCC (l4:00-l5:00h) Day 1; behavior (24h from 00:00h) Day 2; salivary serum amyloid-A (SAA) Day 3. We report, for the first time, that an increase in salivary SAA occurs during subclinical mastitis; SAA was higher in SCM cows and demonstrated a positive (weak) correlation with SCC. The behavioral comparisons revealed that SCM cows had reductions in activity (behavioral transitions and distance moved), social exploration, social reactivity (here: likelihood to move away/be displaced following receipt of agonism), performance of social grooming and head butts, and the receipt of agonistic non-contact challenges. In addition, SCM cows received more head swipes, and spent a greater proportion of time lying with their head on their flank than CTRL cows. SCM cows also displayed an altered feeding pattern; they spent a greater proportion of feeding time in direct contact with two conspecifics, and a lower proportion of feeding time at self-locking feed barriers, than CTRL cows. Behavioral measures were found to correlate, albeit loosely, with serum SAA in a direction consistent with predictions for sickness behavior. These included positive correlations with both lying duration and the receipt of all agonistic behavior, and negative correlations with feeding, drinking, the performance of all social and all agonistic behavior, and social reactivity. We conclude that changes in salivary SAA, social behavior, and activity offer potential in the detection of subclinical mastitis and recommend further investigation to substantiate and refine our findings. 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. Global population genomics of two subspecies of Cryptosporidium hominis during 500 years of evolution Cryptosporidiosis is a major global health problem and a primary cause of diarrhoea, particularly in young children in low- and middle-income countries (LMICs). The zoonotic Cryptosporidium parvum and anthroponotic C. hominis cause most human infections. Here, we present a comprehensive whole-genome study of C. hominis, comprising 114 isolates from 16 countries within five continents. We detect two lineages with distinct biology and demography, which diverged circa 500 years ago. We consider these lineages two subspecies and propose the names C. hominis hominis and C. hominis aquapotentis (gp60 subtype IbA10G2). In our study, C. h. hominis is almost exclusively represented by isolates from LMICs in Africa and Asia and appears to have undergone recent population contraction. In contrast, C. h. aquapotentis was found in high-income countries, mainly in Europe, North America and Oceania, and appears to be expanding. Notably, C. h. aquapotentis is associated with high rates of direct human-to-human transmission, which may explain its success in countries with well-developed environmental sanitation infrastructure. Intriguingly, we detected genomic regions of introgression following secondary contact between the subspecies. This resulted in high diversity and divergence in genomic islands of putative virulence genes (GIPVs), including muc5 (CHUDEA2_430) and a hypothetical protein (CHUDEA6_5270). This diversity is maintained by balancing selection, suggesting a coevolutionary arms race with the host. Lastly, we find that recent gene flow from C. h. aquapotentis to C. h. hominis, likely associated with increased human migration, may be driving evolution of more virulent C. hominis variants. A memetic algorithm enables efficient local and global all-atom protein-protein docking with backbone and sidechain flexibility Protein complex formation is encoded by specific interactions at the atomic scale, but the computational cost of modeling proteins at this level often requires the use of simplified energy models and limited conformational flexibility. In particular, the use of all-atom energy functions, backbone and sidechain flexibility results in rugged energy landscapes that are difficult to explore. In this study we develop a protein-protein docking algorithm, EvoDOCK, that combine the strength of a differential evolution algorithm for efficient exploration of the global search space with the benefits of a local optimization method to refine detailed atomic interactions. EvoDOCK enabled accurate and fast local and global protein-protein docking using an all-atom energy function with side-chain flexibility. Comparison with a standard method built on Monte Carlo optimization demonstrated improved accuracy and with increases in computational speed of up to 35 times. The evolutionary algorithm also enabled efficient atomistic docking with backbone flexibility. Muscle Activation Profile During Perturbed Walking is Modulated According to Body State During an unexpected loss of balance, avoiding a fall requires people to readjust their footing rapidly and effectively. We investigated the activation patterns of the ankle and knee muscles, and muscle fiber type recruitment resulting from unannounced, medio-lateral (i.e., right/left) horizontal-surface walking perturbations in twenty healthy adults (27.00{+/-}2.79 years, 10 females). Surface electromyography (sEMG) total spectral power for specific frequency bands (40-60Hz, 60-150Hz, 150-250Hz, 250-400Hz and 400-1000Hz), from tibialis anterior (TA) and vastus lateralis (VL) muscles were analyzed. Compared to non-perturbed walking, we found a significant increase in the total spectral power of lower-extremity muscles during the first 3 seconds after perturbation. When two feet were on the ground in time of perturbation we found a different muscle fiber type recruitment pattern between VL and TA muscles. This was not significant for perturbations implemented when one foot was on the ground. Our findings suggest that muscle operating frequency is modulated in real time to fit body state and somatosensory input in the context of functional goal requirements such as a rapid change of footing in response to unexpected loss of balance in single and double-support phases of gait.

New & NoteworthyTo study muscle spectral profiles in response to loss of balance, we investigated the dynamics of muscle spectral power changes, across different frequency bands after unannounced mechanical perturbations during walking. We showed increased activation of high-frequency motor units of the lower-limb muscles, subside 3 seconds after perturbation. Differences in power increase of specific frequency bands suggest that muscle activation is modulated in real time to fit body state in the context of functional goal requirements. 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-catalyzed aggregation of islet amyloid polypeptide is dominated by secondary nucleation Type II diabetes is characterized by the loss of pancreatic {beta}-cells. This loss is thought to be a consequence of membrane disruption, caused by the aggregation of islet amyloid polypeptide (IAPP) into amyloid fibrils. However, the molecular mechanisms of IAPP aggregation in the presence of membranes have remained unclear. Here, we use kinetic analysis to elucidate the aggregation mechanism of IAPP in the presence of mixed zwitterionic and anionic lipid membranes. The results converge to a model in which aggregation on the membrane is strongly dominated by secondary nucleation, i.e. the formation of new nuclei on the surface of existing fibrils. The critical nucleus consists of a single IAPP molecule, and anionic lipids catalyze both primary and secondary nucleation, but not elongation. The fact that anionic lipids promote secondary nucleation implies that these events take place at the interface between the membrane and existing fibrils, demonstrating that fibril growth occurs at least to some extent on the membrane surface. These new insights into the mechanism of IAPP aggregation on membranes may help to understand IAPP toxicity and will be important for the development of therapeutics to prevent {beta}-cell death in type II diabetes. Co-option of immune effectors by the hormonal signalling system triggering metamorphosis in Drosophila melanogaster Insect metamorphosis is regulated by the production, secretion and degradation of two peripheral hormones: 20-hydroxyecdysone (ecdysone) and juvenile hormone (JH). In addition to their roles in developmental regulation, increasing evidence suggests that these hormones are involved in innate immunity processes, such as phagocytosis and the induction of antimicrobial peptide (AMP) production. AMP regulation includes systemic responses as well as local responses at surface epithelia that contact with the external environment. At pupariation, Drosophila melanogaster increases dramatically the expression of three AMP genes, drosomycin (drs), drosomycin-like 2 (drsl2) and drosomycin-like 5 (drsl5). We show that the systemic action of drs at pupariation is dependent on ecdysone signalling in the fat body and operates via the ecdysone downstream target, Broad-Z4. In parallel, ecdysone also regulates local responses, specifically through the activation of drsl2 expression in the gut. Finally, we confirm the relevance of this ecdysone dependent AMP expression for the control of bacterial persistence by showing that flies lacking drs expression in the fat body have higher bacterial persistence over metamorphosis. Together, our data emphasize the importance of the association between ecdysone signalling and immunity using in vivo studies and establish a new role for ecdysone at pupariation, which impacts developmental success by regulating the immune system in a stage-dependent manner. We speculate that this co-option of immune effectors by the hormonal system may constitute a pre-emptive mechanism to control bacterial numbers in the pupa at the core of metamorphosis evolution. Cell-specific Bioorthogonal Tagging of Glycoproteins Altered glycosylation is an undisputed corollary of cancer development. Understanding these alterations is paramount but hampered by limitations underlying cellular model systems. For instance, the intricate interactions between tumour and host cannot be adequately recapitulated in monoculture of tumour-derived cell lines. More complex co-culture models usually rely on sorting procedures for proteome analyses and rarely capture the details of protein glycosylation. Here, we report a strategy termed Bio-Orthogonal Cell line-specific Tagging of Glycoproteins (BOCTAG). Cells are equipped by transfection with an artificial biosynthetic pathway that transforms bioorthogonally tagged sugars into the corresponding nucleotide-sugars. Only transfected cells incorporate bioorthogonal tags into glycoproteins in the presence of non-transfected cells. We employ BOCTAG as an imaging technique and to annotate cell-specific glycosylation sites in mass spectrometry-glycoproteomics. We demonstrate application in co-culture and mouse models, allowing for profiling of the glycoproteome as an important modulator of cellular function. A protease-activatable luminescent biosensor andreporter cell line for authentic SARS-CoV-2 infection Efforts to define serological correlates of protection against COVID-19 have been hampered by the lack of a simple, scalable, standardised assay for SARS-CoV-2 infection and antibody neutralisation. Plaque assays remain the gold standard, but are impractical for high-throughput screening. In this study, we show that expression of viral proteases may be used to quantitate infected cells. Our assays exploit the cleavage of specific oligopeptide linkers, leading to the activation of cell-based optical biosensors. First, we characterise these biosensors using recombinant SARS-CoV-2 proteases. Next, we confirm their ability to detect viral protease expression during replication of authentic virus. Finally, we generate reporter cells stably expressing an optimised luciferase-based biosensor, enabling viral infection to be measured within 24 h in a 96- or 384-well plate format, including variants of concern. We have therefore developed a luminescent SARS-CoV-2 reporter cell line, and demonstrated its utility for the relative quantitation of infectious virus and titration of neutralising antibodies.

Author summaryTechniques for measuring infection with SARS-CoV-2 in the laboratory are laborious and time-consuming, and different laboratories use different approaches. There is therefore no generally agreed way to quantitate neutralising antibodies against SARS-CoV-2, which block infection with the virus and protect people from COVID-19. In this study, we describe a new way to measure SARS-CoV-2 infection, which is much simpler and faster than existing methods. It relies on the production of a specific protease enzyme by the virus, which is able to cleave and activate an engineered protein biosensor in infected cells. This biosensor emits light in the presence of viral infection, and the amount of light released is used as a readout for the amount of infectious SARS-CoV-2 present. The signal is very sensitive, so the number of infected cells required is very small, and the method can be scaled-up to test many samples at once. In particular, we demonstrate how it can be used to detect different variants of SARS-CoV-2, and quantitate neutralising antibodies against these viruses. 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.