title: Behavioral context shapes vocal sequences in two anuran species with different repertoire sizes

abstract:
Acoustic signals in animals serve to convey context-dependent information to receivers. Birds and mammals combine diverse sounds into complex sequences to communicate, but these sequences largely remain understudied in other taxa. Anuran vocalizations are a prominent feature of their life history, and function in defense of territories and to attract mates. However, despite the spectacular diversity of anurans in tropical regions of the world, vocal diversity and communication strategies remain relatively poorly studied. Specifically, studies of vocal sequences and context-dependent vocal patterns in frogs remain few. Here, we investigated the context-dependent vocal repertoire and the use of vocal sequences by two anuran species belonging to different lineages, both endemic to the hyper-diverse Western Ghats of India. By recording vocal sequences both when frogs were alone and in the presence of a territorial rival, we present evidence that both species modify their vocal repertoire according to context. Specifically, one species appends notes to generate more complex sequences, whereas the other shifts to different note types, resulting in different sequences for different contexts. Thus, despite differences in repertoire size, both frog species are capable of adjusting the temporal sequence of vocalizations to communicate in different contexts. This study highlights the need for further studies of insular frogs, to understand how diversification across these continental islands has influenced the evolution of vocal repertoires, vocal sequence patterns and communication systems.

Lay SummaryAnimals employ complex sequences of acoustic signals to communicate in diverse behavioral contexts. Here, we demonstrate that two frog species with different vocal repertoires both modify the sequence of note emissions in the presence of a territorial rival. These patterns demonstrate that anurans are capable of complex shifts in the patterns of their vocalization, to communicate different messages to different receivers. Our findings demonstrate the value of studying behavioral diversity in tropical regions.


title: No evidence for a special role of language in feature-based categorization

abstract: The relationship between language and human thought is the subject of long-standing debate. One specific claim implicates language in feature-based categorization. According to this view, language resources facilitate object categorization based on a certain feature (e.g., color). Specifically, it is hypothesized that verbal labels help maintain focus on a relevant categorization criterion and reduce interference from irrelevant features. As a result, language impairment is expected to affect categorization of items grouped according to a single feature (low-dimensional categories, e.g.,  Things that are yellow), where many irrelevant features need to be inhibited, more than categorization of items that share many features (high-dimensional categories, e.g.,  Animals), where few irrelevant features need to be inhibited. In two behavioral studies with individuals with aphasia, we failed to find consistent support for the role of language in low-dimensional categorization. We also collected fMRI data from healthy adults and observed little activity in language-responsive brain regions during both low-dimensional and high-dimensional categorization. Combined, these results demonstrate that the language system is not implicated in object categorization. Our work adds to the growing evidence that, although language may assist in accessing task-relevant information (e.g., instructions), many cognitive tasks in adult brains proceed without recruiting the language system.

category: neuroscience


title: Does the Zebra Finch Mating Song Circuit Use Spike Times Efficiently?

abstract: Precise and reliable spike times are thought to subserve multiple possible functions, including improving the accuracy of encoding stimuli or behaviours relative to other coding schemes. Indeed, repeating sequences of spikes with sub-millisecond precision exist in nature, such as the synfire chain of spikes in area HVC of the zebra-finch mating-song circuit. Here, we analyzed what impact precise and reliable spikes have on the encoding accuracy for both the zebra-finch and more generic neural circuits using computational modelling. Our results show that neural circuits can use precisely timed spikes to encode signals with a higher-order accuracy than a conventional rate code. Circuits with precisely timed and reliably emitted spikes increase their encoding accuracy linearly with network size, which is the hallmark signature of an efficient code. This qualitatively differs from circuits that employ a rate code which increase their encoding accuracy with the square-root of network size. However, this improved scaling is dependent on the spikes becoming more accurate and more reliable with larger networks. Finally, we discuss how to test this scaling relationship in the zebra mating song circuit using both neural data and song-spectrogram-based recordings while taking advantage of the natural fluctuation in HVC network size due to neurogenesis. The zebra-finch mating-song circuit may represent the most likely candidate system for the use of spike-timing-based, efficient coding strategies in nature.


title: HIV-1 requires capsid remodelling at the nuclear pore for nuclear entry and integration

abstract:
The capsid (CA) lattice of the HIV-1 core plays a key role during infection. From the moment the core is released into the cytoplasm, it interacts with a range of cellular factors that, ultimately, direct the pre-integration complex to the integration site. For integration to occur, the CA lattice must disassemble. Early uncoating or a failure to do so has detrimental effects on virus infectivity, indicating that an optimal stability of the viral core is crucial for infection. Here, we introduced cysteine residues into HIV-1 CA in order to induce disulphide bond formation and engineer hyper-stable mutants that are slower or unable to uncoat, and then followed their replication. From a panel of mutants, we identified three with increased capsid stability in cells and found that, whilst the M68C/E212C mutant had a 5-fold reduction in reverse transcription, two mutants, A14C/E45C and E180C, were able to reverse transcribe to approximately WT levels. Moreover, these mutants only had a 5-fold reduction in 2-LTR circle production, suggesting that not only could reverse transcription complete in hyper-stable cores, but that the nascent viral cDNA could enter the nuclear compartment. Furthermore, we observed significant levels of A14C/E45C mutant capsid in nuclear and chromatin-associated fractions implying that the hyper-stable cores themselves entered the nucleus. Immunofluorescence studies revealed that although the A14C/E45C mutant capsid reached the nuclear pore with the same kinetics as wild type capsid, it was then retained at the pore in association with Nup153. Crucially, infection with the hyper-stable mutants did not promote CPSF6 re-localisation to nuclear speckles, despite the mutant capsids being competent for CPSF6 binding. These observations suggest that hyper-stable cores are not able to uncoat, or remodel, enough to pass through or dissociate from the nuclear pore and integrate successfully. This, is turn, highlights the importance of capsid lattice flexibility for nuclear entry. In conclusion, we hypothesise that during a productive infection, a capsid remodelling step takes place at the nuclear pore that releases the core complex from Nup153, and relays it to CPSF6, which then localises it to chromatin ready for integration.

AUTHOR SUMMARYThe mature viral core of human immunodeficiency virus (HIV) consists of a highly organised lattice formed by capsid molecules that encloses the viral RNA and viral enzymes. This lattice is crucial during the early stages of viral replication, as it has to break down - uncoat - at the right time and place in order for the viral DNA to integrate successfully. Lentiviruses, like HIV, can infect non-dividing cells and are able to access the host cell DNA by entering the nucleus through nuclear pores. Until recently, uncoating was thought to occur in the cytoplasm as the whole core was thought too large to pass through the nuclear pore. However, lately it has been suggested that uncoating might occur at the nuclear pore or even inside the nucleus and the site of uncoating is currently hotly debated. By investigating HIV mutants with an increased lattice stability, we have shown that lattice flexibility is crucial for nuclear entry. Provocatively, we observed hyper-stable mutant capsid in nuclear and chromatin-associated fractions suggesting that uncoating is not required for nuclear entry. Nonetheless, microscopy experiments suggested that these hyper-stable mutants were retained on the inner side of the nuclear pore, and were impaired for downstream events in the nucleus, leading to a severe infectivity defect. Therefore, we believe that an essential uncoating, or capsid lattice remodelling event normally takes place at the nuclear pore.

category: microbiology


title: Ciliary Hedgehog signaling regulates cell survival to build the facial midline

abstract: Craniofacial defects are among the most common phenotypes caused by ciliopathies, yet the developmental and molecular etiology of these defects is poorly understood. We investigated multiple mouse models of human ciliopathies (including Tctn2, Cc2d2a and Tmem231 mutants) and discovered that each displays hypotelorism, a narrowing of the midface. As early in development as the end of gastrulation, Tctn2 mutants displayed reduced activation of the Hedgehog (HH) pathway in the prechordal plate, the head organizer. This prechordal plate defect preceded a reduction of HH pathway activation and Shh expression in the adjacent neurectoderm. Concomitant with the reduction of HH pathway activity, Tctn2 mutants exhibited increased cell death in the neurectoderm and facial ectoderm, culminating in a collapse of the facial midline. Enhancing HH signaling by decreasing the gene dosage of a negative regulator of the pathway, Ptch1, decreased cell death and rescued the midface defect in both Tctn2 and Cc2d2a mutants. These results reveal that ciliary HH signaling mediates communication between the prechordal plate and the neurectoderm to provide cellular survival cues essential for development of the facial midline.

category: developmental biology


title: The archaeal triphosphate tunnel metalloenzyme SaTTM defines structural determinants for the diverse activities in the CYTH protein family

abstract:
CYTH is a large protein superfamily that is conserved in all three domains of life with its unique triphosphate tunnel metalloenzyme (TTM) fold. Besides phosphatase functions, e.g. as RNA triphosphatase, inorganic polyphosphatase or thiamine triphosphatase, some CYTH orthologs cyclize nucleotide triphosphates to 3,5-cyclic nucleotides. So far, archaeal CYTH proteins are annotated as adenylyl cyclases although experimental evidence is lacking. To address this gap, we characterized a CYTH ortholog, SaTTM, from the crenarchaeote Sulfolobus acidocaldarius. Our initial in silico studies suggested a close relationship between archaeal CYTH enzymes and class IV adenylyl cyclases compared to the other CYTH-subclasses, but biochemical data showed no cyclic nucleotide production. Instead, our structural and functional analyses show a classical TTM behavior. The Ca2+-inhibited Michaelis complex indicates a two-metal ion reaction mechanism analogous to other TTMs. Different co-crystal structures of SaTTM further reveal conformational dynamics in SaTTM, let us to assume feedback inhibition in TTMs due to tunnel closure in the product state. Combining our structural insights with sequence-similarity network based in silico analysis, we further set out a firm molecular basis for distinguishing CYTH orthologs with phosphatase activities from class IV adenylyl cyclases.

Major highlights- CyaB-like class IV adenylyl cyclase homologs in archaea are triphosphatases.
- The co-crystal structure of SaTTM in sulfate and triphosphate bound state revealed conformational transition of the TTM tunnel during catalysis.
- Atomic insights into TTM inhibition by calcium and pyrophosphate.
- In silico and structure-function analysis revealed the molecular determinant for functional diversification among CYTH proteins.

category: biochemistry


title: Elimination of GPI2 suppresses glycosylphosphatidylinositol GlcNAc transferase activity and alters GPI glycan modification in Trypanosoma brucei

abstract: The biosynthesis of glycosylphosphatidylinositol (GPI) membrane protein anchors is initiated in the endoplasmic reticulum by transfer of GlcNAc from the sugar nucleotide UDP-GlcNAc to phosphatidylinositol. The reaction is catalyzed by GPI GlcNAc transferase, a multi-subunit complex comprising the catalytic subunit Gpi3/PIG-A, as well as at least five other subunits including the hydrophobic protein Gpi2 which is essential for activity in yeast and mammals, but whose function is not known. Here we exploited Trypanosoma brucei (Tb), an early diverging eukaryote and important model organism, to investigate the function of Gpi2. We generated trypanosomes that lack TbGPI2 and found that in TbGPI2-null parasites (i) GPI GlcNAc transferase activity is reduced but not lost, in contrast with the situation in yeast and human cells, (ii) the GPI GlcNAc transferase complex persists, but its architecture is affected, with loss of at least the TbGPI1 subunit, and (iii) the GPI anchors of the major surface proteins are underglycosylated when compared with their wild-type counterparts, indicating the importance of TbGPI2 for reactions that are expected to occur in the Golgi apparatus. Additionally, TbGPI2-null parasites were unable to perform social motility, a form of collective migration on agarose plates. Immunofluorescence microscopy localized TbGPI2 to the endoplasmic reticulum as expected, but also to the Golgi apparatus, suggesting that in addition to its expected function as a subunit of the GPI GlcNAc transferase complex, TbGPI2 may have an enigmatic non-canonical role in Golgi-localized GPI anchor modification in trypanosomes.


title: Amelioration of hemophilia B through CRISPR/Cas9 induced homology-independent targeted integration

abstract: Site-specific integration of exogenous gene through genome editing is a promising strategy for gene therapy. However, homology-directed repair (HDR) only occurring in proliferating cells is inefficient especially in vivo. To investigate the efficacy of Cas9-induced homology-independent targeted integration (HITI) strategy for gene therapy, a rat hemophilia B model was generated and employed. Through HITI, a DNA sequence encoding the last exon of rat Albumin (rAlb) gene fused with a high-specific-activity Factor IX variant (R338L) using T2A, was inserted into the last intron of rAlb via recombinant adeno-associated viral (rAAV). The knock-in efficiency reached up to 3.66% determined by ddPCR. The clotting time was reduced to normal level 4 weeks after treatment, and the circulating FIX level was gradually increased up to 52% of normal during 9 months even after partial hepatectomy, demonstrating the amelioration of hemophilia. Through PEM-seq, no significant off-targeting effect was detected. Moreover, this study provides a promising therapeutic approach for hereditary diseases.

category: bioengineering


title: Cell cycle gene regulation dynamics revealed by RNA velocity and deep-learning

abstract: The cell cycle is a fundamental process of life, however, a quantitative understanding of gene regulation dynamics in the context of the cell cycle is still far from complete. Single-cell RNA-sequencing (scRNA-seq) technology gives access to its dynamics without externally perturbing the cell. Here, we build a high-resolution map of the cell cycle transcriptome based on scRNA-seq and deep-learning. By generating scRNA-seq libraries with high depth, in mouse embryonic stem cells and human fibroblasts, we are able to observe cycling patterns in the unspliced-spliced RNA space for single genes. Since existing methods in scRNA-seq are not efficient to measure cycling gene dynamics, we propose a deep learning approach to fit these cycling patterns sorting single cells across the cell cycle. We characterize the cell cycle in asynchronous pluripotent and differentiated cells identifying major waves of transcription during the G1 phase and systematically study the G1-G0 transition where the cells exit the cycle. Our work presents to the scientific community a broader understanding of RNA velocity and cell cycle maps, that we applied to pluripotency and differentiation. Our approach will facilitate the study of the cell cycle in multiple cellular models and different biological contexts, such as cancer and development.

category: systems biology


title: Longitudinal allometry of sulcal morphology in health and schizophrenia

abstract:
Scaling between subcomponents of cortical folding and total brain volume (TBV) in healthy individuals (HI) is allometric, i.e. non-linear. It is unclear whether this is also true in individuals with schizophrenia (SZ) or first-episode psychosis (FEP). The current study first confirmed normative allometric scaling norms in HI using discovery and replication samples. Cross-sectional and longitudinal diagnostic differences in folding subcomponents were then assessed using an allometric analytic framework.

Structural imaging from a longitudinal (sample 1: HI and SZ, nHI Baseline = 298, nSZ Baseline = 169, nHI Follow-up = 293, nSZ Follow-up = 168, a total of 1087 images, all individuals [≥] 2 images, age 16-69 years) and a cross-sectional sample (sample 2: nHI = 61 and nFEP = 89, age 10-30 years) is leveraged to calculate global folding and its nested subcomponents: sulcation index (SI, total sulcal/cortical hull area) and determinants of sulcal area; sulcal length and sulcal depth.

Scaling of the SI, sulcal area, and sulcal length with TBV in SZ and FEP was allometric and did not differ from HI. Longitudinal age trajectories demonstrated steeper loss of SI and sulcal area through adulthood in SZ. Longitudinal allometric analysis revealed that both annual change in SI and sulcal area was significantly stronger related to change in TBV in SZ compared to HI.

Our results detail the first evidence of the disproportionate contribution of changes in SI and sulcal area to TBV changes in SZ. Longitudinal allometric analysis of sulcal morphology provides deeper insight into lifespan trajectories of cortical folding in SZ.


title: NK cell metabolic adaptation to infection promotes survival and viral clearance.

abstract: Natural killer (NK) cells are essential for early protection against virus infection, and must metabolically adapt to the energy demands of activation. Here, we found upregulation of the metabolic adaptor hypoxia inducible factor-1 (HIF-1) is a feature of NK cells during murine cytomegalovirus (MCMV) infection in vivo. HIF-1-deficient NK cells failed to control viral load, causing increased morbidity. No defects were found in effector functions of HIF-1 KO NK cells however, their numbers were significantly reduced. Loss of HIF-1 did not affect NK cell proliferation during in vivo infection and in vitro cytokine stimulation. Instead, we found HIF1-deficient NK cells showed increased expression of the pro-apoptotic protein Bim and glucose metabolism was impaired during cytokine stimulation in vitro. Similarly, during MCMV infection HIF1-deficient NK cells upregulated Bim and had increased caspase activity. Thus, NK cells require HIF-1-dependent metabolic functions to repress Bim expression and sustain cell numbers for an optimal virus response.

category: immunology


title: Plasma Membrane Calcium ATPase (PMCA) Regulates Stoichiometry of CD4+ T-cell Compartments

abstract: Immune responses involve mobilization of T cells within naive and memory compartments. Tightly regulated Ca2+ levels are essential for balanced immune outcomes. How Ca2+ contributes to regulating compartment stoichiometry is unknown. Here, we show that plasma membrane Ca2+ ATPase 4 (PMCA4) is differentially expressed in human CD4+ T compartments yielding distinct store operated Ca2+ entry (SOCE) profiles. Modulation of PMCA4 yielded a more prominent increase of SOCE in memory than in naive CD4+ T cell. Interestingly, downregulation of PMCA4 reduced the effector compartment fraction and led to accumulation of cells in the naive compartment. In silico analysis and chromatin immunoprecipitation unraveled Ying Yang 1 (YY1) as a transcription factor regulating PMCA4 expression. Analyses of PMCA and YY1 expression patterns following activation and of PMCA promoter activity following downregulation of YY1 highlight repressive role of YY1 on PMCA expression. Our findings show that under the transcriptional control by YY1, PMCA4 adapts Ca2+ levels to cellular requirements during effector and quiescent phases and thereby represent a potential target to intervene with the outcome of the immune response.


title: A novel Mycobacterium tuberculosis-specific subunit vaccine provides synergistic immunity upon co-administration with Bacillus Calmette-Guerin

abstract: Tuberculosis (TB) remains a global health crisis. Following encouraging clinical results of subunit vaccination and revaccination with Bacillus Calmette-Guerin (BCG), it has been suggested to combine BCG and subunit vaccines for increased efficacy. Current subunit vaccines are almost exclusively designed as BCG boosters. The goal of this study was to design a subunit vaccine that does not share antigens with BCG and explore the advantages of a BCG+subunit vaccine co-administration strategy, where the two vaccines do not cross-react. Eight protective antigens were selected to create a Mycobacterium tuberculosis (Mtb)-specific subunit vaccine, named H107. Whereas subunit vaccines with BCG-shared antigens displayed cross-reactivity to BCG in vivo in both mice and humans, H107 showed no cross-reactivity and did not inhibit BCG colonization in mice. Encouragingly, co-administering H107 with BCG (BCG+H107) led to increased adaptive immune responses against both H107 and BCG leading to improved BCG-mediated immunity. In contrast to subunit vaccines with BCG-shared antigens,  boosting BCG with H107 led to substantial expansion of clonal diversity in the T cell repertoire, and BCG+H107 co-administration conferred significantly increased Th17 responses and less differentiated CD4 T cells. CD4 T cells induced by BCG+H107 maintained functional superiority after Mtb infection, and BCG+H107 provided significantly increased long-term protection compared to both BCG and H107 alone, as well as, BCG co-administered with a subunit vaccine composed of antigens shared with BCG. Overall, we identify several advantages of combining an Mtb-specific subunit vaccine with BCG and introduce H107 as a BCG-complementing vaccine with distinctive value for co-administration with BCG.


title: Prediction of ligand-receptor pharmacological activities using a combined docking and machine learning approach

abstract: G protein coupled receptors (GPCRs) are valuable therapeutic targets for many diseases. A central question of GPCR drug discovery is to understand what determines the agonism or antagonism of ligands which bind them. Ligands exert their action via the interactions in the ligand binding pocket. We hypothesised that there is a common set of receptor interactions made by ligands of diverse structures that mediate their action and that among a large dataset of different ligands, the functionally important interactions will be over-represented. We computationally docked ~2700 known {beta}2AR ligands to multiple {beta}2AR structures, generating ca 75,000 docking poses and predicted all atomic interactions between the receptor and the ligand. We used machine learning (ML) techniques to identify specific interactions that correlate with the agonist or antagonist activity of these ligands. The interpretation of ML analysis in human understandable form allowed us to construct an exquisitely detailed structure-activity relationship that identifies small changes to the ligands that invert their activity and thus helps to guide the drug discovery process. This approach can be readily applied to any drug target.

category: bioinformatics


title: Combined docking and machine learning identifies key molecular determinants of ligand pharmacological activity on β2 adrenoceptor.


title: Hubness reduction improves clustering and trajectory inference in single-cell transcriptomic data

abstract:
BackgroundSingle-cell RNA-seq datasets are characterized by large ambient dimensionality, and their analyses can be affected by various manifestations of the dimensionality curse. One of these manifestations is the hubness phenomenon, i.e. existence of data points with surprisingly large incoming connectivity degree in the neighbourhood graph. Conventional approach to dampen the unwanted effects of high dimension consists in applying drastic dimensionality reduction. It remains unexplored if this step can be avoided thus retaining more information than contained in the low-dimensional projections, by correcting directly hubness.

ResultsWe investigate the phenomenon of hubness in scRNA-seq data in spaces of increasing dimensionality. We also link increased hubness to increased levels of dropout in sequencing data. We show that hub cells do not represent any visible technical or biological bias. The effect of various hubness reduction methods is investigated with respect to the visualization, clustering and trajectory inference tasks in scRNA-seq datasets. We show that hubness reduction generates neighbourhood graphs with properties more suitable for applying machine learning methods; and that it outperforms other state-of-the-art methods for improving neighbourhood graphs. As a consequence, clustering, trajectory inference and visualisation perform better, especially for datasets characterized by large intrinsic dimensionality.

ConclusionHubness is an important phenomenon in sequencing data. Reducing hubness can be beneficial for the analysis of scRNA-seq data with large intrinsic dimensionality in which case it can be an alternative to drastic dimensionality reduction.


title: Analyzing Assay Redundancy in Metabolomics using Unique Ion Signatures

abstract:
Targeted, untargeted and data-independent acquisition (DIA) metabolomics workflows are often hampered by ambiguous identification based on either MS1 information alone or relatively few MS2 fragment ions. While DIA methods have been popularized in proteomics, it is less clear whether they are suitable for metabolomics workflows due to their large precursor isolation windows and complex co-isolation patterns. Here, we quantitatively investigate the conditions necessary for unique metabolite detection in complex backgrounds using precursor and fragment ion mass-to-charge separation, comparing three benchmarked MS methods (MS1, MRM, DIA). Our simulations show that DIA outperformed MS1-only and MRM-based methods with regards to specificity by a factor of ~2.8-fold and ~1.8-fold, respectively. Additionally, we show that our results are not dependent on the number of transitions used or the complexity of the background matrix. Finally, we show that collision energy is an important factor in unambiguous detection and that a single collision energy setting per compound cannot achieve optimal pairwise differentiation of compounds. Our analysis demonstrates the power of using both high resolution precursor and high resolution fragment ion m/z for unambiguous compound detection. This work also establishes DIA as an emerging MS acquisition method with high selectivity for metabolomics, outperforming both DDA and MRM with regards to unique compound identification potential.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=132 SRC="FIGDIR/small/434579v2_ufig1.gif" ALT="Figure 1">
View larger version (23K):
org.highwire.dtl.DTLVardef@142b100org.highwire.dtl.DTLVardef@1728a1corg.highwire.dtl.DTLVardef@c953c8org.highwire.dtl.DTLVardef@19b856a_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Analyzing Assay Specificity in Metabolomics using Unique Ion Signature Simulations


title: Targeting de novo lipogenesis and the Lands cycle induces ferroptosis in KRAS-mutant lung cancer

abstract: Mutant KRAS (KM) is the most common oncogene in lung cancer (LC). KM regulates several metabolic networks, but their role in tumorigenesis is still not sufficiently characterized to be exploited in cancer therapy. To identify metabolic networks specifically deregulated in KMLC, we characterized the lipidome of genetically engineered LC mice, cell lines, patient derived xenografts and primary human samples. We also determined that KMLC, but not EGFR-mutant (EGFR-MUT) LC, is enriched in triacylglycerides (TAG) and phosphatidylcholines (PC). We also found that KM upregulates fatty acid synthase (FASN), a rate-limiting enzyme in fatty acid (FA) synthesis promoting the synthesis of palmitate and PC. We determined that FASN is specifically required for the viability of KMLC, but not of LC harboring EGFR-MUT or wild type KRAS. Functional experiments revealed that FASN inhibition leads to ferroptosis, a reactive oxygen species (ROS)-and iron-dependent cell death. Consistently, lipidomic analysis demonstrated that FASN inhibition in KMLC leads to accumulation of PC with polyunsaturated FA (PUFA) chains, which are the substrate of ferroptosis. Integrating lipidomic, transcriptome and functional analyses, we demonstrated that FASN provides saturated (SFA) and monounsaturated FA (MUFA) that feed the Lands cycle, the main process remodeling oxidized phospholipids (PL), such as PC. Accordingly, either inhibition of FASN or suppression of the Lands cycle enzymes PLA2 and LPCAT3, promotes the intracellular accumulation of lipid peroxides and ferroptosis in KMLC both in vitro and in vivo. Our work supports a model whereby the high oxidative stress caused by KM dictates a dependency on newly synthesized FA to repair oxidated phospholipids, establishing a targetable vulnerability. These results connect KM oncogenic signaling, FASN induction and ferroptosis, indicating that FASN inhibitors already in clinical trial in KMLC patients (NCT03808558) may be rapidly deployed as therapy for KMLC.

category: cancer biology


title: Micro-region transcriptomics of fixed human tissue using Pick-Seq

abstract: Spatial transcriptomics and multiplexed imaging are complementary methods for studying tissue biology and disease. Recently developed spatial transcriptomic methods use fresh-frozen specimens but most diagnostic specimens, clinical trials, and tissue archives rely on formaldehyde-fixed tissue. Here we describe the Pick-Seq method for deep spatial transcriptional profiling of fixed tissue. Pick-Seq is a form of micro-region sequencing in which small regions of tissue, containing 5-20 cells, are mechanically isolated on a microscope and then sequenced. We demonstrate the use of Pick-Seq with several different fixed and frozen human specimens. Application of Pick-Seq to a human melanoma with complex histology reveals significant differences in transcriptional programs associated with tumor invasion, proliferation, and immuno-editing. Parallel imaging confirms changes in immuno-phenotypes and cancer cell states. This work demonstrates the ability of Pick-Seq to generate deep spatial transcriptomic data from fixed and archival tissue with multiplexed imaging in parallel.


title: Antagonistic effects of chemical mixtures on the oxidative stress response are silenced by heat stress and reversed under dietary restriction

abstract: Chemical agents released into the environment can induce oxidative stress in organisms, which is detrimental for health. Although environmental exposures typically incorporate multiple chemicals, organismal studies on oxidative stress derived from chemical agents commonly study exposures to individual compounds. In this work, we explore how chemical mixtures drive the oxidative stress response under various conditions in the nematode C. elegans, by quantitatively assessing levels of gst-4 expression. Our results indicate that naphthoquinone mixtures drive responses differently than individual components, and that altering environmental conditions, such as increased heat and reduced food availability, result in dramatically different oxidative stress responses mounted by C. elegans. When exposed to heat, the oxidative stress response is diminished. Notably, when exposed to limited food, the oxidative stress response specific to juglone is significantly heightened, while identified antagonistic interactions between some naphthoquinone components in mixtures are abolished. This implies that organismal responses to xenobiotics is confounded by environment and stressor interactions. Given the high number of variables under study, and their potential combinations, a simplex centroid design was used to capture such non-trivial response over the design space. This makes the case for the adoption of Design of Experiments approaches as they can greatly expand the experimental space probed in noisy biological readouts, and in combinatorial experiments. Our results also reveal gaps in our current knowledge of the organismal oxidative stress response, which can be addressed by employing sophisticated design of experiments approaches to identify significant interactions.

category: physiology


title: Electrophysiological signatures of brain aging in autism spectrum disorder.

abstract: Emerging evidence suggests that aging processes may be altered in adults with autism spectrum disorder (ASD). However, it remains unclear if oscillatory slowing, a key neurophysiological change in the aging brain, manifests atypically in this population. This study sought to examine patterns of age-related oscillatory slowing in adults with ASD, captured by reductions in the brains peak alpha frequency. Resting-state EEG data from adults (18-70 years) with ASD (N=93) and age-matched neurotypical (NT) controls (N=87) were pooled from three independent datasets. A robust curve-fitting procedure quantified the peak frequency of alpha oscillations (7-13Hz) across all brain regions. Associations between peak alpha frequency and age were assessed and compared between groups. Consistent with characteristic patterns of oscillatory slowing, peak alpha frequency was negatively associated with age across the entire sample (p<.0001). A significant group by age interaction revealed that this relationship was more pronounced in adults with ASD (p<.01), suggesting that that age-related oscillatory slowing may be accelerated in this population. Scalable EEG measures such as peak alpha frequency could provide insights into neural aging that are crucially needed to inform care plans and preventive interventions that can promote successful aging in ASD.


title: An RNase III processed, antisense RNA pair regulates a Campylobacter jejuni colonization factor

abstract: Bacterial small RNAs (sRNAs) are important post-transcriptional regulators in stress responses and virulence. They can be derived from an expanding list of genomic contexts, such as processing from parental transcripts by RNase E. The role of RNase III in sRNA biogenesis is less well understood despite its well-known roles in rRNA processing, RNA decay, and cleavage of sRNA-mRNA duplexes. Here, we show that RNase III processes a pair of cis-encoded sRNAs (CJnc190 and CJnc180) of the foodborne pathogen Campylobacter jejuni. While CJnc180 processing by RNase III requires CJnc190, In contrast, RNase III processes CJnc190 independent of CJnc180 via cleavage of an intramolecular duplex. We also show that CJnc190 directly represses translation of the colonization factor PtmG by targeting a G-rich ribosome binding site, and uncover that CJnc180 is a cis-acting antagonist of CJnc190, indirectly affecting ptmG regulation. Our study highlights a role for RNase III in sRNA biogenesis and adds cis-encoded RNAs to the expanding diversity of transcripts that antagonize bacterial sRNAs.


title: RNase III-mediated processing of a trans-acting bacterial sRNA and its cis-encoded antagonist


title: Thermal stress triggers productive viral infection of a key coral reef symbiont

abstract: Climate change-driven ocean warming is increasing the frequency and severity of bleaching events, in which corals appear whitened after losing their dinoflagellate endosymbionts (family Symbiodiniaceae). Viral infections of Symbiodiniaceae may contribute to some bleaching signs, but little empirical evidence exists to support this hypothesis. We present the first temporal analysis of a lineage of Symbiodiniaceae-infecting positive-sense single-stranded RNA viruses ( dinoRNAVs) in coral colonies, which were exposed to a 5-day heat treatment (+2.1{degrees}C). A total of 124 dinoRNAV major capsid protein gene  aminotypes (unique amino acid sequences) were detected from five colonies of two closely related Pocillopora-Cladocopium (coral-symbiont) combinations in the experiment; most dinoRNAV aminotypes were shared between the two coral-symbiont combinations (64%) and among multiple colonies (82%). Throughout the experiment, seventeen dinoRNAV aminotypes were found only in heat-treated fragments, and 22 aminotypes were detected at higher relative abundances in heat-treated fragments. DinoRNAVs in fragments of some colonies exhibited higher alpha diversity and dispersion under heat stress. Together, these findings provide the first empirical evidence that exposure to high temperatures triggers some dinoRNAVs to switch from a persistent to a productive infection mode within heat-stressed corals. Over extended time frames, we hypothesize that cumulative dinoRNAV production in the Pocillopora-Cladocopium system could affect colony symbiotic status, for example, by decreasing Symbiodiniaceae densities within corals. This study sets the stage for reef-scale investigations of dinoRNAV dynamics during bleaching events.

category: ecology


title: Structural insights into hormone recognition by the human glucose-dependent insulinotropic polypeptide receptor

abstract: Glucose-dependent insulinotropic polypeptide (GIP) is a peptide hormone that exerts crucial metabolic functions by binding and activating its cognate receptor, GIPR. As an important therapeutic target, GIPR has been subjected to intensive structural studies without success. Here, we report the cryo-EM structure of the human GIPR in complex with GIP and a Gs heterotrimer at a global resolution of 2.9 [A]. GIP adopts a single straight helix with its N terminus dipped into the receptor transmembrane domain (TMD), while the C-terminus is closely associated with the extracellular domain and extracellular loop 1. GIPR employs conserved residues in the lower half of the TMD pocket to recognize the common segments shared by GIP homologous peptides, while uses non-conserved residues in the upper half of the TMD pocket to interact with residues specific for GIP. These results provide a structural framework of hormone recognition and GIPR activation.

category: molecular biology


title: Cell signaling pathways in human mutant PAX6 corneal cells: an in vitro model for aniridia-related keratopathy

abstract:
PURPOSETo establish an in vitro model of aniridia-related keratopathy (ARK) using CRISPR/Cas9 engineered human keratocytes with mutations in the PAX6 gene, and to study the Notch Homolog 1, Translocation-Associated (Notch1), sonic hedgehog (SHH), mammalian target of rapamycin (mTOR), and Wnt/{beta}-catenin signaling pathways in the PAX6 mutant keratocytes.

METHODSPrimary human keratocytes were isolated from healthy corneas. Keratocytes were transduced with Cas9 lentiviral particles in order to create cells stably expressing Cas9 nuclease. Lentiviral particles carrying PAX6 sgRNA were transduced into the Cas9 keratocytes creating mutants. Analysis of signaling pathways was assessed by RT-qPCR for gene expression and western blot for protein expression.

RESULTSHuman keratocytes stably expressing Cas9 nuclease were created. Keratocytes carrying PAX6 gene mutation were successfully generated. PAX6 mutant keratocytes showed modified expression patterns of extracellular matrix components such as collagens and fibrotic markers. Analysis of the Notch1, SHH, mTOR, and Wnt/{beta}-catenin signaling pathways in the PAX6 mutant keratocytes revealed altered gene and protein expression of the key players involved in these pathways.

CONCLUSIONSA properly functioning PAX6 gene in keratocytes is crucial for the regulation of signaling pathways important for cell fate determination, proliferation, and inflammation. Pax6 mutation in the in vitro settings leads to changes in these pathways which resemble those found in corneas of patients with ARK.


title: Conformational specificity of opioid receptors is determined by subcellular location irrespective of agonist

abstract: The prevailing model for the variety in drug responses is that they stabilize distinct active states of their G protein-coupled receptor (GPCR) targets, allowing coupling to different effectors. However, whether the same ligand can produce different GPCR active states based on the environment of receptors in cells is a fundamental unanswered question. Here we address this question using live cell imaging of conformational biosensors that read out distinct active conformations of the {delta}-opioid receptor (DOR), a physiologically relevant GPCR localized to Golgi and the surface in neurons. We show that, although Golgi and surface pools of DOR regulated cAMP, the two pools engaged distinct conformational biosensors in response to the same ligand. Further, DOR recruited arrestin on the plasma membrane but not the Golgi. Our results suggest that the same agonist drives different conformations of a GPCR at different locations, allowing receptor coupling to distinct effectors at different locations.

category: cell biology


title: Copper dependent ERK1/2 phosphorylation is essential for the viability of neurons and not glia

abstract: Intracellular Copper [Cu(I)] has been hypothesized to play role in the differentiation of the neurons. This necessitates understanding the role of Cu(I) not only in the neurons but also in the glia considering their anatomical proximity, contribution towards ion homeostasis, neuronal physiology, and neurodegeneration. In this study we did a systematic investigation of the changes in the cellular copper homeostasis during neuronal and glial differentiation and the pathways triggered by them. Our study demonstrates increased mRNA for the plasma membrane copper transporter CTR1 leading to an increased Cu(I) during neuronal (PC-12) differentiation. ATP7A is retained in the Trans Golgi Network (TGN) despite high Cu(I) demonstrating its utilization in triggering the pathways towards the neuronal differentiation. One of these pathways is ERK1/2 phosphorylation accompanying the differentiation of both PC-12 and human fetal brain derived neuronal progenitor cells. The study demonstrates that the ERK1/2 phosphorylation is essential for the viability of the neurons. In contrast, differentiated C-6 (glia) cells contain low intracellular copper and significant downregulation of the ERK1/2 phosphorylation. Interestingly ATP7A shows vesicular localization despite the low copper in the glia. In addition to the TGN in the perinuclear region, ATP7A localizes into RAB11 positive recycling endosomes in the glial neurites, not observed in the neurons. Our study demonstrates role of the copper dependent ERK1/2 phosphorylation in the neuronal differentiation. Whereas glial differentiation largely involves sequestration of Cu(I) into the endosomes potentially (i) for ready release to the neurons (ii) rendering cytosolic copper unavailable for pathways like the ERK1/2 activation.


title: Optogenetic manipulation of YAP cellular localisation and function

abstract: YAP, an effector of the Hippo signalling pathway, promotes organ growth and regeneration. Prolonged YAP activation results in uncontrolled proliferation and cancer. Therefore, exogenous regulation of YAP activity has potential translational applications. We present a versatile optogenetic construct (optoYAP) for manipulating YAP localisation, and consequently its activity and function. We attached a LOV2 domain that photocages a nuclear localisation signal (NLS) to the N-terminus of YAP. In 488 nm light, the LOV2 domain unfolds, exposing the NLS, which shuttles optoYAP into the nucleus. Nuclear import of optoYAP is reversible and tuneable by light intensity. In cell culture, activated optoYAP promotes YAP target gene expression, cell proliferation, and anchorage-independent growth. Similarly, we can utilise optoYAP in zebrafish embryos to modulate target genes. OptoYAP is functional in both cell culture and in vivo, providing a powerful tool to address basic research questions and therapeutic applications in regeneration and disease.


title: CLEMSite, a software for automated phenotypic screens using light microscopy and FIB-SEM

abstract: Correlative light and electron microscopy (CLEM) combines two imaging modalities, balancing out the limits of one technique with the other. In recent years, Focused Ion Beam Scanning Electron Microscopy (FIB-SEM) has emerged as a flexible method that enables semi-automated volume acquisition at the ultrastructural level. We present a toolset for adherent cultured cells that enables tracking and finding cell regions previously identified in light microscopy, in the FIB-SEM along with automatic acquisition of high-resolution volume datasets. We detect a grid pattern in both modalities (LM and EM), which identifies common reference points. The novel combination of these techniques enables complete automation of the workflow. This includes setting the coincidence point of both ion and electron beams, automated evaluation of the image quality and constantly tracking the sample position with the microscopes field of view reducing or even eliminating operator supervision. We show the ability to target the regions of interest in EM within 5 {micro}m accuracy, while iterating between different targets and implementing unattended data acquisition. Our results demonstrate that executing high throughput volume acquisition in electron microscopy is possible.


title: Nonsense Mediated RNA Decay Is a Unique Vulnerability of Cancer Cells with SF3B1 and U2AF1 Mutations

abstract: Nonsense-mediated RNA decay (NMD) is well recognized as an RNA surveillance pathway that targets aberrant mRNAs with premature translation termination codons (PTCs) for degradation; however, its molecular mechanisms and roles in health and disease remain incompletely understood. In this study, we developed a novel reporter system that can accurately measure NMD activity in individual cells. By carrying out a genome-wide CRISPR/Cas9 knockout screen using this reporter system, we identified novel NMD-promoting factors, including multiple components of the SF3B complex and other U2 spliceosome factors. Interestingly, we also found that cells with mutations in the U2 spliceosome genes SF3B1 and U2AF1--which are commonly found in myelodysplastic syndrome (MDS) and cancers--have overall attenuated NMD activity. Furthermore, we found that compared to wild type cells, SF3B1 and U2AF1 mutant cells are more sensitive to NMD inhibition, a phenotype that is accompanied by elevated DNA replication obstruction, DNA damage and chromosomal instability. Remarkably, the sensitivity of spliceosome mutant cells to NMD inhibition could be rescued by overexpression of RNase H1, which removes R-loops in the genome. Together, our findings shed new light on the functional interplay between NMD and RNA splicing and suggest a novel strategy for the treatment of MDS and cancers with spliceosome mutations.


title: Deep Learning Shows Cellular Senescence Is a Barrier to Cancer Development

abstract: Cellular senescence is a critical component of aging and many age-related diseases, but understanding its role in human health is challenging in part due to the lack of exclusive or universal markers. Using neural networks, we achieve high accuracy in predicting senescence state and type from the nuclear morphology of DAPI-stained human fibroblasts, murine astrocytes and fibroblasts derived from premature aging diseases in vitro. After generalizing this approach, the predictor recognizes an increasing rate of senescent cells with age in H&E-stained murine liver tissue and human dermal biopsies. Evaluating corresponding medical records reveals that individuals with increased senescent cells have a significantly decreased rate of malignant neoplasms, lending support for the protective role of senescence in limiting cancer development. In sum, we introduce a novel predictor of cellular senescence and apply it to diagnostic medical images, indicating cancer occurs more frequently for those with a lower rate of senescence.


title: Mechano-transduction via the pectin-FERONIA complex regulates ROP6 GTPase signaling in Arabidopsis

abstract: During growth and morphogenesis, plant cells respond to mechanical stresses resulting from spatiotemporal changes in the cell wall that bear high internal turgor pressure. Microtubule (MT) arrays are re-organized to align in the direction of maximal tensile stress to guide the synthesis of cellulose, reinforcing the local cell wall. However, how mechanical forces regulate MT re-organization remains largely unknown. Here, we demonstrate that mechanical signaling that is based on the CrRLK1L subfamily receptor kinase FERONIA (FER) regulates the reorganization of cortical MT in cotyledon epidermal pavement cells (PC) in Arabidopsis. Recessive mutations in FER compromised MT response to mechanical perturbations such as single cell ablation, compression and Isoxaben treatment in these pavement cells. These perturbations promoted the activation of ROP6 GTPase that acts directly downstream of FER. Furthermore, defects in the ROP6 signaling pathway negated the reorganization of cortical MTs induced by these stresses. Finally, reduction in highly demethylesterified pectin, which binds the extracellular malectin domain of FER and is required for FER-mediated ROP6 activation, also impacted mechanical induction of cortical MT reorganization. Taken together our results suggest that the FER-pectin complex senses and/or transduce mechanical forces to regulate MT organization through activating the ROP6 signaling pathway in Arabidopsis.


title: Distinct Hypoxia-induced Translational Profiles of Embryonic and Adult-derived Macrophages

abstract: Tissue homeostasis and repair are orchestrated by resident and newly recruited macrophages that alter their gene expression program in response to changes in tissue microenvironment. Embryonic macrophages, such as fetal liver derived macrophages (FLDM) seed the organs, including heart and lung during embryonic development and persist throughout the adult lifetime, while bone marrow-derived macrophages (BMDM) are recruited following an acute perturbation. Transcriptome analyses of FLDM and BMDM identified differences between them at the level of RNA expression, which correlates imperfectly with protein levels. Post-transcriptional regulation by microRNAs (miRNAs) and RNA-binding proteins determines mRNA stability and translation rate and may override transcriptional cues in response to environmental changes, such as hypoxia. To identify distinct features of FLDM and BMDM response to hypoxia at the level of translation, we employed translating ribosome affinity purification (TRAP) to isolate polysomal RNA. RNA-seq profiling of translated RNA identified distinct hypoxia-induced translational signature of BMDM (Ly6e, vimentin and glycolysis-associated enzymes Pgk1, Tpi1, Aldoa, Ldha) and FLDM (chemokines Ccl7 and Ccl2). By translational profiling of BMDM and FLDM with deletion of the RNA-binding protein HuR, we identified transcripts that were dependent on HuR. These findings highlight the importance of HuR and identify its distinct targets for post-transcriptional regulation of gene expression in embryonic vs. adult-derived macrophages.


title: Intraflagellar transport trains can turn around without the ciliary tip complex

abstract: Cilia and flagella are microtubule doublet based organelles found across the eukaryotic tree of life. Their very high aspect ratio and crowded interior are unfavourable to diffusive transport for their assembly and maintenance. Instead, a highly dynamic system of intraflagellar transport (IFT) trains moves rapidly up and down the cilium. However, the mechanism of how these trains turn around upon reaching the ciliary tip has remained elusive. It has been hypothesized that there exists a dedicated calcium-dependent protein-based machinery at the ciliary tip to mediate this conversion. In this work, we use physical and chemical methods to manipulate IFT in the cilia of the unicellular green alga Chlamydomonas reinhardtii to show that no such stationary tip-machinery is required for IFT turnaround. Instead, we demonstrate that the conversion from anterograde to retrograde IFT trains is a calcium independent intrinsic ability of the IFT system.


title: High throughput genome-wide single cell protein:DNA binding site mapping by targeted insertion of promoters (TIP-seq)

abstract: Assessing cell to cell, and importantly, chromosome to chromosome, heterogeneity in cellular phenotypes is a central goal of modern cell biology. However, chromatin profiling in single cells has been extremely challenging, and single chromosome profiling has not been achieved. In cases where single cell methods have shown promise, success has been mainly limited to histone proteins and/or require highly specialized equipment or cell type specific protocols and are relatively low throughput. Here, we have combined the advantages of tagmentation, linear amplification and combinatorial indexing to produce a high throughput single cell DNA binding site mapping method that does not require specialized equipment and is capable of multiplexing several samples/ target proteins in one experiment. Targeted Insertion of Promoters (TIP-seq) uses Tn5 fused to protein A (as with CUT&Tag) to insert a T7 RNA polymerase promoter into sites adjacent to an antibody bound to a chromatin protein of interest, followed by linear amplification of flanking DNA with T7 polymerase, cDNA preparation and PCR indexing. Tip-seq provides [~]10-fold higher unique reads and thus higher coverage per single cell compared to state-of-the-art methods. We apply TIP-seq to map histone modifications, RNA PolII and CTCF binding sites in single human and mouse cells. TIP-seq will also be adaptable for other platforms, such as 10X genomics and ICELL8. In summary, TIP-seq provides a high-throughput, low-cost method for single cell protein mapping, that yields substantially higher coverage per cell and signal to noise than existing methods.

category: genomics


title: Cell-Tak coating may cause mis-normalization of Seahorse metabolic flux data

abstract:
Metabolic flux investigations of cells and tissue samples are a rapidly advancing tool in diverse research areas. Reliable methods of data normalization are crucial for an adequate interpretation of results and to avoid a misinterpretation of experiments and incorrect conclusions. The most common methods for metabolic flux data normalization are to cell number, DNA and protein. Data normalization may be affected by a variety of factors, such as density, healthy state, adherence efficiency, or proportional seeding of cells.

The mussel-derived adhesive Cell-Tak is often used to immobilize poorly adherent cells. Here we demonstrate that this coating may strongly affect the fluorescent detection of DNA leading to an incorrect and highly variable normalization of metabolic flux data. Protein assays are much less affected and cell counting can virtually completely remove the effect of the coating. Cell-Tak coating also affects cell shape in a cell line-specific manner and may change cellular metabolism.

Based on these observations we recommend cell counting as a gold standard normalization method for Seahorse metabolic flux measurements with protein content as a reasonable alternative.


title: Ion dynamics at the energy-deprived tripartite synapse

abstract:
The anatomical and functional organization of neurons and astrocytes at  tripartite synapses is essential for reliable neurotransmission, which critically depends on ATP. In low energy conditions, synaptic transmission fails, accompanied by a breakdown of ion gradients, changes in membrane potentials and cell swelling. The resulting cellular damage and cell death are causal to the often devastating consequences of an ischemic stroke. The severity of ischemic damage depends on the age and the brain region in which a stroke occurs, but the reasons for this differential vulnerability are far from understood. In the present study, we address this question by developing a comprehensive biophysical model of a glutamatergic synapse to identify key determinants of synaptic failure during energy deprivation. Our model is based on fundamental biophysical principles, includes dynamics of the most relevant ions, i.e., Na+, K+, Ca2+, Cl- and glutamate, and is calibrated with experimental data. It confirms the critical role of the Na+/K+-ATPase in maintaining ion gradients, membrane potentials and cell volumes. Our simulations demonstrate that the system exhibits two stable states, one physiological and one pathological. During energy deprivation, the physiological state may disappear, forcing a transit to the pathological state, which can be reverted when blocking voltage-gated Na+ and K+ channels. Our model predicts that the transition to the pathological state is favoured if the extracellular space fraction is small. A reduction in the extracellular space volume fraction, as, e.g. observed with ageing, will thus promote the brains susceptibility to ischemic damage. Our work provides new insights into the brains ability to recover from energy deprivation, with translational relevance for diagnosis and treatment of ischemic strokes.

Author summaryThe brain consumes energy to keep ion concentrations at normal working conditions. In the case of energy deprivation (ED), e.g., during a stroke, synaptic communication fails first. Inspired by our recent experimental work on ED, we formulated a novel computational model to explore initial events during ED. Our model reproduces time courses for several ions from different experimental data. In some cases, the system returns to baseline upon restoring energy supply. In others, we observe that neurons and astrocytes cannot recover accompanied by cell swelling. There is a threshold depending on the depth and duration of ATP depletion differentiating these cases. Also, smaller extracellular spaces hamper recovery more. This result may explain clinical observations of increased vulnerability to stroke as the size of the extracellular space shrinks with ageing.

category: biophysics


title: Scaling down the microbial loop: data-driven modelling of growth interactions in a diatom-bacterium co-culture

abstract: An intricate set of interactions characterizes marine ecosystems. One of the most important is represented by the so-called microbial loop, which includes the exchange of dissolved organic matter (DOM) from phototrophic organisms to heterotrophic bacteria. Here, it can be used as the major carbon and energy source. Arguably, this interaction is one of the foundations of the entire ocean food-web. Carbon fixed by phytoplankton can be redirected to bacterial cells in two main ways; either i) bacteria feed on dead (eventually lysed) phytoplankton cells or ii) DOM is actively released by phytoplankton cells (a widespread process that may result in up to 50% of the fixed carbon leaving the cell). In this work, we have set up a co-culture of the model diatom Phaeodactylum tricornutum and the model chemoheterotrophic bacterium Pseudoalteromonas haloplanktis TAC125 and used this system to study the interactions between these two representatives of the microbial loop. We show that the bacterium can indeed thrive on diatom-derived carbon and that this growth can be sustained by both diatom dead cells and diatom-released compounds. These observations were formalized in a network of putative interactions between P. tricornutum and P. haloplanktis and implemented in a mathematical model that reproduces the observed co-culture dynamics, suggesting that our hypotheses on the interactions occurring in this two-player system can accurately explain the experimental data.


title: The long-time orphan protist Meringosphaera mediterranea Lohmann, 1902  is a centrohelid heliozoan

abstract: Meringosphaera is an enigmatic marine protist without clear phylogenetic affiliation, but it has long been suggested to be a chrysophytes-related autotroph. Microscopy-based reports indicate that it has a worldwide distribution, but no sequence data exists so far. We obtained the first 18S rDNA sequence for M. mediterranea (identified using light and electron microscopy) from the West Coast of Sweden. Observations of living cells revealed granulated axopodia and up to 6 globular photosynthesizing bodies about 2 m in diameter, the nature of which requires further investigation. The ultrastructure of barbed undulating spine scales and patternless plate scales with a central thickening is in agreement with previous reports. Molecular phylogenetic analysis placed M. mediterranea inside the NC5 environmental clade of Centroplasthelida (Haptista) along with additional environmental sequences, together closely related to Choanocystidae. This placement is supported by similar scales in Meringosphaera and Choanocystidae. We searched the Tara Oceans 18S-V9 metabarcoding dataset which revealed four OTUs with 95.5-98.5% similarity, with oceanic distribution similar to that based on morphological observations. The current taxonomic position and species composition of the genus are discussed. The planktonic lifestyle of M. mediterranea contradicts the view of some authors that centrohelids enter the plankton only temporarily.

category: zoology


title: Dopamine and reward-related vigor in younger and older human participants

abstract: Vigor reflects how motivated one is to respond to a stimulus. We previously showed that humans are more vigorous when more reward is available on average, and that this relationship is modulated by the dopamine precursor levodopa. Dopamine signalling and probabilistic reward learning degrade with age, so the relationship between vigor and reward should change with age. We test this and assess whether the relationship between vigor and reward correlates with D1 dopamine receptor availability measured using Positron Emission Tomography. We measured response times of 30 older and 30 younger subjects during an oddball discrimination task where rewards varied systematically between trial. Reward rate had a similar impact on the vigor of both groups. We observed a weak positive association across subjects between ventral striatal dopamine receptor availability and effect of average reward rate on response time, which was in the opposite direction to our prediction. Overall, the effect of reward on response vigor is similar between younger and older humans and is weakly sensitive to dopamine D1 receptor availability.


title: Mechanical bistability enabled by ectodermal compression facilitates Drosophila mesoderm invagination

abstract: Apical constriction driven by non-muscle myosin II ("myosin") provides a well-conserved mechanism to mediate epithelial folding. It remains unclear how contractile forces near the apical surface of a cell sheet drive out-of-the-plane bending of the sheet and whether myosin contractility is required throughout folding. By optogenetic-mediated acute inhibition of myosin, we find that during Drosophila mesoderm invagination, myosin contractility is critical to prevent tissue relaxation during the early, "priming" stage of folding but is dispensable for the actual folding step after the tissue passes through a stereotyped transitional configuration. The binary response suggests that the mesoderm is mechanically bistable during gastrulation. Combined modeling analysis and experimental measurements suggest that the observed mechanical bistability may arise from apicobasal shrinkage of the surrounding ectoderm, which promotes mesoderm invagination by facilitating a buckling transition. Our results suggest that Drosophila mesoderm invagination requires a joint action of local myosin contractility and mechanical bistability of the epithelium to trigger epithelial buckling.


title: Mechanical bistability of the mesoderm epithelium facilitates mesoderm invagination during Drosophila gastrulation


title: Environmental persistence but not per capita transmission rates of a chytrid fungus determines infection prevalence

abstract: AO_SCPLOWBSTRACTC_SCPLOWUnderstanding local-scale variability in disease dynamics can be important for informing strategies for surveillance and management. For example, the amphibian chytrid fungus (Batrachochytrium dendrobatidis; Bd), which is implicated in population declines and species extinctions of amphibians, causes spatially variable epizootics and extirpations of its hosts. Outbreak heterogeneity could be driven by differential survival of zoospores, the free-living infectious life stage of Bd, or the persistence of dead zoospores and/or its metabolites in water, which could induce resistance among hosts. To gain a mechanistic understanding of the potential for variation in local transmission dynamics of Bd, we conducted Bd survival and infection experiments and then fit models to discern how Bd mortality, decomposition, and per-capita transmission rate vary among water sources. We found that infection prevalence differed among water sources, which was driven by differences in mortality rates of Bd zoospores, rather than differences in per-capita transmission rates. Specifically, zoospore mortality rates varied significantly among pond water treatments and were lower in artificial spring water compared to pond water sources. These results suggest that variation in Bd infection dynamics could be a function of differences in exposure of hosts to live Bd. In contrast to the persistence of live zoospores, we found that rates of decomposition of dead zoospores did not vary among water sources. These results may suggest that exposure of hosts to dead Bd or its metabolites, which have been shown to induce acquired resistance, might not commonly vary among nearby sites. Ultimately, a mechanistic understanding of the drivers of variable epizootics of Bd could lead to increases in the effectiveness of surveillance and management strategies.


title: Optimizing smartphone-based canopy hemispherical photography

abstract:
O_LIHemispherical photography (HP) is a long-standing tool for forest canopy characterization. Currently, there are low-cost fisheye lenses to convert smartphones into high-portable HP equipment; however, they cannot be used whenever since HP is sensitive to illumination conditions. To obtain sound results outside diffuse light conditions, a deep-learning-based system needs to be developed. A ready-to-use alternative is the multiscale color-based binarization algorithm, but it can provide moderate-quality results only for open forests. To overcome this limitation, I propose coupling it with the model-based local thresholding algorithm. I call this coupling the MBCB approach.
C_LIO_LIMethods presented here are part of the R package CAnopy IMage ANalysis (caiman), which I am developing. The accuracy assessment of the new MBCB approach was done with data from a pine plantation and a broadleaf native forest.
C_LIO_LIThe coefficient of determination (R2) was greater than 0.7, and the root mean square error (RMSE) lower than 20 %, both for plant area index calculation.
C_LIO_LIResults suggest that the new MBCB approach allows the calculation of unbiased canopy metrics from smartphone-based HP acquired in sunlight conditions, even for closed canopies. This facilitates large-scale and opportunistic sampling with hemispherical photography.
C_LI


title: Distance decay 2.0 - a global synthesis of taxonomic and functional turnover in ecological communities.

abstract: Understanding the variation in community composition and species abundances, i.e., {beta}-diversity, is at the heart of community ecology. A common approach to examine {beta}-diversity is to evaluate directional turnover in community composition by measuring the decay in the similarity among pairs of communities along spatial or environmental distances. We provide the first global synthesis of taxonomic and functional distance decay along spatial and environmental distance by analysing 149 datasets comprising different types of organisms and environments. We modelled an exponential distance decay for each dataset using generalized linear models and extracted r2 and slope to analyse the strength and the rate of the decay. We studied whether taxonomic or functional similarity has stronger decay across the spatial and environmental distances. We also unveiled the factors driving the rate of decay across the datasets, including latitude, spatial extent, realm, and organismal features. Taxonomic distance decay was stronger along spatial and environmental distances compared with functional distance decay. The rate of taxonomic spatial distance decay was the fastest in the datasets from mid-latitudes while the rate of functional decay increased with latitude. Overall, datasets covering larger spatial extents showed a lower rate of decay along spatial distances but a higher rate of decay along environmental distances. Marine ecosystems had the slowest rate of decay. This synthesis is an important step towards a more holistic understanding of patterns and drivers of taxonomic and functional {beta}-diversity.


title: Pooled Single-Molecule transcriptomics identifies a giant gene under balancing selection in sunflower

abstract: Genes under balancing selection control phenotypes such as immunity, color or sex, but are difficult to identify. Self-incompatibility genes are under negative frequency-dependent selection, a special case of balancing selection, with up to 30 to 50 alleles segregating per population. We developed a method based on pooled Single-Molecule transcriptomics to identify balanced polymorphisms expressed in tissues of interest. We searched for multi-allelic, non-recombining genes causing self-incompatibility in wild sunflower (Helianthus annuus). A diversity scan in pistil identified a gene, Ha7650b,that displayed balanced polymorphism and colocalized with a quantitative trait locus for self-incompatibility. Unexpectedly, Ha7650b displayed gigantism (400 kb), which was caused by increase in intron size as a consequence of suppressed recombination. Ha7650b emerged after a whole-genome duplication (29 millions years ago) followed by tandem duplications and neofunctionalisation. Ha7650b shows expression, genetic location, genomic neighbourhood and predicted function that provide strong evidence that it is involved in self-incompatibility. Pooled Single-Molecule transcriptomics is an affordable and powerful new method that makes it possible to identify diversity and structural outliers simultaneously. It will allow a breakthrough in the discovery of self-incompatibility genes and other expressed genes under balancing selection.

category: evolutionary biology


title: Genomic and Chemical Diversity of Commercially Available Industrial Hemp Accessions

abstract: High consumer demand for cannabidiol (CBD) has made industrial hemp (Cannabis sativa) an extremely high-value crop. However, high demand has resulted in the industry developing faster than the research, resulting in the sale of many hemp accessions with inconsistent performance and chemical profiles. To determine the genetic and phenotypic consistency in available CBD hemp varieties, we obtained seed or clones from 22 different named accessions. Genotypes (~48,000 SNPs) and chemical profiles (% CBD and THC by dry weight) were determined for up to 8 plants per accession. Many accessions--including several with the same name--showed little consistency either genetically or chemically. Most seed-grown accessions deviated significantly from their purported levels of CBD and THC based on the supplied certificates of analysis. Several also showed evidence of an active tetrahydrocannabinolic acid (THCa) synthase gene, leading to unacceptably high levels of THC in female flowers. We conclude that the current market for CBD-hemp varieties is highly unreliable, making many purchases risky for growers. We suggest options for addressing these issues, such using unique names and developing seed and plant certification programs to ensure the availability of high-quality, verified planting materials.


title: Fast and accurate genome-wide predictions and structural modeling of protein-protein interactions using Galaxy

abstract:
Protein-protein interactions play a crucial role in almost all cellular processes. Identifying interacting proteins reveals insight into living organisms and yields novel drug targets for disease treatment. Here, we present a publicly available, automated pipeline to predict genome-wide protein-protein interactions and produce high-quality multimeric structural models.

Application of our method to the Human and Yeast genomes yield protein-protein interaction networks similar in quality to common experimental methods. We identified and modeled Human proteins likely to interact with the papain-like protease of SARS-CoV2s non-structural protein 3 (Nsp3). We also produced models of SARS-CoV2s spike protein (S) interacting with myelin-oligodendrocyte glycoprotein receptor (MOG) and dipeptidyl peptidase-4 (DPP4). The presented method is capable of confidently identifying interactions while providing high-quality multimeric structural models for experimental validation.

The interactome modeling pipeline is available at usegalaxy.org and usegalaxy.eu.


title: Nanoparticle Cellular Internalization is Not Required for RNA Delivery to Mature Plant Leaves

abstract: Rapidly growing interest in nanoparticle-mediated delivery of DNA and RNA to plants requires a better understanding of how nanoparticles and their cargoes translocate in plant tissues and into plant cells. However, little is known about how the size and shape of nanoparticles influences transport in plants and use of their cargoes, limiting development and deployment of nanotechnology in plant systems. Here, we employ non-biolistically delivered DNA-modified gold nanoparticles (AuNP) spanning various sizes (5 - 20 nm) and shapes (spheres and rods) to systematically investigate their transport following infiltration into Nicotiana benthamiana (Nb) leaves. Generally, smaller AuNPs demonstrate more rapid, higher, and longer-lasting levels of association with plant cell walls compared to larger AuNPs. We observe internalization of rod-shaped but not spherical AuNPs into plant cells, yet surprisingly, 10 nm spherical AuNP functionalized with small-interfering RNA (siRNA) are most efficient at siRNA delivery and inducing gene silencing in mature plant leaves. These results indicate the importance of nanoparticle size in efficient biomolecule delivery, and, counterintuitively, demonstrate that efficient cargo delivery is possible and potentially optimal in the absence of nanoparticle cellular internalization. Our results highlight nanoparticle features of importance for transport within plant tissues, providing a mechanistic overview of how nanoparticles can be designed to achieve efficacious bio-cargo delivery for future developments in plant nanobiotechnology.


title: SCAPTURE: a deep learning-embedded pipeline that captures polyadenylation information from 3' tag-based RNA-seq of single cells

abstract: Single-cell RNA-seq (scRNA-seq) profiles gene expression with a resolution that empowers depiction of cell atlas in complex systems. Here, we developed a stepwise computational pipeline SCAPTURE to identify, evaluate, and quantify cleavage and polyadenylation sites (PASs) from 3 tag-based scRNA-seq. SCAPTURE detects PASs de novo in single cells with high sensitivity and accuracy, enabling detection of previously unannotated PASs. Quantified alternative PAS transcripts refine cell identities, enriching information extracted from scRNA-seq.


title: Assigning confidence to structural annotations from mass spectra with COSMIC

abstract: Untargeted metabolomics experiments rely on spectral libraries for structure annotation, but these libraries are vastly incomplete; in silico methods search in structure databases but cannot distinguish between correct and incorrect annotations. As biological interpretation relies on accurate structure annotations, the ability to assign confidence to such annotations is a key outstanding problem. We introduce the COSMIC workflow that combines structure database generation, in silico annotation, and a confidence score consisting of kernel density p-value estimation and a Support Vector Machine with enforced directionality of features. In evaluation, COSMIC annotates a substantial number of hits at small false discovery rates, and outperforms spectral library search for this purpose. To demonstrate that COSMIC can annotate structures never reported before, we annotated twelve novel bile acid conjugates; nine structures were confirmed by manual evaluation and two structures using synthetic standards. Second, we annotated and manually evaluated 315 molecular structures in human samples currently absent from the Human Metabolome Database. Third, we applied COSMIC to 17,400 experimental runs and annotated 1,715 structures with high confidence that were absent from spectral libraries.


title: A structural property for reduction of biochemical networks

abstract: Large-scale biochemical models are of increasing sizes due to the consideration of interacting organisms and tissues. Model reduction approaches that preserve the flux phenotypes can simplify the analysis and predictions of steady-state metabolic phenotypes. However, existing approaches either restrict functionality of reduced models or do not lead to significant decreases in the number of modelled metabolites. Here, we introduce an approach for model reduction based on the structural property of balancing of complexes that preserves the steady-state fluxes supported by the network and can be efficiently determined at genome scale. Using two large-scale mass-action kinetic models of Escherichia coli, we show that our approach results in a substantial reduction of 99% of metabolites. Applications to genome-scale metabolic models across kingdoms of life result in up to 55% and 85% reduction in the number of metabolites when arbitrary and mass-action kinetics is assumed, respectively. We also show that predictions of the specific growth rate from the reduced models match those based on the original models. Since steady-state flux phenotypes from the original model are preserved in the reduced, the approach paves the way for analysing other metabolic phenotypes in large-scale biochemical networks.


title: Spectracular: minimizing spectral overlap in multicolor flow cytometry experiments

abstract: Selecting fluorochromes for polychromatic panels for flow cytometry is complex and time-consuming. Poorly designed panels can result in overlap between the emission spectra of the different fluorochromes, making their signals difficult to separate. While assessing all possible panels is simple to do programmatically, the combinatorial complexity of most real world problems renders brute force computation impractical. Here we present a novel complexity optimization algorithm for fast design of minimal spectral overlap fluorochrome combinations. To aid researchers in designing fluorochrome panels we implemented the algorithm in a web server, Spectracular, that also considers instrument laser configuration and allows users to define various constraints related to antibody-fluorochrome availability and marker co-expression.


title: Minimizing spectral overlap in multicolor flow cytometry experiments


title: CHRONIC INFLAMMATION IN ULCERATIVE COLITIS CAUSES LONG TERM CHANGES IN GOBLETCELL FUNCTION

abstract:
ObjectiveOne of the features of ulcerative colitis (UC) is a defect in the protective mucus layer. This has been attributed to a reduced number of goblet cells (GC). However, it is not known whether abnormal GC mucus secretion also contributes to the reduced mucus layer. Our aims were to test the hypothesis that GC secretion was abnormal in UC with the changes persistent in colonoids even in the absence of immune cells.

DesignColonoids were established from intestinal stem cells of healthy subjects (HS) and from patients with UC (inactive and active sites). Colonoids were maintained as undifferentiated (UD) or induced to differentiate (DF) and studied as 3D or monolayers on Transwell filters. Total RNA was extracted for quantitative real-time PCR analysis. Carbachol and PGE2 mediated stimulation followed by examination of mucus layer by MUC2 IF/confocal microscopy and TEM were performed.

ResultsColonoids derived from patients with UC can be propagated over many passages; however, they exhibit a reduced rate of growth and TEER compared with colonoids from HS. Differentiated UC colonoid monolayers form a thin and non-continuous mucus layer. UC colonoids have increased expression of secretory lineage markers: ATOH1 and SPDEF, including MUC2 positive GCs and ChgA positive enteroendocrine cells but failed to secrete mucin when exposed to the cholinergic agonist carbachol and PGE2, which caused increased secretion in HS. Exposure to TNF- (5days), reduced the number of GC with a greater percentage decrease in UC colonoids compared to HS.

ConclusionsAbnormal mucus layer in UC is due to long term changes in epithelial cells that lead to abnormal mucus secretion as well as effects of the inflammatory environment to reduce the number of GC. This continued defect in GC mucus secretion may be involved in UC recurrence.


title: Predicting tomato field-yield using continuous monitoring of young tomato water status

abstract: To address the challenge of predicting tomato yields in the field, we used whole-plant functional phenotyping to evaluate water relations under well-irrigated and drought conditions. The genotypes tested are known to exhibit variability in their yields in wet and dry fields. The examined lines included two lines with recessive mutations that affect carotenoid biosynthesis, zeta z2083 and tangerine t3406, both isogenic to the processing tomato variety M82. The two mutant lines were reciprocally grafted onto M82, and multiple physiological characteristics were measured continuously, before, during and after drought treatment in the greenhouse. A comparative analysis of greenhouse and field yields showed that the whole-canopy stomatal conductance (gsc) in the morning and cumulative transpiration (CT) were strongly correlated with field measurements of total yield (TY: r2 = 0.9 and 0.77, respectively) and plant vegetative weight (PW: r2 = 0.6 and 0.94, respectively). Furthermore, the minimum CT during drought and the rate of recovery when irrigation was resumed were both found to predict resilience.


title: The potential of dynamic physiological traits in young tomato plants to predict field-yield performance


title: Predicting intelligence from fMRI data of the human brain in a few minutes of scan time

abstract: In recent years, the prediction of individual behaviour from the fMRI-based functional connectome has become a major focus of research. The motivation behind this research is to find generalizable neuromarkers of cognitive functions. However, insufficient prediction accuracies and long scan time requirements are still unsolved issues. Here we propose a new machine learning algorithm for predicting intelligence scores of healthy human subjects from resting state (rsfMRI) or task-based fMRI (tfMRI). In a cohort of 390 unrelated test subjects of the Human Connectome Project, we found correlations between the observed and the predicted general intelligence of more than 50 percent in tfMRI, and of around 59 percent when results from two tasks are combined. Surprisingly, we found that the tfMRI data were significantly more predictive of intelligence than rsfMRI even though they were acquired at much shorter scan times (approximately 10 minutes versus 1 hour). Existing methods that we investigated in a benchmark comparison underperformed on tfMRI data and produced prediction accuracies well below our results. Our proposed algorithm differs from existing methods in that it achieves dimensionality reduction via ensemble learning and partial least squares regression rather than via brain parcellations or ICA decompositions. In addition, it introduces Ricci-Forman curvature as a novel type of edge weight.

category: neuroscience


title: Fibrinogen and Complement Factor H are promising CSF protein biomarker(s) for Parkinson's disease with cognitive impairment- A Proteomics and ELISA based study

abstract: Cognitive impairment is a debilitating non-motor symptom of Parkinsons disease (PD). The diagnosis of PD with cognitive impairment (PDCI) is essentially through clinical and neuropsychological examinations. There is an emerging need to identify biomarker(s) to foresee cognitive decline in PD patients, at an early stage. We performed label-free unbiased nontargeted proteomics (Q-TOF LC/MS-MS) in CSF of non-neurological control (NNC); PDCI; PD and normal pressure hydrocephalus (NPH), followed by targeted ELISA for validation. The diagnosis was confirmed by neuropsychological and MRI assessments prior to CSF collection. Of the 282 proteins identified by mass spectrometry, 42 were differentially altered in PD, PDCI and NPH. Further, 28 proteins were altered in PDCI and 25 in NPH. An interesting overlap of certain proteins was noted both in PDCI and NPH. Five significantly upregulated proteins in PDCI were fibrinogen, gelsolin, complement factor-H, apolipoprotein A-IV and apolipoprotein A-I. Whereas carnosine dipeptidase 1, carboxypeptidase E, dickkpof 3 and secretogranin 3 precursor proteins were down-regulated. NPH also had few uniquely altered proteins viz. insulin-like growth factor-binding protein, ceruloplasmin, -1 antitrypsin, VGF nerve growth factor, neural cell adhesion molecule L1 like protein. Interestingly, the ELISA-derived protein concentrations correlated well with the neuropsychological scores of certain cognitive domains. Executive function was affected both in PDCI and NPH. In PD, Wisconsin card sorting test (WCST) percentile correlated positively with ApoA-IV and negatively with the ratio of ApoAI: ApoA-IV. Thus assessment of a battery of proteins like fibrinogen--chain, CFAH and ApoAI: ApoA-IV ratio alongside neuropsychological could be reliable biomarkers to distinguish PDCI and NPH.


title: A new drug candidate for glycogen storage disorders enhances glycogen catabolism: Lessons from Adult Polyglucosan Body Disease models

abstract:
This work employs Adult Polyglucosan Body Disease (APBD) models to explore the efficacy and mechanism of action of 144DG11, a new polyglucosan-reducing lead compound discovered by a high-throughput screen (HTS). APBD is an adult onset glycogen storage disorder (GSD) manifesting as a debilitating progressive axonopathic leukodystrophy. APBD is caused by glycogen branching enzyme (GBE) deficiency leading to poorly branched and insoluble glycogen inclusions, which precipitate as neuropathogenic polyglucosans (PG). 144DG11 led to prolonged survival and improved motor parameters in a GBE knockin (Gbeys/ys) APBD mouse model. Histopathologically, 144DG11 reduced PG and glycogen levels in brain, liver, heart, and peripheral nerve. Indirect calorimetry experiments revealed that 144DG11 increases carbohydrate burn at the expense of fat burn, suggesting metabolic mobilization of pathogenic PG. These results were also reflected at the cellular level by increased glycolytic, mitochondrial and total ATP production. Mechanistically, we show that the molecular target of 144DG11 is the lysosomal membrane protein LAMP1, whose interaction with the compound, similar to LAMP1 knockdown, enhanced autolysosomal degradation of glycogen and lysosomal acidification. Enhanced mitochondrial activity and lysosomal modifications were also the most pronounced effects of 144DG11 in APBD patient fibroblasts as discovered by image-based multiparametric phenotyping analysis and corroborated by proteomics. In summary, this work presents a broad mechanistic and target-based characterization of 144DG11 in in vivo and cell models of the prototypical GSD APBD. This investigation warrants development of 144DG11 into a safe and efficacious GSD therapy.

One Sentence SummaryA new compound, demonstrated to ameliorate APBD in vivo and ex vivo by autophagic catabolism of glycogen, may potentially become a universal drug for glycogen storage disorders.


title: Style Transfer Using Generative Adversarial Networks for Multi-Site MRI Harmonization

abstract: Large data initiatives and high-powered brain imaging analyses require the pooling of MR images acquired across multiple scanners, often using different protocols. Prospective cross-site harmonization often involves the use of a phantom or traveling subjects. However, as more datasets are becoming publicly available, there is a growing need for retrospective harmonization, pooling data from sites not originally coordinated together. Several retrospective harmonization techniques have shown promise in removing cross-site image variation. However, most unsupervised methods cannot distinguish between image-acquisition based variability and cross-site population variability, so they require that datasets contain subjects or patient groups with similar clinical or demographic information. To overcome this limitation, we consider cross-site MRI image harmonization as a style transfer problem rather than a domain transfer problem. Using a fully unsupervised deep-learning framework based on a generative adversarial network (GAN), we show that MR images can be harmonized by inserting the style information encoded from a reference image directly, without knowing their site/scanner labels a priori. We trained our model using data from five large-scale multi-site datasets with varied demographics. Results demonstrated that our styleencoding model can harmonize MR images, and match intensity profiles, successfully, without relying on traveling subjects. This model also avoids the need to control for clinical, diagnostic, or demographic information. Moreover, we further demonstrated that if we included diverse enough images into the training set, our method successfully harmonized MR images collected from unseen scanners and protocols, suggesting a promising novel tool for ongoing collaborative studies.


title: Limited effects of m6A modification on mRNA partitioning into stress granules

abstract: Recent studies have argued that the m6A modification of mRNAs promotes mRNA recruitment to stress granules through the interaction with YTHDF proteins (Anders et al., 2018; Ries et al., 2019). However, mRNAs that contain multiple m6A modified sites partition similarly into stress granules in both wild-type and m6A-deficient cells by single-molecule FISH suggesting m6A modifications play a minor role in mRNA partitioning into stress granules. Moreover, multiple linear regression analysis suggests m6A modification plays a minimal role in stress granule recruitment. Finally, the artificial tethering of 25 YTHDF proteins on reporter mRNAs leads to only a modest increase in mRNA partitioning to stress granules. These results indicate m6A modification makes a small, but measurable, contribution to recruiting specific mRNAs to stress granules.


title: Sponges are highly resistant to radiation exposure and cancer

abstract: There are no reports of cancer in sponges, despite them having somatic cell turnover, long lifespans and no specialized adaptive immune cells. In order to investigate whether sponges are cancer resistant, we exposed a species of sponge, Tethya wilhelma, to X-rays. We found that T. wilhelma can withstand 600 Gy of X-ray radiation. That is approximately 100 times the lethal dose for humans. A single high dose of X-rays did not induce cancer in sponges, providing the first experimental evidence of cancer resistance in the phylum, Porifera. Following X-ray exposure, we found an overexpression of genes involved in DNA repair, signaling transduction pathways and epithelial to mesenchymal transition. Sponges have the highest level of radiation resistance that has yet been observed in animals that have sustained somatic cell turnover. This may make them an excellent model system for studying cancer resistance and developing new approaches for cancer prevention and treatment.


title: Predicting continuous ground reaction forces from accelerometers during uphill and downhill running: A recurrent neural network solution

abstract:
Ground reaction forces (GRFs) are important for understanding human movement, but their measurement is generally limited to a laboratory. Previous studies used neural networks to predict GRF waveforms during running from wearable device data, but these predictions are limited to the stance phase of level-ground running. We sought to develop a recurrent neural network capable of predicting continuous normal (perpendicular to surface) GRFs across a range of running speeds and slopes from accelerometer data.

19 subjects ran on a force-measuring treadmill at five slopes (0{degrees}, {+/-}5{degrees}, {+/-}10{degrees}) and three speeds (2.5, 3.33, 4.17 m/s) per slope with sacral- and shoe-mounted accelerometers. We then trained a recurrent neural network to predict normal GRF waveforms frame-by-frame. The predicted versus measured GRF waveforms had an average {+/-} SD RMSE of 0.16 {+/-} 0.04 BW and relative RMSE of 6.4 {+/-} 1.5% across all conditions and subjects.

The recurrent neural network predicted continuous normal GRF waveforms across a range of running speeds and slopes with greater accuracy than neural networks implemented in previous studies. This approach may facilitate predictions of biomechanical variables outside the laboratory in near real-time and improves the accurately of quantifying and monitoring external loads experienced by the body when running.


title: A Deep Learning Model for Molecular Label Transfer that Enables Cancer Cell Identification from Histopathology Images

abstract: Deep learning cancer classification systems have the potential to improve cancer diagnosis. However, development of these computational approaches depends on prior annotation through a pathologist. This initial step relying on a manual, low-resolution, time-consuming process is highly variable and subject to observer variance. To address this issue, we developed a novel method, H&E Molecular neural network (HEMnet). This two-step process utilises immunohistochemistry as an initial molecular label for cancer cells on a H&E image and then we train a cancer classifier on the overlapping clinical histopathological images. Using this molecular transfer method, we show that HEMnet accurately distinguishes colorectal cancer from normal tissue at high resolution without the need for an initial manual histopathologic evaluation. Our validation study using histopathology images from TCGA samples accurately estimates tumour purity. Overall, our method provides a path towards a fully automated delineation of any type of tumor so long as there is a cancer-oriented molecular stain available for subsequent learning. Software, tutorials and interactive tools are available at: https://github.com/BiomedicalMachineLearning/HEMnet


title: Kawasaki Disease patient stratification and pathway analysis based on host transcriptomic and proteomic profiles

abstract: The aetiology of Kawasaki Disease (KD), an acute inflammatory disorder of childhood, remains unknown despite various triggers of KD having been proposed. Host  omic profiles offer insights into the host response to infection and inflammation, with the interrogation of multiple  omic levels in parallel providing a more comprehensive picture. We used differential abundance analysis, pathway analysis, clustering and classification techniques to explore whether the host response in KD is more similar to the response to bacterial or viral infection at the transcriptomic and proteomic levels through comparison of  omic profiles from children with KD to those with bacterial and viral infections. Pathways activated in patients with KD included those involved in anti-viral and anti-bacterial responses. Unsupervised clustering showed that the majority of KD patients clustered with bacterial patients on both  omic levels, whilst application of diagnostic signatures specific for bacterial and viral infections revealed that many transcriptomic KD samples had low probabilities of having bacterial or viral infections, suggesting that KD may be triggered by a different process not typical of either common bacterial or viral infections. Clustering based on the transcriptomic and proteomic responses during KD revealed three clusters of KD patients on both  omic levels, suggesting heterogeneity within the inflammatory response during KD. The observed heterogeneity may reflect differences in the host response to a common trigger, or variation dependent on different triggers of the condition.