title: Assessing neuromodulation effects of theta burst stimulation to the prefrontal cortex using TMS-evoked potentials

abstract: Theta burst stimulation (TBS), a form of repetitive transcranial magnetic stimulation (TMS), is capable of non-invasively modulating cortical excitability. TBS is gaining popularity as a therapeutic tool for psychiatric disorders such as depression, in which the dorsolateral prefrontal cortex (DLPFC) is the main therapeutic target. However, the neuromodulatory effects of TBS on prefrontal regions remain unclear. An emerging tool to assess neuromodulation in non-motor regions is concurrent transcranial magnetic stimulation and electroencephalography (TMS-EEG) to measure TMS-evoked potentials (TEPs). We assessed twenty-four healthy participants (13 males, mean age 25.2{+/-}9.9 years) following intermittent TBS, continuous TBS, and sham applied to the left DLPFC using a double-blinded crossover design. TEPs were obtained at baseline and 2-, 15-, and 30-min post-stimulation. Four TEP components (N40, P60, N100 and P200) were analysed using mixed effects repeated measures models (MRMM). Results indicate no significant effects for any assessed components (all p>.05). The largest effect size (Cohens d = -0.5) comparing iTBS and sham was obtained for the N100 component at 15 minutes post-stimulation. This result was in the same direction but smaller than found in previous studies, suggesting that the true effect size may be lower than previously reported. Accurate estimates of the effects sizes and inter-individual heterogeneity will critically inform clinical applications using TEPs to assess the neuromodulatory effects of TBS.

category: physiology


title: Regulation of the Molecular Shuttle Service: A Multiscale Model of Aquaporin 2 Recycling

abstract: The response of cells to their environment is driven by a variety of proteins and messenger molecules. In eukaryotes their distribution, which is regulated by a vesicular transport system, is important for a tight cellular response. The recycling of aquaporin 2 between membrane and storage region is a crucial part of the body water homeostasis and its malfunction can lead to Diabetes insipidus. To understand the regulation of this system, we aggregated pathways and mechanisms from literature and derived three models in a hypothesis-driven approach. Furthermore, we combined the models to a single system to gain insight into key regulatory mechanisms. To achieve this we developed a multiscale computational framework for the modeling and simulation of cellular systems. The analysis of the system rationalises that the compartmentalization of cAMP in renal principal cells is a result of the protein kinase A signalosome and can only occur if specific cellular components are observed in conjunction. Endocytotic and exocytotic processes are inherently connected and can be regulated by the same protein kinase A signal.

category: systems biology


title: SPIKE-TRIGGERED INTRASPINAL MICROSTIMULATION IMPROVES MOTOR PERFORANCE IN AN AMBULATORY RAT MODEL OF SPINAL CORD INJURY

abstract: The purpose of this study was to determine if spike-triggered intraspinal microstimulation (ISMS) results in improved motor performance in an ambulatory rat model of spinal cord injury (SCI). Experiments were carried out in adult male Sprague Dawley rats with 175 kdyn moderate T8 contusion injury. Rats were randomly assigned to one of two groups: Control or Activity Dependent Stimulation (ADS) therapy. Four weeks post-SCI, all rats were implanted with a recording electrode in the left hindlimb motor cortex and a fine-wire, custom-made stimulating electrode in the contralateral lumbar spinal cord. Intracortical and intraspinal microstimulation were used to find sites of similar hip representation areas, which were paired together for ADS therapy. In the ADS therapy group, spike-stimulus conditioning was administered for 4 hours/day, 4 days/week, for 4 weeks via a tethered cable in a testing chamber. During therapy sessions, single-unit spikes were discriminated in real time in the hindlimb motor cortex and used to trigger stimulation in the spinal cord ventral horn. The optimal stimulus intensity (50% ISMS movement threshold) and spike-stimulus delay (10ms) determined in preliminary anesthetized preparations were used during ADS. Control rats were similarly implanted with electrodes but did not receive stimulation therapy. Motor performances of each rat were evaluated before SCI contusion, once a week post-SCI for four weeks (prior to electrode implantation), and once a week post-conditioning for four weeks. Behavioral testing included BBB scoring, Ledged Beam walking, Horizontal Ladder walking, treadmill kinematics via the DigiGait and TreadScan system, and open field walking using OptiTrack kinematic analysis. BBB scores were significantly improved in ADS rats compared to Control rats after 1 week of therapy. In the ADS therapy rats, BBB scores were significantly improved after two weeks of ADS therapy when compared to pre-therapy. Foot fault scores on the Horizontal Ladder were significantly lower in ADS rats compared to pre-therapy ADS and Control rats after 1 week of therapy and returned to pre-injury measures after three weeks of ADS therapy. The Ledged Beam test and kinematic analysis using the DigiGait and TreadScan system showed deficits after SCI in both ADS and Control rats but there were no significant differences between groups after 4 weeks of ADS therapy. These results show that activity dependent stimulation after spinal cord injury using spike-triggered ISMS enhances behavioral recovery of locomotor function as measured by the BBB score and the Horizontal Ladder task.

category: neuroscience


title: SPIKE-TRIGGERED INTRASPINAL MICROSTIMULATION IMPROVES MOTOR PERFORMANCE IN AN AMBULATORY RAT MODEL OF SPINAL CORD INJURY


title: Comparison of undirected frequency-domain connectivity measures for cerebro-peripheral analysis

abstract: Analyses of cerebro-peripheral connectivity aim to quantify ongoing coupling between brain activity (measured by MEG/EEG) and peripheral signals such as muscle activity, continuous speech, or physiological rhythms (such as pupil dilation or respiration). Due to the distinct rhythmicity of these signals, undirected connectivity is typically assessed in the frequency domain. This leaves the investigator with two critical choices, namely a) the appropriate measure for spectral estimation (i.e., the transformation into the frequency domain) and b) the actual connectivity measure. As there is no consensus regarding best practice, a wide variety of methods has been applied. Here we systematically compare combinations of six standard spectral estimation methods (comprising fast Fourier and continuous wavelet transformation, bandpass filtering, and short-time Fourier transformation) and six connectivity measures (phase-locking value, Gaussian-Copula mutual information, Rayleigh test, weighted pairwise phase consistency, magnitude squared coherence, and entropy). We provide performance measures of each combination for simulated data (with precise control over true connectivity), a single-subject set of real MEG data, and a full group analysis of real MEG data. Our results show that, overall, wppc and gcmi tend to outperform other connectivity measures, while entropy was the only measure sensitive to bimodal deviations from a uniform phase distribution. For group analysis, choosing the appropriate spectral estimation method appeared to be more critical than the connectivity measure. We discuss practical implications (sampling rate, SNR, computation time, and data length) and aim to provide recommendations tailored to particular research questions.


title: Computational neural network provides naturalistic solution for recovery of finger dexterity after stroke

abstract: Finger dexterity is a fundamental movement skill of humans and the ability to individuate fingers imparts high motor flexibility. Disruption of dexterity due to brain injury reduces quality of life. Thus, understanding the neurological mechanisms responsible for recovery is critical to effective neurorehabilitation. Two neuronal pathways have been proposed to play crucial roles in finger individuation: the corticospinal tract, originating from primary motor cortex and premotor areas, and the subcortical reticulospinal tract, originating from the reticular formation in the brainstem. Finger individuation in patients with lesions to these pathways may recover. However, it remains an open question how the cortical-reticular network reorganizes and contributes to this recovery following a stroke. We hypothesized that interactive connections between cortical and subcortical neurons reflect dynamics appropriate for generating outgoing commands for finger movement. To test this hypothesis, we developed an Artificial Neural Network (ANN) representing a premotor planning input layer, a cortical layer including excitatory and inhibitory neurons and, a reticular layer that control motoneurons eliciting unilateral flexion of two fingers. The ANN was trained to reproduce "normal" activity of finger individuation and strength. Analysis of the trained ANN revealed that the natural dynamical solution was a near-linear relationship between the force of the instructed and uninstructed finger, resembling individuation patterns in humans. A simulated stroke lesion was then applied to the ANN and the resulting finger dexterity was assessed at multiple stages post stroke. Analysis revealed: (1) increased unintended force produced by uninstructed fingers (i.e., enslaving) and (2) weakening of the force in the instructed finger immediately after stroke, (3) improved finger control during recovery that typically occurs early after stroke, and (4) association of this behavior with increased neural plasticity of the residual neurons, as reflected by strengthening of connectivity weights between premotor and focal cortical excitatory and inhibitory neurons, but reduction in connectivity in shared cortical neurons. Interestingly, the network solution predicted that the reticulospinal pathway also contributed to the improved behavior. Lastly, the ANN also predicts the effect of cortical lesion size on finger individuation. Our model provides a framework by which to understand a number of experimental findings. The model solution suggests that a key mechanism of finger individuation is establishment of an interactive relationship between cortical and subcortical regions, appropriate to produce desired finger movement.


title: Fate transitions in Drosophila neural lineages: a single cell road map to mature neurons.

abstract: Neuron specification and maturation are essential for proper central nervous system development. However, the precise mechanisms that govern neuronal maturation remain poorly understood. Here, we use single-cell RNA sequencing, combined with a conditional genetic strategy to analyse neuronal lineages and their new born neurons in the Drosophila larval brain. A focused analysis on the transcriptional alterations that occur right after neuron generation revealed that neuron maturation starts shortly after neuronal birth, with transcription, but no translation, of mature neuronal features such as neurotransmitter (NT) genes. Using NT gene Choline acetyltransferase as an example, we show that the timings of translation initiation are not solely dependent on neuron age but are rather coordinated with the animal developmental stage. This study is the first characterization of the initial phases of neuron maturation, supporting a model where neuron maturation is coordinated with the animal developmental stage through post-transcriptional regulation of terminal effector genes.


title: Fatty acid dysregulation in the anterior cingulate cortex of depressed suicides with a history of child abuse

abstract: Child abuse (CA) strongly increases the lifetime risk of suffering from major depression and predicts an unfavorable course for the illness. Severe CA has been associated with a specific dysregulation of oligodendrocyte function and thinner myelin sheaths in the human anterior cingulate cortex (ACC) white matter. Given that myelin is extremely lipid-rich, it is plausible that these findings may be accompanied by a disruption of the lipid profile that composes the myelin sheath. This is important to explore since the composition of fatty acids (FA) in myelin phospholipids can influence its stability, permeability, and compactness. Therefore, the objective of this study was to quantify and compare FA concentrations in postmortem ACC white matter in the choline glycerophospholipid pool (ChoGpl), a key myelin phospholipid pool, between adult depressed suicides with a history of CA (DS-CA) matched depressed suicides without CA (DS) and healthy non-psychiatric controls (CTRL). Total lipids were extracted according to the Folch method and separated into respective classes using thin-layer chromatography. FA methyl esters from the ChoGpl fraction were quantified using gas chromatography. Our analysis revealed a strong age-related decrease in most FAs, and specific effects of CA in FAs from the arachidonic acid synthesis pathway, which was further validated with RNA-sequencing data. Furthermore, the concentration of most FAs was found to decrease with age. By extending the previous molecular level findings linking CA with altered myelination in the ACC, these results provide further insights regarding white matter alterations associated with early-life adversity.


title: Timing of transcriptomic peripheral blood mononuclear cell responses of sheep to Fasciola hepatica infection differs from those of cattle, reflecting different disease phenotypes

abstract:
Infection with the zoonotic trematode Fasciola hepatica, common in many regions with a temperate climate, leads to delayed growth and loss of productivity in cattle, while infection in sheep can have more severe effects, potentially leading to death. Previous transcriptomic analyses revealed upregulation of TGFB1, cell death and Toll-like receptor signalling, T-cell activation, and inhibition of nitric oxide production in macrophages in response to infection. However, the differences between ovine and bovine responses have not yet been explored. The objective of this study was to further investigate the transcriptomic response of ovine peripheral blood mononuclear cells (PBMC) to F. hepatica infection, and to elucidate the differences between ovine and bovine PBMC responses.

Sixteen male Merino sheep were randomly assigned to infected or control groups (n = 8 per group) and orally infected with 120 F. hepatica metacercariae. Transcriptomic data was generated from PBMC at 0, 2 and 16 weeks post-infection (wpi), and analysed for differentially expressed (DE) genes between infected and control animals at each time point (analysis 1), and for each group relative to time 0 (analysis 2). Analysis 2 was then compared to a similar study performed previously on bovine PBMC.

A total of 453 DE genes were found at 2 wpi, and 2 DE genes at 16 wpi (FDR < 0.1, analysis 1). Significantly overrepresented biological pathways at 2 wpi included role of PKR in interferon induction and anti-viral response, death receptor signalling and RIG-I-like receptor signalling, which suggested that an activation of innate response to intracellular nucleic acids and inhibition of cellular apoptosis were taking place. Comparison of analysis 2 with the previous bovine transcriptomic study revealed that anti-inflammatory response pathways which were significantly overrepresented in the acute phase in cattle, including IL-10 signalling, Th2 pathway, and Th1 and Th2 activation were upregulated only in the chronic phase in sheep. We propose that the earlier activation of anti-inflammatory responses in cattle, as compared with sheep, may be related to the general absence of acute clinical signs in cattle. These findings offer scope for "smart vaccination" strategies for this important livestock parasite.

category: immunology


title: Inflammation induces pro-NETotic neutrophils via TNFR2 signaling

abstract: Cytokines released during chronic inflammatory diseases induce pro-inflammatory properties in polymorphonuclear neutrophils (PMN). Here we show that in vitro cytokine treatment leads to the development of a subgroup of human PMN expressing CCR5, termed CCR5+ cytokine-induced PMN (CCR5+ cPMN). Auto/paracrine TNF signaling increases intracellular neutrophil elastase (ELANE) abundance and induces NETosis in CCR5+ cPMN. Triggering of CCR5 amplifies NETosis. Membranous TNF (mTNF) outside-in signaling induces the formation of reactive oxygen species, a known activator of NETosis. In vivo, we find an increased number of CCR5+ cPMN in the peripheral blood and inflamed lamina propria of patients with ulcerative colitis (UC) but not Crohns disease (CD). Notably, failure of anti-TNF therapy is associated with higher frequencies of CCR5+ cPMN. In conclusion, we identify a phenotype of pro-NETotic, CCR5 positive PMN present in inflamed tissue in vivo and inducible in vitro. These cells may reflect an important component of tissue damage during chronic inflammation and could be of diagnostic value.


title: Linking plant growth promoting arbuscular mycorrhizal colonization with bacterial plant sulfur supply

abstract:
Sulfur (S) exists in organically bound complexes ([~]95%), predominantly as sulfonates, and are not directly plant available. Specific soil bacteria can mobilise sulfonates but very little is known about these bacteria in the hyphosphere. Since mycorrhizal fungi support growth of the majority of land plants, hyphosphere desulfonating bacteria may be of substantial benefit to the plant host. This study analysed the effect of AM inoculation with Rhizophagus irregularis (former G. intraradices, Glomus) and a mix of six AM species (Mixed) on PGP, microbial communities and sulfonate mobilising bacteria with L. perenne, Agrostis stolonifera and Plantago lanceolata as plant hosts in bi-compartmental microcosms and A. stolonifera in PGP pot experiments. AM inoculation significantly increased plant growth, percentage root colonisation and the quantity of cultivable desulfonating bacteria in the hyphosphere over pre-inoculated soil for all plants. Community analysis via PCR-DGGE revealed significantly different bacterial and fungal communities post inoculation. Analysis of the sulfonate mobilising asfA gene revealed a significantly altered community and novel bacterial isolates with this important functional ability post-inoculation. The results demonstrate that AM inoculation increased plant biomass yield, AM root colonisation and altered bacterial and fungal community dynamics in the hyphosphere. AM inoculated microcosms had an increased abundance of desulfonating bacteria that may be beneficial for plant-S supply.

Research highlightsO_LIInoculation with AM fungi was shown to promote plant growth and harbour larger populations of sulfonate mobilising bacteria.
C_LIO_LIPost-inoculation hyphospheric bacterial and saprotrophic fungal communities were shown to differ significantly in composition and abundance.
C_LIO_LIAnalysis of sulfonate mobilising bacteria revealed novel presumptive species in possession of the asfA gene associated with AM hyphae.
C_LIO_LIAM inoculation was shown to significantly impact the asfA positive bacterial community composition.
C_LI

category: microbiology


title: Classical swine fever virus NS5A inhibits NF-κB signaling by targeting NEMO

abstract:
The NS5A non-structural protein of classical swine fever virus (CSFV) is a multifunctional protein involved in viral genomic replication, protein translation and regulation of cellular signaling pathways, and assembly of infectious virus particles. Previous report showed that NS5A inhibited nuclear factor kappa B (NF-{kappa}B) signaling induced by poly(I:C); however, the mechanism was not elucidated. Here, we report that NS5A interacts with NF-{kappa}B essential modulator (NEMO),a regulatory subunit of the I{kappa}B kinase (IKK) complex, and that the zinc finger domain of NEMO essential for NEMO ubiquitination and IKK activation is required for the interaction of NEMO with NS5A. Viral infection or NS5A expression by itself reduced the protein level of NEMO. Mechanistic analysis revealed that NS5A mediated proteasomal degradation of NEMO. Ubiquitination assay showed that NS5A induced K27-but not K48-linked polyubiquitination of NEMO. In addition, NS5A blocked k63-linked polyubiquitination of NEMO, thereby inhibiting activation of IKK and NF-{kappa}B. These findings indicate that NS5A inhibits NF-{kappa}B signaling by mediating proteasomal degradation of NEMO and blocking k63-linked polyubiquitination of NEMO, thereby revealing a novel mechanism by which CSFV inhibits host innate immunity.

ImportanceClassical swine fever (CSF) is a economically important swine viral disease leading to hemorrhagic fever and immuno-suppression. In order to successfully infect and replicate in a host cell, viruses have evolved various strategies to antagonize host innate immunity. It is known that CSFV non-structural protein Npro interacts with interferon regulatory factor 3 (IRF3) and mediates its proteasomal degradation, thereby inhibiting the production of type I interferon. However, no other mechanism by which CSFV antagonizes host innate immunity has so far been reported. Here, we show that NS5A inhibits NF-{kappa}B signaling by mediating proteasomal degradation of NEMO and by blocking k63-linked polyubiquitination of NEMO, thereby revealing a novel mechanism by which CSFV antagonizes host innate immunity.


title: Mucosal immunization with DTaP confers protection against Bordetella pertussis infection and cough in Sprague-Dawley rats

abstract: Pertussis is a respiratory disease caused by the Gram-negative pathogen, Bordetella pertussis (Bp). The transition from a whole cell pertussis vaccine (wP; DTP) to an acellular pertussis vaccine (aP; DTaP; Tdap) correlates with an increase in pertussis cases, despite widespread vaccine implementation and coverage, and it is now appreciated that the protection provided by aP rapidly wanes. To recapitulate the localized immunity observed from natural infection, mucosal vaccination with aP was explored using the coughing rat model of pertussis. Immunity induced by both oral gavage (OG) and intranasal (IN) vaccination of aP in Bp challenged rats over a nine-day infection was compared to intramuscular (IM)-wP and IM-aP immunized rats that were used as positive controls as IM immunization is the current route for wP and aP vaccination. Our data demonstrate that both IN and OG immunization of aP resulted in production of anti-Bp IgG antibody titers similar to IM-wP and IM-aP vaccinated controls post-challenge. IN-aP also induced anti-Bp IgA antibodies in the nasal cavity. Immunization with IM-wP, IM-aP, IN-aP, and OG-aP immunization protected against Bp induced cough, while OG-aP immunization did not protect against respiratory distress. Mucosal immunization (IN-aP and OG-aP) also protected against acute inflammation and decreased bacterial burden in the lung compared to mock vaccinated challenge (MVC) rats. The data presented in this study suggests that mucosal vaccination with aP can induce a mucosal immune response and provide protection against Bp challenge.


title: Impaired succinate oxidation prevents growth and influences drug susceptibility in Mycobacterium tuberculosis

abstract:
Succinate is a major focal point in mycobacterial metabolism and respiration, serving as both an intermediate of the TCA cycle and a direct electron donor for the respiratory chain. Mycobacterium tuberculosis encodes multiple enzymes predicted to be capable of catalyzing the oxidation of succinate to fumarate, including two different succinate dehydrogenases (Sdh1 and Sdh2) and a separate fumarate reductase (Frd) with possible bi-directional behavior. Previous attempts to investigate the essentiality of succinate oxidation in M. tuberculosis have relied on the use of single-gene deletion mutants, raising the possibility that the remaining enzymes could catalyze succinate oxidation in the absence of the other. To address this, we report on the use of mycobacterial CRISPR interference (CRISPRi) to construct single, double, and triple transcriptional knockdowns of sdhA1, sdhA2, and frdA in M. tuberculosis. We show that the simultaneous knockdown of sdhA1 + sdhA2 is required to prevent succinate oxidation and overcome the functional redundancy within these enzymes. Succinate oxidation was demonstrated to be essential for the optimal growth of M. tuberculosis, with the combined knockdown of sdhA1 + sdhA2 significantly impairing the activity of the respiratory chain and preventing growth on a range of carbon sources. Moreover, impaired succinate oxidation was shown to influence the activity of several antitubercular drugs against M. tuberculosis, including potentiating the activity of bioenergetic inhibitors and attenuating the activity of cell wall inhibitors. Together, these data provide fundamental insights into mycobacterial physiology, energy metabolism, and antimicrobial susceptibility.

ImportanceNew drugs are urgently required to combat the tuberculosis epidemic that claims 1.5 million lives annually. Inhibitors of mycobacterial energy metabolism have shown significant promise clinically; however, further advancing this nascent target space requires a more fundamental understanding of the respiratory enzymes and pathways used by Mycobacterium tuberculosis. Succinate is a major focal point in mycobacterial metabolism and respiration; yet the essentiality of succinate oxidation, and the consequences of inhibiting this process, are poorly defined. In this study, we demonstrate that impaired succinate oxidation prevents the optimal growth of M. tuberculosis on a range of carbon sources and significantly reduces the activity of the electron transport chain. Moreover, we show that impaired succinate oxidation both positively and negatively influences the activity of a variety of anti-tuberculosis drugs. Combined, these findings provide fundamental insights into mycobacterial physiology and drug susceptibility that will be useful in the continued development of bioenergetic inhibitors.


title: Pickering stabilization of a dynamic intracellular emulsion

abstract:
Biomolecular condensates are cellular compartments that form by phase separation in the absence of limiting membranes. Studying the P granules of C. elegans, we find that condensate dynamics are regulated by protein clusters that adsorb to the condensate interface. Using in vitro reconstitution, live observations and theory, we demonstrate that localized assembly of P granules is controlled by MEG-3, an intrinsically disordered protein that forms low dynamic assemblies on P granules. Following classic Pickering emulsion theory, MEG-3 clusters lower surface tension and slow down coarsening. During zygote polarization, MEG-3 recruits DYRK/MBK-2 kinase to accelerate localized growth of the P granule emulsion. By tuning condensate-cytoplasm exchange, interfacial clusters regulate the structural integrity of biomolecular condensates, reminiscent of the role of lipid bilayers in membrane-bound organelles.

One Sentence SummaryBiomolecular condensates are stabilized by interfacial nanoscale protein clusters.

category: molecular biology


title: Visualizing sequential compound fusion and kiss-and-run in live excitable cells

abstract: Vesicle fusion is assumed to occur at flat membrane of excitable cells. In live neuroendocrine cells, we visualized vesicle fusion at {Omega}-shape membrane generated by preceding fusion, termed sequential compound fusion, which may be followed by fusion pore closure, termed compound kiss-and-run. These novel fusion modes contribute to vesicle docking, multi-vesicular release, asynchronous release, and endocytosis. We suggest modifying current models of exo-endocytosis to include these new fusion modes.

category: cell biology


title: Inhibition of N-myristoyltransferase Promotes Naive Pluripotency in Mouse and Human Pluripotent Stem Cells

abstract: Naive and primed states are distinct states of pluripotency during early embryonic development that can be captured and converted to each other in vitro. To elucidate the regulatory mechanism of pluripotency, we performed a recessive genetic screen of homozygous mutant mouse embryonic stem cells (mESCs) and found that suppression of N-myristoyltransferase (Nmt) promotes naive pluripotency. Disruption of Nmt1 in mESCs conferred resistance to differentiation. Suppression of Nmt in mouse epiblast stem cells (mEpiSCs) promoted the conversion from the primed to the naive state. This effect was independent of Src, which is a major substrate of Nmt and is known to promote differentiation of mESCs. Suppression of Nmt in naive-state human induced pluripotent stem cells (hiPSCs) increased the expression of the naive-state marker. These results indicate that Nmt is a novel target for the regulation of naive pluripotency conserved between mice and humans.

category: developmental biology


title: Cell adhesions link subcellular actomyosin dynamics to tissue scale force production during vertebrate convergent extension.

abstract: Axis extension is a fundamental biological process that shapes multicellular organisms. The design of an animals body plan is encoded in the genome and execution of this program is a multiscale mechanical progression involving the coordinated movement of proteins, cells, and whole tissues. Thus, a key challenge to understanding axis extension is connecting events that occur across these various length scales. Here, we use approaches from proteomics, cell biology, and tissue biomechanics to describe how a poorly characterized cell adhesion effector, the Armadillo Repeat protein deleted in Velo-Cardio-Facial syndrome (Arvcf) catenin, controls vertebrate head-to-tail axis extension. We find that Arvcf catenin is required for axis extension within the intact organism but is not required for extension of isolated tissues. We then show that the organism scale phenotype is caused by a modest defect in force production at the tissue scale that becomes apparent when the tissue is challenged by external resistance. Finally, we show that the tissue scale force defect results from dampening of the pulsatile recruitment of cell adhesion and cytoskeletal proteins to cell membranes. These results not only provide a comprehensive understanding of Arvcf function during an essential biological process, but also provide insight into how a modest cellular scale defect in cell adhesion results in an organism scale failure of development.


title: H3K9 tri-methylation at Nanog times differentiation commitment and enables the acquisition of primitive endoderm fate

abstract: Mouse Embryonic Stem (ES) cells have an inherent propensity to explore distinct gene-regulatory states associated with either self-renewal or differentiation. This property is largely dependent on ERK activity, which promotes silencing of pluripotency genes, most notably of the transcription factor Nanog. Here, we aimed at identifying repressive histone modifications that would mark the Nanog locus for inactivation in response to ERK activity. We found histone H3 lysine 9 tri-methylation (H3K9me3) focally enriched between the Nanog promoter and its -5kb enhancer. While in undifferentiated ES cells H3K9me3 at Nanog depends on ERK activity, in somatic cells it becomes ERK-independent. Moreover, upon deletion of the region harbouring H3K9me3, ES cells display reduced heterogeneity of NANOG expression, delayed commitment into differentiation and impaired ability to acquire a primitive endoderm fate. We suggest that establishment of irreversible H3K9me3 at specific master regulators allows the acquisition of particular cell fates during differentiation.


title: A mechanistic density functional theory for ecology across scales

abstract: Our ability to predict the properties of a system typically diminishes as the number of its interacting constituents rises. This poses major challenges for understanding natural ecosystems, and humanitys effects on them. How do macroecological patterns emerge from the interplay between species and their environment? What is the impact on complex ecological systems of human interventions, such as extermination of large predators, deforestation, and climate change? The resolution of such questions is hampered in part by the lack of a holistic approach that unifies ecology across temporal and spatial scales. Here we use density functional theory, a computational method for many-body problems in physics, to develop a novel computational framework for ecosystem modelling. Our methods accurately fit experimental and synthetic data of interacting multi-species communities across spatial scales and can project to unseen data. Our mechanistic framework provides a promising new avenue for understanding how ecosystems operate and facilitates quantitative assessment of interventions.

category: ecology


title: Climatic niche conservatism shapes the ecological assembly of Hawaiian arthropod communities

abstract: Spatial variation in climatic conditions along elevation gradients provides an important backdrop by which communities assemble and diversify. Lowland habitats tend to be connected through time, whereas highlands can be continuously or periodically isolated, conditions that have been hypothesized to promote high levels of species endemism. This tendency is expected to be accentuated among taxa that show niche conservatism within a given climatic envelope. While species distribution modeling approaches have allowed extensive exploration of niche conservatism among target taxa, a broad understanding of the phenomenon requires sampling of entire communities. Species-rich groups such as arthropods are ideal case studies for understanding ecological and biodiversity dynamics along elevational gradients given their important functional role in many ecosystems, but community-level studies have been limited due to their tremendous diversity. Here, we develop a novel semi-quantitative metabarcoding approach that combines specimen counts and size-sorting to characterize arthropod community-level diversity patterns along two elevational gradients across two volcanoes on the island of Hawai i. We find that arthropod communities between the two transects become increasingly distinct compositionally at higher elevations. Resistance surface approaches suggest that climatic differences between sampling localities are an important driver in shaping beta-diversity patterns, though the relative importance of climate varies across taxonomic groups. Nevertheless, the climatic niche position of OTUs between transects was highly correlated, suggesting that climatic filters shape the colonization between adjacent volcanoes. Taken together, our results highlight climatic niche conservatism as an important factor shaping ecological assembly along elevational gradients and suggest topographic complexity as an important driver of diversification.


title: Testing concordance and conflict in spatial replication of landscape genetics inferences

abstract: The field of landscape genetics relates habitat features and genetic information to infer dispersal and genetic connectivity between populations or individuals distributed across a landscape. Such studies usually focus on a small portion of a species range, and the degree to which these geographically restricted results can be extrapolated to different areas of a species range remains poorly understood. Studies that have focused on spatial replication in landscape genetics processes either evaluate a small number of sites, are informed by a small set of genetic markers, analyze only a small subset of environmental variables, or implement models that do not fully explore parameter space. Here, we used a broadly distributed ectothermic lizard (Sceloporus occidentalis, Western Fence lizard) as a model species to evaluate the full role of topography, climate, vegetation, and roads on dispersal and genetic differentiation. We conducted landscape genetics analyses in five areas within the Sierra Nevada mountain range, using thousands of ddRAD genetic markers distributed across the genome, implemented in the landscape genetics program ResistanceGA. Across study areas, we found a great deal of consistency in the variables impacting genetic connectivity, but also noted site-specific differences in the factors in each study area. High-elevation colder areas were consistently found to be barriers to gene flow, as were areas of high ruggedness and slope. High temperature seasonality and high precipitation during the winter wet season also presented a substantial barrier to gene flow in a majority of study areas. The effect of other landscape variables on genetic differentiation was more idiosyncratic and depended on specific attributes at each site. Vegetation type was found to substantially affect gene flow only in the southernmost Sequoia site, likely due to a higher proportion of desert habitat here, thereby fragmenting habitats that have lower costs to dispersal. The effect of roads also varied between sites and may be related to differences in road usage and amount of traffic in each area. Across study areas, canyons were always substantially implicated as facilitators to dispersal and key features linking populations and maintaining genetic connectivity across landscapes. We emphasize that spatial data layers are complex and multidimensional, and a careful consideration of associations between variables is vital to form sound conclusions about the critical factors affecting dispersal and genetic connectivity across space.


title: Extracting Phylogenetic Information of Human Mitochondrial DNA by Linear Autoencoder

abstract: We used a linear autoencoder (LAE) and its learning dynamics to analyze the high-order structure of human mitochondrial DNA (mtDNA). A total of 360 complete human mtDNA sequences were collected from the MITOMAP database and transformed into 1024-dimensional vectors of pentanucleotide frequencies. We compressed those into a three-dimensional (3D) coordinates by an LAE at each step of training by gradient descent with respect to the quadratic error function. Along the time axis of training epochs, the compressed 3D coordinates were gradually clustered and separated in accordance with the order of the genetic distance in the phylogenetic tree of human mtDNA haplogroups. This suggests that there is an association between the learning dynamics of LAE and the high-dimensional structure of human mtDNA sequences, similar to that of phylogenetic analysis and evolutionary pathways: the five clusters eventually contained only a single haplogroup of L0, M, N, R, and U, while the L3 cluster contained a small number of M members and The packing was comparable to that realized in learning dynamics similar to genetic classification and evolutionary pathways by LAE in principal component analysis (PCA), but somewhat denser than PCA.

category: evolutionary biology


title: Recovery of deleted deep sequencing data sheds more light on the early Wuhan SARS-CoV-2 epidemic

abstract: The origin and early spread of SARS-CoV-2 remains shrouded in mystery. Here I identify a data set containing SARS-CoV-2 sequences from early in the Wuhan epidemic that has been deleted from the NIHs Sequence Read Archive. I recover the deleted files from the Google Cloud, and reconstruct partial sequences of 13 early epidemic viruses. Phylogenetic analysis of these sequences in the context of carefully annotated existing data further supports the idea that the Huanan Seafood Market sequences are not fully representative of the viruses in Wuhan early in the epidemic. Instead, the progenitor of currently known SARS-CoV-2 sequences likely contained three mutations relative to the market viruses that made it more similar to SARS-CoV-2s bat coronavirus relatives.


title: A conserved trans regulatory loop involving an odorant binding protein controls male mating behavior in flies

abstract: A major goal in evolutionary biology is to understand how natural variation is maintained in sexually selected and sexually dimorphic traits. Hypotheses to explain genetic variation in sexually selected traits include context-dependent fitness effects, epistatic interactions, and pleiotropic constraints. The house fly, Musca domestica, is a promising system to investigate how these factors affect polymorphism in sexually selected traits. Two common Y chromosomes (YM and IIIM) segregate as stable polymorphisms in natural house fly populations, appear to be locally adapted to different thermal habitats, and differentially affect male mating success. Here, we perform a meta-analysis of RNA-seq data which identifies genes encoding odorant binding proteins (in the Obp56h family) as differentially expressed between the heads of males carrying YM and IIIM Differential expression of Obp56h has been associated with variation in male mating behavior in Drosophila melanogaster. We find differences in male mating behavior between house flies carrying the Y chromosomes that are consistent with the relationship between male mating behavior and expression of Obp56h in D. melanogaster. We also find that male mating behaviors in house fly are affected by temperature, and the same temperature differentials further affect the expression of Obp56h genes. However, we show that temperature-dependent effects cannot explain the maintenance of genetic variation for male mating behavior in house fly. Using a network analysis and allele-specific expression measurements, we find evidence that the house fly IIIM chromosome is a trans regulator of Obp56h gene expression. Moreover, we find that Obp56h disproportionately affects the expression of genes on the D. melanogaster chromosome that is homologous to the house fly IIIM chromosome. This provides evidence for a conserved trans regulatory loop involving Obp56h expression that affects male mating behavior in flies. The complex regulatory architecture controlling Obp56h expression suggests that variation in male mating behavior could be maintained by epistasis or pleiotropic constraints.


title: Quantifying the impact of data sharing on variant classification

abstract: 4Healthcare is increasingly leveraging genomic data to inform diagnosis, monitoring, and treatment of certain diseases with genetic predisposition. Associating patient data such as family history and de novo status with a genomic variant helps classify that variant as being pathogenic or benign. Indeed, many variants are already classified by experts, but the majority of variants are very rare, have no associated patient data, and are therefore of uncertain significance. This research models the hypothetical sharing of patient data across institutions in order to accelerate the time it takes to classify a variant. Using conservative assumptions described in the paper, we found that the probability of classifying a pathogenic variant which occurs at the rate of 1 in 100,000 people increases from less than 25% to nearly 80% after just one year when sequencing centers share their clinical data. After 5 years, the probability of classifying such a variant is nearly 100%.

category: genetics


title: Essential role for Ggct in erythrocyte antioxidant defense

abstract: GGCT encodes {gamma}-glutamyl cyclotransferase enzyme activity, and its expression is up-regulated in various human cancers. {gamma}-glutamyl cyclotransferase enzyme activity was originally purified from human red blood cells (RBCs), however physiological function of GGCT in RBCs is still not clear. Here we reported that Ggct deletion in mouse leads to splenomegaly and progressive anemia phenotypes, due to elevated oxidative damage and shortened life span of Ggct-/- RBCs. Ggct-/- RBCs have increased reactive oxygen species (ROS), and are more sensitive to H2O2 induced damage compared to control RBCs. Glutathione (GSH) and GSH synthesis precursor L-cysteine are decreased in Ggct-/- RBCs. Our study suggests a critical function of Ggct in RBC redox balance and life span maintenance through regulating GSH metabolism.


title: BMI-CNV: A Bayesian framework for multiple genotyping platforms detection of copy number variation

abstract:
Whole-exome sequencing (WES) enables detection of Copy number variations (CNVs) with high resolution in protein-coding regions. However, variations in the intergenic or intragenic regions are excluded from studies. Fortunately, samples have been previously sequenced by other genotyping platforms, such as SNP array. Moreover, conventional single sample-based methods suffer from high false discovery rate due to prominent data noise. Therefore, methods for integrating multiple genotyping platforms and samples are highly demanded for improved CNV detection.

We developed BMI-CNV, a Bayesian Multi-sample and Integrative CNV (BMI-CNV) profiling method with data sequenced by both WES and microarray. For the multi-sample integration, we identify the shared CNVs regions across samples using a Bayesian probit stick-breaking process model coupled with a Gaussian Mixture model estimation. With extensive simulations, BMI-CNV outperformed existing methods with remarkably improved accuracy. By applying to the matched 1000 genomes project and HapMap project data, we showed that BMI-CNV accurately detected common variants. We further applied it to The Research of International Cancer of Lung (TRICL) consortium with matched WES and OncoArray data and identified lung cancer risk associated genes in 17q11.2, 1p36.12, 8q23.1 and 5q22.2 regions, which may provide new insights into the etiology of lung cancer.


title: Genomic prediction with allele dosage information in highly polyploid species

abstract: Several studies have shown how to leverage allele dosage information to improve the accuracy of genomic selection models in autotetraploids. In this study we expanded the methodology used for genomic selection in autotetraploids to higher (and mixed) ploidy levels. We adapted the models to build covariance matrices of both additive and digenic dominance effects that are subsequently used in genomic selection models. We applied these models using estimates of ploidy and allele dosage to sugarcane and sweet potato datasets and validated our results by also applying the models in simulated data. For the simulated datasets, including allele dosage information led up to 140% higher mean predictive abilities in comparison to using diploidized markers. Including dominance effects was highly advantageous when using diploidized markers, leading to mean predictive abilities which were up to 115% higher in comparison to only including additive effects. When the frequency of heterozygous genotypes in the population was low, such as in the sugarcane and sweet potato datasets, there was little advantage in including allele dosage information in the models. Overall, we show that including allele dosage can improve genomic selection in highly polyploid species under higher frequency of different heterozygous genotypic classes and high dominance degree levels.


title: Genomic sequencing of Lowe syndrome trios reveal a mechanism for the heterogeneity of neurodevelopmental phenotypes

abstract: Lowe syndrome is an X-linked recessive monogenic disorder resulting from mutations in the OCRL gene that encodes a phosphatidylinositol 4,5 bisphosphate 5-phosphatase. The disease affects three organs-the kidney, brain and eye and clinically manifests as proximal renal tubule dysfunction, neurodevelopmental delay and congenital cataract. Although Lowe syndrome is a monogenic disorder, there is considerable heterogeneity in clinical presentation; some individuals show primarily renal symptoms with minimal neurodevelopmental impact whereas others show neurodevelopmental defect with minimal renal symptoms. However, the molecular and cellular mechanisms underlying this clinical heterogeneity remain unknown. Here we analyze a Lowe syndrome family in whom affected members show clinical heterogeneity with respect to the neurodevelopmental phenotype despite carrying an identical mutation in the OCRL gene. Genome sequencing and variant analysis in this family identified a large number of damaging variants in each patient. Using novel analytical pipelines and segregation analysis we prioritize variants uniquely present in the patient with the severe neurodevelopmental phenotype compared to those with milder clinical features. The identity of genes carrying such variants underscore the role of additional gene products enriched in the brain or highly expressed during brain development that may be determinants of the neurodevelopmental phenotype in Lowe syndrome. We also identify a heterozygous variant in CEP290, previously implicated in ciliopathies that underscores the potential role of OCRL in regulating ciliary function that may impact brain development. More generally, our findings demonstrate analytic approaches to identify high-confidence genetic variants that could underpin the phenotypic heterogeneity observed in monogenic disorders.

category: genomics


title: Genomics analysis of Drosophila sechellia response to Morinda citrifolia fruit diet

abstract: Drosophila sechellia is an island endemic host specialist that has evolved to consume the toxic fruit of Morinda citrifolia, also known as noni fruit. Recent studies by our group and others have examined genome-wide gene expression responses of fruit flies to individual highly abundant compounds found in noni responsible for the fruits unique chemistry and toxicity. In order to relate these reductionist experiments to the gene expression responses to feeding on noni fruit itself, we fed rotten noni fruit to adult female D. sechellia and performed RNA-sequencing. Combining the reductionist and more wholistic approaches, we have identified candidate genes that may contribute to each individual compound and those that play a more general role in response to the fruit as a whole. Using the compound specific and general responses, we used transcription factor prediction analyses to identify the regulatory networks and specific regulators involved in the responses to each compound and the fruit itself. The identified genes and regulators represent the possible genetic mechanisms and biochemical pathways that contribute to toxin resistance and noni specialization in D. sechellia.


title: Mechanisms of DNA repair have evolved to minimise the probability of nonsense mutations.

abstract: Variation in sequence mutability has important implications for evolutionary models and predicting disease occurrence, and is driven in part by evolutionary divergence in mechanisms of DNA repair. The aim of this study was twofold: first, to assess the effect of local sequence context on substitution rates in the mouse lineage; second, to investigate the relationship between sequence mutability and selection. We show that the 7-mer context (i.e three bases either side of the base of interest) explains more variation in substitution rates between chromosomes in the mouse lineage than either the 3-mer, 5-mer, or 9-mer contexts. Furthermore, we also show that 7-mer substitutions with the potential to cause nonsense mutations when they occur in translated sequences occur at a lower rate across the genome than 7-mer substitutions with the potential to cause synonymous mutations. We propose that mechanisms of DNA repair have evolved to prioritise substitutions that are more likely to be deleterious to fitness.


title: Reproductive skew, fitness costs, and winner-loser effects in social-dominance evolution

abstract: Social hierarchies can increase reproductive skew in group-living animals. Using game theory we investigate how the opportunity for differently ranked individuals to acquire resources influences reproductive skew, costs of hierarchy formation, and winner and loser effects. Individuals adjust their aggressive and submissive behaviour through reinforcement learning. The learning is based on perceived rewards and penalties, which depend on relative fighting ability. From individualbased simulations we determine evolutionary equilibria of traits that control an individuals learning. We examine situations that differ in the extent of monopolisation of contested resources by dominants and in the amounts of uncontested resources that are distributed independently of rank. With costly fighting, we find that stable dominance hierarchies form, such that reproductive skew mirrors the distribution of resources over ranks. Individuals pay substantial costs of interacting, in particular in high-skew situations, with the highest costs paid by intermediately ranked individuals. For cases where dominants monopolise contested resources there are notable winner and loser effects, with winner effects for high ranks and very pronounced loser effects for lower ranks. The effects are instead weak when acquired resources increase linearly with rank. We compare our results on contest costs and winner-loser effects with field and experimental observations.

category: animal behavior and cognition


title: Online Reference Trajectory Adaptation: A Personalized Control Strategy for Lower Limb Exoskeletons

abstract: This paper presents a novel method for reference trajectory adaptation in lower limb rehabilitation exoskeletons during walking. Our adaptation rule is extracted from a cost function that penalizes both interaction force and trajectory modification. By adding trajectory modification term into the cost function, we restrict the boundaries of the reference trajectory adaptation according to the patients motor capacity. The performance of the proposed adaptation method is studied analytically in terms of convergence and optimality. We also developed a realistic dynamic walking simulator and utilized it in performance analysis of the presented method. The proposed trajectory adaptation technique guarantees convergence to a stable, reliable, and rhythmic reference trajectory with no prior knowledge about the human intended motion. Our simulations demonstrate the convergence of exoskeleton trajectories to those of simulated healthy subjects while the exoskeleton trajectories adapt less to the trajectories of patients with reduced motor capacity (less reliable trajectories). Furthermore, the gait stability and spatiotemporal parameters such as step time symmetry and minimum toe off clearance enhanced by the adaptation in all subjects. The presented mathematical analysis and simulation results show the applicability and effectiveness of the proposed method and its potential to be applied for trajectory adaptation in lower limb rehabilitation exoskeletons.

category: bioengineering


title: Evaluating the cost of pharmaceutical purification for a long-duration space exploration medical foundry

abstract: There are medical treatment vulnerabilities in longer-duration space missions present in the current International Space Station crew health care system with risks, arising from spaceflight-accelerated pharmaceutical degradation and resupply lag times. Bioregenerative life support systems may be a way to close this risk gap by leveraging in situ resource utilization (ISRU) to perform pharmaceutical synthesis and purification. Recent literature has begun to consider biological ISRU using microbes and plants as the basis for pharmaceutical life support technologies. However, there has not yet been a rigorous analysis of the processing and quality systems required to implement biologically-produced pharmaceuticals for human medical treatment. In this work, we use the equivalent system mass (ESM) metric to evaluate pharmaceutical purification processing strategies for longer-duration space exploration missions. Monoclonal antibodies, representing a diverse therapeutic platform capable of treating multiple space-relevant disease states, were selected as the target products for this analysis. We investigate the ESM resource costs (mass, volume, power, cooling, and crew time) of an affinity-based capture step for monoclonal antibody purification as a test case within a manned Mars mission architecture. We compare six technologies (three biotic capture methods and three abiotic capture methods), optimize scheduling to minimize ESM for each technology, and perform scenario analysis to consider a range of input stream compositions and pharmaceutical demand. We also compare the base case ESM to scenarios of alternative mission configuration, equipment models, and technology reusability. Throughout the analyses, we identify key areas for development of pharmaceutical life support technology and improvement of the ESM framework for assessment of bioregenerative life support technologies.


title: TedSim: temporal dynamics simulation of single cell RNA-sequencing data and cell division history

abstract:
Recently, the combined scRNA-seq and CRISPR/Cas9 genome editing technologies have enabled simultaneous readouts of gene expressions and lineage barcodes, which allows for the reconstruction of the cell division tree, and makes it possible to trace the origin of each cell type. Computational methods are emerging to take advantage of the jointly profiled scRNA-seq and lineage barcode data to better reconstruct the cell division history or to infer the cell state trajectories.

Here, we present TedSim (single cell Temporal dynamics Simulator), a simulator that simulates the cell division events from the root cell to present-day cells, simultaneously generating the lineage barcodes and scRNA-seq data. In particular, TedSim generates cells from multiple cell types through cell division events. TedSim can be used to benchmark and investigate computational methods which use either or both of the two types of data, scRNA-seq and lineage barcodes, to study cell lineages, ancestral cells or cell trajectories. TedSim is available at: https://github.com/Galaxeee/TedSim.

category: bioinformatics


title: A combined test for feature selection on sparse metaproteomics data - alternative to missing value imputation

abstract: One of the difficulties encountered in the statistical analysis of metaproteomics data is the high proportion of missing values, which are usually treated by imputation. Nevertheless, imputation methods are based on restrictive assumptions regarding missingness mechanisms, namely "at random" or "not at random". To circumvent these limitations in the context of feature selection in a multi-class comparison, we propose a univariate selection method that combines a test of association between missingness and classes, and a test for difference of observed intensities between classes. This approach implicitly handles both missingness mechanisms. We performed a quantitative and qualitative comparison of our procedure with imputation-based feature selection methods on two experimental data sets, as well as simulated data with various scenarios regarding the missingness mechanisms and the nature of the difference of expression (differential intensity or differential missingness). Whereas we observed similar performances in terms of prediction on the experimental data set, the feature ranking and selection from various imputation-based methods were strongly divergent. We showed that the combined test reaches a compromise by correlating reasonably with other methods, and remains efficient in all simulated scenarios unlike imputation-based feature selection methods.


title: A combined test for feature selection on sparse metaproteomics data - an alternative to missing value imputation


title: ROPIUS0: A deep learning-based protocol for protein structure prediction and model selection and its performance in CASP14

abstract: Protein structure prediction has recently been revolutionized when AlphaFold2 [1] predicted protein structures with near-experimental accuracy in the latest CASP14 season of critical assessment of methods of protein structure prediction (CASP). Among numerous implications, this breakthrough has led to a rapidly growing number of high-quality structural models [2]. We present a protocol ROPIUS0 for protein structure prediction and model selection and discuss its benefits in the new era of structure prediction. At the core of the ROPIUS0 protocol is the deep learning module developed for the selection of protein structural models. It is shown that the direct use of predicted inter-residue distances may be sufficient to discriminate between correct and incorrect protein folds, considering only a small fraction of predicted distances. We extensively tested the protocol: In the latest CASP14 prediction season, a ROPIUS0 variant based on model selection ranked 13th in the category of tertiary structure prediction. Its performance is on par with top-performing automated prediction servers when tested on the CASP13 dataset, and it performs similarly on a CAMEO dataset. The results suggest ways to improve searching for structurally similar and homologous proteins without considerably increasing speed. Our new open-source threading tool based on comparing a subset of inter-residue distances demonstrates the effectiveness and application of the deep learning module of the ROPIUS0 protocol.


title: Efficient Generative Modelling of Protein Structure Fragments using a Deep Markov Model

abstract: Fragment libraries are often used in protein structure prediction, simulation and design as a means to significantly reduce the vast conformational search space. Current state-of-the-art methods for fragment library generation do not properly account for aleatory and epistemic uncertainty, respectively due to the dynamic nature of proteins and experimental errors in protein structures. Additionally, they typically rely on information that is not generally or readily available, such as homologous sequences, related protein structures and other complementary information. To address these issues, we developed BIFROST, a novel take on the fragment library problem based on a Deep Markov Model architecture combined with directional statistics for angular degrees of freedom, implemented in the deep probabilistic programming language Pyro. BIFROST is a probabilistic, generative model of the protein backbone dihedral angles conditioned solely on the amino acid sequence. BIFROST generates fragment libraries with a quality on par with current state-of-the-art methods at a fraction of the run-time, while requiring considerably less information and allowing efficient evaluation of probabilities.


title: Contamination detection and microbiome exploration with GRIMER

abstract: Exploring microbiome data is a time-consuming task that can be only partially automated due to the specific requirements and goals of each project. Visualizations and analysis platforms are crucial to better guide this step. Best practices in the field are constantly evolving and many pitfalls can lead to biased outcomes. Compositionality of data and sample contamination are two important points that should be carefully considered in early stages of microbiome studies. Detecting contamination can be a challenging task, especially in low-biomass samples or in studies lacking proper controls by design. However, external evidences and commonly identified contaminant taxa can be used to discover and mitigate contamination. We propose GRIMER, a tool that automates analysis, generates plots and runs external tools to create a portable dashboard integrating annotation, taxonomy and metadata. It unifies several sources of evidence towards contamination detection. GRIMER is independent of quantification methods and directly analyses contingency tables to create an interactive and offline report. GRIMER reports can be created in seconds and are accessible for non-specialists, providing an intuitive set of charts to explore data distribution among observations and samples and its connections with external sources. Further, we compiled an extensive list of common contaminants and possible external contaminant taxa reported in the literature and use it to annotate data. GRIMER is open-source and available at: https://gitlab.com/dacs-hpi/grimer.


title: Deep Learning Decodes Multi-layered Roles of Polycomb Repressor Complex 2 During Mouse CNS Development

abstract: A prominent aspect of most, if not all, central nervous systems (CNSs) is that anterior regions (brain) are larger than posterior ones (spinal cord). Studies in Drosophila and mouse have revealed that the Polycomb Repressor Complex 2 (PRC2), a protein complex responsible for applying key repressive histone modifications, acts by several mechanisms to promote anterior CNS expansion. However, it is unclear what the full spectrum of PRC2 action is during embryonic CNS development and how PRC2 integrates with the epigenetic landscape. We removed PRC2 function from the developing mouse CNS, by mutating the key gene Eed, and generated spatio-temporal transcriptomic data. To decode the role of PRC2, we developed a method that incorporates standard statistical analyses with probabilistic deep learning to integrate the transcriptomic response to PRC2 inactivation with epigenetic information from ENCODE. This multi-variate analysis corroborates the central involvement of PRC2 in anterior CNS expansion, and reveals layered regulation via PRC2. These findings uncover a differential logic for the role of PRC2 upon functionally distinct gene categories that drive CNS anterior expansion. To support the analysis of emerging multi-modal datasets, we provide a novel bioinformatics package that integrates transcriptomic and epigenetic datasets to identify regulatory underpinnings of heterogeneous biological processes.


title: Integrated gene landscapes uncover multi-layered roles of repressive histone marks during mouse CNS development


title: Annotating precision for integrative structural models using deep learning

abstract:
SummaryA single precision value is currently reported for an integrative model. However, precision may vary for different regions of an integrative model owing to varying amounts of input information. We develop PrISM (Precision for Integrative Structural Models), to efficiently identify high and low-precision regions for integrative models.

Availability and ImplementationPrISM is available at https://github.com/isblab/prism; benchmark data used in this paper is available at doi:10.5281/zenodo.6241200

Contactshruthiv@ncbs.res.in

Supplementary informationThis article contains Supplementary data online.


title: PrISM: precision for integrative structural models


title: TKGWV2: An ancient DNA relatedness pipeline for ultra-low coverage whole genome shotgun data

abstract: Estimation of genetically related individuals is playing an increasingly important role in the ancient DNA field. In recent years, the numbers of sequenced individuals from single sites have been increasing, reflecting a growing interest in understanding the familial and social organisation of ancient populations. Although a few different methods have been specifically developed for ancient DNA, namely to tackle issues such as low-coverage homozygous data, they require a 0.1 - 1x minimum average genomic coverage per analysed pair of individuals between. Here we present an updated version of a method that enables estimates of 1st and 2nd-degrees of relatedness with as little as 0.026x average coverage, or around 1.3 million aligned reads per sample - 4 times less data than 0.1x. By using simulated data to estimate false positive error rates, we further show that a threshold even as low as 0.012x, or around 600,000 reads, will always show 1st-degree relationships as related. Lastly, by applying this method to published data, we are able to identify previously undocumented relationships using individuals previously excluded from kinship analysis due to their very low coverage. This methodological improvement has the potential to enable relatedness estimation on ancient whole genome shotgun data during routine low-coverage screening, and therefore improve project management when decisions need to be made on which individuals are to be further sequenced.


title: Role of mitochondrial genetic interactions in determining adaptation to high altitude in human population around the globe

abstract: Physiological and haplogroup studies performed to understand high-altitude adaptation in humans are limited to individual genes and polymorphic sites. Due to stochastic evolutionary forces, the frequency of a polymorphism is affected by changes in the frequency of a near-by polymorphism on the same DNA sample making them connected in terms of evolution. Here, first we provide a method to model these mitochondrial polymorphisms as "co-mutation networks" for three high-altitude populations, Tibetan, Ethiopian and Andean. Then, by transforming these co-mutation networks into weighted and undirected gene-gene interaction (GGI) networks, we were able to identify functionally enriched genetic interactions of CYB and CO3 genes in Tibetan and Andean populations, while NADH dehydrogenase genes in Ethiopian population playing a significant role in high altitude adaptation. These co-mutation based genetic networks provide insights into the role of different set of genes in high-altitude adaptation human sub-populations.


title: Stochastic Resonance Behavior of DNA Translocation with an Oscillatory Electric Field

abstract: Stochastic resonance (SR) describes the synchronization between noise of a system and an applied oscillating field to achieve an optimized response signal. In this work, we use simulations to investigate the phenomenon of SR of a single stranded DNA driven through a nanopore when an oscillating electric field is added. The system is comprised of a MspA protein nanopore embedded in a membrane and different lengths of DNA is driven from one end of the pore to the other via a constant potential difference. We superimposed an oscillating electric field on top of the existing electric field. The source of noise is due to thermal fluctuations, since the system is immersed in solution at room temperature. Here, the signal optimization we seek is the increase in translocation time of DNA through the protein nanopore. Normally, translocation time scales linearly with DNA length and inversely with driving force in a drift dominated regime. We found a non-monotonic dependence of the mean translocation time with the frequency of the oscillating field. This non-monotonic behavior of the translocation time is observed for all lengths of DNA, but SR occurs only for longer DNA. Furthermore, we also see evidence of DNA extension being influenced by the oscillating field while moving through the nanopore.

category: biophysics


title: Kinetic pathways of water exchange in the first hydration shell of magnesium: Influence of water model and ionic force field

abstract: Water exchange between the first and second hydration shell is essential for the role of Mg2+ in biochemical processes. In order to provide microscopic insights into the exchange mechanism, we resolve the exchange pathways by all-atom molecular dynamics simulations and transition path sampling. Since the exchange kinetics relies on the choice of the water model and the ionic force field, we systematically investigate the influence of seven different polarizable and non-polarizable water and three different Mg2+ models. In all cases, water exchange can occur either via an indirect or direct mechanism (exchanging molecules occupy different/same position on water octahedron). In addition, the results reveal a crossover from an interchange dissociative (Id) to an associative (Ia) reaction mechanism dependent on the range of the Mg2+-water interaction potential of the respective force field. Standard non-polarizable force fields follow the Id mechanism in agreement with experimental results. By contrast, polarizable and long-ranged non-polarizable force fields follow the Ia mechanism. Our results provide a comprehensive view on the influence of the water model and ionic force field on the exchange dynamics and the foundation to assess the choice of the force field in biomolecular simulations.


title: Coupled dynamics of charged macromolecules and counterions mediated by binding enzymes

abstract: We investigate the role of active coupling on the transport properties of the macromolecules. The active coupling comes due to bound enzymes with a segment of the macromolecule wherein the enzyme exerts an electrostatic force on the segment of the macromolecule, and eventually, it gets unbound due to the thermal fluctuations. This binding and unbinding process generates active fluctuations in the dynamics of the macromolecule. Starting with segment dynamics and correlations for three dynamical models with active coupling, we obtain the cooperative diffusivity for the realistic charged macromolecules with hydrodynamics. First, we construct the three models by incorporating the features of a real polymer systematically, starting from simple Rouse dynamics with active coupling. We further include segment-segment interactions and in addition, hydrodynamic interactions with active coupling. Our obtained scaling form for segment-segment correlations for the models in terms of the size exponent of the polymer indicating that hydrodynamic and segment-segment interactions along with the active coupling lead to new scaling regimes. We finally study the dynamics of a homogeneously charged flexible polymer in an infinitely dilute solution where enzymes and counterions affect the dynamics of the polymers. We analytically investigate how these active fluctuations affect the coupled dynamics of the polymer and counterions. It turns out that these active fluctuations enhance the effective diffusivity of the polymer. The derived closed-form expression for diffusivity is pertinent to accurate interpretation of light scattering data on multi-component systems with binding-unbinding equilibria.


title: Nucleotide Dynamics During Flossing of Polycation-DNA-Polycation through a Nanopore using Molecular Dynamics

abstract: The full potential of solid state nanopores is yet to be realized for genome sequencing. Due to its robustness it can handle strong voltage amplitude and frequency. The effect of strong alternating voltage on the dynamics of nucleotides during translocation has been explored. We proposed a setup consisting of single stranded DNA covalently linked with symmetric polycations at both ends fashioned after the proposal of Kasianowicz. 1 Such a setup allows for repeated back and forth motion of the DNA along the nanopore (1.45 nm diameter and 1.53 nm thick) by simply switching the voltage polarity if the polycation tail is sufficiently long ([≥] 10) and the applied voltage is below 0.72 volts, but the average residence time of the nucleotides are too small to be of any practical use (6-30 ns). When alternating voltage of higher frequency is applied, it enhances the average residence time of the nucleotides by an order of magnitude to [~] 0.1 {micro}s relative to direct voltage but the individual trajectories are too stochastic. Since, we are able to collect repeated read on the dynamics of individual nucleotides, we obtained the most probable time of appearance of a nucleotide within the nanopore. With such construct we were able to get almost linear dependence of most probable time versus nucleotide index, after gaussian fitting.


title: Unfolding of RNA via Translocation Through a Nanopore

abstract: RNA unfolding and refolding are important biological phenomena, which occur during the transfer of genetic information from DNA to RNA to proteins. During these processes, RNA is found in single stranded, secondary and tertiary structures, including secondary conformations like hairpins and pseudoknots. Understanding the diverse conformations of RNA and how these influence the dynamics of unfolding and refolding is crucial to gain insight to fundamental biological processes. In this work, we employ coarse-grained Langevin dynamics simulations of a simple model of different RNA hairpins passing through a geometric nanopore to find the influence of structural changes on the translocation dynamics. The threshold voltage of unfolding depends on the length of the hairpin attached to the tail. The lag time to unfold is longer for smaller applied voltages and for the architectures containing a longer hairpin attached to the tail. Chain translocation dynamics for different architectures are largely collapsed by the threshold. A distinct signature for the base unfolding time was observed for the bases around the unpaired bases in all the RNA hairpin models. These results can motivate future technologies or experiments that use translocation to predict secondary structures of polynucleotides.


title: Langevin dynamics simulation of protein dynamics in nanopores at microsecond timescales

abstract: With rapid advancement in the fields of nanopore analysis of protein, it has become imperative to develop modeling framework for understanding the protein dynamics in nanopores. Such modeling framework should include the effects of electro-osmosis, as it plays significant role during protein translocation in confinement. Currently, the molecular dynamics simulations that include the hydrodynamic effects are limited to a timescale of few 100 ns. These simulations give insight about important events like protein unfolding which occurs in this timescale. But many electrophoresis experiments are limited by a detector resolution of ~ 2.5 s. Analytical theory has been used to interpret protein dynamics at such large timescale. There is a need for molecular modeling of more complex environment and protein shapes which cannot be accounted for by analytical theory. We have developed a framework to study globular protein dynamics in nanopores by using langevin dynamics on a rigid body model of protein and the hydrodynamics is accounted by analytical theory for simple cylindrical nanopore geometry. This framework has been applied to study the dynamics of Ubiquitin translocation in SiNx nanopore by Nir et al26. They have reported 7 times decrease in average dwell time of the protein inside the nanopore in response to a small change in pH from 7.0 to 7.2 and the modification of protein charge was attributed for such drastic change. Closer examination using our simulation revealed that the electro-osmotic effects originating due to very small change in the surface electrostatic potential of the nanopore could lead to such a drastic change in protein dynamics.


title: Voltage-driven polyelectrolyte complexation inside a nanopore

abstract: We have investigated how a pair of oppositely charged macromolecules can be driven by an electric field to form a polyelectrolyte complex inside a nanopore. To observe and isolate an individual complex pair, a model protein nanopore, embedded in artificial phospholipid membrane, allowing compartmentalization (cis/trans) is employed. A polyanion in the cis and a polycation in the trans compartments are subjected to electrophoretic capture by the pore. We find that the measured ionic current across the pore has a distinguishable signature of complex formation, which is different from the signature of the passage of individual molecules through the pore. The ionic current signature allows us to detect the interaction between the two oppositely charged macromolecules and thus, enables us to measure the lifetime of the complex inside the nanopore. After showing that we can isolate a complex pair in the nanopore, we studied the effects of molecular identity on the nature of interaction in different complex pairs. In contrast to the irreversible conductance state of the alpha-hemolysin (HL) channel in the complexation of poly-styrene-sulfonate (PSS) and poly-L-lysine (PLL), a reversible conductance state is observed during complexation between single stranded DNA (ssDNA) and PLL. This suggests that there is a weak interaction between ssDNA and PLL, when compared to the interaction in a PSS-PLL complex. Analysis of the PSS-PLL complexation events and its lifetime inside the nanopore supports a four step-mechanism: (i) The polyanion is captured by the pore, (ii) the polyanion starts threading through the pore. (iii) The polycation is captured, a complex pair is formed in the pore, and the polyanion slides along the polycation. (iv) The complex pair can be pulled through the pore into the trans compartment or it can dissociate. Additionally, we have developed a simple theoretical model, which describes the lifetime of the complex inside the pore. The observed reversible two-state conductance across HL channel during ssDNA-PLL complexation, is described as the binding/unbinding of PLL during the translocation of ssDNA. This enables us to evaluate the apparent rate constants for association/dissociation and equilibrium dissociation constants for the interaction of PLL with ssDNA. This work throws light on the behavior of polyelectrolyte complexes in an electric field and enhances our understanding of the electrical aspects of inter-macromolecular interactions, which plays an extremely important role in the organization of macromolecules in the crowded and confined cellular environment.


title: Inhibiting adenine synthesis attenuates glioblastoma cell stemness and temozolomide resistance

abstract: Glioblastoma (GBM) is a lethal brain cancer exhibiting high levels of drug resistance, a feature partially imparted by tumor cell stemness. Recent work shows that homozygous MTAP deletion, a genetic alteration occurring in about half of all GBMs, promotes stemness in GBM cells. Exploiting MTAP loss-conferred deficiency in adenine salvage, we demonstrate that transient adenine blockade via treatment with L-Alanosine (ALA), an inhibitor of de novo adenine synthesis, attenuates stemness of MTAP-deficient GBM cells. This ALA-induced reduction in stemness is accompanied by compromised mitochondrial function, highlighted by diminished spare respiratory capacity. Direct pharmacological inhibition of mitochondrial respiration recapitulates the effect of ALA on GBM cell stemness, suggesting ALA targets stemness partially via affecting mitochondrial function. Finally, in agreement with diminished stemness and compromised mitochondrial function, we show that ALA sensitizes GBM cells to temozolomide (TMZ) in vitro and in an orthotopic GBM model. Collectively, these results identify critical roles of adenine supply in maintaining mitochondrial function and stemness of GBM cells, highlight a critical role of mitochondrial function in sustaining GBM stemness, and implicate adenine synthesis inhibition as a complementary approach for treating MTAP-deleted GBMs.

category: cancer biology


title: Immunosuppressive traits of the hybrid epithelial/mesenchymal phenotype

abstract: Recent preclinical and clinical data suggests enhanced metastatic fitness of hybrid epithelial/ mesenchymal (E/M) phenotypes, but mechanistic details regarding their survival strategies during metastasis remain unclear. Here, we investigate immune-evasive strategies of hybrid E/M states. We construct and simulate the dynamics of a minimalistic regulatory network encompassing the known associations among regulators of EMT (epithelial-mesenchymal transition) and PD-L1, an established immune-suppressor. Our model simulations, integrated with single-cell and bulk RNA-seq data analysis, elucidate that hybrid E/M cells can have high levels of PD-L1, similar to those seen in cells with a full EMT phenotype, thus obviating the need for cancer cells to undergo a full EMT to be immune-evasive. Specifically, in breast cancer, we show the co-existence of hybrid E/M phenotypes, enhanced resistance to anti-estrogen therapy and increased PD-L1 levels. Our results underscore how the emergent dynamics of interconnected regulatory networks can coordinate different axes of cellular fitness during metastasis.


title: WEE1 kinase is a therapeutic vulnerability in CIC-DUX4 undifferentiated sarcoma.

abstract: CIC-DUX4 rearrangements define an aggressive and chemotherapy-insensitive subset of undifferentiated sarcomas. The CIC-DUX4 fusion drives oncogenesis through direct transcriptional upregulation of cell cycle and DNA replication genes. Notably, CIC-DUX4- mediated CCNE1 upregulation compromises the G1/S transition, conferring a potential survival dependence on the G2/M cell cycle checkpoint. Through an integrative transcriptional and kinase activity screen using patient-derived specimens, we now show that CIC-DUX4 sarcomas depend on the G2/M checkpoint regulator, WEE1, as an adaptive survival mechanism. Specifically, CIC-DUX4 sarcomas depend on WEE1 activity to limit DNA damage and unscheduled mitotic entry. Consequently, genetic or pharmacologic WEE1 inhibition in vitro and in vivo leads to rapid DNA damage-associated apoptotic induction of patient-derived CIC-DUX4 sarcomas. Thus, we identify WEE1 as an actionable therapeutic vulnerability in CIC-DUX4 sarcomas.


title: Anti-clotting functions of lymphatics form the natural on-off switch for immune recognition by controlling the antigens and immune cells access to the lymph nodes

abstract: The ability of lymph to clot indicates that similar to blood vessels lymphatics must have means to counteract this process. We analyzed their hemostatic properties, tailoring them for the potential therapeutic applications. Inflammatory stimuli induced tissue factor-dependent focal lymph clotting while blocking thrombomodulin lead to widespread but also transient occlusion of collecting vessels. Decellularization of lymphatics resulted in tissue factor-independent lymphatic occlusion by widespread and persistent lymph clots. These decellularized basement membrane remnants of collectors were capable of macromolecules drainage and leukocytes transit only when lymph clotting was inhibited with heparin. In occluded  ghost vessels, fibrin was replaced with transient basement membrane-rich inclusions to be eventually reperfused. During regeneration, ghost vessels were filled with granuloma-like clusters of antigen-presenting and T cells. Despite that, immune response against allografts placed under non-drained skin was not developed as long lymphatics remained occluded, and graft survival was prolonged together with the delay of lymphatic regeneration with anti-lymphangiogenic therapy. Other potential applications of lymphatic hemostatic control functions include blocking pathogen or metastasis spread but also should be considered in the treatment of lymphedema.

category: physiology


title: Effect of a low water salinity on oocyte maturation, oocyte hydration, ovulation and egg quality in rainbow trout.

abstract:
Water salinity is an important environmental factor known to have detrimental effects on salmonid reproduction, mostly when migrating female broodfish are held in sea water. In contrast, data obtained in freshwater are scarce and the impact of low water salinity, and more specifically of low water concentrations in sodium, chloride and potassium, during reproduction in freshwater is currently unknown. For this reason, and because ion and water fluxes are critical for the final steps of the female gamete formation, including oocyte hydration and ovulation, the aim of the present study was to investigate the impact of low salinity water on final oocyte maturation, ovulation and, ultimately, on egg quality, using rainbow trout as a physiological model and relevant aquaculture species.

Fish from the same commercial strain were raised either in a site characterized by low concentrations of Na+, K+, and Cl- ions in the water or in a closely located control site exhibiting higher concentration in these elements. Egg quality and duration of final oocyte maturation were investigated using innovative phenotyping tools such as automatic assessment of egg viability using the VisEgg system and non-invasive echograph-based monitoring of final oocyte maturation duration, respectively. Oocyte hydration during final oocyte maturation and after ovulation was also investigated. Finally, molecular phenotyping was performed using real-time PCR-based monitoring of several key players of final oocyte maturation and ovulation associated with ion and water transport, inflammation, proteolytic activity, and coagulation. Oocyte hydration and gene expression data were analyzed in the light of the duration of final oocyte maturation.

Here we show that low water salinity (i.e., low water concentration in chloride, sodium and potassium) negatively influences final oocyte maturation, ovulation and, ultimately, egg quality. Low water salinity triggered delayed ovulation and lower oocyte viability. When investigating the impact of low water salinity on final oocyte maturation duration, individuals presenting the most severe phenotypes exhibited impaired oocyte hydration and abnormally reduced gene expression levels of several key players of the ovulatory process. While the under expression of water (i.e., aquaporins) and ion (i.e., solute carriers) transporters is consistent with impaired oocyte hydration, our observations also indicate that the overall ovulatory gene expression program is disrupted. Our results raise the question of the mechanisms underlying the negative influence of low salinity water on the dynamics of the preovulatory phase, on the control of the oocyte homeostasis, including hydration, and on the overall success of the maturation-ovulation process.

HighlightsO_LILow water salinity impairs final oocyte maturation and egg quality in rainbow trout
C_LIO_LILow water salinity induces delayed ovulation and impaired oocyte hydration
C_LIO_LILow water salinity induces a dysregulation of several key ovulatory genes
C_LIO_LIMonitoring of final oocyte maturation can be performed using ultrasound staging
C_LI


title: Effect of a low concentration in chloride, sodium and potassium on oocyte maturation, oocyte hydration, ovulation and egg quality in rainbow trout.


title: The RNA-binding protein RBP45D of Arabidopis plays a role in epigenetic control of flowering time and DCL3-independent RNA-directed DNA methylation

abstract: RNA-directed DNA methylation (RdDM) helps to defend plants against invasive nucleic acids. In the canonical form of RdDM, 24-nt small interfering RNAs (siRNAs) are produced by DICER-LIKE 3 (DCL3). Here, we describe the Arabidopsis thaliana prors1 (LUC) transgenic system, in which transcriptional gene silencing (TGS) is independent of DLC3. A forward genetics screen performed with this system identified both known components of RdDM, and the RNA-binding protein RBP45D. RBP45D promotes DNA methylation, and its loss delays flowering, especially at high temperature, presumably mediated by elevated FLC levels. RBP45D is localized to the nucleus, where it is associated with snRNAs and snoRNAs. RBP45D maintains siRNA production originating from the LUC transgene, but does not alter mRNA levels or affect processing of transcripts of known RdDM genes. We suggest that RBPD45 facilitates DCL3-independent siRNA production by stabilising either the precursor RNA or the - as yet unidentified - slicer protein.

category: plant biology


title: Multiplexed Screening of Population-level Synthetic Gene Oscillator Libraries

abstract: Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales.

category: synthetic biology


title: Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics


title: FACT-mediated maintenance of chromatin integrity during transcription is essential for viability of mammalian stem cells

abstract: Preservation of nucleosomes during replication has been extensively studied, while the maintenance of nucleosomes during transcription has gotten less attention. The histone chaperone FACT is involved in transcription elongation, although whether it disassembles or assembles nucleosomes during this process is still unclear. We deleted the FACT subunit in adult mice to clarify the function of FACT in mammals. FACT loss was lethal due to the loss of the earliest progenitors in bone marrow and intestine, while mor differentiated cells were not affected. Using cells isolated from several tissues, we showed that FACT loss was lethal only for stem cells but not cells differentiated in vitro. FACT depletion led to increased chromatin accessibility in a transcription-dependent manner, suggesting that nucleosomes are lost during transcription in the absence of FACT. The most prominent response to the loss of nucleosomes was the activation of interferon signaling and the accumulation of immunocytes in sensitive organs. FACT maintained chromatin integrity during transcription in mammalian adult stem cells, suggesting that chromatin transcription in these cells is different from more differentiated cells.

category: systems biology


title: Oleate metabolism using kinetic 13C dilution strategy deciphered the potential role of global transcription regulator arcA in Escherichia coli

abstract: Microbial metabolism of long-chain fatty acids (LCFA; > C12) is of relevance owing to their presence in various nutrient niches. Microbes have evolved to metabolize LCFA by expressing relevant genes coordinated by various transcriptional regulators. Among the global transcriptional regulators, the metabolic control conferred by arcA (aerobic respiration control) under a LCFA medium is lacking. This work is targeted to unravel the metabolic features of E.coli MG1655 and its knockout strain {Delta}arcA under oleate (C18:1) as a sole carbon source, providing novel insights into the flexibility of the global regulators in maintaining the cellular physiology. Owing to the availability and cost of stable isotope LCFA tracers, we adopted a novel kinetic 13C dilution strategy. This allowed us to quantify the 13C dilution rates in the amino acids that retro-biosynthetically shed light on the central metabolic pathways in actively growing cells. Our data comprehensively mapped oleate oxidization in E.coli via the pathways of {beta}-oxidation, TCA cycle, anaplerotic and gluconeogenesis. Interestingly, arcA knockout showed expeditious growth (~60%) along with an increased oleate utilization rate (~55%) relative to the wild-type. {Delta}arcA also exhibited higher 13C dilution rates (> 20%) in proteinogenic amino acids than the wild-type. Overall, the study established the de-repression effect conferred by {Delta}arcA in E.coli, which resulted in a phenotype with reprogrammed metabolism favouring higher oleate assimilation. The outcomes suggest rational metabolic engineering of regulators as a strategy to develop smart cells for enhanced biotransformation of LCFA. This study also opens an avenue for adopting a kinetic 13C dilution strategy to decipher the cellular metabolism of a plethora of substrates, including other LCFA in microbes.


title: Biorational Management And Mycosis Studies Of Grape Thrips

abstract:
Field experiment was conducted to test the bioefficacy of various biorational insecticides against grapevine thrips at AICRP on Fruits, Department of Horticulture, MPKV, Rahuri, during 2017-18. Results revealed that standard check emamectin benzoate 5% SG consistently proved to be the most promising by recording the least thrips population (3.10/shoot). Among biorational insecticides neem oil 2% (4.09/shoot) proved as best treatment followed by, karanj oil (4.51/shoot) and neemazol (5.08/shoot). While entomopathogenic fungi Lecanicillium lecanii recorded less population (4.24/shoot) and emerged as best treatment as compared to the Metarhizium anisopliae (4.87/shoot) and Beauveria bassiana (5.34/shoot). However chlli methanolic extract (6.29/shoot), garlic methanolic extract (6.78/shoot), chilli water extract (6.85/shoot) and garlic water extract (7.08/shoot) are the least effective treatments.

Incremental Cost Benefit Ratio (ICBR) in respect of different treatments ranged between 1.30 to 7.92. The highest ICBR of 1:7.92 was recorded in the treatment with emamectin benzoate 5 SG, and it was followed by Lecanicillium lecanii (1:6.34) and Metarhizium anisopliae (1:5.32). Although neem oil and karanj oil has great reduction of thrips population, but has less incremental cost benefit ratio i.e. 2.81 and 3.04, respectively, due to high dose and its cost.

The mycosis test of three entomopathogenic fungi viz. Beauveria bassiana, Metarhizium anisopliae, Lecanicillium lecanii were studied on grape thrips. Mycosis by Beauveria bassiana was confirmed the pathogenicity of entomopathogenic fungi on grape thrips. Highly pronounced mycosis was observed by Metarhizium anisopliae on the dead bodies of thrips. Mycosis test of Lecanicillium lecanii was also proved on grape thrips (plate - 1, 2 and 3).

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=86 SRC="FIGDIR/small/449356v1_figP1.gif" ALT="Figure 1">
View larger version (51K):
org.highwire.dtl.DTLVardef@ce72dorg.highwire.dtl.DTLVardef@11d5d3borg.highwire.dtl.DTLVardef@d3fc3forg.highwire.dtl.DTLVardef@155a62a_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOPlate 1:C_FLOATNO Mycosis of B.eauveria bassiana on grape thrips

C_FIG O_FIG O_LINKSMALLFIG WIDTH=188 HEIGHT=200 SRC="FIGDIR/small/449356v1_figP2.gif" ALT="Figure 2">
View larger version (97K):
org.highwire.dtl.DTLVardef@104ea8org.highwire.dtl.DTLVardef@9c124corg.highwire.dtl.DTLVardef@754939org.highwire.dtl.DTLVardef@1a7d04_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOPlate 2:C_FLOATNO Mycosis of Metarhizium anisopliae on grape thrips

C_FIG O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=186 SRC="FIGDIR/small/449356v1_figP3.gif" ALT="Figure 3">
View larger version (85K):
org.highwire.dtl.DTLVardef@105cf36org.highwire.dtl.DTLVardef@66ebfborg.highwire.dtl.DTLVardef@10ede35org.highwire.dtl.DTLVardef@11e942c_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOPlate 3:C_FLOATNO Mycosis of Lecanicillium lecanii on grape thrips

C_FIG

category: zoology