TRAWLING: a Transcriptome Reference AWare of spLIciNG events. Alternative splicing is critical for human gene expression regulation and plays an important role in multiple human diseases. In this context, RNA sequencing has emerged as powerful approach to detect alternative splicing events.

In parallel, fast alignment-free methods have emerged as a viable alternative to quantify gene and transcript level abundance from RNAseq data. However, the ability to detect differential splicing events is dependent on the annotation of the transcript reference provided by the user.

Here, we introduce a new reference transcriptome aware of splicing events, TRAWLING, which simplifies the detection of aberrant splicing events in a fast and simple way. In addition, we evaluate the performances and the benefits of aligning transcriptome data to TRAWLING using three different RNA sequencing datasets: whole transcriptome sequencing, single cell RNA sequencing and Digital RNA with pertUrbation of Genes.

Collectively, our comprehensive evaluation underlines the value of using TRAWLING in transcriptomic data analysis.

Availability and implementationOur code is available at https://github.com/Novartis/TRAWLING
Estimating cell-type-specific gene co-expression networks from bulk gene expression data with an application to Alzheimer's disease Inferring and characterizing gene co-expression networks has led to important insights on the molecular mechanisms of complex diseases. Most co-expression analyses to date have been performed on gene expression data collected from bulk tissues with different cell type compositions across samples. As a result, the co-expression estimates only offer an aggregate view of the underlying gene regulations and can be confounded by heterogeneity in cell type compositions, failing to reveal gene coordination that may be distinct across different cell types. In this paper, we describe a flexible framework for estimating cell-type-specific gene co-expression networks from bulk sample data, without making specific assumptions on the distributions of gene expression profiles in different cell types. We develop a novel sparse least squares estimator, referred to as CSNet, that is efficient to implement and has good theoretical properties. Using CSNet, we analyzed the bulk gene expression data from a cohort study on Alzheimers disease and identified previously unknown cell-type-specific co-expressions among Alzheimers disease risk genes, suggesting cell-type-specific disease pathology for Alzheimers disease.
Autophagosomes fuse to phagosomes and facilitate the degradation of apoptotic cells in Caenorhabditis elegans Autophagosomes are double-membrane intracellular vesicles that degrade protein aggregates, intracellular organelles, and other cellular components. During the development of the nematode Caenorhabditis elegans, many somatic and germ cells undergo apoptosis. These cells are engulfed and degraded by their neighboring cells. We discovered a novel role of autophagosomes in facilitating the degradation of apoptotic cells using a real-time imaging technique. Specifically, the double-membrane autophagosomes in engulfing cells are recruited to the surfaces of phagosomes containing apoptotic cells and subsequently fuse to phagosomes, allowing the inner vesicle to enter the phagosomal lumen. Mutants defective in the production of autophagosomes display significant defects in the degradation of apoptotic cells, demonstrating the importance of autophagosomes to this process. The signaling pathway led by the phagocytic receptor CED-1, the adaptor protein CED-6, and the large GTPase dynamin (DYN-1) promotes the recruitment of autophagosomes to phagosomes. Moreover, the subsequent fusion of autophagosomes with phagosomes requires the functions of the small GTPase RAB-7 and the HOPS complex components. Further observations suggest that autophagosomes provide apoptotic cell-degradation activities in addition to and in parallel of lysosomes. Our findings reveal that, unlike the single-membrane, LC3-associated phagocytosis (LAP) vesicles reported for mammalian phagocytes, the canonical double-membrane autophagosomes facilitate the clearance of C. elegans apoptotic cells. These findings add autophagosomes to the collection of intracellular organelles that contribute to phagosome maturation, identify novel crosstalk between the autophagy and phagosome maturation pathways, and discover the upstream signaling molecules that initiate this crosstalk. Flavobacterium johnsoniae Tyrosine Ammonia Lyase (FjTAL) in-silico Structure Prediction and Molecular Docking to L-Tyrosine, p-Coumaric Acid (pHCA) and Caffeic Acid Responsible for the conversion of L-tyrosine to p-coumaric acid in Flavobacterium johnsoniae, FjTAL has drawn the attention of many biochemical engineers who wish to carry out a sustainable biosynthetic scheme for the production of aromatic compounds. In this study, with the aid of various computational tools, the secondary and tertiary structures of FjTAL have been predicted. The results suggest that FjTAL forms a homo-tetramer when active as a cytosolic enzyme and it is mostly consisted of alpha helices. With the aid of molecular docking, one can hypothesize that FjTAL is likely to bind to L-tyrosine, p-coumaric acid, and caffeic acid with a similar molecular mechanism and thus, p-coumaric acid and caffeic acid may exhibit a negative feedback response toward the enzyme and inhibit its activity competitively. Two distinct binding pockets have been discovered, one of which contains highly conserved residues among several species. The residues which form the prosthetic group 3,5-dihydro-5-methylidene-4H-imidazol-4-one (MIO) also emerge in the evolutionary conserved binding pocket. The other discovered cavity, could either be a second binding site for the ligands or simply an artifact of the molecular docking task. Nuclear Progestin Receptor Mediated Linkage of Blood Coagulation and Ovulation Ovulation is a dramatic remodeling process that includes rupture of blood capillaries and clotting, but coagulation is not thought to directly regulate this process. Herein, we report remarkable increases of coagulation factors V (f5, [~]3145-fold) and tissue factor (f3a, [~]120-fold) in zebrafish ovarian follicle cells during ovulation. This increase was mediated through the nuclear progestin receptor (Pgr), which is essential for ovulation in zebrafish, and was totally abolished in ovarian follicular cells from pgr-/- mutants. In addition, promoter activities of f5 and f3a were significantly enhanced by progestin (DHP) via Pgr. Similar regulation of human F5 promoter activity was induced via human PGRB, suggesting a conserved mechanism. Site-directed mutagenesis of the zebrafish f5 promoter further demonstrated a direct regulation of coagulation factors via progestin response elements. Moreover, a stark increase of erythrocytes occurred in capillaries meshed in wildtype preovulatory follicles but was absent in pgr-/- mutants. Interestingly, anticoagulants significantly inhibited ovulation both in vitro and in vivo, respectively. Furthermore, reduced fecundity was observed in f5+/- female zebrafish. Taken together, our study provides plausible evidence for steroid regulation of coagulation factors, and a new hypothesis for blood clotting triggered ovulation in vertebrates. 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. Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

Experimental ApproachFunctional analysis of OMA demonstrates that the slope factor implemented by Black et al. affects relationships among parameters of the OMA, leading to wrong estimates of operational efficacy and affinity of an agonist. In contrast, functional analysis of OMA with Hill coefficient as a slope factor has not this flaw.

Key ResultsFitting of the OMA with Black et al. slope factor to concentration-response curves of several models (homotropic allosteric modulation, substrate inhibition and non-competitive auto-inhibition) gives wrong estimates of operational efficacy and affinity. In contrast, fitting of the OMA modified by the Hill coefficient to concentration-response curves of these several models gives correct estimates of operational efficacy and affinity.

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. Activation of the type 3 secretion system of enteropathogenic E. coli leads to remodeling of its membrane composition and function The cell envelope of gram-negative bacteria is a complex structure, essential for bacterial survival and for resistance to many antibiotics. Channels that cross the bacterial envelope and the host cell membrane form secretion systems that are activated upon attachment to host, enabling bacteria to inject effector molecules into the host cell, required for bacterial-host interaction. The type III secretion system (T3SS) is critical for the virulence of several pathogenic bacteria, including enteropathogenic E. coli (EPEC). The EPEC T3SS activation is associated with repression of carbon storage regulator (CsrA), resulting in gene expression remodeling, which is known to affect EPEC central carbon metabolism and contributes to the adaptation to a cell-adherent lifestyle in a poorly understood manner. We reasoned that the changes in bacterial envelope upon attachment to host and the activation of a secretion system may involve a modification of the lipid composition of bacterial envelope. Accordingly, we performed a lipidomics analysis on mutant strains that simulate T3SS activation. We saw a shift in glycerophospholipid metabolism towards the formation of lysophospholipids, attributed to corresponding upregulation of the phospholipase pldA and the acyltransferase ygiH upon T3SS activation in EPEC. We also detected a shift from menaquinones and ubiquinones to undecaprenyl lipids, concomitant to abnormal synthesis of O-antigen. The remodeling of lipid metabolism is mediated by CsrA and associated with increased bacteria cell size and Zeta potential, and a corresponding alteration in EPEC permeability to vancomycin, increasing the sensitivity of T3SS-activated strains and of adherent wild type EPEC to the antibiotic.

ImportanceThe characterization of EPEC membrane lipid metabolism upon attachment to host is an important step towards a better understanding the shift of EPEC, a notable human pathogen, from a planktonic to adherent life style. It may also apply to other pathogenic bacteria that use this secretion system. We predict that upon attachment to host cells the lipid remodeling upon T3SS activation contributes to bacterial fitness and promotes host colonization, and show that it is associated with increased cell permeability and higher sensitivity to vancomycin.

To the best of our knowledge, this is the first demonstration of a bacterial lipid remodeling due to activation of a secretion system. Reconstitution of minimal motility system based on Spiroplasma swimming by expressing two bacterial actins in synthetic minimal bacterium Motility is one of the most important features of life, but its evolutionary origin remains unknown. In this study, we focused on Spiroplasma, commensal, or parasitic bacteria. They swim by switching the helicity of a ribbon-like cytoskeleton that comprises six proteins, each of which evolved from a nucleosidase and bacterial actin called MreB. We expressed these proteins in a synthetic, non-motile minimal bacterium, JCVI-syn3.0B, whose reduced genome was computer-designed and chemically synthesized. The synthetic bacterium exhibited swimming motility with features characteristic of Spiroplasma swimming. Moreover, some combinations of the two proteins produced a helical cell shape and swimming, suggesting that the swimming originated from the differentiation and coupling of bacterial actins, and we also obtained a minimal system for motility of the synthetic bacterium.

One-Sentence SummaryThe minimal system comprised two bacterial actins that provided cell helicity and swimming to the synthetic minimal bacterium. Ubiquity and origins of structural maintenance of chromosomes (SMC) proteins in eukaryotes Structural maintenance of chromosomes (SMC) protein complexes are common in Bacteria, Archaea, and Eukaryota. SMC proteins, together with the proteins related to SMC (SMC-related proteins), constitute a superfamily of ATPases. Bacteria/Archaea and Eukaryotes are distinctive from one another in terms of the repertory of SMC proteins. A single type of SMC protein is dimerized in the bacterial and archaeal complexes, whereas eukaryotes possess six distinct SMC subfamilies (SMC1-6), constituting three heterodimeric complexes, namely cohesin, condensin, and SMC5/6 complex. Thus, to bridge the homodimeric SMC complexes in Bacteria and Archaea to the heterodimeric SMC complexes in Eukaryota, we need to invoke multiple duplications of an SMC gene followed by functional divergence. However, to our knowledge, the evolution of the SMC proteins in Eukaryota had not been examined for more than a decade. In this study, we reexamined the ubiquity of SMC1-6 in phylogenetically diverse eukaryotes that cover the major eukaryotic taxonomic groups recognized to date and provide two novel insights into the SMC evolution in eukaryotes. First, multiple secondary losses of SMC5 and SMC6 occurred in the eukaryotic evolution. Second, the SMC proteins constituting cohesin and condensin (i.e., SMC1-4), and SMC5 and SMC6 were derived from closely related but distinct ancestral proteins. Based on the above-mentioned findings, we discuss how SMC1-6 have diverged from the archaeal homologs. A platform of patient-derived microtumors identifies therapeutic vulnerabilities in ovarian cancer. BackgroundIn light of the frequent development of therapeutic resistance in cancer treatment, there is a strong need for personalized model systems representing patient tumor heterogeneity, while enabling parallel drug testing and identification of appropriate treatment responses in individual patients. Using ovarian cancer as a prime example of a heterogeneous tumor disease with complex microenvironment and high recurrence rates, we developed a 3D preclinical ovarian cancer model comprised of patient-derived microtumors (PDM) and autologous tumor-infiltrating lymphocytes (TILs) for identification of treatment vulnerabilities and validation of anti-cancer drug efficacy using immunohistochemistry, immune cell phenotyping, functional assays and protein profiling analyses.

MethodsPDM and TILs were isolated from fresh primary ovarian cancer tissue specimen using mechanical disruption and limited enzymatic digestion and were subsequently cultured in suspension in defined media in the absence of serum. The heterogeneous cellular composition of isolated PDM as well as autologous TILs was analyzed by FFPE immunohistochemistry and multi-color flow cytometry, respectively. For in-depth protein profiling of PDM we established Reverse Phase Protein Array (RPPA) analyses of >110 total and phospho-proteins. Treatment efficacy in response to chemotherapeutics as well as immunotherapeutic compounds was assessed in PDM and PDM-TIL co-cultures using a functional viability assay in microplate format.

ResultsThe enzymatic digestion of primary ovarian cancer tissue and suspension culture in defined serum-free media allowed fast and efficient recovery of patient-derived microtumors (PDM). Immunohistochemical analyses demonstrated histopathological comparability of ovarian cancer PDM with corresponding patient tumor tissue. Reverse Phase Protein Array (RPPA)-based analyses of >110 total and phospho-proteins enabled the identification of patient-specific sensitivities to standard, platinum-based therapy and thereby the prediction of potential treatment-responders. Finally, combining PDM and autologous TILs for individual efficacy testing of immune checkpoint inhibitor treatment demonstrated the potential for patient-specific enhancement of cytotoxic TIL activity by this therapeutic approach.

Conclusion3D patient-derived ovarian cancer microtumors represent a preclinical, ex vivo tumor model that reflects intertumoral heterogeneity and represent the cellular complexity of individual patient tumors. Combining protein pathway analysis and anti-cancer drug efficacy testing of PDM enables drug mode-of-action analyses and therapeutic sensitivity prediction within a clinically relevant time frame after surgery. Follow-up studies in larger cohorts are currently under way to further evaluate the applicability of this platform to support clinical decision-making and personalizing cancer treatment. Identification of CD4+ Sub-population of Resident Cardiac Fibroblasts Linked to Myocardial Fibrosis Infiltration with inflammatory T-cells and accumulation of cardiac myofibroblasts are hallmarks of cardiac fibrosis and maladaptive remodeling. The origin, identity, and functions of the resident cardiac cells involved in this process are, however, unclear. To determine the identity of cells contained in regions exhibiting fibrosis, mass cytometry profiling was performed using resident human ventricular cardiac fibroblasts and right ventricle autopsy tissues from individuals diagnosed with pulmonary hypertension and SUGEN/hypoxia rats. Results showed that a subpopulation of resident myocardial fibroblasts expresses increased levels of CD4+, a helper T-cell surface marker, in addition to mesenchymal markers in humans and rats. Characterization of the resident cardiac fibroblast subpopulation, both structurally and functionally, using transcriptome and secretome analysis of the secreted cytokines, chemokines, proteins, and metabolites, evidenced that IL-1{beta} induces a phenotypic switch of human cardiac fibroblasts from mesenchymal to CD4+ lymphoidal lineage in vitro. RNA sequencing (RNA-seq) analysis of FACS-sorted CD4-expressing cardiac fibroblasts further revealed that the transcriptome of such IL-1{beta}-induced CD4+ fibroblast population exhibited classical lymphoidal and stem cell-like signatures. Lastly, reversal of cell clustering, phosphorylation of MAPK p38 and NF-{kappa}B p65, and phenotypic switching was achieved with the administration of an IL-1R antagonist. In conclusion, we have identified a subpopulation of cardiac fibroblasts which exhibits structural and functional attributes of both mesenchymal and lymphoid cells which is induced by IL-1{beta}-IL-1R-NFkB pathway for differentiation of cardiac fibroblast cells. These data suggest that cardiac fibroblast transdifferentiation during inflammation may form the basis for maladaptive remodeling during myocardial fibrosis. Chemomechanically voxelated niches for programmable histogenesis Stem-cell derived tissue models are commonly cultured under globally-delivered stimuli that trigger histogenesis via self-organizing activity. However, the culture of such tissue models is prone to stochastic behavior, limiting the reproducibility of cellular composition and resulting in non-physiological architectures. To overcome these shortcomings, we developed a method for printing cell niche microenvironments with microstructured cues that mediate local histogenic processes, including mechanosensing and differentiation of selected cell types. Microstructured cues include independently tunable mechano-chemical properties, with conjugated peptides, proteins, and morphogens across a range of Youngs moduli. By rationally designing niches, we mediate the structure of tissues derived from stem-cell-progenitor sources, including a bone-fat assembly from stromal mesenchyme, and embryonic tissues derived from hiPSC. We show that microstructured cues can recapitulate mechano-chemical signals resembling early embryonic histogenesis. This outcome includes a role for niche mechanics in human embryonic organization, where soft niche mechanics bias markers of mesendodermal differentiation and epithelial-to-mesenchymal-transition (EMT), as well as a demonstration of a material-mediated morphogen signaling centers able to induce foci of mesenchymal and EMT differentiation. Thus, microstructured materials can mediate local histogenic processes to enhance the structure and composition of tissue models. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. Sgs1BLM independent role of Dna2DNA2 nuclease at DNA double strand break is inhibited by Nej1XLF The two major pathways of DNA double strand break (DSB) repair, non-homologous end-joining (NHEJ) and homologous recombination (HR), are highly conserved from yeast to mammals. The regulation of 5 DNA resection controls repair pathway choice and influences repair outcomes. Nej1 was first identified as a canonical NHEJ factor involved in stimulating the ligation of broken DNA ends and more recently, it was shown to be important for DNA end-bridging and inhibiting 5 resection mediated by Dna2-Sgs1. Nej1 interacts with Sae2 and this impacts DSB repair in three ways. First, Nej1 inhibits MRX-Sae2 interactions and Sae2 localization to a DSB. Second, Nej1 inhibits Sae2-dependent recruitment of Dna2 in the absence of Sgs1. Third, NEJ1 and SAE2 showed an epistatic relationship for DNA end-bridging, an event that restrains the broken ends and reduces the frequency of genomic deletions from developing at the DSB. Deletion of NEJ1 suppressed the synthetic lethality of sae2{Delta} sgs1{Delta} and was dependent on the nuclease activity of Dna2. These Nej1 functions promote end-joining DSB repair, but could also be relevant for controlling resection initiation during HR repair.

HighlightsO_LINej1 physically interacts with Sae2 and inhibits end-resection at a DSB.
C_LIO_LINej1 inhibits Sae2 interactions with the MRX complex.
C_LIO_LINej1 inhibits Sae2-dependent recruitment of Dna2 to a DSB.
C_LIO_LINEJ1 and SAE2 are epistatic for DNA end-bridging.
C_LIO_LIDeletion of NEJ1 suppresses the synthetic lethality of sae2{Delta} sgs1{Delta}, which is dependent on Dna2 nuclease activity.
C_LI Evidence for spatiotemporal shift in demersal fisheries management priority areas in the western Mediterranean Marine Protected Areas (MPAs) are a promising management tool for the conservation and recovery of marine ecosystems, as well as fisheries management. MPAs are generally established as permanent closures but marine systems are dynamic, which has generated debate in favour of more dynamic designs. As a consequence, the identification of priority areas should assess their persistence in space and time. Here, we develop a step-by-step approach to assess the spatiotemporal dynamics of fisheries management priority areas using standard fishery-independent survey data. To do so, we fit Bayesian hierarchical spatiotemporal SDM models to different commercially important demersal species and use the resulting maps to fit different spatial prioritisation configurations. We use these results to assess the spatiotemporal dynamics of fisheries priority areas. The proposed method is illustrated through a western Mediterranean case study using fishery-independent trawl survey data on six commercially important species collected over 17 years. We identified two fisheries priority area patterns in the study area, each predominant during a different time-period of the study, asserting the importance of regularly re-assessing MPA designs. 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. Multi-regional module-based signal transmission in mouse visual cortex The visual cortex is organized hierarchically, but the presence of extensive recurrent and parallel pathways make it challenging to decipher how signals flow between neuronal populations. Here, we tracked the flow of spiking activity recorded from six interconnected levels of the mouse visual hierarchy. By analyzing leading and lagging spike-timing relationships among pairs of simultaneously recorded neurons, we created a cellular-scale directed network graph. Using a module-detection algorithm to cluster neurons based on shared connectivity patterns, we uncovered two multi-regional communication modules distributed across the hierarchy. The direction of signal flow between and within these modules, differences in layer and area distributions, and distinct temporal dynamics suggest that one module is positioned to transmit feedforward sensory signals, whereas the other integrates inputs for recurrent processing. These results suggest that multi-regional functional modules may be a fundamental feature of organization beyond cortical areas that supports signal propagation across hierarchical recurrent networks. 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. The OPAQUE1/DISCORDIA2 myosin XI is required for phragmoplast guidance during asymmetric cell division in maize Asymmetric divisions produce cells with different fates and are critical for development. Here we show that a maize myosin XI protein, OPAQUE1 (O1), is necessary for asymmetric divisions during maize stomatal development. We analyzed stomatal precursor cells prior to and during asymmetric division to determine why o1 mutants have abnormal division planes. Cell polarization and nuclear positioning occur normally in the o1 mutant, and the future site of division is correctly specified. The defect in o1 occurs during late cytokinesis, when the plant-specific phragmoplast -- made of microtubules, actin and other proteins - forms the nascent cell plate. The phragmoplast becomes misguided and does not meet the previously established division site. Initial phragmoplast guidance is correct in o1. However, as phragmoplast expansion continues, phragmoplasts in o1 stomatal precursor cells become misguided and do not meet the cortex at the established division site. To understand how this myosin protein contributes to phragmoplast guidance, we identified O1-interacting proteins. Other myosins, specific actin-binding proteins, and maize kinesins related to the Arabidopsis thaliana division site markers PHRAGMOPLAST ORIENTING KINESINs (POKs) interact with O1. We propose that different myosins are important at multiple steps of phragmoplast expansion, and the O1 actin motor and POK-like microtubule motors work together to ensure correct late-stage phragmoplast guidance. Orientation dependent proton transverse relaxation in the human brain white matter: The magic angle effect on a cylindrical helix PurposeTo overcome some limitations of prior orientation-dependent proton transverse relaxation formalisms in white matter (WM) with a novel framework based on the generalized magic angle effect function.

MethodsA cylindrical helix model was developed embracing anisotropic rotational and translational diffusion of restricted molecules in human brain WM, with the former characterized by an axially symmetric system. Transverse relaxation rates R2 and [Formula] were divided into isotropic [Formula] and anisotropic parts, [Formula], with  denoting an open angle and{varepsilon} 0 an orientation ({Phi}) offset from DTI-derived primary diffusivity direction. The proposed framework (Fit A) was compared with prior models without{varepsilon} 0 on previously published water and methylene proton transverse relaxation rates from developing, healthy, and pathological WM at 3T. Goodness of fit was represented by root-mean-square error (RMSE). F-test and linear correlation were used with statistical significance set to P [&le;] 0.05.

ResultsFit A significantly (P<0.01) outperformed prior models as demonstrated by reduced RMSEs, e.g., 0.349 vs. 0.724 in myelin water. Fitted{varepsilon} 0 was in good agreement with calculated{varepsilon} 0 from directional diffusivities. Compared with those from healthy adult, the fitted [Formula], and  from neonates were substantially reduced but{varepsilon} 0 increased, consistent with incomplete myelination. Significant positive [Formula] and negative ( and [Formula]) correlations were found with aging (demyelination) in elderly.

ConclusionThe developed framework can better characterize orientation dependences from a wide range of proton transverse relaxation measurements in human brain WM, shedding new light on myelin microstructural alterations at the molecular level. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Experimental evolution reveals post-transcriptional regulation as a novel driver of Leishmania fitness gain The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network driving parasite fitness, with genome instability causing highly reproducible, gene dosage-dependent changes in protein abundance linked to post-transcriptional regulation. These in turn were associated with a gene dosage-independent reduction in abundance of flagellar transcripts and a coordinated increase in abundance of coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of fitness-adapted ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings. Comparative Analysis of non-coding small-RNAs in Pseudomonas aeruginosa Keratitis Strains with Different Antibiotic Susceptibility Pseudomonas aeruginosa, is a gram-negative bacterium causes opportunistic or nosocomial infections in immunocompromised individuals. In recent years, a steady increase in human corneal infections of P. aeruginosa has been reported with increased multi-drug resistance (MDR) or extensively drug resistance (XDR). Several non-coding sRNAs, has been identified to regulate various physiological processes in P. aeruginosa, including biofilm formation, quorum sensing. However, the regulatory mechanism of sRNAs in MDR/XDR pathways of P. aeruginosa keratitis strains is not yet studied. In this study, we identified bacterial sRNAs in publicly available P. aeruginosa keratitis genomes and investigated their regulatory role in MDR/XDR pathways using bioinformatic analysis. Totally, 46 P. aeruginosa keratitis strains from different geographical regions were included. Of 46, Eight (30%) out of Twenty-seven and Nine (52) out of Nineteen P. aeruginosa strains from India and Australia were identified as not-MDR. Whereas, 10 (38%) Indian and 9 (47%) Australian strains were identified as MDR. Eight Indian strains were identified as XDR. Out of 46 strains, 23 (50%) carried ExoU, 21(45%) carried ExoS and two (5%) strains carried both ExoU and ExoS, exotoxins for their virulence. The sRNA, SPA0021 was identified in 18 MDR/XDR and 6 not-MDR strains along with UCBPP-PA14. Interestingly, majority of the imipenem resistant P. aeruginosa keratitis strains from the present study was found to be carried SPA0023 sRNA (18 out of 30 strains). The outer membrane porin protein OprD, identified as binding target of SPA0023. Negative regulation or inactivation of OprD, reported in increased imipenem resistance in P. aeruginosa. Mutation analysis revealed that SPA0023 carrying P. aeruginosa keratitis strains contains a lesser number of amino acid changes in OprD protein than other strains. These findings indicate, imipenem resistance in SPA0023 carried strains might arose from the negative regulation or inhibition of OprD by SPA0023. However, functional studies are warranted with large number of P. aeruginosa keratitis strains to confirm the negative regulation of OprD by SPA0023 and imipenem resistance. Chromatin-contact atlas reveals disorder-mediated protein interactions and moonlighting chromatin-associated RBPs RNA-binding proteins (RBPs) play diverse roles in regulating co-transcriptional RNA-processing and chromatin functions, but our knowledge of the repertoire of chromatin-associated RBPs (caRBPs) and their interactions with chromatin remains limited. Here, we developed SPACE (Silica Particle Assisted Chromatin Enrichment) to isolate global and regional chromatin components with high specificity and sensitivity, and SPACEmap to identify the chromatin-contact regions in proteins. Applied to mouse embryonic stem cells, SPACE identified 1,459 chromatin-associated proteins, [~]48% of which are annotated as RBPs, indicating their dual roles in chromatin and RNA-binding. Additionally, SPACEmap stringently verified chromatin-binding of 404 RBPs and identified their chromatin-contact regions. Notably, SPACEmap showed that about half of the caRBPs bind to chromatin by intrinsically disordered regions (IDRs). Studying SPACE and total proteome dynamics from mES cells grown in 2iL and serum medium indicates significant correlation (R = 0.62). One of the most dynamic caRBPs is Dazl, which we find co-localized with PRC2 at transcription start sites of genes that are distinct from Dazl mRNA binding. Dazl and other PRC2-colocalised caRBPs are rich in intrinsically disordered regions (IDRs), which could contribute to the formation and regulation of phase-separated PRC condensates. Together, our approach provides an unprecedented insight into IDR-mediated interactions and caRBPs with moonlighting functions in native chromatin. Risky choice: probability weighting explains Independence Axiom violations in monkeys Expected Utility Theory (EUT) provides axioms for maximizing utility in risky choice. The Independence Axiom (IA) is its most demanding axiom: preferences between two options should not change when altering both options equally by mixing them with a common gamble. We tested common consequence (CC) and common ratio (CR) violations of the IA over several months in thousands of stochastic choices using a large variety of binary option sets. Three monkeys showed consistently few outright Preference Reversals (8%) but substantial graded Preference Changes (46%) between the initial preferred gamble and the corresponding altered gamble. Linear Discriminant Analysis (LDA) indicated that gamble probabilities predicted most Preference Changes in CC (72%) and CR (88%) tests. The Akaike Information Criterion indicated that probability weighting within Cumulative Prospect Theory (CPT) explained choices better than models using Expected Value (EV) or EUT. Fitting by utility and probability weighting functions of CPT resulted in nonlinear and non-parallel indifference curves (IC) in the Marschak-Machina triangle and suggested IA non-compliance of models using EV or EUT. Indeed, CPT models predicted Preference Changes better than EV and EUT models. Indifference points in out-of-sample tests were closer to CPT-estimated ICs than EV and EUT ICs. Finally, while the few outright Preference Reversals may reflect the long experience of our monkeys, their more graded Preference Changes corresponded to those reported for humans. In benefitting from the wide testing possibilities in monkeys, our stringent axiomatic tests contribute critical information about risky decision-making and serve as basis for investigating neuronal decision mechanisms.