Dynamic trajectories of connectome state transitions are heritable The brains functional connectome is dynamic, constantly reconfiguring in an individual-specific manner. However, which characteristics of such reconfigurations are subject to genetic effects, and to what extent, is largely unknown. Here, we identified heritable dynamic features, quantified their heritability, and determined their association with cognitive phenotypes. In resting-state fMRI, we obtained multivariate features, each describing a temporal or spatial characteristic of connectome dynamics jointly over a set of connectome states. We found strong evidence for heritability of temporal features, particularly fractional occupancy (FO) and transition probability (TP), describing the trajectory of state transitions. Genetic effects explained a substantial proportion of phenotypic variance of these features (h2=.39, 95% CI= [.24.54] for FO; h2=.43, 95% CI=[.29.57] for TP). Moreover, these temporal phenotypes were associated with cognitive performance. Contrarily, we found no robust evidence for heritability of spatial features of the dynamic states (i.e., states Modularity and connectivity pattern). Genetic effects may therefore primarily contribute to how the connectome transitions across states, rather than the precise spatial instantiation of the states in individuals. In sum, genetic effects impact the duration spent in each connectivity configuration and the frequency of shifting between configurations, and such temporal features may act as endophenotypes for cognitive abilities.
State-dependent network interactions differentially gate sensory input at the motor and command neuron level in Caenorhabditis elegans Neural responses are influenced by both external stimuli and internal network states. While network states have been linked to behavioral and stimulus states, little is known about how sensory inputs are filtered by whole-brain activity to affect motor and command neurons. Here, we recorded whole-brain activity of Caenorhabditis elegans experiencing bacterial food stimuli, and modeled how sensory inputs affect motor and command neurons in a network state-dependent manner. First, we classified active neurons into six functional clusters: two sensory neuron clusters (ON, OFF), and four motor/command neuron clusters (AVA, RME, SMDD, SMDV). Using encoding models, we found that ON and OFF sensory neurons that respond to onset and removal of bacteria, respectively, employ different adaptation strategies. Next, we used decoding models to show that bacterial onset and removal differentially drive AVA and RME cluster activity. To explore state-dependent effects on AVA and RME clusters, we developed a model that identified network states and fitted submodels for each state to predict how each of the six functional clusters drive AVA and RME cluster activity. We also identified network states in which AVA and RME clusters were either largely unperturbed by or receptive to bacterial sensory input. Furthermore, this model allowed us to disentangle the state-dependent contributions of stimulus timescales and bacterial content to neural activity. Collectively, we present an interpretable approach for modeling network dynamics that goes beyond implication of neurons in particular states, and moves toward explicitly dissecting how neural populations work together to produce state dependence.

Significance StatementA major function of the brain is to transform sensory information into behavior. As the first receiver of sensory input, sensory neuron activity is often most correlated with stimulus features. However, this high-fidelity representation of sensory input becomes diluted as it travels to downstream neurons, where sensory information is integrated with network activity. By the time sensory information reaches motor neurons, it is often difficult to dissociate the influence of sensory input from the influence of network activity. Here, we describe a method that is fully interpretable such that we can show how neural populations on a whole-brain scale interact to produce network states. From there, we can attribute motor neuron activity to network history and sensory input.
Intramolecular 13C signals inform about carbon allocation in plants Stable isotopes at natural abundance are key tools to study physiological processes occurring outside the temporal scope of manipulation and monitoring experiments. Whole-molecule carbon isotope ratios (13C/12C) enable assessments of plant carbon uptake yet conceal information about carbon allocation. Here, we identify an intramolecular 13C/12C signal at treering glucose C-5 and C-6 and develop experimentally testable theories on its origin. More specifically, we assess the potential of processes within C3 metabolism for signal introduction based (inter alia) on constraints on signal propagation posed by metabolic networks. We propose that the intramolecular signal reports carbon allocation into major metabolic pathways in actively photosynthesising leaf cells including the anaplerotic, shikimate, and non-mevalonate pathway. We support our theoretical framework by linking it to previously reported whole-molecule 13C/12C increases in cellulose of ozone-treated Betula pendula and a highly significant relationship between the intramolecular signal and tropospheric ozone concentration. Our theory postulates a pronounced preference of leaf-cytosolic triose-phosphate isomerase to catalyse the forward reaction in vivo (dihydroxyacetone phosphate to glyceraldehyde 3-phosphate). In conclusion, intramolecular 13C/12C analysis resolves information about carbon uptake and allocation enabling more comprehensive assessments of carbon metabolism than whole-molecule 13C/12C analysis.

HighlightIntramolecular 13C/12C analysis resolves information about carbon uptake and allocation (and associated environmental controls) enabling more comprehensive assessments of carbon metabolism, plant-environment interactions, and environmental variability than whole-molecule 13C/12C analysis. RSC and GRFs confer promoter directionality by limiting divergent noncoding transcription The directionality of gene promoters - the ratio of protein-coding over divergent noncoding transcription - is highly variable and regulated. How promoter directionality is controlled remains poorly understood. Here, we show that the chromatin remodelling complex RSC and general regulatory factors (GRFs) dictate promoter directionality by attenuating divergent transcription. At gene promoters that are highly directional, depletion of RSC leads to a relative increase in divergent noncoding transcription and thus a decrease in promoter directionality. We find that RSC facilitates nucleosome positioning upstream in promoters at the sites of divergent transcription. These highly directional promoters are also enriched for the binding of GRFs such as Reb1 and Abf1. Ectopic targeting of divergent transcription initiation sites with GRFs or the dCas9 protein suppresses divergent transcription. Our data suggest that RSC and GRFs play a pervasive role in limiting divergent transcription. We propose that any DNA binding factor, when stably associated with cryptic transcription start sites, form barriers for repressing divergent transcription. Our study provides an explanation as to why certain promoters are more directional than others. Creation of a single cell RNA seq meta-atlas to define human liver immune homeostasis The liver is unique in both its ability to maintain immune homeostasis and in its potential for immune tolerance following solid organ transplantation. Single-cell RNA sequencing (scRNA seq) is a powerful approach to generate highly dimensional transcriptome data to understand cellular phenotypes. However, when scRNA data is produced by different groups, with different data models, different standards, and samples processed in different ways, it can be challenging to draw meaningful conclusions from the aggregated data. The goal of this study was to establish a method to combine  human liver scRNA seq datasets by 1) characterizing the heterogeneity between studies and 2) using the meta-atlas to define the dominant phenotypes across immune cell subpopulations in healthy human liver. Publicly available scRNA seq data generated from liver samples obtained from a combined total of 17 patients and [~]32,000 cells were analyzed. Liver-specific immune cells (CD45+) were extracted from each dataset, and immune cell subpopulations (myeloid cells, NK and T cells, plasma cells, and B cells) were examined using dimensionality reduction (UMAP), differential gene expression, and ingenuity pathway analysis. All datasets co-clustered, but cell proportions differed between studies. Gene expression correlation demonstrated similarity across all studies, and canonical pathways that differed between datasets were related to cell stress and oxidative phosphorylation rather than immune-related function. Next, a meta-atlas was generated via data integration and compared against PBMC data to define gene signatures for each hepatic immune subpopulation. This analysis defined key features of hepatic immune homeostasis, with decreased expression across immunologic pathways and enhancement of pathways involved with cell death. This method for meta-analysis of scRNA seq data provides a novel approach to broadly define the features of human liver immune homeostasis. Specific pathways and cellular phenotypes described in this human liver immune meta-atlas provide a critical reference point for further study of immune mediated disease processes within the liver. A nesprin-4/kinesin-1 cargo model for nuclear positioning in cochlear outer hair cells Nuclear positioning is important for the functionality of many cell types and is mediated by interactions of cytoskeletal elements and nucleoskeleton proteins. Nesprin proteins, part of the linker of nucleoskeleton and cytoskeleton (LINC) complex, have been shown to participate in nuclear positioning in multiple cell types. Outer hair cells (OHCs) in the inner ear are specialized sensory epithelial cells that utilize somatic electromotility to amplify auditory signals in the cochlea. Recently, nesprin-4 (encoded by Syne4) was shown to play a crucial role in nuclear positioning in OHCs. Syne4 deficiency in humans and mice leads to mislocalization of the OHC nuclei and cell death resulting in deafness. However, it is unknown how nesprin-4 mediates the position of the nucleus, and which other molecular components are involved in this process. Here, we show that the interaction of nesprin-4 and the microtubule motor kinesin-1 is mediated by a conserved 4 amino-acid motif. Using in-vivo AAV gene delivery, we show that this interaction is critical for nuclear positioning and hearing in mice. Nuclear mislocalization and cell death of OHCs coincide with the onset of hearing and electromotility and are solely restricted to outer, but not inner, hair cells. Likewise, the C. elegans functional homolog of nesprin-4, UNC-83, uses a similar motif to mediate interactions between migrating nuclei and kinesin-1. Overall, our results suggest that OHCs require unique cellular machinery for proper nuclear positioning at the onset of electromotility. This machinery relies on the interaction between nesprin-4 and kinesin-1 motors supporting a microtubule cargo model for nuclear positioning. Thermal stimulation of pressure phosphenes To investigate effect on the intensity of pressure phosphenes (PP) of various methods of heating hands, as well as manual influence on cervical spine. This study included subjective assessments of the PP intensity in 10 healthy men, and chronometry of bioelectrical activity of brain and heart using electroencephalography (EEG) and electrocardiography (ECG). EEG and ECG frequency spectra respond synchronously to pressure, heat and light. The time of excitation of potential oscillations in the visual area of the cortex of both hemispheres is consistent with the delay in the onset of PP flaring. The effect of PP stimulation is enhanced when the hands and body are heated in a humid environment and at temperatures above 42 {degrees}C. Stimulation of PP by heating indicates the convergence of impulses from neurons of the lateral geniculate body (LGB) and nuclei of the thalamus, which are responsible for thermoreceptors in the skin of the palms and fingers. The thermal stimulation effect of PP is enhanced at temperatures above 42 {degrees}C due to the fact that thermoreceptors work as pain receptors. The mechanism of PP generation is dominated by the processes of redistribution and recombination of charges in the layers of the retina and LGB. Symbiodiniaceae conduct under natural bleaching stress during advanced gametogenesis stages of a mesophotic coral The mesophotic coral Alveopora allingi from the northern Gulf of Eilat/Aqaba, Red Sea, is affected by year-round partial coral-bleaching events. During these events, the migration of Symbiodiniaceae takes place from the coral-host mesoglea to the developed oocytes in bleached parts of colonies of A. allingi but not in the non-bleached parts. Additionally, these oocytes are abnormal, missing part of the structural material of the peripheral areas and are also significantly larger in the bleached areas of the colonies. Hence, we suggest a parasitic behavior of the symbionts or a commensalism relationship which enhance symbionts needs during bleaching periods and may boost the gametogenesis development in these corals. We propose that evolutionarily, this behavior may greatly contribute to the symbiont community survival throughout the bleaching period, and it can also be beneficial for the hosts persistence and adaptation to bleaching through the acquisition of a specific symbiont community following the bleaching event. The pH Dependence of Niclosamide Solubility, Dissolution, and Morphology Motivates Potentially Universal Mucin-Penetrating Nasal and Throat Sprays for COVID19, its Contagious Variants, and Other Respiratory Viral Infections MotivationWith the coronavirus pandemic still raging, prophylactic nasal and early treatment throat sprays could help prevent infection and reduce viral load. Niclosamide has the potential to treat a broad range of viral infections if local bioavailability is optimized as mucin-penetrating solutions as opposed to more traditional microparticle-based sprays that cannot penetrate the mucin.

ExperimentalpH-dependence of supernatant concentrations and dissolution rates of niclosamide were measured in buffered solutions by Nanodrop-UV/Vis-spectroscopy for niclosamide from different suppliers, as precipitated material, and as cosolvates. Data was compared to predictions from Henderson Hasselbalch and precipitation pH models. Optimal microscopy was used to observe the morphologies of precipitated and converted niclosamide.

ResultsSupernatant-concentrations of niclosamide increased with increasing pH, from 1.77uM at pH 3.66 to 30uM at pH 8, and more rapidly from 90uM at pH8.5 to 300uM at pH 9.1, reaching 641uM at pH 9.5. Logarithmic rates for dissolution increased by [~]3x for pHs 8.62 to 9.44. However, when precipitated from supersaturated solution, niclosamide equilibrated to much lower final supernatant concentrations, reflective of more stable polymorphs at each pH that were also apparent for niclosamide from other suppliers and cosolvates.

ConclusionsGiven niclosamides activity against COVID19, its more contagious variants, and other respiratory viral infections, these niclosamide solutions, that put the virus in lockdown, could represent universal prophylactic nasal and early treatment throat sprays. As solutions they would be the simplest and potentially most effective formulations from both an efficacy standpoint as well as manufacturing and distribution, with no cold chain. They now just need testing. Barcode aware adaptive sampling for GridION and PromethION Oxford Nanopore sequencers. Adaptive sampling enables selection of individual DNA molecules from sequencing libraries, a unique property of nanopore sequencing. Here we develop our adaptive sampling tool readfish to become "barcode-aware" enabling selection of different targets within barcoded samples or filtering out individual barcodes. We show that multiple human genomes can be assessed for copy number and structural variation on a single sequencing flow cell using sample specific customised target panels on both GridION and PromethION devices. Wolbachia load variation in Drosophila is more likely caused by drift than by host genetic factors Symbiosis is a continuum of long-term interactions ranging from mutualism to parasitism, according to the balance between costs and benefits for the protagonists. The density of endosymbionts is, in both cases, a key factor that determines both the transmission of symbionts and the host extended phenotype, and is thus tightly regulated within hosts. However, the evolutionary and molecular mechanisms underlying bacterial density regulation are currently poorly understood. In this context, the symbiosis between the fruit fly and its intracellular bacteria Wolbachia (wMelPop strain) is particularly interesting to study. Although vertically transmitted, the symbiont is pathogenic, and a positive correlation between virulence and wMelPop density is observed. In addition, the number of repeats of a bacterial genomic region -Octomom-is positively correlated with Wolbachia density, underlying a potential genetic mechanism that controls bacterial density. Interestingly, the number of repeats varies between host individuals, but most likely also within them. Such genetic heterogeneity within the host could promote conflicts between bacteria themselves and with the host, notably by increasing within-host competition between symbiont genotypes through a process analogous to the tragedy of the commons. To characterize the determinisms at play in the regulation of bacterial density, we first introgressed wMelPop in different genetic backgrounds of D. melanogaster, and found different density levels and Octomom copy numbers in each host lineage. To determine whether such variations reflect a host genetic determinism on density regulation through Octomom copy number selection, we replicated the introgressions and performed reciprocal crosses on the two Drosophila populations with the most extreme density levels. In both experiments, we detected an absence of directionality in the patterns of infection, associated with a strong instability of these patterns across generations. Given that bacterial density was highly correlated with Octomom copy numbers in all experiments, these results rather suggest a strong influence of drift and a random increase in the frequency of certain bacterial variants. We then discuss how drift, both on the symbiont population during transmission and on the host population, could limit the efficiency of selection in such a symbiotic system, and the consequences of drift on the regulation of density and composition of bacterial populations. Task-independent metrics of computational hardness predict human cognitive performance The survival of human organisms depends on our ability to solve complex tasks in the face of limited cognitive resources. However, little is known about the factors that drive the complexity of those tasks. Here, building on insights from computational complexity theory, we quantify the computational hardness of cognitive tasks using a set of task-independent metrics related to the computational resource requirements of individual instances of a task. We then examine the relation between those metrics and human behavior and find that they predict both time spent on a task as well as accuracy in three canonical cognitive tasks. Our findings demonstrate that performance in cognitive tasks can be predicted based on generic metrics of their inherent computational hardness.

TeaserThe ability of humans to solve cognitive tasks is affected by generic mathematical properties of problems related to their computational complexity. Integrative analysis reveals aged clonal B cells, microenvironment and c-Myc activation in the origin of age-related lymphoma While cancer is an age-related disease, most studies focus on genetically engineered younger mouse models. Here, we uncover how cancer develops as a consequence of the naturally aged immune system in mice. B-cell lymphoma frequently occurs in aged mice and is associated with increased cell size, splenomegaly, and a novel clonal B-cell population. Age-emergent B cells clonally expand outside of germinal centers driven by somatic mutations, activated c-Myc and hypermethylated promoters, and both genetically and epigenetically recapitulate human follicular and diffuse-large B-cell lymphomas. Mechanistically, mouse cancerous B cells originate from age-associated B cells, which are atypical memory B cells. Age-associated B cells secrete a spectrum of proinflammatory cytokines and activate paracrinally the expression of c-Myc in surrounding B cells. Although clonal B cells are a product of an aging microenvironment, they evolve being self-sufficient and support malignancy when transferred into young mice. Inhibition of mTOR and c-Myc attenuates premalignant changes in B cells during aging and emerges as a therapeutic strategy to delay the onset of age-related lymphoma. Together, we uncover how aging generates cancerous B cells, characterize a model that captures the origin of spontaneous cancer during aging and identify interventions that may postpone age-associated lymphoma. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. A Histoplasma capsulatum lipid metabolic map identifies antifungal targets Lipids play a fundamental role in fungal cell biology, being essential cell membrane components and major targets of antifungal drugs. A deeper knowledge of lipid metabolism is key for developing new drugs and a better understanding of fungal pathogenesis. Here we built a comprehensive map of the Histoplasma capsulatum lipid metabolic pathway by incorporating proteomic and lipidomic analyses. We performed genetic complementation and overexpression of H. capsulatum genes in Saccharomyces cerevisiae to validate reactions identified in the map and to determine enzymes responsible for catalyzing orphan reactions. The map led to the identification of both the fatty acid desaturation and the sphingolipid biosynthesis pathways as targets for drug development. We found that the sphingolipid biosynthesis inhibitor myriocin, the fatty acid desaturase inhibitor thiocarlide and the fatty acid analog 10-thiastearic acid inhibit H. capsulatum growth in nanomolar to low micromolar concentrations. These compounds also reduced the intracellular infection in an alveolar macrophage cell line. Overall, this lipid metabolic map revealed pathways that can be targeted for drug development. The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. Supervised Learning Model Predicts Protein Adsorption to Carbon Nanotubes Engineered nanoparticles are advantageous for numerous biotechnology applications, including biomolecular sensing and delivery. However, testing the compatibility and function of nanotechnologies in biological systems requires a heuristic approach, where unpredictable biofouling via protein corona formation often prevents effective implementation. Moreover, rational design of biomolecule-nanoparticle conjugates requires prior knowledge of such interactions or extensive experimental testing. Toward better applying engineered nanoparticles in biological systems, herein, we develop a random forest classifier (RFC) trained with proteomic mass spectrometry data that identifies proteins that adsorb to nanoparticles, based solely on the proteins amino acid sequence. We model proteins that populate the corona of a single-walled carbon nanotube (SWCNT)-based optical nanosensor and study whether there is a relationship between the proteins amino acid-based properties and the proteins adsorption to SWCNTs. We optimize the classifier and characterize the classifier performance against other models. To evaluate the predictive power of our model, we apply the classifier to rapidly identify proteins with high binding affinity to SWCNTs, followed by experimental validation. We further determine protein features associated with increased likelihood of SWCNT binding: high content of solvent-exposed glycine residues and non-secondary structure-associated amino acids. Conversely, proteins with high content of leucine residues and beta-sheet-associated amino acids are less likely to form the SWCNT protein corona. The classifier presented herein provides a step toward undertaking the otherwise intractable problem of predicting protein-nanoparticle interactions, which is needed for more rapid and effective translation of nanobiotechnologies from in vitro synthesis to in vivo use.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=66 SRC="FIGDIR/small/449132v2_ufig1.gif" ALT="Figure 1">
View larger version (22K):
org.highwire.dtl.DTLVardef@bbc175org.highwire.dtl.DTLVardef@9a0b94org.highwire.dtl.DTLVardef@16e22e1org.highwire.dtl.DTLVardef@1b40a84_HPS_FORMAT_FIGEXP  M_FIG C_FIG Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less-functional Misfolded protein conformations with decreased functionality can bypass the proteostasis machinery and remain soluble in vivo. This is an unexpected phenomenon as several cellular quality control mechanisms have evolved to rid cells of misfolded proteins. Three questions, then, are: how is it structurally possible for long-lived, soluble, misfolded proteins to bypass the proteostasis machinery and processes? How widespread are these soluble, misfolded states across the proteome? And how long do they persist for? Here, we address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of all proteins exhibit subpopulations of misfolded conformations that are likely to bypass molecular chaperones, avoid aggregation, and not be rapidly degraded. These misfolded states may persist for months or longer for some proteins. Structurally characterizing these misfolded states, we observe they have a large amount of native structure, but also contain localized misfolded regions from non-native changes in entanglement, in which a protein segment threads through a loop formed by another portion of the protein that is not found in the native state. The surface properties of these misfolded states are native like, suggesting they may bypass the proteostasis machinery and its regulatory processes to remain soluble, while their entanglements make these states long-lived kinetic traps, as disentanglement requires unfolding of already folded portions of the protein. In terms of function, we predict that one-third of proteins have subpopulations that misfold into less-functional states that have structurally perturbed functional sites yet remain soluble. Data from limited-proteolysis mass spectrometry experiments, which interrogate the misfolded conformations populated by proteins upon unfolding and refolding, are consistent with the structural changes seen in the entangled states of glycerol-3-phosphate dehydrogenase upon misfolding. These results provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass cellular quality controls, and indicate, unexpectedly, this may be a wide-spread phenomenon in proteomes. Such entanglements are observed in many native structures, suggesting the non-native entanglements we observe are plausible. More broadly, these near-native entangled structures suggest a hypothesis for how synonymous mutations can modulate downstream protein structure and function, with these mutations partitioning nascent proteins between these kinetically trapped states. Bayesian Additive Regression Trees for Genotype by Environment Interaction Models We propose a new class of models for the estimation of genotype by environment (GxE) interactions in plant-based genetics. Our approach, named AMBARTI, uses semi-parametric Bayesian additive regression trees to accurately capture marginal genotypic and environment effects along with their interaction in a fully Bayesian model. We demonstrate that our approach is competitive or superior to similar models widely used in the literature via both simulation and a real world data set. Furthermore, we introduce new types of visualisation to properly assess both the marginal and interactive predictions from the model. An R package that implements our approach is available at https://github.com/ebprado/ambarti. 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. wnt16 exerts pleiotropic effects on bone and lean mass in zebrafish Osteoporosis, a disease of bone fragility, and sarcopenia, a condition of reduced muscle mass and strength, commonly occur in the same individual--a condition termed osteosarcopenia. Genetic variants at the CPED1-WNT16 locus have pleiotropic effects on human bone mineral density (BMD) and lean tissue mass. However, our current understanding of Wnt16 and its influence on bone fails to explain effects on lean mass at this locus. In this study, we sought to better understand the biological mechanism underlying pleiotropy at the CPED1-WNT16 locus on BMD and lean mass. Using single cell analysis, microCT imaging, and genetic approaches in zebrafish, we show that wnt16 exerts dual influence on muscle and bone, and provide evidence that this occurs through a function in somite development, a period in which precursors for both of these tissues are specified. We show that mutations in wnt16 alter vertebral size and lean mass during axial skeletogenesis. We show that wnt16+ cells transcriptionally resemble cells within the amniote dermomyotome, and provide evidence they function as myogenic precursors during embryonic development. We also provide evidence that WNT16 is a gene of major effect on lean mass at the CPED-WNT16 locus. Our studies support the potential for wnt16 to exert pleiotropic effects on bone and lean mass through a function in myogenic precursors. Notch signaling via Hey1 and Id2b regulates Muller glia's regenerative response to retinal injury Zebrafish Muller glia (MG) respond to retinal injury by suppressing Notch signaling and producing progenitors for retinal repair. A certain threshold of injury-derived signal must be exceeded in order to engage MG in a regenerative response (MGs injury-response threshold). Pan-retinal Notch inhibition expands the zone of injury-responsive MG at the site of focal injury, suggesting that Notch signaling regulates MGs injury-response threshold. We found that Notch signaling enhanced chromatin accessibility and gene expression at a subset of regeneration-associated genes in the uninjured retina. Two Notch effector genes, hey1 and id2b, were identified that reflect bifurcation of the Notch signaling pathway, and differentially regulate MGs injury-response threshold and proliferation of MG-derived progenitors. Furthermore, Notch signaling component gene repression in the injured retina suggests a role for Dll4, Dlb, and Notch3 in regulating Notch signaling in MG and epistasis experiments confirm that the Dll4/Dlb-Notch3-Hey1/Id2b signaling pathway regulates MGs injury-response threshold and proliferation. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. Behavioral and postural analyses establish sleep-like states for mosquitoes that can impact host landing and blood feeding Sleep is an evolutionarily conserved process that has been described in different animal systems. For insects, sleep characterization has been primarily achieved using behavioral and electrophysiological correlates in a few systems. Sleep in mosquitoes, which are important vectors of disease-causing pathogens, has not been directly examined. This is surprising as circadian rhythms, which have been well studied in mosquitoes, influence sleep in other systems. In this study, we characterized sleep in mosquitoes using body posture analysis and behavioral correlates, and quantified the effect of sleep deprivation on sleep rebound, host landing and blood-feeding propensity. Body and appendage position metrics revealed a clear distinction between the posture of mosquitoes in their putative sleep and awake states for multiple species, which correlate with a reduction in responsiveness to host cues. Sleep assessment informed by these posture analyses indicated significantly more sleep during periods of low activity. Nighttime and daytime sleep deprivation resulting from the delivery of vibration stimuli induced sleep rebound in the subsequent phase in day and night active mosquitoes, respectively. Lastly, sleep deprivation suppressed host landing in both laboratory and field settings and also impaired blood feeding of a human host when mosquitoes would normally be active. These results suggest that quantifiable sleep states occur in mosquitoes, and highlight the potential epidemiological importance of mosquito sleep.