title: Designing anti-predator training to maximize learning and efficacy assessments

abstract: Despite the growing need to use conservation breeding and translocations in species recovery, many attempts to reintroduce animals to the wild fail due to predation post-release. Released animals often lack appropriate behaviours for survival, including anti-predator responses. Anti-predator training--a method for encouraging animals to exhibit wariness and defensive responses to predators--has been used to help address this challenge with varying degrees of success. The efficacy of anti-predator training hinges on animals learning to recognize and respond to predators, but learning is rarely assessed, or interventions miss key experimental controls to document learning. An accurate measure of learning serves as a diagnostic tool for improving training if it otherwise fails to reduce predation. Here we present an experimental framework for designing anti-predator training that incorporates suitable controls to infer predator-specific learning and illustrate their use with the critically endangered Hawaiian crow,  alal[a] (Corvus hawaiiensis). We conducted anti-predator training within a conservation breeding facility to increase anti-predator behaviour towards a natural predator, the Hawaiian hawk,  io (Buteo solitaries). In addition to running live-predator training trials, we included two control groups, aimed at determining if responses could otherwise be due to accumulated stress and agitation, or to generalized increases in fear of movement. We found that without these control groups we may have wrongly concluded that predator-specific learning occurred. Additionally, despite generations in human care that can erode anti-predator responses,  alal[a] showed unexpectedly high levels of predatory wariness during baseline assessments. We discuss the implications of a learning-focused approach to training for managing endangered species that require improved behavioural competence for dealing with predatory threats, and the importance of understanding learning mechanisms in diagnosing behavioural problems.


title: Scaphoideus titanus Ball feeding behaviour on three grapevine cultivars with different susceptibilities to Flavescence doree

abstract: Scaphoideus titanus (Ball) is a grapevine-feeder leafhopper, and the most important vector of Flavescence doree of grapevine (FD), a disease associated with phytoplasmas belonging to ribosomal subgroups 16Sr-V-C and -D. FD is a major constraint to viticulture in several European countries and, so far, its control has relied on roguing of infected plants and insecticide applications against the vector. Detailed knowledge on different levels of the multifaceted phytoplasma-plant-vector relationship is required to envisage and explore more sustainable ways to control the disease spread. In the present work, S. titanus feeding behaviour was described on three grapevine cultivars: Barbera (susceptible to FD), Brachetto, and Moscato (tolerant to FD) using the Electrical Penetration Graph (EPG) technique. Interestingly, no differences were highlighted in the non-phloem feeding phases, thus suggesting that the tested cultivars have no major differences in the biochemical composition or structure of the leaf cuticle, epidermis or mesophyll, that can affect the first feeding activities. On the contrary, the results showed significant differences in leafhopper feeding behaviour in terms of the duration of the phloem feeding phase, longer on Barbera and shorter on Brachetto and Moscato, and of the frequency of interruption-salivation events inside the phloem, higher on Brachetto and Moscato. These findings indicate a possible preference for the Barbera cultivar, a better host for the leafhopper. Scaphoideus titanus feeding behaviour on Barbera correlates with an enhanced FDp transmission efficiency, thus explaining, at least in part, the higher susceptibility of this cultivar to FD. The mechanisms for the possible non-preference for Brachetto and Moscato are discussed, and an antixenosis is hypothesized. We propose that breeding for resistance against FD should take into account both plant traits associated with the response to the phytoplasmas and to the vector.

category: animal behavior and cognition


title: Molecular chaperone ability to inhibit amyloid-derived neurotoxicity, but not amorphous protein aggregation, depends on a conserved pH-sensitive Asp residue

abstract: Proteins can self-assemble into amyloid fibrils or amorphous aggregates and thereby cause disease. Molecular chaperones can prevent both these types of protein aggregation, but the respective mechanisms are not fully understood. The BRICHOS domain constitutes a disease-associated small heat shock protein-like chaperone family, with activities against both amyloid toxicity and amorphous protein aggregation. Here, we show that the activity of two BRICHOS domain families against Alzheimers disease associated amyloid-{beta} neurotoxicity to mouse hippocampi in vitro depends on a conserved aspartate residue, while the ability to suppress amorphous protein aggregation is unchanged by Asp to Asn mutations. The conserved Asp in its ionized state promotes structural flexibility of the BRICHOS domain and has a pKa value between pH 6.0-7.0, suggesting that chaperone effects against amyloid toxicity can be affected by physiological pH variations. Finally, the Asp is evolutionarily highly conserved in >3000 analysed BRICHOS domains but is replaced by Asn in some BRICHOS families and animal species, indicating independent evolution of molecular chaperone activities against amyloid fibril formation and non-fibrillar amorphous protein aggregation.

category: biochemistry


title: Rapid, Reference-Free Human Genotype Imputation with Denoising Autoencoders

abstract:
Genotype imputation is a foundational tool for population genetics. Standard statistical imputation approaches rely on the co-location of large whole-genome sequencing-based reference panels, powerful computing environments, and potentially sensitive genetic study data. This results in computational resource and privacy-risk barriers to access to cutting-edge imputation techniques. Moreover, the accuracy of current statistical approaches is known to degrade in regions of low and complex linkage disequilibrium.

Artificial neural network-based imputation approaches may overcome these limitations by encoding complex genotype relationships in easily portable inference models. Here we demonstrate an autoencoder-based approach for genotype imputation, using a large, commonly-used reference panel, and spanning the entirety of human chromosome 22. Our autoencoder-based genotype imputation strategy achieved superior imputation accuracy across the allele-frequency spectrum and across genomes of diverse ancestry, while delivering at least 4-fold faster inference run time relative to standard imputation tools.

category: bioinformatics


title: Physics Augmented U-Net: A High-Frequency Aware Generative Prior for Microscopy

abstract: A key challenge in optical microscopy is to image fast at high-resolution. To address this problem, we propose "Physics Augmented U-Net", which combines deep learning and structured illumination microscopy (SIM). In SIM, the structured illumination aliases out-of-band high-frequencies to the passband of the microscope; thus SIM captures some high-frequencies even when the image is sampled at low-resolution. To utilize these features, we propose a three-element method: 1) a modified U-Net model, 2) a physics-based forward model of SIM 3) an inference algorithm combining the two models. The modified U-Net architecture is similar to the seminal work, but the bottleneck is modified by concatenating two latent vectors, one encoding low-frequencies (LFLV), and the other encoding high-frequencies (HFLV). LFLV is learned by U-Net contracting path, and HFLV is learned by a second encoding path. In the inference mode, the high-frequency encoder is removed; HFLV is then optimized to fit the measured microscopy images to the output of the forward model for the generated image by the U-Net. We validated our method on two different datasets under different experimental conditions. Since a latent vector is optimized instead of a 2D image, the inference mode is less computationally complex. The proposed model is also more stable compared to other generative prior-based methods. Finally, as the forward model is independent of the U-Net, Physics Augmented U-Net can enhance resolution on any variation of SIM without further retraining.


title: StabilitySort: assessment of protein stability changes on a genome-wide scale to prioritise potentially pathogenic genetic variation

abstract:
SummaryMissense mutations that change protein stability are strongly associated with human inherited genetic disease. With the recent availability of predicted structures for all human proteins generated using the AlphaFold2 prediction model, genome-wide assessment of the stability effects of genetic variation can, for the first time, be easily performed. This facilitates the interrogation of personal genetic variation for potentially pathogenic effects through the application of stability metrics. Here, we present a novel algorithm to prioritise variants predicted to strongly destabilise essential proteins, available as both a standalone software package and a web-based tool. We demonstrate the utility of this tool by showing that at values of the Stability Sort Z-score above 1.6, pathogenic, protein-destabilising variants from ClinVar are detected at a 58% enrichment, over and above the destabilising (but presumably non-pathogenic) variation already present in the HapMap NA12878 genome.

Availability and ImplementationStabilitySort is available as both a web service (http://130.56.244.113/StabilitySort/) and can be deployed as a standalone system (https://gitlab.com/baaron/StabilitySort).

ContactDan.Andrews@anu.edu.au


title: A single-cell transcriptomic landscape of innate and adaptive intratumoral immunity in triple negative breast cancer during chemo- and immunotherapies

abstract:
Breast cancer (BC) constitutes a major health problem worldwide, making it the most common malignancy in women. Current treatment options for BC depend primarily on histological type, molecular markers, clinical aggressiveness and stage of disease. Immunotherapy, such as anti-PD-1, have shown combinatorial clinical activity with chemotherapy in triple negative breast cancer (TNBC) delineating some therapeutic combinations as more effective than others. However, a clear overview of the main immune cell populations involved in these treatments has never been provided.

Here, an assessment of the immune landscape in the tumour microenvironment (TME) of two TNBC mouse models (4T1 and EMT6 cell lines) has been performed using single-cell RNA sequencing (scRNA-seq) technology. Specifically, immune cells were evaluated in untreated conditions and after being treated with chemotherapy or immunotherapy used as single agents or in combination. A decrease of regulatory T cells, compared to the untreated TME, was found in treatments with in vivo efficacy as well as {gamma}{delta} T cells, which have a pro-tumoral activity in mice. Focusing on Cd8 T cells, across all the conditions, a general increase of exhausted-like Cd8 T cells was confirmed in pre-clinical treatments with low efficacy; on the other hand, an opposite trend was found for the proliferative Cd8 T cells. Regarding macrophages, M2-like cells were found enriched in treatments with low efficacy while opposite behaviour was associated with M1-like macrophages. For both cell lines, similar proportions of B cells were detected with an increase of proliferative B cells in treatments that involved cisplatin in combination with anti-PD-1. The fine-scale characterization of the immune TME in this work can lead to new insights on the diagnosis and treatment of TNBC for a possible application at the clinical level.

category: cancer biology


title: GBS and a newly developed mRNA-GBS approach to link population genetic and transcriptome analyses reveal pattern differences between sites and treatments in red clover (Trifolium pratense L.)

abstract: The important worldwide forage crop red clover (Trifolium pratense L.) is widely cultivated as cattle feed and for soil improvement. Wild populations and landraces have great natural diversity that could be used to improve cultivated red clover. However, to date, there is still insufficient knowledge about the natural genetic and phenotypic diversity of the species. Here, we developed a low-cost transcriptome analysis (mRNA-GBS) with reduced complexity and compared the results with population genetic (GBS) and previously published mRNA-Seq data, to assess whether analysis of intraspecific variation within and between populations and transcriptome responses is possible simultaneously. The mRNA-GBS approach was successful. SNP analyses from the mRNA-GBS approach revealed comparable patterns to the GBS results, but it was not possible to link transcriptome analyses with reduced complexity and sequencing depth to previously published greenhouse and field expression studies. The use of short sequences upstream of the poly(A) tail of mRNA to reduce complexity are promising approaches that combine population genetics and expression profiling to analyze many individuals with trait differences simultaneously and cost-effectively, even in non-model species. Our mRNA-GBS approach revealed too many additional short mRNA sequences, hampering sequence alignment depth and SNP recovery. Optimizations are being discussed. Nevertheless, our study design across different regions in Germany was also challenging as the use of differential expression analyses with reduced complexity, in which mRNA is fragmented at specific sites rather than randomly, is most likely counteracted under natural conditions by highly complex plant reactions at low sequencing depth.

category: plant biology


title: A new protocol for speed vernalisation of winter cereals

abstract: There are many challenges facing the development of high-yielding, nutritious crops for future environments. One limiting factor is generation time, which prolongs research and plant breeding timelines. Recent advances in speed breeding protocols have dramatically reduced generation time for many short-day and long-day species by optimising light and temperature conditions during plant growth. However, winter crops with a vernalisation requirement still require up to 6-10 weeks in low-temperature conditions before transition to reproductive development. Here, we tested a suite of environmental conditions and protocols to investigate if vernalisation can be satisfied more efficiently. We identified a vernalisation method consisting of exposing seeds at the soil surface to an extended photoperiod of 22 h day:2 h night at 10{degrees}C with transfer to speed breeding conditions that dramatically reduces generation time in both winter wheat (Triticum aestivum) and winter barley (Hordeum vulgare). Implementation of this protocol achieved up to five generations per year for winter wheat or barley, instead of the two typically obtained under standard vernalisation and plant growth conditions. The protocol has great potential to enhance training and to accelerate research, pre-breeding, and breeding outcomes focussed on winter crop improvement.


title: Mutated neuronal voltage-gated CaV2.1 channels causing familial hemiplegic migraine 1 increase the susceptibility for cortical spreading depolarization and seizures and worsen outcome after experimental traumatic brain injury

abstract: Patients suffering from familial hemiplegic migraine type 1 (FHM1) may have a disproportionally severe outcome after head trauma, but the underlying mechanisms are unclear. Hence, we subjected knock-in mice carrying the severer S218L or milder R192Q FHM1 gain-of-function missense mutation in the CACNA1A gene that encodes the 1A subunit of neuronal voltage-gated CaV2.1 (P/Q-type) calcium channels and their wild-type (WT) littermates to experimental traumatic brain injury (TBI) by controlled cortical impact (CCI) and investigated cortical spreading depolarizations (CSDs), lesion volume, brain edema formation, and functional outcome. After TBI, all mutant mice displayed considerably more CSDs and seizures than WT mice, while S218L mutant mice had a substantially higher mortality. Brain edema formation and the resulting increase in intracranial pressure was more pronounced in mutant mice, while only S218L mutant mice had larger lesion volumes and worse functional outcome. Here we show that gain of CaV2.1 channel function worsens histopathological and functional outcome after TBI in mice. This phenotype was associated with a higher number of CSDs, increased seizure activity, and more pronounced brain edema formation. Hence, our results suggest increased susceptibility for CSDs and seizures as potential mechanisms for bad outcome after TBI in FHM1 mutation carriers.

category: neuroscience


title: Rapid volumetric brain changes after acute psychosocial stress

abstract:
Rapid structural brain plasticity after acute stress has been shown in animals. It is unknown whether such stress-related brain changes also occur in humans, in which they have been found, using structural magnetic resonance imaging (MRI), after motor learning and visual stimulation. We here investigated grey matter volume (GMV) changes after acute stress in humans and tested their relation to psychophysiological stress measures.

Sixty-seven healthy men (25.8{+/-}2.7 years) completed a standardized psychosocial laboratory stressor (Trier Social Stress Test) or a control version while blood, saliva, heart rate, and psychometrics were sampled. T1-weighted MP2RAGE images at 3T MRI were acquired 45 min before and 90 min after intervention onset. GMV changes were analysed using voxel-based morphometry. Associations with endocrine, autonomic, and subjective stress measures were tested with linear models.

We found significant group-by-time interactions in several brain clusters including anterior/mid-cingulate cortices and bilateral insula: GMV was increased in the stress group relative to the control group, in which several clusters showed a GMV decrease. We found no significant group-by-time interaction for other MRI parameters, including cerebral blood flow, but a significant association of GMV changes with state anxiety and heart rate variability changes.

In summary, we show rapid GMV changes following acute psychosocial stress in humans. The results suggest that endogenous circadian brain changes are counteracted by acute stress and generally emphasize the influence of stress on the brain.


title: Correspondence between gene expression and neurotransmitter receptor and transporter density in the human brain

abstract: Neurotransmitter receptors modulate the signaling between neurons. Thus, neurotransmitter receptors and transporters play a key role in shaping brain function. Due to the lack of comprehensive neurotransmitter receptor/transporter density datasets, microarray gene expression is often used as a proxy for receptor densities. In the present report, we comprehensively test the expression-density association for a total of 27 neurotransmitter receptors, receptor binding-sites, and transporters across 9 different neurotransmitter systems, using both PET and autoradiography imaging modalities. We find poor spatial correspondences between gene expression and density for all neurotransmitter receptors and transporters except four single-protein metabotropic receptors (5-HT1A, D2, CB1, and MOR). These expression-density associations are related to population variance and change across different classes of laminar differentiation. Altogether, we recommend using direct measures of receptor and transporter density when relating neurotransmitter systems to brain structure and function.


title: Emotion Dynamics as Hierarchical Bayesian Inference in Time

abstract: What fundamental property of our environment would be most valuable and optimal in characterizing the emotional dynamics we experience in our daily life? Empirical work has shown that an accurate estimation of uncertainty is necessary for our optimal perception, learning, and decision-making. However, the role of this uncertainty in governing our affective dynamics remains unexplored. Using Bayesian encoding, decoding and computational modelling, we show that emotional experiences naturally arise due to ongoing uncertainty estimations in a hierarchical neural architecture. This hierarchical organization involves a number of prefrontal sub-regions, with the lateral orbitofrontal cortex having the highest representational complexity of uncertainty. Crucially, this representational complexity, was sensitive to temporal fluctuations in uncertainty and was predictive of participants predisposition to anxiety. Furthermore, the temporal dynamics of uncertainty revealed a distinct functional double dissociation within the OFC. Specifically, the medial OFC showed higher connectivity with the DMN, while the lateral OFC with that of the FPN in response to the evolving affect. Finally, we uncovered a temporally predictive code updating individuals beliefs swiftly in the face of fluctuating uncertainty in the lateral OFC. A biologically relevant and computationally crucial parameter in theories of brain function, we extend uncertainty to be a defining component of complex emotions.


title: NAc-DBS corrects depression-like behaviors in CUMS mouse model via disinhibition of DA neurons in the VTA

abstract: Major depressive disorder (MDD) is characterized by diverse debilitating symptoms that include loss of motivation and anhedonia. If multiple medications, psychotherapy, and electroconvulsive therapy fail in some patients with MDD, their condition is then termed treatment - resistant depression (TRD). MDD can be associated with abnormalities in the reward-system-dopaminergic mesolimbic pathway, in which the nucleus accumbens (NAc) and ventral tegmental area (VTA) play major roles. Deep brain stimulation (DBS) applied to the NAc alleviates the depressive symptoms of MDD. However, the mechanism underlying the effects of this DBS has remained elusive. In this study, using the chronic unpredictable mild stress (CUMS) mouse model, we investigated the behavioral and neurobiological effects of NAc - DBS on the multidimensional depression - like phenotypes induced by CUMS by integrating behavioral, in vivo microdialysis coupled with high-performance liquid chromatography - electrochemical detector (HPLC - ECD), calcium imaging, pharmacological, and genetic manipulation methods in freely moving mice. We found that long-term and repeated, but not single, NAc-DBS induced robust antidepressant responses in CUMS mice. Moreover, even a single trial NAc-DBS led to the elevation of the {gamma}-aminobutyric acid (GABA) neurotransmitter, accompanied by the increase in dopamine (DA) neuron activity in the VTA. Both the inhibition of the GABAA receptor activity and knockdown of the GABAA-1 gene in VTA-GABA neurons blocked the antidepressant effect of NAc-DBS in CUMS mice. Our results showed that NAc- DBS could disinhibit VTA-DA neurons by regulating the level of GABA and the activity of VTA-GABA in the VTA and could finally correct the depression-like behaviors in the CUMS mouse model.


title: ANXA2 enhances virus replication through negatively regulating cGAS-STING and RLRs2 mediated signal pathway

abstract:
Host nucleic acid receptors can recognize the viral DNA or RNA upon virus infection, which further triggers multiple signal pathways to promote the translocation of the interferon regulatory factor 3 (IRF3) into nucleus and produce type I interferon (IFN), leading to the host antiviral response. Here, we report a novel negative regulator Annexin A2 (ANXA2) that regulates type I IFN production through multiple mechanisms. Ectopic expression of ANXA2 inhibited the production of type I IFN induced by DNA- and RNA viruses and enhanced virus replication, while knockout of ANXA2 expression enhanced the production of type I IFN and inhibited virus replication. Mechanistically, ANXA2 not only disrupted MDA5 recruiting MAVS, but also inhibited the interaction between MAVS and TRAF3 upon RNA virus infection. In addition, ANXA2 impacted the translocation of STING from endoplasmic reticulum to Golgi apparatus upon DNA virus infection. Interestingly, ANXA2 also inhibited IRF3 phosphorylation and nuclear translocation through competing with TANK-binds kinase 1 (TBK1) and inhibitor-{kappa}B kinase {varepsilon} (IKK{varepsilon}) for binding to IRF3. Anxa2 deficiency in vivo increased the production of type I IFN, which resulted in suppression of encephalomyocarditis virus (EMCV) replication. Our findings reveal that ANXA2, as a negative regulator of type I IFN production, plays an important role in regulating the host antiviral responses.

Author summaryAnnexin is a family of evolutionarily conserved multi-gene proteins, which are widely distributed in various tissues and cells of plants and animals. These proteins can reversibly bind to phospholipid membranes and to calcium ions (Ca2+). To date, several studies have confirmed that some members of the Annexin family regulate the antiviral innate immune response. Until now, regulation of the production of type I IFN by ANXA2 is not reported. In this study, ANXA2 were found to strongly inhibit the production of type I IFN, leading to increased virus replication while knockout of ANXA2 expression inhibited virus replication by increasing the amount of IFN. Compared with wild-type littermates, ANXA2 deficiency mice produced more type I IFN to inhibit virus replication. Our results provide methanistic insights into the novel role of ANXA2 in the antiviral innate immune responses.

category: immunology


title: FIRST REPORT OF HALOBACTERIA DOMINANCE IN A TROPICAL CAVE MICROBIOME

abstract: Scarce studies on microbial diversity in tropical caves have been published, a subterranean system still neglected from a microbiological point of view. Although most published studies are about temperate caves, usually archaeas and fungi have less attention than bacterial communities. Here, the microbiome structure and composition in a tropical cave system, as well the main environmental drivers, were studied during the wet and dry season. Soil and sediments from three different habitats at the cave (surface, entrance cave and dark zone) were sampled. Samples were characterized (temperature, air and substrate humidity, salinity, pH, nitrogen and organic carbon content, and chemical composition) and the microbiome was assessed by high-throughput sequencing, using amplicon sequencing (16S and ITS). Prokaryotic communities were dominated by Halobacteria, Actinobacteria and Bacilli, while fungal communities showed high abundance of Sordariomycetes. Microbiomes from the cave entrance, where a significantly elevated salinity levels were found, supported up to 63% of Haloarchaea compared to the other habitats studied. Differences in community structure were significant between habitats, but no influence of the season was observed. Main environmental drivers of community assembly included nitrogen and organic carbon content, temperature, and salinity. This is the first report of Halobacteria dominance in cave habitats, where they likely play important roles in nitrogen and phosphorus cycles. The cave entrance had lower diversity, but higher degree of microbial endemism, which characterize it as an important cave ecotone. The prevalence of heterotrophic microbial groups implies trophic structure based on detritivores, particularly in the dark zones. Our study brings new insights on microbiome composition in the underexplored tropical cave habitats.

category: microbiology


title: Reproductive developmental dynamics during the early and intermediate prepuberal periods in Nelore breed (Bos indicus) calves

abstract: The demand for calves as oocyte donors for in vitro embryo production is growing. However, Bos indicus have a late puberty, and some aspects of the reproductive physiology during the prepubertal period remain unclear. We characterized endocrine and morpho-functional reproductive features in Nelore calves (n=8) at 2-5 (early prepubertal period, EPP) and from 8-11 months old (mo., intermediate prepubertal period, IPP). The calves ovaries and uterus were B-mode transrectal ultrasonography examined, and blood samples were collected every second week. The antral follicles number and size, and ovarian and uterine horn diameters, were recorded, and plasma FSH and LH concentrations were measured (RIA). Non-pregnant, non-lactating cyclic Nelore cows (n=8) were used as controls for endocrine endpoints. Somatic development was monitored by monthly weighing, and 3D scanning of the rump area. The somatic and endocrine endpoints were compared within and between EPP and IPP, and between each period and control cows. Associations were determined by the Spearman correlation method, and the developmental rates were determined by non-linear regression. All morphological endpoints, except antral follicle count, increased (P < 0.001) from the EPP to the IPP. However, within each period differences occurred only at EPP. During the EPP LH and FSH plasma concentrations were similar (P > 0.05), whereas during the IPP LH was lower (P < 0.05) and FSH was higher (P < 0.001) than control cows. The EPP calves showed moderate to high positive correlations among ovarian, uterine, and somatic endpoints. Conversely, the IPP such correlations were mostly weak. In summary, distinct ovarian activity and development patterns of primary and secondary sexual characteristics occurred in Nelore calves at 2-5 mo compared to 8-11 mo.

category: developmental biology


title: On eco-evolutionary dynamics of phenologies in competitive communities and their robustness to climate change

abstract: Global changes currently cause temporal shifts in the favourable conditions for different phases of species life cycles. Phenologies characterizing temporal presence, may adapt through heritable evolution in response to these changes. Given a community context, this evolution may cause a change in the phenology overlap and thus a change of interspecific interactions such as competition. Using a model in which phenologies compete and coevolve, we study the conditions under which diversity emerges, as well as their annual distribution. We find that the environment richness (food quantity, light, pollinators...) and competition constrain the diversity and spread of phenologies. A robust pattern of phenologies distribution emerges consistent with Swedish flowering observations. Once a stable community is reached, we apply a progressive change in environmental conditions. We found that adaptation eventually restored diversity, but that the simulated change often led to numerous extinctions due to increased competition. The percentage of diversity lost depends on the speed of change and on the initial diversity. Phenologies already pre-adapted to the new environmental conditions drive the restoration of diversity after the change. We finally study a spatial version of the model in which local communities are organized along an environmental gradient. Pre-change, allowing dispersal decreases the local adaptation of phenologies to their local fixed environmental conditions. Dispersal however largely enhances the maintenance of biodiversity in changing environments, though its benefits are not homogeneous in space. Evolution remains the only rescue mechanism for southern phenotypes.

category: evolutionary biology


title: Extreme sampling for genetic rare variant association analysis of dichotomous traits with focus on infectious disease susceptibility

abstract: As genomic sequencing becomes more accurate and less costly, large cohorts and consortiums of cohorts are providing high power for rare variant association studies for many conditions. When large sample sizes are not attainable and the phenotype under study is continuous, an extreme phenotypes design can provide high statistical power with a small to moderate sample size. We extend the extreme phenotypes design to the dichotomous infectious disease outcome by sampling on extremes of the pathogenic exposure instead of sampling on extremes of phenotype. We use a likelihood ratio test (LRT) to test the significance of association between infection status and presence of susceptibility rare variants. More than 10 billion simulations are studied to assess the method. The method results in high sample enrichment for rare variants affecting susceptibility. Greater than 90% power to detect rare variant associations is attained in reasonable scenarios. The ordinary case-control design requires orders of magnitude more samples to achieve the same power. The Type I error rate of the LRT is accurate even for p-values < 10-7. We find that erroroneous exposure assessment can lead to power loss more severe than excluding the observations with errors. Nevertheless, careful sampling on exposure extremes can make a study feasible by providing adequate statistical power. Limitations of this method are not unique to this design, and the power is never less than that of the ordinary case-control design. The method applies without modification to other dichotomous outcomes that have strong association with a continuous covariate.

category: genomics


title: T cell receptor repertoire signatures associated with COVID-19 severity

abstract: T cell receptor (TCR) repertoires are critical for antiviral immunity. Determining the TCR repertoires composition, diversity, and dynamics and how they change during viral infection can inform the molecular specificity of viral infection such as SARS-CoV-2. To determine signatures associated with COVID-19 disease severity, here we performed a large-scale analysis of over 4.7 billion sequences across 2,130 TCR repertoires from COVID-19 patients and healthy donors. TCR repertoire analyses from these data identified and characterized convergent COVID-19 associated CDR3 gene usages, specificity groups, and sequence patterns. T cell clonal expansion was found to be associated with upregulation of T cell effector function, TCR signaling, NF-kB signaling, and Interferon-gamma signaling pathways. Machine learning approaches accurately predicted disease severity for patients based on TCR sequence features, with certain high-power models reaching near-perfect AUROC scores across various predictor permutations. These analyses provided an integrative, systems immunology view of T cell adaptive immune responses to COVID-19.


title: Pupillometry and the vigilance decrement: Task-evoked but not baseline pupil measures reflect declining performance in visual vigilance tasks

abstract: Baseline and task-evoked pupil measures are known to reflect the activity of the nervous systems central arousal mechanisms. With the increasing availability, affordability and flexibility of video-based eye tracking hardware, these measures may one day find practical application in real-time biobehavioral monitoring systems to assess performance or fitness for duty in tasks requiring vigilant attention. But real-world vigilance tasks are predominantly visual in their nature and most research in this area has taken place in the auditory domain. Here we explore the relationship between pupil size--both baseline and task-evoked--and behavioral performance measures in two novel vigilance tasks requiring visual target detection: 1) a traditional vigilance task involving prolonged, continuous, and uninterrupted performance (n = 28), and 2) a psychomotor vigilance task (n = 25). In both tasks, behavioral performance and task- evoked pupil responses declined as time spent on task increased, corroborating previous reports in the literature of a vigilance decrement with a corresponding reduction in task-evoked pupil measures. Also in line with previous findings, baseline pupil size did not show a consistent relationship with performance measures. We discuss our findings considering the adaptive gain theory of locus coeruleus function and question the validity of the assumption that baseline (prestimulus) pupil size and task-evoked (poststimulus) pupil measures correspond to the tonic and phasic firing modes of the LC.


title: Variance in Variants: Propagating Genome Sequence Uncertainty into Phylogenetic Lineage Assignment

abstract:
Genetic sequencing is subject to many different types of errors, but most analyses treat the resultant sequences as if they are known without error. Next generation sequencing methods rely on significantly larger numbers of reads than previous sequencing methods in exchange for a loss of accuracy in each individual read. Still, the coverage of such machines is imperfect and leaves uncertainty in many of the base calls. On top of this machine-level uncertainty, there is uncertainty induced by human error, such as errors in data entry or incorrect parameter settings. In this work, we demonstrate that the uncertainty in sequencing techniques will affect downstream analysis and propose a straightforward method to propagate the uncertainty.

Our method uses a probabilistic matrix representation of individual sequences which incorporates base quality scores as a measure of uncertainty that naturally lead to resampling and replication as a framework for uncertainty propagation. With the matrix representation, resampling possible base calls according to quality scores provides a bootstrap- or prior distribution-like first step towards genetic analysis. Analyses based on these re-sampled sequences will include a more complete evaluation of the error involved in such analyses.

We demonstrate our resampling method on SARS-CoV-2 data. The resampling procedures adds a linear computational cost to the analyses, but the large impact on the variance in downstream estimates makes it clear that ignoring this uncertainty may lead to overly confident conclusions. We show that SARS-CoV-2 lineage designations via Pangolin are much less certain than the bootstrap support reported by Pangolin would imply and the clock rate estimates for SARS-CoV-2 are much more variable than reported.


title: Mapping lesion, structural disconnection, and functional disconnection to symptoms in semantic aphasia

abstract: Patients with semantic aphasia have impaired control of semantic retrieval, often accompanied by executive dysfunction following left hemisphere stroke. Many but not all of these patients have damage to the left inferior frontal gyrus, important for semantic and cognitive control. Yet semantic and cognitive control networks are highly distributed, including posterior as well as anterior components. Accordingly, semantic aphasia might not only reflect local damage but also white matter structural and functional disconnection. Here we characterise the lesions and predicted patterns of structural and functional disconnection in individuals with semantic aphasia and relate these effects to semantic and executive impairment. Impaired semantic cognition was associated with infarction in distributed left- hemisphere regions, including in the left anterior inferior frontal and posterior temporal cortex. Lesions were associated with executive dysfunction within a set of adjacent but distinct left frontoparietal clusters. Performance on executive tasks was also associated with interhemispheric structural disconnection across the corpus callosum. In contrast, poor semantic cognition was associated with small left-lateralized structurally disconnected clusters, including in the left posterior temporal cortex. Little insight was gained from functional disconnection symptom mapping. These results demonstrate that while left- lateralized semantic and executive control regions are often damaged together in stroke aphasia, these deficits are associated with distinct patterns of structural disconnection, consistent with the bilateral nature of executive control and the left-lateralized yet distributed semantic control network.


title: Targeted memory reactivation in human REM sleep elicits recurrent, temporally compressed reactivation

abstract: Memories are reactivated during non-rapid eye movement (NREM) sleep, but the question of whether equivalent reactivation also occurs in rapid eye movement (REM) sleep is hotly debated. To examine this, we used a technique called targeted memory reactivation (TMR) in which sounds are paired with learned material in wake, and then re-presented in subsequent sleep to trigger reactivation. We then used machine learning classifiers to identify TMR-induced reactivation in REM. The reactivation we measured was temporally compressed by approximately five times during REM compared to wakeful performance of the task, and often occurred twice after a single TMR cue. Reactivation strength positively predicted overnight performance improvement and was only apparent in trials with high theta activity. These findings provide strong evidence for memory reactivation in human REM sleep after TMR as well as an initial characterisation of this reactivation.


title: Benchmarking Geometric Deep Learning for Cortical Segmentation and Neurodevelopmental Phenotype Prediction

abstract: The emerging field of geometric deep learning extends the application of convolutional neural networks to irregular domains such as graphs, meshes and surfaces. Several recent studies have explored the potential for using these techniques to analyse and segment the cortical surface. However, there has been no comprehensive comparison of these approaches to one another, nor to existing Euclidean methods, to date. This paper benchmarks a collection of geometric and traditional deep learning models on phenotype prediction and segmentation of sphericalised neonatal cortical surface data, from the publicly available Developing Human Connectome Project (dHCP). Tasks include prediction of postmenstrual age at scan, gestational age at birth and segmentation of the cortical surface into anatomical regions defined by the M-CRIB-S atlas. Performance was assessed not only in terms of model precision, but also in terms of network dependence on image registration, and model interpretation via occlusion. Networks were trained both on sphericalised and anatomical cortical meshes. Findings suggest that the utility of geometric deep learning over traditional deep learning is highly task-specific, which has implications for the design of future deep learning models on the cortical surface. The code, and instructions for data access, are available from https://github.com/Abdulah-Fawaz/Benchmarking-Surface-DL.


title: Clinical implementation of single-cell RNA sequencing using liver fine needle aspirate tissue sampling and centralized processing captures compartment specific immuno-diversity

abstract: Blood samples are frequently collected in human studies of the immune system but poorly represent tissue-resident immunity. Understanding the immunopathogenesis of tissue-restricted diseases, such as chronic hepatitis B, necessitates direct investigation of local immune responses. We developed a workflow that enables frequent, minimally invasive collection of liver fine-needle aspirates in multi-site international studies and centralized single-cell RNA sequencing data generation using the Seq-Well S3 picowell-based technology. All immunological cell types were captured, including liver macrophages, and showed distinct compartmentalization and transcriptional profiles, providing a systematic assessment of the capabilities and limitations of peripheral blood samples when investigating tissue-restricted diseases. The ability to electively sample the liver of chronic viral hepatitis patients and generate high-resolution data will enable multi-site clinical studies to power fundamental and therapeutic discovery.


title: Draft genome of the Brazilian railroad worm Phrixothrix hirtus E.Olivier (Phengodidae: Coleoptera)

abstract:
The Neotropical region is the richest in bioluminescent Coleoptera species, however, its bioluminescence megadiversity is still underexplored in terms of genomic organization and evolution, mainly within the Phengodidae family. The railroad worm Phrixothrix hirtus is an important biological model and symbolic species due to its bicolor bioluminescence, being the only organism that produces true red light among bioluminescent terrestrial species. Here, we performed the partial genome assembly of P. hirtus, combining short and long reads generated with Illumina sequencing, providing an important source of genomic information and a framework for comparative genomic analyses for the evaluation of the bioluminescent system in Elateroidea. The estimated genome size has [~]3.4Gb, 32% of GC content, and 67% of repetitive elements, being the largest genome described in the Elateroidea superfamily. Several events of gene family expansions associated with anatomical development and morphogenesis, as well as distinct odorant-binding receptors and retrotransposable elements were found in this genome. Similar molecular functions and biological processes are shared with other studied species of Elateriformia. Common genes putatively associated with bioluminescence production and control, including two luciferase genes that displayed 7 exons and 6 introns, and genes that could be involved in luciferin biosynthesis were found, indicating that there are no clear differences about the presence or absence of gene families associated with bioluminescence in Elateroidea. In P. hirtus the conversion of L- to D-luciferin seems to involve additional steps using a Palmitoyl-CoA thioesterase instead of an Acyl-CoA synthetase, which was found in Lampyridae species.

HighlightsO_LIFirst draft genome assembly of Phengodidae, the largest one described in Coleoptera;
C_LIO_LIGene family expansions associated with anatomical development and morphogenesis;
C_LIO_LIBioluminescent control and luciferin biosynthesis genes are common within Elateroidea;
C_LIO_LIDespite similar bioluminescent system, metabolic routes may have evolved independently;
C_LI


title: Structural basis of dynamic P5CS filaments

abstract: The bifunctional enzyme {Delta}1-pyrroline-5-carboxylate synthase (P5CS) is central to the synthesis of proline and ornithine. Pathogenic mutations in P5CS gene (ALDH18A1) lead to neurocutaneous syndrome and skin relaxation connective tissue disease in humans, and P5CS deficiency seriously damages the ability to resist adversity in plants, which has an essential role in agriculture and human health. Recently, P5CS has been demonstrated forming the cytoophidium in vivo and filaments in vitro. However, the underlying mechanism for the function of P5CS filamentation and catalyze the synthesis of P5C is hardly accessible without structural basis. Here, we have succeeded in determining the full-length structures of Drosophila P5CS filament in three states at resolution from 3.1 to 4.3 [A] under cryo-electron microscopy, we observed the distinct ligand-binding states and conformational changes for GK and GPR domain separately. These structures show the distinctive spiral filament is assembled by P5CS tetramers and stabilized by multiple interfaces. Point mutations that deplete such interactions disturb P5CS filamentation and greatly reduce the activity. Our findings reveal a previously undescribed mechanism that filamentation is crucial for the coordination between GK and GPR domains, and provide insights into structural basis for catalysis function of P5CS filament.


title: Initiation of ensemble kinesin-3 motility is regulated by the rigidity of cargo-motor attachment

abstract: Intracellular cargo transport is powered by molecular motors that move on their respective filamentous tracks. A key component in this process is the tether between cargo and motor, which is often connected by long slender coiled-coils. Several studies have identified mechanisms that regulate cargo transport and can be broadly categorized into regulation of the motor ATPase activity by autoinhibition, cargo adapters and modifications in the cytoskeletal tracks. The regulatory effects of cargo-motor linkers have been described in kinesin-3 subfamily motors. However, the effects of cargo-motor linker rigidity on ensemble cargo transport has not been explored. Here we have built a DNA origami scaffold, which can be tethered with multiple kinesin-3 motors using either single or double-stranded DNA linkages, mimicking rigid versus flexible cargo-motor linkages. Using this system, we show that regardless of the motor numbers attached to the cargo, only linkers with a lesser degree of freedom allow motors to engage with microtubule tracks. Together, our work identifies that the rigidity of cargo-motor linkages influences motor motility. This opens up the possibilities to identify new factors that can influence the rigidity of cargo-motor linkages that in turn can regulate intracellular cargo transport.


title: A versatile new ubiquitin detection and purification tool derived from a bacterial deubiquitylase

abstract: Protein ubiquitylation is an important post-translational modification affecting an wide range of cellular processes. Due to the low abundance of ubiquitylated species in biological samples, considerable effort has been spent on developing methods to purify and detect ubiquitylated proteins. We have developed and characterized a novel tool for ubiquitin detection and purification based on OtUBD, a high-affinity ubiquitin-binding domain derived from an Orientia tsutsugamushi deubiquitylase. We demonstrate that OtUBD can be used to purify both monoubiquitylated and polyubiquitylated substrates from yeast and human tissue culture samples and compare their performance with existing methods. Importantly, we found conditions for either selective purification of covalently ubiquitylated proteins or co-isolation of both ubiquitylated proteins and their interacting proteins. As a proof-of-principle for these newly developed methods, we profiled the ubiquitylome and ubiquitin-associated proteome of the yeast Saccharomyces cerevisiae. Combining OtUBD affinity purification with quantitative proteomics, we identified potential substrates for E3 ligases Bre1 and Pib1. OtUBD provides a versatile, efficient, and economical tool for ubiquitin researchers with specific advantages over other methods, such as in detecting monoubiquitylation or ubiquitin linkages to noncanonical sites.


title: miR2Trait: an integrated resource for investigating miRNA-disease associations

abstract: MicroRNAs are key components of cellular regulatory networks, and breakdown in miRNA function could lead to cascading effects culminating in pathophenotypes. A better understanding of the role of miRNAs in diseases would aid human health. Here, we have devised a method for comprehensively mapping the associations between miRNAs and diseases by merging on a common key between two curated omics databases. The resulting bidirectional resource, miR2Trait, is more detailed than earlier catalogues, uncovers new relationships, and includes analytical utilities to interrogate and extract knowledge from these datasets. The resource could aid in identifying the disease enrichment of a user-given set of miRNAs and analyzing the miRNA profile of a specified diseasome. miR2Trait is available as both a web-server (https://sas.sastra.edu/pymir18) and an open-source command-line interface (https://github.com/miR2Trait) under MIT license for both commercial and non-commercial use. The datasets are available for download at: https://doi.org/10.6084/m9.figshare.8288825.


title: Polycomb Group Ring Finger Protein 6 suppresses Myc-induced lymphomagenesis

abstract: Max is an obligate dimerization partner for the Myc transcription factors and for several repressors, such as Mnt, Mxd1-4 and Mga, collectively thought to antagonize Myc function in transcription and oncogenesis. Mga, in particular, is part of the variant Polycomb group repressive complex PRC1.6. Here, we show that ablation of the distinct PRC1.6 subunit Pcgf6 - but not Mga - accelerates Myc-induced lymphomagenesis in E{micro}-myc transgenic mice. Unexpectedly, however, Pcgf6 loss shows no significant impact on transcriptional profiles, in neither pre-tumoral B-cells, nor lymphomas. Altogether, these data unravel an unforeseen, Mga- and PRC1.6-independent tumor suppressor activity of Pcgf6.


title: Quasi-Entropy Closure: A Fast and Reliable Approach to Close the Moment Equations of the Chemical Master Equation

abstract:
MotivationThe Chemical Master Equation is the most comprehensive stochastic approach to describe the evolution of a (bio-)chemical reaction system. Its solution is a time-dependent probability distribution on all possible configurations of the system. As the number of possible configurations is typically very large, the Master Equation is often practically unsolvable. The Method of Moments reduces the system to the evolution of a few moments of this distribution, which are described by a system of ordinary differential equations. Those equations are not closed, since lower order moments generally depend on higher order moments. Various closure schemes have been suggested to solve this problem, with different advantages and limitations. Two major problems with these approaches are first that they are open loop systems, which can diverge from the true solution, and second, some of them are computationally expensive.

ResultsHere we introduce Quasi-Entropy Closure, a moment closure scheme for the Method of Moments which estimates higher order moments by reconstructing the distribution that minimizes the distance to a uniform distribution subject to lower order moment constraints. Quasi-Entropy closure is similar to Zero-Information closure, which maximizes the information entropy. Results show that both approaches outperform truncation schemes. Moreover, Quasi-Entropy Closure is computationally much faster than Zero-Information Closure. Finally, our scheme includes a plausibility check for the existence of a distribution satisfying a given set of moments on the feasible set of configurations. Results are evaluated on different benchmark problems.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/470753v1_ufig1.gif" ALT="Figure 1">
View larger version (31K):
org.highwire.dtl.DTLVardef@8c2beforg.highwire.dtl.DTLVardef@9876f7org.highwire.dtl.DTLVardef@1aa5db7org.highwire.dtl.DTLVardef@1eed4b4_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: systems biology


title: Crowding changes appearance systematically in peripheral, amblyopic, and developing vision

abstract:
Visual crowding is the disruptive effect of clutter on object recognition. Although most prominent in adult peripheral vision, crowding also disrupts foveal vision in typically-developing children and those with strabismic amblyopia. Do these crowding effects share the same mechanism? Here we exploit observations that crowded errors in peripheral vision are not random: target objects appear either averaged with the flankers (assimilation), or replaced by them (substitution). If amblyopic and developmental crowding share the same mechanism then their errors should be similarly systematic. We tested foveal vision in children aged 3-8 years with typical vision or strabismic amblyopia, and peripheral vision in typical adults. The perceptual effects of crowding were measured by requiring observers to adjust a reference stimulus to match the perceived orientation of a target  Vac-Man element. When the target was surrounded by flankers that differed by {+/-}30{degrees}, all three groups (adults and children with typical or amblyopic vision) reported orientations between the target and flankers (assimilation). Errors were reduced with {+/-}90{degrees} differences, but primarily matched the flanker orientation (substitution) when they did occur. A population pooling model of crowding successfully simulated this pattern of errors in all three groups. We conclude that the perceptual effects of amblyopic and developing crowding are systematic and resemble the near periphery in adults, suggesting a common underlying mechanism.

PrecisCrowding strongly disrupts peripheral vision, as well as the foveal vision of children with typical vision and amblyopia. We show that typically developing and amblyopic children make the same crowded errors as adults in the visual periphery, consistent with a common mechanism in all three cases.


title: NamZ1 and NamZ2 from the oral pathogen Tannerella forsythia are peptidoglycan processing exo-β-N-acetylmuramidases with distinct substrate specificity

abstract:
The Gram-negative periodontal pathogen Tannerella forsythia is inherently auxotrophic for N-acetylmuramic acid (MurNAc), which is an essential carbohydrate constituent of the peptidoglycan (PGN) of the bacterial cell wall. Thus, to build up its cell wall, T. forsythia strictly depends on the salvage of exogenous MurNAc or sources of MurNAc, such as polymeric or fragmentary PGN, derived from cohabiting bacteria within the oral microbiome. In our effort to elucidate how T. forsythia satisfies its demand for MurNAc, we recognized that the organism possesses three putative orthologs of the exo-{beta}-N-acetylmuramidase BsNamZ from Bacillus subtilis, which cleaves non-reducing end, terminal MurNAc entities from the artificial substrate pNP-MurNAc and the naturally-occurring disaccharide substrate MurNAc-{beta}-1,4-N-acetylglucosamine (GlcNAc). TfNamZ1 and TfNamZ2 were successfully purified as soluble, pure recombinant His6-fusions and characterized as exo-lytic {beta}-N-acetylmuramidases with distinct substrate specificities. The activity of TfNamZ1 was considerably lower compared to TfNamZ2 and BsNamZ, in the cleavage of pNP-MurNAc and MurNAc-GlcNAc. When peptide-free PGN glycans were used as substrates, we revealed striking differences in the specificity and mode of action of these enzymes, as analyzed by mass spectrometry. TfNamZ1, but not TfNamZ2 or BsNamZ, released GlcNAc-MurNAc disaccharides from these glycans. In addition, glucosamine (GlcN)-MurNAc disaccharides were generated when partially N-deacetylated PGN glycans from B. subtilis 168 were applied. This characterizes TfNamZ1 as a unique disaccharide-forming exo-lytic {beta}-N-acetylmuramidase (exo-disaccharidase), and, TfNamZ2 and BsNamZ as sole MurNAc monosaccharide-lytic exo-{beta}-N-acetylmuramidases.

IMPORTANCETwo exo-{beta}-N-acetylmuramidases from T. forsythia belonging to glycosidase family GH171 (www.cazy.org) were shown to differ in their activities, thus revealing a functional diversity within this family: NamZ1 releases disaccharides (GlcNAc-MurNAc/GlcN-MurNAc) from the non-reducing ends of PGN glycans, whereas NamZ2 releases terminal MurNAc monosaccharides. This work provides a better understanding of how T. forsythia may acquire the essential growth factor MurNAc by the salvage of PGN from cohabiting bacteria in the oral microbiome, which may pave avenues for the development of anti-periodontal drugs. On a broad scale, our study indicates that the utilization of PGN as a nutrient source, involving exo-lytic N-acetylmuramidases with different modes of action, appears to be a general feature of bacteria, particularly among the phylum Bacteroidetes.


title: Glucose-1,6-bisphosphate, a key metabolic regulator, is synthesized by a distinct family of α-D-phosphohexomutases widely distributed in prokaryotes

abstract: The reactions of -D-phosphohexomutases (PHM) are ubiquitous, key to primary metabolism and essential for several processes in all domains of life. The functionality of these enzymes relies on an initial auto-phosphorylation step which requires the presence of -D-glucose-1,6-bisphosphate (Glc-1,6-BP). While well investigated in vertebrates, the origin of this activator compound in bacteria is unknown. Here we show that the Slr1334 protein from the unicellular cyanobacterium Synechocysitis sp. PCC 6803 is a Glc-1,6-BP-synthase. Biochemical analysis revealed that Slr1334 efficiently converts fructose-1,6-bisphosphate (Frc-1,6-BP) and -D-glucose-1-phosphate/-D-glucose-6-phosphate into Glc-1,6-BP and also catalyzes the reverse reaction. As inferred from phylogenetic analysis, the slr1334 product belongs to a primordial subfamily of PHMs that is present especially in deeply branching bacteria and also includes human commensals and pathogens. Remarkably, the homologue of Slr1334 in the human gut bacterium Bacteroides salyersiae catalyzes the same reaction, suggesting a conserved and essential role for the members of this PHM subfamily.


title: FoQDE2-dependent milRNA promotes Fusarium oxysporum f. sp. cubense virulence by targeting a glycosyl hydrolase coding gene at transcriptional level

abstract:
MicroRNAs (miRNAs) are small non-coding RNAs that regulate protein-coding gene expression primarily found in plants and animals. Fungi produce microRNA-like RNAs (milRNAs) that are structurally similar to miRNAs and functionally important in various biological processes. The fungus Fusarium oxysporum f. sp. cubense (Foc) is the causal agent of Panama disease that threatens global banana production. It remains uncharacterized about the biosynthesis and functions of milRNAs in Foc. In this study, we investigated the biological function of milRNAs contributing to Foc pathogenesis. Within 24 hours post infecting the host, the Argonaute coding gene FoQDE2, and two Dicer coding genes FoDCL1 and FoDCL2, all of which are involved in milRNA biosynthesis, were significantly induced. FoQDE2 deletion mutant exhibited decreased virulence and hypersensitivity to hydrogen peroxide (H2O2). These results indicate that milRNA biosynthesis is crucial for Foc pathogenesis. By small RNA sequencing, we identified 364 small RNA-producing loci in the Foc genome, 25 of which were significantly downregulated in the FoQDE2 deletion mutant, from which milR-87 was verified as a FoQDE2-depedent milRNA based on qRT-PCR analysis. Through deletion and overexpression of milR-87 in the wild-type Foc strain, functions of milR-87 were studied. The results showed that milR-87 is crucial for Foc virulence in infection process. We furthermore identified a glycosyl hydrolase-coding gene, FOIG_15013, as the direct target of milR-87. The FOIG_15013 deletion mutant displayed a dramatic increase in the growth, conidiation and virulence. Transient expression of FOIG_15013 in Nicotiana benthamiana leaves activates the host defense responses. Collectively, this study documents the involvement of milRNAs in the manifestation of the devastating fungal disease in banana, and demonstrates the importance of milRNAs in the pathogenesis and other biological processes. Further analyses of the biosynthesis and expression regulation of fungal milRNAs may offer a novel strategy to combat devastating fungal diseases.

Author summaryThe fungus Fusarium oxysporum f. sp. cubense (Foc) is the causal agent of Panama disease that threatens global banana production. As a typical representative of F. oxysporum species complex, the pathogen has been widely concerned. However, pathogenesis of Foc is not fully elucidated. In particular, pathogenic regulatory mechanism of the microRNA like small RNAs (milRNAs) found in Foc is unknown. Here, we found that FoQDE2, one Argonaute coding gene, and two Dicer coding genes FoDCL1 and FoDCL2, which are involved in milRNA biosynthesis, are significantly induced during the early infection stage of Foc. The results suggested that the milRNAs biosynthesis mediated by these genes may play an active role in the infection process of Foc. Based on this assumption, we subsequently found a FoQDE2-dependent milRNA (milR-87) and identified its target gene. Functional analysis showed that FoQDE2, miR-87 and its target gene were involved in the pathogenicity of Foc in different degree. The studies help us gain insight into the pathogenesis with FoQDE2, milR-87, and its target gene as central axis in Foc. The identified pathogenicity-involved milRNA provides an active target for developing novel and efficient biocontrol agents against Panama disease.


title: Redox-active secondary metabolites act as interspecies modulators of antibiotic resilience

abstract: Bacterial opportunistic pathogens make a wide range of secondary metabolites both in the natural environment and when causing infections, yet how these molecules mediate microbial interactions and their consequences for antibiotic treatment are still poorly understood. Here, we explore the role of two redox-active secondary metabolites, pyocyanin and toxoflavin, as interspecies modulators of antibiotic resilience. We find that these molecules dramatically change susceptibility levels of diverse bacteria to clinical antibiotics. Pyocyanin is made by Pseudomonas aeruginosa, while toxoflavin is made by Burkholderia gladioli, organisms that infect cystic fibrosis and other immunocompromised patients. Both molecules alter the susceptibility profile of pathogenic species within the "Burkholderia cepacia complex" to different antibiotics, either antagonizing or potentiating their effects, depending on the drugs class. Defense responses regulated by the redox-sensitive transcription factor SoxR potentiate the antagonistic effects these metabolites have against fluoroquinolones, and the presence of genes encoding SoxR and the efflux systems it regulates can be used to predict how these metabolites will affect antibiotic susceptibility of different bacteria. Finally, we demonstrate that inclusion of secondary metabolites in standard protocols used to assess antibiotic resistance can dramatically alter the results, motivating the development of new tests for more accurate clinical assessment.


title: MgrB dependent colistin resistance in Klebsiella pneumoniae is associated with an increase in host-to-host transmission.

abstract:
Due to its high transmissibility, Klebsiella pneumoniae (Kpn) is one of the leading causes of nosocomial infections. Here, we studied the biological cost of colistin resistance, an antibiotic of last resort, of this opportunistic pathogen using a murine model of gut colonization and transmission. Colistin resistance in Kpn is commonly the result of inactivation of the small regulatory protein MgrB. Without a functional MgrB, the two-component system PhoPQ is constitutively active, leading to increased lipid A modifications and subsequent colistin resistance. Using an engineered MgrB mutant, we observed that MgrB-dependent colistin resistance is not associated with a fitness defect during in vitro growth conditions. However, colistin-resistant Kpn colonizes the murine gut poorly, which may be due to the decreased production of capsular polysaccharide by the mutant. The colistin-resistant mutant of Kpn had increased survival outside the host when compared to the parental colistin-sensitive strain. We attribute this enhanced survivability to dysregulation of the PhoPQ two-component system and accumulation of the master stress regulator RpoS. The enhanced survival of the colistin resistant strain may be a key factor in the observed rapid host-to-host transmission in our model. Together, our data demonstrate that colistin-resistant Kpn experiences a biological cost in gastrointestinal colonization. However, this cost is mitigated by enhanced survival outside the host, increasing the risk of transmission. Additionally, it underscores the importance of considering the entire life cycle of a pathogen to truly determine the biological cost associated with antibiotic resistance.

ImportanceThe biological cost associated with colistin resistance in Klebsiella pneumoniae (Kpn) was examined using a murine model of Kpn gut colonization and fecal-oral transmission. A common mutation resulting in colistin resistance in Kpn is a loss-of-function mutation of the small regulatory protein MgrB that regulates the two-component system PhoPQ. Even though colistin resistance in Kpn comes with a fitness defect in gut colonization, it increases bacterial survival outside the host enabling it to more effectively transmit to a new host. The enhanced survival is dependent upon the accumulation of RpoS and dysregulation of the PhoPQ. Hence, our study expands our understanding of the underlying molecular mechanism contributing to the transmission of colistin-resistant Kpn.


title: An ecological perspective on microbial genes of unknown function in soil

abstract: Genes that remain hypothetical, uncharacterized, and unannotated comprise a substantial portion of metagenomic datasets and are likely to be particularly prevalent in soils where poorly characterized taxa predominate. Documenting the prevalence, distribution, and potential roles of these genes of unknown function is an important first step to understanding their functional contributions in soil communities. We identified genes of unknown function from 50 soil metagenomes and analyzed their environmental distributions and ecological associations. We found that genes of unknown function are prevalent in soils, particularly fine-textured, higher pH soils that harbor greater abundances of Crenarchaeota, Gemmatimonadota, Nitrospirota, and Methylomirabilota. We identified 43 dominant (abundant and ubiquitous) gene clusters of unknown function and determined their associations with soil microbial phyla and other "known" genes. We found that these dominant unknown genes were commonly associated with microbial phyla that are relatively uncharacterized, with the majority of these dominant unknown genes associated with mobile genetic elements. This work demonstrates a strategy for investigating genes of unknown function in soils, emphasizes the biological insights that can be learned by adopting this strategy, and highlights specific hypotheses that warrant further investigation regarding the functional roles of abundant and ubiquitous genes of unknown function in soil metagenomes.


title: Global inflammatory response in in vitro organ cultured testes using single-cell RNA-sequencing

abstract: In vitro functional sperm production is important for understanding spermatogenesis and for the treatment of male infertility. Here, we describe similarities and differences between testis tissues in vivo and in vitro and clarify abnormalities in the early stage of in vitro spermatogenesis at single-cell resolution. While in vitro spermatogenesis progressed similarly to in vivo spermatogenesis until the early pachytene spermatocyte stage, a noticeable acute inflammatory response occurred in immune cells and non-immune testicular somatic cells immediately after cultivation. Inhibitor treatment revealed that NLRP3 inflammasome signaling is key to the inflammation. We observed damaged/dead germ cell accumulation in cultured testis, which may be due to dysfunctional phagocytosis by Sertoli cells. Our data revealed abnormal testicular milieu of in vitro cultured testes caused by tissue-wide sterile inflammation, in which the danger-associated molecular pattern-NLRP3 inflammasome axis may be a key element.

category: molecular biology


title: Differential mobility and self-association of Arc/Arg3.1 in the cytoplasm and nucleus of living cells

abstract: Arc, also known as Arg3.1, is an activity-dependent immediate-early gene product that plays essential roles in memory consolidation. A pool of Arc is located in the postsynaptic cytoplasm, where it promotes AMPA receptor endocytosis and cytoskeletal remodeling. However, Arc is also found in the nucleus, a major portion being associated with promyelocytic leukemia nuclear bodies (PML-NBs). Nuclear Arc has been implicated in epigenetic control of gene transcription associated with learning and memory. In this study, we use a battery of fluorescence nanoimaging approaches to characterize the behavior of Arc in living cells. Our results indicate that in the cytoplasm, Arc exists predominantly as monomers and dimers associated with slowly diffusing particles. In contrast, nuclear Arc is almost exclusively monomeric and displays a higher diffusivity than cytoplasmic Arc. We further show that Arc moves freely and rapidly between PML-NBs and the nucleoplasm, and that its movement within PML-NBs is relatively unobstructed.

category: cell biology


title: The HIV-1 capsid core is an opportunistic nuclear import receptor

abstract: The movement of viruses and other large macromolecular cargo through nuclear pore complexes (NPCs) is poorly understood. The human immunodeficiency virus type 1 (HIV-1) provides an attractive model to interrogate this process due to the genetic and cell biological assays to score virus nuclear entry in living cells. Although initial studies of HIV-1 infection of nondividing cells focused on karyophilic virion proteins, subsequent work revealed the viral capsid (CA), the chief structural component of the pre-integration complex (PIC), to be a critical determinant in nuclear transport1. In support of this model, HIV-1 interactions with NPCs can be altered through CA mutation2, which makes direct contact with nucleoporins (Nups)3-5. Here we identify Nup35, Nup153, and POM121 to coordinately support HIV-1 nuclear entry. For Nup35 and POM121, this dependence was strongly dependent cyclophilin A (CypA) interaction with CA. Mutation of CA or removal of soluble host factors changed the interaction with the NPC. Collectively, these findings implicate the HIV-1 CA hexameric lattice that encapsulates the viral genome as a macromolecular nuclear transport receptor (NTR) that exploits soluble host factors to modulate NPC requirements during nuclear invasion.


title: Cationic inhalable particles for enhanced drug delivery to M. tuberculosis infected macrophages

abstract: Inhalable microparticle-based drug delivery platforms are being investigated extensively for Tuberculosis (TB) treatment as they offer efficient deposition in lungs and improved pharmacokinetics of the encapsulated cargo. However, the effect of physical parameters of microcarriers on interaction with Mycobacterium tuberculosis (Mtb) infected mammalian cells is underexplored. In this study, we report that Mtb-infected macrophages are highly phagocytic and microparticle surface charge plays a major role in particle internalization by infected cells. Microparticles of different sizes (0.5 - 2 m) were internalized in large numbers by Mtb-infected THP-1 macrophages and murine primary Bone Marrow Derived Macrophages in vitro. Drastic improvement in particle uptake was observed with cationic particles in vitro and in mice lungs. Rapid uptake of rifampicin-loaded cationic microparticles allowed high intracellular accumulation of the drug and lead to enhanced anti-bacterial function when compared to non-modified rifampicin-loaded microparticles. Cytocompatibility assay and histological analysis in vivo confirmed that the formulations were safe and did not elicit any adverse reaction. Additionally, pulmonary delivery of cationic particles in mice resulted in two-fold higher uptake in resident alveolar macrophages compared to non-modified particles. This study provides a framework for future design of drug carriers to improve delivery of anti-TB drugs inside Mtb-infected cells.

category: bioengineering


title: Accurate sequence-dependent coarse-grained model for conformational and elastic properties of double-stranded DNA

abstract: We introduce MADna, a sequence-dependent coarse-grained model of double-stranded DNA (dsDNA), where each nucleotide is described by three beads localized at the sugar and base moieties, and at the phosphate group. The sequence dependence is included by considering a step-dependent parameterization of the bonded interactions, which are tuned in order to reproduce the values of key observables obtained from exhaustive atomistic simulations from literature. The predictions of the model are benchmarked against an independent set of all-atom simulations, showing that it captures with high fidelity the sequence dependence of conformational and elastic features beyond the single step considered in its formulation. A remarkably good agreement with experiments is found for both sequence-averaged and sequence-dependent conformational and elastic features, including the stretching and torsion moduli, the twist-stretch and twist-bend couplings, the persistence length and the helical pitch. Overall, for the inspected quantities, the model has a precision comparable to atomistic simulations, hence providing a reliable coarse-grained description for the rationalization of singlemolecule experiments and the study of cellular processes involving dsDNA. Owing to the simplicity of its formulation, MADna can be straightforwardly included in common simulation engines.

category: biophysics


title: The SPICA Coarse-Grained Force Field for Proteins and Peptides

abstract: A coarse-grained (CG) model for peptides and proteins was developed as an extension of the SPICA (Surface Property fItting Coarse grAined) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulation with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. The side chain analogs reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable dimerization free energies in water compared to those from AA-MD. A series of peptides/proteins adsorbed or penetrated into the membrane simulated by the CG-MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins comparable to those in the orientation of protein in membrane (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated a good performance of the SPICA FF.


title: AKT1 is Required for a Complete Palbociclib-induced Senescence Phenotype in BRAF-V600E-Driven Human Melanoma

abstract: Cellular senescence is a carefully regulated process of proliferative arrest accompanied by numerous functional and morphologic changes. Senescence allows damaged cells to avoid neoplastic proliferation, however induction of the senescence-associated secretory phenotype (SASP) can promote tumor growth. The complexity of the senescence response may limit the efficacy of anti-neoplastic agents, such as CDK4/6 inhibitors (Cdk4/6i), that induce a senescence-like, non-proliferative state in tumor cells. The AKT kinase family plays an important role in cellular growth and division, and is commonly hyperactive in many cancers including melanoma. AKT activity has also been implicated in regulation of senescence. The three AKT isoforms play both redundant and unique roles in tumorigenesis and cancer progression. To interrogate the role of AKT isoforms in the induction of cellular senescence by Cdk4/6i, we generated isoform specific AKT knockout human BRAF-V600E mutated melanoma cell lines. We found that the CDK4/6i Palbociclib induced a form of senescence in these cells that was dependent on AKT1. As a potential mechanism, we evaluated the activity of the cGAS-STING pathway, recently implicated in cellular senescence. While we showed cGAS-STING function to be dependent on AKT1, pharmacologic inhibition of either cGAS or STING had little effect on senescence. However, we found SASP factors to require NF-kB function, in part dependent on a stimulatory phosphorylation of IKK by AKT1 previously reported in other models. In summary, we provide the first evidence of a novel, isoform specific role for AKT1 in therapy-induced senescence in human melanoma cells acting through NF-kB but independent of cGAS-STING.


title: MDR1 Drug Efflux Pump Promotes Intrinsic and Acquired Resistance to PROTACs in Cancer Cells

abstract: PROTACs (Proteolysis-Targeting Chimeras) represent a promising new class of drugs that selectively degrade proteins of interest from cells. PROTACs targeting oncogenes are avidly being explored for cancer therapies, with several currently in clinical trials. Drug resistance represents a significant challenge in cancer therapies, and the mechanism by which cancer cells acquire resistance to PROTACs remains poorly understood. Using proteomics, we discovered acquired and intrinsic resistance to PROTACs in cancer cells can be mediated by upregulation of the drug efflux pump MDR1. PROTAC-resistant cells could be re-sensitized to PROTACs through co-administering MDR1 inhibitors. Notably, co-treatment of MDR1-overexpressing colorectal cancer cells with MEK1/2 or KRASG12C degraders and the dual ErbB receptor/MDR1 inhibitor lapatinib exhibited potent drug synergy due to simultaneous blockade of MDR1 and ErbB receptor activity. Together, our findings suggest that concurrent blockade of MDR1 will likely be required in combination with PROTACs to achieve durable protein degradation and therapeutic response in cancer.


title: A Novel Oncogenic Function of PRC2 Heterogeneity in Medulloblastoma

abstract: Intratumor epigenetic heterogeneity is emerging as a key mechanism underlying tumor evolution and drug resistance. Medulloblastomas, the most common childhood malignant brain tumor, are classified into four subtypes including SHH medulloblastomas, which are characterized by elevated SHH signaling and a cerebellum granule neuron precursor (CGNP) cell-of-origin. Medulloblastomas are highly associated with epigenetic abnormalities. We observed that the histone H3K27 methyltransferase polycomb repressor complex 2 (PRC2) is often heterogeneous within individual SHH medulloblastoma tumors. Using mouse models, we showed that while a complete deletion of the PRC2 core subunit EED inhibited medulloblastoma growth, a mosaic deletion of EED significantly enhanced tumor growth. EED is intrinsically required for CGNP maintenance by inhibiting both neural differentiation and cell death. Complete EED deletion led to CGNP depletion and reduced occurrence of medulloblastoma. Surprisingly, we found that medulloblastomas with mosaic EED levels grew faster than did control wildtype tumors and expressed increased levels of oncogenes such as Igf2. Igf2 is directly repressed by PRC2 and has been demonstrated to be both necessary and sufficient for SHH medulloblastoma progression. We showed that IGF2 mediated the oncogenic effects of PRC2 heterogeneity in tumor growth. Using a human medulloblastoma cell line, we generated clones with different EED levels and confirmed that EEDlow cells could stimulate the growth of EEDhigh cells through derepressed IGF2 signals. Thus, PRC2 heterogeneity controls medulloblastoma growth through both intrinsic growth competence and non-cell autonomous mechanisms in distinct tumor subclones. We reveal a novel oncogenic function of PRC2 heterogeneity in tumor development.


title: Microbial ecology of sulfur biogeochemical cycling at a mesothermic hot spring atop Northern Himalayas, India

abstract:
O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=146 SRC="FIGDIR/small/470874v1_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@1dfd0d5org.highwire.dtl.DTLVardef@107811forg.highwire.dtl.DTLVardef@1ae4af6org.highwire.dtl.DTLVardef@1bb90a1_HPS_FORMAT_FIGEXP  M_FIG Graphical Abstract C_FIG Sulfur Related Prokaryotes (SRP) residing in hot spring present good opportunity for exploring the limitless possibilities of integral ecosystem processes. Metagenomic analysis further expand the phylogenetic breadth of these extraordinary sulfur metabolizing microorganisms, as well a their complex metabolic networks and syntrophic interactions in environmental biosystems. Through this study, we explored and expanded the microbial genetic repertoire with focus on sulfur cycling genes through metagenomic analysis of sulfur (S) contaminated hot spring, located at the Northern Himalayas. The analysis revealed rich diversity of microbial consortia with established roles in S cycling such as Pseudomonas, Thioalkalivibrio, Desulfovibrio and Desulfobulbaceae (Proteobacteria). The major gene families inferred to be abundant across microbial mat, sediment and water were assigned to Proteobacteria as reflected from the RPKs (reads per kilobase) categorized into translation and ribosomal structure and biogenesis. Analysis of sequence similarity showed conserved pattern of both dsrAB genes (n=178) retrieved from all metagenomes while other sulfur disproportionation proteins were diverged due to different structural and chemical substrates. The diversity of sulfur oxidizing bacteria (SOB) and sulfate reducing bacteria (SRB) with conserved (r)dsrAB suggests for it to be an important adaptation for microbial fitness at this site. Here, we confirm that (i) SRBs belongs to{delta} -Proteobacteria occurring independent LGT of dsr genes to different and few novel lineages (ii) also, the oxidative and reductive dsr evolutionary time scale phylogeny, proved that the earliest (not first) dsrAB proteins belong to anaerobic Thiobacillus with other (rdsr) oxidizers. Further, the structural prediction of unassigned DsrAB proteins confirmed their relatedness with species of Desulfovibrio (TM score= 0.86; 0.98; 0.96) and Archaeoglobus fulgidus (TM score= 0.97; 0.98). We proposed that the genetic repertoire might provide the basis of studying time scale evolution and horizontal gene transfer of these genes in biogeochemical S cycling and the complementary genes could be implemented in biotechnology and bioremediation applications.


title: Bacterial community response to species overrepresentation or omission is strongly influenced by life in spatially structured habitats

abstract:
Variations in type and strength of interspecific interactions in natural bacterial communities (e.g., synergistic to inhibitory) affect species composition and community functioning. The extent of interspecific interactions is often modulated by environmental factors that constrain diffusion pathways and cell mobility and limit community spatial arrangement. We studied how spatially structured habitats affect interspecific interactions and influence the resulting bacterial community composition. We used a bacterial community made of 11 well-characterized species that grew in porous habitats (comprised of glass beads) under controlled hydration conditions or in liquid habitats. We manipulated the initial community composition by overrepresenting or removing selected members, and observed community composition over time. Life in porous media reduced the number and strength of interspecific interactions compared to mixed liquid culture, likely due to spatial niche partitioning in porous habitats. The community converged to similar species composition irrespective of the initial species mix, however, the dominant bacterial species was markedly different between liquid culture and structured porous habitats. Moreover, differences in water saturation levels of the porous medium affected community assembly highlighting the need to account for habitat structure and physical conditions to better understand and interpret assembly of bacterial communities. We point at the modulation of bacterial interactions due to spatial structuring as a potential mechanism for promoting community stability and species coexistence, as observed in various natural environments such as soil or human gut.

ImportanceBacteria live as complex multispecies communities essential for healthy and functioning ecosystems ranging from soil to the human gut. The bacterial species that form these communities can have positive or negative impact on each other, promoting or inhibiting each others growth. Yet, the factors controlling the balance of such interactions in nature, and how these influence the community, are not fully understood. Here, we show that bacterial interactions are modified by life in spatially structured bacterial habitats. These conditions exert important control over the resulting bacterial community regardless of initial species composition. The study demonstrates limitations of inferences from bacterial communities grown in liquid culture relative to behaviour in structured natural habitats such as soil.


title: Dynamic utilization of low-molecular-weight organic substrates across a microbial growth rate gradient

abstract:
Constantly in flux, low-molecular-weight organic substances (LMWOSs) are at the nexus between microorganisms, plant roots, detritus, and the soil mineral matrix. Nominal oxidation state of carbon (NOSC) has been put forward as one way to parameterize microbial uptake rates of LMWOSs and efficiency of carbon incorporation into new biomass. In this study, we employed an ecophysiological approach to test these proposed relationships using targeted exometabolomics (1H-NMR, HR-LCMS) coupled with stable isotope (13C) probing. We assessed the role of compound class and oxidation state on uptake kinetics and substrate-specific carbon use efficiency (SUE) during the growth of three model soil microorganisms (Penicillium spinulosum, Paraburkholderia solitsugae, and Ralstonia pickettii) in media containing 34 common LMWOSs. Microbial isolates were chosen to span a gradient in growth rate (0.046-0.316 hr-1) and differ phylogenetically (a fungal isolate and two bacterial isolates). Clustered, co-utilization of LMWOSs occured for all three organisms, but temporal cluster separation was most apparent for P. solitsugae. Potential trends (p <0.05) for early utilization of more oxidized substrates were present for the two bacterial isolates (P. solitsugae and R. pickettii), but high variability (R2 > 0.15) and a small effect of NOSC indicate these are not useful relationships for prediction. The SUEs ranged from 0.16-0.99 and the hypothesized inverse relationship between NOSC and SUE was not observed. Thus, our results do not provide compelling support for NOSC as a predictive tool, implying that metabolic strategies of organisms may be more important than chemical identity in determining LMWOS cycling in soils.

ImportanceCommunity-level observations from soils indicate that low-molecular-weight compounds of higher oxidation state tend to be depleted from soil solution faster and incorporated less efficiently into microbial biomass under oxic conditions. Here, we tested hypothetical relationships between substrate chemical characteristics and the order of substrate utilization by aerobic heterotrophs at the population-level in culture, using two bacterial isolates (Ralstonia pickettii and Paraburkholderia solitsugae) and one fungal isolate from soil (Penicillium spinulosum). We found weak relationships indicating earlier uptake of more oxidized substrates by the two bacterial isolates but no relationship for the fungal isolate. We found no relationship between substrate identity and substrate use efficiency. Our findings indicate that substrate chemical characteristics have limited utility for modeling the depletion of low-molecular-weight organics from soil solution and incorporation into biomass over broader phylogenetic gradients.


title: CRISPR gene-drive systems based on Cas9 nickases promote super-Mendelian inheritance in Drosophila

abstract: CRISPR-based gene-drive systems have been proposed for managing insect populations, including disease-transmitting mosquitoes, due to their ability to bias their inheritance towards super-Mendelian rates (>50%). Current technologies employ a Cas9 that introduces DNA double-strand breaks into the opposing wildtype allele to replace it with a copy of the gene drive allele via DNA homology-directed repair. Yet, the use of different Cas9s versions is unexplored, and alternative approaches could increase the available toolkit for gene-drive designs. Here, we report a novel gene-drive approach that relies on Cas9 nickases that generate staggered paired nicks in DNA to propagate the engineered gene-drive cassette. We show that generating 5 overhangs in the system yields efficient allelic conversion. The nickase gene-drive arrangement produces large, stereotyped deletions that are advantageous for targeting essential genes. Indeed, our nickase approach should expand the repertoire for gene-drive designs aimed at applications in mosquitoes and beyond.


title: A nickase Cas9 gene-drive system promotes super-Mendelian inheritance in Drosophila


title: MAL suppresses OSCC tumorigenesis by maintaining epithelial cell differentiation

abstract:
Oral squamous cell carcinoma (OSCC) is widely recognized as an optimal model for precise medicine guided molecular biomarkers of cancer, however, few clinical practices were applied till now. Based on the data from our own studies and published papers, it was found that the expression of MAL was significantly decreased in epithelial cancer as compared with normal tissues, and exhibited a opposite association with pathological grade. To study the molecular events related to deficiency of MAL during carcinogenesis, occurrence and development, a Mal knockout mouse model was constructed and consistently reproduced and bred. The Mal knockout mice are highly vulnerable to tumor induction by carcinogen of 4NQO, evidenced by their extremely earlier carcinogenesis, higher incidence, and more aggressive growth. Analysis of scRNA-seq data indicated that Mal knockout mice lost the ability in maintaining epithelial cell differentiation and get more prone to carcinogen with a remarkably higher incidence of epithelial malignancy. Further analyses identified putative co-functional genes of MAL, including DSG1, AQP3 and S100A8, which are key factors in maintaining epithelial cell differentiation. To conclude, the current study exhibits the clinical significance and explains the tumor suppressing function of MAL. The results also suggest the potential of MAL and its co-functional genes being biomarkers for designing the prevention and/or differentiation therapy strategies in OSCC.

SignificanceMAL is found to be strongly opposite with tumor pathological grade from clinical and in vivo studies in OSCC. We propose MAL and its co-functional genes, including DSG1, AQP3 and S100A8, as key factors in maintaining epithelial cell differentiation and are valuable targets for designing prevention and differentiation therapy strategies in OSCC.


title: HSC-independent definitive hematopoietic cells persist into adult life

abstract: The stem cell theory that all blood cells are derived from hematopoietic stem cell (HSC) is a central dogma in hematology. However, various types of blood cells are already produced from hemogenic endothelial cells (HECs) before the first HSCs appear at embryonic day (E)11 in the mouse embryo. This early blood cell production from HECs, called HSC-independent hematopoiesis, includes primitive and definitive erythromyeloid progenitors that transiently support fetal blood homeostasis until HSC-derived hematopoiesis is established. Lymphoid potential has traditionally been detected in the extra-embryonic yolk sac (YS) and/or embryos before HSC emergence, but the actual presence of lymphoid progenitors at this stage remains unknown. In addition, whether HSCs in the fetal liver are the main source of innate-like B-1a cells has been controversial. Here, using complementary lineage tracing mouse models, we show that HSC-independent multipotent progenitors (MPPs) and HSC-independent adoptive B-lymphoid progenitors persist into adult life. Furthermore, HSCs minimally contribute to the peritoneal B-1a cell pool; most B-1a cells are originated directly from ECs in the YS and embryo and HSC-independent for life. Our discovery of extensive HSC-independent MPP and B-lymphoid progenitors in adults attests to the complex blood developmental dynamics through embryo to adult that underpin the immune system and challenges the paradigm of HSC theory in hematology.


title: Structural Basis for an Unprecedented Enzymatic Alkylation in Cylindrocyclophane Biosynthesis

abstract: The cyanobacterial enzyme CylK assembles the cylindrocyclophane natural products by performing two unusual alkylation reactions, forming new carbon-carbon bonds between aromatic rings and secondary alkyl halide substrates. This transformation is unprecedented in biology and the structure and mechanism of CylK are unknown. Here, we report x-ray crystal structures of CylK, revealing a distinctive fusion of a Ca2+ binding domain and a {beta}-propeller fold. We use a mutagenic screening approach to locate CylKs active site at its domain interface, identifying two residues, Arg105 and Tyr473, that are required for catalysis. Anomalous diffraction datasets collected with bound bromide ions, a product analog, suggest these residues interact with the alkyl halide electrophile. Additional mutagenesis and molecular dynamics simulations implicates Asp440 and Glu374 in activating the nucleophilic aromatic ring. Bioinformatic analysis of CylK homologs from other cyanobacteria establishes that they conserve these key catalytic amino acids but they are likely associated with divergent reactivity and altered secondary metabolism. By gaining a molecular understanding of this unusual biosynthetic transformation, this work fills a gap in our understanding of how alkyl halides are activated and used by enzymes as biosynthetic intermediates, informing enzyme engineering, catalyst design, and natural product discovery.


title: Single-molecule biophysics experiments in silico: Towards a physical model of a replisome

abstract: The interpretation of single-molecule experiments is frequently aided by computational modeling of biomolecular dynamics. The growth of computing power and ongoing validation of computational models suggest that it soon may be possible to replace some experiments out-right with computational mimics. Here we offer a blueprint for performing single-molecule studies in silico using a DNA binding protein as a test bed. We demonstrate how atomistic simulations, typically limited to sub-millisecond durations and zeptoliter volumes, can guide development of a coarse-grained model for use in simulations that mimic experimental assays. We show that, after initially correcting excess attraction between the DNA and protein, qualitative consistency between several experiments and their computational equivalents is achieved, while additionally providing a detailed portrait of the underlying mechanics. Finally the model is used to simulate the trombone loop of a replication fork, a large complex of proteins and DNA.