Loss of developmentally derived Irf8+ macrophages promotes hyperinnervation and arrhythmia in the adult zebrafish heart.

Recent developments in cardiac macrophage biology have broadened our understanding of the critical functions of macrophages in the heart. As a result, there is further interest in understanding the independent contributions of distinct subsets of macrophage to cardiac development and function. Here, we demonstrate that genetic loss of interferon regulatory factor 8 (Irf8)-positive embryonic-derived macrophages significantly disrupts cardiac conduction, chamber function, and innervation in adult zebrafish. At 4 months post-fertilization (mpf), homozygous irf8st96/st96 mutants have significantly shortened atrial action potential duration and significant differential expression of genes involved in cardiac contraction. Functional in vivo assessments via electro- and echocardiograms at 12 mpf reveal that irf8 mutants are arrhythmogenic and exhibit diastolic dysfunction and ventricular stiffening. To identify the molecular drivers of the functional disturbances in irf8 null zebrafish, we perform single cell RNA sequencing and immunohistochemistry, which reveal increased leukocyte infiltration, epicardial activation, mesenchymal gene expression, and fibrosis. Irf8 null hearts are also hyperinnervated and have aberrant axonal patterning, a phenotype not previously assessed in the context of cardiac macrophage loss. Gene ontology analysis supports a novel role for activated epicardial-derived cells (EPDCs) in promoting neurogenesis and neuronal remodeling in vivo. Together, these data uncover significant cardiac abnormalities following embryonic macrophage loss and expand our knowledge of critical macrophage functions in heart physiology and governing homeostatic heart health.

Conumee 2.0: enhanced copy-number variation analysis from DNA methylation arrays for humans and mice.

MOTIVATION: Copy-number variations (CNVs) are common genetic alterations in cancer and their detection may impact tumor classification and therapeutic decisions. However, detection of clinically relevant large and focal CNVs remains challenging when sample material or resources are limited. This has motivated us to create a software tool to infer CNVs from DNA methylation arrays which are often generated as part of clinical routines and in research settings. RESULTS: We present our R package, conumee 2.0, that combines tangent normalization, an adjustable genomic binning heuristic, and weighted circular binary segmentation to utilize DNA methylation arrays for CNV analysis and mitigate technical biases and batch effects. Segmentation results were validated in a lung squamous cell carcinoma dataset from TCGA (n = 367 samples) by comparison to segmentations derived from genotyping arrays (Pearson's correlation coefficient of 0.91). We further introduce a segmented block bootstrapping approach to detect focal alternations that achieved 60.9% sensitivity and 98.6% specificity for deletions affecting CDKN2A/B (60.0% and 96.9% for RB1, respectively) in a low-grade glioma cohort from TCGA (n = 239 samples). Finally, our tool provides functionality to detect and summarize CNVs across large sample cohorts. AVAILABILITY AND IMPLEMENTATION: Conumee 2.0 is available under open-source license at: https://github.com/hovestadtlab/conumee2.

Near-lifespan longitudinal tracking of brain microvascular morphology, topology, and flow in male mice.

In age-related neurodegenerative diseases, pathology often develops slowly across the lifespan. As one example, in diseases such as Alzheimer's, vascular decline is believed to onset decades ahead of symptomology. However, challenges inherent in current microscopic methods make longitudinal tracking of such vascular decline difficult. Here, we describe a suite of methods for measuring brain vascular dynamics and anatomy in mice for over seven months in the same field of view. This approach is enabled by advances in optical coherence tomography (OCT) and image processing algorithms including deep learning. These integrated methods enabled us to simultaneously monitor distinct vascular properties spanning morphology, topology, and function of the microvasculature across all scales: large pial vessels, penetrating cortical vessels, and capillaries. We have demonstrated this technical capability in wild-type and 3xTg male mice. The capability will allow comprehensive and longitudinal study of a broad range of progressive vascular diseases, and normal aging, in key model systems.

Investigating the correlation between early vascular alterations and cognitive impairment in Alzheimer's disease in mice with SD-OCT.

Vascular alterations have recently gained some attention with their strong association with Alzheimer's disease (AD). We conducted a label-free in vivo optical coherence tomography (OCT) longitudinal imaging using an AD mouse model. We achieved the tracking of the same individual vessels over time and conducted an in-depth analysis of temporal dynamics in vasculature and vasodynamics using OCT angiography and Doppler-OCT. The AD group showed an exponential decay in both vessel diameter and blood flow change with the critical timepoint before 20 weeks of age, which precedes cognitive decline observed at 40 weeks of age. Interestingly, for the AD group, the diameter change showed the dominance in arterioles over venules, but no such influence was found in blood flow change. Conversely, three mice groups with early vasodilatory intervention did not show any significant change in both vascular integrity and cognitive function compared to the wild-type group. We found early vascular alterations and confirmed their correlation with cognitive impairment in AD.

Deep Learning and Simulation for the Estimation of Red Blood Cell Flux With Optical Coherence Tomography.

We present a deep learning and simulation-based method to measure cortical capillary red blood cell (RBC) flux using Optical Coherence Tomography (OCT). This method is more accurate than the traditional peak-counting method and avoids any user parametrization, such as a threshold choice. We used data that was simultaneously acquired using OCT and two-photon microscopy to uncover the distribution of parameters governing the height, width, and inter-peak time of peaks in OCT intensity associated with the passage of RBCs. This allowed us to simulate thousands of time-series examples for different flux values and signal-to-noise ratios, which we then used to train a 1D convolutional neural network (CNN). The trained CNN enabled robust measurement of RBC flux across the entire network of hundreds of capillaries.

Viral genome sequence datasets display pervasive evidence of strand-specific substitution biases that are best described using non-reversible nucleotide substitution models.

Background: The vast majority of phylogenetic trees are inferred from molecular sequence data (nucleotides or amino acids) using time-reversible evolutionary models which assume that, for any pair of nucleotide or amino acid characters, the relative rate of X to Y substitution is the same as the relative rate of Y to X substitution. However, this reversibility assumption is unlikely to accurately reflect the actual underlying biochemical and/or evolutionary processes that lead to the fixation of substitutions. Here, we use empirical viral genome sequence data to reveal that evolutionary non-reversibility is pervasive among most groups of viruses. Specifically, we consider two non-reversible nucleotide substitution models: (1) a 6-rate non-reversible model (NREV6) in which Watson-Crick complementary substitutions occur at identical relative rates and which might therefor be most applicable to analyzing the evolution of genomes where both complementary strands are subject to the same mutational processes (such as might be expected for double-stranded (ds) RNA or dsDNA genomes); and (2) a 12-rate non-reversible model (NREV12) in which all relative substitution types are free to occur at different rates and which might therefore be applicable to analyzing the evolution of genomes where the complementary genome strands are subject to different mutational processes (such as might be expected for viruses with single-stranded (ss) RNA or ssDNA genomes). Results: Using likelihood ratio and Akaike Information Criterion-based model tests, we show that, surprisingly, NREV12 provided a significantly better fit to 21/31 dsRNA and 20/30 dsDNA datasets than did the general time reversible (GTR) and NREV6 models with NREV6 providing a better fit than NREV12 and GTR in only 5/30 dsDNA and 2/31 dsRNA datasets. As expected, NREV12 provided a significantly better fit to 24/33 ssDNA and 40/47 ssRNA datasets. Next, we used simulations to show that increasing degrees of strand-specific substitution bias decrease the accuracy of phylogenetic inference irrespective of whether GTR or NREV12 is used to describe mutational processes. However, in cases where strand-specific substitution biases are extreme (such as in SARS-CoV-2 and Torque teno sus virus datasets) NREV12 tends to yield more accurate phylogenetic trees than those obtained using GTR. Conclusion: We show that NREV12 should, be seriously considered during the model selection phase of phylogenetic analyses involving viral genomic sequences.

Deep learning toolbox for automated enhancement, segmentation, and graphing of cortical optical coherence tomography microangiograms.

Optical coherence tomography angiography (OCTA) is becoming increasingly popular for neuroscientific study, but it remains challenging to objectively quantify angioarchitectural properties from 3D OCTA images. This is mainly due to projection artifacts or "tails" underneath vessels caused by multiple-scattering, as well as the relatively low signal-to-noise ratio compared to fluorescence-based imaging modalities. Here, we propose a set of deep learning approaches based on convolutional neural networks (CNNs) to automated enhancement, segmentation and gap-correction of OCTA images, especially of those obtained from the rodent cortex. Additionally, we present a strategy for skeletonizing the segmented OCTA and extracting the underlying vascular graph, which enables the quantitative assessment of various angioarchitectural properties, including individual vessel lengths and tortuosity. These tools, including the trained CNNs, are made publicly available as a user-friendly toolbox for researchers to input their OCTA images and subsequently receive the underlying vascular network graph with the associated angioarchitectural properties.

Determination of confocal profile and curved focal plane for OCT mapping of the attenuation coefficient.

The attenuation coefficient has proven to be a useful tool in numerous biological applications, but accurate calculation is dependent on the characterization of the confocal effect. This study presents a method to precisely determine the confocal effect and its focal plane within a sample by examining the ratio of two optical coherence tomography (OCT) images. The method can be employed to produce a single-value estimate, or a 2D map of the focal plane accounting for the curvature or tilt within the sample. Furthermore, this method is applicable to data obtained with both high numerical aperture (NA) and low-NA lenses, thereby furthering the applicability of the attenuation coefficient to high-NA OCT data. We test and validate this method using standard samples of Intralipid 20% and 5%, improving the accuracy to 99% from 65% compared to the traditional method and preliminarily show applicability to real biological data of glioblastoma acquired in vivo in a murine model.

Consciousness Indexing and Outcome Prediction with Resting-State EEG in Severe Disorders of Consciousness.

We applied the following methods to resting-state EEG data from patients with disorders of consciousness (DOC) for consciousness indexing and outcome prediction: microstates, entropy (i.e. approximate, permutation), power in alpha and delta frequency bands, and connectivity (i.e. weighted symbolic mutual information, symbolic transfer entropy, complex network analysis). Patients with unresponsive wakefulness syndrome (UWS) and patients in a minimally conscious state (MCS) were classified into these two categories by fitting and testing a generalised linear model. We aimed subsequently to develop an automated system for outcome prediction in severe DOC by selecting an optimal subset of features using sequential floating forward selection (SFFS). The two outcome categories were defined as UWS or dead, and MCS or emerged from MCS. Percentage of time spent in microstate D in the alpha frequency band performed best at distinguishing MCS from UWS patients. The average clustering coefficient obtained from thresholding beta coherence performed best at predicting outcome. The optimal subset of features selected with SFFS consisted of the frequency of microstate A in the 2-20 Hz frequency band, path length obtained from thresholding alpha coherence, and average path length obtained from thresholding alpha coherence. Combining these features seemed to afford high prediction power. Python and MATLAB toolboxes for the above calculations are freely available under the GNU public license for non-commercial use ( https://qeeg.wordpress.com ).