title: Autocrine regulation of adult neurogenesis by the endocannabinoid 2-arachidonoylglycerol (2-AG)

abstract:
The endocannabinoid system modulates adult hippocampal neurogenesis by promoting the proliferation and survival of neural stem and progenitor cells (NSPCs). Specifically, deleting cannabinoid receptors on NSPCs or the constitutive deletion of the endocannabinoid 2-arachidonoylglycerol (2-AG) producing enzyme diacylglycerol lipase alpha (DAGLa) disrupts adult neurogenesis. However, it is not known which cells are the producers of 2-AG relevant to neurogenesis. In this paper, we investigated the cellular source of endocannabinoids in the subgranular zone (SGZ) of the hippocampus, an important neurogenic niche. For this purpose, we used two complementary Cre-deleter mouse strains to delete DAGLa either in neurons or astroglia and NSPCs. Surprisingly, neurogenesis was not altered in mice with a deletion of Dagla in neurons (Syn-Dagla KO), although they are the main source for the endocannabinoids in the brain. In contrast, mice with a specific inducible deletion of Dagla in NPSCs and astrocytes (GLAST-CreERT2-Dagla KO) showed a strongly impaired neurogenesis with significantly reduced proliferation and survival of newborn cells. These results identify Dagla in NSPCs in the SGZ of dentate gyrus or in astrocytes, as the cellular source for 2-AG in adult hippocampal neurogenesis. In summary, 2-AG produced by progenitor cells or astrocytes in the SGZ regulates adult hippocampal neurogenesis.

SummaryDAGLa in neuronal progenitor cells in the SGZ of dentate gyrus is identified as the cellular source for 2-AG in adult hippocampal neurogenesis.


title: REAP: A platform to identify autoantibodies that target the human exoproteome

abstract: Autoantibodies that recognize extracellular proteins (the "exoproteome") exert potent biological effects but have proven challenging to detect with existing screening technologies. Here, we developed Rapid Extracellular Antigen Profiling (REAP) as a technique for comprehensive, high-throughput discovery of exoproteome-targeting autoantibodies. With REAP, patient samples are applied to a genetically-barcoded library containing 2,688 human extracellular proteins displayed on the surface of yeast. Antibody-coated cells are isolated by magnetic selection and deep sequencing of their barcodes is used to identify the displayed antigens. To benchmark the performance of REAP, we screened 77 patients with autoimmune polyendocrinopathy candidiasis ectodermal dystrophy (APECED). REAP sensitively and specifically detected known autoantibody reactivities in APECED in addition to numerous previously unidentified reactivities. We further screened 106 patients with systemic lupus erythematosus (SLE) and identified novel autoantibody reactivities against a diverse set of antigens including growth factors, extracellular matrix components, cytokines, and immunomodulatory proteins. Several of these responses were associated with disease severity and specific clinical manifestations of SLE and exerted potent functional effects on cell signaling ex vivo. These findings demonstrate the utility of REAP to atlas the expansive landscape of exoproteome-targeting autoantibodies and their impacts on patient health outcomes.

category: immunology


title: Cooperative virulence can emerge via horizontal gene transfer but is stabilized by transmission

abstract: Intestinal inflammation fuels Salmonella Typhimurium (S.Tm) transmission despite a fitness cost associated with the expression of virulence. Cheater mutants can emerge that profit from inflammation without enduring this cost. Intestinal virulence in S.Tm is therefore a cooperative trait, and its evolution a conundrum. Horizontal gene transfer (HGT) of cooperative alleles may facilitate the emergence of cooperative virulence, despite its instability. To test this hypothesis, we cloned hilD, coding for a master regulator of virulence, into a conjugative plasmid that is highly transferrable during intestinal colonization. We demonstrate that virulence can emerge by hilD transfer between avirulent strains in vivo. However, this was indeed unstable and hilD mutant cheaters arose within a few days. The timing of cheater emergence depended on the cost. We further show that stabilization of cooperative virulence in S.Tm is dependent on transmission dynamics, strengthened by population bottlenecks, leading cheaters to extinction and allowing cooperators to thrive.

category: evolutionary biology


title: Genetic determinants facilitating the evolution of resistance to carbapenem antibiotics

abstract: In this era of rising antibiotic resistance, in contrast to our increasing understanding of mechanisms that cause resistance, our understanding of mechanisms that influence the propensity to evolve resistance remains limited. Here, we identified genetic factors that facilitate the evolution of resistance to carbapenems, the antibiotic of "last resort," in Klebsiella pneumoniae, the major carbapenem resistant species. In clinical isolates, we found that high-level transposon insertional mutagenesis plays an important role in contributing to high-level resistance frequencies in several major and emerging carbapenem-resistant lineages. A broader spectrum of resistance-conferring mutations for select carbapenems such as ertapenem also enables higher resistance frequencies and importantly, creates stepping-stones to achieve high-level resistance to all carbapenems. These mutational mechanisms can contribute to the evolution of resistance, in conjunction with the loss of systems that restrict horizontal resistance gene uptake, such as the CRISPR-Cas system. Given the need for greater antibiotic stewardship, these findings argue that in addition to considering the current efficacy of an antibiotic for a clinical isolate in antibiotic selection, considerations of future efficacy are also important. The genetic background of a clinical isolate and the exact antibiotic identity can and should also be considered as it is a determinant of a strains propensity to become resistant. Together, these findings thus provide a molecular framework for understanding acquisition of carbapenem resistance in K. pneumoniae with important implications for diagnosing and treating this important class of pathogens.


title: Identification of a rapidly-spreading triple mutant for high-level metabolic insecticide resistance in Anopheles gambiae provides a real-time molecular diagnostic for anti-malarial intervention deployment.

abstract: Insecticide resistance provides both an increasingly pressing threat to the control of vector-borne diseases and insights into the remarkable capacity of natural populations to show rapid evolutionary responses to contemporary selection. Malaria control remains heavily dependent on deployment of pyrethroid insecticides, primarily in long lasting insecticidal nets (LLINs), but resistance in the major malaria vectors has increased over the last 15 years in concert with dramatic expansion of LLIN distributions. Identifying genetic mechanisms underlying high-level resistance in mosquitoes, which may almost entirely overcome pyrethroid efficacy, is crucial for the development and deployment of potentially resistance-breaking tools. Using the Anopheles gambiae 1000 genomes (Ag1000g) data we identified a very recent selective sweep in mosquitoes from Uganda which localized to a cluster of cytochrome P450 genes, including some commonly implicated in resistance. Further interrogation revealed a haplotype involving a trio of mutations, a nonsynonymous point mutation in Cyp6p4 (I236M), an upstream insertion of a partial Zanzibar-like transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae from the Kenya-Uganda border region around Lake Victoria, with stepwise replacement of the double-mutant (Zanzibar-like TE and Cyp6p4-236M) with the triple-mutant haplotype (including Cyp6aa1 duplication), which has spread into the Democratic Republic of Congo and Tanzania. The triple-mutant haplotype is strongly associated with increased expression of genes able to metabolise pyrethroids and is strongly predictive of resistance to pyrethroids most notably deltamethrin, a commonly-used LLIN insecticide. Importantly, there was increased mortality in mosquitoes carrying the triple-mutation when exposed to nets co-treated with the synergist piperonyl butoxide (PBO). Frequencies of the triple-mutant haplotype remain spatially variable within countries, suggesting an effective marker system to guide deployment decisions for limited supplies of PBO-pyrethroid co-treated LLINs across African countries. Duplications of the Cyp6aa1 gene are common in An. gambiae across Africa and, given the enzymes metabolic activity, are likely to be a useful diagnostic for high levels of pyrethroid resistance.


title: A benchmarking of human mitochondrial DNA haplogroup classifiers from whole-genome and whole-exome sequence data

abstract: The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short-read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

category: genomics


title: Ribovore: ribosomal RNA sequence analysis for GenBank submissions and database curation

abstract:
BackgroundThe DNA sequences encoding ribosomal RNA genes (rRNAs) are commonly used as markers to identify species, including in metagenomics samples that may combine many organismal communities. The 16S small subunit ribosomal RNA (SSU rRNA) gene is typically used to identify bacterial and archaeal species. The nuclear 18S SSU rRNA gene, and 28S large subunit (LSU) rRNA gene have been used as DNA barcodes and for phylogenetic studies in different eukaryote taxonomic groups. Because of their popularity, the National Center for Biotechnology Information (NCBI) receives a disproportionate number of rRNA sequence submissions and BLAST queries. These sequences vary in quality, length, origin (nuclear, mitochondria, plastid), and organism source and can represent any region of the ribosomal cistron.

ResultsTo improve the timely verification of quality, origin and loci boundaries, we developed Ribovore, a software package for sequence analysis of rRNA sequences. The ribotyper and ribosensor programs are used to validate incoming sequences of bacterial and archaeal SSU rRNA. The ribodbmaker program is used to create high-quality datasets of rRNAs from different taxonomic groups. Key algorithmic steps include comparing candidate sequences against rRNA sequence profile hidden Markov models (HMMs) and covariance models of rRNA sequence and secondary-structure conservation, as well as other tests. Nine freely available blastn rRNA databases created and maintained with Ribovore are used for checking incoming GenBank submissions and used by the blastn browser interface at NCBI. Since 2018, Ribovore has been used to analyze more than 50 million prokaryotic SSU rRNA sequences submitted to GenBank, and to select at least 10,435 fungal rRNA RefSeq records from type material of 8,350 taxa.

ConclusionRibovore combines single-sequence and profile-based methods to improve GenBank processing and analysis of rRNA sequences. It is a standalone, portable, and extensible software package for the alignment, classification and validation of rRNA sequences. Researchers planning on submitting SSU rRNA sequences to GenBank are encouraged to download and use Ribovore to analyze their sequences prior to submission to determine which sequences are likely to be automatically accepted into GenBank.

category: bioinformatics


title: Substrates modulate charge-reorganization allosteric effects in protein-protein association

abstract:
Protein function may be modulated by an event occurring far away from the functional site, a phenomenon termed allostery. While classically allostery involves conformational changes, we recently observed that charge redistribution within an antibody can also lead to an allosteric effect, modulating the kinetics of binding to target antigen. In the present study, we study the association of a poly-histidine tagged enzyme (phosphoglycerate kinase, PGK) to surface-immobilized anti-His antibodies, finding a significant Charge-Reorganization Allostery (CRA) effect. We further observe that PGKs negatively charged nucleotide substrates modulate CRA substantially, even though they bind far away from the His-tag-antibody interaction interface. In particular, binding of ATP reduces CRA by more than 50%. The results indicate that CRA may be affected by charged substrates bound to a protein and provide further insight into the role of charge redistribution in protein function.

TOC GRAPHIC O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

category: biophysics


title: siRNA enrichment in Argonaute 2-depleted Blattella germanica

abstract:
BackgroundRNA interference (RNAi) is a cellular mechanism used to fight various threats, including transposons, aberrant RNA, and some types of viruses. It relies on dsRNA detection and, through a mechanism involving Dicer-2 (Dcr-2) and Argonaute 2 (AGO2), together with a small RNA fragment (siRNA) as a complementary guide, binds to an RNA and cleaves it.

MethodsUsing the cockroach Blattella germanica as a model, we examined AGO2 activity by depleting its mRNA levels using RNAi and analyzing the phenotypes produced.

ResultsSilencing AGO2 expression had no remarkable effect on nymphal development or reproduction. dsRNA treatment triggered an immediate and transitory increase in AGO2 expression, independently of Dcr-2 action. In addition, we analyzed the siRNAs generated after injecting a heterologous dsRNA in control and AGO2-depleted animals. The results revealed that siRNAs were distributed non-uniformly along the dsRNA sequence. In AGO2-depleted animals, the proportion of 22 nucleotide reads was higher and accumulations of reads appeared in areas less well-represented in the controls. We also detected a slight preference for cytosine as the first nucleotide in controls to be lost in AGO2-depleted individuals.

Conclusions/General significanceThe siRNAs produced from a dsRNA are heterogeneously distributed along the length of the dsRNA and this arrangement depends on the sequence. AGO2 exerts its role as nuclease on the siRNA duplexes independently of its action on the homologous mRNA. This study sheds light on an extremely useful process for reverse genetics in laboratories, in addition to the design of more effective, specific, and eco-friendly pest-control strategies.

Highlights- RNA interference is used to fight cellular threats including transposons and viruses
- Silencing AGO2 in Blattella germanica does not affect development or reproduction
- siRNA distribution along the dsRNA sequence is non-uniform and sequence dependent
- AGO2 depletion increased the proportion of 22 nt reads in certain areas

category: physiology


title: Cone mitochondria act as microlenses to enhance light delivery and confer Stiles-Crawford-like direction sensitivity

abstract: Evolution endeavors to maximize the function of biological structures in organisms, and the vertebrate eye is no exception. Cone photoreceptors in the retina are among the most energy-demanding cells in our body, necessitating numerous mitochondria. Intriguingly, these mitochondria adopt a peculiar spatial aggregation immediately beneath the cone outer segment (OS) that houses light-sensitive opsin molecules. Here we demonstrate, via direct live imaging and computational modeling of ground squirrel cones, that such mitochondria bundles concentrate light to enter the OS for detection. This "microlens"-like feature of cone mitochondria produces an angular dependence of light intensity quantitively consistent with the Stiles-Crawford effect, a psychophysical phenomenon believed to improve visual resolution. Thus, in addition to their function as a necessary powerhouse, cone mitochondria play a critical optical role.

category: neuroscience


title: Regulation of PD1 signaling is associated with prognosis in glioblastoma multiforme

abstract: Glioblastoma is an aggressive cancer of the brain and spine. While analysis of glioblastoma  omics data has somewhat improved our understanding of the disease, it has not led to direct improvement in patient survival. Cancer survival is often characterized by differences in expression of particular genes, but the mechanisms that drive these differences are generally unknown. We therefore set out to model the regulatory mechanisms that associate with glioblastoma survival. We inferred individual patient gene regulatory networks using data from two different expression platforms from The Cancer Genome Atlas (n=522 and 431). We performed a comparative network analysis between patients with long- and short-term survival, correcting for patient age, sex, and neoadjuvant treatment status. We identified seven pathways associated with survival, all of which were involved in immune system signaling. Differential regulation of PD1 signaling was validated in an independent dataset from the German Glioma Network (n=70). We found that transcriptional repression of genes in this pathway--for which treatment options are available--was lost in short-term survivors and that this was independent of mutation burden and only weakly associated with T-cell infiltrate. These results provide a new way to stratify glioblastoma patients that uses network features as biomarkers to predict survival, and identify new potential therapeutic interventions, thus underscoring the value of analyzing gene regulatory networks in individual cancer patients.

category: systems biology


title: A single-cell atlas of lymphocyte adaptive immune repertoires and transcriptomes reveals age-related differences in convalescent COVID-19 patients

abstract: COVID-19 disease outcome is highly dependent on adaptive immunity from T and B lymphocytes, which play a critical role in the control, clearance and long-term protection against SARS-CoV-2. To date, there is limited knowledge on the composition of the T and B cell immune receptor repertoires [T cell receptors (TCRs) and B cell receptors (BCRs)] and transcriptomes in convalescent COVID-19 patients of different age groups. Here, we utilize single-cell sequencing (scSeq) of lymphocyte immune repertoires and transcriptomes to quantitatively profile the adaptive immune response in COVID-19 patients of varying age. We discovered highly expanded T and B cells in multiple patients, with the most expanded clonotypes coming from the effector CD8+ T cell population. Highly expanded CD8+ and CD4+ T cell clones show elevated markers of cytotoxicity (CD8: PRF1, GZMH, GNLY; CD4: GZMA), whereas clonally expanded B cells show markers of transition into the plasma cell state and activation across patients. By comparing young and old convalescent COVID-19 patients (mean ages = 31 and 66.8 years, respectively), we found that clonally expanded B cells in young patients were predominantly of the IgA isotype and their BCRs had incurred higher levels of somatic hypermutation than elderly patients. In conclusion, our scSeq analysis defines the adaptive immune repertoire and transcriptome in convalescent COVID-19 patients and shows important age-related differences implicated in immunity against SARS-CoV-2.


title: Rapidly Increasing SARS-CoV-2 Neutralization by Intravenous Immunoglobulins Produced from Plasma Collected During the 2020 Pandemic

abstract: Immunoglobulin (IG) lots (N=176) released since March 2020 were tested for SARS-CoV-2 neutralizing antibodies, with first positive results for September 2020 lots, mean = 1.8 IU/ml, 46% of lots positive. From there, values steadily increased, in correlation with the cumulative COVID-19 incidence, to reach a mean of 36.7 IU/ml and 93% of lots positive by January 2021. Extrapolating the correlation, IGs could reach an anti-SARS-CoV-2 potency of ~400 IU/ml by July 2021. At that stage, prophylactic IG treatment for primary/secondary immunodeficiency could contain similar doses of anti-SARS-CoV-2 as convalescent plasma which is used for treatment of COVID-19.


title: A combination of cross-neutralizing antibodies synergizes to prevent SARS-CoV-2 and SARS-CoV pseudovirus infection

abstract: Coronaviruses have caused several epidemics and pandemics including the ongoing coronavirus disease 2019 (COVID-19). Some prophylactic vaccines and therapeutic antibodies have already showed striking effectiveness against COVID-19. Nevertheless, concerns remain about antigenic drift in SARS-CoV-2 as well as threats from other sarbecoviruses. Cross-neutralizing antibodies to SARS-related viruses provide opportunities to address such concerns. Here, we report on crystal structures of a cross-neutralizing antibody CV38-142 in complex with the receptor binding domains from SARS-CoV-2 and SARS-CoV. Our structural findings provide mechanistic insights into how this antibody can accommodate antigenic variation in these viruses. CV38-142 synergizes with other cross-neutralizing antibodies, in particular COVA1-16, to enhance neutralization of SARS-CoV-2 and SARS-CoV. Overall, this study provides valuable information for vaccine and therapeutic design to address current and future antigenic drift in SARS-CoV-2 and to protect against zoonotic coronaviruses.

category: microbiology


title: Surface-aerosol stability and pathogenicity of diverse MERS-CoV strains from 2012 - 2018

abstract: Middle East Respiratory Syndrome coronavirus (MERS-CoV) is a coronavirus that infects both humans and dromedary camels and is responsible for an ongoing outbreak of severe respiratory illness in humans in the Middle East. While some mutations found in camel-derived MERS-CoV strains have been characterized, the majority of natural variation found across MERS-CoV isolates remains unstudied. Here we report on the environmental stability, replication kinetics and pathogenicity of several diverse isolates of MERS-CoV as well as SARS-CoV-2 to serve as a basis of comparison with other stability studies. While most of the MERS-CoV isolates exhibited similar stability and pathogenicity in our experiments, the camel derived isolate, C/KSA/13, exhibited reduced surface stability while another camel isolate, C/BF/15, had reduced pathogenicity in a small animal model. These results suggest that while betacoronaviruses may have similar environmental stability profiles, individual variation can influence this phenotype, underscoring the importance of continual, global viral surveillance.


title: Reproductive plasticity in both sexes interacts to determine mating behaviour and fecundity

abstract:
Organisms alter their phenotype in response to variation in their environment by expressing phenotypic plasticity. Both sexes exhibit such plasticity in response to contrasting environmental and social cues, and this can reflect the influence of sexual conflict. However, theory predicts that plasticity expressed by both sexes may either maximise the sex-specific fitness of both, or of one sex at the expense of the other. Hence empirical tests of the predictions are sorely needed. Here we conducted novel tests of the fitness effects of interacting reproductive plasticity in Drosophila melanogaster. First, prior to mating, males were kept alone, or with same sex rivals, and females were kept alone, in same sex, or mixed sex groups. Second, we conducted matings between individuals from all these social treatments under  choice and  no choice scenarios. The results showed that males and females can both plastically respond to these socio-sexual environments to influence the expression of mating duration, mating latency, and fecundity. These plastic responses interacted significantly to determine mating latency and fecundity. Effects on mating latency were also observed under both choice and no-choice conditions, but in opposing directions. Variation in the outcome of interacting plasticity pivoted around the outcomes observed with focal females that had been maintained in same-sex environments prior to mating. However, not all fitness-related traits examined responded in the same way. Mating duration was determined largely by the social environment of the male. Our results show that the expression of some, but not all fitness-related reproductive traits can be determined by the outcome of interacting behavioural plasticity expressed by both sexes. This highlights the need for new predictive theory informed by these empirically-derived parameters. Overall, we conclude that variation in the expression of shared traits due to interacting plasticity represents an important and novel facet of sexual interactions.

Impact SummaryAnimals and plants are able to respond to variation in their environment by expressing phenotypic plasticity. In sexual organisms, both males and females can exhibit such plasticity but the cues they respond to and the fitness consequences of these actions may be different between the sexes, and even conflicting. For example, males may respond to the presence of competitors by altering their mating behaviour or ejaculate transfer to increase their own, but not necessarily their mates reproductive output. However, females may also express phenotypic plasticity in response to their social and sexual environment to maximise their own fitness. Theory suggests that plasticity expressed by both sexes may either maximise the sex-specific fitness of both, or of one sex at the expense of the other. So far, little experimental work has been conducted to explore such interacting plasticity. Here we conducted novel tests of the fitness effects of interacting plasticity in the fruit fly Drosophila melanogaster. In doing so, we provide novel experimental evidence for interacting behavioural plasticity. We show that males and females can plastically respond to their socio-sexual environment to influence the expression of mating duration, mating latency, and fecundity. These plastic responses, while induced to increase the fitness interests of each sex, interact in the case of mating latency and fecundity and may reflect the outcome of sexual conflict. Our findings suggest that studies of reproductive behaviour should carefully consider the socio-sexual environment of both males and females and highlight the need for new predictive theory informed by empirically-derived parameters. Overall, we show that interacting plasticity between sexes represents an important and novel facet of sexual interactions.


title: Evidence for Multiple Teosinte Hybrid Zones in Central Mexico

abstract:
1O_LIHybrid zones provide an excellent opportunity for studying population dynamics and whether hybrid genetic architectures are locally adaptive. The genus Zea contains many diverse wild taxa collectively called teosinte. Zea mays ssp. parviglumis, the lowland progenitor of maize (Zea mays ssp. mays), and its highland relative Zea mays ssp. mexicana live parapatrically and, while putative hybrids have been identified in regions of range overlap, these have never been deeply explored.
C_LIO_LIHere we use a broadly sampled SNP data set to identify and confirm 112 hybrids between Zea mays ssp. parviglumis and Zea mays ssp. mexicana, mostly clustered in three genetically and geographically distinct hybrid groups in Central Mexico.
C_LIO_LIThese hybrid groups inhabit intermediate environments relative to parental taxa. We demonstrate that these individuals are true hybrids and not products of isolation by distance or ancestral to parviglumis and mexicana. This work expands on previous studies, clearly identifying hybrid zones in Zea, genetically characterizing hybrid groups, and showing what appear to be unique genetic architectures of hybridization in distinct hybrid groups.
C_LIO_LIWith the potential for local adaptation, variable hybrid zone dynamics, and differential architectures of hybridization, we present these teosinte hybrids and parental taxa as a promising model system for studying hybridization and hybrid zones.
C_LI


title: Learning Sparse Log-Ratios for High-Throughput Sequencing Data

abstract: The automatic discovery of sparse biomarkers that are associated with an outcome of interest is a central goal of bioinformatics. In the context of high-throughput sequencing (HTS) data, and compositional data (CoDa) more generally, an important class of biomarkers are the log-ratios between the input variables. However, identifying predictive log-ratio biomarkers from HTS data is a combinatorial optimization problem, which is computationally challenging. Existing methods are slow to run and scale poorly with the dimension of the input, which has limited their application to low- and moderate-dimensional metagenomic datasets. Building on recent advances from the field of deep learning, we present CoDaCoRe, a novel learning algorithm that identifies sparse, interpretable, and predictive log-ratio biomarkers. Our algorithm exploits a continuous relaxation to approximate the underlying combinatorial optimization problem. This relaxation can then be optimized efficiently using the modern ML toolbox, in particular, gradient descent. As a result, CoDaCoRe runs several orders of magnitude faster than competing methods, all while achieving state-of-the-art performance in terms of predictive accuracy and sparsity. We verify the outperformance of CoDaCoRe across a wide range of microbiome, metabolite, and microRNA benchmark datasets, as well as a particularly high-dimensional dataset that is outright computationally intractable for existing sparse log-ratio selection methods.1


title: BIAPSS - BioInformatic Analysis of liquid-liquid Phase-Separating protein Sequences

abstract:
Liquid-liquid phase separation (LLPS) has recently emerged as a cornerstone mechanism underlying the biogenesis of membraneless organelles (MLOs). However, a quantitative molecular grammar of protein sequences that controls the LLPS remains poorly understood. The progress in this field is hampered by the insufficiency of comprehensive databases and associated computational infrastructure for targeting biophysical and statistical analysis of phase separating biopolymers. Therefore, we have created a novel open-source web platform named BIAPSS (BioInformatic Analysis of liquid-liquid Phase-Separating protein Sequences) which contains interactive data analytic tools in combination with a comprehensive repository of bioinformatic data for on-the-fly exploration of sequence-dependent properties of proteins with known LLPS behavior. BIAPSS includes a residue-resolution biophysical analyzer for interrogating individual protein sequences (SingleSEQ tab). The latter allows users to correlate regions prone to phase separation with a large array of physicochemical attributes and various short linear motifs. BIAPSS also includes global statistics derived over the universe of most of the known LLPS-driver protein sequences (MultiSEQ tab) for revealing the regularities and sequence-specific signals driving phase separation. Finally, BIAPSS incorporates an extensive cross-reference section that links all entries to primary LLPS databases and other external resources thereby serving as a central navigation hub for the phase separation community. All of the data used by BIAPSS is freely available for download as well-formatted pre-processed data with detailed descriptions, facilitating rapid implementation in user-defined computational protocols.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=70 SRC="FIGDIR/small/430806v2_ufig1.gif" ALT="Figure 1">
View larger version (24K):
org.highwire.dtl.DTLVardef@1d90ddorg.highwire.dtl.DTLVardef@111d084org.highwire.dtl.DTLVardef@90b3b0org.highwire.dtl.DTLVardef@51d71_HPS_FORMAT_FIGEXP  M_FIG TOC - graphical abstract

C_FIG Author summaryProteins, especially those with low complexity and intrinsically disordered regions, have recently come into the limelight because of mounting evidence showing that these regions can drive the formation of membraneless organelles (MLOs) in cells. The underlying physical mechanism for forming MLOs is liquid-liquid phase separation (LLPS); a thermodynamically driven process whereby a cellular milieu with a relatively well-mixed distribution of biomolecules gets decomposed into liquid droplets where the concentration of selected biomolecules is higher. Deciphering molecular sequence grammar of phase separation has turned out to be challenging because of the complexity of this process in cells and the vastness of sequence space of LLPS-driver proteins. While the field is still in its infancy the growth of experimental data has already spurred the creation of several major databases which collect and annotate bimolecular systems with confirmed LLPS behavior. What is currently missing is a framework that would leverage the existing databases by integrating them with deep biophysical and bioinformatic analysis for identifying statistically significant features of protein sequences implicated in LLPS. In this work, we have addressed this challenge by creating an open-source web platform named BIAPSS (BioInformatic Analysis of liquid-liquid Phase-Separating protein Sequences) which integrates a comprehensive repository of pre-processed bioinformatic data for LLPS-driver protein sequences with interactive analytic applications for on-the-fly analysis of biophysical features relevant for LLPS behavior. BIAPSS empowers users with novel and effective tools for exploring LLPS-related sequence signals for individual proteins (SingleSEQ tab) and globally by integrating common regularities across subgroups or the entire LLPS sequence superset (MultiSEQ). The long-term plan for BIAPSS is to serve as a unifying hub for the experimental and computational community with a comprehensive set of analytic tools, biophysically featured data, and standardized protocols facilitating the identification of sequence hot spots driving the LLPS, which all can support applications for designing new sequences of biomedical interest.


title: Physical phenotype of blood cells is altered in COVID-19

abstract: Clinical syndrome coronavirus disease 2019 (COVID-19) induced by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is characterized by rapid spreading and high mortality worldwide. While the pathology is not yet fully understood, hyper-inflammatory response and coagulation disorders leading to congestions of microvessels are considered to be key drivers of the still increasing death toll. Until now, physical changes of blood cells have not been considered to play a role in COVID-19 related vascular occlusion and organ damage. Here we report an evaluation of multiple physical parameters including the mechanical features of five frequent blood cell types, namely erythrocytes, lymphocytes, monocytes, neutrophils, and eosinophils. More than 4 million blood cells of 17 COVID-19 patients at different levels of severity, 24 volunteers free from infectious or inflammatory diseases, and 14 recovered COVID-19 patients were analyzed. We found significant changes in erythrocyte deformability, lymphocyte stiffness, monocyte size, and neutrophil size and deformability. While some of these changes recovered to normal values after hospitalization, others persisted for months after hospital discharge, evidencing the long-term imprint of COVID-19 on the body.


title: Lack of effect of different pain-related manipulations on opioid self-administration, reinstatement of opioid seeking, and opioid choice in rats

abstract: Rationale and Objective: Pain-related factors increase risk for opioid addiction, and opioid-induced pain relief may function as a negative reinforcer to increase opioid taking and seeking. However, experimental pain-related manipulations generally do not increase opioid self-administration in rodents. This discrepancy may reflect insufficient learning of pain-relief contingencies or confounding effects of pain-related behavioral impairments. Here we determined if pairing noxious stimuli with opioid self-administration would promote pain-related reinstatement of opioid seeking or increase opioid choice over food. Methods: In Experiment 1, rats self-administered fentanyl in the presence or absence of repeated intraplantar capsaicin injections in distinct contexts to model context-specific exposure to cutaneous nociception. After capsaicin-free extinction in both contexts, we tested if capsaicin would reinstate fentanyl seeking. In Experiment 2, rats self-administered heroin after intraperitoneal (i.p.) lactic acid injections to model acute visceral inflammatory pain. After lactic acid-free extinction, we tested if lactic acid would reinstate heroin seeking. In Experiment 3, we tested if repeated i.p. lactic acid or intraplantar Complete Freund Adjuvant (CFA; to model sustained inflammatory pain) would increase fentanyl choice over food. Results: In Experiments 1-2, neither capsaicin nor lactic acid reinstated opioid seeking after extinction, and lactic acid did not increase heroin-induced reinstatement. In Experiment 3, lactic acid and CFA decreased reinforcement rate without affecting fentanyl choice. Conclusions: Results extend the range of conditions across which pain-related manipulations fail to increase opioid seeking in rats and suggest that enhanced opioid-addiction risk in humans with chronic pain involves factors other than enhanced opioid reinforcement and relapse.


title: Post mortem validation of DSI compressed sensing with optical imaging

abstract: While many useful microstructural indices, as well as orientation distribution functions, can be obtained from multi-shell dMRI data, there is growing interest in exploring the richer set of microstructural features that can be extracted from the full ensemble average propagator (EAP). The EAP can be readily computed from diffusion spectrum imaging (DSI) data, at the cost of a very lengthy acquisition. Compressed sensing (CS) has been used to make DSI more practical by reducing its acquisition time. CS applied to DSI (CS-DSI) attempts to reconstruct the EAP from significantly undersampled q-space data. We present a post mortem validation study where we evaluate the ability of CS-DSI to approximate not only fully sampled DSI but also multi-shell acquisitions with high fidelity. Human brain samples are imaged with high-resolution DSI at 9.4T and with polarization-sensitive optical coherence tomography (PSOCT). The latter provides direct measurements of axonal orientations at microscopic resolutions, allowing us to evaluate the mesoscopic orientation estimates obtained from diffusion MRI, in terms of their angular error and the presence of spurious peaks. We test two fast, dictionary-based, L2-regularized algorithms for CS-DSI reconstruction. We find that, for a CS acceleration factor of R=3, i.e., an acquisition with 171 gradient directions, one of these methods is able to achieve both low angular error and low number of spurious peaks. With a scan length similar to that of high angular resolution multi-shell acquisition schemes, this CS-DSI approach is able to approximate both fully sampled DSI and multi-shell data with high accuracy. Thus it is suitable for orientation reconstruction and microstructural modeling techniques that require either grid- or shell-based acquisitions. We find that the signal-to-noise ratio (SNR) of the training data used to construct the dictionary can have an impact on the accuracy of CS-DSI, but that there is substantial robustness to loss of SNR in the test data. Finally, we show that, as the CS acceleration factor increases beyond R=3, the accuracy of these reconstruction methods degrade, either in terms of the angular error, or in terms of the number of spurious peaks. Our results provide useful benchmarks for the future development of even more efficient q-space acceleration techniques.


title: High-fidelity approximation of grid- and shell-based sampling schemes from undersampled DSI using compressed sensing: Post mortem validation


title: miR-34a is not a candidate tumor suppressor

abstract: While several microRNAs (miRNAs) have been proposed to act as tumor suppressors, a consensual definition of tumor suppressing miRNAs is still missing. Similarly to coding genes, we propose that tumor suppressor miRNAs must show evidence of genetic or epigenetic inactivation in cancers, and exhibit an anti-tumorigenic (e.g., anti-proliferative) activity under endogenous expression levels. Here we observe that this definition excludes the most extensively studied tumor suppressor candidate miRNA, miR-34a. In analyzable cancer types, miR-34a does not appear to be down-regulated in primary tumors relatively to normal adjacent tissues. Deletion of miR-34a is occasionally found in human cancers, but it does not seem to be driven by an anti-tumorigenic activity of the miRNA, since it is not observed upon smaller, miR-34a-specific alterations. Its anti-proliferative action was observed upon large, supra-physiological transfection of synthetic miR-34a in cultured cells, and our data indicates that endogenous miR-34a levels do not have such an effect. Our results therefore argue against a general tumor suppressive function for miR-34a, providing an explanation to the lack of efficiency of synthetic miR-34a administration against solid tumors.

category: genetics


title: A rationalized definition of tumor suppressor microRNAsexcludes miR-34a


title: A rationalized definition of general tumor suppressor microRNAs excludes miR-34a


title: Dynamics and modulation of the human snoRNome

abstract:
BackgroundSmall nucleolar RNAs (snoRNAs) are mid-size non-coding RNAs required for ribosomal RNA modification, implying a ubiquitous tissue distribution linked to ribosome synthesis. However, increasing numbers of studies identify extra-ribosomal roles of snoRNAs in modulating gene expression, suggesting more complex snoRNA abundance patterns. Therefore, there is a great need for mapping the snoRNome in different human tissues as the blueprint for snoRNA functions.

ResultsWe used a low structure bias RNA-Seq approach to accurately quantify snoRNAs and compare them to the entire transcriptome in seven healthy human tissues (breast, ovary, prostate, testis, skeletal muscle, liver and brain). We identified 475 expressed snoRNAs categorized in two abundance classes that differ significantly in their function, conservation level and correlation with their host gene: 390 snoRNAs are uniformly expressed and 85 are enriched in the brain or reproductive tissues. Most tissue-enriched snoRNAs are embedded in lncRNAs and display strong correlation of abundance with them, whereas uniformly expressed snoRNAs are mostly embedded in protein-coding host genes and are mainly non- or anticorrelated with them. 59% of the non-correlated or anticorrelated protein-coding host gene/snoRNA pairs feature dual-initiation promoters, as opposed to only 16% of the correlated non-coding host gene/snoRNA pairs.

ConclusionsOur results demonstrate that snoRNAs are not a single homogeneous group of housekeeping genes but include highly regulated tissue-enriched RNAs. Indeed, our work indicates that the architecture of snoRNA host genes varies to uncouple the host and snoRNA expressions in order to meet the different snoRNA abundance levels and functional needs of human tissues.


title: Annotation of snoRNA abundance across human tissues reveals complex snoRNA-host gene relationships


title: Elucidating vimentin interaction with zinc ions and its interplay with oxidative modifications through crosslinking assays and molecular dynamics simulations

abstract: The intermediate filament protein vimentin is involved in essential cellular processes, including cell division and stress responses. Vimentin oxidative modifications impact network reorganization and its single cysteine residue, Cys328, acts as a redox sensor. Vimentin binds zinc, which influences its assembly by undefined mechanisms. Here, results from combined biochemical and molecular dynamics studies support that zinc ions interact with Cys328 in its thiolate form, whereas Glu329 and Asp331 stabilize zinc coordination. Vimentin oxidation can induce disulfide crosslinking, implying a close proximity of cysteine residues in certain vimentin associations, validated by our computational models. Notably, micromolar zinc concentrations selectively prevent Cys328 alkylation and crosslinking. These effects are not mimicked by magnesium, consistent with the fewer magnesium ions hosted at the cysteine region. Altogether, our results pinpoint the region surrounding Cys328, highly conserved in type III intermediate filaments, as a hot spot for zinc binding, which modulates Cys328 reactivity and vimentin assembly.

category: biochemistry


title: Selective inhibition reveals the regulatory function of DYRK2 in protein synthesis and calcium entry

abstract: The dual-specificity tyrosine phosphorylation-regulated kinase DYRK2 has emerged as a key regulator of cellular processes such as proteasome-mediated protein degradation. To gain further insights into its function, we took a chemical biology approach and developed C17, a potent small-molecule DYRK2 inhibitor, through multiple rounds of structure-based optimization guided by a number of co-crystallized structures. C17 displayed an effect on DYRK2 at a single-digit nanomolar IC50 and showed outstanding selectivity for the human kinome containing 467 other human kinases. Using C17 as a chemical probe, we further performed quantitative phosphoproteomic assays and identified several novel DYRK2 targets, including eukaryotic translation initiation factor 4E-binding protein 1 (4E-BP1) and stromal interaction molecule 1 (STIM1). DYRK2 phosphorylated 4E-BP1 at multiple sites, and the combined treatment of C17 with AKT and MEK inhibitors showed synergistic 4E-BP1 phosphorylation suppression. The phosphorylation of STIM1 by DYRK2 substantially increased the interaction of STIM1 with the ORAI1 channel, and C17 impeded the store-operated calcium entry process. Collectively, these studies further expand our understanding of DYRK2 and provide a valuable tool to further pinpoint its biological function.


title: Nse5/6 is a negative regulator of the ATPase activity of the Smc5/6 complex

abstract: The multi-component Smc5/6 complex plays a critical role in the resolution of recombination intermediates formed during mitosis and meiosis, and in the cellular response to replication stress. Using recombinant proteins, we have reconstituted a series of defined S. cerevisiae SMC5/6 complexes, visualised them by negative stain electron microscopy, and tested their ability to function as an ATPase. We find that only the six protein  holo-complex is capable of turning over ATP and that its activity is significantly increased by the addition of double-stranded DNA to reaction mixes. Furthermore, stimulation is wholly dependent on functional ATP-binding pockets in both Smc5 and Smc6. Importantly, we demonstrate that budding yeast Nse5/6 acts as a negative regulator of Smc5/6 ATPase activity, binding to the head-end of the complex to suppress turnover, irrespective of the DNA-bound status of the complex.


title: The bacterial iron sensor IdeR recognizes its DNA targets by indirect readout

abstract: IdeR is the main transcriptional regulator controlling iron homeostasis genes in Actinobacteria, including species from the Corynebacterium, Mycobacterium, and Streptomyces genera, as well as the erythromycin-producing bacterium Saccharopolyspora erythraea. Despite being a well-studied transcription factor since the identification of the Diphtheria toxin repressor DtxR three decades ago, the details of how IdeR proteins recognize their highly conserved 19-bp DNA target remain to be elucidated. The results of our structural and mutational studies support a model wherein IdeR uses an indirect readout mechanism, identifying its targets via the sequence-specific DNA backbone structure rather than through direct contacts with the DNA bases. Furthermore, we show that IdeR efficiently recognizes a shorter palindromic sequence corresponding to a half binding site as compared to the full 19-bp target previously reported, expanding the number of potential target genes controlled by IdeR proteins.


title: Local thermodynamics governs the formation and dissolution of protein condensates in living cells

abstract:
Membraneless compartments, also known as condensates, provide chemically distinct environments and thus spatially organize the cell. A well-studied example of condensates is P granules in the roundworm C. elegans which play an important role in the development of the germline. P granules are RNA-rich protein condensates that share the key properties of liquid droplets such as a spherical shape, the ability to fuse, and fast diffusion of their molecular components. An outstanding question is to what extent phase separation at thermodynamic equilibrium is appropriate to describe the formation of condensates in an active cellular environment. To address this question, we investigate the response of P granule condensates in living cells to temperature changes. We observe that P granules dissolve upon increasing the temperature and recondense upon lowering the temperature in a reversible manner. Strikingly, this temperature response can be captured by in vivo phase diagrams which are well described by a Flory-Huggins model at thermodynamic equilibrium. This finding is surprising due to active processes in a living cell. To address the impact of such active processes on intra-cellular phase separation, we discuss temperature heterogeneities. We show that, for typical estimates of the density of active processes, temperature represents a well-defined variable and that mesoscopic volume elements are at local thermodynamic equilibrium. Our findings provide strong evidence that P granule assembly and disassembly are governed by phase separation based on local thermal equilibria where the non-equilibrium nature of the cytoplasm is manifested on larger scales.

SIGNIFICANCE STATEMENTLiving cells rely on a continuous flux of energy to spatially organize biochemical processes. It remained unclear whether cells can achieve this spatial organization via thermodynamic principles. Here, we report the striking behavior of a cold-blooded organism that reacts to environmental temperature changes similar to a thermodynamic system at local equilibrium. Our key finding is that protein-rich droplets form and dissolve reversibly with temperature due to changes in the organism?s entropy. We show that the organism uses a specific molecule to extend droplet stability to the natural temperature range of the organisms habitat. Due to the relevance of such protein droplets for the organism?s fertility, our works shed light on how molecular components could facilitate biological functions via thermodynamic principles.


title: Development of an exosomal biomarker signature to detect minimal residual disease in dogs with osteosarcoma using a novel xenograft platform and machine learning

abstract:
Osteosarcoma has a guarded prognosis. A major hurdle in developing more effective osteosarcoma therapies is the lack of disease-specific biomarkers to predict risk, prognosis, or therapeutic response. Exosomes are secreted extracellular microvesicles emerging as powerful diagnostic tools. However, their clinical application is precluded by challenges in identifying disease-associated cargo from the vastly larger background of normal exosome cargo. We developed a method using canine osteosarcoma in mouse xenografts to distinguish tumor-derived from host-response exosomal mRNAs. The model allows for the identification of canine osteosarcoma-specific gene signatures by RNA sequencing and a species-differentiating bioinformatics pipeline. An osteosarcoma-associated signature consisting of five gene transcripts (SKA2, NEU1, PAF1, PSMG2, and NOB1) was validated in dogs with spontaneous osteosarcoma by qRT-PCR, while a machine learning model assigned dogs into healthy or disease groups. Serum/plasma exosomes were isolated from 53 dogs in distinct clinical groups ("healthy", "osteosarcoma", "other bone tumor", or "non-neoplastic disease"). Pre-treatment samples from osteosarcoma cases were used as the training set and a validation set from post-treatment samples was used for testing, classifying as "osteosarcoma-detected" or "osteosarcoma-NOT detected". Dogs in a validation set whose post-treatment samples were classified as "osteosarcoma-NOT detected" had longer remissions, up to 15 months after treatment. In conclusion, we identified a gene signature predictive of molecular remissions with potential applications in the early detection and minimal residual disease settings. These results provide proof-of-concept for our discovery platform and its utilization in future studies to inform cancer risk, diagnosis, prognosis, and therapeutic response.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=147 SRC="FIGDIR/small/429432v2_ufig1.gif" ALT="Figure 1">
View larger version (41K):
org.highwire.dtl.DTLVardef@1b612c4org.highwire.dtl.DTLVardef@8a1847org.highwire.dtl.DTLVardef@3ffd64org.highwire.dtl.DTLVardef@114df8c_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: cancer biology


title: Development of an exosomal gene signature to detect residual disease in dogs with osteosarcoma using a novel xenograft platform and machine learning


title: Identification of phenotype-specific networks from paired gene expression-cell shape imaging data

abstract: The morphology of breast cancer cells is often used as an indicator of tumour severity and prognosis. Additionally, morphology can be used to identify more fine-grained, molecular developments within a cancer cell, such as transcriptomic changes and signaling pathway activity. Delineating the interface between morphology and signaling is important to understand the mechanical cues that a cell processes in order to undergo epithelial-to-mesenchymal transition and consequently metastasize. However, the exact regulatory systems that define these changes remain poorly characterised. In this study, we employ a network-systems approach to integrate imaging data and RNA-seq expression data. Our workflow allows the discovery of unbiased and context-specific gene expression signatures and cell signaling sub-networks relevant to the regulation of cell shape, rather than focusing on the identification of previously known, but not always representative, pathways. By constructing a cell-shape signaling network from shape-correlated gene expression modules and their upstream regulators, we found central roles for developmental pathways such as WNT and Notch as well as evidence for the fine control of NFkB signaling by numerous kinase and transcriptional regulators. Further analysis of our network implicates a gene expression module enriched in the RAP1 signaling pathway as a mediator between the sensing of mechanical stimuli and regulation of NFkB activity, with specific relevance to cell shape in breast cancer.


title: The genetic basis of adaptation through the evolution of mating systems

abstract:
Although adaptation can be realized through the fixation of beneficial alleles that increase viability, many plant populations may adapt through the evolution of self-fertilization, especially when pollination becomes inefficient. However, the genetic basis of adaptation through the evolution of selfing remains unclear. Using population genetic models, I study adaptation through the fixation of alleles that increase the selfing rate (selfing modifiers) from new mutations or/and standing variation. For adaptive alleles unrelated to selfing, it is known that selfing promotes adaptation from a new mutation only when the beneficial alleles are recessive, and the probability of adaptation from standing variation is nearly independent of dominance, and always decreases with the selfing rate. In contrast, for adaptation through the evolution of selfing, when it occurs by fixation of a newly arisen mutation, a population that already has a high selfing rate may be more (less) likely to adapt than outcrossers even when the modifier is dominant (recessive) if the modifier is weakly (strongly) selected. Also, adaptation from standing variation is more likely through recessive modifier alleles, with the highest fixation probability found in partially selfing populations, but fixation is fastest when dominance is intermediate. When there are multiple modifiers, adaptation through new mutations is more likely when selfing is controlled by few large-effect rather than many slight-effect modifiers. This study suggests that to understand the genetic basis of adaptation, it is necessary to determine the ecological and genetic advantages of adaptive alleles.

Significance statementThis study, by deriving the selective coefficient and effective population size, investigated the genetic basis of adaptation through fixation of modifier alleles that increase the selfing rate, which is shown to differ in several aspects from that through evolution of mating-unrelated alleles. Specifically, when adaptation is from new mutations, the dominance of a selfing modifier allele below which selfing increases the fixation probability depends on the strength of pollen limitation and pollen discounting. Adaptation from standing variation is more likely through recessive modifier alleles and in populations with an intermediate selfing rate. This work suggests it is important to have a mechanistic understanding of how adaptive alleles increase individual fitness in environment that is unfavorable to the population.


title: Extensive introgression among North American wild grapes (Vitis) fuels biotic and abiotic adaptation

abstract: Introgressive hybridization can introduce adaptive genetic variation into a species or population. To evaluate the evolutionary forces that contribute to introgression, we studied six Vitis species that are native to the Southwestern United States and potentially useful for breeding grapevine (V. vinifera) rootstocks. By creating a reference genome from one wild species, V. arizonica, and by resequencing 130 accessions, we focused on identifying putatively introgressed regions (pIRs) between species. We found that up to ~8% of extant genome is attributable to introgression between species. The pIRs tended to be gene poor, located in regions of high recombination and enriched for genes implicated in disease resistance functions. To assess potential pIR function, we explored SNP associations to bioclimatic variables and to bacterial levels after infection with the causative agent of Pierces Disease. pIRs were enriched for SNPs associated with both climate and bacterial levels, suggesting potential drivers of adaptive events. Altogether, this study yields insights into the genomic extent of introgression, potential pressures that shape adaptive introgression, and the history of economically important wild relatives of a critical crop.


title: Alternative female and male developmental strategies in the dynamic balance of human visual perception

abstract: The numerous multistable phenomena in vision, hearing and touch attest that the inner workings of perception are prone to instability. We investigated a visual example - binocular rivalry - with an accurate no-report paradigm, and uncovered developmental and maturational lifespan trajectories that were specific for age and sex. To interpret these trajectories, we hypothesized that conflicting objectives of visual perception - such as stability of appearance, sensitivity to visual detail, and exploration of fundamental alternatives - change in relative importance over the lifespan. Computational modelling of our empirical results allowed us to estimate this putative development of stability, sensitivity, and exploration over the lifespan. Our results confirmed prior findings of developmental psychology and appear to quantify important aspects of neurocognitive phenotype. Additionally, we report atypical function of binocular rivalry in autism spectrum disorder and borderline personality disorder. Our computational approach offers new ways of quantifying neurocognitive phenotypes both in development and in dysfunction.


title: Alternative female and male developmental trajectories in the dynamic balance of human visual perception


title: Extensive N4 cytosine methylation is essential for Marchantia sperm function

abstract: 4-methylcytosine (4mC) is an important DNA modification in prokaryotes, but its relevance, and even presence in eukaryotes have been mysterious. Here we show that spermatogenesis in the liverwort Marchantia polymorpha involves two waves of extensive DNA methylation reprogramming. First, 5-methylcytosine (5mC), a well-known eukaryotic DNA modification, expands from transposons to the entire genome. Notably, the second wave installs 4mC throughout genic regions, covering over 50% of CG sites in sperm. 4mC is catalyzed by two novel methyltransferases (MpDN4MT1a and MpDN4MT1b) specifically expressed during late spermiogenesis. Deletion of MpDN4MT1s eliminates 4mC, alters the sperm transcriptome, and produces sperm with swimming defects. Our results reveal extensive 4mC in a eukaryote and define a new family of eukaryotic methyltransferases, thereby expanding the repertoire of functional eukaryotic DNA modifications.

category: developmental biology


title: The making of cauliflowers: the story of unsuccessful flowers

abstract:
The arrangement of plant organs, called phyllotaxis, produce remarkable spiral or whorled patterns. Cauliflowers present a unique phyllotaxis with a multitude of spirals over a wide range of scales. How such a self-similar fractal organization emerges from developmental mechanisms has remained elusive. Combining experimental assays with modeling, we found that cauliflowers arise due to the hysteresis of the bistable floral network that generates inflorescences imprinted by a transient floral state. We further show how additional mutations affecting meristem growth dynamics can induce the production of conical phyllotactic structures reminiscent of the conspicuous fractal Romanesco shape. This study reveals how the spectacular morphological modification of the inflorescences in cauliflower and Romanesco shape arises from the hysteresis of the genetic programs controlling inflorescence development.

One Sentence SummaryThe molecular making of cauliflowers

category: plant biology


title: The role of auxin and sugar signaling in dominance inhibition of inflorescence growth by fruit load

abstract:
Dominance inhibition of shoot growth by fruit load is a major factor that regulates shoot architecture and limits yield in agriculture and horticulture crops. In annual plants, the inhibition of inflorescence growth by fruit load occurs at a late stage of inflorescence development termed the end of flowering transition. Physiological studies show that this transition is mediated by production and export of auxin from developing fruits in close proximity to the inflorescence apex. In the meristem, cessation of inflorescence growth is controlled in part by the age dependent pathway, which regulates the timing of arrest. Here, results show that the end of flowering transition is a two-step process in which the first stage is characterized by a cessation of inflorescence growth, while immature fruit continue to develop. At this stage, dominance inhibition of inflorescence growth by fruit load correlates with a selective dampening of auxin transport in the apical region of the stem. Subsequently, an increase in auxin response in the vascular tissues of the apical stem where developing fruits are attached marks the second stage for the end of flowering transition. Similar to the vegetative and floral transition, the end of flowering transition correlates with a change in sugar signaling and metabolism in the inflorescence apex. Taken together, our results suggest that during the end of flowering transition, dominance inhibition of inflorescence shoot growth by fruit load is mediated by auxin and sugar signaling.

One-sentence summaryDominance inhibition of inflorescence shoot growth by fruit load is involves auxin and sugar signaling during the end of flowering transition.


title: An improved ChEC-seq method accurately maps the genome-wide binding of transcription coactivators and sequence-specific transcription factors

abstract: Mittal and colleagues have raised questions about mapping transcription factor locations on DNA using the MNase-based ChEC-seq method (Mittal et al., 2021). Partly due to this concern, we modified the experimental conditions of the MNase cleavage step and subsequent computational analyses, resulting in more stringent conditions for mapping protein-DNA interactions (Donczew et al., 2020). The revised method (dx.doi.org/10.17504/protocols.io.bizgkf3w) answers questions raised by Mittal et al. and, without changing earlier conclusions, identified widespread promoter binding of the transcription coactivators TFIID and SAGA at active genes. The revised method is also suitable for accurately mapping the genome-wide locations of DNA sequence-specific transcription factors.


title: Integration of electric, calcium, reactive oxygen species and hydraulic signals during rapid systemic signaling in plants

abstract:
The sensing of abiotic stress, mechanical injury, or pathogen attack by a single plant tissue results in the activation of systemic signals that travel from the affected tissue to the entire plant, alerting it of an impending stress or pathogen attack. This process is essential for plant survival during stress and is termed systemic signaling. Among the different signals triggered during this process are calcium, electric, reactive oxygen species (ROS) and hydraulic signals. These are thought to propagate at rapid rates through the plant vascular bundles and to regulate many of the systemic processes essential for plant survival. Although the different signals activated during systemic signaling are thought to be interlinked, their coordination and hierarchy remain to be determined. Here, using a combination of advanced whole-plant imaging and hydraulic pressure measurements, we studied the activation of all four systemic signals in wild type and different Arabidopsis thaliana mutants subjected to a local high light (HL) stress or wounding. Our findings reveal that in response to wounding systemic changes in membrane potential, calcium, ROS, and hydraulic pressure are coordinated by glutamate receptor-like (GLR) proteins 3.3 and 3.6, while in response to HL the respiratory burst oxidase homolog D-driven systemic ROS signal could be separated from systemic changes in membrane potential and calcium levels. We further determine that plasmodesmata functions are required for systemic changes in membrane potential, calcium, and ROS during systemic signaling. Our findings shed new light on the different mechanisms that integrate different systemic signals in plants during stress.

Significance statementThe ability of plants to transmit a signal from a stressed or wounded tissue to the entire plant, termed systemic signaling, is key to plant survival during conditions of environmental stress. At least four different systemic signals are thought to be involved in this process: electric, calcium, reactive oxygen and hydraulic. However, how are they coordinated and whether they can be stress-specific is mostly unknown. Here we report that different types of stimuli can induce different types of systemic signals that may or may not be linked with each other. We further reveal that hydraulic waves can be actively regulated in plants in response to wounding, and that proteins that regulate plasmodesmata pores play a key role in systemic signaling.


title: Musicianship and melodic predictability enhance neural gain in auditory cortex during pitch deviance detection

abstract:
When listening to music, pitch deviations are more salient and elicit stronger prediction error responses when the melodic context is predictable and when the listener is a musician. Yet, the neuronal dynamics and changes in synaptic efficacy underlying such effects remain unclear. Here, we employed dynamic causal modeling (DCM) to investigate whether the magnetic mismatch negativity response (MMNm)--and its modulation by context predictability and musical expertise--are associated with enhanced neural gain of auditory areas, as a plausible mechanism for encoding precision-weighted prediction errors. Using Bayesian model comparison, we asked whether models with intrinsic connections within primary auditory cortex (A1) and superior temporal gyrus (STG)--typically related to gain control--or extrinsic connections between A1 and STG--typically related to propagation of prediction and error signals--better explained magnetoencephalography (MEG) responses. We found that, compared to regular sounds, out-of-tune pitch deviations were associated with lower intrinsic (inhibitory) connectivity in A1 and STG, and lower backward (inhibitory) connectivity from STG to A1, consistent with disinhibition and enhanced neural gain in these auditory areas. More predictable melodies were associated with disinhibition in right A1, while musicianship was associated with disinhibition in left A1 and reduced connectivity from STG to left A1. These results indicate that musicianship and melodic predictability, as well as pitch deviations themselves, enhance neural gain in auditory cortex during deviance detection. Our findings are consistent with predictive processing theories suggesting that precise and informative error signals are selected by the brain for subsequent hierarchical processing.

Significance statementIn complex auditory contexts, being able to identify informative signals is of paramount importance. Such is the case of music listening, where surprising sounds play a fundamental role in its perceptual, aesthetical, and emotional experience. Crucially, surprising sounds in the pitch dimension are more easily detected and generate stronger cortical responses when melodies are predictable and when the listener is a musician. Using Dynamic Causal Modelling, here we show that such effects arise from a local increase in neural gain within auditory areas, rather than from passing of prediction and error signals between brain regions. Consistent with predictive processing theories, this suggests that the enhanced precision of auditory predictive models--through melodic predictability and musical training--up-regulates the processing of informative error signals in the brain.


title: Abrupt remapping in human CA3/dentate gyrus signals resolution of memory interference

abstract: Remapping refers to a decorrelation of hippocampal representations of similar spatial environments. While it has been speculated that remapping may contribute to the resolution of episodic memory interference in humans, direct evidence is surprisingly limited. Here, we tested this idea using high-resolution, pattern-based fMRI analyses. We show that activity patterns in human CA3/dentate gyrus exhibit an abrupt, temporally-specific decorrelation of highly similar memory representations that is precisely coupled with behavioral expressions of successful learning. Strikingly, the magnitude of this learning-related decorrelation was predicted by the amount of pattern overlap during initial stages of learning, with greater initial overlap leading to stronger decorrelation. Finally, we show that remapped activity patterns carry relatively more information about learned episodic associations compared to competing associations, further validating the learning-related significance of remapping. Collectively, these findings establish a critical link between hippocampal remapping and episodic memory interference and provide novel insight into why remapping occurs.


title: Postpartum corticosterone and fluoxetine shift the tryptophan-kynurenine pathway in dams

abstract:
Perinatal depression (PND) affects 15% of mothers. Selective serotonin reuptake inhibitors (SSRIs) are currently the first-line of treatment for PND but are not always efficacious. Previously, we found significant reductions in plasma tryptophan concentrations and higher hippocampal proinflammatory cytokine, IL-1b levels, due to maternal SSRI treatment. Both inflammation and tryptophan-kynurenine metabolic pathway (TKP) are associated with SSRI efficacy in individuals with major depressive disorder (MDD). TKP is divided into neuroprotective and neurotoxic pathways. Higher metabolite concentrations of the neurotoxic pathway are associated with depression onset and implicated in SSRI efficacy. Metabolites in TKP were investigated in a rodent model of de novo postpartum depression (PPD) given treatment with the SSRI, fluoxetine (FLX). Dams were administered corticosterone (CORT) (40mg/kg, s.c.), and treated with the SSRI, fluoxetine (FLX) (10mg/kg, s.c.), during the postpartum for 22 days after parturition. Plasma TKP metabolite concentrations were quantified on the last day of treatment. Maternal postpartum CORT increased neurotoxic metabolites and co-enzyme/cofactors in dams (3-hydroxykynurenine, 3-hydroxyanthranilic acid, vitamin B2, flavin adenine dinucleotide). The combination of both CORT and FLX shifted the neuroprotective-to-neurotoxic ratio towards neurotoxicity. Postpartum FLX decreased plasma xanthurenic acid concentrations. Together, our data indicate higher neurotoxic TKP expression due to maternal postpartum CORT treatment, similar to clinical presentation of MDD. Moreover, maternal FLX treatment showed limited efficacy to influence TKP metabolites, which may correspond to its limited efficacy to treat depressive-like endophenotypes. Overall suggesting changes in TKP may be used as a biomarker of de novo PPD and antidepressant efficacy and targeting this pathway may serve as a potential therapeutic target.

HighlightsO_LITryptophan-kynurenine pathway (TKP) is altered by postpartum corticosterone (CORT)
C_LIO_LIPostpartum CORT upregulated neurotoxic more metabolites (3HK, 3HAA)
C_LIO_LIPostpartum fluoxetine (FLX) increased xanthurenic acid concentrations
C_LIO_LIPostpartum CORT and FLX together shifted the TKP balance towards neurotoxicity
C_LI


title: Phenotypic analysis of catastrophic childhood epilepsy genes: The Epilepsy Zebrafish Project

abstract: Genetic engineering techniques have contributed to the now widespread use of zebrafish to investigate gene function, but zebrafish-based human disease studies, and particularly for neurological disorders, are limited. Here we used CRISPR-Cas9 to generate 40 single-gene mutant zebrafish lines representing catastrophic childhood epilepsies. We evaluated larval phenotypes using electrophysiological, behavioral, neuro-anatomical, survival and pharmacological assays. Phenotypes with unprovoked electrographic seizure activity (i.e., epilepsy) were identified in zebrafish lines for 8 genes; ARX, EEF1A, GABRB3, GRIN1, PNPO, SCN1A, STRADA and STXBP1. A unifying epilepsy classification scheme was developed based on local field potential recordings and blinded scoring from ~3300 larvae. We also created an open-source database containing sequencing information, survival curves, behavioral profiles and representative electrophysiology data. We offer all zebrafish lines as a resource to the neuroscience community and envision them as a starting point for further functional analysis and/or identification of new therapies.


title: Parameter estimation for correlated Ornstein-Uhlenbeck time-series

abstract: In this article, we develop a Maximum likelihood (ML) approach to estimate parameters from correlated time traces that originate from coupled Ornstein-Uhlenbeck processes. The most common technique to characterize the correlation between time-series is to calculate the Pearson correlation coefficient. Here we show that for time series with memory (or a characteristic relaxation time), our method gives more reliable results, but also results in coupling coefficients and their uncertainties given the data. We investigate how these uncertainties depend on the number of samples, the relaxation times and sampling time. To validate our analytic results, we performed simulations over a wide range of correlation coefficients both using our maximum likelihood solutions and Markov-Chain Monte-Carlo (MCMC) simulations. We found that both ML and MCMC result in the same parameter estimations. We also found that when analyzing the same data, the ML and MCMC uncertainties are strongly correlated, while ML underestimates the uncertainties by a factor of 1.5 to 3 over a large range of parameters. For large datasets, we can therfore use the less computationally expensive maximum likelihood method to run over the whole dataset, and then we can use MCMC on a few samples to determine the factor by which the ML method underestimates the uncertainties. To illustrate the application of our method, we apply it to time series of brain activation using fMRI measurements of the human default mode network. We show that our method significantly improves the interpretation of multi-subject measurements of correlations between brain regions by providing parameter confidence intervals for individual measurements, which allows for distinguishing between the variance from differences between subjects from variance due to measurement error.


title: The C. elegans Notch proteins LIN-12 and GLP-1 are tuned to lower force thresholds for activation than Drosophila Notch

abstract: The conserved transmembrane receptor Notch mediates cell fate decisions in all animals. In the absence of ligand, a Negative Regulatory Region (NRR) in the Notch ectodomain adopts an autoinhibited confirmation, masking an ADAM protease cleavage site [1, 2]; ligand binding makes the cleavage site accessible, leading to shedding of the Notch ectodomain as the first step of signal transduction [3, 4]. In Drosophila and vertebrates, the ligands are all single-pass transmembrane Delta/Serrate/LAG-2 (DSL) proteins; the endocytic adaptor Epsin binds to the ubiquitinated intracellular domain, and the resulting Clathrin-mediated endocytosis exerts a "pulling force" that exposes the cleavage site in the NRR [4-6]. However, in C. elegans, the presence of natural secreted DSL proteins [7] and other observations suggested that Epsin-mediated endocytosis may not be required to activate the Notch proteins LIN-12 and GLP-1. Here, we confirm that neither Epsin nor the cytosolic domains of DSL proteins are required for Notch signaling in C. elegans. Furthermore, we provide evidence that the NRRs of LIN-12 and GLP-1 are tuned to a lower force level than the NRR of Drosophila Notch. Finally, we show that adding a Leucine "plug" that occludes the cleavage site in vertebrate and Drosophila Notch proteins but is absent in the C. elegans Notch proteins [1, 2] renders the LIN-12 and GLP-1 NRRs dependent on Epsin-mediated ligand endocytosis, indicating that greater force is now required to expose the cleavage site. Thus, the NRRs of LIN-12 and GLP-1 appear to be tuned to a lower force threshold, accounting for the different requirements for signaling in C. elegans.


title: Different trajectories of polyploidization shape the genomic landscape of the Brettanomyces bruxellensis yeast species

abstract: Polyploidization events are observed across the tree of life and occurred in many fungi, plant and animal species. Polyploidy is thought to be an important source of speciation and tumorigenesis. However, the origins of polyploid populations are not always clear and little is known about the precise nature and structure of their complex genome. Using a long-read sequencing strategy, we sequenced a large number of isolates from the Brettanomyces bruxellensis yeast species, which is found in anthropized environments (e.g. beer, contaminant of wine, kombucha and ethanol production) and characterized by several polyploid subpopulations. To reconstruct the polyploid genomes, we phased them by using different strategies and we found that each subpopulation had a unique polyploidization history with distinct trajectories. The polyploid genomes contain either genetically closely related (with a genetic divergence < 1%) or diverged copies (> 3%), indicating auto- as well as allopolyploidization events. These latest events have occurred independently with a specific and unique donor in each of the polyploid subpopulations, and exclude the known Brettanomyces sister species as possible donors. Finally, loss of heterozygosity events have shaped the structure of these polyploid genomes and underline their dynamic. Overall, our study highlights the multiplicity of the trajectories leading to polyploid genomes within a same species.


title: Cancer cells induce hepatocytes apoptosis in co-opted colorectal cancer liver metastatic lesions

abstract: Vessel co-option in colorectal cancer liver metastases (CRCLM) has been recognized as one of the mechanistic pathways that contribute to resistance against anti-angiogenic therapy. In vessel co-opted CRCLM lesions, the cancer cells are highly motile that move toward and along the pre-existing sinusoidal vessels and hijack them to gain access to nutrient. The movement of cancer cells is accompanied by replacement of the hepatocytes. However, the molecular mechanisms by which this replacement occurs are unclear yet. To examine the involvement of apoptosis in hepatocytes replacement by cancer cells in co-opted lesions, we conducted immunohistochemical staining for chemonaive CRCLM specimens using pro-apoptotic markers antibody, such as cleaved caspase-3 and cleaved poly (ADP-ribose) polymerase-1 (PARP-1). The results suggested overexpression of pro-apoptotic markers in liver parenchyma of co-opted lesions compared to angiogenic lesions, specifically the hepatocytes that are in close proximity to the cancer cells. Importantly, co-culturing hepatocytes with colorectal cancer cells induced overexpression of pro-apoptotic markers in the hepatocytes. Altogether, these results propose that cancer cells could exploit apoptosis to replace the hepatocytes and establish vessel co-option in CRCLM.


title: Beyond massive univariate tests: Covariance regression reveals complex patterns of functional connectivity related to attention-deficit/hyperactivity disorder, age, sex, and response control

abstract: We applied a novel Covariate Assisted Principal (CAP) whole-matrix regression approach to identify resting-state functional connectivity (FC) brain networks associated with attention-deficit/hyperactivity disorder (ADHD) and response control. Participants included 8-12 year-old children with ADHD (n=115, 29 girls) and typically developing controls (n=102, 35 girls) with a resting-state fMRI scan and go/no-go task behavioral data. We modeled three sets of covariates to identify resting-state networks associated with ADHD, age, sex, and response control. Four networks were identified across models revealing complex interactions between subregions of cognitive control, default mode, subcortical, visual, and somatomotor networks that relate to age, response control, and a diagnosis of ADHD among girls and boys. Unique networks were also identified in each of the three models suggesting some specificity to the covariates of interest. These findings demonstrate the utility of our novel covariance regression approach to studying functional brain networks relevant for development, behavior, and psychopathology.


title: Get rid of the beat in mobile EEG applications: A framework towards automated cardiogenic artifact detection and removal in single-channel EEG

abstract: Brain activity recordings outside clinical or laboratory settings using mobile EEG systems have recently gained popular interest allowing for realistic long-term monitoring and eventually leading to identification of possible biomarkers for diseases. The less obtrusive, minimized systems (e.g. single-channel EEG, no ECG reference) have the drawback of artifact contamination with varying intensity that are particularly difficult to identify and remove. We developed brMEGA, the first algorithm for automated detection and removal of cardiogenic artifacts using non-linear time-frequency analysis and machine learning to (1) detect whether and where cardiogenic artifacts exist, and (2) remove those artifacts. We compare our algorithm against visual artifact identification and a previously established approach and validate it in one real and semi-real datasets. We demonstrated that brMEGA successfully identifies and substantially removes cardiogenic artifacts in single-channel EEG recordings. Moreover, recovery of cardiogenic artifacts gives the opportunity for future extraction of heart rate features without ECG measurement.


title: Schizophrenia-linked protein tSNARE1 regulates endolysosomal trafficking in cortical neurons

abstract: TSNARE1, which encodes the protein tSNARE1, is a high-confidence gene candidate for schizophrenia risk, but nothing is known about its cellular or physiological function. We identified the major gene products of TSNARE1 and their cytoplasmic localization and function in endolysosomal trafficking in murine cortical neurons. We validated four primary isoforms of TSNARE1 expressed in human brain, all of which encode a syntaxin-like Qa SNARE domain. RNA-sequencing data from adult and fetal human brain suggested that the majority of tSNARE1 lacks a transmembrane domain that is thought to be necessary for membrane fusion. Live-cell imaging demonstrated that brain tSNARE1 isoforms localized to compartments of the endosomal network. The most abundant brain isoform, tSNARE1c, localized most frequently to Rab7+ late endosomal compartments. Expression of either tSNARE1b or tSNARE1c, which differ only in their inclusion or exclusion of a Myb-like domain, delayed the trafficking of known endosomal cargo, Neep21, into late endosomal and lysosomal compartments. These data suggest that tSNARE1 regulates endosomal trafficking in cortical neurons, likely through negatively regulating endocytic trafficking or maturation of early endosomes to late endosomes.


title: Unexpected organellar locations of ESCRT machinery in Giardia intestinalis and complex evolutionary dynamics spanning the transition to parasitism in the lineage Fornicata

abstract:
Comparing a parasitic lineage to its free-living relatives is a powerful way to understand how the evolutionary transition to parasitism occurred. Giardia intestinalis (Fornicata) is a leading cause of gastrointestinal disease world-wide and is famous for its unusual complement of cellular compartments, such as having peripheral vacuoles instead of typical endosomal compartments. Endocytosis plays an important role in Giardias pathogenesis. Endosomal sorting complexes required for transport (ESCRT) are membrane-deforming proteins associated with the late endosome/multivesicular body (MVB). MVBs are ill-defined in G. intestinalis and roles for identified ESCRT-related proteins are not fully understood in the context of its unique endocytic system. Furthermore, components thought to be required for full ESCRT functionality have not yet been documented in this species.

We used genomic and transcriptomic data from several Fornicata species to clarify the evolutionary genome streamlining observed in Giardia, as well as to detect any divergent orthologs of the Fornicata ESCRT subunits. We observed differences in the ESCRT machinery complement between Giardia strains. Microscopy-based investigations of key components of ESCRT machinery such as GiVPS36and GiVPS25 link them to peripheral vacuoles, highlighting these organelles as simplified MVB equivalents. Unexpectedly, we show ESCRT components associated with the Endoplasmic Reticulum, and for the first time, mitosomes. Finally, we identified the rare ESCRT component CHMP7 in several fornicate representatives, including Giardia, and show that contrary to current understanding, CHMP7 evolved from a gene fusion of VPS25 and SNF7 domains, prior to the last eukaryotic common ancestor, over 1.5 billion years ago. Our findings show that ESCRT machinery in G. intestinalis is far more varied and complete than previously thought, and associating to multiple cellular locations and presenting changes in ESCRT complement which pre-date adoption of a parasitic lifestyle.


title: The clinical and functional effects of TERT variants in myelodysplastic syndrome

abstract: Germline pathogenic TERT variants are associated with short telomeres and an increased risk of developing myelodysplastic syndrome (MDS) among patients with a telomere biology disorder. We identified TERT rare variants in 41 of 1514 MDS patients (2.7%) without a clinical diagnosis of telomere biology disorder who underwent allogeneic transplantation. Patients with TERT rare variants had shorter telomere length (p<0.001) and younger age at MDS diagnosis (52 vs. 59 years, p=0.03) than patients without a TERT rare variant. In multivariable analyses, TERT rare variants were associated with inferior overall survival (p=0.034) driven by an increased incidence of non-relapse mortality (NRM) (p=0.015). Death from a non-infectious pulmonary cause was more frequent among patients with a TERT rare variant. According to ACMG/AMP guidelines and Sherloc criteria, 39 TERT rare variants were classified as VUS and one as likely pathogenic. Therefore, we cloned all rare missense variants and quantified their impact on telomere elongation in a cell-based assay. We found that 36 of 40 variants had severe or intermediate impairment in their capacity to elongate telomeres. Using a homology model of human TERT bound to the shelterin protein TPP1, we inferred that TERT rare variants disrupt domain-specific functions, including catalysis, protein-RNA interactions, and recruitment to telomeres. Our results indicate that the contribution of TERT rare variants to MDS pathogenesis and NRM risk is underrecognized and routine screening for TERT rare variants in MDS patients regardless of age or clinical suspicion could identify clinically inapparent telomere biology disorders and improve transplant outcomes through risk-adapted approaches.


title: Understanding conformational dynamics from macromolecular crystal diffuse scattering

abstract: All macromolecular crystals contain some extent of disorder. The diffraction from such crystals contains diffuse scattering in addition to Bragg peaks and this scattering contains information about correlated displacements in the constituent molecules. While much work has been performed recently in decoding the dynamics of the crystalline ordering, the goal of understanding the internal dynamics of the molecules within a unit cell has been out-of-reach. In this article, we propose a general framework to extract the internal conformational modes of a macromolecule from diffuse scattering data. We combine insights on the distribution of diffuse scattering from short- and long-range disorder with a Bayesian global optimization algorithm to obtain the best fitting internal motion modes to the data. To illustrate the efficacy of the method, we apply it to a publicly available dataset from triclinic lysozyme. Our mostly parameter-free approach can enable the recovery of a much richer, dynamic structure from macromolecular crystallography.


title: Targeted enzyme assisted chemotherapy (TEAC) - a novel microRNA-guided and selenium-based regimen to specifically eradicate hepatocellular carcinoma

abstract: Despite progress in the treatment of non-visceral malignancies, the prognosis remains poor for malignancies of visceral organs and novel therapeutic approaches are urgently required. Here we introduce a novel therapeutic regimen by treatment with Se-methylselenocysteine (MSC) and concomitant tumor-specific induction of Kynurenine aminotransferase 1 (KYAT1) in hepatocellular carcinoma (HCC) cell lines, using either vector-based and/or lipid nanoparticle-mediated delivery of mRNA. Supplementation of MSC in KYAT1 overexpressed cells resulted in significantly increased cytotoxicity as compared to MSC alone. Furthermore, microRNA antisense targeted sites for miR122, known to be widely expressed in normal hepatocytes whilst downregulated in hepatocellular carcinoma, were added to specifically limit cytotoxicity in HCC cells, thereby limiting off-target effects. KYAT1 expression was significantly reduced in cells with high levels of miR122 supporting the concept of miR-guided induction of tumor-specific cytotoxicity. The addition of alpha-ketoacid favored the production of methylselenol, enhancing the cytotoxic efficacy of MSC in HCC cells, with no effects on primary human hepatocytes. Altogether, the proposed regimen offers great potential to safely and specifically target hepatic tumors that are currently untreatable.


title: Locomotor performance and CNS responses to hypoxia in a cyclic nucleotide-gated channel mutant of adult Drosophila

abstract: Drosophila provides an excellent opportunity to explore the genetic basis for behavioral and CNS responses to hypoxia. Cyclic guanosine monophosphate (cGMP) modulates the speed of recovery from anoxia in adults and mediates hypoxia-related behaviors in larvae. Cyclic nucleotide-gated channels (CNG) and cGMP-activated protein kinase (PKG) are two cGMP downstream targets. PKG is involved in behavioral tolerance to hypoxia and anoxia in adults, however little is known about CNG channels. We used a CNGL mutant with reduced CNGL transcripts to investigate the contribution of CNGL to the hypoxia response. In control flies (w1118), hypoxia immediately reduced path length per minute in a locomotor assay. Flies took 30-40 mins in air to recover from 15 mins hypoxia. CNGL mutants had reduced locomotion under normoxia and impaired recovery from hypoxia, similar to the effects of pan-neural CNGL knockdown. In the CNGL mutants hypoxia caused an acute increase in path length per minute followed by a gradual increase during hypoxia. Basal levels of CNS extracellular K+ concentrations were reduced in the mutants. In response to hypoxia, the mutants had an increased extracellular K+ concentration change, reduced time to reach the K+ concentration peak, and delayed recovery time. Genetic manipulation to increase cGMP in the CNGL mutants eliminated the impairment of recovery from hypoxia and partially compensated for the effects of hypoxia on CNS K+. Although the neural mechanisms have yet to be determined, CNGL channels and cGMP signaling are involved in the hypoxia response of adult Drosophila.


title: Thresholded Partial Least Squares: Fast Construction of Interpretable Whole-brain Decoders

abstract: Recent neuroimaging research has shown that it is possible to decode mental states and predict future consumer behavior from brain activity data (a time-series of images). However, the unique characteristics (and high dimensionality) of neuroimaging data, coupled with a need for neuroscientifically interpretable models, has largely discouraged the use of the entire brains data as predictors. Instead, most neuroscientific research uses "regionalized" (partial-brain) data to reduce the computational burden and to improve interpretability (i.e., localizability of signal), at the cost of losing potential information. Here we propose a novel approach that can build whole-brain neural decoders (using the entire data set and capitalizing on the full correlational structure) that are both interpretable and computationally efficient. We exploit analytical properties of the partial least squares algorithm to build a regularized regression model with variable selection that boasts (in contrast to most statistical methods) a unique  fit-once-tune-later approach where users need to fit the model only once and can choose the best tuning parameters post-hoc. We demonstrate its efficacy in a large neuroimaging dataset against off-the-shelf prediction methods and show that our new method scales exceptionally with increasing data size, yields more interpretable results, and uses less computational memory, while retaining high predictive power.