title: Predicting trait phenotypes from knowledge of the topology of gene networks

abstract: In many fields there is interest in manipulating genes and gene networks to realize improved trait phenotypes. The practicality of doing so, however, requires accepted theory on the properties of gene networks that is well-tested by empirical results. The extension of quantitative genetics to include models that incorporate properties of gene networks expands the long tradition of studying epistasis resulting from gene-gene interactions. Here we consider NK models of gene networks by applying concepts from graph theory and Boolean logic theory, motivated by a desire to model the parameters that influence predictive skill for trait phenotypes under the control of gene networks; N defines the number of graph nodes, the number of genes in the network, and K defines the number of edges per node in the graph, representing the gene-gene interactions. We define and consider the attractor period of an NK network as an emergent trait phenotype for our purposes. A long-standing theoretical treatment of the dynamical properties of random Boolean networks suggests a transition from long to short attractor periods as a function of the average node degree K and the bias probability P in the applied Boolean rules. In this paper we investigate the appropriateness of this theory for predicting trait phenotypes on random and real microorganism networks through numerical simulation. We show that: (i) the transition zone between long and short attractor periods depends on the number of network nodes for random networks; (ii) networks derived from metabolic reaction data on microorganisms also show a transition from long to short attractor periods, but at higher values of the bias probability than in random networks with similar numbers of network nodes and average node degree; (iii) the distribution of phenotypes measured on microorganism networks shows more variation than random networks when the bias probability in the Boolean rules is above 0.75; and (iv) the topological structure of networks built from metabolic reaction data is not random, being best approximated, in a statistical sense, by a lognormal distribution. The implications of these results for predicting trait phenotypes where the genetic architecture of a trait is a gene network are discussed.


title: Epigenetic regulation clocks the multigenerational olfactory imprinting in C. elegans

abstract: Imprinting is an early sensory life experience that induces adult behaviours, such as mother recognition or homing. In a previous study, we demonstrated a striking olfactory imprinting in C. elegans that can be inherited over generations. When exposed to specific odorants during a timely controlled post-hatch period, C. elegans worms display during adulthood an enhanced migration towards these molecules. In order to unveil some of the genetic and epigenetic factors that are responsible for such a behavioural plasticity, we assessed the role of heterochronic genes using a candidate gene approach. We report here that translation of the Hunchback-Like 1 (HBL1) transcription factor in the sensory processing interneuron AIY, is a determining factor for olfactory plasticity timing in C.elegans. HBL1 may associate to the SPR1/CoREST co-repressor, the lysine demethylase SPR5/LSD1 and the histone deacetylase HDA3 to lengthen the plasticity period, whereas the translation initiation factor IFE-4 and the histone deacetylase HDA2 abridge it. We also observed that lengthened plasticity periods allow proportionally faster stable behavioral adaptation of C. elegans populations. We conclude that plasticity timing is a key factor, not only to transiently adapt individuals but also to stably adapt animal populations via multigenerational accumulation of experience.

category: animal behavior and cognition


title: Glutathione synthetase overexpression in Acidithiobacillus ferrooxidans improves halotolerance of iron oxidation

abstract:
Acidithiobacillus ferrooxidans are well-studied iron- and sulfur-oxidizing acidophilic chemolithoautotrophs that are exploited for their ability to participate in the bioleaching of metal sulfides. Here, we overexpressed the endogenous glutamate--cysteine ligase and glutathione synthetase genes in separate strains and found that glutathione synthetase overexpression increased intracellular glutathione levels. We explored the impact of pH on the halotolerance of iron oxidation in wild type and engineered cultures. The increase in glutathione allowed the modified cells to grow under salt concentrations and pH conditions that are fully inhibitory to wild type cells. These results indicate that glutathione overexpression can be used to increase halotolerance in A. ferrooxidans and would likely be a useful strategy on other acidophilic bacteria.

ImportanceThe use of acidophilic bacteria in the hydrometallurgical processing of sulfide ores can enable many benefits including the potential reduction of environmental impacts. The cells involved in bioleaching tend to have limited halotolerance, and increased halotolerance could enable several benefits, including a reduction in the need for fresh water resources. We show that the genetic modification of A. ferrooxidans for the overproduction of glutathione is a promising strategy to enable cells to resist the oxidative stress that can occur during growth in the presence of salt.

category: bioengineering


title: High-throughput single-cell TCR - pMHC dissociation rate measurements performed by an autonomous microfluidic cellular processing unit

abstract: We developed an integrated microfluidic cellular processing unit (mCPU) capable of autonomously isolating single cells, perform, measure, and on-the-fly analyze cell-surface dissociation rates, followed by recovery of selected cells. We performed proof-of-concept, high-throughput single-cell experiments characterizing pMHC-TCR interactions on live CD8+ T cells. The mCPU platform analyzed TCR-pMHC dissociation rates with a throughput of 50 cells per hour and hundreds of cells per run, and we demonstrate that cells can be selected, enriched, and easily recovered from the device.


title: miRNAFinder: A Comprehensive Web Resource for Plant Pre-microRNA Classification

abstract: microRNAs (miRNAs) are known as one of the small non-coding RNA molecules that control the expression of genes at the RNA level, while some operate at the DNA level. They typically range from 20-24 nucleotides in length and can be found in the plant and animal kingdoms as well as in some viruses. Computational approaches have overcome the limitations of the experimental methods and have performed well in identifying miRNAs. Compared to mature miRNAs, precursor miRNAs (pre-miRNAs) are long and have a hairpin loop structure with structural features. Therefore, most in-silico tools are implemented for pre-miRNA identification. This study presents a multilayer perceptron (MLP) based classifier implemented using 180 features under sequential, structural, and thermodynamic feature categories for plant pre-miRNA identification. This classifier has a 92% accuracy, a 94% specificity, and a 90% sensitivity. We have further tested this model with other small non-coding RNA types and obtained 78% accuracy. Furthermore, we introduce a novel dataset to train and test machine learning models, addressing the overlapping data issue in the positive training and testing datasets presented in PlantMiRNAPred for the classification of real and pseudo-plant pre-miRNAs. The new dataset and the classifier that can be used with any plant species are deployed on a web server freely accessible at http://mirnafinder.shyaman.me/.

category: bioinformatics


title: Capturing hidden regulation based on noise change of gene expression level from single cell RNA-seq in yeast

abstract: Recent progress in high throughput single cell RNA-seq (scRNA-seq) has activated the development of data-driven inferring methods of gene regulatory networks. Most network estimations assume that perturbations produce downstream effects. However, the effects of gene perturbations are sometimes compensated by a gene with redundant functionality (functional compensation). In order to avoid functional compensation, previous studies constructed double gene deletions, but its vast nature of gene combinations was not suitable for comprehensive network estimation. We hypothesized that functional compensation may emerge as a noise change without mean change (noise-only change) due to varying physical properties and strong compensation effects. Here, we show compensated interactions, which are not detected by mean change, are captured by noise-only change quantified from scRNA-seq. We investigated whether noise-only change genes caused by a single deletion of STP1 and STP2, which have strong functional compensation, are enriched in redundantly regulated genes. As a result, noise-only change genes are enriched in their redundantly regulated genes. Furthermore, novel downstream genes detected from noise change are enriched in "transport", which is related to known downstream genes. Herein, we suggest the noise difference comparison has the potential to be applied as a new strategy for network estimation that capture even compensated interaction.


title: PCaDB - a comprehensive and interactive database for transcriptomes from prostate cancer population cohorts

abstract: Prostate cancer (PCa) is a heterogeneous disease with highly variable clinical outcomes which presents enormous challenges in the clinical management. A vast amount of transcriptomics data from large PCa cohorts have been generated, providing extraordinary opportunities for the molecular characterization of the PCa disease and the development of diagnostic and prognostic signatures. The lack of an inclusive collection and harmonization of the scattered public datasets constrains the extensive use of the valuable resources. In this study, we present a user-friendly database, PCaDB, for a comprehensive and interactive analysis and visualization of gene expression profiles from 77 transcriptomics datasets with 9,068 patient samples. PCaDB also includes a single-cell RNA-sequencing (scRNAseq) dataset for normal human prostates and 30 published PCa prognostic signatures. The comprehensive data resources and advanced analytical methods equipped in PCaDB would greatly facilitate data mining to understand the heterogeneity of PCa and to develop machine learning models for accurate PCa diagnosis and prognosis to assist on clinical decision-making. PCaDB is publicly available at http://bioinfo.jialab-ucr.org/PCaDB/.


title: The Genome Sequence Archive Family: Towards Explosive Data Growth and Diverse Data Types

abstract: The Genome Sequence Archive (GSA) is a data repository for archiving raw sequence data, which provides data storing and sharing services for worldwide scientific communities. Considering explosive data growth with diverse data types, here we present the GSA family by expanding into a set of resources for raw data archive with different purposes, namely, GSA (https://ngdc.cncb.ac.cn/gsa/), GSA for Human (GSA-Human, https://ngdc.cncb.ac.cn/gsa-human/), and Open Archive for Miscellaneous Data (OMIX, https://ngdc.cncb.ac.cn/omix/). Compared with the 2017 version, GSA has been significantly updated in data model, online functionalities, and web interfaces. GSA-Human, as a new partner of GSA, is a data repository specialized in human genetics-related data with controlled access and security. OMIX, as a critical complement to the two resources mentioned above, is an open archive for miscellaneous data. Together, all these resources form a family of resources dedicated to archiving explosive data with diverse types, accept data submissions from all over the world and provide free open access to all publicly available data in support of worldwide research activities.


title: spiralize: an R Package for Visualizing Data on Spirals

abstract:
SummarySpiral layout has two major advantages for data visualization. First, it is able to visualize data with long axes, which greatly improves the resolution of visualization. Second, it is efficient for time series data to reveal periodic patterns. Here we present the R package spiralize that provides a general solution for visualizing data on spirals. spiralize implements numerous graphics functions so that self-defined high-level graphics can be easily implemented by users. The flexibility and power of spiralize are demonstrated by five examples from real-world datasets.

Availability and implementationThe spiralize package and documentations are freely available at the Comprehensive R Archive Network (CRAN) https://CRAN.R-project.org/package=spiralize.

Contactz.gu@dkfz.de or d.huebschmann@dkfz.de

Supplementary informationSupplementary data are available at Bioinformatics online.


title: Improved estimation of cell type-specific gene expression through deconvolution of bulk tissues with matrix completion

abstract: Single cell RNA-seq (scRNA-seq) has been widely used to uncover cellular heterogeneity, however, the constraints of cost make it impractical as a routine on large patient cohorts. Here we present ENIGMA, a method that accurately deconvolute bulk tissue RNA-seq into single cell-type resolution given the knowledge gained from scRNA-seq. ENIGMA applies a matrix completion strategy to minimize the distance between mixture transcriptome and weighted combination of cell type-specific expression, allowing quantification of cell type proportions and reconstruction of cell type-specific transcriptome. The superior performance of ENIGMA was validated in simulated and realistic datasets, including disease-related tissues, demonstrating its ability in novel biological findings.


title: Accurate estimation of cell-type resolution transcriptome in bulk tissue through matrix completion


title: Evaluation of FRET X for Single-Molecule Protein Fingerprinting

abstract: Single-molecule protein identification is a novel, as of yet unrealized concept with potentially groundbreaking applications in biological research. We propose a method called FRET X (Forster Resonance Energy Transfer via DNA eXchange) fingerprinting, in which the FRET efficiency is read out between exchangeable dyes on protein-bound DNA docking strands, and accumulated FRET efficiency values constitute the fingerprint for a protein. To evaluate the feasibility of this approach, we simulated fingerprints for hundreds of proteins using a coarse-grained lattice model and experimentally demonstrated FRET X fingerprinting on a system of model peptides. Measured fingerprints are in agreement with our simulations, corroborating the validity of our modeling approach. In a simulated complex mixture of >300 human proteins of which only cysteines, lysines and arginines were labeled, a support vector machine was able to identify constituents with 95% accuracy. We anticipate that our FRET X fingerprinting approach will form the basis of an analysis tool for targeted proteomics.

category: biophysics


title: Intratumoral Heterogeneity Promotes Collective Cancer Invasion Through NOTCH1 Variation

abstract:
Cellular and molecular heterogeneity within tumors has long been associated with the progression of cancer to an aggressive phenotype and a poor prognosis. However, how such intratumoral heterogeneity contributes to the invasiveness of cancer is largely unknown. Here, using a multidisciplinary approach, we investigate the interaction between molecular subtypes within bladder microtumors and the corresponding effects on their invasiveness. Our results reveal heterogeneous microtumors formed by multiple molecular subtypes possess enhanced invasiveness compared to individual cells, even when both cells are not invasive individually. To examine the molecular mechanism of intratumoral heterogeneity mediated invasiveness, live single cell biosensing, RNA interference, and CRISPR-Cas9 gene editing approaches were applied to investigate and control the composition of the microtumors. An agent-based computational model was also developed to evaluate the influence of NOTCH1 variation on DLL4 expression within a microtumor. The data indicate that variation in NOTCH1 expression can lead to upregulation of DLL4 expression within the microtumor and enhancement of microtumor invasiveness. Overall, our results reveal a novel mechanism of heterogeneity mediated invasiveness through intratumoral variation of gene expression.

Summary statementThis study reveals a mechanism that Notch1 variation, instead of the average value, promotes the invasiveness of microtumor, providing a link between intratumoral heterogeneity and collective cancer invasion.

category: cancer biology


title: Conversion of cancer-associated fibroblasts from pro- to antitumor improves the sensitivity of pancreatic cancer to chemotherapeutics

abstract: Previous therapeutic attempts to deplete cancer-associated fibroblasts (CAFs) or inhibit their proliferation in pancreatic ductal adenocarcinoma (PDAC) were not successful in mice or patients. Thus, CAFs may be tumor suppressive or heterogeneous, with distinct cancer-restraining and -promoting CAFs (rCAFs and pCAFs, respectively). Here, we show that induced expression of the glycosylphosphatidylinositol-anchored protein Meflin, a rCAF-specific marker, in CAFs by genetic and pharmacological approaches improved the chemosensitivity of mouse PDAC. A chemical library screen identified Am80, a synthetic, non-natural retinoid, as a reagent that effectively induced Meflin expression in CAFs. Am80 administration improved the sensitivity of PDAC to chemotherapeutics, accompanied by increases in tumor vessel area and intratumoral drug delivery. Mechanistically, Meflin was involved in the suppression of tissue stiffening by interacting with lysyl oxidase to inhibit its collagen crosslinking activity. These data suggested that modulation of CAF heterogeneity may represent a strategy for PDAC treatment.


title: Toxicological profile of umbilical cord blood-derived small extracellular vesicles

abstract: The development and adoption of cell therapies has been largely limited by difficulties associated with their safety, handling and storage. Extracellular vesicles (EV) have recently emerged as a likely mediator for the therapeutic effect of cells, offering several advantages over cell therapies. Due to their small size and inability to expand and metastasize, EV are generally considered safer than cell transplantation. Nevertheless, few studies have scrutinized the toxicity profile of EV, particularly after repeated high dose administration. The present study aimed to evaluate a preparation of small EV obtained from umbilical cord blood mononuclear cells (UCB-MNC-sEV) for its cytotoxicity in different cell lines, as well as its differential accumulation, distribution and toxicity following repeated intravenous (IV) administrations in a rodent model. In vitro, repeated sEV exposure in concentrations up to 1x1011 particles/ml had no deleterious impact on the viability or metabolic activity of peripheral blood mononuclear cells, THP-1 monocytes, THP-1-derived macrophages, normal dermal human fibroblasts or human umbilical vein endothelial cells. DiR-labeled sEV, injected IV for four weeks in healthy rats, were detected in clearance organs, particularly kidneys, spleen and liver, similarly to control dye. Moreover, repeated administrations during six and twelve weeks of up to 1x1010 total particles of sEV-dye were well tolerated, with no changes in general hematological cell counts, or kidney and liver toxicity markers. Importantly, unlabeled sEV likewise did not induce significant alterations in cellular and biochemical blood parameters, nor any morphological changes in heart, kidney, lung, spleen, or liver tissue. In sum, our data shows that UCB-MNC-sEV have no significant toxicity in vitro or in vivo, even when administered repeatedly at high concentrations, therefore confirming their safety profile and potential suitability for future clinical use.

category: pharmacology and toxicology


title: Susceptibility to diet induced obesity at thermoneutral conditions is independent of UCP1

abstract:
ObjectiveActivation of uncoupling protein 1 (UCP1) in brown adipose tissue (BAT) upon cold stimulation leads to substantial increase in energy expenditure to defend body temperature. Increases in energy expenditure after a high caloric food intake, termed diet-induced thermogenesis, are also attributed to BAT. These properties render BAT a potential target to combat diet-induced obesity. However, studies investigating the role of UCP1 to protect against diet-induced obesity are controversial and rely on the phenotyping of a single constitutive UCP1-knockout model.

To address this issue, we generated a novel UCP1-knockout model by Cre-mediated deletion of Exon 2 in the UCP1 gene. We studied the effect of constitutive UCP1 knockout on metabolism and the development of diet-induced obesity.

MethodsUCP1 knockout and wildtype mice were housed at 30{degrees}C and fed a control diet for 4-weeks followed by 8-weeks of high-fat diet. Body weight and food intake were monitored continuously over the course of the study and indirect calorimetry was used to determine energy expenditure during both feeding periods.

ResultsBased on Western blot analysis, thermal imaging and noradrenaline test, we confirmed the lack of functional UCP1 in knockout mice. However, body weight gain, food intake and energy expenditure were not affected by deletion of UCP1 gene function during both feeding periods.

ConclusionConclusively, we show that UCP1 does not protect against diet-induced obesity at thermoneutrality. Further we introduce a novel UCP1-KO mouse enabling the generation of conditional UCP1-knockout mice to scrutinize the contribution of UCP1 to energy metabolism in different cell types or life stages.

category: physiology


title: Integrative analysis of the shikonin metabolic network identifies new gene connections and reveals evolutionary insight into shikonin biosynthesis

abstract: Plant specialized 1,4-naphthoquinones present a remarkable case of convergent evolution. Species across multiple discrete orders of vascular plants produce diverse 1,4-naphthoquinones via one of several pathways using different metabolic precursors. Evolution of these pathways was preceded by events of metabolic innovation and many appear to share connections with biosynthesis of photosynthetic or respiratory quinones. Here, we sought to shed light on the metabolic connections linking shikonin biosynthesis with its precursor pathways and on the origins of shiknoin metabolic genes. Downregulation of Lithospermum erythrorhizon geranyl diphosphate synthase (LeGPPS), recently shown to have been recruited from a cytoplasmic farnesyl diphosphate synthase (FPPS), resulted in reduced shikonin production and a decrease in expression of mevalonic acid and phenylpropanoid pathway genes. Next, we used LeGPPS and other known shikonin pathway genes to build a coexpression network model for identifying new gene connections to shikonin metabolism. Integrative in silico analyses of network genes revealed candidates for biochemical steps in the shikonin pathway arising from Boraginales-specific gene family expansion. Multiple genes in the shikonin coexpression network were also discovered to have originated from duplication of ubiquinone pathway genes. Taken together, our study provides evidence for transcriptional crosstalk between shikonin biosynthesis and its precursor pathways, identifies several shikonin pathway gene candidates and their evolutionary histories, and establishes additional evolutionary links between shikonin and ubiquinone metabolism. Moreover, we demonstrate that global coexpression analysis using limited transcriptomic data obtained from targeted experiments is effective for identifying gene connections within a defined metabolic network.

category: plant biology


title: Determination of the Optical Properties of Thylakoid Membranes from Whole Leaf Reflectance Measurements Using a Capacitor and Unbound Electron Model

abstract: This paper demonstrates that a capacitor equivalent along with unbound electrons can be used to model thylakoid membranes in grana stacks. From whole leaf reflectance measurements at normal incidences at 660nm wavelength and taken from the literature (Hu et al. 2020), refractive indices are obtained from the Fresnels equation for transverse electric (TE) and transverse magnetic (TM) polarization. The TE and TM polarizations for external reflectance depict the Brewster angle at which the magnitude of the reflected electric vector is zero; the internal reflections show that theres a narrow angle window of about 10 degrees before the internally refracted light goes into critical angle. The clustering and separation of reflection measurements with angle of incidence is explained using Fresnels equation; the cross-over angle beyond the Brewsters angle for internal reflection at which reflectance measurements separate with incidence angle is given by [Formula]. The predicted relaxation times from a capacitor and unbound electron model gave favorable comparisons against commonly reported fluorescence times in the 0.1 to 1 ns range (our results gave 0.5-0.8 ns). The di-electric constant for the membrane is estimated to be 5. The stacking number (number of grana layers) is consistent with the light penetration depth (skin depth). The magnetic permeability was shown to be close to that of vacuum and therefore the thylakoid lacks any magnetic properties as would be expected for such a transparent media. An in-vivo estimate based on thermal equilibrium of molecules for the permanent dipole moment of the chlorophyll molecule gave 2,025D (Debye).


title: A 3D gene expression atlas of the floral meristem based on spatial reconstruction of single nucleus RNA sequencing data

abstract: Identity and functions of plant cells are influenced by their precise cellular location within the plant body. Cellular heterogeneity in growth and differentiation trajectories results in organ patterning. Therefore, assessing this heterogeneity at molecular scale is a major question in developmental biology. Single-cell transcriptomics (scRNA-seq) allows to characterize and quantify gene expression heterogeneity in developing organs at unprecedented resolution. However, the original physical location of the cell is lost during the scRNA-seq procedure. To recover the original location of cells is essential to link gene activity with cellular function and morphology. Here, we reconstruct genome-wide gene expression patterns of individual cells in a floral meristem by combining single-nuclei RNA-seq with 3D spatial reconstruction. By this, gene expression differences among meristematic domains giving rise to different tissue and organ types can be determined. As a proof of principle, the data are used to trace the initiation of vascular identity within the floral meristem. Our work demonstrates the power of spatially reconstructed single cell transcriptome atlases to understand plant morphogenesis. The floral meristem 3D gene expression atlas can be accessed at http://threed-flower-meristem.herokuapp.com


title: Characterization of the Brassica napus flavonol synthase gene family reveals bifunctional flavonol synthases

abstract: Flavonol synthase (FLS) is a key enzyme for the formation of flavonols, which are a subclass of the flavonoids. FLS catalyses the conversion of dihydroflavonols to flavonols. The enzyme belongs to the 2-oxoglutarate-dependent dioxygenases (2-ODD) superfamily. We characterized the FLS gene family of Brassica napus that covers 13 genes, based on the genome sequence of the B. napus cultivar Express 617. The goal was to unravel which BnaFLS genes are relevant for seed flavonol accumulation in the amphidiploid species B. napus. Two BnaFLS1 homeologs were identified and shown to encode bifunctional enzymes. Both exhibit FLS activity as well as flavanone 3-hydroxylase (F3H) activity, which was demonstrated in vivo and in planta. BnaFLS1-1 and -2 are capable of converting flavanones into dihydroflavonols and further into flavonols. Analysis of spatio-temporal transcription patterns revealed similar expression profiles of BnaFLS1 genes. Both are mainly expressed in reproductive organs and co-expressed with the genes encoding early steps of flavonoid biosynthesis. Our results provide novel insights into flavonol biosynthesis in B. napus and contribute information for breeding targets with the aim to modify the flavonol content in rapeseed.


title: Intrinsic Network Activity Reflects the Fluctuating Experience of Tonic Pain

abstract:
Although we know sensation is continuous, research on long-lasting and continuously changing stimuli is scarce and the dynamic nature of ongoing cortical processing is largely neglected.

In a longitudinal study with 152 fMRI sessions, participants were asked to continuously rate the intensity of applied tonic heat pain for 20 minutes. Using group independent component analysis and dual-regression, we extracted the subjects time courses of intrinsic network activity. The relationship between the dynamic fluctuation of network activity with the varying time courses of three pain processing entities was computed: pain intensity, the direction of pain intensity changes and temperature.

We were able to dissociate the spatio-temporal patterns of objective (temperature) and subjective (pain intensity/changes of pain intensity) aspects of pain processing in the human brain. We found two somatosensory networks with distinct functions: one network which encodes the small fluctuations in temperature and consists mainly of bilateral SI. A second right-lateralised network that encodes the intensity of the subjective experience of pain consists of SI, SII, the PCC, and the thalamus.

We revealed the somatosensory dynamics that build up towards a current subjective percept of pain. The timing suggests a cascade of subsequent processing steps towards the current pain percept.

category: neuroscience


title: A model of rapid homeostatic plasticity accounts for hidden, long-lasting changes in a neuronal circuit after exposure to high potassium.

abstract: Neural circuits must both function reliably and flexibly adapt to changes in their environment. We studied how both biological neurons and computational models respond to high potassium concentrations. Pyloric neurons of the crab stomatogastric ganglion (STG) initially become quiescent, then recover spiking activity in high potassium saline. The neurons retain this adaptation and recover more rapidly in subsequent high potassium applications, even after hours in control saline. We constructed a novel activity-dependent computational model that qualitatively captures these results. In this model, regulation of conductances is gated on and off depending on how far the neuron is from its target activity. This allows the model neuron to retain a trace of past perturbations even after it returns to its target activity in control conditions. Thus, perturbation, followed by recovery of normal activity, can hide cryptic changes in neuronal properties that are only revealed by subsequent perturbations.


title: BOLD fMRI signals of visual white matter encode visuotopic information and predict effective connectivity between visual areas

abstract:
The functional significance of BOLD signals in white matter (WM) remains unclear. The current study investigated whether 7T BOLD fMRI signal from visual WM tracts contains high fidelity retinotopic information and whether it correlates with the effective connectivity between visual areas. Population receptive field (pRF) analysis of the 7T retinotopy dataset from the Human Connectome Project revealed clear contralateral retinotopic representations from two visual WM bundles: optic radiation (OR) and vertical occipital fasciculus (VOF). The retinotopic organization of OR is consistent with post-mortem studies. The pRF size of WM voxels also increases with eccentricity. Based on the retinotopic maps of OR, we investigated whether BOLD signals in OR during visual stimulation are related to the resting-state effective connectivity between the lateral geniculate nucleus (LGN) and the primary visual cortex (V1). Results show that visually-evoked BOLD responses in OR correlate with the feedforward and feedback connectivity between the LGN and V1 during resting state. These findings demonstrate that WM BOLD signals contain high fidelity information such as visual field maps, and also predict the functional connectivity of brain areas.

Significance statementWhite matter (WM) tracks conduct spiking activity between distant neurons. Weak fluctuations of BOLD signals in the WM can be detected with fMRI, but their functional relevance remains largely unknown. Here we characterized the visual field map properties of two major visual WM bundles: the optic radiation (OR) and vertical occipital fasciculus (VOF). Population receptive field analysis of the WM BOLD signals revealed clear visual field maps in both WM tracts. Effective connectivity analysis further showed that visually evoked BOLD responses in OR can predict the resting thalamo-cortical functional connectivity. These findings demonstrate that WM BOLD signals contain highly specific functional information and could directly index the functional connectivity between brain areas.


title: Individually Unique Dynamics of Cortical Connectivity Reflect the Ongoing Intensity of Chronic Pain

abstract:
BackgroundChronic pain diseases are characterised by an ongoing and fluctuating endogenous pain, yet it remains to be elucidated how this is reflected by the dynamics of ongoing functional cortical connections. The present study addresses this disparity by taking the individual perspective of pain patients into account, which is the varying intensity of endogenous pain.

MethodsTo this end, we investigated the cortical encoding of 20 chronic back pain patients and 20 chronic migraineurs in four repeated fMRI sessions. During the recording, the patients were asked to continuously rate their pain intensity. A brain parcellation approach subdivided the whole brain into 408 regions. A 10 s sliding-window connectivity analysis computed the pair-wise and time-varying connectivity between all brain regions across the entire recording period. Linear mixed effects models were fitted for each pair of brain regions to explore the relationship between cortical connectivity and the observed trajectory of the patients fluctuating endogenous pain.

ResultsTwo pain processing entities were taken into account: pain intensity (high, middle, low pain) and the direction of pain intensity changes (rising vs. falling pain). Overall, we found that periods of high and increasing pain were predominantly related to low cortical connectivity. For chronic back pain this applies to the pain intensity-related connectivity for limbic and cingulate areas, and for the precuneus. The change of pain intensity was subserved by connections in left parietal opercular regions, right insular regions, as well as large parts of the parietal, cingular and motor cortices. The change of pain intensity direction in chronic migraine was reflected by decreasing connectivity between the anterior insular cortex and orbitofrontal areas, as well as between the PCC and frontal and ACC regions.

ConclusionsInterestingly, the group results were not mirrored by the individual patterns of pain-related connectivity, which is suggested to deny the idea of a common neuronal core problem for chronic pain diseases. In a similar vein, our findings are supported by the experience of clinicians, who encounter patients with a unique composition of characteristics: personality traits, various combinations of symptoms, and a wide range of individual responses to treatment. The diversity of the individual cortical signatures of chronic pain encoding results adds to the understanding of chronic pain as a complex and multifaceted disease. The present findings support recent developments for more personalised medicine.


title: Deficiency of the ywhaz gene, involved in neurodevelopmental disorders, alters brain activity and behaviour in zebrafish

abstract: Genetic variants in YWHAZ contribute to psychiatric disorders such as autism spectrum disorder and schizophrenia, and have been related to an impaired neurodevelopment in humans and mice. Here, we used zebrafish to further understand the mechanisms by which YWHAZ contributes to neurodevelopmental disorders. We observed that ywhaz expression was panneuronal during developmental stages and restricted to Purkinje cells in the adult cerebellum, cells that are described to be reduced in number in autistic patients. We then performed whole-brain imaging in wild-type and ywhaz CRISPR/Cas9 knockout (KO) larvae and found altered neuronal activity and connectivity in the hindbrain. Adult ywhaz KO fish display decreased levels of monoamines in the hindbrain and freeze when exposed to novel stimuli, a phenotype that can be reversed with drugs that target monoamine neurotransmission. These findings suggest an important role for ywhaz in establishing neuronal connectivity during development and modulating both neurotransmission and behaviour in adults.


title: Genomic Modules and Intramodular Network Convergency of Susceptibility and Resilience in Multimodeled Stress in Male Mice

abstract: The multifactorial etiology of stress-related disorders is a challenge in developing synchronized medical standards for treatment and diagnosis. It is largely unknown whether there exists molecular convergence in preclinical models of stress generated using disparate construct validity. Using RNA-sequencing (RNA-seq), we investigated the genomic signatures in the ventral hippocampus, which mostly regulates affective behavior, in mouse models that recapitulate the hallmarks of anxiety and depression. Chronic oral corticosterone (CORT), a model that causes a blunted endocrine response to stress, induced anxiety- and depression-like behavior in wildtype mice and mice heterozygous for the gene coding for brain-derived neurotrophic factor (BDNF) Val66Met, a variant associated with genetic susceptibility to stress. In a separate set of mice, chronic social defeat stress led to a susceptible or a resilient population, whose proportion was dependent on housing conditions, standard housing or enriched environment. A rank-rank-hypergeometric (RRHO) analysis of the RNA-seq data set across models demonstrated that in mice treated with CORT and susceptible mice raised in standard housing differentially expressed genes (DEGs) converged toward gene networks involved in similar biological functions. Weighted gene co-expression analysis generated 54 unique modules of interconnected gene hubs, two of which included a combination of all experimental groups and were significantly enriched in DEGs, whose function was consistent with that predicted in the RRHO GO analysis. This multimodel approach showed transcriptional synchrony between models of stress with hormonal, environmental or genetic construct validity shedding light on common genomic drivers that embody the multifaceted nature of stress-related disorders.


title: Early neuroadaptations to an obesogenic diet identify the schizophrenia-related ErbB4 receptor in obesity-induced hippocampal abnormalities

abstract:
Childhood obesity leads to hippocampal atrophy and altered cognition. However, the molecular mechanisms underlying these impairments are poorly understood. The neurotrophic factor neuregulin-1 (NRG1) and its cognate ErbB4 receptor play critical roles in hippocampal maturation and function. This study aimed to determine whether altered NRG1-ErbB4 activities may partly explain hippocampal abnormalities in rats exposed to an obesogenic Western-like diet (WD). Lewis rats were randomly divided into four groups (12 rats/group): 1) control diet+vehicle (CDV); 2) CD+NRG1 (CDN) (daily intraperitoneal injections: 5 g/kg/day; between postnatal day, PND 21-PND 41); 3) WD+VEH (WDV); 4) WD+NRG1 (WDN). Neurobehavioral assessments were performed at PND 43-49. Brains were harvested for MRI and molecular analyses at PND 49. We found that NRG1 administration reduced hippocampal volume (7%) and attenuated hippocampal-dependent cued fear conditioning in CD rats (56%). NRG1 administration reduced PSD-95 protein expression (30%) and selectively reduced hippocampal cytokine levels (IL-33, GM-CSF, CCL-2, IFN-{gamma}) while significantly impacting microglia morphology (increased span ratio and reduced circularity). WD rats exhibited reduced right hippocampal volume (7%), altered microglia morphology (reduced density and increased lacunarity), and increased levels of cytokines implicated in neuroinflammation (IL-1, TNF-, IL-6). Notably, NRG1 synergized with the WD to increase hippocampal ErbB4 phosphorylation and the tumor necrosis alpha converting enzyme (TACE/ADAM17) protein levels. Together, these data suggest a novel interaction between obesogenic diet exposure and TACE/ADAM17-NRG1-ErbB4 signaling during hippocampal maturation. Our results indicate that supraoptimal ErbB4 activities may contribute to the abnormal hippocampal structure and cognitive vulnerabilities observed in obese individuals.

HighlightsO_LIObesogenic diet consumption during adolescence induces anxiety-like behaviors before the onset of obesity-related changes in metabolism.
C_LIO_LIObesogenic diet-driven abnormal behaviors co-occurred with alterations in hippocampal pro-inflammatory cytokine profiles.
C_LIO_LIObesogenic diet consumption attenuates exogenous NRG1 effects on hippocampal-related behaviors and structure.
C_LIO_LIExogenous NRG1 administration during adolescence resulted in reduced hippocampal volumes and domain-specific cognitive impairments.
C_LIO_LIExogenous NRG1 administration has potent immunomodulatory actions and alters hippocampal microglia morphology.
C_LI


title: Prenatal alcohol exposure disrupts hippocampal sharp-wave ripple-associated spike dynamics

abstract: Prenatal alcohol exposure (PAE) is among the most common developmental insults to the nervous system and is characterized by memory disruption. There is a pressing need to identify physiological alterations that help explain this memory impairment. Hippocampal sharp-wave ripples (SPW-Rs) are a compelling candidate for this purpose as they are the electrophysiological signatures of memory consolidation. We report that rats exposed to moderate prenatal alcohol display abnormalities restricted to SPW-R episodes that manifest as decreased recruitment of CA1 pyramidal cells and interneurons to SPW-R events, altered excitation during SPW-Rs, and decreased cell assembly activation rate. These differences observed at the single neuron and the population level may limit the ability of memory trace reactivation during SPW-Rs through the disruption of the intrinsic structure of cell sequences. Together, our results suggest that alterations in hippocampal SPW-R spike dynamics may underlie alcohol exposure-related memory deficits.


title: ICER controls HIV-1 infection and replication in elite controllers

abstract: A rare subset of HIV-infected individuals, termed elite controllers (ECs), can maintain long-term control over HIV replication in the absence of antiretroviral therapy (ART). To elucidate the biological mechanism of resistance to HIV replication at the molecular and cellular levels, we performed RNA sequencing and identified alternative splicing variants from ECs, HIV-infected individuals undergoing ART, ART-naive HIV-infected individuals, and healthy controls. Differential gene expression patterns that are specific to ECs and may influence HIV resistance were identified, including alternative RNA splicing and exon usage variants of the CREM/ICER gene (cAMP-responsive element modulator/inducible cAMP early repressors). The knockout and knockdown of specific ICER exons that were found to be upregulated in ECs resulted in significantly increased HIV infection in CD4+ T cell line and primary CD4+ T cells. Overexpression of ICER isoforms decreased HIV infection in primary CD4+ T cells. Furthermore, ICER regulated HIV-1 LTR promoter activity in a Tat-dependent manner. Together, these results suggest that ICER is an HIV host factor that may contribute to HIV resistance of ECs. These findings will help elucidate the mechanisms of HIV control by ECs and may yield a new approach for treatment of HIV.

category: immunology


title: Synthetic dual-function RNA reveals features necessary for target regulation

abstract:
Small base pairing RNAs (sRNAs) and small proteins comprise two classes of regulators that allow bacterial cells to adapt to a wide variety of growth conditions. A limited number of transcripts encoding both of these activities, regulation of mRNA expression by base pairing and a small regulatory protein, have been identified. Given that few have been characterized, little is known about the interplay between the two regulatory functions. To investigate the competition between the two activities, we constructed synthetic dual-function RNAs, hereafter referred to as MgtSR or MgtRS, comprised of the Escherichia coli sRNA MgrR and the open reading frame encoding the small protein MgtS. MgrR is a 98 nt base pairing sRNA that negatively regulates eptB encoding phosphoethanolamine transferase. MgtS is a 31 aa small inner membrane protein that is required for the accumulation of MgtA, a magnesium (Mg2+) importer. Expression of the separate genes encoding MgrR and MgtS is normally induced in response to low Mg2+ by the PhoQP two-component system. By generating various versions of this synthetic dual-function RNA, we probed how the organization of components and the distance between the coding and base pairing sequences contribute to the proper function of both activities of a dual-function RNA. By understanding the features of natural and synthetic dual-function RNAs, future synthetic molecules can be constructed to maximize their regulatory impact.

IMPORTANCEDual-function RNAs encode a small protein and also base pair with mRNAs to act as small, regulatory RNAs. Given that only a limited number of dual-function RNAs have been characterized, further study of these regulators is needed to increase understanding of their features. This study demonstrates that a functional synthetic dual-regulator can be constructed from separate components and used to study the functional organization of dual-function RNAs, with the goal of exploiting these regulators.

category: microbiology


title: A transferable IncC-IncX3 hybrid plasmid co-carrying blaNDM-4, tet(X4), and tmexCD3-toprJ3 confers resistance to carbapenem and tigecycline

abstract: Tigecycline is a last-resort antimicrobial that exhibits promising activity against carbapenemase-producing Enterobacterales (CPE). However, mobile tigecycline resistance genes, tet(X) and tmexCD-toprJ, have emerged in China and have spread possibly worldwide. Tet(X) family proteins, Tet(X3) to Tet(X14), function as tigecycline-inactivating enzymes, and TMexCD-TOprJ complexes function as efflux pumps for tigecycline. Here, we report a CPE isolate co-harboring both emerging tigecycline resistance factors for the first time. A carbapenem- and tigecycline-resistant Klebsiella aerogenes NUITM-VK5 was isolated from an urban drainage in Vietnam in 2021 and a plasmid pNUITM-VK5_mdr co-carrying tet(X4) and tmexCD3-toprJ3 along with the carbapenemase gene blaNDM-4 was identified in NUITM-VK5. pNUITM-VK5_mdr was transferred to Escherichia coli by conjugation and simultaneously conferred high-level resistance against multiple antimicrobials, including carbapenems and tigecycline. An efflux pump inhibitor canceled TMexCD3-TOprJ3-mediated tigecycline resistance, suggesting that both tigecycline resistance factors independently and additively contribute to the high-level resistance. The plasmid had the IncX3 and IncC replicons and was estimated to be a hybrid of plasmids with different origins. Unlike IncX3 plasmids, IncC plasmids are stably maintained in an extremely broad range of bacterial hosts in humans, animals, and environment. Thus, future global spread of multidrug-resistance plasmids such as pNUITM-VK5_mdr poses a public health crisis.


title: Mitotic chromosome condensation requires phosphorylation of the centromeric protein KNL-2 in C. elegans

abstract:
Centromeres are chromosomal regions that serve as sites for kinetochore formation and microtubule attachment, processes that are essential for chromosome segregation during mitosis. Centromeres are almost universally defined by the histone variant CENP-A. In the holocentric nematode C. elegans, CENP-A deposition depends on the loading factor KNL-2. Depletion of either CENP-A or KNL-2 results in defects in centromere maintenance, chromosome condensation and kinetochore formation, leading to chromosome segregation failure. Here, we show that KNL-2 is phosphorylated by CDK-1, and that mutation of three C-terminal phosphorylation sites causes chromosome segregation defects and an increase in embryonic lethality. In strains expressing phosphodeficient KNL-2, CENP-A and kinetochore proteins are properly localised, indicating that the role of KNL-2 in centromere maintenance is not affected. Instead, the mutant embryos exhibit reduced mitotic levels of condensin II on chromosomes and significant chromosome condensation impairment. Our findings separate the functions of KNL-2 in CENP-A loading and chromosome condensation and demonstrate that KNL-2 phosphorylation regulates the cooperation between centromeric regions and the condensation machinery in C. elegans.

SUMMARY STATEMENTPhosphorylation of the essential centromere protein KNL-2 is required for mitotic chromosome condensation, but not for the role of KNL-2 in centromere maintenance and kinetochore formation.

category: cell biology


title: Life history mediates the trade-offs among different components of demographic resilience

abstract: Summary paragraphAccelerating rates of biodiversity loss underscore the need to understand how species achieve resilience -their ability to resist and recover from a/biotic disturbances. Yet, the factors determining the resilience of species remain poorly understood, due to disagreements on its definition and the lack of large-scale analyses. Here, we investigate how the life history of 785 natural populations of animals and plants predict their intrinsic ability to be resilient. We show that demographic resilience can be achieved through different combinations of compensation, resistance, and recovery after a disturbance. We demonstrate that these resilience components are highly correlated with life history traits related to the species pace of life and reproductive strategy. Species with longer generation times require longer recovery times post-disturbance, while those with greater reproductive capacity have greater resistance and compensation. Our findings highlight the key role of life history traits to understand species resilience, improving our ability to predict how natural populations cope with disturbance regimes.

category: ecology


title: Mega-fire in Redwood Tanoak Forest Reduces Bacterial and Fungal Richness and Selects for Pyrophilous Taxa and Traits that are Phylogenetically Conserved

abstract: Mega-fires of unprecedented size, intensity, and socio-economic impacts have surged globally due to climate change, fire suppression, and development. Soil microbiomes are critical for post-fire plant regeneration and nutrient cycling, yet how mega-fires impact the soil microbiome remains unclear. We had a serendipitous opportunity to obtain pre- and post-fire soils from the same sampling locations because the 2016 Soberanes Fire, a mega-fire burning >500 Km2, burned with high severity throughout several of our established redwood-tanoak plots. This makes our study the first to examine microbial fire response in redwood-tanoak forests. We re-sampled soils immediately post-fire from two burned plots and one unburned plot to elucidate the effect of mega-fire on soil microbiomes. We used Illumina MiSeq sequencing of 16S and ITS1 to determine that both bacterial and fungal richness were reduced by 38-70% in burned plots, with richness unchanged in the unburned plot. Fire altered composition by 27% for bacteria and 24% for fungi, whereas the unburned plots experienced no change in fungal and negligible change in bacterial composition. We observed several pyrophilous taxa previously observed in Pinaceae forests, indicating that these microbes are likely general fire-responders across forest types. Further, the pyrophilous taxa that positively responded to fire were phylogenetically conserved, suggesting shared evolutionary traits. For bacteria, fire selected for increased Firmicutes and Actinobacteria. For fungi, fire selected for the Ascomycota classes Pezizomycetes and Eurotiomycetes and for a Basidiomycota class of heat-resistant Geminibasidiomycete yeasts. We hypothesize that microbes share analogous fire response to plants and propose a trait-based conceptual model of microbial response to fire that builds from Grimes Competitor-Stress tolerator-Ruderal framework (C-S-R) and its recent applications to microbes. Using this framework and established literature on several microbial species, we hypothesize some generalizable principals to predict which microbial taxa will respond to fire.


title: Microbial community composition of translocated ancient woodland soil: a case study

abstract:
Soil translocation is an ecological habitat restoration technique which consists of moving the entire topsoil from a donor site to a chosen receptor site. We investigated changes in soil chemistry and microbiology three years following the salvage of semi-ancient woodland soil and materials (0.94 ha) to a nearby receptor pasture due to road widening works (Kent, UK). We sampled i) intact woodland soils adjacent to the area of soils that was translocated to represent the lost donor site, ii) the soil three years after it had been translocated, and iii) grassland soils adjacent to the translocated soil to represent the original receptor site. The intention was to ascertain if shifts in soil chemistry and microbial community composition (Phospholipid Fatty-acid analysis: PLFA) occurred due to soil translocation. PLFA signature biomarkers demonstrated the overall microbial community profile of the translocated and woodland soils were similar; however, salvaged soils had a 40% increase in the Arbuscular Mycorrhizal Fungi (AMF) bioindicator fatty acid 16:1{omega}5, a 10% decrease in the Gram-positive bacterial fatty acids, and increased pH (5.01-5.77) compared to the original donor woodland soil. The AMF bioindicator and the first Principal Component (PC1) of the PCA of PLFA data positively correlated with soil pH (r2=0.94 and r2=0.88 respectively) across all three experimental groups. Considering that soil pH increases with depth in this location, it is likely that mixing of soil horizons during translocation increased the topsoil pH causing changes in the soil microbial communities. We concluded that after three years, the chemical and microbial properties of the salvaged soil were characteristic of a woodland soil but showed signs of disturbance.

Implications for practiceO_LITranslocated semi-ancient woodland soils may still exhibit similar microbiological characteristics to the donor site three years after being moved to a pasture receptor site.
C_LIO_LIHowever, disturbances expressed as significant changes in chemistry and microbial communities can affect soils following their translocation.
C_LIO_LIWhere soil pH varies with depth, translocation by loose-tipping may redistribute acidity within the soil profile, with important cascading impacts on the soil microbial communities.
C_LI


title: Does Communal Breeding Promote an Increase in Social Immunity in Burying Beetles? A Test Case with Nicrophorus Defodiens

abstract: Communal breeding is a reproductive system in which more than a single pair of individuals share parental care duties. Burying beetles (genus Nicrophorus) breed on small vertebrate carcasses, which is used as a food source for their young. On larger carcasses, burying beetles will breed communally, forming multiple male-female associations. A significant and costly component of parental investment by burying beetles is the preservation of the carcass with secretions containing immune molecules. Because this immune investment is for the benefit of the offspring, the behavior is a form of social immunity. We test the hypothesis that communal breeding in burying beetles evolved as a mechanism to increase the social immune investment on larger carcasses, which are more difficult to preserve. We used N. defodiens, a communal breeding burying beetle species to test the hypothesis. There were two experimental treatments wherein, the females either bred communally or non-communally. Our results show that the combined immune activity in the secretions were higher in communally breeding pairs than in the immune contribution of single male-female pairs. However, subordinate females were rarely observed on the carcass, and the level of social immune activity of dominant females was lower than females breeding singly. These data suggest that communal breeding in N. defodiens decreases the level of investment in social immunity. Our results demonstrate that the presence of multiple females, which is common under natural conditions, can greatly complicate patterns of social immunity investment in burying beetles.


title: Twisting to freedom: The evolution of copulation termination techniques across 48 species of sepsids (Diptera, Sepsidae)

abstract: Studies of insect mating behaviour usually focus on what happens before and during copulation. Few pay close attention to the actions needed to end copulation. However, genital separation after copulation is likely to be an important cause of mechanical stress and injuries because it often involves the withdrawal of heavily armoured male intromittent organs from membranous female reproductive tracts. Difficult and/or slow separations can also reduce male and female fitness by increasing their exposure to predation. We here report the results of a comparative study of separation behaviour in 48 species of Sepsidae (Diptera) and one outgroup. We find a surprising amount of qualitative and quantitative behavioural variability within and between species. We characterize and reconstruct three types of behaviours: 1) The sepsid ancestor likely used  back-off; a gentle separation technique that does not involve any pulling or twisting (https://youtu.be/EbkJvOaubZ0). 2) This potentially gave rise to the most common  pull technique where the male turns 180 degrees and pulls in an opposite direction from the female (https://youtu.be/oLf4xGpkk1s). This separation can be quick and straightforward, but in some species the  pull is slow and protracted and we routinely find dead males and/or females attached to their living partners in the latter (difficult: https://youtu.be/MbYPbXN6jr0; failure: https://youtu.be/leTiXefFzCc). 3) Finally, several species use  twist, a new technique where the male rotates >360 degrees from the initial mounting position (https://youtu.be/WMUXbIPyLbk). We document that species capable of using  twist have shorter and less variable separation times than those using "pull". However, many species capable of  twist also retain the ability to use  pull ( back-off/pull= 8%;  pull only= 41%;  twist/  pull= 24%;  twist only = 27%). Overall, our study suggests that separation behaviour can vary among closely related species and highlights the significance of studying variable behavioural traits in a phylogenetic context.

category: evolutionary biology


title: A supergene in seaweed flies modulates male traits and female perception

abstract: Supergenes, tightly linked allelic combinations that underlie complex adaptive phenotypes represent a critical mechanism protecting intra-specific polymorphism [1, 2]. Supergenes represent some of the best examples of balancing selection in nature and there is increasing evidence that disassortative mating, when individuals preferentially mate with dissimilar phenotypes, is a key force stabilizing supergene polymorphisms. Yet, the underlying biological mechanisms and genetic basis of disassortative mating remain poorly known. Here, we examine a possible mechanism of disassortative mating driven by female mate choice in relation to the overdominant Cf-Inv(1) supergene in the seaweed fly Coelopa frigida by investigating chemical communication and its genomic architecture. We show that Cf-Inv(1) strongly affects chemical signaling; cuticular hydrocarbon (CHC) composition differed between genotypes in males but not females across two continents. In tandem, Cf-Inv(1) affected female perception of these compounds; females are able to sense 36 compounds from the male CHC cocktail but show differential perception between genotypes for almost half of them. This indicates that the genetic underpinnings of male traits and female perceptions are tightly linked within Cf-Inv(1) which likely facilitates disassortative mating [3]. A differential expression approach based on candidate genes for CHC biosynthesis and odorant detection revealed differential expression for CHC biosynthesis in males alone but broad changes in odorant receptors across both sexes. Furthermore, odorant genes clustered together within Cf-Inv(1), with some of them differing between arrangements by 8.3% at the protein level, suggesting evolution via tandem duplication then divergence. We propose that the tight linkage between overdominant loci, male traits, and female perception has helped to maintain the Cf-Inv(1) polymorphism across its range in the face of supergene degeneration.


title: Fin Shape, Asymmetry, and Evolutionary Ecomorphology in Triggerfishes and Filefishes (Superfamily: Balistoidea)

abstract: Triggerfishes and filefishes exhibit a wide range of fin and body morphologies, inhabit many marine habitats, and feed on a variety of benthic and pelagic organisms. Particular morphologies are predicted to provide functional advantages for swimming behaviors that facilitate life in diverse habitats and feeding guilds. Ecomorphological relationships can, in turn, inform evolutionary patterns of morphological convergence. We quantified morphological diversity of 80 balistoid species using geometric morphometrics and assigned each species a primary habitat and feeding mode. Results revealed strong evidence for evolutionary integration among body and fin shapes as well as widespread convergence of both high and low aspect ratio (AR) dorsal and anal fins, the fins that power steady locomotion in these fishes. Dorsal and anal fins were determined to be moderately to highly asymmetrical in most species. Families exhibited considerable overlap in fin and body shapes, but triggerfishes generally exhibited higher AR and more asymmetrical fins than filefishes. Fin asymmetry was not strongly associated with ecology. Planktivorous and offshore-pelagic species exhibited high AR dorsal and anal fins suitable for high endurance swimming performance, while benthic grazing and structured reef species exhibited convergence on low AR median fins more suitable for facilitating maneuverability.


title: daf-16/FOXO blocks adult cell fate in Caenorhabditis elegans dauer larvae via lin-41/TRIM71

abstract: Many tissue-specific stem cells maintain the ability to produce multiple cell types during long periods of non-division, or quiescence. FOXO transcription factors promote quiescence and stem cell maintenance, but the mechanisms by which FOXO proteins promote multipotency during quiescence are still emerging. The single FOXO ortholog in C. elegans, daf-16, promotes entry into a quiescent and stress-resistant larval stage called dauer in response to adverse environmental cues. During dauer, stem and progenitor cells maintain or re-establish multipotency to allow normal development to resume after dauer. We find that during dauer, daf-16/FOXO prevents epidermal stem cells (seam cells) from prematurely adopting differentiated, adult characteristics. In particular, dauer larvae that lack daf-16 misexpress collagens that are normally adult-enriched. Using col-19p::gfp as an adult cell fate marker, we find that all major daf-16 isoforms contribute to opposing col-19p::gfp expression during dauer. By contrast, daf-16(0) larvae that undergo non-dauer development do not misexpress col-19p::gfp. Adult cell fate and the timing of col-19p::gfp expression are regulated by the heterochronic gene network, including lin-41 and lin-29. lin-41 encodes an RNA-binding protein orthologous to LIN41/TRIM71 in mammals, and lin-29 encodes a conserved zinc finger transcription factor. In non-dauer development lin-41 opposes adult cell fate by inhibiting the translation of lin-29, which directly activates col-19 transcription and promotes adult cell fate. We find that during dauer, lin-41 blocks col-19p::gfp expression, but surprisingly, lin-29 is not required in this context. Additionally, daf-16 promotes the expression of lin-41 in dauer larvae. The col-19p::gfp misexpression phenotype observed in dauer larvae with reduced daf-16 requires the downregulation of lin-41, but does not require lin-29. Taken together, this work demonstrates a novel role for daf-16/FOXO as a heterochronic gene that promotes expression of lin-41/TRIM71 to contribute to multipotent cell fate in a quiescent stem cell model.

category: genetics


title: The histone chaperone FACT facilitates heterochromatin spreading through regulation of histone turnover and H3K9 methylation states

abstract: Heterochromatin formation requires three distinct steps: nucleation, self-propagation (spreading) along the chromosome, and faithful maintenance after each replication cycle. Impeding any of those steps induces heterochromatin defects and improper gene expression. The essential histone chaperone FACT has been implicated in heterochromatin silencing, however, the mechanisms by which FACT engages in this process remain opaque. Here, we pin-pointed its function to the heterochromatin spreading process. FACT impairment reduces nucleation-distal H3K9me3 and HP1/Swi6 accumulation at subtelomeres and derepresses genes in the vicinity of heterochromatin boundaries. FACT promotes spreading by repressing heterochromatic histone turnover, which is crucial for the H3K9me2 to me3 transition that enables spreading. FACT mutant spreading defects are suppressed by removal of the H3K9 methylation antagonist Epe1 via nucleosome stabilization. Together, our study identifies FACT as a histone chaperone that specifically promotes heterochromatin spreading and lends support to the model that regulated histone turnover controls the propagation of epigenetic marks.


title: High-throughput variant detection using a color-mixing strategy in real-time PCR reactions

abstract: Many diseases are related to multiple genetic alterations along a single gene. Probing for highly multiple (>10) variants in a single qPCR tube is not possible due to a limited number of fluorescence channels and one variant per channel, so many more tubes are needed. Here, we experimentally validate our novel color-mixing strategy that uses fluorescence combinations as digital color codes to probe multiple variants simultaneously. The color-mixing strategy relies on a simple intra-tube assay that can probe for 15 variants as part of an inter-tube assay that can probe for an exponentially increased number of variants. The color-mixing strategy is achieved using multiplex double-stranded toehold probes modified with fluorophores and quenchers; the probes are designed to be quenched or luminous after binding to wildtype or variant templates. We used the color-mixing strategy to probe for 21 pathogenic mutations in thalassemia and to distinguish between heterozygous and homozygous variants in 6 tubes, with a specificity of 99% and a sensitivity of 94%. To support tuberculosis diagnosis, we used the same strategy to simultaneously probe in Mycobacterium tuberculosis for rifampicin-resistance mutations occurring within one 81-bp region and one 48-bp region in rpoB gene, plus five isoniazid-resistance mutations in inhA and katG genes.


title: High-throughput variant detection using a color-mixing strategy


title: A role for worm cutl-24 in background- and parent-of-origin-dependent ER stress resistance

abstract: Organisms in the wild can acquire disease- and stress-resistance traits that outstrip the programs endogenous to humans. Finding the molecular basis of such natural resistance characters is a key goal of evolutionary genetics. Standard statistical-genetic methods toward this end can perform poorly in organismal systems that lack high rates of meiotic recombination, like Caenorhabditis worms. Here we discovered unique ER stress resistance in a wild Kenyan C. elegans isolate, which in inter-strain crosses was passed by hermaphrodite mothers to hybrid offspring. We developed an unbiased version of the reciprocal hemizygosity test, RH-seq, to explore the genetics of this parent-of-origin-dependent phenotype. Among top-scoring gene candidates from a partial-coverage RH-seq screen, we focused on the neuronally-expressed, cuticlin-like gene cutl-24 for validation. In gene disruption and controlled crossing experiments, we found that cutl-24 was required in Kenyan hermaphrodite mothers for ER stress tolerance in their inter-strain hybrid offspring; cutl-24 was also a contributor to the trait in purebred backgrounds. These data establish the Kenyan strain allele of cutl-24 as a determinant of a natural stress-resistant state, and they set a precedent for the dissection of natural trait diversity in invertebrate animals without the need for a panel of meiotic recombinants.


title: Complex Mutation Profiles in Mismatch Repair and Ribonucleotide Reductase Mutations Reveal Novel Repair Substrate Specificity of MutS Homolog (MSH) Complexes

abstract: Determining mutation signatures is standard for understanding the etiology of human tumors and informing cancer treatment. Multiple determinants of DNA replication fidelity prevent mutagenesis that leads to carcinogenesis, including the regulation of free deoxyribonucleoside triphosphate (dNTP) pools by ribonucleotide reductase (RNR) and repair of replication errors by the mismatch repair (MMR) system. We identified genetic interactions between rnr1 alleles that elevate dNTP levels and MMR. We then utilized a targeted deep-sequencing approach to determine mutational signatures associated with MMR pathway defects. By combining rnr1 and msh mutations to increase dNTP levels and alter the mutational load, we uncovered previously unreported specificities of Msh2-Msh3 and Msh2-Msh6. Msh2-Msh3 is uniquely able to direct repair of G/C single base deletions in GC runs, while Msh2-Msh6 specifically directs repair of substitutions at G/C dinucleotides. We also identified broader sequence contexts that influence variant profiles in different genetic backgrounds. Finally, we observed that the mutation profiles in double mutants were not necessarily an additive relationship of mutation profiles in single mutants. Our results have implications for interpreting mutation signatures from human tumors, particularly when MMR is defective.


title: Identification of transcription termination defects at DNA hypomethylated transcription termination sites in DNA methyltransferase 3a-deficient vertebrates.

abstract: CpG methylation in genomic DNA is well known as a repressive epigenetic marker in eukaryotic transcription, and DNA methylation of the promoter regions is correlated with silencing of gene expression. In contrast to the promoter regions, the function of DNA methylation during transcription termination remains to be elucidated. A recent study has revealed that mouse DNA methyltransferase 3a (Dnmt3a) mainly functions in de novo methylation in the promoter and gene body regions (including transcription termination sites (TTSs)) during development. To investigate the relationship between DNA methylation overlapping the TTSs and transcription termination, we employed two strategies: informatic analysis using already deposited datasets of Dnmt3a-/- mouse cells and the zebrafish model system. Bioinformatic analysis using methylome and transcriptome data showed that hypomethylated differentially methylated regions overlapping the TTSs were associated with increased read counts and chimeric transcripts downstream of TTSs in Dnmt3a-/- Agouti-related protein neurons, but not in Dnmt3a-/- embryonic stem cells and mouse embryonic fibroblasts. We experimentally detected increased read-through and chimeric transcripts downstream of hypomethylated TTSs in zebrafish maternal-zygotic dnmt3aa-/- mutant. This study is the first to identify transcription termination defects in DNA hypomethylated TTSs in Dnmt3a-/- vertebrates.

category: genomics


title: Probing molecular specificity with deep sequencing and biophysically interpretable machine learning

abstract: Quantifying sequence-specific protein-ligand interactions is critical for understanding and exploiting numerous cellular processes, including gene regulation and signal transduction. Next-generation sequencing (NGS) based assays are increasingly being used to profile these interactions with high-throughput. However, these assays do not provide the biophysical parameters that have long been used to uncover the quantitative rules underlying sequence recognition. We developed a highly flexible machine learning framework, called ProBound, to define sequence recognition in terms of biophysical parameters based on NGS data. ProBound quantifies transcription factor (TF) behavior with models that accurately predict binding affinity over a range exceeding that of previous resources, captures the impact of DNA modifications and conformational flexibility of multi-TF complexes, and infers specificity directly from in vivo data such as ChIP-seq without peak calling. When coupled with a new assay called Kd-seq, it determines the absolute affinity of protein-ligand interactions. It can also profile the kinetics of kinase-substrate interactions. By constructing a biophysically robust foundation for profiling sequence recognition, ProBound opens up new avenues for decoding biological networks and rationally engineering protein-ligand interactions.


title: Genome of the parasitoid wasp Dinocampus coccinellae reveals extensive duplications, accelerated evolution, and independent origins of thelytokous parthenogeny and solitary behavior

abstract: Dinocampus coccinellae (Hymenoptera: Braconidae) is a generalist parasitoid wasp that parasitizes >50 species of predatory lady beetles (Coleoptera: Coccinellidae), with thelytokous parthenogeny as its primary mode of reproduction. Here we present the first high quality genome of D. coccinellae using a combination of short and long read sequencing technologies, followed by assembly and scaffolding of chromosomal segments using Chicago+ HiC technologies. We also present a first-pass ab initio genome annotation, and resolve timings of divergence and evolution of (1) solitary behavior vs eusociality, (2) arrhenotokous vs thelytokous parthenogenesis, and (3) rates of gene loss and gain among Hymenopteran lineages. Our study finds (1) at least two independent origins of eusociality and solitary behavior among Hymenoptera, (2) two independent origins of thelytokous parthenogenesis from ancestral arrhenotoky, and (3) accelerated rates of gene duplications, loss, and gain along the lineages leading to D. coccinellae. Our work both affirms the ancient divergence of Braconid wasps from ancestral Hymenopterans and accelerated rates of evolution in response to adaptations to novel hosts, including polyDNA viral co-evolution.


title: Evaluation of persistence and fate of ex vivo edited-HSC modified with donor template and Its role in correcting Sickle Cell Disease

abstract: Sickle cell disease (SCD) is caused by a single nucleotide transversion in exon 1 of the HBB gene that changes the hydrophobicity of adult globin ({beta}A), leading to substantial morbidity and reduced lifespan. Ex vivo autologous gene editing utilizing co-delivery of a designer nuclease along with a DNA donor template allows for precise homology-directed repair (HDR). These gene corrected cells when engrafted into the bone marrow (BM) can prove to be therapeutic and serves as an alternative to HLA-matched BM transplantation. In the current study, we extensively explored the role of single stranded oligonucleotide (ssODN) and recombinant adeno-associated 6 (rAAV6) donor template delivery to introduce a codon-optimized change (E6optE) or a sickle mutation (E6V) change following Crispr/Cas9-mediated cleavage of HBB in healthy human mobilized peripheral blood stem cells (mPBSCs). We achieved efficient HDR in vitro in edited cells and observed robust human CD45+ engraftment in the BM of NBSGW mice at 16-17 weeks. Notably, recipients of ssODN-modified HSC exhibited a significantly higher proportion of HDR-modified cells within individual BM, CD34+ and CD235+ compartments of both E6optE and E6V cohorts. We further assessed key functional outcomes including RNA transcripts analysis and globin sub-type expression. Our combined findings demonstrate the capacity to achieve clinically relevant HDR in vitro and in vivo using both donor template delivery method. The use of ssODN donor template-delivery is consistently associated with higher levels of gene correction in vivo as demonstrated by sustained engraftment of HDR-modified HSC and erythroid progeny. Finally, the HDR-based globin protein expression was significantly higher in the E6V ssODN-modified animals compared to the rAAV6-modified animals confirming that the ssODN donor template delivery outperforms rAAV6-donor template delivery.


title: Malignancy and NF-kB signalling strengthen coordination between the expression of mitochondrial and nuclear-encoded oxidative phosphorylation genes

abstract: Mitochondria are ancient endosymbiotic organelles crucial to eukaryotic growth and metabolism. Mammalian mitochondria carry a small genome containing thirteen protein-coding genes with the remaining mitochondrial proteins encoded by the nuclear genome. Little is known about how coordination between the two sets of genes is achieved. Correlation analysis of RNA-seq expression data from large publicly-available datasets is a common method to leverage genetic diversity to infer gene co-expression modules. Here we use this method to investigate nuclear-mitochondrial gene expression coordination. We identify a pitfall in correlation analysis that results from the large variation in the proportion of transcripts from the mitochondrial genome in RNA-seq data. Commonly used normalization techniques based on total read count (such as FPKM or TPM) produce artefactual negative correlations between mitochondrial- and nuclear-encoded transcripts. This also results in artefactual correlations between pairs of nuclear-encoded genes, thus having important consequences for inferring co-expression modules beyond mitochondria. We show that these effects can be overcome by normalizing using the median-ratio normalization (MRN) or trimmed mean of M values (TMM) methods. Using these normalizations, we find only weak and inconsistent correlations between mitochondrial and nuclear-encoded mitochondrial genes in the majority of healthy human tissues from the GTEx database. However, a subset of healthy tissues with high expression of NF-{kappa}B show significant coordination supporting a role for NF-{kappa}B in retrograde signalling. Contrastingly, most cancer types show robust coordination of nuclear and mitochondrial OXPHOS gene expression, identifying this as a feature of gene regulation in cancer.


title: Network-Based Discovery of Opioid Use Vulnerability in Rats Using the Bayesian Stochastic Block Model

abstract: Opioid use disorder is a psychological condition that affects over 200,000 people per year in the U.S., causing the Centers for Disease Control and Prevention to label the crisis as a rapidly spreading public health epidemic. It has been found that the behavioral relationship between opioid exposure and development of opioid use disorder varies greatly between individuals, implying existence of sup-populations with varying degrees of opioid vulnerability. In this study, we assessed several behavioral variables across heroin taking, refraining and seeking to establish how these factors interact with one another resulting in a heroin dependent, resilient, or vulnerable behavioral phenotype. Over 400 (male and female) heterogeneous stock rats were used in these two studies, and data were collected from two geographically distinct locations. Rats underwent heroin self-administration training, followed by a progressive ratio and heroin-primed reinstatement test. Next, rats underwent extinction training and a cue-induced reinstatement test. To assess how these variables contribute to heroin addiction vulnerability, we developed a network-based data analysis workflow. Specifically, we integrated different cohorts of rats, remove possible batch effects, and constructed a rat-rat similarity network based on their behavioral patterns. We then implemented community detection on this similarity network using a Bayesian degree-corrected stochastic block model to uncover sub-populations of rats with differing levels of opioid vulnerability. We discovered three distinct behavioral sub-populations, each with significantly different behavioral outcomes that allowed for unique characterization of each cluster in terms of vulnerability to opioid use and seeking. We implement this analysis workflow as an open source R package, named mlsbm.


title: In Situ Volumetric Imaging and Analysis of FRESH 3D Bioprinted Constructs Using Optical Coherence Tomography

abstract:
As 3D bioprinting has grown as a fabrication technology, so too has the need for improved analytical methods to characterize these engineered constructs. This is especially challenging for soft tissues composed of hydrogels and cells as these materials readily deform, posing a barrier when trying to assess print fidelity and other properties non- destructively. Indeed, given the importance of structure-function relationships in most tissue types, establishing that the 3D architecture of the bioprinted construct matches its intended anatomic design is critical. Here we report development of a multimaterial bioprinting platform with integrated optical coherence tomography (OCT) for in situ volumetric imaging, error detection, and 3D reconstruction. While generally applicable to extrusion-based 3D bioprinting, we also report improvements to the Freeform Reversible Embedding of Suspended Hydrogels (FRESH) bioprinting process through new collagen bioink compositions, support bath optical clearing, and machine pathing. This enables high-fidelity 3D volumetric imaging with micron scale resolution over centimeter length scales, the ability to detect a range of print defect types within a 3D volume, and real-time imaging of the printing process at each print layer. These advances provide FRESH and other extrusion-based 3D bioprinting approaches with a comprehensive methodology for quality assessment that has been absent in the field to date, paving the way for translation of these engineered tissues to the clinic and ultimately to achieving regulatory approval.

TeaserTransparent FRESH support bath enables in situ 3D volumetric imaging and validation of patient-derived tissue constructs.


title: Monitoring the effects of chemical stimuli on live cells with metasurface-enhanced infrared reflection spectroscopy

abstract: Infrared spectroscopy has found wide applications in the analysis of biological materials. A more recent development is the use of engineered nanostructures - plasmonic metasurfaces - as substrates for metasurface-enhanced infrared reflection spectroscopy (MEIRS). Here, we demonstrate that strong field enhancement from plasmonic metasurfaces enables the use of MEIRS as a highly informative analytic technique for real-time monitoring of cells. By exposing live cells cultured on a plasmonic metasurface to chemical compounds, we show that MEIRS can be used as a label-free phenotypic assay for detecting multiple cellular responses to external stimuli: changes in cell morphology, adhesion, lipid composition of the cellular membrane, as well as intracellular signaling. Using a focal plane array detection system, we show that MEIRS also enables spectro-chemical imaging at the single-cell level. The described metasurface-based all-optical sensor opens the way to a scalable, high-throughput spectroscopic assay for live cells.


title: Shortening of 3' UTRs in most cell types composing tumor tissues implicates alternative polyadenylation in protein metabolism

abstract:
During pre-mRNA maturation 3 end processing can occur at different polyadenylation sites in the 3 untranslated region (3 UTR) to give rise to transcript isoforms that differ in the length of their 3UTRs. Longer 3 UTRs contain additional cis-regulatory elements that impact the fate of the transcript and/or of the resulting protein.

Extensive alternative polyadenylation (APA) has been observed in cancers, but the mechanisms and roles remain elusive. In particular, it is unclear whether the APA occurs in the malignant cells or in other cell types that infiltrate the tumor. To resolve this, we developed a computational method, called SCUREL, that quantifies changes in 3UTR length between groups of cells, including cells of the same type originating from tumor and control tissue. We used this method to study APA in human lung adenocarcinoma (LUAD).

SCUREL relies solely on annotated 3UTRs and on control systems, such as T cell activation and spermatogenesis gives qualitatively similar results at much greater sensitivity compared to the previously published scAPA method.

In the LUAD samples, we find a general trend towards 3UTR shortening not only in cancer cells compared to the cell type of origin, but also when comparing other cell types from the tumor vs. the control tissue environment. However, we also find high variability in the individual targets between patients. The findings help to understand the extent and impact of APA in LUAD, which may support improvements in diagnosis and treatment.


title: Optimal balancing of clinical factors in large scale clinical RNA-Seq studies

abstract: Omics technologies are ubiquitous in biomedical research. However, improper sample selection is an often-overlooked complication with large omics studies, resulting in confounding effects that can disrupt the internal validity of a study and lead to false conclusions. Here, we present a method called BalanceIT, which uses a genetic algorithm to identify an optimal set of samples with balanced clinical factors for large-scale omics experiments. We apply our approach to two large RNA-Seq studies in autism (1) to find a post-hoc balanced sample set among an imbalanced study, and (2) to design an optimal study that allows for efficient batch correction. Our approach leads to near-perfect estimates of differential gene expression, superior performance of pathway-level enrichment analysis, and consistent network dysregulation patterns of autism symptom severity. These results provide empirical support for the importance of balanced experimental design, and BalanceIT will be invaluable for large-scale study design and batch effect correction.


title: Benchmarking supervised signature-scoring methods for single-cell RNA sequencing data in cancer

abstract: Quantifying the activity of gene expression signatures is common in analyses of single-cell RNA sequencing data. Methods originally developed for bulk samples are often used for this purpose without accounting for contextual differences between bulk and single-cell data. More broadly, these methods have not been benchmarked. Here we benchmark four such supervised methods, including single sample gene set enrichment analysis (ssGSEA), AUCell, Single Cell Signature Explorer (SCSE), and a new method we developed, Jointly Assessing Signature Mean and Inferring Enrichment (JASMINE). Using cancer as an example, we show cancer cells consistently express more genes than normal cells. This imbalance leads to bias in performance by bulk-sample-based ssGSEA in gold standard tests and down sampling experiments. In contrast, single-cell-based methods are less susceptible. Our results suggest caution should be exercised when using bulk-sample-based methods in single-cell data analyses, and cellular contexts should be taken into consideration when designing benchmarking strategies.


title: The flying insect thoracic cuticle is heterogenous in structure and in thickness-dependent modulus gradation

abstract:
The thorax is a specialized structure central to an insects ability to fly. In the thorax, flight muscles are surrounded by a thin layer of cuticle. The structure, composition, and material properties of this chitinous structure may influence the efficiency of the thorax in flight. However, these properties, as well as their variation throughout anatomical regions of the thorax or between insect taxa, are not known. In this work, we provide a multi-faceted assessment of thorax cuticle for fliers with asynchronous (honey bee; Apis mellifera) and synchronous (hawkmoth; Manduca sexta) muscles. We investigated cuticle structure using histology, material composition through confocal laser scanning microscopy, and modulus gradation with nanoindentation. Our results suggest that cuticle properties of the thorax are highly dependent on anatomical region and species. Modulus gradation, but not mean modulus, differed between the two types of fliers. In some regions, A. mellifera had a positive linear modulus gradient from cuticle interior to exterior of about 2 GPa. In M. sexta, the modulus gradients were variable and were not well represented by linear fits with respect to cuticle thickness. We utilized finite element modeling to assess how measured modulus gradients influenced maximum stress in cuticle. Stress was reduced when cuticle with a linear gradient was compressed from the high modulus side. These results support the protective role of the A. mellifera thorax cuticle. Our multi-faceted assessment advances our understanding of thorax cuticle structural and material heterogeneity and the potential benefit of material gradation to flying insects.

Graphical Abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

Statement of SignificanceThe insect thorax is essential for efficient flight but questions remain about the contribution of exoskeletal cuticle. We assessed the thorax cuticle using a high resolution multi-faceted approach to determine how cuticle properties vary within thorax anatomical regions and between fliers with asynchronous (honey bee; Apis mellifera) and synchronous (hawkmoth; Manduca sexta) muscles. We examined structure using histological staining, modulus using nanoindentation, and material composition using confocal scanning light microscopy. We further utilized finite element modeling to understand the effect of the modulus gradations observed experimentally on stress accumulation. Cuticle properties vary through cuticle thickness, by thorax region, and between flight lineages.


title: Proximity-dependent biotinylation to elucidate the interactome of TNK2 non-receptor tyrosine kinase

abstract: Non-receptor tyrosine kinases represent an important class of signaling molecules which are involved in driving diverse cellular pathways. Although, the large majority have been well-studied in terms of their protein-binding partners, the interactomes of some of the key non-receptor tyrosine kinases such as TNK2 (also known as activated Cdc42-associated kinase 1 or ACK1) have not been systematically investigated. Aberrant expression and hyperphosphorylation of TNK2 has been implicated in a number of cancers. However, the exact proteins and cellular events that mediate phenotypic changes downstream of TNK2 are unclear. Biological systems that employ proximity-dependent biotinylation methods, such as BioID, are being increasingly used to map protein-protein interactions as they provide increased sensitivity in discovering interaction partners. In this study, we employed BioID coupled to the biotinylation site identification technology (BioSITe) method that we recently developed to perform molecular mapping of intracellular proteins associated with TNK2. We also employed stable isotope labeling with amino acids in cell culture (SILAC) to quantitatively explore the interactome of TNK2. By performing a controlled comparative analysis between full-length TNK2 and its truncated counterpart, we were not only able to confidently identify site-level biotinylation of previously well-established TNK2 binders and substrates such as NCK1, NCK2, CTTN, STAT3, but also discover several novel TNK2 interacting partners. We validated TNK2 interaction with one of the novel TNK2 interacting protein, clathrin interactor 1 (CLINT1), using immunoblot analysis. Overall, this work reveals the power of the BioSITe method coupled to BioID and highlights several molecules that warrant further exploration to assess their functional significance in TNK2-mediated signaling.


title: Small Vessel Diseases: 3D Characteristics of the Vasculature and White Matter

abstract: Cerebral small vessel disease (SVD) is associated with white matter hyperintensities (WMHs), thereby contributing to vascular dementia and movement disorder. However, the pathomechanisms responsible for WMH-related small vessel degeneration remain poorly understood due to the technical limitations of current methods. The aim of this study was to clarify the 2- and 3-dimensional (2D and 3D) pathological features of small vessels and white matter (WM) in the brains of patients with cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL), HTRA1-autosomal dominant disease (HTRA1-AD) and sporadic SVD (sSVD). From a cohort of 86 consecutive autopsied patients with SVDs, we retrieved those with genetically confirmed CADASIL and HTRA1-AD (three and three, respectively), and four with sSVD. We quantitatively evaluated WM and vascular changes in the frontal portion of the centrum semiovale and temporal lobe using conventional 2D and chemically cleared 3D analytical methods with light-sheet fluorescence microscopy. Quantitatively, the WM pathology, including the density of myelin, axons and gliosis, was most severe in CADASIL, but unexpectedly sSVD was second in order of severity, followed by HTRA1-AD. The density of clasmatodendrocytes, known to be irreversibly injured astrocytes, was considerably highest in HTRA1-AD. The vascular pathology, including arteriole and capillary sclerosis and the extent of the perivascular space, was most severe in CADASIL, whereas the density of smooth muscle actin (SMA) positivity was most decreased in HTRA1-AD. 3D immunohistochemistry for SMA demonstrated two distinct patterns of SMA loss within the vessels: (1) CADASIL and sSVD: diffuse loss, being prominent in small branches, (2) HTRA1-AD: selective loss in main branches. Overall, the extent of WM and vascular degeneration is most severe in CADASIL, whereas SMA loss is most evident in HTRA1-AD. These differences in the size and distribution of affected vessels may be related to the heterogeneous WM pathology and underlying pathomechanisms of SVD.

category: pathology


title: Structural studies on cell wall biogenesis-I: A method to produce high-quality protoplasts from alga Micrasterias denticulata, an emerging model in plant biology

abstract: Micrasterias denticulate is a freshwater unicellular green alga emerging as a model system in plant cell biology. This is an algae that has been examined in the context of cell wall research from early 1970s. Protoplast production from such a model system is important for many downstream physiological and cell biological studies. The algae produce intact protoplast in a straight two-step protocol involving 5% mannitol, 2% cellulysin, 4mM calcium chloride under a  temperature ramping strategy. The process of protoplast induction and behavior of protoplast was examined by light microscopy and reported in this study.


title: Mechanistic model-driven exometabolomic characterisation of human dopaminergic neuronal metabolism

abstract: Patient-derived cellular models are a powerful approach to study human disease, especially neurode-generative diseases, such as Parkinsons disease, where affected primary neurons, e.g., substantia nigra dopaminergic neurons, are almost inaccessible. Starting with a comprehensive generic reconstruction of human metabolism, Recon3D, we generated a high-quality, constraint-based, genome-scale, in silico model of human dopaminergic neuronal metabolism (iDopaNeuro1). It is a synthesis of extensive manual curation of the biochemical literature on neuronal metabolism, together with novel, quantitative, transcriptomic and targeted exometabolomic data from human stem cell-derived, midbrain-specific, dopaminergic neurons in vitro. Thermodynamic constraint-based modelling with iDopaNeuro1 is qualitatively accurate (92% correct) and quantitatively accurate (Spearman rank 0.7) at predicting metabolite secretion or uptake, given quantitative exometabolomic constraints on uptakes, or secretions, respectively. iDopaNeuro1 is also qualitatively accurate at predicting the consequences of metabolic perturbations, e.g., complex I inhibition (Spearman rank 0.69) in a manner consistent with literature on monogenic mitochondrial Parkinsons disease. The iDopaNeuro1 model provides a foundation for a quantitative systems biochemistry approach to metabolic dysfunction in Parkinsons disease. Moreover, the plethora of novel mathematical and computational approaches required to develop it are generalisable to study any other disease associated with metabolic dysfunction.

category: systems biology