title: Identification and quantification of chimeric sequencing reads in a highly multiplexed RAD-seq protocol

abstract: Highly multiplexed approaches have become a common practice in genomic studies. They have improved the cost-effectiveness of genotyping hundreds of individuals by using combinatorially-barcoded adapters. These strategies, however, can potentially misassign reads to incorrect samples. Here we used a modified quaddRAD protocol to analyse the occurrence of index hopping and PCR chimeras in a series of experiments with up to a 100 multiplexed samples per sequencing lane (total n = 639). We created two types of sequencing libraries: four libraries of Type A, where PCR reactions were run on individual samples before multiplexing, and three libraries of Type B, where PCRs were run on pooled samples. We used fixed pairs of inner barcodes to identify chimeric reads. Type B libraries show a higher percentage of misassigned reads (1.15%) compared to Type A libraries (0.65%). We also quantify the commonly undetectable chimeric sequences that occur whenever multiplexed groups of samples with different outer barcodes are sequenced together on a single flow cell. Our results suggest that these types of chimeric sequences represent up to 1.56% and 1.29% of reads in Type A and B libraries, respectively. We review the source of such errors, provide recommendations for developing highly-multiplexed RAD-seq protocols and analysing the resulting data to minimise the generation of chimeric sequences, allow their quantification, and provide finer control over the number of PCR cycles necessary to generate enough input DNA for library preparation.

category: bioinformatics


title: GBRAP: a tool to retrieve, parse and analyze GenBank files of viral and bacterial species

abstract:
SummaryGenBank files contain genomic data of sequenced living organisms. Here, we present GBRAP (GenBank Retrieving, Analyzing and Parsing software), a tool written in Python 3 that can be used to easily download, parse and analyze viral and bacterial GenBank files, even when contain more than one genomic sequence for each species. GBRAP can analyze more files simultaneously through single command-line parameters that give as output a single table showing the genomic characteristics of each organism. It is also able to calculate Shannon, LZSS (Lempel-Ziv-Storer-Szymanski) and topological entropy for both the entire genome and its constitutive elements such as genes, rRNAs, tRNAs, tmRNAs and ncRNAs together with Chargaffs second parity rule scores obtained using different mathematical methods. Moreover, GBRAP can calculate, the number, the length and the nucleotides abundance of genomic components for each DNA strand and for the overlapping regions among the two complementary helixes. To our knowledge, this is the only software capable of providing this type of genomic analyses all together in a single tool, that, therefore can be used by the scientists interested in both genomics and evolutionary research.

Availability and implementationThe data underlying this article are available from the corresponding author on reasonable request.


title: piRNAdb: A piwi-interacting RNA database

abstract:
MotivationFound in several metazoan species, piwi-interacting RNAs (piRNAs) regulate the expression of a wide variety of transposable elements and genes associated with cellular development, differentiation, and growth. Despite their importance, piRNAs are not well known and are still underexplored. To facilitate piRNA research, a comprehensive and easy-to-use piRNA database is still needed.

ResultsHere, we present piRNAdb, an integrative and user-friendly database designed to encompass several aspects of piRNAs. We selected piRNAs from four reliable human RNA-sequencing datasets to start our database. After data processing, we displayed these sequences, their genomic location, clustering information, putative targets on known genes and transposable elements, as well as direct links to other databases. In this first release, piRNAdb catalogues 27,329 piRNAs, as well as 23,380 genes that are putative targets and 47,060 associated gene ontology terms, both of which are organized and linked to each respective piRNA. Finally, to improve information exchange and increase the confidence of sequences, a feedback system is provided to users of piRNAdb.

ConclusionThe inclusion of new features to facilitate piRNA analyses, data visualization, and integration is the major pillar of piRNAdb. Our main goal was to make this database an easy interface between the data and the user. We believe that this web tool achieves this objective by providing a streamlined and well-organized data repository for piRNAs and that it will be extremely useful to those already studying piRNAs and to the broader community.

AvailabilitypiRNAdb is available freely and is compatible with smartphones and tablets: https://www.pirnadb.org/


title: The genes controlling normal function of citrate and spermine secretion is lost in aggressive prostate cancer and prostate model systems

abstract:
BackgroundSecretion of the metabolites citrate and spermine into prostate lumen is a unique hallmark for normal prostate epithelial cells. However, the identity of the genes controlling citrate and spermine secretion remains mostly unknown despite their obvious relevance for progression to aggressive prostate cancer.

Materials & MethodsIn this study, we have correlated simultaneous measurement of citrate/spermine and transcriptomics data. We have refined these gene correlations in 12 prostate cancer cohorts containing 2915 tissue samples to create a novel gene signature of 150 genes connected with citrate and spermine secretion. We further explored the signature in public data, interrogating over 18 000 samples from various tissues and model systems, including 3826 samples from prostate and prostate cancer.

ResultsIn prostate cancer, the expression of this gene signature is gradually lost in tissue from normal epithelial cells through PIN, low grade (Gleason <= 7), high grade cancer (Gleason >= 8) and metastatic lesions. The accuracy of the signature is validated by its unique enrichment in prostate compared to other tissues, and its strong enrichment in epithelial tissue compartments compared to stroma. Several zinc-binding proteins that are not previously investigated in the prostate are present in the gene signature, suggesting new mechanisms for controlling zinc homeostasis in citrate/spermine secretion. However, the absence of the gene signature in all common prostate normal and cancer cell-lines, as well as prostate organoids, underlines the challenge to study the role of these genes during prostate cancer progression in model systems.

ConclusionsA large collection of transcriptomics data integrated with metabolomics identifies the genes related to citrate and spermine secretion in the prostate, and show that the expression of these genes gradually decreases on the path towards aggressive prostate cancer. In addition, the study questions the relevance of currently available model systems to study metabolism in prostate cancer development.


title: A waveform-independent measure of recurrent neural activity

abstract: Rhythmic neural activity, so called oscillations, play a key role for neural information transmission, processing and storage. Neural oscillations in distinct frequency bands are central to physiological brain function and alterations thereof have been associated with several neurological and psychiatric disorders. The most common methods to analyse neural oscillations, e.g. short-term Fourier transform or wavelet analysis, assume that measured neural activity is composed of a series of symmetric prototypical waveforms, e.g. sinusoids. However, usually the models generating the signal, including waveform shapes of experimentally measured neural activity are unknown. Decomposing asymmetric waveforms of nonlinear origin using these classic methods may result in spurious harmonics visible in the estimated frequency spectra. Here, we introduce a new method for capturing rhythmic brain activity based on recurrences of similar states in phase-space. This method allows for a time-resolved estimation of amplitude fluctuations of recurrent activity irrespective of or specific to waveform-shapes. The algorithm is derived from the well-established field of recurrence analysis, which has rarely been adopted in neuroscience. In this paper, we show its advantages and limitations in comparison to short-time Fourier transform and wavelet convolution using periodic signals of different waveform shapes. Further, we demonstrate its application using experimental data, i.e. intracranial electrophysiological recordings from the human motor cortex of one epilepsy patient.


title: Objective method reveals minimum time required to quantify net cost of transport across slow to fast walking speeds

abstract:
The U-shaped net cost of transport (COT) curve of walking has helped scientists understand the biomechanical basis that underlies energy minimization during walking. However, to produce an individuals net COT curve, data must be analyzed during periods of steady-rate metabolism. Traditionally, studies analyze the last few minutes of a 6-10 min trial, assuming that steady-rate metabolism has been achieved. Yet, it is possible that an individual achieves steady rates of metabolism much earlier. However, there is no consensus on how to objectively quantify steady-rate metabolism across a range of walking speeds. Therefore, we developed an objective method to determine the minimum time needed for humans to achieve steady rates of metabolism across slow to fast walking speeds. We hypothesized that a shorter time window could be used to produce a net COT curve that is comparable to the net COT curve created using traditional methods. We analyzed metabolic data from twenty-one subjects who completed several 7-min walking trials ranging from 0.50-2.00 m/s. We partitioned the metabolic data for each trial into moving 1-min, 2-min, and 3 min intervals and calculated their slopes. We statistically compared these slope values to values derived from the last 3-min of the 7-min trial, our  gold standard comparison. We found that a minimum of 2 min is required to achieve steady-rate metabolism and that data from 2-4 min yields a net COT curve that is not statistically different from the one derived from experimental protocols that are generally accepted in the field.

New and NoteworthyA simple slope method reveals that 4-min of walking at speeds between 0.50-1.75 m/s is the minimum time required to produce an accurate net cost of transport curve.

category: physiology


title: Transiently increased intercommunity regulation characterizes concerted cell phenotypic transition

abstract: Phenotype transition takes place in many biological processes such as differentiation and reprogramming. A fundamental question is how cells coordinate switching of expressions of clusters of genes. Through analyzing single cell RNA sequencing data in the framework of transition path theory, we studied how such a genome-wide expression program switching proceeds in three different cell transition processes. For each process we reconstructed a reaction coordinate describing the transition progression, and inferred the gene regulation network (GRN) along the reaction coordinate. In all three processes we observed common pattern that the effective number and strength of regulation between different communities increase first and then decrease. The change accompanies with similar change of the GRN frustration, defined as overall confliction between the regulation received by genes and their expression states, and GRN heterogeneity. While studies suggest that biological networks are modularized to contain perturbation effects locally, our analyses reveal a general principle that during a cell phenotypic transition intercommunity interactions increase to concertedly coordinate global gene expression reprogramming, and canalize to specific cell phenotype as Waddington visioned.

category: systems biology


title: The HDAC inhibitor CI-994 acts as a molecular memory aid by facilitating synaptic and intra-cellular communication after learning

abstract:
Long-term memory formation relies on synaptic plasticity, activity-dependent transcription and epigenetic modifications. Multiple studies have shown that HDAC inhibitor (HDACi) treatments can enhance individual aspects of these processes, and thereby act as putative cognitive enhancers. However, their mode of action is not fully understood. In particular, it is unclear how systemic application of HDACis, which are devoid of substrate specificity, can target pathways that promote memory formation. In this study, we explore the electrophysiological, transcriptional and epigenetic responses that are induced by CI-994, a class I HDAC inhibitor, combined with contextual fear conditioning (CFC) in mice. We show that CI-994-mediated improvement of memory formation is accompanied by enhanced long-term potentiation in the hippocampus, a brain region recruited by CFC, but not in the striatum, a brain region not primarily implicated in contextual memory formation. Furthermore, using a combination of bulk and single cell RNA sequencing, we find that synaptic plasticity-promoting gene expression cascades are more strongly engaged in the hippocampus than in the striatum, but only when HDACi treatment co-occurred with CFC, and not by either treatment alone. Lastly, using ChIP-sequencing, we show that the combined action of HDACi application and conditioning is required to elicit enhancer histone acetylation in pathways that may underlie improved memory performance. Together, our results indicate that systemic HDACi administration amplifies brain-region specific processes that are naturally induced by learning. These findings shed light onto the mode of action of HDACis as cognitive enhancers.

Significance StatementMemory formation relies on a plethora of functions, including epigenetic modifications. Over the past years, multiple studies have indicated the potential of HDAC inhibitors (HDACi) to act as cognitive enhancers, but their mode of action is not fully understood. Here, we tested whether HDACi treatment improves memory formation via "cognitive epigenetic priming", stipulating that HDACis - without inherent target specificity - specifically enhance plasticity-related processes. We found that combining HDACi with fear learning, but not either treatment alone, enhances synaptic plasticity as well as memory-promoting transcriptional signaling in the hippocampus, a brain area known to be recruited by fear learning, but not in others. These results lend experimental support to the theory of "cognitive epigenetic priming".

category: neuroscience


title: Neural signatures of compression in the retelling of past events

abstract: When we retell our past experiences, we aim to reproduce some version of the original events; this reproduced version is often temporally compressed relative to the original. How does such compression of memories manifest in brain activity? One possibility is that a compressed retrieved memory manifests as a neural pattern which is more dissimilar to the original, relative to a more detailed or vivid memory. However, we argue that measuring raw dissimilarity alone is insufficient, as it confuses a variety of interesting and uninteresting changes. To address this problem, we examine brain pattern changes that are consistent across people. We show that temporal compression in individuals retelling of past events predicts systematic encoding-to-recall transformations in a number of higher associative regions. These findings elucidate how neural representations are not simply reactivated, but can also be transformed due to temporal compression during a universal form of human memory expression: verbal retelling.


title: Mitophagy enhancers against phosphorylated Tau-induced mitochondrial and synaptic toxicities in Alzheimer disease

abstract: The purpose of our study is to determine the protective effects of mitophagy enhancers against phosphorylated tau (P-tau)-induced mitochondrial and synaptic toxicities in Alzheimers disease (AD). Mitochondrial abnormalities, including defective mitochondrial dynamics, biogenesis, axonal transport and impaired clearance of dead mitochondria are linked to P-tau in AD. Mitophagy enhancers are potential therapeutic candidates to clear dead mitochondria and improve synaptic and cognitive functions in AD. We recently optimized the doses of mitophagy enhancers urolithin A, actinonin, tomatidine, nicotinamide riboside in immortalized mouse primary hippocampal (HT22) neurons. In the current study, we treated mutant Tau expressed in HT22 (mTau-HT22) cells with mitophagy enhancers and assessed mRNA and protein levels of mitochondrial/synaptic genes, cell survival and mitochondrial respiration. We also assessed mitochondrial morphology in mTau-HT22 cells treated and untreated with mitophagy enhancers. Mutant Tau-HT22 cells showed increased fission, decreased fusion, synaptic & mitophagy genes, reduced cell survival and defective mitochondrial respiration. However, these events were reversed in mitophagy enhancers treated mTau-HT22 cells. Cell survival was increased, mRNA and protein levels of mitochondrial fusion, synaptic and mitophagy genes were increased, and mitochondrial fragmentation is reduced in mitophagy enhancers treated mTau-HT22 cells. Further, urolithin A showed strongest protective effects among all enhancers tested in AD. Our combination treatments of urolithin A + EGCG, addition to urolithin A and EGCG individual treatment revealed that combination treatments approach is even stronger than urolithin A treatment. Based on these findings, we cautiously propose that mitophagy enhancers are promising therapeutic drugs to treat mitophagy in patients with AD.


title: High permeability of the CSF flow in the juvenile mouse brain by brain-wide DCE-MRI

abstract: The vasculature system with a lymphatic function in the brain is manifested as meningeal lymphatic vessels and the glymphatic system, which drain waste products from cerebrospinal fluid (CSF) and produce interstitial fluid. Invasion of stimulated immune cells or inflammatory cytokines during the maturation is regarded as a sign of the disruption of brain functions through excessive immune stresses. However, it is unclear when the lymphatic system is functionally completed and in which parts of the organ the brain immunity is privileged (i.e., the blood-brain barrier is fully established). Here, we visualized the whole-head signal of gadolinium (Gd) contrast agents in mouse brains to investigate the CSF flow from the signals of Gd-contrast agents using magnetic resonance imaging (MRI) in the course of developmental stages. We found higher Gd-signals in the olfactory bulbs, prefrontal cortex, parietal surface regions of the neocortex along with dorsal meningeal lymphatic vessels, ventral midbrain, and a part of cerebellum, as well as in the basal brain regions, at the immature stage of postnatal (P) 4 weeks, compared to the P8-12 weeks. Our results suggest that the barrier of the vasculature system in mouse brains is still permeable until P1-month and small molecules can leak into the parenchyma.


title: Individual differences in behavioral flexibility predict future volitional ethanol consumption in mice

abstract:
Cognitive control is key to regulating alcohol intake and preventing relapse. Behavioral inflexibility can prevent adaptive strategies such as mindfulness or other relapse-prevention behaviors. In a mouse model we investigated whether individual variability in behavioral flexibility (using attentional set-shifting task; ASST) predicts future alcohol intake. Adult male and female C57BL/6J mice were subjected to ASST using a bowl digging paradigm where mice identify a baited bowl based on compound odor and textural cues. This was completed prior to any alcohol exposure. Individual performance across mice varied within the group. We integrated several metrics, specifically ASST stage completed, trials to completion and errors performed to produce an individual performance index measure of behavioral flexibility. After, ASST mice were trained to drink ethanol (15%, v/v, 1hr/day) for 3-4 weeks until intake stabilized. Using this prospective approach, we identified an inverse relationship between behavioral flexibility and drinking - less flexible mice had a propensity to consume more alcohol. Similar relationships have been identified previously in non-human primates and rats. Our results show that the relationship between alcohol and cognitive flexibility is a robust trait that is conserved across species and can be used in mice to study neural substrates underlying these behaviors.

HIGHLIGHTS- ASST can be used to examine individual differences in cognitive function in mice
- Behavioral inflexibility is related to future higher alcohol consumption
- Executive function can be used as a predictive risk factor for alcohol intake
- This relationship in mice supports previous findings across species


title: Preparing for a second attack: a lesion simulation study on network resilience after stroke

abstract:
Background and PurposeDoes the brain become more resilient after a first stroke to reduce the consequences of a new lesion? Although recurrent strokes are a major clinical issue, whether and how the brain prepares for a second attack is unknown. This is due to the difficulties to obtain an appropriate dataset of stroke patients with comparable lesions, imaged at the same interval after onset. Furthermore, timing of the recurrent event remains unpredictable.

MethodsHere we used a novel clinical lesion simulation approach to test the hypothesis that resilience in brain networks increases during stroke recovery. 16 patients with a lesion restricted to the primary motor cortex were recruited. At 3 time points of the index event (10 days, 3 weeks, 3 months), we mimicked recurrent infarcts by deletion of nodes in brain networks (resting-state fMRI). Graph measures were applied to determine resilience (global efficiency) and wiring cost (mean degree) of the network.

ResultsAt 10 days and 3 weeks after stroke, resilience was similar in patients and controls. However, at 3 months, while motor function had fully recovered, resilience to clinically representative simulated lesions was higher compared to controls (cortical lesion p=0.012; subcortical: p=0.009; cortico-subcortical: p=0.009). Similar results were found after random (p=0.012) and targeted (p=0.015) attacks.

ConclusionOur results suggest that, after a lesion, brain networks reconfigure to increase resilience to future insults. Lesion simulation is an innovative approach, which may have major implications for stroke therapy. Individualized neuromodulation strategies could be developed to foster resilient network reconfigurations after a first stroke to limit the consequences of future attacks.


title: Explicit expression of mesophyll conductance in the traditional leaf photosynthesis-transpiration coupled model and its physiological significances

abstract:
O_LIAlmost all terrestrial biosphere models (TBMs) still assume infinite mesophyll conductance (gm) to estimate photosynthesis and transpiration. This assumption has caused low accuracy of TBMs to predict leaf gas exchange under certain conditions.
C_LIO_LIIn this study, we developed a photosynthesis-transpiration coupled model that explicitly considers gm and designed an optimized parameterization solution through evaluating four different gm estimation methods in 19 C3 species at 31 experimental treatments.
C_LIO_LIResults indicated that temperature responses of the maximum carboxylation rate (Fcmax) and the electron transport rate (Jmax) estimated by fusing the Bayesian retrieval algorithm and the Sharkey online calculator together with gm temperature response estimated by fusing the chlorophyll fluorescence-gas exchange method and anatomy method predicted leaf gas exchange more accurately. The gm temperature response exhibited activation energy ({Delta}Ha) of 63.13 {+/-} 36.89 kJ mol-1 and entropy ({Delta}S) of 654.49 {+/-} 11.36 J K-1 mol-1. The gm optimal temperature (Topt_gm) explained 58% of variations in photosynthesis optimal temperature (ToptA). The gm explicit expression has equally important effects on photosynthesis and transpiration estimations.
C_LIO_LIResults advanced understandings of better representation of plant photosynthesis and transpiration in TBMs.
C_LI

category: ecology


title: Origin of the type I antifreeze gene in flounders in response to Cenozoic climate change

abstract: Antifreeze proteins (AFPs) inhibit ice growth within fish and protect them from freezing in icy seawater. Alanine-rich, alpha-helical AFPs (type I) have independently (convergently) evolved in four branches of fishes, one of which is a subsection of the righteye flounders. The origin of this gene family has been elucidated by sequencing two loci from a starry flounder, Platichthys stellatus, collected off Vancouver Island, British Columbia. The first locus had two alleles that demonstrated the plasticity of the AFP gene family, one encoding 33 AFPs and the other allele only four. In the closely related Pacific halibut, this locus encodes multiple Gig2 (antiviral) proteins, but in the starry flounder, the Gig2 genes were found at a second locus due to a lineage-specific duplication event. An ancestral Gig2 gave rise to a 3-kDa "skin" AFP isoform, encoding three Ala-rich 11-a.a. repeats, that is expressed in skin and other peripheral tissues. Subsequent gene duplications, followed by internal duplications of the 11 a.a. repeat and the gain of a signal sequence, gave rise to circulating AFP isoforms. One of these, the "hyperactive" 32-kDa Maxi likely underwent a contraction to a shorter 3.3-kDa "liver" isoform. Present day starry flounders found in Pacific Rim coastal waters from California to Alaska show a positive correlation between latitude and AFP gene dosage, with the shorter allele being more prevalent at lower latitudes. This study conclusively demonstrates that the flounder AFP arose from the Gig2 gene, so it is evolutionarily unrelated to the three other classes of type I AFPs from non-flounders. Additionally, this gene arose and underwent amplification coincident with the onset of ocean cooling during the Cenozoic ice ages.

category: evolutionary biology


title: Variation in intraspecific demography drives localised concordance but species-wide discordance in responses to Plio-Pleistocene climatic change

abstract: Understanding how species biology may facilitate resilience to climate change remains a critical factor in detecting and protecting species at risk of extinction. Many studies have focused on the role of particular ecological traits in driving species responses, but less so on demographic history and levels of standing genetic variation. We used environmental and genomic datasets to reconstruct the phylogeographic histories of two ecologically similar and largely co-distributed freshwater fishes to assess the degree of concordance in their responses to Plio-Pleistocene climatic changes. Although several co-occurring populations demonstrated concordant demographic histories, idiosyncratic population size changes were found at the range edges of the more spatially restricted species. Discordant responses between species were associated with low standing genetic variation in peripheral populations. This might have hindered adaptive potential, as documented in recent population declines and extinctions of the two species. Our results highlight both the role of spatial scale in the degree of concordance in species responses to climate change, and the importance of standing genetic variation in facilitating range shifts. Even when ecological traits are similar between species, long-term genetic diversity and historical population demography may lead to discordant responses to ongoing and future climate change.


title: Comparison of breeding strategies for the creation of a synthetic pig line

abstract: Creating a new synthetic line by crossbreeding means complementary traits from pure breeds can be combined in the new population. Although diversity is generated during the crossbreeding stage, in this study, we analyze diversity management before selection starts. Using genomic and phenotypic data from animals belonging to the first generation (G0) of a new line, different simulations were run to evaluate diversity management during the first generations of a new line and to test the effects of starting selection at two alternative times, G3 and G4. Genetic diversity was characterized by allele frequency, inbreeding coefficients based on genomic and pedigree data, and expected heterozygosity. Breeding values were extracted at each generation to evaluate differences in starting selection at G3 or G4. All simulations were run for ten generations. A scenario with genomic data to manage diversity during the first generations of a new line was compared with a random and a selection scenario. As expected, loss of diversity was higher in the selection scenario, while the scenario with diversity control preserved diversity. We also combined the diversity management strategy with different selection scenarios involving different degrees of diversity control. Our simulation results show that a diversity management strategy combining genomic data with selection starting at G4 and a moderate degree of diversity control generates genetic progress and preserves diversity.

category: genetics


title: Automated 3D reconstruction of the fetal thorax in the standard atlas space from motion-corrupted MRI stacks for 21-36 weeks GA range

abstract:
Slice-to-volume registration (SVR) methods allow reconstruction of high-resolution 3D images from multiple motion-corrupted stacks. SVR-based pipelines have been increasingly used for motion correction for T2-weighted fetal MRI since they allow more informed and detailed diagnosis of brain and body anomalies including congenital heart defects (Lloyd et al., 2019). Recently, fully automated rigid SVR reconstruction of the fetal brain in the atlas space was achieved in (Salehi et al., 2019) that used convolutional neural networks (CNNs) for segmentation and pose estimation. However, these CNN-based methods have not yet been applied to the fetal trunk region. Meanwhile, the existing rigid and deformable SVR (DSVR) solutions (Uus et al., 2020) for the fetal trunk region are limited by the requirement of manual input as well the narrow capture range of the classical gradient descent based registration methods that cannot resolve severe fetal motion frequently occurring at the early gestational age (GA). Furthermore, in our experience, the conventional 2D slice-wise CNN-based brain masking solutions are reportedly prone to errors that require manual corrections when applied on a wide range of acquisition protocols or abnormal cases in clinical setting.

In this work, we propose a fully automated pipeline for reconstruction of the fetal thorax region for 21-36 weeks GA range T2-weighted MRI datasets. It includes 3D CNN-based intra-uterine localisation of the fetal trunk and landmark-guided pose estimation steps that allow automated DSVR reconstruction in the standard radiological space irrespective of the fetal trunk position or the regional stack coverage. The additional step for generation of the common template space and rejection of outliers provides the means for automated exclusion of stacks affected by low image quality or extreme motion. The pipeline was evaluated on a series of experiments including fetal MRI datasets and simulated rotation motion. Furthermore, we performed a qualitative assessment of the image reconstruction quality in terms of the definition of vascular structures on 100 early (median 23.14 weeks) and late (median 31.79 weeks) GA group MRI datasets covering 21 to 36 weeks GA range.


title: Mitochondrial stress-induced GDF15-GFRAL axis promotes anxiety-like behavior and CRH-dependent anorexia

abstract: Growth differentiation factor 15 (GDF15) is a stress-induced cytokine that modulates food intake and energy metabolism. Until now, most mechanistic studies on GDF15 rely on pharmacological interventions using exogenous GDF15, but little is known about its mode of action when induced both chronically and endogenously. Mitochondrial stress is one of the most described physiological conditions that induces GDF151, and therefore an important model to study the underlying mechanisms of endogenous GDF15s action. Here, using a mouse model of mitochondrial dysfunction via elevated respiratory uncoupling in skeletal muscle, we show a circadian oscillation of muscle-derived GDF15 to promote a daytime-restricted anorexia without signs of nausea or reduced physical activity, contrary to findings using recombinant GDF152-5. We find that mitochondrial stress-induced GDF15 associates with increased anxiety and hypothalamic corticotropin releasing hormone (CRH) induction, without further activation of the hypothalamic-pituitary-adrenal (HPA) axis and corticosterone response. Strikingly, the daytime-restricted anorexia, lean phenotype, systemic shift in substrate metabolism and anxiety-like behavior are completely abolished in conditions of mitochondrial stress coupled with genetic ablation of the GDF15 receptor GDNF receptor alpha-like (GFRAL), which is predominantly expressed in the hindbrain. Finally, we demonstrate that stress-induced GDF15-GFRAL signaling is required for hypothalamic CRH induction to control diurnal food intake in a CRH-receptor 1 (CRHR1)-dependent manner. With this, we uncover for the first time a molecular target of the GDF15-GFRAL axis that links anxiolytic and anorectic behavior as downstream effects of the chronic activation of this pathway by mitochondrial stress.


title: Urethane Anesthesia Exhibits Neurophysiological Correlates of Unconsciousness and is Distinct from Sleep

abstract: Urethane is a general anesthetic widely used in animal research. It is unique among anesthetics because urethane anesthesia alternates between macroscopically distinct electrographic states: a slow-wave state that resembles NREM sleep (NREMure), and an activated state with features of both REM sleep and wakefulness (REMure). However, the relationship between urethane anesthesia and physiological sleep is still unclear. In this study, electroencephalography (EEG) and electromyography were recorded in chronically prepared rats during natural sleep-wake states and during urethane anesthesia. We subsequently analyzed the EEG signatures associated with the loss of consciousness and found that, in comparison to natural sleep-wake states, the power, coherence, directed connectivity and complexity of brain oscillations are distinct during urethane. We also demonstrate that both urethane states have clear EEG signatures of general anesthesia. Thus, despite superficial similarities that have led others to conclude that urethane is a model of sleep, the electrocortical traits of depressed and activated states during urethane anesthesia differ from physiological sleep states.


title: Long- and short-term history effects in a spiking network model of statistical learning

abstract: The statistical structure of the environment is often important when making decisions. There are multiple theories of how the brain represents statistical structure. One such theory states that neural activity spontaneously samples from probability distributions. In other words, the network spends more time in states which encode high-probability stimuli. Existing spiking network models implementing sampling lack the ability to learn the statistical structure from observed stimuli and instead often hard-code a dynamics. Here, we focus on how arbitrary prior knowledge about the external world can both be learned and spontaneously recollected. We present a model based upon learning the inverse of the cumulative distribution function. Learning is entirely unsupervised using biophysical neurons and biologically plausible learning rules. We show how this prior knowledge can then be accessed to compute expectations and signal surprise in downstream networks. Sensory history effects emerge from the model as a consequence of ongoing learning.


title: Chicken Auditory Supporting Cells Express Interferon Response Genes during Regeneration towards Nascent Sensory Hair Cells In Vivo

abstract: The avian hearing organ is the basilar papilla that, in sharp contrast to the mammalian cochlea, can regenerate sensory hair cells and thereby recover from complete deafness within weeks. The mechanisms that trigger, sustain, and terminate the regenerative response in vivo are largely unknown. Here, we profile the changes in gene expression in the chicken basilar papilla after aminoglycoside antibiotic-induced hair cell loss using RNA-sequencing. The most prominent changes in gene expression were linked to the upregulation of interferon response genes which occurred in supporting cells, confirmed by single-cell RNA-sequencing and in situ hybridization. We determined that the JAK/STAT signaling pathway is essential for the interferon gene response in supporting cells, set in motion by hair cell loss. Four days after ototoxic damage, we identified newly regenerated, nascent auditory hair cells that express genes linked to termination of the interferon response. These cells are incipient modified neurons that represent a population of hair cells en route towards obtaining their location-specific and fully functional cell identity. The robust, transient expression of immune-related genes in supporting cells suggests a potential functional involvement of JAK/STAT signaling and interferon in sensory hair cell regeneration.


title: Modulation of sensory perception by hydrogen peroxide enables Caenorhabditis elegans to find a niche that provides both food and protection from hydrogen peroxide

abstract: Hydrogen peroxide (H2O2) is the most common chemical threat that organisms face. Here, we show that H2O2 alters the bacterial food preference of Caenorhabditis elegans, enabling the nematodes to find a safe environment with food. H2O2 induces the nematodes to leave food patches of laboratory and microbiome bacteria when those bacterial communities have insufficient H2O2-degrading capacity. The nematodes behavior is directed by H2O2-sensing neurons that promote escape from H2O2 and by bacteria-sensing neurons that promote attraction to bacteria. However, the input for H2O2-sensing neurons is removed by bacterial H2O2-degrading enzymes and the bacteria-sensing neurons perception of bacteria is prevented by H2O2. The resulting cross-attenuation provides a general mechanism that ensures the nematodes behavior is faithful to the lethal threat of hydrogen peroxide, increasing the nematodes chances of finding a niche that provides both food and protection from hydrogen peroxide.


title: Transcranial Focused Ultrasound Neuromodulation of the Thalamic Visual Pathway in a Large Animal Model

abstract:
BackgroundNeuromodulation of deep brain structures via transcranial ultrasound stimulation (TUS) is a promising, but still elusive approach to non-invasive treatment of brain disorders.

ObjectiveThe purpose of this study was to determine whether MR-guided TUS of the lateral geniculate nucleus (LGN) can modulate visual evoked potentials (VEPs) and cortical brain oscillations in a large animal model.

MethodsThe lateral geniculate nucleus (LGN) on one side was identified with T2-weighted MRI in sheep (all male, n=9). MR acoustic radiation force imaging (MR-ARFI) was used to confirm a tight sonication focus. Electroencephalographic (EEG) signals were recorded, and the visual evoked potential (VEP) peak-to-peak amplitude (N70 and P100) was calculated for each trial. Time-frequency spectral analysis was performed to elucidate the effect of TUS on cortical brain dynamics.

ResultsThe VEP peak-to-peak amplitude was reversibly suppressed relative to baseline during TUS. Dynamic spectral analysis demonstrated a change in cortical oscillations when TUS is paired with visual sensory input, but not when TUS is applied by itself.

ConclusionTUS non-invasively delivered to LGN can neuromodulate visual activity and oscillatory dynamics in large mammalian brains.


title: Periodic attention operates faster during more complex visual search

abstract: Attention has been found to sample visual information periodically, in a wide range of frequencies below 20 Hz. This periodicity may be supported by brain oscillations at corresponding frequencies. We propose that part of the discrepancy in periodic frequencies observed in the literature is due to differences in attentional demands, resulting from heterogeneity in tasks performed. To test this hypothesis, we used visual search and manipulated task complexity, i.e., target discriminability (high, medium, low) and number of distractors (set size), while electro-encephalography was simultaneously recorded. We replicated previous results showing that the phase of pre-stimulus low-frequency oscillations predicts search performance. Crucially, such effects were observed at increasing frequencies within the theta-alpha range (6-18 Hz) for decreasing target discriminability. In medium and low discriminability conditions, correct responses were further associated with higher post-stimulus phase-locking than incorrect ones, in increasing frequency and latency. Finally, the larger the set size, the later the post-stimulus effect peaked. Together, these results suggest that increased complexity (lower discriminability or larger set size) requires more attentional cycles to perform the task, partially explaining discrepancies between reports of attentional sampling. Low-frequency oscillations structure the temporal dynamics of neural activity and aid top-down, attentional control for efficient visual processing.


title: Laryngeal and swallow dysregulation following acute cervical spinal cord injury

abstract: Proper function of the larynx is vital to airway protection, including swallow. While the swallow reflex is controlled by the brainstem, patients with cervical spinal cord injuries (cSCI) are likely at increased risk of disordered swallow (dysphagia) and pneumonia, and the underlying mechanisms are unknown. We aimed to determine if acute spinal cord injury would disrupt swallow function in animal models. We hypothesized that 1) loss of descending efferent information to the diaphragm would affect swallow and breathing differently, and that 2) loss of ascending spinal afferent information would alter central swallow regulation to change motor drive to the upper airway. We recorded amplitudes of laryngeal and inspiratory muscle electromyograms (EMGs), submental and pharyngeal muscle EMGs, and cardiorespiratory measures in freely breathing pentobarbital-anesthetized cats and rats. First, we assessed the effect of a lateral hemisection at the second cervical level (C2) in cats during breathing. Posterior cricoarytenoid (laryngeal abductor) EMG activity during inspiration increased nearly two-fold, indicating that inspiratory laryngeal drive increased following cSCI. Ipsilateral to the injury, the crural diaphragm EMG was significantly reduced during breathing (62 {+/-} 25 percent change post-injury), but no animal had a complete termination of all activity; 75% of animals had an increase in contralateral diaphragm recruitment after cSCI, but this did not reach significance. Next, we assessed the effect of C2 lateral hemisection in cats during swallow. The thyroarytenoid (laryngeal adductor) and thyropharyngeus (pharyngeal constrictor) both increased EMG activity during swallow, indicating increased upper airway drive during swallow following cSCI. There was no change in the number of swallows stimulated per trial. We also found that diaphragm activity during swallow (schluckatmung) was bilaterally suppressed after lateral C2 hemisection, which was unexpected because this injury did not suppress contralateral diaphragm activity during breathing. Swallow-breathing coordination was also affected by cSCI, with more post-injury swallows occurring during early expiration. Finally, because we wanted to determine if the chest wall is a major source of feedback for laryngeal regulation, we performed T1 total transections in rats. As in the cat C2 lateral hemisection, a similar increase in inspiratory laryngeal activity (posterior cricoarytenoid) was the first feature noted after rat T1 complete spinal cord transection. In contrast to the cat C2 lateral hemisection, diaphragmatic respiratory drive increased after T1 transection in every rat (215 {+/-} 63 percent change), and this effect was significant. Overall, we found that spinal cord injury alters laryngeal drive during swallow and breathing, and alters swallow-related diaphragm activity. Our results show behavior-specific effects, suggesting that swallow may be more affected than breathing is by cSCI, and emphasizing the need for additional studies on laryngeal function during breathing and swallow after spinal cord injury.


title: Extrasynaptic signaling enables an asymmetric juvenile motor circuit to produce a symmetric mature gait.

abstract: Bilaterians generate motor patterns with symmetries that correspond to their body plans. This is thought to arise from wiring symmetries in their motor circuits. We show that juvenile C. elegans larva has an asymmetrically wired motor circuit, but still generates a bending pattern with dorsal-ventral symmetry. In the juvenile circuit, wiring between excitatory and inhibitory motor neurons coordinates contraction of dorsal muscles with relaxation of ventral muscles, producing dorsal bends. Ventral bending is not driven by analogous wiring. Instead, ventral muscles are excited uniformly by premotor interneurons through extrasynaptic signaling. Ventral bends occur in anti-phasic entrainment to activity of the same motor neurons that drive dorsal bends. During maturation, the juvenile motor circuit is replaced by two homologous motor circuits that separately drive dorsal and ventral bending. Modeling reveals that the juveniles immature motor circuit provides an adequate solution to generate an adult-like gait long before the animal matures. Circuit degeneracy minimizes functional disruption during development.


title: Extrasynaptic signaling enables an asymmetric juvenile motor circuit to produce a symmetric gait


title: No evidence for aberrant expression of amyloid β and phosphorylated tau proteins in herpes simplex virus-infected human neurons in vivo

abstract: Increasing evidence implicates the neurotropic herpes simplex virus 1 (HSV-1) in the pathogenesis of Alzheimers Disease (AD). However, it is unclear whether previously reported findings in HSV-1 cell culture and animal models can be translated to humans. Here, we analyzed clinical specimens from latently HSV-1 infected individuals and individuals with lytic HSV infection of the brain (herpes simplex encephalitis; HSE). Latent HSV-1 DNA load in trigeminal ganglia was identical between AD patients and controls, and latently HSV-infected neurons did not express amyloid {beta} (A{beta}) or hyperphosphorylated tau (pTau). Some HSE patient brains presented with ageing-related intraneuronal A{beta} accumulations, neurofibrillary tangles (NFT) or extracellular A{beta} plaques, but these were neither restricted to HSV-infected neurons nor brain regions containing virus-infected cells. Analysis of unique brain material from an AD patient with concurrent HSE showed that HSV-infected cells frequently localized close to A{beta} plaques and NFT, but did not exacerbate AD-related pathology in relation to HSV infection. HSE-associated neuroinflammation was not associated with specific A{beta} or pTau phenotypes. Collectively, the data indicate that neither latent nor lytic HSV infection of human neurons in vivo is directly associated with aberrant A{beta} or pTau protein expression.

category: microbiology


title: High-throughput functionalization of the Toxoplasma kinome uncovers a novel regulator of invasion and egress

abstract: Protein kinases regulate fundamental aspects of cell biology in all eukaryotes, making them attractive chemotherapeutic targets in Apicomplexan parasites such as the causative agents of malaria (Plasmodium spp.) and toxoplasmosis (Toxoplasma gondii). However, the precise roles of individual parasite kinases cannot be inferred simply from sequence identity, due to rewiring of signaling pathways and the shifting repertoire of kinases across species. To systematically examine the parasite kinome, we developed a high-throughput (HiT) CRISPR-mediated tagging strategy to endogenously label all predicted cytosolic protein kinases with a synthetic sequence encoding the minimal auxin-inducible degron (mAID) linked to a fluorophore and epitope tag. The system enables the assembly of thousands of tagging vectors from synthetic sequences in a single reaction and the pooled generation of mutants to examine kinase localization and function. We examined the phenotypes associated with kinase knock-down in 1,160 arrayed clones by replica-plating in the presence or absence of auxin and found broad defects across the lytic cycle for 109 clonal isolates, assigning localizations to 39 proteins, and associating 15 kinases within 6 distinct morphological phenotypes. The relative fitness of tagged alleles was also examined by tracking the relative abundance of individual guide RNAs as parasite populations progressed through the lytic cycle, in the presence or absence of auxin. Pooled screening had a high predictive value and differentiated between delayed and acute death. Demonstrating the value of this resource, we identified a novel kinase associated with delayed death as a novel regulator of invasion and egress. We call the previously unstudied kinase Store Potentiating/Activating Regulatory Kinase (SPARK), based on its impact on intracellular Ca2+ stores at key moments during the lytic cycle. Despite having a similar kinase domain to the mammalian PDK1, SPARK lacks the canonical lipid-binding domain and we find no indication SPARK positively regulates other AGC kinases, suggesting a rewiring of signaling pathways to accommodate parasite adaptations. The HiT vector screening system extends the applications of genome-wide screens into complex cellular phenotypes, providing a scalable and versatile platform for the dissection of apicomplexan cell biology.


title: Transcriptional dynamics of methane-cycling microbiomes are linked to seasonal CH4 fluxes in two hydromorphic and organic-rich grassland soils

abstract: Soil CH4 fluxes are driven by CH4-producing and -consuming microorganisms that determine whether soils are sources or sinks of this potent greenhouse gas. Using quantitative metatranscriptomics, we linked CH4-cycling microbiomes to net surface CH4 fluxes throughout a year in two drained peatland soils differing in grassland land-use intensity and physicochemical properties. CH4 fluxes were highly dynamic; both soils were net CH4 sources in autumn and winter and sinks in spring and summer. Despite similar net CH4 emissions, methanogen and methanotroph loads, as determined by small subunit rRNA transcripts per gram soil, differed strongly between sites. In contrast, mRNA transcript abundances were similar in both soils and correlated well with CH4 fluxes. The methane monooxygenase to methanogenesis mRNA ratio was higher in spring and summer, when the soils were net CH4 sinks. CH4 uptake was linked to an increased proportion of USC and {gamma} and pmoA2 pmoA transcripts. We assume that methanogen transcript abundance may be useful to approximate changes in net surface CH4 emissions from drained peat soils; high methanotroph to methanogen ratios would indicate CH4 sink properties. Our study shows the strength of quantitative metatranscriptomics; mRNA transcript abundance holds promising indicator to link soil microbiome functions to ecosystem-level processes.


title: Maturation of human intestinal epithelium from pluripotency in vitro

abstract: Methods to generate human intestinal tissue from pluripotent stem cells (PSCs) open new inroads into modeling intestine development and disease. However, current protocols require organoid transplantation into an immunocompromised mouse to achieve matured and differentiated epithelial cell states. Inspired by developmental reconstructions from primary tissues, we establish a regimen of inductive cues that enable stem cell maturation and epithelial differentiation entirely in vitro. We show that the niche factor Neuregulin1 (NRG1) promotes morphological change from proliferative epithelial cysts to matured epithelial tissue in three-dimensional cultures. Single-cell transcriptome analyses reveal differentiated epithelial cell populations, including diverse secretory and absorptive lineages. Comparison to multi-organ developmental and adult intestinal cell atlases confirm the specificity and maturation state of cell populations. Altogether, this work opens a new direction to use in vitro matured epithelium from human PSCs to study human intestinal epithelium development, disease, and evolution in controlled culture environments.

category: cell biology


title: How the naked mole-rat resists senescence: a constraints-based theory

abstract: Here, I present a theory describing how the stabilization of constraints imposed on chromatin dynamics by the naked mole-rats histone H1.0 protein--which in terminally differentiated cells constrains the accessibility of the nucleosome core particle for histone-modifying enzymes and chromatin remodeling factors--explains its resistance to both senescence and cancer. Further, this theory predicts that a mutant house mouse displaying such stabilization will be similarly resistant to both senescence and cancer. A proof-of-concept computational analysis is presented and two predictions for the direct testing of the theory are provided. These experiments comprise, as test subjects, mutant naked mole-rats synthesizing a house mouse (Mus musculus)-like histone H1.0, and mutant house mice synthesizing a naked mole-rat-like histone H1.0. The predictions are that the constraints on chromatin dynamics embodied by the respective mutant histone H1.0 proteins will negate the otherwise significant resistance to both senescence and cancer of the naked mole-rats and, conversely, confer such resistance to the house mice. A verification of these predictions will imply that constraints on chromatin dynamics embodied by naked mole-rat-like histone H1.0 proteins may confer significant resistance to both senescence and age-related cancer to otherwise senescence-prone and/or cancer-susceptible multicellular species, including humans.

category: developmental biology


title: Prescribed versus wildfire impacts on exotic plants and soil microbes in California grasslands

abstract:
O_LIPrescribed fire is often used as a management tool to decrease exotic plant cover and increase native plant cover in grasslands. These changes may also be mediated by fire impacts on soil microbial communities, which drive plant productivity and function. Yet, the ecological effects of prescribed burns compared to wildfires on either plant or soil microbial composition remain unclear.
C_LIO_LIHere, we investigated the impacts of a spring prescribed fire versus a fall wildfire on plant cover and community composition and bacterial and fungal richness, abundance, and composition in a California grassland. We used qPCR of 16S and 18S to assess impacts on bacterial and fungal abundance and Illumina MiSeq of 16S and ITS2 to assess impacts on bacterial and fungal richness and composition.
C_LIO_LIWildfire had stronger impacts than prescribed fire on microbial communities and both fires had similar impacts on plants with both prescribed and wildfire reducing exotic plant cover but neither reducing exotic plant richness. Fungal richness declined after the wildfire but not prescribed fire, but bacterial richness was unaffected by either. Yet increasing char levels in both fire types resulted in reduced bacterial and fungal biomass, and both fire types slightly altered bacterial and fungal composition.
C_LIO_LIExotic and native plant diversity differentially affected soil microbial diversity, with native plant diversity leading to increased arbuscular mycorrhizal fungal richness while exotic plant diversity better predicted bacterial richness. However, the remainder of the soil microbial communities were more related to aspects of soil chemistry including cation exchange capacity, organic matter, pH and phosphorous.
C_LIO_LISynthesis and applications. Understanding the different ecological effects of prescribed fires and wildfires on plant and soil communities are key to enhancing a prevalent management action and to guide potential management opportunities post wildfires. Our coupled plant and soil community sampling allowed us to capture the sensitivity of the fungal community to fire and highlights the importance of potentially incorporating management actions such as soil or fungal amendments to promote this critical community that mediates native plant performance.
C_LI


title: On the Origin of Feces: Fungal diversity, distribution, and conservation implications from feces of small mammals

abstract: Fungi are extremely diverse, but only a small fraction of the total estimated species have been characterized. Often, the extent of diversity and distribution of fungal communities is difficult or near impossible to assess. This is due to the fact that many Fungi are cryptic and persist predominantly hidden within substrates such as soil or plant material. This is particularly true for hypogeous sporocarps, including truffles, which are extremely difficult to survey in a systematic manner. However, fungi with hypogeous sporocarps have evolved traits that make them highly attractive to animals, such as small mammals, which ingest and disperse fungal spores through defecation. Here, samples of feces from 138 small mammals collected in the western United States were assessed for total fungal diversity using a dual-index metabarcoding, high-throughput Illumina sequencing approach. Our findings exhibit many identifications within Agaricomycetidae, with 65 of the 138 samples containing sequences belonging to several species of the hypogeous truffle genus Rhizopogon. Metadata, such as geospatial coordinates, for each captured small mammal can be used as a proxy for the presence or absence of Rhizopogon species identified in their feces. Utilizing these proxy data, along with publicly available data on observation and occurrence from data repositories such as GBIF and MycoPortal, it is possible to augment our current knowledge of the extent of occurrence and area of occupancy of cryptic hypogeous fungi without direct observation, further enhancing our ability to assess their conservation status.


title: Elevated alpha diversity in disturbed sites obscures regional decline and homogenization of amphibian diversity

abstract: Loss of natural habitat is one of the major threats for biodiversity worldwide. Habitat conversion not only changes diversity and species composition locally (alpha diversity) but might also lead to large-scale homogenization of species communities and decrease in regional species richness (gamma diversity). We investigated the effect of farmland conversion on amphibian communities in Rwanda and compared local and regional (country-wide) taxonomic, functional and phylogenetic diversity between natural and farmland sites (agricultural marais). Alpha diversity was higher in the disturbed farmland than in natural sites. However, species turnover among farmland sites was much lower than among natural sites, resulting in highly homogenized amphibian communities and much lower country-wide taxonomic, functional and phylogenetic gamma diversity in farmland compared to natural sites. The few frog species found in farmland were mostly disturbance-tolerant species that are widespread in Eastern Africa and beyond. In contrast, most of the regionally endemic frog species that make this region a continent-scale hotspot of amphibian diversity were found only in the natural habitats. Ongoing farmland conversion might lead to a loss of regional endemism and a widespread homogenization of species communities across sub-Saharan Africa.


title: The Influence of Model Violation on Phylogenetic Inference: A Simulation Study

abstract: AO_SCPLOWBSTRACTC_SCPLOWPhylogenetic inference typically assumes that the data has evolved under Stationary, Reversible and Homogeneous (SRH) conditions. Many empirical and simulation studies have shown that assuming SRH conditions can lead to significant errors in phylogenetic inference when the data violates these assumptions. Yet, many simulation studies focused on extreme non-SRH conditions that represent worst-case scenarios and not the average empirical dataset. In this study, we simulate datasets under various degrees of non-SRH conditions using empirically derived parameters to mimic real data and examine the effects of incorrectly assuming SRH conditions on inferring phylogenies. Our results show that maximum likelihood inference is generally quite robust to a wide range of SRH model violations but is inaccurate under extreme convergent evolution.


title: Idiosyncratic epistasis leads to global fitness-correlated trends

abstract:
Epistasis can dramatically affect evolutionary trajectories. In recent decades, protein-level fitness landscapes have revealed extensive idiosyncratic epistasis among specific mutations. In contrast, other work has found ubiquitous and apparently non-specific patterns of global diminishing-returns and increasing-costs epistasis among mutations across the genome. Here, we use a hierarchical CRISPR gene drive system to construct all combinations of 10 missense mutations from across the genome in budding yeast, and measure their fitness in six environments. We show that the resulting fitness landscapes exhibit global fitness-correlated trends, but that these trends emerge from specific idiosyncratic interactions. This provides the first experimental validation of recent theoretical work that has argued that fitness-correlated trends can emerge as the generic consequence of idiosyncratic epistasis.

One-Sentence SummaryA genome-spanning fitness landscape reveals how idiosyncratic genetic interactions lead to global epistatic patterns.


title: Deleterious Mutations Accumulate Faster in Allopolyploid than Diploid Cotton (Gossypium) and Unequally Between Subgenomes

abstract: Whole genome duplication (polyploidization) is among the most dramatic mutational processes in nature, so understanding how natural selection differs in polyploids relative to diploids is an important goal. Population genetics theory predicts that recessive deleterious mutations accumulate faster in allopolyploids than diploids due to the masking effect of redundant gene copies, but this prediction is hitherto unconfirmed. Here, we use the cotton genus (Gossypium), which contains seven allopolyploids derived from a single polyploidization event 1-2 million years ago, to investigate deleterious mutation accumulation. We use two methods of identifying deleterious mutations at the nucleotide and amino acid level, along with whole-genome resequencing of 43 individuals spanning six allopolyploid species and their two diploid progenitors, to demonstrate that deleterious mutations accumulate faster in allopolyploids than in their diploid progenitors. We find that, unlike what would be expected under models of demographic changes alone, strongly deleterious mutations show the biggest difference between ploidy levels, and this effect diminishes for moderately and mildly deleterious mutations. We further show that the proportion of nonsynonymous mutations that are deleterious differs between the two co-resident subgenomes in the allopolyploids, suggesting that homoeologous masking acts unequally between subgenomes. Our results provide a genome-wide perspective on classic notions of the significance of gene duplication that likely are broadly applicable to allopolyploids, with implications for our understanding of the evolutionary fate of deleterious mutations. Finally, we note that some measures of selection (e.g. dN/dS,{pi} N/{pi}S) may be biased when species of different ploidy levels are compared.


title: Zebrafish Information Network, the knowledgebase for Danio rerio research

abstract: The Zebrafish Information Network (ZFIN, zfin.org) is the central repository for zebrafish genetic and genomic data. ZFIN expertly curates, integrates, and displays zebrafish data including genes, alleles, human disease models, gene expression, phenotype, gene function, orthology, morpholino, CRISPR, TALEN, and antibodies. ZFIN makes zebrafish research data Findable, Accessible, Interoperable, and Reusable (FAIR) through nomenclature, curatorial and annotation activities, web interfaces, and data downloads. ZFIN is a founding member of the Alliance of Genome Resources, providing zebrafish data for integration into the cross species platform as well as contributing to model organism data harmonization efforts.


title: Expedited gene delivery for osteochondral defect repair in a rabbit knee model: a one-year investigation

abstract:
ObjectiveTo evaluate a single-step, gene-based procedure for repairing osteochondral lesions.

DesignOsteochondral lesions were created in the patellar groove of skeletally mature rabbits. Autologous bone marrow aspirates were mixed with adenovirus vectors carrying cDNA encoding green fluorescent protein (Ad.GFP) or transforming growth factor-{beta}1 (Ad.TGF-{beta}1) and allowed to clot. The clotted marrow was press-fit into the defects. Animals receiving Ad.GFP were euthanized at 2 weeks and intra-articular expression of GFP examined by fluorescence microscopy. Animals receiving Ad.TGF-{beta}1 were euthanized at 3 months and 12 months; repair was compared to empty defects using histology and immunohistochemistry. Complementary in vitro experiments assessed transgene expression and chondrogenesis in marrow clots and fibrin gels. In a subsequent pilot study, repair at 3 months using a fibrin gel to encapsulate Ad.TGF-{beta}1 was evaluated.

ResultsAt 2 weeks, GFP expression was seen at variable levels within the cartilaginous lesion. At 3 months, there was a statistically significant improvement in healing of lesions receiving Ad.TGF-{beta}1, although variability was high. At 12 months, there was no difference between the empty defects and those receiving Ad.TGF-{beta}1 in overall score and cartilage score, but the bone healing score remained higher. Variability was again high. In vitro experiments suggested that variability reflected variable transduction efficiency and chondrogenic activity of the marrow clots; using fibrin gels instead of marrow provided more uniformity in healing.

ConclusionsThis approach to improving the repair of osteochondral lesions holds promise but needs further refinement to reduce variability and provide a more robust outcome.


title: AMPK regulates germline stem cell integrity and quiescence through a mir-1/tbc-7/rab-7 pathway in C. elegans.

abstract: During periods of energetic stress, Caenorhabditis elegans can undergo a global quiescent stage known as "dauer". During this stage, all germline stem cells undergo G2 cell cycle arrest through an AMPK-dependent mechanism. In animals that lack AMPK signalling, the germ cells fail to arrest, undergo uncontrolled proliferation and lose their reproductive capacity. These germline defects are accompanied by an altered chromatin landscape and gene expression program. We identified an allele of tbc-7, a RabGAP protein that functions in the neurons, which when compromised, suppresses the germline hyperplasia in the dauer larvae, as well as the post-dauer sterility and somatic defects characteristic of AMPK mutants. This mutation also corrects the abundance and aberrant distribution of transcriptionally activating and repressive chromatin marks in animals that otherwise lack all AMPK signalling. We identified RAB-7 as one of the potential RAB proteins that is regulated by tbc-7 and show that the activity of RAB-7 is critical for the maintenance of germ cell integrity during the dauer stage. A singular small RNA, mir-1, was identified as a direct negative regulator of tbc-7 expression through the analysis of seed sequences on the 3UTR of tbc-7. Animals lacking mir-1 are post-dauer sterile, displaying a similar phenotype to AMPK mutants. Altogether, our findings describe a novel mir-1/tbc-7/rab-7 pathway occurring in the neurons that regulates the germ line in a cell non-autonomous manner.


title: The genome of an apodid holothuroid (Chiridota heheva) provides insights into its adaptation to deep-sea reducing environment

abstract: Cold seeps and hydrothermal vents are deep-sea reducing environments that are characterized by a lack of oxygen, photosynthesis-derived nutrients and a high concentration of reducing chemicals. Apodida is an order of deep-sea echinoderms lacking tube feet and complex respiratory trees, which are commonly found in holothurians. Chiridota heheva Pawson & Vance, 2004 (Apodida: Chiridotidae) is one of the few echinoderms that resides in deep-sea reducing environments. Unlike most cold seep and hydrothermal vent-dwelling animals, C. heheva does not survive by maintaining an epi- or endosymbiotic relationship with chemosynthetic microorganisms. The species acquires nutrients by extracting organic components from sediment detritus and suspended material. Here, we report a high-quality genome of C. heheva as a genomic reference for echinoderm adaptation to reducing environments. Chiridota heheva likely colonized its current habitats in the early Miocene. The expansion of the aerolysin-like protein family in C. heheva compared with other echinoderms might be involved in the disintegration of microbes during digestion, which in turn facilitates the species adaptation to cold seep environments. Moreover, several hypoxia-related genes were subject to positive selection in the genome of C. heheva, which contributes to their adaptation to hypoxic environments.

category: genomics


title: A population based expression atlas provides insights into disease resistance and other physiological traits in cassava (Manihot esculenta Crantz)

abstract: Cassava, a food security crop in Africa, is grown throughout the tropics and subtropics. Although cassava can provide high productivity in suboptimal conditions, the yield in Africa is substantially lower than in other geographies. The yield gap is attributable to many challenges faced by cassava in Africa, including susceptibility to diseases and poor soil conditions. In this study, we carried out 3RNA sequencing on 150 accessions from the National Crops Resources Research Institute, Ugandan for 5 tissue types, providing population-based transcriptomics resources to the research community in a web-based queryable cassava expression atlas. Differential expression and weighted gene co-expression network analysis were performed to detect 8,820 significantly differentially expressed genes (DEGs), revealing similarity in expression patterns between tissue types and the clustering of detected DEGs into 18 gene modules. As a confirmation of data quality, differential expression and pathway analysis targeting cassava mosaic disease (CMD) identified 27 genes observed in the plant-pathogen interaction pathway, several previously identified CMD resistance genes and two peroxidase family proteins different from the CMD2 gene. Present research work represents a novel resource towards understanding complex traits at expression and molecular levels for the development of resistant and high-yielding cassava varieties, as exemplified with CMD.


title: Sexual dimorphism and the effect of wild introgressions on recombination in cassava (Manihot esculenta Crantz) breeding germplasm

abstract: Recombination has essential functions in meiosis, evolution, and breeding. The frequency and distribution of crossovers dictate the generation of new allele combinations and can vary across species and between sexes. Here, we examine recombination landscapes across the 18 chromosomes of cassava (Manihot esculenta Crantz) with respect to male and female meioses and known introgressions from the wild relative Manihot glaziovii. We used SHAPEIT2 and duoHMM to infer crossovers from genotyping-by-sequencing (GBS) data and a validated multi-generational pedigree from the International Institute of Tropical Agriculture (IITA) cassava breeding germplasm consisting of 7,020 informative meioses. We then constructed new genetic maps and compared them to an existing map previously constructed by the International Cassava Genetic Map Consortium (ICGMC). We observed higher recombination rates in females compared to males, and lower recombination rates in M. glaziovii introgression segments on chromosomes 1 and 4, with suppressed recombination along the entire length of the chromosome in the case of the chromosome 4 introgression. Finally, we discuss hypothesized mechanisms underlying our observations of heterochiasmy and crossover suppression and discuss the broader implications for plant breeding.


title: Prolactin and the regulation of parental care and helping behavior in cooperatively breeding white-browed sparrow weaver societies

abstract: In many cooperatively breeding societies non-breeding individuals help to rear the offspring of breeders. The physiological mechanisms that regulate such cooperative helping behavior are poorly understood, but may have been co-opted, during the evolution of cooperative breeding, from pre-existing mechanisms that regulated parental care. Key among these may be a role for prolactin. Here we investigate whether natural variation in circulating prolactin levels predicts both parental and helper contributions to nestling provisioning in cooperatively breeding white-browed sparrow weavers, Plocepasser mahali. In sparrow weaver groups, a single dominant pair monopolize reproduction and non-breeding subordinates help with nestling feeding. We show that: (i) among parents, dominant females feed nestlings at higher rates, make longer provisioning visits and have higher prolactin levels than dominant males; and (ii) among subordinates, engaged in cooperative helping behavior, those within their natal groups feed nestlings at higher rates and have higher prolactin levels than immigrants. Moreover, continuous variation in prolactin levels positively predicts nestling-provisioning rates and mean provisioning visit durations when all bird classes are combined. These relationships are principally driven by differences among bird classes in both circulating prolactin levels and provisioning traits; the more limited within-class variation in prolactin and provisioning traits were not evidently correlated, highlighting a likely role for additional mechanisms in the fine-scale regulation of care. Our findings broadly support the hypothesis that parental care and cooperative helping behavior are regulated by a common underlying mechanism and highlight the need for experimentation to now establish the causality of any role for prolactin.

category: animal behavior and cognition


title: Microbial Biotextiles for a Circular Materials Economy

abstract: A circular materials economy can disrupt linear production models and mitigate its environmental threats. We harness extracellular biosynthesis of microbial nanocellulose (MC), coupled to a protocol that represents a serendipitous intersection between ancient textile techniques and green processing as a route to high-performance bioleather with a circular life cycle. Specifically, a plant-based lecithin phosphocholine treatment modulates cellulose cross-linking through phosphate and methylene groups, to yield a biodegradable material with superior mechanical and flame-retardant properties. A broad color palette and controlled color modulation was achieved with cultural heritage dyes and waste-to-resource strategies. Life cycle impact assessment shows up to an order of magnitude reduction in environmental impacts relative to conventional textiles, making MC relevant for widespread application in fashion, interiors, construction and insulating materials. The translational potential of this approach is tremendous, as use of microbes and green chemistry to engineer regenerative, high-performance products can accelerate a paradigm shift to a circular bioeconomy.

category: bioengineering


title: Advanced Functional Biotextiles for a Circular Materials Economy


title: Anisotropic Rod-Shaped Particles Influence Injectable Granular Hydrogel Properties and Cell Invasion

abstract: Granular hydrogels have emerged as a new class of injectable and porous biomaterials that improve integration with host tissue when compared to solid hydrogels. Granular hydrogels are typically prepared using spherical particles and this study considers whether particle shape (i.e., isotropic spheres versus anisotropic rods) influences granular hydrogel properties and cellular invasion. Simulations predict that anisotropic rods influence pore shape and interconnectivity, as well as bead transport through granular assemblies. Photocrosslinkable norbornene-modified hyaluronic acid is used to produce spherical and rod-shaped particles using microfluidic droplet generators and formed into shear-thinning and self-healing granular hydrogels at low and high particle packing. Rod-shaped particles form granular hydrogels that have anisotropic and interconnected pores, with pore number and size, storage moduli, and extrusion forces influenced by particle shape and packing. Robust in vitro sprouting of endothelial cells from embedded cellular spheroids is observed with rod-shaped particles, including higher sprouting densities and sprout lengths when compared to hydrogels with spherical particles. Cellular invasion into granular hydrogels when injected subcutaneously in vivo is significantly greater with rod-shaped particles, whereas a gradient of cellularity is observed with spherical particles. Overall, this work demonstrates potentially superior functional properties of granular hydrogels with rod-shaped particles for tissue repair.


title: PcGCE is a potent elicitor of defense responses in aspen

abstract: Using microbial enzymes in transgenesis is a powerful means to introduce new functionalities in plants. Glucuronoyl esterase (GCE) is a microbial enzyme hydrolyzing the ester bond between lignin and 4-O-methyl--D-glucuronic acid present as a side chain of glucuronoxylan. This bond mediates lignin-carbohydrate complex (LCC) formation, considered as crucial factor of lignocellulose recalcitrance to saccharification. Previous studies showed that hybrid aspen (Populus tremula L. x tremuloides Michx.) constitutively expressing Phanerochaete carnosa Burt GCE (PcGCE) had better efficiency of cellulose-to-glucose conversion but were stunned and had lower cellulose content indicating that more studies are needed to design strategy for deployment of this enzyme in planta. Here we report that the transgenic plants exhibit premature leaf senescence, increased accumulation of calcium oxalate crystals, tyloses and necrotic lesions and have strongly activated immune defense reactions as revealed by their altered profiles of transcriptomes, metabolomes and hormones in the leaves. To elucidate if these effects are triggered by damage-associated molecular patterns (DAMPs) or by PcGCE protein perceived as a pathogen-associated molecular pattern (PAMP), we ectopically expressed in aspen an enzymatically inactive PcGCES217A. The mutated PcGCE induced similar growth retardation, leaf necrosis and premature senescence as the active one, providing evidence that PcGCE protein is recognized as PAMP. Transcriptomics analysis of young expanding leaves of 35S:PcGCE plants identified several candidates for receptors of PcGCE, which were not expressed in developing wood tissues. Grafting experiments showed that PcGCE transcripts are not cell-to-cell mobile and that PcGCE expressing leaves augment systemic responses. In agreement, expressing PcGCE in developing wood by using the wood-specific promoter (WP), avoided all off-target effects. Moreover, WP:PcGCE lines grew better than control plants providing evidence that this strategy can be used in transgenic crops dedicated for biorefinery.


title: Computationally-guided tuning of ligand sensitivity in a GPCR-based sensor

abstract: Genetically-encoded fluorescent sensors for neuromodulators are increasingly used molecular tools in neuroscience. However, these protein-based biosensors are often limited by the sensitivity of the protein scaffold towards endogenous ligands. Here, we explored the possibility of applying computational design approaches for enhancing sensor sensitivity. Using the dopamine sensor dLight1 as proof of concept, we designed two variants that boost the sensors potency (EC50) for dopamine and norepinephrine by up to 5- and 15-fold, respectively. Interestingly, the largest effects were obtained through improved designed allosteric transmission in the transmembrane region of the sensor. Our approach should prove generally useful for enhancing sensing capabilities of a large variety of neuromodulator sensors.


title: Engineered neuron-targeting, placental mesenchymal stromal cell-derived extracellular vesicles for in utero treatment of myelomeningocele

abstract: This study investigated the feasibility and efficiency of neuron-targeting hybrid placental mesenchymal stromal cell-derived extracellular vesicles (PMSC-EVs), engineered by membrane fusion with Targeted Axonal Import (TAxI) peptide modified, TrkB agonist 7,8-DHF-loaded liposomes for treatment of myelomeningocele (MMC) via intra-amniotic cavity administration. The prepared TAxI modified liposomes with 7,8-DHF were used to fuse with PMSC-EVs. Different fusion approaches were investigated and freeze-thaw-extrude method was found to be the optimal. The engineered PMSC-EVs had a uniform particle size and efficiently loaded 7,8-DHF. It also had typical markers of native EVs. Freeze-thaw-extrude process did not change the release profile of 7,8-DHF from engineered EVs compared to TAxI modified, 7,8-DHF loaded liposomes. The engineered EVs could elicit TrkB phosphorylation depending on the incorporation of 7,8-DHF while native EVs did not. The engineered EVs increased neurite outgrowth of apoptotic cortical neurons induced by staurosporine, suggesting that they exhibited neuroprotective function. In a rodent model of MMC, neuron-targeting, engineered EVs became an active targeting delivery system to MMC defect sites. Pups treated with engineered EVs had the lowest density of apoptotic cells and displayed a therapeutic outcome. The study suggests the potential use of engineered hybrid, active neuron-targeting EVs for the in utero treatment of MMC.


title: JaponicusDB: Rapid deployment of a model organism database for an emerging model species

abstract: The fission yeast Schizosaccharomyces japonicus has recently emerged as a powerful system for studying the evolution of essential cellular processes, drawing on similarities as well as key differences between S. japonicus and the related, well-established model Schizosaccharomyces pombe. We have deployed the open-source, modular code and tools originally developed for PomBase, the S. pombe model organism database (MOD), to create JaponicusDB (www.japonicusdb.org), a new MOD dedicated to S. japonicus. By providing a central resource with ready access to a growing body of experimental data, ontology-based curation, seamless browsing and querying, and the ability to integrate new data with existing knowledge, JaponicusDB supports fission yeast biologists to a far greater extent than any other source of S. japonicus data. JaponicusDB thus enables S. japonicus researchers to realise the full potential of studying a newly emerging model species, and illustrates the widely applicable power and utility of harnessing reusable PomBase code to build a comprehensive, community-maintainable repository of species-relevant knowledge.


title: Protein sequence profile prediction using ProtAlbert transformer1

abstract: Protein profiles have many applications in bioinformatics. To construct the profile from a protein sequence, the sequence is aligned with database. However, sometimes there are no similar sequences with the query. This paper proposes a method based on pre-trained ProtAlbert transformer to predict the profile for a single protein sequence without alignment. The performance of transformers on natural languages is impressive. Protein sequences can be viewed as a language; therefore, we can benefit from using these models. We analyze the attention heads in different layers of ProtAlbert to show that the transformer can capture five essential protein characteristics of the family from a single protein sequence. These assessments are performed on the CASP13 dataset to find representative heads for each of five protein characteristics. Then, these heads are investigated on one thermophilic and two mesophilic proteins as case studies. The results show the significant attention heads for protein family properties extracted from a single protein sequence. This analysis led us to propose an algorithm called PA_SPP for profile prediction using only a single protein sequence as input. In our algorithm, we apply the masked language modeling method of ProtAlbert. The results display high similarity between the predicted profiles and HSSP profiles.


title: Protein sequence profile prediction using ProtAlbert transformer


title: Biclique extension as an effective approach to predict novel interaction partners in metabolic compound-protein interaction networks.

abstract:
MotivationMetabolic networks are complex systems of connected chemical reactions with physical interactions between metabolites and proteins playing a critical role for both metabolic conversion and regulation. In this study, we aimed to predict previously unknown compound-protein interactions (CPI) by transferring the concept of biclique extension, which was developed in the context of drug-target interaction prediction and that is based on the rationale that interactions that readily extend an existing biclique are real, to metabolic CPI networks.

ResultsWe developed and tested a workflow to predict CPIs based on the concept of extending existing bicliques and applied it to E. coli and human using their respective known CPI network as input. Depending on the chosen biclique size, for the E. coli network we reached a sensitivity of 39% with an associated precision of 59%. For the larger human CPI network, a sensitivity of 78% with a false-positive rate of less than 5% and an associated precision of 75% was obtained. At more stringent settings, a precision as high as 95% was attainable at the expense of a lowered recall. Prediction performance significantly exceeded that obtained using randomized networks as input. Predicted novel interactions were tested for biomolecular function involvement, with TCA-cycle and ribosomal processes found associated with particularly pronounced statistical significance. As we demonstrate, our approach holds great potential to increase efficiency of experimental testing of CPIs and can readily be transferred to other species of interest.

Availability and implementationThe R code and datasets are available at https://github.com/SandraThieme/BiPredict.


title: IgMAT: immunoglobulin sequence multi-species annotation tool for any species including those with incomplete antibody annotation or unusual characteristics.

abstract:
The advent and continual improvement of high-throughput sequencing technologies has made immunoglobulin repertoire sequencing accessible and informative regardless of study species. However, to fully map changes in polyclonal dynamics, precise annotation of these constantly rearranging genes is pivotal. For this reason, data agnostic tools able to learn from presented data are required. Most sequence annotation tools are designed primarily for use with human and mouse antibody sequences which use databases with fixed species lists, applying very specific assumptions which select against unique structural characteristics. We present IgMAT, which utilises a reduced amino acid alphabet, incorporates multiple HMM alignments into a single consensus and enables the incorporation of user defined databases to better represent their species of interest.

Availability and implementationIgMAT has been developed as a python module, and is available on GitHub (https://github.com/TPI-Immunogenetics/igmat) for download under GPLv3 license.

Supplementary informationModel Breakdowns


title: Nonparametric Interrogation of Transcriptional Regulation in Single-Cell RNA and Chromatin Accessibility Multiomic Data

abstract: Epigenetic control of gene expression is highly cell-type- and context-specific. Yet, despite its complexity, gene regulatory logic can be broken down into modular components consisting of a transcription factor (TF) activating or repressing the expression of a target gene through its binding to a cis-regulatory region. Recent advances in joint profiling of transcription and chromatin accessibility with single-cell resolution offer unprecedented opportunities to interrogate such regulatory logic. Here, we propose a nonparametric approach, TRIPOD, to detect and characterize three-way relationships between a TF, its target gene, and the accessibility of the TFs binding site, using single-cell RNA and ATAC multiomic data. We apply TRIPOD to interrogate cell-type-specific regulatory logic in peripheral blood mononuclear cells and contrast our results to detections from enhancer databases, cis-eQTL studies, ChIP-seq experiments, and TF knockdown/knockout studies. We then apply TRIPOD to mouse embryonic brain data during neurogenesis and gliogenesis and identified known and novel putative regulatory relationships, validated by ChIP-seq and PLAC-seq. Finally, we demonstrate TRIPOD on SHARE-seq data of differentiating mouse hair follicle cells and identify lineage-specific regulation supported by histone marks for gene activation and super-enhancer annotations.


title: Deriving time-concordant event cascades from gene expression data: A case study for Drug-Induced Liver Injury (DILI)

abstract:
Adverse event pathogenesis is often a complex process which compromises multiple events ranging from the molecular to the phenotypic level. In toxicology, Adverse Outcome Pathways (AOPs) aim to formalize this as temporal sequences of events, in which event relationships should be supported by causal evidence according to the tailored Bradford-Hill criteria. One of the criteria is whether events are consistently observed in a certain temporal order and, in this work, we study this time concordance using the concept of "first activation" as data-driven means to generate hypotheses on potentially causal mechanisms. As a case study, we analysed liver data from repeat- dose studies in rats from the TG-GATEs database which comprises measurements across eight timepoints, ranging from 3 hours to 4 weeks post-treatment. We identified time concordant gene expression-derived events preceding adverse histopathology, which serves as surrogate readout for Drug-Induced Liver Injury (DILI). We find known mechanisms in DILI to be time-concordant, and show further that significance, frequency and log fold change (logFC) of differential expression are metrics which can additionally prioritize events although not necessary to be mechanistically relevant. Moreover, we used the temporal order of transcription factor (TF) expression and regulon activity to identify transcriptionally regulated TFs and subsequently combined this with prior knowledge on functional interactions to derive detailed gene-regulatory mechanisms, such as reduced Hnf4a activity leading to decreased expression and activity of Cebpa. At the same time, also potentially novel events are identified such as Sox13 which is highly significantly time- concordant and shows sustained activation over time. Overall, we demonstrate how time-resolved transcriptomics can derive and support mechanistic hypotheses by quantifying time concordance and how this can be combined with prior causal knowledge, with the aim of both understanding mechanisms of toxicity, as well as potential applications to the AOP framework. We make our results available in the form of a Shiny app (https://github.com/anikaliu/DILICascades_App), which allows users to query events of interest in more detail.

Author SummaryUnderstanding mechanisms from systems-scale biological data is of great relevance in toxicology as well as drug discovery; however how to generate causal hypotheses instead of correlations is by no means clear. In this work, we study the conserved temporal order of events and present an automatable framework to quantify and characterize time concordance across a large set of time- series. We apply this concept to events derived from time-resolved gene expression and histopathology from the TG-GATEs in vivo liver data as a case study. We were able to recover known events involved in the pathogenesis of Drug-Induced Liver Injury (DILI), and identify potentially novel pathway and transcription factors (TFs) which precede adverse histopathology. As complementary sources of evidence for causality, we additionally show how time concordance and prior knowledge on plausible interactions between TFs can be combined to derive causal hypotheses on the TFs mode of regulation and interaction partners. Overall, the results derived in our case study can serve as valuable hypothesis-free starting points for the development of Adverse Outcome Pathways for DILI, and demonstrate that our approach provides a novel angle to prioritize mechanistically relevant events.


title: Literature-Based Discovery beyond the ABC paradigm: a contrastive approach

abstract:
Literature-Based Discovery (LBD) aims to help researchers to identify relations between concepts which are worthy of further investigation by text-mining the biomedical literature. The vast majority of the LBD research follows the ABC model: a relation (A,C) is a candidate for discovery if there is some intermediate concept B which is related to both A and C. The ABC model has been successful in applications where the search space is strongly constrained, but there is limited evidence about its usefulness when applied in a broader context.

Through a case study of 8 recent discoveries related to neurodegenerative diseases (NDs), we show the limitations of the ABC model in an open-ended context. The study emphasizes the impact of the choice of source data and extraction method on the resulting knowledge base: different "views" of the biomedical literature offer different levels of accuracy and coverage. We propose a novel contrastive approach which leverages these differences between "views" in order to target relations between concepts of interest. We explore various parameters and demonstrate the relevance of our approach through quantitative evaluation on the 8 target discoveries.

The source data used in this article are publicly available. The different parts of the software used to process the data are published under open-source license and provided with detailed instructions. The main code for this paper is available at https://github.com/erwanm/lbd-contrast (required dependencies are detailed in the documentation). A prototype of the system is also provided as an online exploration tool at brainmend.adaptcentre.ie.


title: Revisiting the Use of Structural Similarity Index in Hi-C

abstract: Identification of dynamic changes in chromatin conformation is a fundamental task in genetics. In 2020, Galan et al.1 presented CHESS (Comparison of Hi-C Experiments using Structural Similarity), a novel computational algorithm designed for systematic identification of structural differences in chromatin-contact maps. Using CHESS, the same group recently reported that chromatin organization is largely maintained across tissues during dorsoventral patterning of fruit fly embryos despite tissue-specific chromatin states and gene expression2. However, here we show that the primary outputs of CHESS-namely, the structural similarity index (SSIM) profiles-are nearly identical regardless of the input matrices, even when query and reference reads were shuffled to destroy any significant differences. This issue stems from the dominance of the regional counting noise arising from stochastic sampling in chromatin-contact maps, reflecting a fundamentally incorrect assumption of the CHESS algorithm. Therefore, biological interpretation of SSIM profiles generated by CHESS requires considerable caution.


title: pChem: a modification-centric assessment tool for the performance of chemoproteomic probes

abstract: Chemoproteomics has emerged as a key technology to expand the functional space in complex proteomes for probing fundamental biology and for discovering new small molecule-based therapies. Here we report a modification-centric computational tool termed pChem to provide a streamlined pipeline for unbiased performance assessment of chemoproteomic probes. The pipeline starts with an experimental setting for isotopically coding probe-derived modifications (PDMs) that can be automatically recognized by pChem, with masses accurately calculated and sites precisely localized. Further, pChem exports on-demand reports by scoring the profiling efficiency, modification-homogeneity and proteome-wide residue selectivity of a tested probe. The performance and robustness of pChem were benchmarked by applying it to eighteen bioorthogonal probes. Of note, the analyses reveal that the formation of unexpected PDMs can be driven by endogenous reactive metabolites (e.g., bioactive aldehydes and glutathione). Together, pChem is a powerful and user-friendly tool that aims to facilitate the development of probes for the ever-growing field of chemoproteomics.