Three-dimensional Epigenome Statistical Model: Genome-wide Chromatin Looping Prediction

Machine and Deep Learning Methods for Predicting 3D Genome Organization

Three-dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, topologically associating domains (TADs), and A/B compartments, play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers and transcription factor binding site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, and TAD boundaries) and analyze their pros and cons. We also point out obstacles to the computational prediction of 3D interactions and suggest future research directions.

A systematic analyses of different bioinformatics pipelines for genomic data and its impact on deep learning models for chromatin loop prediction

Genomic data analysis has witnessed a surge in complexity and volume, primarily driven by the advent of high-throughput technologies. In particular, studying chromatin loops and structures has become pivotal in understanding gene regulation and genome organization. This systematic investigation explores the realm of specialized bioinformatics pipelines designed specifically for the analysis of chromatin loops and structures. Our investigation incorporates two protein (CTCF and Cohesin) factor-specific loop interaction datasets from six distinct pipelines, amassing a comprehensive collection of 36 diverse datasets. Through a meticulous review of existing literature, we offer a holistic perspective on the methodologies, tools and algorithms underpinning the analysis of this multifaceted genomic feature. We illuminate the vast array of approaches deployed, encompassing pivotal aspects such as data preparation pipeline, preprocessing, statistical features and modelling techniques. Beyond this, we rigorously assess the strengths and limitations inherent in these bioinformatics pipelines, shedding light on the interplay between data quality and the performance of deep learning models, ultimately advancing our comprehension of genomic intricacies.

Developmental dysplasia of the hip caused by homozygous TRIM33 pathogenic variant affecting downstream BMP pathway

BACKGROUND

Developmental dysplasia of the hip (DDH), formerly termed congenital dislocation of the hip, is the most common congenital disease of the musculoskeletal system in newborns. While familial predilection to DDH has been well documented, the molecular genetics/pathways of this common disorder are poorly understood.

METHODS

Linkage analysis and whole exome sequencing; real-time PCR studies of skin fibroblasts.

RESULTS

Consanguineous Bedouin kindred presented with DDH with apparent autosomal recessive heredity. Linkage analysis and whole exome sequencing delineated a single 3.2 Mbp disease-associated chromosome 1 locus (maximal multipoint Logarithm of the Odds score 2.3), containing a single homozygous variant with a relevant expression pattern: addition of threonine in TRIM33 (NM_015906.4); c.1648_1650dup. TRIM33 encodes a protein that acts both in the TGF-β and the BMP pathways; however, it has been mostly studied regarding its function in the TGF-β pathway. Since BMPs are known to act in bone formation, we focused on the BMP pathway, in which TRIM33 functions as a transcription factor, both an activator and repressor. Skin fibroblasts of two affected girls and a heterozygous TRIM33 variant carrier were assayed through reverse-transcription PCR for expression of genes known to be downstream of TRIM33 in the BMP pathway: fibroblasts of affected individuals showed significantly reduced expression of DLX5, significantly increased expression of BGLAP, increased expression of ALPL and no change in expression of RUNX2 or of TRIM33 itself.

CONCLUSIONS

DDH can be caused by a biallelic variant in TRIM33, affecting the BMP pathway.

ELUCIDATING THE DIFFERENTIAL IMPACTS OF EQUIVALENT GATING-CHARGE MUTATIONS IN VOLTAGE-GATED SODIUM CHANNELS

Voltage-gated sodium (Nav) channels are pivotal for cellular signaling and mutations in Nav channels can lead to excitability disorders in cardiac, muscular, and neural tissues. A major cluster of pathological mutations localizes in the voltage-sensing domains (VSDs), resulting in either gain-of-function (GoF), loss-of-function (LoF) effects, or both. However, the mechanism behind this functional divergence of mutations at equivalent positions remains elusive. Through hotspot analysis, we identified three gating charges (R1, R2, and R3) as major mutational hotspots in VSDs. The same amino-acid substitutions at equivalent gating-charge positions in VSDI and VSDII of the cardiac sodium channel Nav1.5 show differential gating-property impacts in electrophysiology measurements. We conducted 120 μs molecular dynamics (MD) simulations on wild-type and six mutants to elucidate the structural basis of their differential impacts. Our μs-scale MD simulations with applied external electric fields captured VSD state transitions and revealed the differential structural dynamics between equivalent R-to-Q mutants. Notably, we observed transient leaky conformations in some mutants during structural transitions, offering a detailed structural explanation for gating-pore currents. Our salt-bridge network analysis uncovered VSD-specific and state-dependent interactions among gating charges, countercharges, and lipids. This detailed analysis elucidated how mutations disrupt critical electrostatic interactions, thereby altering VSD permeability and modulating gating properties. By demonstrating the crucial importance of considering the specific structural context of each mutation, our study represents a significant leap forward in understanding structure-function relationships in Nav channels. Our work establishes a robust framework for future investigations into the molecular basis of ion channel-related disorders.

Satellite glial contact enhances differentiation and maturation of human iPSC-derived sensory neurons

Sensory neurons generated from induced pluripotent stem cells (iSNs) are used to model human peripheral neuropathies, however current differentiation protocols produce sensory neurons with an embryonic phenotype. Peripheral glial cells contact sensory neurons early in development and contribute to formation of the canonical pseudounipolar morphology, but these signals are not encompassed in current iSN differentiation protocols. Here, we show that terminal differentiation of iSNs in co-culture with rodent Dorsal Root Ganglion satellite glia (rSG) advances their differentiation and maturation. Co-cultured iSNs develop a pseudounipolar morphology through contact with rSGs. This transition depends on semaphorin-plexin guidance cues and on glial gap junction signaling. In addition to morphological changes, iSNs terminally differentiated in co-culture exhibit enhanced spontaneous action potential firing, more mature gene expression, and increased susceptibility to paclitaxel induced axonal degeneration. Thus, iSNs differentiated in coculture with rSGs provide a better model for investigating human peripheral neuropathies.

Single cell autofluorescence imaging reveals immediate metabolic shifts of neutrophils with activation across biological systems

Neutrophils, the most abundant leukocytes in human peripheral circulation, are crucial for the innate immune response. They are typically quiescent but rapidly activate in response to infection and inflammation, performing diverse functions such as oxidative burst, phagocytosis, and NETosis, which require significant metabolic adaptation. Deeper insights into such metabolic changes will help identify regulation of neutrophil functions in health and diseases. Due to their short lifespan and associated technical challenges, the metabolic processes of neutrophils are not completely understood. This study uses optical metabolic imaging (OMI), which entails optical redox ratio and fluorescence lifetime imaging microscopy of intrinsic metabolic coenzymes NAD(P)H and FAD to assess the metabolic state of single neutrophils. Primary human neutrophils were imaged in vitro under a variety of activation conditions and metabolic pathway inhibitors, while metabolic and functional changes were confirmed with mass spectrometry, oxidative burst, and NETosis measurements. Our findings show that neutrophils undergo rapid metabolic remodeling to a reduced redox state, followed by a shift to an oxidized redox state during activation. Additionally, single cell analysis reveals a heterogeneous metabolic response across neutrophils and donors to live pathogen infection ( Pseudomonas aeruginosa and Toxoplasma gondii ). Finally, consistent metabolic changes were validated with neutrophils in vivo using zebrafish larvae. This study demonstrates the potential of OMI as a versatile tool for studying neutrophil metabolism and underscores its use across different biological systems, offering insights into neutrophil metabolic activity and function at a single cell level.

Sleep induced by mechanosensory stimulation provides cognitive and health benefits in Drosophila

Study Objectives

Sleep is a complex phenomenon regulated by various factors, including sensory input. Anecdotal observations have suggested that gentle rocking helps babies fall asleep, and experimental studies have verified that rocking promotes sleep in both humans and mice. Recent studies have expanded this understanding, demonstrating that gentle vibration also induces sleep in Drosophila. Natural sleep serves multiple functions, including learning and memory, synaptic downscaling, and clearance of harmful substances associated with neurodegenerative diseases. Here, we investigated whether vibration-induced sleep provides similar cognitive and health benefits in Drosophila.

Methods

We administered gentle vibration to flies that slept very little due to a forced activation of wake-promoting neurons and investigated how the vibration influenced learning and memory in the courtship conditioning paradigm. Additionally, we examined the effects of VIS on synaptic downscaling by counting synapse numbers of select neurons. Finally, we determined whether vibration could induce sleep in Drosophila models of Alzheimer's disease (AD) and promote the clearance of Amyloid b (Ab) and Tubulin Associated Unit (TAU).

Results

Vibration-induced sleep enhanced performance in a courtship conditioning paradigm and reduced the number of synapses in select neurons. Moreover, vibration improved sleep in Drosophila models of AD, promoting the clearance of Ab and TAU.

Conclusions

Mechanosensory stimulation offers a promising non-invasive avenue for enhancing sleep, potentially providing associated cognitive and health benefits.

Microglia Morphological Response to Mesenchymal Stromal Cell Extracellular Vesicles Demonstrates EV Therapeutic Potential for Modulating Neuroinflammation

Background

Mesenchymal stromal cell derived extracellular vesicles (MSC-EVs) are a promising therapeutic for neuroinflammation. MSC-EVs can interact with microglia, the resident immune cells of the brain, to exert their immunomodulatory effects. In response to inflammatory cues, such as cytokines, microglia undergo phenotypic changes indicative of their function e.g. morphology and secretion. However, these changes in response to MSC-EVs are not well understood. Additionally, no disease-relevant screening tools to assess MSC-EV bioactivity exist, which has further impeded clinical translation. Here, we developed a quantitative, high throughput morphological profiling approach to assess the response of microglia to neuroinflammation-relevant signals and whether this morphological response can be used to indicate the bioactivity of MSC-EVs.

Results

Using an immortalized human microglia cell-line, we observed increased size (perimeter, major axis length) and complexity (form factor) upon stimulation with interferon-gamma (IFN-γ) and tumor necrosis factor-alpha (TNF-α). Upon treatment with MSC-EVs, the overall morphological score (determined using principal component analysis) shifted towards the unstimulated morphology, indicating that MSC-EVs are bioactive and modulate microglia. The morphological effects of MSC-EVs in TNF-γ/IFN-α stimulated cells were concomitant with reduced secretion of 14 chemokines/cytokines (e.g. CXCL6, CXCL9) and increased secretion of 12 chemokines/cytokines (e.g. CXCL8, CXCL10). Proteomic analysis of cell lysates revealed significant increases in 192 proteins (e.g. HIBADH, MEAK7, LAMC1) and decreases in 257 proteins (e.g. PTEN, TOM1, MFF) with MSC-EV treatment. Of note, many of these proteins are involved in regulation of cell morphology and migration. Gene Set Variation Analysis revealed upregulation of pathways associated with immune response, such as regulation of cytokine production, immune cell infiltration (e.g. T cells, NK cells) and morphological changes (e.g. Semaphorin, RHO/Rac signaling). Additionally, changes in microglia mitochondrial morphology were measured suggesting that MSC-EV modulate mitochondrial metabolism.

Conclusion

This study comprehensively demonstrates the effects of MSC-EVs on human microglial morphology, cytokine secretion, cellular proteome, and mitochondrial content. Our high-throughput, rapid, low-cost morphological approach enables screening of MSC-EV batches and manufacturing conditions to enhance EV function and mitigate EV functional heterogeneity in a disease relevant manner. This approach is highly generalizable and can be further adapted and refined based on selection of the disease-relevant signal, target cell, and therapeutic product.

The Surgical Method of Craniectomy Differentially Affects Acute Seizures, Brain Deformation and Behavior in a TBI Animal Model

Traumatic brain injury (TBI) is the leading cause of morbidity and mortality worldwide. Multiple injury models have been developed to study this neurological disorder. One such model is the lateral fluid-percussion injury (LFPI) rodent model. The LFPI model can be generated with different surgical procedures that could affect the injury and be reflected in neurobehavioral dysfunction and acute EEG changes. A craniectomy was performed either with a trephine hand drill or with a trephine electric drill that was centered over the left hemisphere of adult, male Sprague Dawley rats. Sham craniectomy groups were assessed by hand-drilled (ShamHMRI) and electric-drilled (ShamEMRI) to evaluate by MRI. Then, TBI was induced in separate groups (TBIH) and (TBIE) using a fluid-percussion device. Sham-injured rats (ShamH/ShamE) underwent the same surgical procedures as the TBI rats. During the same surgery session, rats were implanted with screw and microwire electrodes positioned in the neocortex and hippocampus and the EEG activity was recorded 24 hours for the first 7 days after TBI for assessing the acute EEG seizure and Gamma Event Coupling (GEC). The electric drilling craniectomy induced greater tissue damage and sensorimotor deficits compared to the hand drill. Analysis of the EEG revealed acute seizures in at least one animal from each group after the procedure. Both TBI and Sham rats from the electric drill groups had a significant greater total number of seizures than the animals that were craniectomized manually (p<0.05). Similarly, EEG functional connectivity was lower in ShamE compared to ShamH rats. These results suggest that electrical versus hand drilling craniectomies produce cortical injury in addition to the LFPI which increases the likelihood for acute post-traumatic seizures. Differences in the surgical approach could be one reason for the variability in the injury that makes it difficult to replicate results between preclinical TBI studies.

Current Status of Clinical Trials Design and Outcomes in Retinal Gene Therapy

With the rapid expansion of methods encompassed by the term gene therapy, new trials exploring the safety and efficacy of these methods are initiated more frequently. As a result, important questions arise pertaining the design of these trials and patient participation. One of the most important aspects of any clinical trial is the ability to measure the trial's outcome in a manner that will reflect the effect of the treatment and allow its quantification, whether the trial is aimed at preservation or restoration of retinal cells (photoreceptors and others), vision, or both. Here we will review the existing methods for quantification of trial outcomes, stressing the importance of assessing the participant's visual function and not just visual acuity. We will also describe the key considerations in trial design. Finally, as patient safety remains the primary concern in any trial participation, we will outline the key principles in that regard.

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Sitetack: A Deep Learning Model that Improves PTM Prediction by Using Known PTMs

Post-translational modifications (PTMs) increase the diversity of the proteome and are vital to organismal life and therapeutic strategies. Deep learning has been used to predict PTM locations. Still, limitations in datasets and their analyses compromise success. Here we evaluate the use of known PTM sites in prediction via sequence-based deep learning algorithms. Specifically, PTM locations were encoded as a separate amino acid before sequences were encoded via word embedding and passed into a convolutional neural network that predicts the probability of a modification at a given site. Without labeling known PTMs, our model is on par with others. With labeling, however, we improved significantly upon extant models. Moreover, knowing PTM locations can increase the predictability of a different PTM. Our findings highlight the importance of PTMs for the installation of additional PTMs. We anticipate that including known PTM locations will enhance the performance of other proteomic machine learning algorithms.

Neuroanatomical Predictors of Transcranial Direct Current Stimulation (tDCS)-Induced Modifications in Neurocognitive Task Performance in Typically Developing Individuals

Transcranial direct current stimulation (tDCS) is a noninvasive neuromodulation technique gaining more attention in neurodevelopmental disorders (NDDs). Due to the phenotypic heterogeneity of NDDs, tDCS is unlikely to be equally effective in all individuals. The present study aimed to establish neuroanatomical markers in typically developing (TD) individuals that may be used for the prediction of individual responses to tDCS. Fifty-seven male and female children received 2 mA anodal and sham tDCS, targeting the left dorsolateral prefrontal cortex (DLPFCleft), right inferior frontal gyrus, and bilateral temporoparietal junction. Response to tDCS was assessed based on task performance differences between anodal and sham tDCS in different neurocognitive tasks (N-back, flanker, Mooney faces detection, attentional emotional recognition task). Measures of cortical thickness (CT) and surface area (SA) were derived from 3 Tesla structural MRI scans. Associations between neuroanatomy and task performance were assessed using general linear models (GLM). Machine learning (ML) algorithms were employed to predict responses to tDCS. Vertex-wise estimates of SA were more closely linked to differences in task performance than measures of CT. Across ML algorithms, highest accuracies were observed for the prediction of N-back task performance differences following stimulation of the DLPFCleft, where 65% of behavioral variance was explained by variability in SA. Lower accuracies were observed for all other tasks and stimulated regions. This suggests that it may be possible to predict individual responses to tDCS for some behavioral measures and target regions. In the future, these models might be extended to predict treatment outcome in individuals with NDDs.

Molecular phenotypes segregate missense mutations in SLC13A5 Epilepsy

The sodium-coupled citrate transporter (NaCT, SLC13A5) mediates citrate uptake across the plasma membrane via an inward Na + gradient. Mutations in SLC13A5 cause early infantile epileptic encephalopathy type-25 (EIEE25, SLC13A5 Epilepsy) due to impaired citrate uptake in neurons. Despite clinical identification of disease-causing mutations, underlying mechanisms and cures remain elusive. We mechanistically classify the molecular phenotypes of six mutations. C50R, T142M, and T227M exhibit impaired citrate transport despite normal expression at the cell surface. G219R, S427L, and L488P are hampered by low protein expression, ER retention, and reduced transport. Mutants' mRNA levels resemble wildtype, suggesting post-translational defects. Class II mutations display immature core-glycosylation and shortened half-lives, indicating protein folding defects. These experiments provide a comprehensive understanding of the mutation's defects in SLC13A5 Epilepsy at the biochemical and molecular level and shed light into the trafficking pathway(s) of NaCT. The two classes of mutations will require fundamentally different treatment approaches to either restore transport function, or enable correction of protein folding defects.

Summary

Loss-of-function mutations in the SLC13A5 causes SLC13A5-Epilepsy, a devastating disease characterized by neonatal epilepsy. Currently no cure is available. We clarify the molecular-level defects to guide future developments for phenotype-specific treatment of disease-causing mutations.

Predicting A/B compartments from histone modifications using deep learning

The three-dimensional organization of genomes plays a crucial role in essential biological processes. The segregation of chromatin into A and B compartments highlights regions of activity and inactivity, providing a window into the genomic activities specific to each cell type. Yet, the steep costs associated with acquiring Hi-C data, necessary for studying this compartmentalization across various cell types, pose a significant barrier in studying cell type specific genome organization. To address this, we present a prediction tool called compartment prediction using recurrent neural networks (CoRNN), which predicts compartmentalization of 3D genome using histone modification enrichment. CoRNN demonstrates robust cross-cell-type prediction of A/B compartments with an average AuROC of 90.9%. Cell-type-specific predictions align well with known functional elements, with H3K27ac and H3K36me3 identified as highly predictive histone marks. We further investigate our mispredictions and found that they are located in regions with ambiguous compartmental status. Furthermore, our model's generalizability is validated by predicting compartments in independent tissue samples, which underscores its broad applicability.

Trigeminal innervation and tactile responses in mouse tongue

The mammalian tongue is richly innervated with somatosensory, gustatory and motor fibers. These form the basis of many ethologically important functions such as eating, speaking and social grooming. Despite its high tactile acuity and sensitivity, the neural basis of tongue mechanosensation remains largely mysterious. Here we explored the organization of mechanosensory afferents in the tongue and found that each lingual papilla is innervated by Piezo2 + trigeminal neurons. Notably, each fungiform papilla contained highly specialized ring-like sensory neuron terminations that asymmetrically circumscribe the taste buds. Myelinated lingual afferents in the mouse lingual papillae did not form corpuscular sensory end organs but rather had only free nerve endings. In vivo single-unit recordings from the trigeminal ganglion revealed lingual low-threshold mechanoreceptors (LTMRs) with conduction velocities in the Aδ range or above and distinct adaptation properties ranging from intermediately adapting (IA) to rapidly adapting (RA). IA units were sensitive to both static indentation and stroking, while RA units had a preference for tangential forces applied by stroking. Lingual LTMRs were not directly responsive to rapid cooling or chemicals that can induce astringent or numbing sensations. Sparse labeling of lingual afferents in the tongue revealed distinct terminal morphologies and innervation patterns in fungiform and filiform papillae. Together, our results indicate that fungiform papillae are mechanosensory structures, while suggesting a simple model that links the functional and anatomical properties of tactile sensory neurons in the tongue.

Cerebral perivascular spaces as predictors of dementia risk and accelerated brain atrophy

Cerebral small vessel disease, an important risk factor for dementia, lacks robust, in vivo measurement methods. Perivascular spaces (PVS) on brain MRI are surrogates for small parenchymal blood vessels and their perivascular compartment, and may relate to brain health. We developed a novel, robust algorithm to automatically assess PVS count and size on MRI, and investigated their relationship with dementia risk and brain atrophy. We analyzed 46,478 clinical measurements of cognitive functioning and 20,845 brain MRI scans from 10,004 participants (71.1±9.7 years-old, 56.6% women). Fewer PVS and larger PVS diameter at baseline were associated with higher dementia risk and accelerated brain atrophy. Longitudinal trajectories of PVS markers were significantly different in non-demented individuals who converted to dementia compared with non-converters. In simulated placebo-controlled trials for treatments targeting cognitive decline, screening out participants less likely to develop dementia based on our PVS markers enhanced the power of the trial. These novel radiographic cerebrovascular markers may improve risk-stratification of individuals, potentially reducing cost and increasing throughput of clinical trials to combat dementia.

Chemical manipulation of m1A mediates its detection in human tRNA

N 1-methyl adenosine (m1A) is a widespread RNA modification present in tRNA, rRNA, and mRNA. m1A modification sites in tRNAs are evolutionarily conserved and its formation on tRNA is catalyzed by methyltransferase TRMT61A and TRMT6 complex. m1A promotes translation initiation and elongation. Due to its positive charge under physiological conditions, m1A can notably modulate RNA structure. It also blocks Watson-Crick-Franklin base-pairing and causes mutation and truncation during reverse transcription. Several misincorporation-based high-throughput sequencing methods have been developed to sequence m1A. In this study, we introduce a reduction-based m1A sequencing (red-m1A-seq). We report that NaBH4 reduction of m1A can improve the mutation and readthrough rates using commercially available RT enzymes to give a better positive signature, while alkaline-catalyzed Dimroth rearrangement can efficiently convert m1A to m6A to provide good controls, allowing the detection of m1A with higher sensitivity and accuracy. We applied red-m1A-seq to sequence human small RNA, and we not only detected all the previously reported tRNA m1A sites, but also new m1A sites in mt-tRNAAsn-GTT and 5.8S rRNA.

The sodium-proton exchangers sNHE and NHE1 control plasma membrane hyperpolarization in mouse sperm

Sperm capacitation, crucial for fertilization, occurs in the female reproductive tract and can be replicated in vitro using a medium rich in bicarbonate, calcium, and albumin. These components trigger the cAMP-PKA signaling cascade, proposed to promote hyperpolarization of the mouse sperm plasma membrane through activation of SLO3 K+ channel. Hyperpolarization is a hallmark of capacitation: proper membrane hyperpolarization renders higher in vitro fertilizing ability, while Slo3 KO mice are infertile. However, the precise regulation of SLO3 opening remains elusive. Our study challenges the involvement of PKA in this event and reveals the role of Na+/H+ exchangers. During capacitation, calcium increase through CatSper channels activates NHE1, while cAMP directly stimulates the sperm-specific NHE, collectively promoting the alkalinization threshold needed for SLO3 opening. Hyperpolarization then feeds back Na+/H+ activity. Our work is supported by pharmacology, and a plethora of KO mouse models, and proposes a novel pathway leading to hyperpolarization.

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Pre-challenge gut microbial signature predicts RhCMV/SIV vaccine efficacy in rhesus macaques

Background

RhCMV/SIV vaccines protect ∼59% of vaccinated rhesus macaques against repeated limiting-dose intra-rectal exposure with highly pathogenic SIVmac239M, but the exact mechanism responsible for the vaccine efficacy is not known. It is becoming evident that complex interactions exist between gut microbiota and the host immune system. Here we aimed to investigate if the rhesus gut microbiome impacts RhCMV/SIV vaccine-induced protection.

Methods

Three groups of 15 rhesus macaques naturally pre-exposed to RhCMV were vaccinated with RhCMV/SIV vaccines. Rectal swabs were collected longitudinally both before SIV challenge (after vaccination) and post challenge and were profiled using 16S rRNA based microbiome analysis.

Results

We identified ∼2,400 16S rRNA amplicon sequence variants (ASVs), representing potential bacterial species/strains. Global gut microbial profiles were strongly associated with each of the three vaccination groups, and all animals tended to maintain consistent profiles throughout the pre-challenge phase. Despite vaccination group differences, using newly developed compositional data analysis techniques we identified a common gut microbial signature predictive of vaccine protection outcome across the three vaccination groups. Part of this microbial signature persisted even after SIV challenge. We also observed a strong correlation between this microbial signature and an early signature derived from whole blood transcriptomes in the same animals.

Conclusions

Our findings indicate that changes in gut microbiomes are associated with RhCMV/SIV vaccine-induced protection and early host response to vaccination in rhesus macaques.

Evaluating Large Language Models in Extracting Cognitive Exam Dates and Scores

Importance

Large language models (LLMs) are crucial for medical tasks. Ensuring their reliability is vital to avoid false results. Our study assesses two state-of-the-art LLMs (ChatGPT and LlaMA-2) for extracting clinical information, focusing on cognitive tests like MMSE and CDR.

Objective

Evaluate ChatGPT and LlaMA-2 performance in extracting MMSE and CDR scores, including their associated dates.

Methods

Our data consisted of 135,307 clinical notes (Jan 12th, 2010 to May 24th, 2023) mentioning MMSE, CDR, or MoCA. After applying inclusion criteria 34,465 notes remained, of which 765 underwent ChatGPT (GPT-4) and LlaMA-2, and 22 experts reviewed the responses. ChatGPT successfully extracted MMSE and CDR instances with dates from 742 notes. We used 20 notes for fine-tuning and training the reviewers. The remaining 722 were assigned to reviewers, with 309 each assigned to two reviewers simultaneously. Inter-rater-agreement (Fleiss' Kappa), precision, recall, true/false negative rates, and accuracy were calculated. Our study follows TRIPOD reporting guidelines for model validation.

Results

For MMSE information extraction, ChatGPT (vs. LlaMA-2) achieved accuracy of 83% (vs. 66.4%), sensitivity of 89.7% (vs. 69.9%), true-negative rates of 96% (vs 60.0%), and precision of 82.7% (vs 62.2%). For CDR the results were lower overall, with accuracy of 87.1% (vs. 74.5%), sensitivity of 84.3% (vs. 39.7%), true-negative rates of 99.8% (98.4%), and precision of 48.3% (vs. 16.1%). We qualitatively evaluated the MMSE errors of ChatGPT and LlaMA-2 on double-reviewed notes. LlaMA-2 errors included 27 cases of total hallucination, 19 cases of reporting other scores instead of MMSE, 25 missed scores, and 23 cases of reporting only the wrong date. In comparison, ChatGPT's errors included only 3 cases of total hallucination, 17 cases of wrong test reported instead of MMSE, and 19 cases of reporting a wrong date.

Conclusions

In this diagnostic/prognostic study of ChatGPT and LlaMA-2 for extracting cognitive exam dates and scores from clinical notes, ChatGPT exhibited high accuracy, with better performance compared to LlaMA-2. The use of LLMs could benefit dementia research and clinical care, by identifying eligible patients for treatments initialization or clinical trial enrollments. Rigorous evaluation of LLMs is crucial to understanding their capabilities and limitations.

A Semi-supervised Pipeline for Accurate Neuron Segmentation with Fewer Ground Truth Labels

Recent advancements in two-photon calcium imaging have enabled scientists to record the activity of thousands of neurons with cellular resolution. This scope of data collection is crucial to understanding the next generation of neuroscience questions, but analyzing these large recordings requires automated methods for neuron segmentation. Supervised methods for neuron segmentation achieve state of-the-art accuracy and speed but currently require large amounts of manually generated ground truth training labels. We reduced the required number of training labels by designing a semi-supervised pipeline. Our pipeline used neural network ensembling to generate pseudolabels to train a single shallow U-Net. We tested our method on three publicly available datasets and compared our performance to three widely used segmentation methods. Our method outperformed other methods when trained on a small number of ground truth labels and could achieve state-of-the-art accuracy after training on approximately a quarter of the number of ground truth labels as supervised methods. When trained on many ground truth labels, our pipeline attained higher accuracy than that of state-of-the-art methods. Overall, our work will help researchers accurately process large neural recordings while minimizing the time and effort needed to generate manual labels.

Short-chain fatty acids propionate and butyrate control growth and differentiation linked to cellular metabolism

The short-chain fatty acids (SCFA) propionate and butyrate are produced in large amounts by microbial metabolism and have been identified as unique acyl lysine histone marks. In order to better understand the function of these modifications we used ChIP-seq to map the genome-wide location of four short-chain acyl histone marks H3K18pr/bu and H4K12pr/bu in treated and untreated colorectal cancer (CRC) and normal cells, as well as in mouse intestines in vivo . We correlate these marks with open chromatin regions along with gene expression to access the function of the target regions. Our data demonstrate that propionate and butyrate act as promoters of growth, differentiation as well as ion transport. We propose a mechanism involving direct modification of specific genomic regions, resulting in increased chromatin accessibility, and in case of butyrate, opposing effects on the proliferation of normal versus CRC cells.

Identification of extracellular membrane protein ENPP3 as a major cGAMP hydrolase, cementing cGAMP's role as an immunotransmitter

cGAMP is a second messenger that is synthesized in the cytosol upon detection of cytosolic dsDNA and passed between cells to facilitate downstream immune signaling. ENPP1, an extracellular enzyme, was the only metazoan cGAMP hydrolase known to regulate cGAMP levels to dampen anti-cancer immunity. Here, we uncover ENPP3 as the second and only other metazoan cGAMP hydrolase under homeostatic conditions. ENPP3 has a tissue expression pattern distinct from that of ENPP1 and accounts for all remaining cGAMP hydrolysis activity in mice lacking ENPP1. Importantly, we also show that as with ENPP1, selectively abolishing ENPP3's cGAMP hydrolase activity results in diminished cancer growth and metastasis of certain tumor types. Both ENPP1 and ENPP3 are extracellular enzymes, suggesting the dominant role that extracellular cGAMP must play as a mediator of cell-cell innate immune communication. Our work clearly shows that ENPP1 and ENPP3 non-redundantly dampen extracellular cGAMP-STING signaling, pointing to ENPP3 as a new target for cancer immunotherapy.

Transient RNA structures underlie highly pathogenic avian influenza virus genesis

Highly pathogenic avian influenza viruses (HPAIVs) cause severe disease and high fatality in poultry1. They emerge exclusively from H5 and H7 low pathogenic avian influenza viruses (LPAIVs)2. Although insertion of a furin-cleavable multibasic cleavage site (MBCS) in the hemagglutinin gene was identified decades ago as the genetic basis for LPAIV-to-HPAIV transition3,4, the exact mechanisms underlying said insertion have remained unknown. Here we used an innovative combination of bioinformatic models to predict RNA structures forming around the influenza virus RNA polymerase during replication, and circular sequencing5 to reliably detect nucleotide insertions. We show that transient H5 hemagglutinin RNA structures predicted to trap the polymerase on purine-rich sequences drive nucleotide insertions characteristic of MBCSs, providing the first strong empirical evidence of RNA structure involvement in MBCS acquisition. Insertion frequencies at the H5 cleavage site were strongly affected by substitutions in flanking genomic regions altering predicted transient RNA structures. Introduction of H5-like cleavage site sequences and structures into an H6 hemagglutinin resulted in MBCS-yielding insertions never observed before in H6 viruses. Our results demonstrate that nucleotide insertions that underlie H5 HPAIV emergence result from a previously unknown RNA-structure-driven diversity-generating mechanism, which could be shared with other RNA viruses.

Cortical Reorganization after Limb Loss: Bridging the Gap between Basic Science and Clinical Recovery

Despite the increasing incidence and prevalence of amputation across the globe, individuals with acquired limb loss continue to struggle with functional recovery and chronic pain. A more complete understanding of the motor and sensory remodeling of the peripheral and central nervous system that occurs postamputation may help advance clinical interventions to improve the quality of life for individuals with acquired limb loss. The purpose of this article is to first provide background clinical context on individuals with acquired limb loss and then to provide a comprehensive review of the known motor and sensory neural adaptations from both animal models and human clinical trials. Finally, the article bridges the gap between basic science researchers and clinicians that treat individuals with limb loss by explaining how current clinical treatments may restore function and modulate phantom limb pain using the underlying neural adaptations described above. This review should encourage the further development of novel treatments with known neurological targets to improve the recovery of individuals postamputation.Significance Statement In the United States, 1.6 million people live with limb loss; this number is expected to more than double by 2050. Improved surgical procedures enhance recovery, and new prosthetics and neural interfaces can replace missing limbs with those that communicate bidirectionally with the brain. These advances have been fairly successful, but still most patients experience persistent problems like phantom limb pain, and others discontinue prostheses instead of learning to use them daily. These problematic patient outcomes may be due in part to the lack of consensus among basic and clinical researchers regarding the plasticity mechanisms that occur in the brain after amputation injuries. Here we review results from clinical and animal model studies to bridge this clinical-basic science gap.

Phosphorylation regulates tau's phase separation behavior and interactions with chromatin

Tau is a microtubule-associated protein often found in neurofibrillary tangles (NFTs) in the brains of patients with Alzheimer's disease (AD). Beyond this context, mounting evidence suggests that tau localizes into the nucleus, where it may play a role in DNA protection and heterochromatin regulation. Models of tau depletion or pathology show loss of genetically silent heterochromatin, aberrant expression of heterochromatic genes, and transposable element activation. The molecular mechanisms behind these observations are currently unclear. Using in vitro biophysical experiments, here we demonstrate that tau can undergo liquid-liquid phase separation (LLPS) with DNA, mononucleosomes, and reconstituted nucleosome arrays under low salt conditions. Low concentrations of tau promote chromatin compaction and protect DNA from digestion. While the material state of samples at physiological salt is dominated by chromatin oligomerization, tau can still associate strongly and reversibly with nucleosome arrays. These properties are driven by tau's strong interactions with linker and nucleosomal DNA, while magic angle spinning (MAS) solid-state NMR experiments show that tau binding does not drastically alter nucleosome structure and dynamics. In addition, tau co-localizes into droplets formed by nucleosome arrays and phosphorylated HP1α, a key heterochromatin constituent thought to function through an LLPS mechanism. Importantly, LLPS and chromatin interactions are disrupted by aberrant tau hyperphosphorylation. These biophysical properties suggest that tau may directly impact DNA and chromatin accessibility and that loss of these interactions could contribute to the aberrant nuclear effects seen in tau pathology.

Network Profile: Improving Response to Malaria in the Amazon through Identification of Inter-Community Networks and Human Mobility in Border Regions of Ecuador, Peru, and Brazil

Objectives

Understanding human mobility's role on malaria transmission is critical to successful control and elimination. However, common approaches to measuring mobility are ill-equipped for remote regions such as the Amazon. This study develops a network survey to quantify the effect of community connectivity and mobility on malaria transmission.

Design

A community-level network survey.

Setting

We collect data on community connectivity along three river systems in the Amazon basin: the Pastaza river corridor spanning the Ecuador-Peru border; and the Amazon and Javari river corridors spanning the Brazil-Peru border.

Participants

We interviewed key informants in Brazil, Ecuador, and Peru, including from indigenous communities: Shuar, Achuar, Shiwiar, Kichwa, Ticuna, and Yagua. Key informants are at least 18 years of age and are considered community leaders.

Primary outcome

Weekly, community-level malaria incidence during the study period.

Methods

We measure community connectivity across the study area using a respondent driven sampling design. Forty-five communities were initially selected: 10 in Brazil, 10 in Ecuador, and 25 in Peru. Participants were recruited in each initial node and administered a survey to obtain data on each community's mobility patterns. Survey responses were ranked and the 2-3 most connected communities were then selected and surveyed. This process was repeated for a third round of data collection. Community network matrices will be linked with eadch country's malaria surveillance system to test the effects of mobility on disease risk.

Findings

To date, 586 key informants were surveyed from 126 communities along the Pastaza river corridor. Data collection along the Amazon and Javari river corridors is ongoing. Initial results indicate that network sampling is a superior method to delineate migration flows between communities.

Conclusions

Our study provides measures of mobility and connectivity in rural settings where traditional approaches are insufficient, and will allow us to understand mobility's effect on malaria transmission.

Tyrosine transfer RNA levels and modifications during blood-feeding and vitellogenesis in the mosquito, Aedes aegypti

Mosquitoes such as Aedes aegypti must consume a blood meal for the nutrients necessary for egg production. Several transcriptome and proteome changes occur post blood meal that likely corresponds with codon usage alterations. Transfer RNA (tRNA) is the adapter molecule that reads messenger RNA (mRNA) codons to add the appropriate amino acid during protein synthesis. Chemical modifications to tRNA enhance codons' decoding, improving the accuracy and efficiency of protein synthesis. Here, we examined tRNA modifications and transcripts associated with the blood meal and subsequent periods of vitellogenesis in A. aegypti. More specifically, we assessed tRNA transcript abundance and modification levels in the fat body at critical times post blood-feeding. Based on a combination of alternative codon usage and identification of particular modifications, we identified that increased transcription of tyrosine tRNAs is likely critical during the synthesis of egg yolk proteins in the fat body following a blood meal. Altogether, changes in both the abundance and modification of tRNA are essential factors in the process of vitellogenin production after blood-feeding in mosquitoes.

Longitudinal Profiling of the Intestinal Microbiome in Children with Cystic Fibrosis Treated with Elexacaftor-Tezacaftor-Ivacaftor

The intestinal microbiome influences growth and disease progression in children with cystic fibrosis (CF). Elexacaftor-tezacaftor-ivacaftor (ELX/TEZ/IVA), the newest pharmaceutical modulator for CF, restores function of the pathogenic mutated CFTR channel. We performed a single-center longitudinal analysis of the effect of ELX/TEZ/IVA on the intestinal microbiome, intestinal inflammation, and clinical parameters in children with CF. Following ELX/TEZ/IVA, children with CF had significant improvements in BMI, ppFEV1 and required fewer antibiotics for respiratory infections. Intestinal microbiome diversity increased following ELX/TEZ/IVA coupled with a decrease in the intestinal carriage of Staphylococcus aureus, the predominant respiratory pathogen in children with CF. There was a reduced abundance of microbiome-encoded antibiotic-resistance genes. Microbial pathways for aerobic respiration were reduced after ELX/TEZ/IVA. The abundance of microbial acid tolerance genes was reduced, indicating microbial adaptation to increased CFTR function. In all, this study represents the first comprehensive analysis of the intestinal microbiome in children with CF receiving ELX/TEZ/IVA.

OprF functions as a latch to direct Outer Membrane Vesicle release in Pseudomonas aeruginosa

Bacterial Outer Membrane Vesicles (OMVs) contribute to virulence, competition, immune avoidance and communication. This has led to great interest in how they are formed. To date, investigation has focused almost exclusively on what controls the initiation of OMV biogenesis. Regardless of the mechanism of initiation, all species face a similar challenge before an OMV can be released: How does the OM detach from the underlying peptidoglycan (PG) in regions that will ultimately bulge and then vesiculate? The OmpA family of OM proteins (OprF in P. aeruginosa) is widely conserved and unusually abundant in OMVs across species considering their major role in PG attachment. OmpA homologs also have the interesting ability to adopt both PG-bound (two-domain) and PG-released (one-domain) conformations. Using targeted deletion of the PG-binding domain we showed that loss of cell wall association, and not general membrane destabilization, is responsible for hypervesiculation in OprF-modified strains. We therefore propose that OprF functions as a 'latch', capable of releasing PG in regions destined to become OMVs. To test this hypothesis, we developed a protocol to assess OprF conformation in live cells and purified OMVs. While >90% of OprF proteins exist in the two-domain conformation in the OM of cells, we show that the majority of OprF in OMVs is present in the one-domain conformation. With this work, we take some of the first steps in characterizing late-stage OMV biogenesis and identify a family of proteins whose critical role can be explained by their unique ability to fold into two distinct conformations.

Neural Circuit Transitions Supporting Developmentally Specific Social Behavior

Environmentally appropriate social behavior is critical for survival across the lifespan. To support this flexible behavior, the brain must rapidly perform numerous computations taking into account sensation, memory, motor-control, and many other systems. Further complicating this process, individuals must perform distinct social behaviors adapted to the unique demands of each developmental stage; indeed, the social behaviors of the newborn would not be appropriate in adulthood and vice versa. However, our understanding of the neural circuit transitions supporting these behavioral transitions has been limited. Recent advances in neural circuit dissection tools, as well as adaptation of these tools for use at early time points, has helped uncover several novel mechanisms supporting developmentally appropriate social behavior. This review, and associated Minisymposium, bring together social neuroscience research across numerous model organisms and ages. Together, this work highlights developmentally regulated neural mechanisms and functional transitions in the roles of the sensory cortex, prefrontal cortex, amygdala, habenula, and the thalamus to support social interaction from infancy to adulthood. These studies underscore the need for synthesis across varied model organisms and across ages to advance our understanding of flexible social behavior.

Evolutionary origin and structural ligand mimicry by the inserted domain of alpha-integrin proteins

Heterodimeric integrin proteins transmit signals through conformational changes upon ligand binding between their alpha (α) and beta (β) subunits. Early in chordate evolution, some α subunits acquired an "inserted" (I) domain, which expanded their ligand binding capacity but simultaneously obstructed the ancestral ligand-binding pocket. While this would seemingly impede conventional ligand-mediated integrin activation, it was proposed that the I domain itself could serve both as a ligand replacement and an activation trigger. Here, we provide compelling evidence in support of this longstanding hypothesis using high-resolution cryo-electron microscopy structures of two distinct integrin complexes: the ligand-free and E-cadherin-bound states of the αEβ7 integrin with the I domain, as well as the α4β7 integrin lacking the I domain in both a ligand-free state and bound to MadCAM-1. We trace the evolutionary origin of the I domain to an ancestral collagen-collagen interaction domain. Our analyses illuminate how the I domain intrinsically mimics an extrinsic ligand, enabling integrins to undergo the canonical allosteric cascade of conformational activation and dramatically expanding the range of cellular communication mechanisms in vertebrates.

Newborn DNA methylation age differentiates long-term weight trajectory: The Boston Birth Cohort

Background

Gestational age (GEAA) estimated by newborn DNA methylation (GAmAge) is associated with maternal prenatal exposures and immediate birth outcomes. However, the association of GAmAge with long-term overweight or obesity (OWO) trajectories is yet to be determined.

Methods

GAmAge was calculated for 831 children from a US predominantly urban, low-income, multi-ethnic birth cohort using Illumina EPIC array and cord-blood DNA samples. Repeated anthropometric measurements aligned with pediatric primary care schedule allowed us to calculate body-mass-index percentiles (BMIPCT) at specific age and to define long-term weight trajectories from birth to 18 years.

Results

Four BMIPCT trajectory groups described the long-term weight trajectories: stable (consistent OWO: "early OWO"; constant normal weight: "NW") or non-stable (OWO by year 1 of follow-up: "late OWO"; OWO by year 6 of follow-up: "NW to very late OWO") BMIPCT. were used GAmAge was a predictor of long-term obesity, differentiating between group with consistently high BMIPCT and group with normal BMIPCT patterns and groups with late OWO development. Such differentiation can be observed in the age periods of birth to 1year, 3years, 6years, 10years, and 14years (p<0.05 for all; multivariate models adjusted for GEAA, maternal smoking, delivery method, and child's sex). Birth weight was a mediator for the GAmAge effect on OWO status for specific groups at multiple age periods.

Conclusions

GAmAge is associated with BMI trajectories from birth to age 18 years, independent of GEAA and birth weight. If further confirmed, GAmAge may serve as an early biomarker for future OWO risk.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: Immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

Background

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

Results

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~944.2 Mb (6732 fragments, N50=1.067 Mb), comprising 23,598 genes (BUSCO=93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes were seen in African compared to South American lineages.

Conclusions

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Simultaneous multimodal three-photon and optical coherence microscopy of the mouse brain in the 1700 nm optical window in vivo

Multimodal microscopy combining various imaging approaches can provide complementary information about tissue in a single imaging session. Here, we demonstrate a multimodal approach combining three-photon microscopy (3PM) and spectral-domain optical coherence microscopy (SD-OCM). We show that an optical parametric chirped-pulse amplification (OPCPA) laser source, which is the standard source for three-photon fluorescence excitation and third harmonic generation (THG), can be used for simultaneous OCM, 3-photon (3P) fluorescence and THG imaging. We validated the system performance in deep mouse brains in vivo with an OPCPA source operating at 1620 nm center wavelength. We visualized small structures such as myelinated axons, neurons, and large fiber tracts in white matter with high spatial resolution non-invasively using linear and nonlinear contrast at >1 mm depth in intact adult mouse brain. Our results showed that simultaneous OCM and 3PM at the long wavelength window can be conveniently combined for deep tissue imaging in vivo.

Interfacial morphodynamics of proliferating microbial communities

In microbial communities, various cell types often coexist by occupying distinct spatial domains. What determines the shape of the interface between such domains-which in turn influences the interactions between cells and overall community function? Here, we address this question by developing a continuum model of a 2D spatially-structured microbial community with two distinct cell types. We find that, depending on the balance of the different cell proliferation rates and substrate friction coefficients, the interface between domains is either stable and smooth, or unstable and develops finger-like protrusions. We establish quantitative principles describing when these different interfacial behaviors arise, and find good agreement both with the results of previous experimental reports as well as new experiments performed here. Our work thus helps to provide a biophysical basis for understanding the interfacial morphodynamics of proliferating microbial communities, as well as a broader range of proliferating active systems.

ECM-Focused Proteomic Analysis of Ear Punch Regeneration in Acomys Cahirinus

In mammals, significant injury is generally followed by the formation of a fibrotic scar which provides structural integrity but fails to functionally restore damaged tissue. Spiny mice of the genus Acomys represent the first example of full skin autotomy in mammals. Acomys cahirinus has evolved extremely weak skin as a strategy to avoid predation and is able to repeatedly regenerate healthy tissue without scar after severe skin injury or full-thickness ear punches. Extracellular matrix (ECM) composition is a critical regulator of wound repair and scar formation and previous studies have suggested that alterations in its expression may be responsible for the differences in regenerative capacity observed between Mus musculus and A. cahirinus , yet analysis of this critical tissue component has been limited in previous studies by its insolubility and resistance to extraction. Here, we utilize a 2-step ECM-optimized extraction to perform proteomic analysis of tissue composition during wound repair after full-thickness ear punches in A. cahirinus and M. musculus from weeks 1 to 4 post-injury. We observe changes in a wide range of ECM proteins which have been previously implicated in wound regeneration and scar formation, including collagens, coagulation and provisional matrix proteins, and matricryptic signaling peptides. We additionally report differences in crosslinking enzyme activity and ECM protein solubility between Mus and Acomys. Furthermore, we observed rapid and sustained increases in CD206, a marker of pro-regenerative M2 macrophages, in Acomys, whereas little or no increase in CD206 was detected in Mus. Together, these findings contribute to a comprehensive understanding of tissue cues which drive the regenerative capacity of Acomys and identify a number of potential targets for future pro-regenerative therapies.

Label-free, real-time monitoring of membrane binding events at zeptomolar concentrations using frequency-locked optical microresonators

Binding events to elements of the cell membrane act as receptors which regulate cellular function and communication and are the targets of many small molecule drug discovery efforts for agonists and antagonists. Conventional techniques to probe these interactions generally require labels and large amounts of receptor to achieve satisfactory sensitivity. Whispering gallery mode microtoroid optical resonators have demonstrated sensitivity to detect single-molecule binding events. Here, we demonstrate the use of frequency-locked optical microtoroids for characterization of membrane interactions in vitro at zeptomolar concentrations using a supported biomimetic membrane. Arrays of microtoroids were produced using photolithography and subsequently modified with a biomimetic membrane, providing high quality (Q) factors (> 10 6 ) in aqueous environments. Fluorescent recovery after photobleaching (FRAP) experiments confirmed the retained fluidity of the microtoroid supported-lipid membrane with a diffusion coefficient of 3.38 ± 0.26   μm 2 ⋅ s - 1 . Utilizing this frequency-locked membrane-on-a-chip model combined with auto-balanced detection and non-linear post-processing techniques, we demonstrate zeptomolar detection levels The binding of Cholera Toxin B- monosialotetrahexosyl ganglioside (GM1) was monitored in real-time, with an apparent equilibrium dissociation constant k d = 1.53  nM . The measured affiny of the agonist dynorphin A 1-13 to the κ -opioid receptor revealed a k d = 3.1  nM using the same approach. Radioligand binding competition with dynorphin A 1-13 revealed a k d in agreement (1.1 nM) with the unlabeled method. The biosensing platform reported herein provides a highly sensitive real-time characterization of membrane embedded protein binding kinetics, that is rapid and label-free, for toxin screening and drug discovery, among other applications.

Long working hours and risk of hypertensive intracerebral haemorrhage among Japanese workers claiming compensation for overwork-related intracerebral haemorrhage: an unmatched case-control study

OBJECTIVES

Overwork-related cerebrovascular and cardiovascular diseases (CCVDs) constitute a major occupational and public health issue worldwide. However, to our knowledge, few studies have reported the underlying pathophysiological mechanisms. We aimed to determine whether patients with extreme workload have a greater risk of developing hypertensive intracerebral haemorrhage (ICH) located in the deep brain areas than patients without extreme workload. We also determined the association between the number of hours of overtime work and the risk of developing hypertensive ICH.

DESIGN

Unmatched case-control study.

SETTING

Database of patients claiming compensation for overwork-related CCVDs in Japan.

PARTICIPANTS

A total of 1215 patients who claimed overwork-related ICH in Japan, of whom 621 had their compensation claim approved (patients with extreme workload) and 622 did not.

PRIMARY AND SECONDARY OUTCOME MEASURES

Logistic regression analysis was performed to calculate the risk of developing hypertensive ICH in patients with extreme workload compared with those without extreme workload, adjusted for confounders such as age, occupation, smoking status, alcohol consumption and medical history. We also calculated the risk of developing hypertensive ICH in compensated patients by average monthly overtime working hours.

RESULTS

Patients with extreme workloads had a significantly higher odds ratio (OR) for developing hypertensive ICH (1.44, 95% CI: 1.10 to 1.88) than those without extreme workloads. ORs for developing hypertensive ICH according to overtime working hours showed a dose-response relationship; an overtime of 100 hours/month was associated with a significantly higher OR (1.31, 95% CI: 0.89 to 1.91; 1.41, 95% CI: 0.95 to 2.11; and 1.50, 95% CI: 1.01 to 2.22 for 60-79.9, 80-99.9 and≥100 hours/month, respectively) than that for workloads of less than 60 hours/month.

CONCLUSIONS

Regarding Japanese workers, overtime work could be associated with the risk of developing hypertensive ICH, and hypertension may play an important role in overwork-related ICH.

Whole blood mitochondrial copy number in clinical populations with mood disorders: a meta-analysis

Background

Major depressive disorder (MDD) and bipolar disorder (BD), are globally prevalent, contributing to significant disease burden and adverse health outcomes. These mood disorders are associated with changes in many aspects of brain reward pathways, yet cellular and molecular changes in the brain are not readily available in clinical populations. Therefore, the use of biomarkers as proxies for changes in the brain are necessary. The proliferation of mitochondria in blood has emerged as a potentially useful biomarker, yet a clear consensus on how these mood disorders impact mitochondrial DNA copy number (mtDNAcn) has not been reached.

Methods

Following PRISMA guidelines for a systematic search, 22 papers met inclusion criteria for meta-analysis (10 MDD, 10 BD, 2 both MDD and BD). We extracted demographic, disorder, and methodological information with mtDNAcn. Using the metafor package for R, calculated effect sizes were used in random effects or meta regression models for MDD and BD.

Results

Our results show a trending increase in mtDNAcn in patients with MDD, which reaches significance when one study with outlying demographic characteristics is excluded. Overall, there was no effect of BD on mtDNAcn, however, further subgroup and meta-regression analysis indicated the effects on mtDNAcn are dependent on BD type.

Conclusions

Together our data suggest whole blood/leukocyte mtDNAcn may be a useful biomarker for mood disorders, with MDD and BD Type II associated with higher mtDNAcn, and BD Type I associated with lower mtDNAcn. Further study of blood mtDNAcn could predict downstream health outcomes or treatment responsivity in individuals with mood disorders.

Efficient ultrasound-mediated drug delivery to orthotopic liver tumors - Direct comparison of doxorubicin-loaded nanobubbles and microbubbles

Liver metastasis is a major obstacle in treating aggressive cancers, and current therapeutic options often prove insufficient. To overcome these challenges, there has been growing interest in ultrasound-mediated drug delivery using lipid-shelled microbubbles (MBs) and nanobubbles (NBs) as promising strategies for enhancing drug delivery to tumors. Our previous work demonstrated the potential of Doxorubicin-loaded C3F8 NBs (hDox-NB, 280 ± 123 nm) in improving cancer treatment in vitro using low-frequency ultrasound. In this study, we investigated the pharmacokinetics and biodistribution of sonicated hDox-NBs in orthotopic rat liver tumors. We compared their delivery and therapeutic efficiency with size-isolated MBs (hDox-MB, 1104 ± 373 nm). Results showed a similar accumulation of hDox in tumors treated with hDox-MBs and unfocused therapeutic ultrasound (hDox-MB+TUS) and hDox-NB+TUS. However, significantly increased apoptotic cell death in the tumor and fewer off-target apoptotic cells in the normal liver were found upon the treatment with hDox-NB+TUS. The tumor-to-liver apoptotic ratio was elevated 9.4-fold following treatment with hDox-NB+TUS compared to hDox-MB+TUS, suggesting that the therapeutic efficacy and specificity are significantly increased when using hDox-NB+TUS. These findings highlight the potential of this approach as a viable treatment modality for liver tumors. By elucidating the behavior of drug-loaded bubbles in vivo, we aim to contribute to developing more effective liver cancer treatments that could ultimately improve patient outcomes and decrease off-target side effects.

Evidence of virulence and antimicrobial resistance in Streptococcus pneumoniae serotype 16F lineages

Introduction

Due to the emergence of non-vaccine serotypes in vaccinated populations, Streptococcus pneumoniae remains a major global health challenge despite advances in vaccine development. Serotype 16F is among the predominant non-vaccine serotypes identified among vaccinated infants in South Africa (SA).

Aim

To characterise lineages and antimicrobial resistance in 16F isolates obtained from South Africa and placed the local findings in a global context.

Methodology

We analysed 10923 S. pneumoniae carriage isolates obtained from infants recruited as part of a broader SA birth cohort. We inferred serotype, resistance profile for penicillin, chloramphenicol, cotrimoxazole, erythromycin and tetracycline, and Global Pneumococcal Sequence Clusters (GPSCs) from genomic data. To ensure global representation, we also included S. pneumoniae carriage and disease isolates from the Global Pneumococcal Sequencing (GPS) project database (n=19,607, collected from 49 countries across five continents, years covered (1995 - 2018), accessed on 17 th March 2022).

Results

Nine percent (934/10923) of isolates obtained from infants in the Drakenstein community in SA and 2% (419/19607) of genomes in the GPS dataset were serotype 16F. Serotype 16F isolates were from 28 different lineages of S. pneumoniae, with GPSC33 and GPSC46 having the highest proportion of serotype 16F isolates at 26% (346/1353) and 53% (716/1353), respectively. Serotype 16F isolates were identified globally, however, most isolates were collected from Africa. GPSC33 was associated with carriage [OR (95% CI) 0.24 (0.09 - 0.66); p=0.003], while GPSC46 was associated with disease [OR (95% CI) 19.9 (2.56 - 906.50); p=0.0004]. 10% (37/346) and 15% (53/346) of isolates within GPSC33 had genes associated with resistance to penicillin and co-trimoxazole, respectively, and 18% (128/716) of isolates within GPSC46 had genes associated with resistance to co-trimoxazole. Resistant isolates formed genetic clusters which may suggest emerging resistant lineages.

Discussion

Serotype 16F lineages are common in Southern Africa. Some of these lineages are associated with disease, and resistance to penicillin and cotrimoxazole. We recommend continuous genomic surveillance to determine long term impact of serotype 16F lineages on vaccine efficacy and antimicrobial therapy globally. Investing in vaccine strategies that offer protection over a wide range of serotypes/lineages remains essential.

DATA SUMMARY

The sequencing reads for the genomes analysed have been deposited in the European Nucleotide Archive and the accession numbers for each isolate are listed in Supplementary Table1 . Phylogenetic tree of serotype 16F pneumococcal genomes and associated metadata are available for download and visualisation on the Microreact website: Phylogenies of seotype 16F, GPSC33 and GPSC46 are available on the Microreact serotype-16F , GPSC33 and GPSC46 , respectively.

IMPACT STATEMENT

This study shows that serotype 16F lineages are predominant in Southern Africa and are associated with disease and antimicrobial resistance. Although serotype 16F has been included in the newer formulation of the upcoming vaccine formulations of PCV21 and IVT-25, continuous surveillance to determine long term impact of serotype 16F lineages on vaccines and antimicrobial therapy remains essential.

Coxsackievirus infection induces direct pancreatic β-cell killing but poor anti-viral CD8+ T-cell responses

Coxsackievirus B (CVB) infection of pancreatic β cells is associated with β-cell autoimmunity. We investigated how CVB impacts human β cells and anti-CVB T-cell responses. β cells were efficiently infected by CVB in vitro, downregulated HLA Class I and presented few, selected HLA-bound viral peptides. Circulating CD8+ T cells from CVB-seropositive individuals recognized only a fraction of these peptides, and only another sub-fraction was targeted by effector/memory T cells that expressed the exhaustion marker PD-1. T cells recognizing a CVB epitope cross-reacted with the β-cell antigen GAD. Infected β cells, which formed filopodia to propagate infection, were more efficiently killed by CVB than by CVB-reactive T cells. Thus, our in-vitro and ex-vivo data highlight limited T-cell responses to CVB, supporting the rationale for CVB vaccination trials for type 1 diabetes prevention. CD8+ T cells recognizing structural and non-structural CVB epitopes provide biomarkers to differentially follow response to infection and vaccination.

Mapping Human Tissues with Highly Multiplexed RNA in situ Hybridization

In situ transcriptomic techniques promise a holistic view of tissue organization and cell-cell interactions. Recently there has been a surge of multiplexed RNA in situ techniques but their application to human tissues and clinical biopsies has been limited due to their large size, general lower tissue quality and high background autofluorescence. Here we report DART-FISH, a versatile padlock probe-based technology capable of profiling hundreds to thousands of genes in centimeter-sized human tissue sections at cellular resolution. We introduced an omni-cell type cytoplasmic stain, dubbed RiboSoma that substantially improves the segmentation of cell bodies. We developed a computational decoding-by-deconvolution workflow to extract gene spots even in the presence of optical crowding. Our enzyme-free isothermal decoding procedure allowed us to image 121 genes in a large section from the human neocortex in less than 10 hours, where we successfully recapitulated the cytoarchitecture of 20 neuronal and non-neuronal subclasses. Additionally, we demonstrated the detection of transcripts as short as 461 nucleotides, including neuropeptides and discovered new cortical layer markers. We further performed in situ mapping of 300 genes on a diseased human kidney, profiled >20 healthy and pathological cell states, and identified diseased niches enriched in transcriptionally altered epithelial cells and myofibroblasts.

Local interneurons in the murine visual thalamus have diverse receptive fields and can provide feature selective inhibition to relay cells

By influencing the type and quality of information that relay cells transmit, local interneurons in thalamus have a powerful impact on cortex. To define the sensory features that these inhibitory neurons encode, we mapped receptive fields of optogenetically identified cells in the murine dorsolateral geniculate nucleus. Although few in number, local interneurons had diverse types of receptive fields, like their counterpart relay cells. This result differs markedly from visual cortex, where inhibitory cells are typically less selective than excitatory cells. To explore how thalamic interneurons might converge on relay cells, we took a computational approach. Using an evolutionary algorithm to search through a library of interneuron models generated from our results, we show that aggregated output from different groups of local interneurons can simulate the inhibitory component of the relay cell's receptive field. Thus, our work provides proof-of-concept that groups of diverse interneurons can supply feature-specific inhibition to relay cells.

3D reconstructions of parasite development and the intracellular niche of the microsporidian pathogen E. intestinalis

Microsporidia are an early-diverging group of fungal pathogens that infect a wide range of hosts. Several microsporidian species infect humans, and infections can lead to fatal disease in immunocompromised individuals. As obligate intracellular parasites with highly reduced genomes, microsporidia are dependent on metabolites from their hosts for successful replication and development. Our knowledge of how microsporidian parasites develop inside the host remains rudimentary, and our understanding of the intracellular niche occupied by microsporidia has thus far relied largely on 2D TEM images and light microscopy. Here, we use serial block face scanning electron microscopy (SBF-SEM) to capture 3D snapshots of the human-infecting microsporidian, Encephalitozoon intestinalis , within host cells. We track the development of E. intestinalis through its life cycle, which allows us to propose a model for how its infection organelle, the polar tube, is assembled de novo in each developing spore. 3D reconstructions of parasite-infected cells provide insights into the physical interactions between host cell organelles and parasitophorous vacuoles, which contain the developing parasites. The host cell mitochondrial network is substantially remodeled during E. intestinalis infection, leading to mitochondrial fragmentation. SBF-SEM analysis shows changes in mitochondrial morphology in infected cells, and live-cell imaging provides insights into mitochondrial dynamics during infection. Together, our data provide insights into parasite development, polar tube assembly, and microsporidia-induced mitochondrial remodeling in the host cell.

A class of secreted mammalian peptides with potential to expand cell-cell communication

Peptide hormones and neuropeptides are fundamental signaling molecules that control diverse aspects of mammalian homeostasis and physiology. Here we demonstrate the endogenous presence of a sequence diverse class of orphan, blood-borne peptides that we call "capped peptides." Capped peptides are fragments of secreted proteins and defined by the presence of two post-translational modifications - N-terminal pyroglutamylation and C-terminal amidation - which function as chemical "caps" of the intervening sequence. Capped peptides share many regulatory characteristics in common with that of other signaling peptides, including dynamic regulation in blood plasma by diverse environmental and physiologic stimuli. One capped peptide, CAP-TAC1, is a tachykinin neuropeptide-like molecule and a nanomolar agonist of multiple mammalian tachykinin receptors. A second capped peptide, CAP-GDF15, is a 12-mer peptide that reduces food intake and body weight. Capped peptides therefore define a largely unexplored class of circulating molecules with potential to regulate cell-cell communication in mammalian physiology.

Alcohol-associated bowel disease: new insights into pathogenesis

Excessive alcohol drinking can cause pathological changes including carcinogenesis in the digestive tract from mouth to large intestine, but the underlying mechanisms are not fully understood. In this review, we discuss the effects of alcohol on small and large intestinal functions, such as leaky gut, dysbiosis and alterations of intestinal epithelium and gut immune dysfunctions, commonly referred to as alcohol-associated bowel disease (ABD). To date, detailed mechanistic insights into ABD are lacking. Accumulating evidence suggests a pathogenic role of ethanol metabolism in dysfunctions of the intestinal tract. Ethanol metabolism generates acetaldehyde and acetate, which could potentially promote functional disruptions of microbial and host components of the intestinal barrier along the gastrointestinal tract. The potential involvement of acetaldehyde and acetate in the pathogenesis of the underlying ABD, including cancer, is discussed. We also highlight some gaps in knowledge existing in the field of ABD. Finally, we discuss future directions in exploring the role of acetaldehyde and acetate generated during chronic alcohol intake in various pathologies affecting different sites of the intestinal tract.