Scalability of large neural network simulations via activity tracking with time asynchrony and procedural connectivity We present a new algorithm to efficiently simulate random models of large neural networks satisfying the property of time asynchrony. The model parameters (average firing rate, number of neurons, synaptic connection probability, and postsynaptic duration) are of the order of magnitude of a small mammalian brain, or of human brain areas. Through the use of activity tracking and procedural connectivity (dynamical regeneration of synapses), both computational and memory complexities of this algorithm are proved to be theoretically linear with the number of neurons. These results are experimentally validated by sequential simulations of millions of neurons and billions of synapses running in few minutes using a single thread of an equivalent desktop computer.
Disruption of miR-18a alters proliferation, photoreceptor replacement kinetics, inflammatory signaling and microglia/macrophage numbers during retinal regeneration in zebrafish In mammals, photoreceptor loss causes permanent blindness, but in zebrafish (Danio rerio), photoreceptor loss reprograms Muller glia to function as stem cells, producing progenitors that regenerate photoreceptors. MicroRNAs (miRNAs) regulate CNS neurogenesis, but the roles of miRNAs in injury-induced neuronal regeneration are largely unknown. In the embryonic zebrafish retina, miR-18a regulates photoreceptor differentiation. The purpose of the current study was to determine, in zebrafish, the function of miR-18a during injury-induced photoreceptor regeneration. RT-qPCR, in-situ hybridization and immunohistochemistry showed that miR-18a expression increases throughout the retina between 1 and 5 days post-injury (dpi). To test miR-18a function during photoreceptor regeneration, we used homozygous miR-18a mutants (miR-18ami5012), and knocked down miR-18a with morpholino oligonucleotides. During photoreceptor regeneration, miR-18ami5012 retinas have fewer mature photoreceptors than WT at 7 and 10 dpi, but there is no difference at 14 dpi, indicating that photoreceptor regeneration is delayed. Labeling dividing cells with 5-bromo-2'-deoxyuridine (BrdU) showed that at 7 and 10 dpi, there are excess dividing progenitors in both mutants and morphants, indicating that miR-18a negatively regulates injury-induced proliferation. Tracing 5-ethynyl-2'-deoxyuridine (EdU) and BrdU-labeled cells showed that in miR-18ami5012 retinas excess progenitors migrate to other retinal layers in addition to the photoreceptor layer. Inflammation is critical for photoreceptor regeneration, and RT-qPCR showed that in miR-18ami5012 retinas, inflammatory gene expression and microglia activation are prolonged. Suppressing inflammation with dexamethasone rescues the miR-18ami5012 phenotype. Together, these data show that in the injured zebrafish retina, disruption of miR-18a alters proliferation, inflammation, the microglia/macrophage response, and the timing of photoreceptor regeneration.
Oct1 recruits the histone lysine demethylase Utx to promote lineage specification The pathways used by cells to transition between undifferentiated, pluripotent state and tissue-specific states are incompletely understood. Here we show that the widely-expressed transcription factor Oct1/Pou2f1 activates silent, developmental lineage-appropriate genes to "canalize" developmental progression. Using Oct1 inducible knockout embryonic stem cells, we show that that Oct1 deficiency impairs mesodermal and terminal muscle differentiation in a manner that can be rescued by Oct1 retroviral expression. Using bulk RNA-seq, we show that mesoderm-specific genes are not correctly induced early in the differentiation timecourse. Single-cell gene expression profiling reveals that Oct1-deficient cells lose coherence in temporal induction of lineage programs, and show inappropriate developmental lineage branching resulting in poorly differentiated cells state with epithelial characteristics and hallmarks of oxidative stress. In embryonic stem cells, Oct1 co-binds with Oct4 to genes critical for mesoderm induction. The Utx/Kdm6a histone lysine demethylase also binds to many of these genes, and using a prototypic Pax3 gene we show that Oct1 recruits Utx to remove inhibitory H3K27me3 marks and activate expression. The specificity of the ubiquitous Oct1 protein for mesodermal genes can be explained by cooperative interactions with lineage-driving Smad transcription factors, as we show that Smad and Oct binding sites frequently coexist mesoderm-specific genes, that Oct1 and Smad3 interact, and that the sites and factors act cooperatively at the Myog enhancer. Overall, these results identify Oct1 as a key mediator of the induction of mesoderm lineage-specific genes. Bacteria swarm-mediated phage transportation disrupts a biofilm inherently protected from phage penetration The treatment of chronic bacterial infections by phages has shown promise in combating antimicrobial resistance. A typical phage particle is at least an order of magnitude larger than an antibiotic molecule. Hence, phages diffuse slower than antibiotics, and can also get trapped in the polymeric mesh of biofilm matrix. By tracking fluorescently labeled lambda phages that do not infect Capnocytophaga gingivalis, a bacterium abundant in the human oral microbiota, we demonstrate active phage transportation by a C. gingivalis swarm. As a result, the rate of disruption of the prey of lambda phage i.e., an Escherichia coli colony, increases 10 times. C. gingivalis drills tunnels within a curli fiber containing E. coli biofilm and increase the efficiency of phage penetration. This provides evidence for a novel mechanism of phage-bacterial warfare. Using machine learning to design adeno-associated virus capsids with high likelihood of viral assembly We study the application of machine learning in designing adeno-associated virus (AAV) capsid sequences with high likelihood of viral assembly, i.e. capsid viability. Specifically, we design and implement Origami, a model-based optimization algorithm, to identify highly viable capsid sequences within the vast space of 2033 possibilities. Our evaluation shows that Origami performs well in terms of optimality and diversity of model-designed sequences. Moreover, these sequences are ranked according to their viability score. This helps designing experiments given budget constraint. A potent human monoclonal antibody with pan-neutralizing activities directly dislocates S trimer of SARS-CoV-2 through binding both up and down forms of RBD The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a global pandemic of novel corona virus disease (COVID-19). The neutralizing monoclonal antibodies (mAbs) targeting the receptor binding domain (RBD) of SARS-CoV-2 are among the most promising strategies to prevent and treat COVID-19. However, SARS-CoV-2 variants of concern (VOCs) profoundly reduced the efficacies of most of mAbs and vaccines approved for clinical use. Herein, we demonstrated mAb 35B5 efficiently neutralizes both wild-type (WT) SARS-CoV-2 and VOCs, including B.1.617.2 (delta) variant, in vitro and in vivo. Cryo-electron microscopy (cryo-EM) revealed that 35B5 neutralizes SARS-CoV-2 by targeting a unique epitope that avoids the prevailing mutation sites on RBD identified in circulating VOCs, providing the molecular basis for its pan-neutralizing efficacy. The 35B5-binding epitope could also be exploited for the rational design of a universal SARS-CoV-2 vaccine. Raman Signal Denoising using Fully Convolutional Encoder Decoder Network Raman spectroscopy is a vibrational method that gives molecular information rapidly and non-invasively. Despite its advantages, the weak intensity of Raman spectroscopy leads to low-quality signals, particularly with tissue samples. The requirement of high exposure times makes Raman a time-consuming process and diminishes its non-invasive property while studying living tissues. Novel denoising techniques using convolutional neural networks (CNN) have achieved remarkable results in image processing. Here, we propose a similar approach for noise reduction for the Raman spectra acquired with 10x lower exposure times. In this work, we developed fully convolutional encoder-decoder architecture (FCED) and trained them with noisy Raman signals. The results demonstrate that our model is superior (p-value < 0.0001) to the conventional denoising techniques such as the Savitzky-Golay filter and wavelet denoising. Improvement in the signal-to-noise ratio values ranges from 20% to 80%, depending on the initial signal-to-noise ratio. Thus, we proved that tissue analysis could be done in a shorter time without any need for instrumental enhancement. Spatial regulation of the anaphase-promoting complex/cyclosome by a microtubule adaptor The anaphase-promoting complex/cyclosome (APC/C) coordinates advancement through mitosis via temporally controlled polyubiquitination of effector proteins. Despite the long-appreciated spatial organization of key events in mitosis mediated largely by cytoskeletal networks, the spatial regulation of APC/C, the major mitotic E3 ligase, is poorly understood. Here, we describe a microtubule-resident protein, PLEKHA5, as an interactor of APC/C and spatial regulator of its activity in mitosis. PLEKHA5 knockdown delayed mitotic progression, causing accumulation of APC/C substrates dependent upon the PLEKHA5-APC/C interaction. To assess the spatial importance of this interaction, we developed microtubule-localized proximity biotinylation tools, which revealed that depletion of PLEKHA5 decreased the extent of APC/C association with the microtubule cytoskeleton. This decreased APC/C microtubule-localization in turn prevented efficient loading of APC/C with its co-activator CDC20, leading to defects in E3 ligase catalytic activity. We propose that PLEKHA5 functions as an adaptor of APC/C that promotes its subcellular localization to microtubules and facilitates its activation by CDC20, thus ensuring the timely turnover of key mitotic APC/C substrates and proper progression through mitosis. A theoretical model of mitochondrial ATP synthase deficiencies. The role of mitochondrial carriers. The m.8993T>G mutation of the mitochondrial MT-ATP6 gene is associated with NARP syndrome (Neuropathy, Ataxia and Retinitis Pigmentosa). The equivalent point mutation introduced in yeast Saccharomyces cerevisiae mitochondrial DNA considerably reduced the activity of ATP synthase and of cytochrome-c-oxidase preventing yeast growth on oxidative substrates. The over-expression of the mitochondrial oxodicarboxylate carrier (Odc1p) is able to rescue the growth on oxidative substrate by increasing the substrate-level phosphorylation of ADP coupled to conversion of -ketoglutarate (AKG) into succinate with an increase in Complex IV activity. Previous studies showed that equivalent point mutations in ATP synthase behave similarly and can be rescued by Odc1p overexpression and/or uncoupler of OXPHOS from ATP synthesis. In order to better understand the mechanism of ATP synthase mutation bypass, we developed a core model of mitochondrial metabolism based on AKG as respiratory substrate. We describe the different possible metabolite output and the ATP/O ratio values as a function of ATP synthase inhibition. Targeting MLL methyltransferases enhance the anti-tumor effects of PI3K inhibition in hormone receptor-positive breast cancer The high frequency of aberrant PI3K pathway activation in hormone receptor-positive (HR+) breast cancer has led to the development, clinical testing, and approval of the p110-selective PI3K inhibitor alpelisib. The limited clinical efficacy of alpelisib and other PI3K inhibitors is partially attributed to the functional antagonism between PI3K and estrogen receptor (ER) signaling, which is mitigated via combined PI3K inhibition and endocrine therapy. We and others have previously demonstrated a chromatin-associated mechanism by which PI3K supports cancer development and antagonizes ER signaling through the modulation of H3K4 methylation. Here we show that inhibition of the H3K4 histone methyltransferase MLL1 in combination with PI3K inhibition impairs HR+ breast cancer clonogenicity and cell proliferation. While combined PI3K/MLL1 inhibition reduces AKT effector signaling and H3K4 methylation, MLL1 inhibition increases PI3K effector signaling and upregulates the expression of receptor tyrosine kinase signaling cascades upstream of and including AKT. These data reveal a feedback loop between MLL1 and AKT in which MLL1 inhibition reactivates AKT. We additionally show that combined PI3K and MLL1 inhibition synergizes to cause cell death in in vitro and in vivo models of HR+ breast cancer, which is enhanced by the additional genetic ablation of the H3K4 methyltransferase and AKT target MLL4. Together, our data provide evidence of a feedback mechanism connecting histone methylation with AKT and may support the preclinical development and testing of pan-MLL inhibitors.

Significance StatementPharmacological inhibition of PI3K provides limited efficacy in PIK3CA-mutated, HR+ breast cancers. Here the authors leverage PI3K/AKT-driven chromatin modification to identify MLL histone methyltransferases as a therapeutic target. Dual PI3K and MLL inhibition synergize to reduce clonogenicity and cell proliferation while enhancing apoptosis in in vitro models, and induces tumor regression in xenograft models of PI3K-activated, HR+ breast cancer. Furthermore, MLL1 inhibition reveals a feedback loop leading to AKT hyperactivation, which is relieved with combined PI3K/MLL inhibition. These findings demonstrate the utility of MLL inhibitors for the treatment of some solid cancers, as patients with HR+ breast cancer characterized by PIK3CA mutation may derive clinical benefit from combined PI3K/MLL inhibition. 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. XRN2-mediated regulation of let-7 microRNAs is of critical pathophysiological significance in Humans MicroRNAs are critical regulators of diverse developmental and physiological processes, and they themselves get regulated both at the level of biogenesis and turnover. We demonstrate that the ribonuclease XRN2 can degrade the mature forms of certain let-7 family members in multiple human cancer cell lines, without affecting their precursors. XRN2 depletion results in a reduction in the expression of a number of oncogenes, and diminishes the proliferative and metastatic potential of cancer cells. The clinical relevance of these observations is also verified in tumour transcriptomics data from public RNA-sequencing datasets, where XRN2 mRNA expression is inversely correlated with the levels of a large number of miRNAs, including let-7 members, and high XRN2 mRNA levels are associated with poor survival in hepatocellular carcinoma, lung adenocarcinoma, and glioblastoma. We demonstrate that the miRNA is released by an as-yet unidentified proteinaceous  miRNA release factor from the grasp of Argonaute before its degradation. Our analyses of the patient-derived transcriptomics data also show that XRN2, via its regulation of let-7, affects multiple pathways in a consistent manner across epithelial and glial cell lineages, and thus, is of critical pathophysiological significance. 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. The initial response of females towards congeneric males matches the propensity to hybridise in Ophthalmotilapia. AO_SCPLOWBSTRACTC_SCPLOWCichlid radiations often harbour closely related species with overlapping niches and distribution ranges. Such species sometimes hybridise in nature, which raises the question how can they coexist. This also holds for the Tanganyika mouthbrooders Ophthalmotilapia ventralis and O. nasuta. Earlier studies found indications of asymmetrical hybridisation with females of O. ventralis accepting males of O. nasuta, but not the other way around. We hypothesised that this was due to differences in the capacity for species recognition. Given the higher propensity of O. ventralis females towards hybridisation, we expect a reduced ability for species recognition in O. ventralis females, compared to O. nasuta females. We staged two experiments, one focusing on 22 female O. nasuta and one on 21 female O. ventralis. These fish were placed in one half of a tank and briefly exposed to a conspecific or a heterospecific male, a conspecific female, or nothing (control). Female response was evaluated by scoring six tracking parameters and by noting the occurrence of ten discrete behaviours before and during the encounter. Females always responded to the presence of another fish by approaching it. Remarkably, for both O. nasuta and O. ventralis, we did not find a different response between encounters with conspecific males and females. However, in agreement with our hypothesis, O. nasuta females behaved differently towards conspecific or heterospecific males, whereas O. ventralis females did not. When presented with a heterospecific male, O. nasuta females performed a lower number of  ram behaviours. Additionally, they never displayed the  flee behaviour, a component of the species mating repertoire that was seen in all but one of the presentations with a conspecific male. Our findings show that differences in species recognition at first encounter predict to a large degree the outcome of the mating process, even in the absence of mating behaviour. Multi-omic analysis of the Arabidopsis clock activator mutant rve 4 6 8 reveals connections to carbohydrate metabolism and proteasome regulation Plants are able to sense changes in their light environments, such as the onset of day and night, as well as anticipate these changes in order to adapt and survive. Central to this ability is the plant circadian clock, a molecular circuit that precisely orchestrates plant cell processes over the course of a day. REVEILLE proteins (RVEs) are recently discovered members of the plant circadian circuitry that activate the evening complex and PRR genes to maintain regular circadian oscillation. The RVE 8 protein and its two homologs, RVE 4 and 6, have been shown to limit the length of the circadian period, with rve 4 6 8 triple-knockout plants possessing an elongated period along with increased leaf surface area, biomass, cell size and delayed flowering relative to wild-type Col-0 plants. Here, using a multi-omics approach consisting of phenomics, transcriptomics, proteomics, and metabolomics we draw novel connections between RVE8-like proteins and a number of core plant cell processes. In particular, we reveal that loss of RVE8-like proteins results in altered carbohydrate, organic acid and lipid metabolism, including a starch excess phenotype at dawn. We further demonstrate that rve 4 6 8 plants have lower levels of 20S proteasome subunits and possess significantly reduced proteasome activity, potentially explaining the increase in cell-size observed in RVE8-like mutants. Overall, this robust, multi-omic dataset, provides substantial new insights into the far reaching impact RVE8-like proteins have on the diel plant cell environment. Structure-Functional-Selectivity Relationship Studies on A-86929 Analogs and Small biaryl Fragments toward Discovery of biased D1 agonists Dopamine regulates normal functions such as movement, reinforcement learning, and cognition, and its dysfunction has been implicated in multiple psychiatric and neurological disorders. Dopamine acts through the D1- (D1R and D5R) and D2-class (D2R, D3R and D4R) of seven transmembrane receptors, and activates both G-protein- and {beta}-arrestin-dependent signaling pathways, to mediate its physiological effects. Current dopamine receptor-based therapies are used to ameliorate motor deficits in Parkinsons disease, or as antipsychotic medications for schizophrenia. These drugs show efficacy for ameliorating only some symptoms caused by dopamine dysfunction and are plagued by debilitating side-effects. Studies in primates and rodents have shown that shifting the balance of dopamine receptor signaling towards the arrestin pathway can be beneficial for inducing normal movement, while reducing motor side-effects such as dyskinesias, and can be efficacious at enhancing cognitive function compared to balanced agonists. Several structure-activity-relationship (SAR) studies have embarked on discovering {beta}-arrestin-biased dopamine agonists, focused on D2 partial agonists, non-catechol D1 agonists, and mixed D1/D2R dopamine receptor agonists. Here, we describe an SAR study to identify novel D1R {beta}-arrestin biased ligands using A-86929, a high-affinity D1R catechol agonist, as a core scaffold. Previously described and novel analogs of A-86929 were synthesized and screened in vitro for structure-functional-selectivity relationships (SFSR) studies to identify chemical motifs responsible for {beta}-arrestin biased activity at both D1 and D2Rs. Most of the A-86929 analogs screened were G protein biased but none of them were exclusively arrestin-biased. Additionally, various catechol aryl fragments were designed and synthesized. Other compounds surveyed included hydroxyl and chloro analogs of dopamine to test for hydrogen bonding and ionic interactions. Some of these small molecular probes displayed weak bias towards the {beta}-arrestin pathway. Continued in-depth SFSR studies informed by structure determination, molecular modeling, and mutagenesis studies will facilitate discovery of potent and efficacious arrestin-biased dopamine receptor ligands. Elucidation of SARS-CoV-2 budding mechanisms through molecular dynamics simulations of M and E protein complexes SARS-CoV-2 and other coronaviruses pose major threats to global health, yet computational efforts to understand them have largely overlooked the process of budding, a key part of the coronavirus life cycle. When expressed together, coronavirus M and E proteins are sufficient to facilitate budding into the ER-Golgi intermediate compartment (ERGIC). To help elucidate budding, we ran atomistic molecular dynamics (MD) simulations using the Feig laboratorys refined structural models of the SARS-CoV-2 M protein dimer and E protein pentamer. Our MD simulations consisted of M protein dimers and E protein pentamers in patches of membrane. By examining where these proteins induced membrane curvature in silico, we obtained insights around how the budding process may occur. Multiple M protein dimers acted together to induce global membrane curvature through protein-lipid interactions while E protein pentamers kept the membrane planar. These results could eventually help guide development of antiviral therapeutics which inhibit coronavirus budding. BlenderPhotonics - an integrated open-source software environment for 3-D meshing and photon simulations in complex tissues SignificanceRapid advances in biophotonics techniques require quantitative, model-based computational approaches to obtain functional and structural information from increasingly complex and multi-scaled anatomies. The lack of efficient tools to accurately model tissue structures and subsequently perform quantitative multi-physics modeling greatly impedes the clinical translation of these modalities.

AimWhile the mesh-based Monte Carlo (MMC) method expands our capabilities in simulating complex tissues by using tetrahedral meshes, the generation of such domains often requires specialized meshing tools such as Iso2Mesh. Creating a simplified and intuitive interface for tissue anatomical modeling and optical simulations is essential towards making these advanced modeling techniques broadly accessible to the user community.

ApproachWe responded to the above challenge by combining the powerful, open-source 3-D modeling software, Blender, with state-of-the-art 3-D mesh generation and MC simulation tools, utilizing the interactive graphical user interface (GUI) in Blender as the front-end to allow users to create complex tissue mesh models, and subsequently launch MMC light simulations.

ResultsHere we present a tutorial to our newly developed Python-based Blender add-on - BlenderPhotonics - to interface with Iso2Mesh and MMC, allowing users to create, configure and refine complex simulation domains and run hardware-accelerated 3-D light simulations with only a few clicks. We provide a comprehensive introduction to this new tool and walk readers through 5 examples, ranging from simple shapes to sophisticated realistic tissue models.

ConclusionBlenderPhotonics is user-friendly and open-source, leveraging the vastly rich ecosystem of Blender. It wraps advanced modeling capabilities within an easy-to-use and interactive interface. The latest software can be downloaded at http://mcx.space/bp. 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. Re-convolving the compositional landscape of primary and recurrent glioblastoma reveals prognostic andtargetable tissue states Glioblastoma is an aggressive diffusely infiltrating neoplasm that spreads beyond surgical resection margins, where it intermingles with non-neoplastic brain cells. This complex microenvironment harbouring infiltrating glioma and non-neoplastic brain cells is the origin of tumor recurrence. Thus, understanding the cellular and molecular features of the glioma microenvironment is therapeutically and prognostically important. Here, we used single-nucleus RNA sequencing (snRNAseq) of primary and recurrent glioma to define three compositional  tissue-states rooted in cohabitation of cell-types and transcriptional states. These comprise a state featuring A) abundant normal brain microenvironment cells, B) reactive/inflamed brain microenvironment, and C) cellular/proliferative tumor. All these states exhibited variable degrees of infiltration by glioma cells. Spatial transcriptomics confirmed that the cell-types and transcriptomics states which compositionally cohabitate tissue states also colocalize in space. Tissue states are clinically significant because they correlate with radiographic, histopathologic, and prognostic features. We found that enrichment of tissue state B signature correlated with shorter survival. Importantly, we found that our compositionally defined tissue states are enriched in distinct metabolic pathways. One such pathway is fatty acid biosynthesis, which was enriched in tissue state B - a state enriched in recurrent glioblastoma and composed of astrocyte-like/mesenchymal glioma cells, reactive astrocytes resembling those seen in neurodegeneration, and monocyte-like myeloid cells. We showed that treating acute slices of GBM with a fatty acid synthesis inhibitor is sufficient to deplete the transcriptional signature of tissue state B. Our findings define a novel compositional approach to the glioma-infiltrated tissue which allows us to discover prognostic and targetable features, paving the way to new mechanistic and therapeutic discoveries. Admixture of evolutionary rates across a hybrid zone Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii, we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates, while genomic regions with greater reproductive isolation evolve at divergent rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and that maintenance of divergent substitution rates can be predicted quantitatively from relative sequence divergence (FST) between the hybridizing species at equilibrium. A corollary is that divergent rates will be maintained in regions linked to barrier loci. Overall, we demonstrate that reproductive isolation maintains not only the final outcome of genomic differentiation, but also the rate at which differentiation accumulates. This new information also suggests that the separation of evolutionary rates co-localizes with the separation of gene pools between genomes of incipient species. Mutations in coral soma and sperm imply lifelong stem cell renewal and cell lineage selection In many animals, the germline differentiates early in embryogenesis, so only mutations that accumulate in germ cells are inherited by offspring1. Exceptions to this developmental process may indicate that other mechanisms have evolved to limit the effects of deleterious mutation accumulation2. Stony corals are animals that can live for hundreds of years3 and have long been thought to produce gametes from somatic tissue4. To clarify conflicting evidence about germline-soma distinction in corals, we sequenced high coverage, full genomes with technical replicates for parent coral branches and their sperm pools. We identified post-embryonic single nucleotide variants (SNVs) unique to each parent branch, then checked if each SNV was shared by the respective sperm pool: 26% of post-embryonic SNVs were shared by the sperm but 74% were not. We also identified germline SNVs, those that were present in the sperm but not in the parent. These data suggest that self-renewing stem cells in corals differentiate into germ and soma throughout the adult life of the colony, with SNV rates and patterns differing markedly in stem, soma, and germ lineages. In addition to informing the important place in the evolutionary spectrum from non-Weismannian to Weismmanian animals that corals occupy, these insights inform how corals may generate adaptive diversity necessary in the face of global climate change. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Sex determination gene transformer regulates the male-female difference in Drosophila fat storage via the adipokinetic hormone pathway Sex differences in whole-body fat storage exist in many species. For example, Drosophila females store more fat than males. Yet, the mechanisms underlying this sex difference in fat storage remain incompletely understood. Here, we identify a key role for sex determination gene transformer (tra) in regulating the male-female difference in fat storage. Normally, a functional Tra protein is present only in females, where it promotes female sexual development. We show that loss of Tra in females reduced whole-body fat storage, whereas gain of Tra in males augmented fat storage. Tras role in promoting fat storage was largely due to its function in neurons, specifically the Adipokinetic hormone (Akh)-producing cells (APCs). Our analysis of Akh pathway regulation revealed a male bias in APC activity and Akh pathway function, where this sex-biased regulation influenced the sex difference in fat storage by limiting triglyceride accumulation in males. Importantly, Tra loss in females increased Akh pathway activity, and genetically manipulating the Akh pathway rescued Tra-dependent effects on fat storage. This identifies sex-specific regulation of Akh as one mechanism underlying the male-female difference in whole-body triglyceride levels, and provides important insight into the conserved mechanisms underlying sexual dimorphism in whole-body fat storage. Migrators Within Migrators: Exploring Transposable Element Dynamics in the Monarch Butterfly, Danaus plexippus BackgroundLepidoptera (butterflies and moths) are an important model system in ecology and evolution. A high-quality chromosomal genome assembly is available for the monarch butterfly (Danaus plexippus), but it lacks an in-depth transposable element (TE) annotation, presenting an opportunity to explore monarch TE dynamics and the impact of TEs on shaping the monarch genome.

ResultsWe find 6.21% of the monarch genome is comprised of TEs, a reduction of 6.85% compared to the original TE annotation performed on the draft genome assembly. Monarch TE content is low compared to two closely related species with available genomes, Danaus chrysippus (33.97% TE) and Danaus melanippus (11.87% TE). The biggest TE contributions to genome size in the monarch are LINEs and Penelope-like elements, and three newly identified families, r2-hero_dPle (LINE), penelope-1_dPle (Penelope-like), and hase2-1_dPle (SINE), collectively contribute 34.92% of total TE content. We find evidence of recent TE activity, with two novel Tc1 families rapidly expanding over recent timescales (tc1-1_dPle, tc1- 2_dPle). LINE fragments show signatures of genomic deletions indicating a high rate of TE turnover. We investigate associations between TEs and wing colouration and immune genes and identify a three-fold increase in TE content around immune genes compared to other host genes.

ConclusionsWe provide a detailed TE annotation and analysis for the monarch genome, revealing a considerably smaller TE contribution to genome content compared to two closely related Danaus species with available genome assemblies. We identify highly successful novel DNA TE families rapidly expanding over recent timescales, and ongoing signatures of both TE expansion and removal highlight the dynamic nature of repeat content in the monarch genome. Our findings also suggest that insect immune genes are promising candidates for future interrogation of TE-mediated host adaptation. Quantifying isolation-by-resistance and connectivity in dendritic ecological networks Statistical methods that quantify individual movements have limited applicability within dendritic ecological networks (DENs), to include riverscapes. We promote an approach herein by deriving an analytical framework (=DO_SCPLOWISTC_SCPLOWNO_SCPLOWETC_SCPLOW and RO_SCPLOWESISTC_SCPLOWNO_SCPLOWETC_SCPLOW) that more appropriately quantifies differentiation by meshing individual DEN segments with fitted distances and annotated environmental features. We first explored the spatial arrangement of diversity in a river network by fitting pairwise distances against graph edges. We then derived effective resistance models (e.g., as a composite of an array of possible environmental covariates) by using a genetic algorithm for heuristic model selection, invoking circuit theory (CO_SCPLOWIRCUITC_SCPLOWSO_SCPLOWCAPEC_SCPLOW) to form a graph-based complement to RO_SCPLOWESISTANCEC_SCPLOWGA (a landscape ecology package). Although numerous pairwise distance metrics can be employed, we utilized genetic distances encompassing 13,218 ddRAD loci from N=762 Speckled Dace (Cyprinidae: Rhinichthys osculus) sampled at 78 Colorado River localities ({micro}=9.77/site). When distances for each locus were fitted into the network, we found intraspecific divergences due to stream reaches were greater than expected when distances alone. In fitting model-averaged  effective resistance models to the same network, we showed the manner by which a host of anthropogenic and environmental variables contribute to patterns of connectivity underlying riverscape genetic differentiation. Our framework represents a contemporary approach to deriving metapopulation/ metacommunity structure within DENs, in that it allows: (a) The extension of projections into unsampled temporal/spatial components; (b) Comparisons to be made among species and/or drainages, both of which benefit multi-species management; and (c) the quantification of locus-specific patterns, from which adaptive hypotheses could then be derived. O_QD"The intimate relation of hydrological conditions to fish life is proverbial; to be hopelessly or dangerously maladjusted to ones environment is to be a  fish out of water" (Hubbs 1941).

C_QD