title: A hybrid capture RNA bait set for resolving genetic and evolutionary relationships in angiosperms from deep phylogeny to intraspecific lineage hybridization

abstract: Novel multi-gene targeted capture probes have been developed with the objective of obtaining multi-locus high quality sequence reads across any angiosperm lineage. Using existing genomic and transcriptomic data, two independent single assay probe/bait sets have been developed, the first targeting conserved exons from 20 low copy nuclear genes (OzBaits_NR V1.0) and the second, 19 plastid gene regions (OZBaits_CP V1.0). These  universal bait sets can efficiently generate DNA sequence data that are suitable for systematics and evolutionary studies of flowering plants. The bait sets can be ordered as Daicel-Arbor Sciences custom myBaits. We demonstrate the utility of the bait set in consistently recovering the targeted genomic regions across an evolutionarily broad range of angiosperm taxa.


title: Open-Access STRS Database Of Populations From The 1000 Genomes Project Using High Coverage Phase 3 Data

abstract: Accurate STR genotyping from next-generation sequencing (NGS) data has been challenging. Haplotype inference and phasing for STRs (HipSTR) was specifically developed to deal with genotyping errors and obtain reliable STR genotypes from whole-genome sequencing datasets. The objective of this investigation was to perform a comprehensive genotyping analysis of a set of STRs of broad forensic interest from the 1000 Genomes populations and release a reliable open-access STR database to the forensic genetics community. A set of 22 STR markers were analyzed using the CRAM files of the 1000 Genomes Project Phase 3 high-coverage (30x) dataset generated by the New York Genome Center (NYGC). HipSTR was used to call genotypes from 2,504 samples from 26 populations organized into five groups: African, East Asian, European, South Asian, and admixed American. The D21S11 marker could not be detected in the present study. Moreover, the Hardy-Weinberg equilibrium analysis, coupled with a comprehensive analysis of allele frequencies, revealed that HipSTR could not identify longer Penta E (and Penta D at a lesser extent) alleles. This issue is probably due to the limited length of sequencing reads available for genotype calling, resulting in heterozygote deficiency. Notwithstanding that, AMOVA, a clustering analysis using STRUCTURE, and a Principal Coordinates Analysis revealed a clear-cut separation between the four major ancestries sampled by the 1000 Genomes Consortium (AFR, EUR, EAS, SAS). Meanwhile, the AMOVA results corroborated previous reports that most of the variance is (97.12%) observed within populations. This set of analyses revealed that except for larger Penta D and Penta E alleles, allele frequencies and genotypes defined by HipSTR from the 1000 Genomes Project phase 3 data and offered as an open-access database are consistent and highly reliable.

category: genetics


title: Fuzzy Integral With Particle Swarm Optimization for CNN-Based Motor-Imagery EEG Classification

abstract: Recently, decoding human electroencephalographic (EEG) data using convolutional neural network (CNN) has driven the state-of-the-art recognition of motor-imagery EEG patterns for brain-computer interfacing (BCI). While a variety of CNN models have been used to classify motor-imagery EEG data, it is unclear if aggregating an ensemble of heterogeneous CNN models could further enhance the classification performance. To integrate the outputs of ensemble classifiers, this work utilizes fuzzy integral with particle swarm optimization (PSO) to estimate optimal confidence levels assigned to classifiers. The proposed frame-work aggregates CNN classifiers and fuzzy integral with PSO, achieving robust performance in single-trial classification of motor-imagery EEG data across various CNN model training schemes depending on the scenarios of BCI usage. This proof-of-concept study demonstrates the feasibility of applying fuzzy fusion techniques to enhance CNN-based EEG decoding and benefit practical applications of BCI.

category: bioengineering


title: Effective expression analysis using gene interaction matrices and convolutional neural networks

abstract: Artificial intelligence recently experienced a renaissance with the advancement of convolutional neural networks (CNNs). CNNs require spatially meaningful matrices (e.g., image data) with recurring patterns, limiting its applicability to high-throughput omics data. We present GIM, a simple, CNN-ready framework for omics data to detect both individual and network-level entities of biological importance. Using gene expression data, we show that GIM-CNNs can outperform comparable neural networks in performance and their design facilitates network-level interpretability. GIM-CNNs provide a means to discover novel disease-relevant factors beyond individual genes and their expression, factors that are likely missed by standard differential gene expression approaches.

category: bioinformatics


title: GMQN: A reference-based method for correcting batch effects as well as probes bias in HumanMethylation BeadChip

abstract: Illumina HumanMethylation BeadChip is one of the most cost-effective ways to quantify DNA methylation levels at the single-base level across the human genome, which makes it a routine platform for epigenome-wide association studies. It has accumulated tens of thousands of DNA methylation array samples in public databases, thus provide great support for data integration and further analysis. However, majority of public DNA methylation data are deposited as processed data without background probes which are widely used in data normalization. Here we present Gaussian mixture quantile normalization (GMQN), a reference based method for correcting batch effects as well as probes bias in HumanMethylation BeadChip. Availability and implementation: https://github.com/MengweiLi-project/gmqn.


title: Hi-C Contacts Encode Heterogeneity in Sub-diffusive Motion of E. coli Chromosomal Loci

abstract: Underneath its apparently simple architecture, the circular chromosome of E. coli is known for displaying complex dynamics in its cytoplasm. Recent experiments have hinted at an inherently heterogeneous dynamics of chromosomal loci, the origin of which has largely been elusive. In this regard, here we investigate the loci dynamics of E. coli chromosome in a minimally growing condition at 30{degrees}C by integrating the experimentally derived Hi-C interaction matrix within a computer model. Our quantitative analysis demonstrates that, while the dynamics of the chromosome is sub-diffusive in a viscoelastic media in general, the diffusion constants and the diffusive exponents are strongly dependent on the spatial coordinates of chromosomal loci. In particular, the loci in Ter Macro-domain display slower mobility compared to the others. The result is found to be robust even in the presence of active noise. Interestingly, a series of control investigations reveal that the absence of Hi-C interactions in the model would have abolished the heterogeneity in loci diffusion, indicating that the observed coordinate-dependent chromosome dynamics is heavily dictated via Hi-C-guided longrange inter-loci communications. Overall, the study underscores the key role of Hi-C interactions in guiding the inter-loci encounter and in modulating the underlying heterogeneity of the loci diffusion.

category: biophysics


title: A G-lectin Receptor Kinase is a Negative Regulator of Arabidopsis Immunity Against Root-Knot Nematode Meloidogyne incognita

abstract:
Root-knot nematodes (Meloidogyne spp., RKN) are responsible for extensive crop losses worldwide. For infection, they penetrate plant roots, migrate between plant cells, and establish feeding sites, known as giant cells, in the root pericycle. Previously, we found that nematode perception and early plant responses were similar to those for microbial pathogens and require the BAK1 co-receptor in Arabidopsis thaliana and tomato. To identify additional receptors involved in this process, we implemented a reverse genetic screen for resistance or sensitivity to RKN using Arabidopsis T-DNA alleles of genes encoding transmembrane receptor-like kinases. This screen identified a pair of allelic mutations with enhanced resistance to RKN in a gene we named ENHANCED RESISTANCE TO NEMATODES 1 (ERN1). ERN1 encodes a G-type lectin receptor kinase (G-LecRK) with a single pass transmembrane domain. Further characterization showed that ern1 mutants displayed stronger activation of MAP kinases, elevated levels of the defense marker MYB51, and enhanced H202 accumulation in roots upon RKN elicitor treatments. Elevated MYB51 expression and ROS burst were also observed in leaves of ern1 mutants upon flg22 treatment. Complementation of ern1.1 with 35S- or native promotor-driven ERN1 rescued the RKN infection and enhanced defense phenotypes. Taken together, our results indicate that ERN1 is an important negative regulator of immunity.

One sentence summaryA plasma membrane localized G-lectin receptor kinase acts as a negative immune regulator by interfering with defense responses activated by nematode and microbial elicitors.

category: plant biology


title: A New Perspective on Individual Reliability beyond Group Effects for Event-related Potentials: A Multisensory Investigation and Computational Modeling

abstract: The dominant approach in investigating the individual reliability for event-related potentials (ERPs) is to extract peak-related features at electrodes showing the strongest group effects. Such a peak-based approach implicitly assumes ERP components showing a stronger group effect are also more reliable, but this assumption has not been substantially validated and few studies have investigated the reliability of ERPs beyond peaks. In this study, we performed a rigorous evaluation of the test-retest reliability of ERPs collected in a multisensory and cognitive experiment from 82 healthy adolescents, each having two sessions. By comparing group effects and individual reliability, we found that a stronger group-level response in ERPs did not guarantee a higher reliability. Further, by simulating ERPs with a computational model, we found that the consistency between group-level ERP responses and individual reliability was modulated by inter-subject latency jitter and inter-trial variability. The current findings suggest that the conventional peak-based approach may underestimate the individual reliability in ERPs. Hence, a comprehensive evaluation of the reliability of ERP measurements should be considered in individual-level neurophysiological trait evaluation and psychiatric disorder diagnosis.

category: neuroscience


title: Linear superposition of responses evoked by individual glottal pulses explain over 80% of the frequency following response to human speech in the macaque monkey

abstract: The frequency-following response (FFR) is a scalp-recorded electrophysiological potential that closely follows the periodicity of complex sounds such as speech. It has been suggested that FFRs reflect the linear superposition of responses that are triggered by the glottal pulse in each cycle of the fundamental frequency (F0 responses) and sequentially propagate through auditory processing stages in brainstem, midbrain, and cortex. However, this conceptualization of the FFR is debated, and it remains unclear if and how well a simple linear superposition can capture the spectro-temporal complexity of FFRs that are generated within the highly recurrent and non-linear auditory system. To address this question, we used a deconvolution approach to compute the hypothetical F0 responses that best explain the FFRs in rhesus monkeys to human speech and click trains with time-varying pitch patterns. The linear superposition of F0 responses explained well over 90% of the variance of click train steady state FFRs and well over 80% of mandarin tone steady state FFRs. The F0 responses could be measured with high signal-to-noise ratio and featured several spectro-temporally and topographically distinct components that likely reflect the activation of brainstem (<5ms; 200-1000 Hz), midbrain (5-15 ms; 100-250 Hz) and cortex (15-35 ms; ~90 Hz). In summary, our results in the monkey support the notion that FFRs arise as the superposition of F0 responses by showing for the first time that they can capture the bulk of the variance and spectro-temporal complexity of FFRs to human speech with time-varying pitch. These findings identify F0 responses as a potential diagnostic tool that may be useful to reliably link altered FFRs in speech and language disorders to altered F0 responses and thus to specific latencies, frequency bands and ultimately processing stages.


title: Social belonging: Brain structure and function is linked to membership in sports teams, religious groups and social clubs

abstract:
Human behaviour across the life span is driven by the psychological need to belong, from kindergarten to bingo nights. Being part of social groups constitutes a backbone for communal life, and confers many benefits for physical and mental health. Capitalizing on neuroimaging and behavioural data from ~40.000 participants from the UK Biobank population cohort, we used structural and functional analyses to explore how social participation is reflected in the human brain. Across three different types of social groups, structural analyses point towards variance in ventromedial prefrontal cortex, fusiform gyrus and anterior cingulate cortex as structural substrates tightly linked to social participation. Functional connectivity analyses emphasized the importance of default mode and limbic network, but also showed differences for sports teams and religious groups as compared to social clubs.

Taken together, our findings establish the structural and functional integrity of the default mode network as a neural signature of social belonging.


title: Home alone: A population neuroscience investigation of brain morphology substrates

abstract: As a social species, ready exchange with peers is a pivotal asset - our "social capital". Yet, single-person households have come to pervade metropolitan cities worldwide, with unknown consequences in the long run. Here, we systematically explore the morphological manifestations associated with singular living in [~]40,000 UK Biobank participants. The uncovered population-level signature spotlights the highly associative default mode network, in addition to findings such as in the amygdala central, cortical and corticoamygdaloid nuclei groups, as well as the hippocampal fimbria and dentate gyrus. Sex-stratified analyses revealed male-specific neural substrates, including somatomotor, saliency and visual systems, while female-specific neural substrates centred on the dorsomedial prefrontal cortex. In line with our demographic profiling results, the discovered neural imprint of living alone is potentially linked to alcohol and tobacco consumption, anxiety, sleep quality as well as daily TV watching. The secular trend for solitary living will require new answers from public-health decision makers.


title: The birth-and-death evolution of cytochrome P450 genes in bees

abstract:
The birth-and-death model of multigene family evolution describes how families can expand by duplication and contract by gene deletion and formation of pseudogenes. The phylogenetic stability of a gene is thought to be related to the degree of functional importance. However, it is unclear how much evolution of a gene in a gene family is driven by adaptive versus neutral processes. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Bees have a high toxin exposure due to their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some bees. Here, we describe the P450s of the orchid bee Euglossa dilemma. Orchid bees are a neotropical clade in which males form perfume bouquets used in courtship displays by collecting a diverse set of volatile compounds, resulting in high chemical compound exposure. We conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We do not find a relationship between the ecology of a bee species and its P450 repertoire. Our analyses reveal that P450 clades can be classified into stable and unstable clades, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures, with signals of adaptive evolution detected, suggesting that both gene duplication and positive selection driving sequence divergence have played a role in the diversification of bee P450s. Our works highlights the complexity of multigene family evolution which does not always follow generalized predictions.

Significance statementGene family evolution is characterized by deletion and duplication, but it is unclear if this is driven by adaptive or neutral evolutionary processes. To investigate the dynamics of gene family evolution we analyze the P450s of ten bee species, a well-studied gene family involved in both physiological and detoxification roles. We find that genes involved in detoxification are more likely to belong to unstable clades, exhibit dynamic evolutionary pressures, and show signals of adaptive evolution. However, we do not find evidence for a relationship between resin collection or sociality, and the P450 repertoire of a bee species. Our findings highlight the complexity of multigene family evolution which does not always follow generalized predictions.

category: molecular biology


title: Guidelines for reliable and reproducible functional enrichment analysis

abstract: Gene set enrichment tests (a.k.a. functional enrichment analysis) are among the most frequently used methods in computational biology. Despite this popularity, there are concerns that these methods are being applied incorrectly and the results of some peer-reviewed publications are unreliable. These problems include the use of inappropriate background gene lists, lack of false discovery rate correction and lack of methodological detail. An example analysis of public RNA-seq reveals that these methodological errors alter enrichment results dramatically. To ascertain the frequency of these errors in the literature, we performed a screen of 186 open access research articles describing functional enrichment results. We find that 95% of analyses using over-representation tests did not implement an appropriate background gene list or did not describe this in the methods. Failure to perform p-value correction for multiple tests was identified in 43% of analyses. Many studies lacked detail in the methods section about the tools and gene sets used. Only 15% of studies avoided major flaws, which highlights the poor state of functional enrichment rigour and reporting in the contemporary literature. We provide a set of minimum standards that should act as a checklist for researchers and peer-reviewers.


title: Motivations, concerns and selection biases when posting preprints: a survey of bioRxiv authors

abstract: Since 2013, the usage of preprints as a means of sharing research in biology has rapidly grown, in particular via the preprint server bioRxiv. Recent studies have found that journal articles that were previously posted to bioRxiv received a higher number of citations or mentions/shares on other online platforms compared to articles in the same journals that were not posted. However, the exact causal mechanism for this effect has not been established, and may in part be related to authors biases in the selection of articles that are chosen to be posted as preprints. We aimed to investigate this mechanism by conducting a mixed-methods survey of 1,444 authors of bioRxiv preprints, to investigate the reasons that they post or do not post certain articles as preprints, and to make comparisons between articles they choose to post and not post as preprints. We find that authors are most strongly motivated to post preprints to increase awareness of their work and increase the speed of its dissemination; conversely, the strongest reasons for not posting preprints centre around a lack of awareness of preprints and reluctance to publicly post work that has not undergone a peer review process. We additionally find weak evidence that authors preferentially select their highest quality, most novel or most significant research to post as preprints, however, authors retain an expectation that articles they post as preprints will receive more citations or be shared more widely online than articles not posted.

category: scientific communication and education


title: Naturalistic hyperscanning with wearable magnetoencephalography

abstract: The evolution of human cognitive function is reliant on complex social interactions which form the behavioural foundation of who we are. These social capacities are subject to dramatic change in disease and injury; yet their supporting neural substrates remain poorly understood. Hyperscanning employs functional neuroimaging to simultaneously assess brain activity in two individuals and offers the best means to understand the neural basis of social interaction. However, present technologies are limited, either by poor performance (low spatial/temporal precision) or unnatural scanning environment (claustrophobic scanners, with interactions via video). Here, we solve this problem by developing a new form of hyperscanning using wearable magnetoencephalography (MEG). This approach exploits quantum sensors for MEG signal detection, in combination with high-fidelity magnetic field control - afforded by a novel "matrix coil" system - to enable simultaneous scanning of two freely moving participants. We demonstrate our approach in a somatosensory task and an interactive ball game. Despite large and unpredictable subject motion, sensorimotor brain activity was delineated clearly in space and time, and correlation of the envelope of neuronal oscillations between people was demonstrated. In sum, unlike existing modalities, wearable-MEG combines high fidelity data acquisition and a naturalistic setting, which will facilitate a new generation of hyperscanning.


title: Impact of pH and removed filtrate on E. coli regrowth and microbial community during storage of electro-dewatered biosolids

abstract: Residual biosolids can be land applied if they meet microbiological requirements at the time of application. Electro-dewatering technology is shown to reduce biosolids bacterial counts to detection limits with little potential for bacterial regrowth during incubations. Here, we investigated the impacts on Escherichia coli regrowth and microbial communities of biosolids pH, removed nutrients via the filtrate, and inhibitory compounds produced in electro-dewatered biosolids. Findings suggest pH as the primary mechanism impacting E. coli regrowth in electro-dewatered biosolids. Propidium monoazide treatments were effective at removing DNA from dead cells, based on the removal of obligate anaerobes observed after anaerobic incubation. Analyses of high throughput sequenced data showed lower alpha-diversities associated with electro-dewatering treatment and incubation time. Moreover, biosolids pH and incubation period were the main factors contributing to the variations in microbial community compositions after incubation. Results highlight the role of electro-dewatered biosolids low pH on inhibiting the regrowth of culturable bacteria as well as reducing the microbial community variance.

category: microbiology


title: Microbial habitability of the early Mars lacustrine environment sustained by iron redox cycling

abstract:
Several studies have reported new data on the estimated compositions of chemical components at Gale crater; however, there is still a lack of information regarding potential past support of biomass and detectable biomarkers of ancient life. In this study we evaluate microbial habitability of early Mars constrained by the recently reconstructed water chemistry at Gale. The modeled community is based on Fe-metabolizing bacteria with the ability to utilize solid-phase iron oxides (e.g., magnetite) as an electron source or sink. Our results illustrate the plausibility of a sustained community in Gale Lake and provides suggestions for future modelled and laboratory-based studies to further evaluate the past habitability of Mars, biosignatures and their preservation potential, and hidden metabolic potential.

One Sentence SummaryThis work provides an existence proof of habitability on early Mars and demonstrates modeling processes by which the habitability of extraterrestrial environments can be explored quantitatively.


title: Functional characterization of 5' UTR cis-acting sequence elements that modulate translational efficiency in P. falciparum and humans

abstract:
BackgroundThe eukaryotic parasite Plasmodium falciparum causes millions of malarial infections annually while drug resistance to common antimalarials is further confounding eradication efforts. Translation is an attractive therapeutic target that will benefit from a deeper mechanistic understanding. As the rate limiting step of translation, initiation is a primary driver of translational efficiency. It is a complex process regulated by both cis and trans acting factors, providing numerous potential targets. Relative to model organisms and humans, P. falciparum mRNAs feature unusual 5 untranslated regions suggesting cis-acting sequence complexity in this parasite may act to tune levels of protein synthesis through their effects on translational efficiency.

MethodsHere, we deployed in vitro translation to compare the role of cis-acting regulatory sequences in P. falciparum and humans. Using parasite mRNAs with high or low translational efficiency, the presence, position, and termination status of upstream "AUG"s, in addition to the base composition of the 5 untranslated regions, were characterized.

ResultsThe density of upstream "AUG"s differed significantly among the most and least efficiently translated genes in P. falciparum, as did the average "GC" content of the 5 untranslated regions. Using exemplars from highly translated and poorly translated mRNAs, multiple putative upstream elements were interrogated for impact on translational efficiency. Upstream "AUG"s were found to repress translation to varying degrees, depending on their position and context, while combinations of upstream "AUG"s had nonadditive effects. The base composition of the 5 untranslated regions also impacted translation, but to a lesser degree. Surprisingly, the effects of cis-acting sequences were remarkably conserved between P. falciparum and humans.

ConclusionWhile translational regulation is inherently complex, this work contributes toward a more comprehensive understanding of parasite and human translational regulation by examining the impact of discrete cis-acting features, acting alone or in context.


title: In vitro discovery and optimization of a human monoclonal antibody that neutralizes neurotoxicity and lethality of cobra snake venom

abstract:
The monocled cobra (Naja kaouthia) is one of the most feared snakes in Southeast Asia. It is a highly dangerous species with a potent venom deriving its toxicity predominantly from abundant long-chain -neurotoxins. The only specific treatment for snakebite envenoming is antivenom, which is based on animal-derived polyclonal antibodies. Despite the lifesaving importance of these medicines over the past 120 years, and their ongoing role in combating snakebite disease, major limitations in safety, supply consistency, and efficacy creates a need for a new generation of improved treatments based on modern biotechnological techniques. Here, we describe the initial discovery and subsequent optimization of a recombinant human monoclonal immunoglobin G (IgG) antibody against -cobratoxin using phage display technology. Affinity maturation of the parental antibody by light chain-shuffling resulted in an 8-fold increase in affinity, translating to a significant increase in in vitro neutralization potency and in vivo efficacy. While the parental antibody prolonged survival of mice challenged with purified -cobratoxin, the optimized antibody prevented lethality when incubated with N. kaouthia whole venom prior to intravenous injection. This study is the first to demonstrate neutralization of whole snake venom by a single recombinant monoclonal antibody. Importantly, this suggests that for venoms whose toxicity relies on a single predominant toxin group, such as that of N. kaouthia, as little as one monoclonal antibody may be sufficient to prevent lethality, thus providing a tantalizing prospect of bringing recombinant antivenoms based on human monoclonal or oligoclonal antibodies to the clinic.

One Sentence SummaryA recombinant human monoclonal immunoglobulin G antibody, discovered and optimized using in vitro methods, was demonstrated to neutralize the lethal effect of whole venom from the monocled cobra in mice via abrogation of -neurotoxin-mediated neurotoxicity.


title: SUMOylation of Dorsal attenuates Toll/NF-κB signalling

abstract:
In Drosophila, Toll/NF-{kappa}B signalling plays key roles in both animal development and in host defence. The activation, intensity and kinetics of Toll signalling is regulated by post-translational modifications such as phosphorylation, SUMOylation or ubiquitination that target multiple proteins in the Toll/NF-{kappa}B cascade.

Here, we have generated a CRISPR-Cas9 edited Dorsal (DL) variant that is SUMO conjugation resistant (SCR). Intriguingly, embryos laid by dlSCR mothers overcome dl haploinsufficiency and complete the developmental program. This ability appears to be a result of higher transcriptional activation by DLSCR. In contrast, SUMOylation dampens DL transcriptional activation, ultimately conferring robustness to the dorso-ventral program. In the larval immune response, dlSCR animals show increase in crystal cell numbers, stronger activation of humoral defence genes, high cactus levels and cytoplasmic stabilization of DL:Cactus complexes. A mathematical model that evaluates the contribution of the small fraction of SUMOylated DL (<5%) suggests that it acts to block transcriptional activation, driven primarily by DL that is not SUMO conjugated.

Our findings define SUMO conjugation as an important regulator of the Toll signalling cascade, in both development and in host defense. Our results broadly indicate that SUMO attenuates DL at the level of transcriptional activation. Further, we hypothesize that SUMO conjugation of DL may be part of a Ubc9 dependant feedback circuit that restrains Toll/NF-{kappa}B signalling.

category: developmental biology


title: Freeze-thaw Cycles Enable a Prebiotically Plausible and Continuous Pathway from Nucleotide Activation to Nonenzymatic RNA Copying

abstract:
Nonenzymatic template-directed RNA copying using chemically activated nucleotides is thought to have played a key role in the emergence of genetic information on the early Earth. A longstanding question concerns the number and nature of different environments that might have been necessary to enable all of the steps from nucleotide synthesis to RNA replication. Here we explore three sequential steps from this overall pathway: nucleotide activation, synthesis of imidazolium-bridged dinucleotides, and template-directed primer extension. We find that all three steps can take place in one reaction mixture, under conditions of multiple freeze-thaw cycles. Recent experiments have demonstrated a potentially prebiotic methyl isocyanide-based nucleotide activation chemistry. Unfortunately, the original version of this approach is incompatible with nonenzymatic RNA copying because the high required concentration of the imidazole activating group prevents the accumulation of the essential imidazolium-bridged dinucleotide needed for primer extension. Here we report that ice eutectic phase conditions facilitate not only the methyl isocyanide-based activation of ribonucleotide 5'-monophosphates with stoichiometric 2-aminoimidazole, but also the subsequent conversion of these activated mononucleotides into imidazolium-bridged dinucleotides. Furthermore, this one pot approach is compatible with template-directed primer extension in the same reaction mixture. Our results suggest that the simple and common environmental fluctuation of freeze-thaw cycles could have played an important role in prebiotic nucleotide activation and nonenzymatic RNA copying.

Significance StatementThe replication of RNA without the aid of evolved enzymes may have enabled the inheritance of useful molecular functions during the origin of life. Several key steps on the path to RNA replication have been studied in isolation, including chemical nucleotide activation, synthesis of a key reactive intermediate, and nonenzymatic RNA copying. Here we report a prebiotically plausible scenario under which these reactions can happen together under mutually compatible conditions. Thus, this pathway could potentially have operated in nature without the complicating requirement for exchange of materials between distinct environments.

category: biochemistry


title: High-resolution structures of the actomyosin-V complex in three nucleotide states provide insights into the force generation mechanism

abstract: The molecular motor myosin undergoes a series of major structural transitions during its force-producing motor cycle. The underlying mechanism and its coupling to ATP hydrolysis and actin binding is only partially understood, mostly due to sparse structural data on actin-bound states of myosin. Here, we report 26 high-resolution cryo-EM structures of the actomyosin-V complex in the strong-ADP, rigor, and a previously unseen post-rigor transition state that binds the ATP analog AppNHp. The structures reveal a high flexibility of myosin in each state and provide valuable insights into the structural transitions of myosin-V upon ADP release and binding of AppNHp, as well as the actomyosin interface. In addition, they show how myosin is able to specifically alter the structure of F-actin. The unprecedented number of high-resolution structures of a single myosin finally enabled us to assemble a nearly complete structural model of the myosin-V motor cycle and describe the molecular principles of force production.


title: Comparative metabolomics with Metaboseek reveals functions of a conserved fat metabolism pathway in C. elegans

abstract: Untargeted metabolomics via high-resolution mass spectrometry (HRMS) can reveal more than 100,000 molecular features in a single sample, many of which may represent unidentified metabolites, posing significant challenges to data analysis. We here introduce Metaboseek, an open-source analysis platform designed for untargeted comparative metabolomics and demonstrate its utility by uncovering biosynthetic functions of a conserved fat metabolism pathway, -oxidation, using C. elegans as a model. Metaboseek integrates modules for molecular feature detection, statistics, molecular formula prediction, and fragmentation (MS/MS) analysis, which uncovered more than 200 previously uncharacterized -oxidation-dependent metabolites in an untargeted comparison of wildtype and -oxidation-defective hacl-1 mutants. The identified metabolites support the predicted enzymatic function of HACL-1 and revealed that -oxidation participates in metabolism of endogenous {beta}-methyl-branched fatty acids and food-derived cyclopropane lipids. Our results showcase compound discovery and feature annotation at scale via untargeted comparative metabolomics applied to a conserved primary metabolic pathway and suggest a model for the metabolism of cyclopropane lipids.


title: Targeting metabolic adaptations in the breast cancer liver metastatic niche using dietary approaches to improve endocrine therapy efficacy

abstract:
Estrogen receptor-positive (ER+) metastatic tumors contribute to nearly 70% of breast cancer-related deaths. Most patients with ER+ metastatic breast cancer (MBC) undergo treatment with the estrogen receptor antagonist fulvestrant (Fulv) as standard-of-care. Yet, among such patients, metastasis in liver is associated with reduced overall survival compared to other metastasis sites. The factors underlying the reduced responsiveness of liver metastases to ER-targeting agents remain unknown, impeding the development of more effective treatment approaches to improve outcomes for patients with ER+ liver metastases. We therefore evaluated site-specific changes in MBC cells and determined the mechanisms through which the liver metastatic niche specifically influences ER+ tumor metabolism and drug resistance. We characterized ER activity of MBC cells both in vitro, using a novel system of tissue-specific extracellular matrix hydrogels representing the stroma of ER+ tumor metastatic sites (liver, lung and bone), and in vivo, in liver and lung metastasis mouse models. ER+ metastatic liver tumors and MBC cells grown in liver hydrogels displayed upregulated expression of glucose metabolism enzymes in response to Fulv. Furthermore, differential ER activity, but not expression, was detected in liver hydrogels. In vivo, increased glucose metabolism led to increased glycogen deposition in liver metastatic tumors, while a fasting-mimicking diet increased efficacy of Fulv treatment to reduce the metastatic burden.

ImplicationsOur findings identify a novel mechanism of endocrine resistance driven by the liver tumor microenvironment. These results may guide the development of dietary strategies to circumvent drug resistance in liver metastasis, with potential applicability in other metastatic diseases.

category: cancer biology


title: Sapogenin based self-assembly structures activating a non-apoptotic cell death via multiple pathways

abstract: Induction of distinct cell death pathways is critical to deal with tumor heterogeneity and therapeutic resistance. In our previous study, we reported a promising saponin analog (AG-08) for cancer therapy inducing non-canonical necrotic cell death. Here, we describe that AG-08 forms unique supramolecular structures responsible for its biological activity. After internalization via non-canonical endocytosis pathway, these structures affect several cell signaling pathways including unfolded protein response, immune response, oxidative stress and heat stress. Moreover, we prepared 18 analogs to reveal the role of residues on the formation of supramolecular structures and biological activities. The results have demonstrated that unique structural features are required for particulate structures and unprecedented cell death mechanism. Although small molecule based supramolecular assemblies have widely been accepted as nuisance for drug discovery studies, our results indicate that they may provide a new research field for anti-cancer drug development studies.

category: pharmacology and toxicology


title: Dual transcriptional analysis of Peronospora belbahrii and Ocimum basilicum in susceptible interactions

abstract: Basil downy mildew, caused by the pathogen Peronospora belbahrii, is a major problem for sweet basil growers worldwide. The genome sequences of both Ocimum basilicum and P. belbahrii were recently completed but extensive transcriptome analysis of this pathosystem has not been completed. RNA sequencing was performed using basil leaf samples collected three and six days after inoculation with sporangia from an Illinois isolate of P. belbahrii and differentially expressed genes were identified. Gene enrichment analysis identified 22 genes that were upregulated at day three, in comparison to mock inoculated leaf samples, that were classified as  defense response to oomycetes; among this group were five orthologues of Arabidopsis thaliana DOWNY MILDEW RESISTANCE 6. During the same time interval, many genes contributing to photosynthesis in the infected leaves were downregulated in comparison to control leaf samples. Many more genes were differentially expressed in the inoculated basil leaves at day six, compared to mock inoculated leaves, as the pathogen began to produce sporangiophores. On days three and six, the pathogen produced high transcript levels of secreted glycoside hydrolases, which likely release sugars from the plant cell walls needed for the growth of the pathogen. These results contribute to a better understanding of the infection process of downy mildew and will aid the development of more effective measures for reducing the severity of the disease.


title: Effects of bacterial lipopolysaccharide and Shiga Toxin on induced Pluripotent Stem Cell-derived Mesenchymal Stem Cells

abstract:
BackgroundMesenchymal Stem Cells can be activated and respond to different bacterial toxins. Lipopolysaccharides (LPS) and Shiga Toxin (Stx) are the two main bacterial toxins present in Hemolytic Uremic Syndrome (HUS) that cause endothelial damage. In this work we aimed to study the response of iPSC-MSC to LPS and/or Stx and its effect on the restoration of injured endothelial cells.

MethodsiPSC-MSC were used as a source of mesenchymal stem cells (MSC) and Human Microvascular Endothelial Cells-1 (HMEC-1) as a source of endothelial cells. iPSC-MSC were treated or not with LPS and/or Stx. For some experiments, Conditioned Media (CM) were collected from each plate and incubated with an anti-Stx antibody to block the direct effect of Stx, or Polymyxin to block the direct effect of LPS. In CM from both treatments, anti-Stx and Polymyxin were used. Results are expressed as mean {+/-} S.E.M. Significant differences (p<0.05) were identified using one way analysis of variance (ANOVA) and Bonferronis Multiple comparison test.

ResultsThe results obtained showed that LPS induced a pro-inflammatory profile on iPSC-MSC, but not Stx, even though they expressed Gb3 receptor. Moreover, LPS induced on iPSC-MSC an increment in migration and adhesion to gelatin substrate. Also, the addition of CM of iPSC-MSC treated with LPS+Stx, decreased the capacity of HMEC-1 to close a wound, and did not favor the formation of new tubes. Proteomic analysis of iPSC-MSC treated with LPS and/or Stx revealed specific protein secretion patterns that support many of the functional results described here.

ConclusionsIn conclusion, these results suggest that iPSC-MSC activated by LPS acquired a pro-inflammatory profile that induces migration and adhesion to extracellular matrix proteins (ECM), but the combination LPS+Stx decreased the repair of endothelial damage. The importance of this work is that it provides knowledge to understand the context in which iPSC-MSC could benefit or not the restoration of tissue injury, taking into account that the inflammatory context in response to a particular bacterial toxin is relevant for iPSC-MSC immunomodulation.

category: immunology


title: Epigenetic reprogramming of human macrophages by an intestinal pathogen

abstract:
Various pathogens systematically reprogram gene expression in macrophages, but the underlying mechanisms are largely unknown. We investigated whether the enteropathogen Yersinia enterocolitica alters chromatin states to reprogram gene expression in primary human macrophages. Genome-wide chromatin immunoprecipitation (ChIP) seq analyses showed that pathogen-associated molecular patterns (PAMPs) induced up- or down-regulation of histone modifications (HMod) at approximately 14500 loci in promoters and enhancers. Effectors of Y. enterocolitica reorganized about half of these dynamic HMod, with the effector YopP being responsible for about half of these modulatory activities. The reorganized HMod were associated with genes involved in immune response and metabolism. Remarkably, the altered HMod also associated with 61 % of all 534 known Rho GTPase pathway genes, revealing a new level in Rho GTPase regulation and a new aspect of bacterial pathogenicity. Changes in HMod were associated to varying degrees with corresponding gene expression, e. g. depending on chromatin localization and cooperation of the HMod. In summary, infection with Y. enterocolitica remodels HMod in human macrophages to modulate key gene expression programs of the innate immune response.

Author SummaryHuman pathogenic bacteria can affect epigenetic histone modifications to modulate gene expression in host cells. However, a systems biology analysis of this bacterial virulence mechanism in immune cells has not been performed. Here we analyzed genome-wide epigenetic histone modifications and associated gene expression changes in primary human macrophages infected with enteropathogenic Yersinia enterocolitica. We demonstrate that Yersinia virulence factors extensively modulate histone modifications and associated gene expression triggered by the pathogen-associated molecular patterns (PAMPs) of the bacteria. The epigenetically modulated genes are involved in several key pathways of the macrophage immune response, including the Rho GTPase pathway, revealing a novel level of Rho GTPase regulation by a bacterial pathogen. Overall, our findings provide an in-depth view of epigenetic and gene expression changes during host-pathogen interaction and might have further implications for understanding of the innate immune memory in macrophages.


title: Yersinia remodels epigenetic histone modifications in human macrophages


title: Phosphorus is a critical factor of the in vitro monoxenic culture method for awide range of arbuscular mycorrhizal fungi culture collections

abstract: Establishing an effective way to propagate a wide range of arbuscular mycorrhizal (AM) fungi species is desirable for mycorrhizal research and agricultural applications. Although the success of mycorrhizal formation is required for spore production of AM fungi, the critical factors for its construction in the in vitro monoxenic culture protocol remain to be identified. In this study, we evaluated the growth of hairy roots from carrot, flax, and chicory, and investigated the effects of the phosphorus (P) concentration in the mother plate, as well as the levels of P, sucrose, and macronutrients in a cocultivation plate with a hairy root, amount of medium of the cocultivation plate, and location of spore inoculation, by utilizing the Bayesian information criterion model selection with greater than 800 units of data. We found that the flax hairy root was suitable for in vitro monoxenic culture, and that the concentration of P in the cocultivation plate was a critical factor for mycorrhizal formation. We showed that an extremely low concentration of P (3 M) significantly improved mycorrhizal formation for AM fungi belonging to the Glomerales order, while a high concentration of P (30 M) was suitable for Diversisporales fungi. Therefore, we anticipate that the refining the P concentration will contribute to future culture collections of a wide range of AM fungi.


title: Characterization of highly ferulate-tolerant Acinetobacter baylyi ADP1 isolates by a rapid reverse-engineering method

abstract:
Adaptive laboratory evolution (ALE) is a powerful approach for improving phenotypes of microbial hosts. Evolved strains typically contain numerous mutations that can be revealed by whole-genome sequencing. However, determining the contribution of specific mutations to new phenotypes is typically challenging and laborious. This task is complicated by factors such as the mutation type, the genomic context, and the interplay between different mutations. Here, a novel approach was developed to identify the significance of mutations in strains derived from Acinetobacter baylyi ADP1. This method, termed Rapid Advantageous Mutation ScrEening and Selection (RAMSES), was used to analyze mutants that emerged from stepwise adaptation to, and consumption of, high levels of ferulate, a common lignin-derived aromatic compound. After whole-genome sequence analysis, RAMSES allowed both rapid determination of effective mutations and seamless introduction of the beneficial mutations into the chromosomes of new strains with different genetic backgrounds. This simple approach to reverse-engineering exploits the natural competence and high recombination efficiency of ADP1. The growth advantage of transformants under selective pressure revealed key mutations in genes related to aromatic transport, including hcaE, hcaK, and vanK, and a gene, ACIAD0482, which is associated with lipopolysaccharide synthesis. This study provides insights into enhanced utilization of industrially relevant aromatic substrates and demonstrates the use of A. baylyi ADP1 as a convenient platform for strain development and evolution studies.

ImportanceMicrobial conversion of lignin-enriched streams is a promising approach for lignin valorization. However, the lignin-derived aromatic compounds are toxic to cells at relevant concentrations. Adaptive laboratory evolution is a powerful approach to develop more tolerant strains, but revealing the underlying mechanisms behind phenotypic improvement typically involves laborious processes. We employed Acinetobacter baylyi ADP1, an aromatic compound degrading strain that may be useful for biotechnology. The natural competence and high recombination efficiency of strain ADP1 can be exploited for critical applications such as the breakdown of lignin and plastics, abundant polymers composed of aromatic subunits. The natural transformability of this bacterium enabled us to develop a novel approach that allows rapid screening of advantageous mutations from ALE-derived aromatic-tolerant ADP1 strains. We clarified the mechanisms and genetic targets for improved tolerance towards common lignin-derived aromatic compounds. This study facilitates metabolic engineering for lignin valorization.


title: Analysis of Ustilago maydis Pot1 reveals context-dependent, dichotomous regulation of homology-directed DNA repair factors by a telomere protein

abstract: The telomere G-strand binding protein Pot1 plays multifaceted roles in telomere maintenance and protection. We examined the structure and activities of Pot1 in Ustilago maydis, a fungal model that recapitulates key features of mammalian telomere regulation. Compared to the well-characterized primate and fission yeast Pot1 orthologs, UmPot1 harbors an extra N-terminal OB-fold domain (OB-N), which was recently shown to be present in most metazoans. UmPot1 binds directly to Rad51 and regulates the latters strand exchange activity. Deleting the OB-N domain, which is implicated in Rad51-binding, caused telomere shortening, suggesting that Pot1-Rad51 interaction facilitates telomere replication. Depleting Pot1 through transcriptional repression triggered growth arrest as well as rampant recombination, leading to multiple telomere aberrations. In addition, telomere repeat RNAs transcribed from both the G- and C-strand were dramatically up-regulated, and this was accompanied by elevated levels of telomere RNA-DNA hybrids. Telomere abnormalities of pot1-deficient cells were suppressed, and cell viability was restored by the deletion of genes encoding Rad51 or Brh2 (the BRCA2 ortholog), indicating that homology-directed repair (HDR) proteins are key mediators of telomere aberrations and cellular toxicity. Together, these observations underscore the complex physical and functional interactions between Pot1 and DNA repair factors, leading to context-dependent and dichotomous effects of HDR proteins on telomere maintenance and protection.


title: Ustilago maydis telomere protein Pot1 harbors an extra N-terminal OB fold and regulates homology-directed DNA repair factors in a dichotomous and context-dependent manner


title: OCT4 activates a Suv39h1-repressive antisense lncRNA to couple histone H3 Lysine 9 methylation to pluripotency

abstract: Histone H3 Lysine 9 (H3K9) methylation, a characteristic mark of heterochromatin, is progressively implemented during development to contribute to cell fate restriction as differentiation proceeds. For instance, in pluripotent mouse Embryonic Stem (ES) cells the global levels of H3K9 methylation are rather low and increase only upon differentiation. Conversely, H3K9 methylation represents an epigenetic barrier for reprogramming somatic cells back to pluripotency. How global H3K9 methylation levels are coupled with the acquisition and loss of pluripotency remains largely unknown. Here, we identify SUV39H1, a major H3K9 di- and tri-methylase, as an indirect target of the pluripotency network of Transcription Factors (TFs). We find that pluripotency TFs, principally OCT4, activate the expression of an uncharacterized antisense long non-coding RNA to Suv39h1, which we name Suv39h1as. In turn, Suv39h1as downregulates Suv39h1 transcription in cis via a mechanism involving the modulation of the chromatin status of the locus. The targeted deletion of the Suv39h1as promoter region triggers increased SUV39H1 expression and H3K9me2 and H3K9me3 levels, leading to accelerated and more efficient commitment into differentiation. We report, therefore, a simple genetic circuitry coupling the global levels of H3K9 methylation to pluripotency in mouse ES cells.


title: The impact of fast radiation on the phylogeny of Bactrocera fruit flies

abstract: True fruit flies (Tephritidae) include several species that cause extensive damage to agriculture worldwide. Among them, species of the genus Bactrocera are widely studied to understand the traits associated to their invasiveness and ecology. Comparative approaches based on a reliable phylogenetic framework are particularly effective, but, to date, molecular phylogenies of Bactrocera are still controversial. Here, we employed a comprehensive genomic dataset to infer a robust backbone phylogeny of eleven representative Bactrocera species and two outgroups. We further provide the first genome scaled inference of their divergence using calibrated relaxed clock. The results of our analyses support a closer relationship of B. dorsalis to B. latifrons than to B. tryoni, in contrast to all mitochondrial-based phylogenies. By comparing different evolutionary models, we show that this incongruence likely derives from the fast and recent radiation of these species that occurred around 2 million years ago, which may be associated with incomplete lineage sorting and possibly (ongoing) hybridization. These results can serve as basis for future comparative analyses and highlight the utility of using large datasets and efficient phylogenetic approaches to study the evolutionary history of species of economic importance.

category: evolutionary biology


title: Empirical estimates of the mutation rate for an alphabaculovirus

abstract:
Mutation rates are of key importance for understanding evolutionary processes and predicting their outcomes. Empirical estimates of mutation rate are available for a number of RNA viruses, but few are available for DNA viruses, which tend to have larger genomes. Whilst some viruses have very high mutation rates, lower mutation rates are expected for viruses with large genomes to ensure genome integrity. Alphabaculoviruses are insect viruses with large genomes and often have high levels of polymorphism, suggesting high mutation rates despite evidence of proofreading activity by the replication machinery. Here, we report an empirical estimate of the mutation rate per base per strand copying (s/n/r) of Autographa californica multiple nucleopolyhedrovirus (AcMNPV). To avoid biases due to selection, we analyzed mutations that occurred in a stable, non-functional genomic insert after five serial passages in Spodoptera exigua larvae. Population bottlenecks, viral mode of replication and thresholds for mutation detection likely affect mutation rate estimates, and we therefore used population genetic models that account for these processes to infer the mutation rate. We estimated a mutation rate of 1x10-7 s/n/r. This estimate was not sensitive to different model assumptions or including whole genome data. The rates at which different classes of mutations accumulate provide good evidence for neutrality of mutations occurring within the inserted region. We therefore present a robust approach for mutation rate estimation for viruses with stable genomes, and strong evidence of a much lower alphabaculovirus mutation rate than supposed based on the high levels of polymorphism observed.

Author SummaryVirus populations can evolve rapidly, driven by the large number of mutations that occur during virus replication. It is challenging to measure mutation rates because selection will affect which mutations are observed: beneficial mutations are overrepresented in virus populations, while deleterious mutations are selected against and therefore underrepresented. Few mutation rates have been estimated for viruses with large DNA genomes, and there are no estimates for any insect virus. Here, we estimate the mutation rate for an alphabaculovirus, a virus that infects caterpillars and has a large, 134 kilobase pair DNA genome. To ensure that selection did not bias our estimate of mutation rate, we studied which mutations occurred in a large artificial region inserted into the virus genome, where mutations did not affect viral fitness. We deep sequenced evolved virus populations, and compared the distribution of observed mutants to predictions from a simulation model to estimate mutation rate. We found evidence for a relatively low mutation rate, of one mutation in every 10 million bases replicated. This estimate is in line with expectations for a virus with self-correcting replication machinery and a large genome.


title: Disruption of c-MYC binding and chromosomal looping involving genetic variants associated with ankylosing spondylitis upstream of RUNX3 promoter

abstract:
BackgroundAnkylosing Spondylitis (AS) is a common form of inflammatory spinal arthritis with a complex aetiology and high heritability, involving more than 100 genetic associations. These include several AS-associated single nucleotide polymorphisms (SNPs) upstream of RUNX3, which encodes the multifunctional RUNT-related transcription factor (TF) 3. The lead associated SNP rs6600247 (p= 2.6 x 10-15) lies ~13kb upstream of the RUNX3 promoter adjacent to a c-MYC TF binding-site. The effect of rs6600247 genotype on DNA binding and chromosome looping were investigated by electrophoretic mobility gel shift assays (EMSA), Western blotting-EMSA (WEMSA) and Chromosome Conformation Capture (3C).

ResultsInterrogation of ENCODE published data showed open chromatin in the region overlapping rs6600247 in primary human CD14+ monocytes in contrast to Jurkat T cell line or primary T-cells. The rs6600247 AS-risk allele is predicted to specifically disrupt a c-MYC binding-site. Using a 50bp DNA probe spanning rs6600247 there was consistently less binding to the AS-risk "C" allele of both purified c-MYC protein and nuclear extracts (NE) from monocyte-like U937 cells. WEMSA on U937 NE and purified c-MYC protein confirmed these differences (n=2; p<0.05). 3C experiments demonstrated negligible interaction between the region encompassing rs6600247 and the RUNX3 promoter. A stronger interaction frequency was demonstrated between the RUNX3 promoter and the previously characterised AS-associated SNP rs4648889.

ConclusionsThe lead SNP rs6600247, located in an enhancer-like region upstream of the RUNX3 promoter, modulates c-MYC binding. However, the region encompassing rs6600247 has rather limited physical interaction with the promoter of RUNX3. In contrast a clear chromatin looping event between the region encompassing rs4648889 and the RUNX3 promoter was observed. These data provide further evidence for complexity in the regulatory elements upstream of the RUNX3 promoter and the involvement of RUNX3 transcriptional regulation in AS.

category: genomics


title: Improved Docking of Protein Models by a Combination of Alphafold2 and ClusPro

abstract: It has been demonstrated earlier that the neural network based program AlphaFold2 can be used to dock proteins given the two sequences separated by a gap as the input. The protocol presented here combines AlphaFold2 with the physics based docking program ClusPro. The monomers of the model generated by AlphaFold2 are separated, re-docked using ClusPro, and the resulting 10 models are refined by AlphaFold2. Finally, the five original AlphaFold2 models are added to the 10 AlphaFold2 refined ClusPro models, and the 15 models are ranked by their predicted aligned error (PAE) values obtained by AlphaFold2. The protocol is applied to two benchmark sets of complexes, the first based on the established protein-protein docking benchmark, and the second consisting of only structures released after May 2018, the cut-off date for training AlphaFold2. It is shown that the quality of the initial AlphaFold2 models improves with each additional step of the protocol. In particular, adding the AlphaFold2 refined ClusPro models to the AlphaFold2 models increases the success rate by 23% in the top 5 predictions, whereas considering the 10 models obtained by the combined protocol increases the success rate to close to 40%. The improvement is similar for the second benchmark that includes only complexes distinct from the proteins used for training the neural network.


title: R-loops between nascent pri-miRNAs and the encoding loci promote co-transcriptional processing of miRNAs in plants

abstract: In most organisms, the maturation of nascent RNAs is coupled to transcription, undergoing many processing steps co-transcriptionally. Unlike in animals, the RNA polymerase II (RNAPII) transcribes microRNAs (miRNAs) as long and structurally variable pri-miRNAs in plants. Current evidence suggests that the miRNA biogenesis complex assembly initiates early during the transcription of pri-miRNAs in plants. However, it is unknown whether miRNA processing occurs co-transcriptionally. Here, we show that plant miRNA biogenesis is coupled to transcription in a process that relies on the formation of DNA:RNA hybrids (R-loops) between the nascent transcript and the encoding loci. We used native elongating transcript sequencing data and imaging techniques to demonstrate that plant miRNA biogenesis occurs co-transcriptionally. We found that the entire biogenesis occurs coupled to transcription for pri-miRNAs processed from the loop but requires a second nucleoplasmic step for those processed from the base of the hairpin. Furthermore, we found that co- and post-transcriptional miRNA processing mechanisms co-exist for most miRNAs in a dynamic balance. Notably, we discovered that R-loops between the 5-end single-stranded arm of the pri-miRNAs and the encoding loci anchor the transcript, promoting co-transcriptional processing. Our data demonstrate the coupling of transcription and miRNA processing in plants and discovered an unexpected function for R-loops promoting RNA processing. Furthermore, our results suggest the neo-functionalization of co-transcriptionally processed miRNAs, boosting countless regulatory scenarios.


title: Deep conservation of histone variants in Thermococcales archaea

abstract: Histones are ubiquitous in eukaryotes where they assemble into nucleosomes, binding and wrapping DNA to form chromatin. One process to modify chromatin and regulate DNA accessibility is the replacement of histones in the nucleosome with paralogous variants. Histones are also present in archaea but whether and how histone variants contribute to the generation of different physiologically relevant chromatin states in these organisms remains largely unknown. Conservation of paralogs with distinct properties can provide prima facie evidence for defined functional roles. We recently revealed deep conservation of histone paralogs with different properties in the Methanobacteriales, but little is known experimentally about these histones. In contrast, the two histones of the model archaeon Thermococcus kodakarensis, HTkA and HTkB, have been examined in some depth, both in vitro and in vivo. HTkA and HTkB exhibit distinct DNA-binding behaviours and elicit unique transcriptional responses when deleted. Here, we consider the evolution of HTkA/B and their orthologs across the order Thermococcales. We find histones with signature HTkA- and HTkB-like properties to be present in almost all Thermococcales genomes. Phylogenetic analysis indicates the presence of one HTkA- and one HTkB-like histone in the ancestor of Thermococcales and long-term maintenance of these two paralogs throughout Thermococcales diversification. Our results support the notion that archaea and eukaryotes have convergently evolved histone variants that carry out distinct adaptive functions. Intriguingly, we also detect more highly diverged histone-fold proteins, related to those found in some bacteria, in several Thermococcales genomes. The functions of these bacteria-type histones remain entirely unknown, but structural modelling suggests that they can form heterodimers with HTkA/B-like histones.


title: Symbiotic bacterial community of Drosophila melanogaster changes with nutritional modifications of the diet but can alleviate negative effects on larval phenotypes

abstract: Obesity is an increasing pandemic and is caused by multiple factors including genotype, psychological stress, and gut microbiota. Our project investigated the effects produced by high fat and high sugar dietary modifications on microbiota and metabolic phenotype of Drosophila melanogaster. Larvae raised on the high fat and high sugar diets exhibited bacterial communities that were compositionally and phylogenetically different from bacterial communities of the larvae raised on normal diets, especially if parental microbiota were removed. Several of the dominant bacteria taxa that are commonly associated with high fat and high sugar diets across model organisms and even human populations showed similar pattern in our results. Corynebacteriaceae and Erysipelotrichaceae were connected with high fat food, while Enterobacteriaceae and Lactobacillaceae were associated with high sugar diets. In addition, we observed that presence of symbiotic microbiota often mitigated the effect that harmful dietary modifications produced on larvae, including elevated triglyceride concentrations and was crucial for Drosophila survival, especially on high sugar peach diets.


title: Non-canonical Staphylococcus aureus pathogenicity island repression

abstract:
Mobile genetic elements (MGEs) control their life cycles by the expression of a master repressor, whose function must be disabled to allow the spread of these elements in nature. Here we describe an unprecedented repression-derepression mechanism involved in the transfer of the Staphylococcus aureus pathogenicity islands (SaPIs). Contrary to the classical phage and SaPI repressors, which are dimers, the SaPI1 repressor StlSaPI1 presents a unique tetrameric conformation, never seen before. Importantly, not just one but two tetramers are required for SaPI1 repression, which increases the novelty of the system. To derepress SaPI1, the phage-encoded protein Sri binds to and induces a conformational change in the DNA binding domains of StlSaPI1, preventing the binding of the repressor to its cognate StlSaPI1 sites. Finally, our findings demonstrate that this system is not exclusive to SaPI1 but widespread in nature. Overall, our results characterise a novel repression-induction system involved in the transfer of MGE-encoded virulence factors in nature.

SignificanceWhile most repressors controlling the transfer of mobile genetic elements are dimers, we demonstrate here that the Staphylococcal pathogenicity island 1 (SaPI1) is repressed by two tetramers, which have a novel structural fold in their body that has never been seen before in other proteins. Moreover, by solving the structure of the SaPI1 repressor in complex with its inducing protein Sri, we have demonstrated that Sri forces the SaPI1 repressor to adopt a conformation that is incompatible with DNA binding, explaining how SaPI1 is induced. Finally, our results demonstrate that this repression system is not exclusive of the SaPIs but widespread in nature. Our studies provide important insights understanding how SaPIs spread in nature.


title: Arg-tRNA synthetase links inflammatory metabolism to RNA splicing via nuclear condensates

abstract: Cells respond to perturbations like inflammation by sensing changes in metabolite levels. Especially prominent is arginine, which has known connections to the inflammatory response. Here, we found that depletion of arginine during inflammation decreased levels of a nuclear form of arginyl-tRNA synthetase (ArgRS). Surprisingly, we found that nuclear ArgRS interacts with serine/arginine repetitive matrix protein 2 (SRRM2), a spliceosomal protein and nuclear speckle component and that arginine depletion impacted both condensate-like nuclear trafficking of SRRM2 and splice-site usage in certain genes. These splice-site usage changes cumulated in synthesis of different protein isoforms that altered cellular metabolism and peptide presentation to immune cells. Our findings uncover a novel mechanism whereby a tRNA synthetase cognate to a key amino acid that is metabolically controlled during inflammation modulates the splicing machinery.


title: Arg-tRNA synthetase links inflammatory metabolism to RNA splicing and nuclear trafficking via SRRM2


title: Human generation times across the past 250,000 years

abstract: The generation times of our recent ancestors can tell us about both the biology and social organization of prehistoric humans, placing human evolution on an absolute timescale. We present a method for predicting historic male and female generation times based on changes in the mutation spectrum. Our analyses of whole-genome data reveal an average generation time of 26.9 years across the past 250,000 years, with fathers consistently older (30.7 years) than mothers (23.2 years). Shifts in sex-averaged generation times have been driven primarily by changes to the age of paternity rather than maternity, though we report a disproportionate increase in female generation times over the past several thousand years. We also find a large difference in generation times among populations, with samples from current African populations showing longer ancestral generation times than non-Africans for over a hundred thousand years, reaching back to a time when all humans occupied Africa.


title: Multiscale Simulation of an Influenza A M2 Channel Mutant Reveals Key Features of Its Markedly Different Proton Transport Behavior

abstract: The influenza A M2 channel, a prototype for the viroporin class of viral channels, is an acid-activated viroporin that conducts protons across the viral membrane, a critical step in the viral life cycle. As the protons enter from the viral exterior, four central His37 residues control the channel activation by binding subsequent protons, which opens the Trp41 gate and allows proton flux to the viral interior. Asp44 is essential for maintaining the Trp41 gate in a closed state at high pH, which results in asymmetric conduction. The prevalent D44N mutant disrupts this gate and opens the C-terminal end of the channel, resulting in overall increased conduction in the physiologically relevant pH range and a loss of this asymmetric conduction. Here, we use extensive Multiscale Reactive Molecular Dynamics (MS-RMD) and Quantum Mechanics/Molecular mechanics (QM/MM) simulations with an explicit, reactive excess proton to calculate the free energy of proton transport in the M2 mutant and to study the dynamic molecular-level behavior of D44N M2. We find that this mutation significantly lowers the barrier of His37 deprotonation in the activated state and shifts the barrier for entry up to the Val27 tetrad. These free energy changes are reflected in structural shifts. Additionally, we show that the increased hydration around the His37 tetrad diminishes the effect of the His37 charge on the channels water structure, facilitating proton transport and enabling activation from the viral interior. Altogether, this work provides key insight into the fundamental characteristics of PT in WT M2 and how the D44N mutation alters this PT mechanism, and it expands our understanding of the role of emergent mutations in viroporins.


title: Differentially optimized cell-free buffer enables robust expression from unprotected linear DNA in exonuclease-deficient extracts

abstract: The use of linear DNA templates in cell-free systems promises to accelerate the prototyping and engineering of synthetic gene circuits. A key challenge is that linear templates are rapidly degraded by exonucleases present in cell extracts. Current approaches tackle the problem by adding exonuclease inhibitors and DNA-binding proteins to protect the linear DNA, requiring additional time- and resource-intensive steps. Here, we delete the recBCD exonuclease gene cluster from the Escherichia coli BL21 genome. We show that the resulting cell-free systems, with buffers optimized specifically for linear DNA, enable near-plasmid levels of expression from {sigma}70 promoters in linear DNA templates without employing additional protection strategies. When using linear or plasmid DNA templates at the buffer calibration step, the optimal potassium glutamate concentrations obtained when using linear DNA were consistently lower than those obtained when using plasmid DNA for the same extract. We demonstrate the robustness of the exonuclease deficient extracts across seven different batches and a wide range of experimental conditions across two different laboratories. Finally, we illustrate the use of the {Delta}recBCD extracts for two applications: toehold switch characterization and enzyme screening. Our work provides a simple, efficient, and cost-effective solution for using linear DNA templates in cell-free systems and highlights the importance of specifically tailoring buffer composition for the final experimental setup. Our data also suggest that similar exonuclease deletion strategies can be applied to other species suitable for cell-free synthetic biology.

category: synthetic biology


title: The response of three-dimensional pancreatic alpha and beta cell co-cultures to oxidative stress

abstract: The pancreatic islets of Langerhans have low endogenous antioxidant levels and are thus especially sensitive to oxidative stress, which is known to influence cell survival and behaviour. As bioengineered islets are gaining interest for therapeutic purposes, it is important to understand how their composition can be optimized to diminish oxidative stress. We investigated how the ratio of the two main islet cell types (alpha and beta cells) and their culture in three-dimensional aggregates could protect against oxidative stress. Monolayer and aggregate cultures were established by seeding the alphaTC1 (alpha) and INS1E (beta) cell lines in varying ratios, and hydrogen peroxide was applied to induce oxidative stress. Viability, oxidative stress, and the level of the antioxidant glutathione were measured. Both aggregation and an increasing prevalence of INS1E cells in the co-cultures conferred greater resistance to cell death induced by oxidative stress. Increasing the prevalence of INS1E cells also decreased the number of alphaTC1 cells experiencing oxidative stress in the monolayer culture. In 3D aggregates, culturing the alphaTC1 and INS1E cells in a ratio of 50:50 prevented oxidative stress in both cell types. Together, the results of this study lead to new insight into how modulating the composition and dimensionality of a co-culture can influence the oxidative stress levels experienced by the cells.

category: cell biology


title: High temporal resolution measurements of movement reveal novel early-life physiological decline in C. elegans.

abstract: Age-related physiological changes are most notable and best-studied late in life, while the nature of aging in early- or middle-aged individuals has not been explored as thoroughly. In C. elegans, studies of movement vs. age generally delineate three distinct phases: sustained, youthful movement; a discrete onset of rapidly progressing impairment; and gross immobility. We investigated whether this first period of early-life adult movement is simply a sustained "healthy" level of high function followed by a discrete "movement catastrophe" -- or whether there are early-life changes in movement that precede future physiological declines. To determine how movement varies during early adult life, we followed isolated individuals throughout life with a previously unachieved combination of duration and temporal resolution. By tracking individuals across the first six days of adulthood, we observed declines in movement starting as early as the first two days of adult life, as well as high interindividual variability in total daily movement. These findings suggest that movement is a highly dynamic behavior early in life, and that factors driving movement decline may begin acting as early as the first day of adulthood. Using simulation studies based on acquired data, we suggest that too infrequent sampling in common movement assays limits observation of early-adult changes in motility, and we propose feasible alternate strategies and a framework for designing assays with increased sensitivity for early movement declines.

category: animal behavior and cognition


title: Antidepressant Sertraline Hydrochloride Inhibits the Growth of HER2+ AU565 Breast Cancer Cell Line through Induction of Apoptosis, and Arrest of Cell Cycle

abstract: Breast cancer is one the most aggressive cancer worldwide, especially Pakistan due to limited therapeutic options. This study was conducted to repurpose the use of selective serotonin reuptake inhibitors (SSRIs), in the treatment of breast cancers, and merit to pursue drug re-positioning in oncology. Anti-proliferative activity of SSRIs, such as fluoxetine, paroxetine, and sertraline hydrochloride on the growth of AU-565, MCF-7, MDA-MB-231, and BT-474 breast cancer cell lines, along with human fibroblast BJ cells was determined in vitro. Changes in nuclear morphology (DAPI staining), and induction of apoptosis (flow cytometry, and caspase-3 activation) were also studied. Sertraline hydrochloride most effectively inhibited the growth of breast cancer cells in vitro. Therefore, pharmacological mechanism involved in sertraline mediated cell death was investigated in HER2+ AU565 cell line. Enhanced nuclear fragmentation, increased Annexin (+) cells, and caspase-3/7 activation indicated that sertraline-mediated cell death could be a result of BCl2-independent apoptosis as evidenced by expression of Bax, and BCl2 genes. Taken together, our results identified sertraline hydrochloride, as a potential candidate for the treatment of HER2-positive breast cancer. Even though these are in vitro results, this study opens great opportunity in the field of drug repurposing for the development of chemotherapeutic agents.


title: Declining RNA integrity in control autopsy brain tissue is robustly and asymmetrically associated with selective neuronal mRNA signal loss

abstract:
RNA integrity numbers (RINs) are a standardized method for semi-quantification of RNA degradation, and are used in quality control prior to transcriptional profiling analysis. Recent work has demonstrated that RINs are associated with downstream transcriptional profiling, and correction procedures are typically employed in bioinformatic analysis pipelines to attempt to control for RINs influence on gene expression. However, relatively little work has been done to determine whether RINs influence is random, or is specifically targeted to a subset of mRNAs. We tested the hypothesis that RIN would be associated with a robust transcriptional profile seen across multiple studies.

To test this, we downloaded subsets of raw transcriptional data from six published studies. We only included control, non-pathological post-mortem human brain tissue (n = 383 samples) in which independent subjects RIN values were also reported. A robust set of mRNAs consistently and significantly correlated with RIN across multiple studies, appearing to be selectively degraded as RIN declines. Many of the affected gene expression pathways are related to neurons (e.g., vesicle, mRNA transport, synapse, and mitochondria), suggesting that neuronal synaptic mRNA may be particularly vulnerable to degradation. Subsequent analysis of the relationship between RIN and vulnerable mRNA expression revealed most of the decay occurred over a relatively narrow RIN range of 7.2-8.6, with RIN values > 8.6 showing a ceiling effect, and those < 7.2 showing a floor effect on gene expression. Our data suggests that the RIN effect is pathway selective and non-linear, which may be an important consideration for current bioinformatic RIN correcting procedures, particularly in datasets in which declining RIN is confounded with a pathology under study (e.g., in Alzheimers disease).

category: bioinformatics


title: Neurophysiological evidence for cognitive map formation during sequence learning

abstract: Humans deftly parse statistics from sequences. Some theories posit that humans learn these statistics by forming cognitive maps, or underlying representations of the latent space which links items in the sequence. Here, an item in the sequence is a node, and the probability of transitioning between two items is an edge. Sequences can then be generated from walks through the latent space, with different spaces giving rise to different sequence statistics. Individual or group differences in sequence learning can be modeled by changing the time scale over which estimates of transition probabilities are built, or in other words, by changing the amount of temporal discounting. Latent space models with temporal discounting bear a resemblance to models of navigation through Euclidean spaces. However, few explicit links have been made between predictions from Euclidean spatial navigation and neural activity during human sequence learning. Here, we use a combination of behavioral modeling and intracranial encephalography (iEEG) recordings to investigate how neural activity might support the formation of space-like cognitive maps through temporal discounting during sequence learning. Specifically, we acquire human reaction times from a sequential reaction time task, to which we fit a model that formulates the amount of temporal discounting as a single free parameter. From the parameter, we calculate each individuals estimate of the latent space. We find that neural activity reflects these estimates mostly in the temporal lobe, including areas involved in spatial navigation. Similar to spatial navigation, we find that low dimensional representations of neural activity allow for easy separation of important features, such as modules, in the latent space. Lastly, we take advantage of the high temporal resolution of iEEG data to determine the time scale on which latent spaces are learned. We find that learning typically happens within the first 500 trials, and is modulated by the underlying latent space and the amount of temporal discounting characteristic of each participant. Ultimately, this work provides important links between behavioral models of sequence learning and neural activity during the same behavior, and contextualizes these results within a broader framework of domain general cognitive maps.


title: External drivers of BOLD signals non-stationarity

abstract: A fundamental challenge in neuroscience is to uncover the principles governing how the brain interacts with the external environment. However, assumptions about external stimuli fundamentally constrain current computational models. We show in silico that unknown external stimulation can produce error in the estimated linear time-invariant dynamical system. To address these limitations, we propose an approach to retrieve the external (unknown) input parameters and demonstrate that the estimated system parameters during external input quiescence uncover spatiotemporal profiles of external inputs over external stimulation periods more accurately. Finally, we unveil the expected (and unexpected) sensory and task-related extra-cortical input profiles using functional magnetic resonance imaging data acquired from 96 subjects (Human Connectome Project) during the resting-state and task scans. Together, we provide evidence that this embodied brain activity model offers information about the structure and dimensionality of the BOLD signals external drivers and shines light on likely external sources contributing to the BOLD signals non-stationarity.


title: Different encoding of reward location in dorsal and ventral hippocampus

abstract: Hippocampal neurons encode a cognitive map for spatial navigation1. When they fire at specific locations in the environment, they are known as place cells2. In the dorsal hippocampus place cells accumulate at current navigational goals, such as learned reward locations3-6. In the intermediate-to-ventral hippocampus (here collectively referred to as ventral hippocampus), neurons fire across larger place fields7-10 and regulate reward- seeking behavior11-16, but little is known about their involvement in reward-directed navigation. Here, we compared the encoding of learned reward locations in the dorsal and ventral hippocampus during spatial navigation. We used calcium imaging with a head- mounted microscope to track the activity of CA1 cells over multiple days during which mice learned different reward locations. In dorsal CA1 (dCA1), the overall number of active place cells increased in anticipation of reward but the recruited cells changed with the reward location. In ventral CA1 (vCA1), the activity of the same cells anticipated the reward locations. Our results support a model in which the dCA1 cognitive map incorporates a changing population of cells to encode reward proximity through increased population activity, while the vCA1 provides a reward-predictive code in the activity of a specific subpopulation of cells. Both of these location-invariant codes persisted over time, and together they provide a dual hippocampal reward-location code, assisting goal- directed navigation17, 18.


title: A generalized cortical activity pattern at internally-generated mental context boundaries during unguided narrative recall

abstract: Current theory and empirical studies suggest that humans segment continuous experiences into events based on the mismatch between predicted and actual sensory inputs; detection of these "event boundaries" evokes transient neural responses. However, boundaries can also occur at transitions between internal mental states, without relevant external input changes. To what extent do such "internal boundaries" share neural response properties with externally-driven boundaries? We conducted an fMRI experiment where subjects watched a series of short movies and then verbally recalled the movies, unprompted, in the order of their choosing. During recall, transitions between movies thus constituted major boundaries between internal mental contexts, generated purely by subjects unguided thoughts. Following the offset of each recalled movie, we observed stereotyped spatial activation patterns in the default mode network, especially the posterior medial cortex, consistent across different movie contents and even across the different tasks of movie watching and recall. Surprisingly, the between-movie boundary patterns were negatively correlated with patterns at boundaries between events within a movie. Thus, major transitions between mental contexts elicit neural phenomena shared across internal and external modes and distinct from within-context event boundary detection, potentially reflecting a cognitive state related to the flushing and reconfiguration of situation models.


title: Retrospective analysis of enhancer activity and transcriptome history

abstract: Cell state changes in development and disease are controlled by gene regulatory networks, the dynamics of which are difficult to track in real time. Here, we utilize an inducible DCM-RNA-polymerase-subunit-b fusion protein, to label active genes and enhancers with a bacterial methylation mark that does not affect gene transcription and is propagated in S-phase. We applied this DCM-time-machine (DCM-TM) technology to study intestinal homeostasis, following enterocyte differentiation back in time, revealing rapid and simultaneous activation of enhancers and nearby genes during intestinal stem cell (ISC) differentiation. We provide new insights in the absorptive-secretory lineage decision in ISC differentiation, and show that ISCs retain a unique chromatin landscape required to maintain ISC identity and delineate future expression of differentiation associated genes. DCM-TM has wide applicability in tracking cell states, providing new insights in the regulatory networks underlying cell state changes in development and differentiation.


title: Illuminating links between cis-regulators and trans-acting variants in the human prefrontal cortex

abstract: Psychiatric disorders exact immense human and economic tolls in societies globally. Underlying many of these disorders is a complex repertoire of genomic variants that influence the expression of genes involved in pathways and processes in the brain. Identifying such variants and their associated brain functions is thus essential for understanding the molecular underpinnings of psychiatric disorders. Genome-wide association studies (GWASes) have provided many variants associated with these disorders; however, our knowledge of the precise biological mechanisms by which these contribute to disease remains limited. In connection with this, expression quantitative trait loci (eQTLs) have provided useful information linking variants to genes and functions. However, most eQTL studies on human brain have focused exclusively on cis-eQTLs. A complete understanding of disease etiology should also include trans-regulatory mechanisms. Thus, we conduct one of the first genome-wide surveys of trans-eQTLs in the dorsolateral prefrontal cortex (DLPFC) by leveraging the large datasets from the PsychENCODE consortium. We identified [~]80,000 trans-eQTLs. We found that a significant number of these overlap with cis-eQTLs, thereby implicating cis-mediators as key players in trans-acting regulation. We show, furthermore, that trans-regulatory mechanisms provide novel insights into psychiatric disease. Particularly, colocalization analysis between trans-eQTLs and schizophrenia (SCZ) GWAS loci identified 90 novel SCZ risk genes and 23 GWAS loci previously uncharacterized by cis-eQTLs. Moreover, these 90 genes tend to be more central in transcriptome-wide co-expression networks and more susceptible to rare variants than SCZ-risk genes associated by cis-variation.


title: Mutation, selection, and quantitative genetic architecture of susceptibility to bacterial pathogens in C. elegans

abstract: Understanding the evolutionary and genetic underpinnings of susceptibility to pathogens is of fundamental importance across a wide swathe of biology. Much theoretical and empirical effort has focused on genetic variants of large effect, but pathogen susceptibility often appears to be a polygenic complex trait. Here we investigate the quantitative genetics of survival over 120 hours of exposure ("susceptibility") of C. elegans to three bacterial pathogens of varying virulence, along with the non-pathogenic OP50 strain of E. coli. We compare the genetic (co)variance input by spontaneous mutations accumulated under minimal selection to the standing genetic (co)variance in a set of [~]50 wild isolates. Three conclusions emerge. First, with one exception, mutations increase susceptibility to pathogens, and susceptibility is uncorrelated with fitness in the absence of pathogens. Second, the orientation in trait space of the heritable (co)variance of wild isolates is sufficiently explained by mutation. However, pathogen susceptibility is clearly under purifying, apparently directional, selection of magnitude similar to that of competitive fitness in the MA conditions. The results provide no evidence for fitness tradeoffs between pathogen susceptibility and fitness in the absence of pathogens, nor that balancing selection is important in maintaining genetic variation for susceptibility to these bacterial pathogens.


title: Roles of cellular senescence in driving bone marrow adiposity in radiation-and aging-associated bone loss

abstract: Osteoporosis is associated with an increase in marrow adipocytes, collectively termed bone marrow adipose tissue (BMAT). An increase in BMAT is linked with decline of mesenchymal progenitors that give rise to osteoblasts which are responsible for bone accrual. Oxidative stress-induced reactive oxygen species, DNA damage, apoptosis and cellular senescence have been associated with reduced osteoprogenitors in a reciprocal fashion to BMAT; however, a direct (causal) link between cellular senescence and BMAT is still elusive. Accumulation of senescent cells occur in naturally aged and in focally radiated bone tissue, but despite amelioration of age- and radiation-associated bone loss after senescent cell clearance, molecular events that precede BMAT accrual are largely unknown. Using a mouse model here we show by RNA-Sequencing data that BMAT-related genes were the most upregulated gene subset in radiated bones. Using focal radiation as a model to understand age-associated changes in bone, we performed a longitudinal assessment of cellular senescence and BMAT. Using qRT-PCR, RNA in situ hybridization and histological assessment of telomere dysfunction as a marker of senescence, we observed increased p21 transcripts in bone lining cells, osteocytes and bone marrow cells, and elevated dysfunctional telomeres in osteocytes starting from day 1 post-radiation, without the presence of BMAT. BMAT was significantly elevated in radiated bones at day 7, confirming the qRT-PCR data in which most BMAT-related genes were elevated by day 7, and the trend continued until day 42 post-radiation. Similarly, elevation in BMAT-related genes was observed in aged bones. The senolytic cocktail of Dasatinib (D) plus Quercetin (Q) - D+Q, which clears senescent cells, reduced BMAT in aged and radiated bones. MicroRNAs (miRs) linked with senescence marker p21 were downregulated in radiated- and aged-bones, while miR-27a, a miR that is associated with increased BMAT, was elevated both in radiated- and aged-bones. D+Q downregulated miR-27a in radiated bones at 42 days post-radiation. Overall, our study provides evidence that BMAT occurrence in oxidatively stressed bone environments, such as radiation and aging, is induced following a common pathway and is dependent on the presence of senescent cells.


title: Protein reservoirs of seeds are composites of amyloid and amyloid-like structures

abstract: The functional roles of amyloids are known across diverse species, but their natural presence in plants remains underexplored. Based on the structure and function of protein bodies found in other organisms, we hypothesize that the seed storage protein bodies (SSPB) may have amyloid structures for controlling physiological functions. Here we show that monocot and dicot SSPB exhibit a speckled-pattern of amyloids interspersed in an amyloid-like matrix, suggesting their composite nature. This is confirmed by amyloid-specific probes, amyloidogenic region scans, infrared signatures, reconstitution studies, and differential degradation during germination. The role of amyloidic composites in seed germination is proved by amyloid content, degradation inhibition and relevant germination parameters, using in-vitro SSPB, ex-vivo protoplasts and in-situ seed sections. The results open new research questions in understanding molecular details of SSPB assembly-disassembly and utilization, in addition to evolutionary roles of amyloid structure and function.


title: Protein reservoirs of seeds are composites of amyloid and amyloid-like structures facilitating germination


title: Crystal structure of a constitutive active mutant of adenosine A2A receptor

abstract: Previously we reported a common activation pathway of the class A G protein-coupled receptors (GPCRs) in which a series of conserved residues/motifs undergo conformational change during extracellular agonist binding and finally induce the coupling of intracellular G protein, and successfully predicted several novel constitutive active or inactive mutations for A2A adenosine receptor (A2AAR) through this mechanism (Zhou et al., 2019). Here we determined the crystal structure of a typical A2AAR constitutive active mutant I92N in complex with agonist UK-432097 to reveal the molecular mechanism of mutation-induced constitutive activity. The mutated I92N forms a hydrophilic interaction network with nearby residues including W6.48 of the CWxP motif, which is absent in the wild type A2AAR. Although the mutant structure is overall similar to the previously determined intermediate state A2AAR structure (PDB ID: 3QAK), the molecular dynamics simulations suggested that the I92N mutant stabilizes the metastable intermediate state through the hydrophilic interaction network and favors the receptor conformational transition towards the active state. This research provides a structural template toward the special pharmacological outcome triggered by conformational mutation and sheds light on future structural or pharmacological studies among class A GPCRs.