Extracellular Vesicle (EV) Mitochondria and EV/HSP27 Mixtures Protect the Ischemic Brain Endothelium Ischemic stroke causes brain endothelial cell (BEC) death and damages tight junction integrity of the blood-brain barrier (BBB). We harnessed the innate mitochondrial load of endothelial cell-derived extracellular vesicles (EVs) and utilized mixtures of EV/exogenous heat shock protein 27 (HSP27) as a one-two punch strategy to increase BEC survival (via EV mitochondria) and preserve their tight junction integrity (via HSP27 effects). We demonstrated that the medium-to-large (m/lEV) but not small EVs (sEV) transferred their mitochondrial load, which subsequently colocalized with the mitochondrial network of the recipient primary human BECs. BECs treated with m/lEVs increased relative ATP levels and displayed superior mitochondrial function. m/lEVs isolated from rotenone (mitochondrial complex I inhibitor) exposed BECs (RTN-m/lEVs) did not increase BECs ATP levels compared to naive m/lEVs. In contrast, RTN-sEVs functionality was minimally affected compared to naive sEVs. Compared to EV/HSP27 mixtures, PEG-DET/EV/HSP27 ternary mixtures increased HSP27 protein loading. (PEG-DET/HSP27)/EV ternary and EV/HSP27 binary mixtures decreased the paracellular permeability of small and large molecular mass fluorescent tracers in oxygen glucose-deprived primary human BECs. This one-two punch approach to increase BEC metabolic function and tight junction integrity is a promising strategy for BBB protection and prevention of long-term neurological dysfunction post-ischemic stroke.

Highlights[tpltrtarr] Medium-to-large extracellular vesicles (m/lEVs), not small EVs contain mitochondria
[tpltrtarr]m/lEVs increased ATP and mitochondrial function in brain endothelial cells (BECs)
[tpltrtarr]Donor m/lEV mitochondria colocalized with host mitochondria in recipient BECs
[tpltrtarr]m/lEVs from rotenone-exposed BECs did not increase recipient BEC ATP levels
[tpltrtarr]PEG-DET/EV/HSP27 ternary mixtures increase exogenous HSP27 protein loading


Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=133 SRC="FIGDIR/small/466491v4_ufig1.gif" ALT="Figure 1">
View larger version (43K):
org.highwire.dtl.DTLVardef@155f618org.highwire.dtl.DTLVardef@8fa778org.highwire.dtl.DTLVardef@13bd44aorg.highwire.dtl.DTLVardef@ee771e_HPS_FORMAT_FIGEXP  M_FIG C_FIG
Impedance scaling for gold and platinum microelectrodes ObjectiveElectrical measurement of the activity of individual neurons is a primary goal for many invasive neural electrodes. Making these "single unit" measurements requires that we fabricate electrodes small enough so that only a few neurons contribute to the signal, but not so small that the impedance of the electrode creates overwhelming noise or signal attenuation. Thus, neuroelectrode design often must strike a balance between electrode size and electrode impedance, where the impedance is often assumed to scale linearly with electrode area.

Approach and main resultsHere we study how impedance scales with neural electrode area and find that the 1 kHz impedance of Pt electrodes (but not Au electrodes) transitions from scaling with area (r-2) to scaling with perimeter (r-1) when the electrode radius falls below 10 microns. This effect can be explained by the transition from planar to spherical diffusion behavior previously reported for electrochemical microelectrodes. Significance: These results provide important intuition for designing small, single unit recording electrodes. Specifically, for materials where the impedance is dominated by a pseudo-capacitance that is associated with a diffusion limited process, the total impedance will scale with perimeter rather than area when the electrode size becomes comparable with the diffusion layer thickness. For Pt electrodes this transition occurs around 10 um radius electrodes. At even lower frequencies (1 Hz) impedance approaches a constant. This transition to r-1 scaling implies that electrodes with a pseudo-capacitance can be made smaller than one might expect before thermal noise or voltage division limits the ability to acquire high-quality single-unit recordings.
A map of binary SARS-CoV-2 protein interactions implicates host immune regulation and ubiquitination Key steps in viral propagation, immune suppression, and pathology are mediated by direct, binary, physical interactions between viral and host proteins. To understand the biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, we generated an unbiased systematic map of binary interactions between viral and host proteins, complementing previous co-complex association maps by conveying more direct mechanistic understanding and potentially enabling targeted disruption of direct interactions. To this end, we deployed two parallel strategies, identifying 205 virus-host and 27 intraviral binary interactions amongst 171 host and 19 viral proteins, and confirming high quality of these interactions via a calibrated orthogonal assay. Host proteins interacting with SARS-CoV-2 proteins are enriched in various cellular processes, including immune signaling and inflammation, protein ubiquitination, and membrane trafficking. Specific subnetworks provide new hypotheses related to viral modulation of host protein homeostasis and T-cell regulation. The binary virus-host protein interactions we identified can now be prioritized as targets for therapeutic intervention. More generally, we provide a resource of systematic maps describing which SARS-CoV-2 and human proteins interact directly. Building bridges: the hippocampus constructs narrative memories across distant events Lifes events are scattered throughout time, yet we often recall different events in the context of an integrated narrative. Prior research suggests that the hippocampus, which supports memory for past events, can sometimes support integration of overlapping associations or separate events, but the conditions which lead to hippocampus-dependent memory integration are unclear. We used functional brain imaging to test whether the ability to form a larger narrative (narrative coherence) drives hippocampal memory integration. During encoding of fictional stories, the hippocampus supported patterns of activity, including activity at boundaries between events, which were more similar between distant events that formed one coherent narrative, compared with events taken from unrelated narratives. Twenty-four hours later, the hippocampus preferentially supported detailed recall of coherent narrative events, through reinstatement of hippocampal activity patterns from encoding. These findings reveal a key function of the hippocampus: the dynamic integration of events into a narrative structure for memory. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Development of a Novel PROTAC using the Nucleic Acid Aptamer as a Targeting Ligand for Tumor Selective Degradation of Nucleolin PROteolysis TArgeting Chimeras (PROTACs) induce targeted protein degradation by hijacking the intracellular ubiquitin proteasome system, thus emerging as a new strategy for drug development. However, most PROTACs generated lack cell-type selectivity and have poorly soluble in water. To address this drawback, we developed a novel PROTAC ZL216 using aptamer AS1411 as a targeting ligand of nucleolin to conjugate with a small molecule ligand of E3 ligase VHL, which shows high aqueous solubility and serum stability. Based on the differential expression of nucleolin on the cell surface, ZL216 could bind to and internalize into breast cancer cells, but not normal breast cells. Furthermore, we revealed that ZL216 promoted the formation of a nucleolin-ZL216-VHL ternary complex in breast cancer cells and potently induced nucleolin degradation in vitro and in vivo. As a result, ZL216 inhibited the proliferation and migration of breast cancer cells. These studies demonstrate that in addition to peptides and small molecule compounds, nuclei acid aptamers can also be used to generate PROTACs, which broadens the toolbox constructing PROTACs and provides a promising strategy for development of tumor-selective PROTACs. 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. Individual-based multiple-unit dissimilarity: novel indices and null model for assessing temporal variability in community composition Beta-diversity was originally defined spatially, i.e., as variation in community composition among sites in a region. However, the concept of beta-diversity has since been expanded to temporal contexts. This is referred to as "temporal beta-diversity", and most approaches are simply an extension of spatial beta-diversity. The persistence and turnover of individuals over time is a unique feature of temporal beta-diversity. Nakadai (2020) introduced the "individual-based beta-diversity" concept, and provided novel indices to evaluate individual turnover and compositional shift by comparing individual turnover between two periods at a given site. However, the proposed individual-based indices are applicable only to pairwise dissimilarity, not to multiple-temporal (or more generally, multiple-unit) dissimilarity. Here, individual-based beta-diversity indices are extended to multiple-unit cases. In addition, a novel type of random permutation criterion related to these multiple-unit indices for detecting patterns of individual persistence is introduced in the present study. To demonstrate the usage the properties of these indices compared to average pairwise measures, I applied them to a dataset for a permanent 50-ha forest dynamics plot on Barro Colorado Island in Panama. Information regarding "individuals" is generally missing from community ecology and biodiversity studies of temporal dynamics. In this context, the methods proposed here are expected to be useful for addressing a wide range of research questions regarding temporal changes in biodiversity, especially studies using traditional individual-tracked forest monitoring data. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. Development and Characterization of PAR-Trackers: New Tools for Detecting Poly(ADP-ribose) In Vitro and In Vivo AO_SCPLOWBSTRACTC_SCPLOWADP-ribosylation (ADPRylation) is a reversible post-translation modification resulting in the covalent attachment of ADP-ribose (ADPR) moieties on substrate proteins. Naturally-occurring protein motifs and domains, including WWEs, PBZs, and macrodomains, act as "readers" for protein-linked ADPR. Although recombinant, antibody-like ADPR detection reagents containing these readers have facilitated the detection of ADPR, they are limited in their ability to capture the dynamic nature of ADPRylation. Herein, we describe and characterize a set of poly(ADP-ribose) (PAR) Trackers (PAR-Ts) - optimized dimerization-dependent or split-protein reassembly PAR sensors in which a naturally occurring PAR binding domain, WWE, was fused to both halves of dimerization-dependent GFP (ddGFP) or split Nano Luciferase (NanoLuc), respectively. We demonstrate that these new tools allow the detection and quantification of PAR levels in extracts, living cells, and living tissues with greater sensitivity, as well as temporal and spatial precision. Importantly, these sensors detect changes in cellular ADPR levels in response to physiological cues (e.g., hormone-dependent induction of adipogenesis without DNA damage), as well as xenograft tumor tissues in living mice. Our results indicate that PAR Trackers have broad utility for detecting ADPR in many different experimental and biological systems. Absence of cardiolipin from the outer-leaflet of a mitochondrial inner-membrane mimic restricts Opa1-mediated fusion Cardiolipin is a tetra-acylated di-phosphatidylglycerol lipid enriched in the matrix-facing (inner) leaflet of the mitochondrial inner membrane. Cardiolipin plays an important role in regulating mitochondria function and dynamics. Yet, the mechanisms connecting cardiolipin distribution and mitochondrial protein function remain indirect. In our previous work, we established an in vitro system reconstituting mitochondrial inner membrane fusion mediated by Opa1. We found that the long form of Opa1 (l-Opa1) works together with the proteolytically processed short form (s-Opa1) to mediate fast and efficient membrane fusion. Here, we extend our reconstitution system to generate supported lipid bilayers with asymmetric cardiolipin distribution. Using this system, we find the presence of cardiolipin on the inter-membrane space-facing (outer) leaflet is important for membrane tethering and fusion. We discuss how the presence of cardiolipin in this leaflet may influence protein and membrane properties, and future applications for this approach. Phase separation of +TIP-networks regulates microtubule dynamics Regulation of microtubule dynamics is essential for diverse cellular functions, and proteins that bind to dynamic microtubule ends can regulate network dynamics. Here we show, that two conserved microtubule end-binding proteins, CLIP-170 and EB3, undergo phase separation and form dense liquid-networks. When CLIP-170 and EB3 act together the multivalency of the network increases, which synergistically increases the amount of protein in the dense phase. In vitro and in cells these liquid networks can condense tubulin. In vitro in the presence of microtubules, EB3/CLIP-170 phase separation can co-condense tubulin all along the microtubule. At this condition microtubule growth speed increases up to two-fold and depolymerization events are strongly reduced, compared to conditions with phase separation deficient networks. Our data show that phase separated EB3/CLIP-170 networks impact microtubule growth dynamics beyond direct protein-microtubule interactions. Agrobacterium mediated transformation and deciphering SNPs in TaLr67 gene homeologs for gene editing in wheat. Diseases adversely affect grain yield of crop plants. Leaf rust is a major disease of wheat. As race-specific resistance breaks down, introduction of newer sources of resistance often from older accessions is necessary. Linkage drag from the donor accessions adversely impacts grain yield. As a result, CIMMYT breeding effort has shifted to using durable resistance, also referred to as race non-specific host resistance, which allows slow rusting but maintains grain yield. One of the key genes for durable resistance is Lr67 and the resistant form of this gene is absent in CIMMYT elite lines. Hence, we have initiated efforts to convert the susceptible copy of Lr67 into its resistant form directly in these elite lines using gene editing. This would eliminate backcrossing and thus save time as well as eliminate linkage drag that would accompany the resistant copy of the gene if it were to be introgressed from an older accession. As first steps, we have isolated and sequenced the genomic copies from each of the A, B, and D genome of the Lr67 gene from three elite lines and an experimental line. Identification of more than 50 single nucleotide polymorphisms (SNPs) in the open reading frames (ORF) among these lines would be useful in designing the guide RNA molecules with precision. Further, we have streamlined genetic transformation of these elite lines, a prerequisite step for gene editing. Rattus Economicus: A neural model for rational regulation of water-seeking effort in rodents In the laboratory, animals motivation to work tends to be positively correlated with reward magnitude. But in nature, rewards earned by work are essential to survival (e.g., working to find water), and the payoff of that work can vary on long timescales (e.g., seasonally). Under these constraints, the strategy of working less when rewards are small could be fatal. We found that instead, rats in a closed economy did more work for water rewards when the rewards were stably smaller, a phenomenon also observed in human labor supply curves. Like human consumers, rats also showed elasticity of demand, consuming far more water per day when its price in effort was lower. The neural mechanisms underlying such "rational" market behaviors remain largely unexplored. We propose a dynamic utility maximization model that can account for the dependence of rat labor supply (trials/day) on the wage rate (ml/trial), and also predict the temporal dynamics of when rats work. Based on data from mice, we hypothesize that SFOGLUT neurons in lamina terminalis continuously compute the instantaneous marginal utility of voluntary work for water reward, and causally determine the amount and timing of work. Venomics and antivenomics of Indian spectacled cobra (Naja naja) from the Western Ghats Venom proteome profiling of Naja naja from the Western Ghats region in Kerala was achieved through SDS-PAGE and RP-HPLC followed by Q-TOF LC-MS/MS analysis, incorporating PEAKS and Novor assisted de novo sequencing methodologies. A total of 115 proteins distributed across 17 different enzymatic and non-enzymatic venom protein families were identified through conventional and 39 peptides through homology-driven proteomics approaches. Fourteen peptides derived through de novo complements the Mascot data indicating the importance of homology-driven approaches in improving protein sequence information. Among the protein families identified, glutathione peroxidase and endonuclease were reported for the first time in the Indian cobra venom. Immunological cross-reactivity assessed using Indian polyvalent antivenoms suggested that VINS showed better EC50 (2.48 g/mL) value than that of PSAV (6.04 g/mL) and Virchow (6.03 g/mL) antivenoms. Western blotting experiments indicated that all the antivenoms elicited poor binding specificities, especially towards low molecular mass proteins. Second-generation antivenomics studies revealed that VINS antivenom was less efficient to detect many low molecular mass proteins such as three-finger toxins and Kunitz-type serine protease Inhibitors. Taken together, the present study enabled a large-scale characterization of the venom proteome of Naja naja from the Western Ghats and emphasized the need for developing more efficient antivenoms.

HighlightsO_LIProteomics of cobra venom resulted in the identification of 115 proteins representing 17 snake venom protein families.
C_LIO_LIDe novo approaches exclusively yielded 39 peptides harbouring multiple amino acid mutations.
C_LIO_LIGlutathione peroxidase and endonuclease were identified for the first time in Indian cobra venom.
C_LIO_LIIndian polyvalent antivenoms showed varying cross-reactivity towards cobra venom.
C_LIO_LIVINS antivenom was less efficient to detect many low molecular mass proteins (< 20 kDa).
C_LI A protein microarray-based in vitro transglutaminase assay platform for epitope mapping and peptide discovery: implication for transglutaminase-mediated immunodominant determination Transglutaminases (TGs) are a family of crosslinking enzymes catalyzing the formation of intra- and inter-molecular glutamine-lysine isopeptide bonds in a calcium dependent manner. Activation of transglutaminases is pathogenically associated with severe human diseases including neurodegenerations, cardiovascular diseases, and autoimmune diseases. Although continuous efforts determining the enzymes substrate preference have been witnessed, a high-throughput assay platform with the "omic" efficiency is still missing for the global identification of substrate-specific TG modification sites. In this study we report a protein microarray-based in vitro TG assay platform for rapid and large-scale (up to 30000 reactions per chip) determination of the glutamine (Q)-bearing TG modification motifs. With this platform we identified the Q16 in superoxide dismutase 1 and Q109 in alpha-synuclein as the primary modification sites for tissue transglutaminase (TG2), the most ubiquitous member of the enzyme family. Of particular interest, we found a close match between the modification motifs and published vaccine epitope sequences in alpha-synuclein, implying an essential and intrinsic role transglutaminase might play in the determination of immunodominant epitopes. Our data collectively suggest the glutamine and its follow-up five residues on the C terminal of a protein compose a minimal determinant motif for TG2 modification and the TG2 modification motifs determined by our platform could finally correspond to the substrates epitope sequences in antigen processing. To screen for site-specific interfering peptides and assist gene editing-based protein engineering, we optimized through onchip amino-acid scanning the TG2 modification motif QQIV in the extracellular matrix protein fibronectin and obtained 14 variants with significantly higher TG2 reactivity that might serve as the competitive inhibitor peptides and 1 with lower reactivity. We further confirmed this optimization approach with 12-mer peptides, the longest ones that could be synthesized on the chip. Taken together, our synthetic transglutaminase assay platform might be able to deliver a precise epitope blueprint for immunotherapeutic targeting and provide proof-of-concept and directional studies for TG-based peptide discovery and protein design. Recovery of high-quality assembled genomes via metagenome binning guided with single-cell amplified genomes BackgroundObtaining high-quality (HQ) reference genomes from microbial communities is crucial for understanding the phylogeny and function of uncultured microbes in complex microbial ecosystems. Despite the improved bioinformatic approaches to generate curated metagenome-assembled genomes (MAGs), existing metagenomic binners often fail to obtain reliable MAGs, and thus, they are nowhere comparable to genomes sequenced from isolates in terms of strain level resolution. Here, we present a single-cell genome-guided metagenome binning (MetaSAG) to reconstruct the strain-resolved genomes from microbial communities at once.

ResultsMetaSAG employs single-cell amplified genomes (SAGs) generated with microfluidic technology as binning guides to recover improved draft genomes with the metagenomic data. To assess the performance of reconstructing genomes from various microbial communities, we compared MetaSAG with four conventional metagenomic binners using a cell mock community, human gut microbiota, and skin microbiota samples. MetaSAG showed precise contig binning and higher recovery rates (>97%) of rRNA and plasmids compared to conventional binners in genome reconstruction from the cell mock community. In human microbiota samples, MetaSAG recovered the largest number of genomes with a total of 103 gut microbial genomes (21 HQ and 65 showed >90% completeness) and 45 skin microbial genomes (10 HQ and 40 showed >90% completeness), respectively. Conventional binners recovered one Staphylococcus hominis genome, whereas MetaSAG recovered two S. hominis genomes from the identical skin microbiota sample. Single-cell sequencing indicated that these S. hominis genomes clearly derived from two distinct strains harboring specifically different plasmids. We found that all conventional S. hominis MAGs had substantial lack or excess of the genome sequences and contamination of other Staphylococcus bacteria (S. epidermidis).

ConclusionsMetaSAG enabled us to obtain the strain-resolved genomes in the mock community and human microbiota samples by assigning metagenomic sequences correctly and covering both highly conserved genes such as rRNA genes and unique extrachromosomal elements, including plasmids. MetaSAG will provide HQ genomes that are difficult to obtain with metagenomic analyses alone and will facilitate the understanding of microbial ecosystems by elucidating detailed metabolic pathways and horizontal gene transfer networks through microbial genomes. MetaSAG is available at https://github.com/kojiari/metasag. Phenotypic heterogeneity driven by plasticity of the intermediate EMT state governs disease progression and metastasis in breast cancer The Epithelial-to-Mesenchymal Transition (EMT) is a developmental cellular program frequently coopted by cancer cells1 and is a key contributor to both heterogeneity in solid tumors2-4 and later stage chemo-resistance and metastasis5,6. Rather than being a switch from an epithelial to a mesenchymal state, increasing evidence points to the existence of intermediate EMT states, wherein cells co-express both epithelial and mesenchymal traits7-13. Multiple stable intermediate EMT states possessing unique characteristics exist across the EMT spectrum7,8,14,15, contributing to the complex heterogeneity of tumors to promote metastasis16. While much work has been carried out identifying and characterizing EMT-inducing transcription factors17-20, the transcriptional and epigenetic networks responsible for the stability and maintenance of the midpoints along the EMT spectrum are poorly defined. In addition, there are currently no approaches to identifying and quantifying intermediate EMT subpopulations within patient tumors to evaluate their prognostic significance. Using clonally isolated derivatives of the SUM149PT breast cancer cell line, we systematically interrogate how each EMT state independently contributes to heterogeneity and influences metastatic progression, uncovering the role of RUNX2 in stabilizing certain intermediate states. Using SUM149PT-derived tumors as a training set, we develop an entropy-based model to quantify phenotypic heterogeneity and EMT status. Remarkably, the majority of cell states captured in the SUM149PT model are represented in a selection of patient tumors, laying the foundation for quantification of epithelial-mesenchymal heterogeneity and understanding the role of the intermediate EMT state in tumor progression. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Classification of Pacemaker Dynamics in the Mouse Intestine by Field Potential Microimaging The flexible and sophisticated movement of the gastrointestinal (GI) tract implies the involvement of mechanisms other than enteric neural circuits, to coordinate excitation in microregions. We thus performed microimaging of pacemaker dynamics in the small intestine of mice since it contains typical network-forming pacemaker cells. A dialysis membrane-reinforced low-impedance microelectrode array (MEA) enabled field potentials over a wide frequency range to be stably measured in microregions. The pacemaker dynamics were classified into basic patterns despite large variations. In the developmental process, pacemaker activity was categorized as either an  expanding or a  migrating pattern that was initiated in or propagated to the MEA sensing area, respectively. The intercellular current of the volume conductor complicated the waveform of both activities. The existence of  expanding and  migrating patterns was attributable to duplicated pacemaker systems such as intracellular Ca2+ oscillation-activated and voltage-gated mechanisms. Additionally, from the spatio-temporal feature during the period of pacemaker events, the  bumpy/aberrant pattern was defined by aberrant, incoherent propagation, and associated with local impairment of excitability, while the  colliding/converging pattern involved the interaction of multiple activities in the MEA area. Interconversion between the four micro-coordination patterns occurred in the same microregion. 5-Hydroxytryptamine (5-HT) promoted  migrating activity, implying an improvement or restoration of spatial conductivity. These results agree well with the action of 5-HT to change GI movement toward propulsion. In conclusion, our MEA method of microimaging classification enables the quantitative assessment of spatio-temporal electric coordination underlying GI motility, suggesting its application to small model animals. Somatic uniparental disomy mitigates the most damaging EFL1 allele combination in Shwachman-Diamond syndrome Shwachman-Diamond syndrome (SDS; OMIM: #260400) is caused by variants in SBDS (Shwachman-Bodian-Diamond syndrome gene), which encodes a protein that plays an important role in ribosome assembly. Recent reports suggest that recessive variants in EFL1 are also responsible for SDS. However, the precise genetic mechanism that leads to EFL1-induced SDS remains incompletely understood. Here we present three unrelated Korean SDS patients that carry biallelic pathogenic variants in EFL1 with biased allele frequencies, resulting from a bone marrow-specific somatic uniparental disomy (UPD) in chromosome 15. The recombination events generated cells that were homozygous for the relatively milder variant, allowing for the evasion of catastrophic physiological consequences. Still, the milder EFL1 variant was solely able to impair 80S ribosome assembly and induce SDS features in cell line, zebrafish, and mouse models. The loss of EFL1 resulted in a pronounced inhibition of terminal oligo-pyrimidine element-containing ribosomal protein transcript 80S assembly. Therefore, we propose a more accurate pathogenesis mechanism of EFL1 dysfunction that eventually leads to aberrant translational control and ribosomopathy. Full-likelihood analysis of genomic data clarifies a complex history of species divergence and introgression in the erato-sara group of Heliconius butterflies AO_SCPLOWBSTRACTC_SCPLOWIntrogressive hybridization plays a key role in adaptive evolution and species diversification in many groups of species. However, frequent hybridization and gene flow between species makes estimation of the species phylogeny and key population parameters challenging. Here, we show that by accounting for phasing and using full-likelihood analysis methods, introgression histories and population parameters can be estimated reliably from whole-genome sequence data. We employ full-likelihood methods under the multispecies coalescent (MSC) model with and without gene flow to analyze the genomic data from six members of the erato-sara clade of Heliconius butterflies and infer the species phylogeny and cross-species introgression events. The methods naturally accommodate random fluctuations in genealogical history across the genome due to deep coalescence. To avoid heterozygote phasing errors in haploid sequences commonly produced by genome assembly methods, we process and compile unphased diploid sequence alignments and use analytical methods to average over uncertainties in heterozygote phase resolution. There is robust evidence for introgression across the genome, both among distantly related species deep in the phylogeny and between sister species in shallow parts of the tree. We obtain chromosome-specific estimates of key population parameters such as introgression directions, times and probabilities, as well as species divergence times and population sizes for modern and ancestral species. We confirm ancestral gene flow between the sara clade and an ancestral population of H. telesiphe, a likely hybrid speciation origin for H. hecalesia, and gene flow between sister species H. erato and H. himera. Inferred introgression among ancestral species also explains the history of two chromosomal inversions deep in the phylogeny of the group. This study illustrates how a full-likelihood approach based on the multispecies coalescent makes it possible to extract rich historical information of species divergence and gene flow from genomic data. Integrated annotation and analysis of genomic features reveal new types of functional elements and large-scale epigenetic phenomena in the developing zebrafish Zebrafish, a popular model for embryonic development and for modelling human diseases, has so far lacked a systematic functional annotation programme akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created the first central repository to store and process zebrafish developmental functional genomic data. Our Data Coordination Center (https://danio-code.zfin.org) combines a total of 1,802 sets of unpublished and reanalysed published genomics data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements in development, including novel classes with distinct features dependent on their activity in time and space. We delineated the distinction between regulatory elements active during zygotic genome activation and those active during organogenesis, identifying new aspects of how they relate to each other. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predict functional relationships between them beyond sequence similarity, extending the utility of zebrafish developmental genomics to mammals. Arabidopsis AtPME2 has a pH-dependent processivity and control cell wall mechanical properties Pectin methylesterases (PMEs) modify homogalacturonans (HG) chemistry and play a key role in regulating primary cell wall mechanical properties. How PME activity can fine-tune pectin structure in the growing plant has remained elusive. Here we report on the Arabidopsis AtPME2, which we found to be highly expressed during lateral root emergence and dark-grown hypocotyl elongation. We produced the mature active enzyme using heterologous expression in Pichia pastoris and characterized it through the use of a generic plant PME antiserum suitable for detecting recombinant and native enzyme independent of species source. At neutral pH AtPME2 is preferentially active on pectins with a degree of 55-70% methylesterification and can be inhibited by PME inhibitor protein (PMEI). We show that the mode of action for AtPME2 can switch from full processivity (at pH 8), creating large blocks of unmethylated galacturonic acid, to low processivity (at pH 5) and relate these observations to the differences in electrostatic potential of the protein at acidic and alkaline pH. To assess the role of AtPME2 in development, we characterized two knock-out lines. We show that in the context of acidified apoplast, low-processive demethylesterification by AtPME2 can loosen the cell wall, with consequent increase in cell elongation and etiolated hypocotyl length. Our study brings insights into how the pH-dependent regulation by PME activity could affect pectin structure and associated cell wall mechanical properties in expansion.

One sentence summaryThe processivity of AtPME2, a pectin methylesterase that fine-tunes cell wall pectins is modulated by pH in vitro and impacts the mechanical properties of the wall, affecting development in planta. Simpati: patient classifier identifies signature pathways based on similarity networks for the disease prediction BACKGROUNDPathway-based patient classification is a supervised learning task which implies a model learning pathways as features to predict the classes of patients. The counterpart of enrichment tools for the pathway analysis are fundamental methods for clinicians and biomedical scientists. They allow to find signature cellular functions which help to define and annotate a disease phenotype. They provide results which lead human experts to manually classify patients. It is a paradox that pathwaybased classifiers which natively resolve this objective are not strongly developed. They could simulate the human way of thinking, decipher hidden multivariate relationships between the deregulated pathways and the disease phenotype, and provide more information than a probability value. Instead, there are currently only two classifiers of such kind, they require a nontrivial hyperparameter tuning, are difficult to interpret and lack in providing new insights. There is the need of new classifiers which can provide novel perspectives about pathways, be easy to apply with different biological omics and produce new data enabling a further analysis of the patients.

RESULTSWe propose Simpati, an innovative and interpretable patient classifier based on pathway-specific patient similarity networks. The first classifier to adopt ad-hoc novel algorithms for such graph type. It standardizes the biological high-throughput dataset of patients profiles with a propagation algorithm that considers the interconnected nature of the cells molecules for inferring a new activity score. This allows Simpati to classify with dense, sparse, and non-homogenous omic data. Simpati organizes patients molecules in pathways represented by patient similarity networks for being interpretable, handling missing data and preserving the patient privacy. A network represents patients as nodes and a novel similarity measure determines how much every pair act co-ordinately in a pathway. Simpati detects signature biological processes based on how much the topological properties of the related networks separate the patient classes. In this step, it includes a new cohesive subgroup detection algorithm to handle patients not showing the same pathway activity as the other class members. An unknown patient is then classified by a unique recommender system which considers how much is similar to known patients and distant from being an outlier. Simpati outperforms previously published classifiers on five cancer datasets described with two biological omics, classifies well with sparse data, identifies more relevant pathways associated to the patients disease than the competitors and has the lowest computational requirements.

CONCLUSIONSimpati can serve as generic-purpose pathway-based classifier of patient classes. It provides signature pathways to unveil the altered biological mechanisms of a disease phenotype and to classify patients according to the learnt pathway-specific similarities. The signature condition and patient prediction can be deciphered considering the patient similarity networks which must reveal the members of a patient class more cohesive and similar than the non-members. Simpati divides the pathways in up and downinvolved. Upinvolved when the signaling cascades generated by the altered molecules of the disease patients impact stronger the pathway than the ones of the control class. We provide an R implementation, a graphical user interface and a visualization function for the patient similarity networks. The software is available at: https://github.com/LucaGiudice/Simpati Multiple origins of nematode-Wolbachia symbiosis in supergroup F and convergent loss of bacterioferritin in filarial Wolbachia The intracellular endosymbiotic proteobacteria Wolbachia have evolved across the phyla nematoda and arthropoda. In Wolbachia phylogeny, supergroup F is the only clade with members from both arthropod and filarial nematode hosts and therefore can provide unique insights into their evolution and biology. In this study, 4 new supergroup F Wolbachia genomes have been assembled using a metagenomic assembly and binning approach, wMoz and wMpe from the human filarial parasites Mansonella ozzardi and Mansonella perstans, and wOcae and wMoviF from the blue mason bee Osmia caerulescens and the sheep ked Melophagus ovinus respectively. A comprehensive phylogenomic analysis revealed two independent origins of filarial Wolbachia in supergroup F, likely from ancestral arthropod hosts, as well as a convergent pseudogenization and loss of the bacterioferritin gene across all filarial Wolbachia. The new genomes provide a valuable resource for further studies on symbiosis, evolution, and the discovery of new antibiotics to treat mansonellosis.

Data depositionRaw Illumina reads for the M. perstans isolate Mpe1 and the M. ozzardi isolate Moz2 have been deposited in the NCBI SRA database under BioProject accession numbers PRJNA666672 and PRJNA666671 respectively. The assembled genome and PGAP annotations for the Wolbachia wMpe and wMoz are available from the NCBI GenBank database with accession numbers JACZHU000000000 and JADAKK000000000 respectively. The annotated genomes for wOcae and wMoviF are available as GenBank accessions JAHDJT000000000 and JAHRBA000000000 respectively, deposited under BioProject accessions PRJNA244005 and PRJNA739984 respectively.