3T-MRI Artificial Intelligence in Patients with Invasive Breast Cancer to Predict Distant Metastasis Status: A Pilot Study

BACKGROUND

The incidence of breast cancer metastasis has decreased over the years. However, 20-30% of patients with early breast cancer still die from metastases. The purpose of this study is to evaluate the performance of a Deep Learning Convolutional Neural Networks (CNN) model to predict the risk of distant metastasis using 3T-MRI DCE sequences (Dynamic Contrast-Enhanced).

METHODS

A total of 157 breast cancer patients who underwent staging 3T-MRI examinations from January 2011 to July 2022 were retrospectively examined. Patient data, tumor histological and MRI characteristics, and clinical and imaging follow-up examinations of up to 7 years were collected. Of the 157 MRI examinations, 39/157 patients (40 lesions) had distant metastases, while 118/157 patients (120 lesions) were negative for distant metastases (control group). We analyzed the role of the Deep Learning technique using a single variable size bounding box (SVB) option and employed a Voxel Based (VB) NET CNN model. The CNN performance was evaluated in terms of accuracy, sensitivity, specificity, and area under the ROC curve (AUC).

RESULTS

The VB-NET model obtained a sensitivity, specificity, accuracy, and AUC of 52.50%, 80.51%, 73.42%, and 68.56%, respectively. A significant correlation was found between the risk of distant metastasis and tumor size, and the expression of PgR and HER2.

CONCLUSIONS

We demonstrated a currently insufficient ability of the Deep Learning approach in predicting a distant metastasis status in patients with BC using CNNs.

CNN-Based Approaches with Different Tumor Bounding Options for Lymph Node Status Prediction in Breast DCE-MRI

Simple Summary

Breast cancer represents the most frequent cancer in women in the world. The state of the axillary lymph node is considered an independent prognostic factor and is currently evaluated only with invasive methods. Deep learning approaches, especially the ones based on convolutional neural networks, offer a valid, non-invasive alternative, allowing extraction of large amounts of the quantitative data that are used to build predictive models. The aim of our work is to evaluate the influence of the peritumoral parenchyma through different bounding box techniques on the prediction of the axillary lymph node in breast cancer patients using a deep learning artificial intelligence approach.

Abstract

Background: The axillary lymph node status (ALNS) is one of the most important prognostic factors in breast cancer (BC) patients, and it is currently evaluated by invasive procedures. Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), highlights the physiological and morphological characteristics of primary tumor tissue. Deep learning approaches (DL), such as convolutional neural networks (CNNs), are able to autonomously learn the set of features directly from images for a specific task. Materials and Methods: A total of 155 malignant BC lesions evaluated via DCE-MRI were included in the study. For each patient’s clinical data, the tumor histological and MRI characteristics and axillary lymph node status (ALNS) were assessed. LNS was considered to be the final label and dichotomized (LN+ (27 patients) vs. LN− (128 patients)). Based on the concept that peritumoral tissue contains valuable information about tumor aggressiveness, in this work, we analyze the contributions of six different tumor bounding options to predict the LNS using a CNN. These bounding boxes include a single fixed-size box (SFB), a single variable-size box (SVB), a single isotropic-size box (SIB), a single lesion variable-size box (SLVB), a single lesion isotropic-size box (SLIB), and a two-dimensional slice (2DS) option. According to the characteristics of the volumes considered as inputs, three different CNNs were investigated: the SFB-NET (for the SFB), the VB-NET (for the SVB, SIB, SLVB, and SLIB), and the 2DS-NET (for the 2DS). All the experiments were run in 10-fold cross-validation. The performance of each CNN was evaluated in terms of accuracy, sensitivity, specificity, the area under the ROC curve (AUC), and Cohen’s kappa coefficient (K). Results: The best accuracy and AUC are obtained by the 2DS-NET (78.63% and 77.86%, respectively). The 2DS-NET also showed the highest specificity, whilst the highest sensibility was attained by the VB-NET based on the SVB and SIB as bounding options. Conclusion: We have demonstrated that a selective inclusion of the DCE-MRI’s peritumoral tissue increases accuracy in the lymph node status prediction in BC patients using CNNs as a DL approach.

Artificial Intelligence in Bone Metastases: An MRI and CT Imaging Review

(1) Background: The purpose of this review is to study the role of radiomics as a supporting tool in predicting bone disease status, differentiating benign from malignant bone lesions, and characterizing malignant bone lesions. (2) Methods: Two reviewers conducted the literature search independently. Thirteen articles on radiomics as a decision support tool for bone lesions were selected. The quality of the methodology was evaluated according to the radiomics quality score (RQS). (3) Results: All studies were published between 2018 and 2021 and were retrospective in design. Eleven (85%) studies were MRI-based, and two (15%) were CT-based. The sample size was <200 patients for all studies. There is significant heterogeneity in the literature, as evidenced by the relatively low RQS value (average score = 22.6%). There is not a homogeneous protocol used for MRI sequences among the different studies, although the highest predictive ability was always obtained in T2W-FS. Six articles (46%) reported on the potential application of the model in a clinical setting with a decision curve analysis (DCA). (4) Conclusions: Despite the variability in the radiomics method application, the similarity of results and conclusions observed is encouraging. Substantial limits were found; prospective and multicentric studies are needed to affirm the role of radiomics as a supporting tool.

AIforCOVID: Predicting the clinical outcomes in patients with COVID-19 applying AI to chest-X-rays. An Italian multicentre study

Recent epidemiological data report that worldwide more than 53 million people have been infected by SARS-CoV-2, resulting in 1.3 million deaths. The disease has been spreading very rapidly and few months after the identification of the first infected, shortage of hospital resources quickly became a problem. In this work we investigate whether artificial intelligence working with chest X-ray (CXR) scans and clinical data can be used as a possible tool for the early identification of patients at risk of severe outcome, like intensive care or death. Indeed, further to induce lower radiation dose than computed tomography (CT), CXR is a simpler and faster radiological technique, being also more widespread. In this respect, we present three approaches that use features extracted from CXR images, either handcrafted or automatically learnt by convolutional neuronal networks, which are then integrated with the clinical data. As a further contribution, this work introduces a repository that collects data from 820 patients enrolled in six Italian hospitals in spring 2020 during the first COVID-19 emergency. The dataset includes CXR images, several clinical attributes and clinical outcomes. Exhaustive evaluation shows promising performance both in 10-fold and leave-one-centre-out cross-validation, suggesting that clinical data and images have the potential to provide useful information for the management of patients and hospital resources.

On the Additional Information Provided by 3T-MRI ADC in Predicting Tumor Cellularity and Microscopic Behavior

BACKGROUND

to evaluate whether Apparent Diffusion Coefficient (ADC) values of invasive breast cancer, provided by 3T Diffusion Weighted-Images (DWI), may represent a non-invasive predictor of pathophysiologic tumor aggressiveness.

METHODS

100 Patients with histologically proven invasive breast cancers who underwent a 3T-MRI examination were included in the study. All MRI examinations included dynamic contrast-enhanced and DWI/ADC sequences. ADC value were calculated for each lesion. Tumor grade was determined according to the Nottingham Grading System, and immuno-histochemical analysis was performed to assess molecular receptors, cellularity rate, on both biopsy and surgical specimens, and proliferation rate (Ki-67 index). Spearman's Rho test was used to correlate ADC values with histological (grading, Ki-67 index and cellularity) and MRI features. ADC values were compared among the different grading (G1, G2, G3), Ki-67 (<20% and >20%) and cellularity groups (<50%, 50-70% and >70%), using Mann-Whitney and Kruskal-Wallis tests. ROC curves were performed to demonstrate the accuracy of the ADC values in predicting the grading, Ki-67 index and cellularity groups.

RESULTS

ADC values correlated significantly with grading, ER receptor status, Ki-67 index and cellularity rates. ADC values were significantly higher for G1 compared with G2 and for G1 compared with G3 and for Ki-67 < 20% than Ki-67 > 20%. The Kruskal-Wallis test showed that ADC values were significantly different among the three grading groups, the three biopsy cellularity groups and the three surgical cellularity groups. The best ROC curves were obtained for the G3 group (AUC of 0.720), for G2 + G3 (AUC of 0.835), for Ki-67 > 20% (AUC of 0.679) and for surgical cellularity rate > 70% (AUC of 0.805).

CONCLUSIONS

3T-DWI ADC is a direct predictor of cellular aggressiveness and proliferation in invasive breast carcinoma, and can be used as a supporting non-invasive factor to characterize macroscopic lesion behavior especially before surgery.

Deep Reinforcement Learning for Fractionated Radiotherapy in Non-Small Cell Lung Carcinoma

Lung cancer is by far the leading cause of cancer death among both men and women. Radiation therapy is one of the main approaches to lung cancer treatment, and its planning is crucial for the therapy outcome. However, the current practice that uniformly delivers the dose does not take into account the patient-specific tumour features that may affect treatment success. Since radiation therapy is by its very nature a sequential procedure, Deep Reinforcement Learning (DRL) is a well-suited methodology to overcome this limitation. In this respect, in this work we present a DRL controller optimizing the daily dose fraction delivered to the patient on the basis of CT scans collected over time during the therapy, offering a personalized treatment not only for volume adaptation, as currently intended, but also for daily fractionation. Furthermore, this contribution introduces a virtual radiotherapy environment based on a set of ordinary differential equations modelling the tissue radiosensitivity by combining both the effect of the radiotherapy treatment and cell growth. Their parameters are estimated from CT scans routinely collected using the Particle Swarm Optimization algorithm. This permits the DRL to learn the optimal behaviour through an iterative trial and error process with the environment. We performed several experiments considering three rewards functions modelling treatment strategies with different tissue aggressiveness and two exploration strategies for the exploration-exploitation dilemma. The results show that our DRL approach can adapt to radiation therapy treatment, optimizing its behaviour according to the different reward functions and outperforming the current clinical practice.

Visual4DTracker: a tool to interact with 3D + t image stacks

Background

Biological phenomena usually evolves over time and recent advances in high-throughput microscopy have made possible to collect multiple 3D images over time, generating \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3D+t$$\end{document}3D+t (or 4D) datasets. To extract useful information there is the need to extract spatial and temporal data on the particles that are in the images, but particle tracking and feature extraction need some kind of assistance.

Results

This manuscript introduces our new freely downloadable toolbox, the Visual4DTracker. It is a MATLAB package implementing several useful functionalities to navigate, analyse and proof-read the track of each particle detected in any \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3D+t$$\end{document}3D+t stack. Furthermore, it allows users to proof-read and to evaluate the traces with respect to a given gold standard. The Visual4DTracker toolbox permits the users to visualize and save all the generated results through a user-friendly graphical user interface. This tool has been successfully used in three applicative examples. The first processes synthetic data to show all the software functionalities. The second shows how to process a 4D image stack showing the time-lapse growth of Drosophila cells in an embryo. The third example presents the quantitative analysis of insulin granules in living beta-cells, showing that such particles have two main dynamics that coexist inside the cells.

Conclusions

Visual4DTracker is a software package for MATLAB to visualize, handle and manually track \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3D+t$$\end{document}3D+t stacks of microscopy images containing objects such cells, granules, etc.. With its unique set of functions, it remarkably permits the user to analyze and proof-read 4D data in a friendly 3D fashion. The tool is freely available at https://drive.google.com/drive/folders/19AEn0TqP-2B8Z10kOavEAopTUxsKUV73?usp=sharing

Positive tissue transglutaminase antibodies with negative endomysial antibodies: Unresolved issues in diagnosing celiac disease

BACKGROUND

The serological screening for celiac disease (CD) is currently based on the detection of anti-transglutaminase (tTG) IgA antibodies, subsequently confirmed by positive endomysial antibodies (EMA). When an anti-tTG IgA positive/EMA IgA negative result occurs, it can be due either to the lower sensitivity of the EMA test or to the lower specificity of the anti-tTG test. This study aimed at verifying how variation in analytical specificity among different anti-tTG methods could account for this discrepancy.

METHODS

A total of 130 consecutive anti-tTG IgA positive/EMA negative samples were collected from the local screening routine and tested using five anti-tTG IgA commercial assays: two chemiluminescence methods, one fluoroimmunoenzymatic method, one immunoenzymatic method and one multiplex flow immunoassay method.

RESULTS

Twenty three/130 (17.7%) patients were diagnosed with CD. In the other 107 cases a diagnosis of CD was not confirmed. The overall agreement among the five anti-tTG methods ranged from 28.5% to 77.7%. CD condition was more likely linked to the positivity of more than one anti-tTG IgA assay (monopositive = 2.5%, positive with ≥ three methods = 29.5%; p = 0.0004), but it was not related to anti-tTG IgA antibody levels (either positive or borderline; p = 0.5).

CONCLUSIONS

Patients with positive anti-tTG/negative EMA have a low probability of being affected by CD. Given the high variability among methods to measure anti-tTG IgA antibodies, anti-tTG-positive/EMA-negative result must be considered with extreme caution. It is advisable that the laboratory report comments on any discordant results, suggesting to consider the data in the proper clinical context and to refer the patient to a CD reference center for prolonged follow up.

Exploiting Multi-Level Parallelism for Stitching Very Large Microscopy Images

Due to the limited field of view of the microscopes, acquisitions of macroscopic specimens require many parallel image stacks to cover the whole volume of interest. Overlapping regions are introduced among stacks in order to make it possible automatic alignment by means of a 3D stitching tool. Since state-of-the-art microscopes coupled with chemical clearing procedures can generate 3D images whose size exceeds the Terabyte, parallelization is required to keep stitching time within acceptable limits. In the present paper we discuss how multi-level parallelization reduces the execution times of TeraStitcher, a tool designed to deal with very large images. Two algorithms performing dataset partition for efficient parallelization in a transparent way are presented together with experimental results proving the effectiveness of the approach that achieves a speedup close to 300×, when both coarse- and fine-grained parallelism are exploited. Multi-level parallelization of TeraStitcher led to a significant reduction of processing times with no changes in the user interface, and with no additional effort required for the maintenance of code.

Comparison between Bioimpedance Analysis and Dual-Energy X-ray Absorptiometry in assessment of body composition in a cohort of elderly patients aged 65-90 years

We compare bioimpedance analysis (BIA) with dual-energy X-ray absorptiometry (DXA) in the assessment of free fat mass (FFM), fat mass (FM) and percentage of body fat under different conditions in relation to age categories, hydration parameters, body mass index (BMI) and sarcopenia. A cross-sectional analysis of body composition was estimated by BIA and DXA in 379 hospitalized elderly patients. In addition, estimates of FFM, FM and percentage of body fat were investigated across different conditions. Paired t-tests, Bland-Altman plot and intraclass correlation coefficient analysis were used to compare methods. Data showed an underestimation of means (BIA versus DXA) of FFM (women: 0,97 kg, p<0,01; men: 1,99 kg; p<0,01), and an overestimation of both the FM (women: +1,11 kg; p<0,01; men: +1,67 kg; p<0,01) and percentage of body fat (women: +2,07 %, p<0,01; men: +2,82 %, p<0,01). BIA underestimated FFM and overestimated FM and percentage of body fat in patients from the age group of 75 to 85 years, in patients with a total body water content <60%, in underweight and normal weight patients and in patients with sarcopenia (p<0,01). The intraclass coefficient results were indicative of poor reproducibility between BIA and DXA for FFM (women: +0,197; men: +0,250) and FM (women: +0,141; men +0,144). BIA is a good alternative for estimation of FFM and FM only in overweight or obese patients or in patients with good hydration status. BIA, on the other hand, is not an accurate method for assessing FFM in sarcopenic patients.

A radiomic approach for adaptive radiotherapy in non-small cell lung cancer patients

The primary goal of precision medicine is to minimize side effects and optimize efficacy of treatments. Recent advances in medical imaging technology allow the use of more advanced image analysis methods beyond simple measurements of tumor size or radiotracer uptake metrics. The extraction of quantitative features from medical images to characterize tumor pathology or heterogeneity is an interesting process to investigate, in order to provide information that may be useful to guide the therapies and predict survival. This paper discusses the rationale supporting the concept of radiomics and the feasibility of its application to Non-Small Cell Lung Cancer in the field of radiation oncology research. We studied 91 stage III patients treated with concurrent chemoradiation and adaptive approach in case of tumor reduction during treatment. We considered 12 statistics features and 230 textural features extracted from the CT images. In our study, we used an ensemble learning method to classify patients' data into either the adaptive or non-adaptive group during chemoradiation on the basis of the starting CT simulation. Our data supports the hypothesis that a specific signature can be identified (AUC 0.82). In our experience, a radiomic signature mixing semantic and image-based features has shown promising results for personalized adaptive radiotherapy in non-small cell lung cancer.

The inter-observer reading variability in anti-nuclear antibodies indirect (ANA) immunofluorescence test: A multicenter evaluation and a review of the literature

Recently there has been an increase demand for Computer-Aided Diagnosis (CAD) tools to support clinicians in the field of Indirect ImmunoFluorescence (IIF), as the novel digital imaging reading approach can help to overcome the reader subjectivity. Nevertheless, a large multicenter evaluation of the inter-observer reading variability in this field is still missing. This work fills this gap as we evaluated 556 consecutive samples, for a total of 1679 images, collected in three laboratories with IIF expertise using HEp-2 cell substrate (MBL) at 1:80 screening dilution according to conventional procedures. In each laboratory, the images were blindly classified by two experts into three intensity classes: positive, negative, and weak positive. Positive and weak positive ANA-IIF results were categorized by the predominant fluorescence pattern among six main classes. Data were pairwise analyzed and the inter-observer reading variability was measured by Cohen's kappa test, revealing a pairwise agreement little further away than substantial both for fluorescence intensity and for staining pattern recognition (k=0.602 and k=0.627, respectively). We also noticed that the inter-observer reading variability decreases when it is measured with respect to a gold standard classification computed on the basis of labels assigned by the three laboratories. These data show that laboratory agreement improves using digital images and comparing each single human evaluation to potential reference data, suggesting that a solid gold standard is essential to properly make use of CAD systems in routine work lab.

Channel interpolation in TMS-EEG: a quantitative study towards an accurate topographical representation

The co-registration of transcranial magnetic stimulation and electroencephalography (TMS-EEG) is emerging as a successful technique for causally exploring cortical mechanisms and connections. However, various artefacts could affect TMS-EEG signals. Correct artefacted channels reconstruction is crucial to obtain accurate topographical representation and consequently accurate inverse problem solution, in order to map in a proper way the global brain responses after the stimulation of one particular brain region of interest. In this paper, we discuss the problem of artefacted channels interpolation in TMS-EEG signals. Aim of the study was to investigate two different interpolation methods evaluating their performance in two datasets: one constituted by 19 EEG channels montage (low-density spatial resolution) and the other one by 60 EEG channels montage (high-density spatial resolution). In addition, these evaluations took place in two different contexts of application: after the averaging of TMS Evoked Potentials (TEPs) in a time interval to obtain a global information in the considered range, and at fixed latencies 100 ms and 300 ms after the TMS stimulus. The results showed that the global reconstruction error was lower at fixed latencies for the high-density electrodes spatial resolution montage.

Sensorimotor cortex excitability and connectivity in Alzheimer's disease: A TMS-EEG Co-registration study

Several studies have shown that, in spite of the fact that motor symptoms manifest late in the course of Alzheimer's disease (AD), neuropathological progression in the motor cortex parallels that in other brain areas generally considered more specific targets of the neurodegenerative process. It has been suggested that motor cortex excitability is enhanced in AD from the early stages, and that this is related to disease's severity and progression. To investigate the neurophysiological hallmarks of motor cortex functionality in early AD we combined transcranial magnetic stimulation (TMS) with electroencephalography (EEG). We demonstrated that in mild AD the sensorimotor system is hyperexcitable, despite the lack of clinically evident motor manifestations. This phenomenon causes a stronger response to stimulation in a specific time window, possibly due to locally acting reinforcing circuits, while network activity and connectivity is reduced. These changes could be interpreted as a compensatory mechanism allowing for the preservation of sensorimotor programming and execution over a long period of time, regardless of the disease's progression. Hum Brain Mapp 37:2083-2096, 2016. © 2016 Wiley Periodicals, Inc.

Computer-based automatic identification of neurons in gigavoxel-sized 3D human brain images

Achieving a comprehensive knowledge of the human brain cytoarchitecture is a fundamental step to understand how the nervous system works, i.e., one of the greatest challenge of 21(st) century science. The recent development of biological tissue labeling and automated microscopic imaging systems has permitted to acquire images at the micro-resolution, which produce a huge quantity of data that cannot be manually analyzed. In case of mammals brain, automatic methods to extract objective information at the microscale have been applied until now to mice, macaque and cat 3D volume images. Here we report a method to automatically localize neurons in a sample of human brain removed during a surgical procedure for the treatments of drug resistant epilepsy in a child with hemimegalencephaly, whose neurons and neurites were fluorescence labelled and finally imaged using the two-photon fluorescence microscope. The method provides the map of both parvalbuminergic neurons and all other cells nuclei with a satisfactory f-score measured using more than two thousand human labelled soma.

Label-free near-infrared reflectance microscopy as a complimentary tool for two-photon fluorescence brain imaging

In vivo two-photon imaging combined with targeted fluorescent indicators is currently extensively used for attaining critical insights into brain functionality and structural plasticity. Additional information might be gained from back-scattered photons from the near-infrared (NIR) laser without introducing any exogenous labelling. Here, we describe a complimentary and versatile approach that, by collecting the reflected NIR light, provides structural details on axons and blood vessels in the brain, both in fixed samples and in live animals under a cranial window. Indeed, by combining NIR reflectance and two-photon imaging of a slice of hippocampus from a Thy1-GFPm mouse, we show the presence of randomly oriented axons intermingled with sparsely fluorescent neuronal processes. The back-scattered photons guide the contextualization of the fluorescence structure within brain atlas thanks to the recognition of characteristic hippocampal structures. Interestingly, NIR reflectance microscopy allowed the label-free detection of axonal elongations over the superficial layers of mouse cortex under a cranial window in vivo. Finally, blood flow can be measured in live preparations, thus validating label free NIR reflectance as a tool for monitoring hemodynamic fluctuations. The prospective versatility of this label-free technique complimentary to two-photon fluorescence microscopy is demonstrated in a mouse model of photothrombotic stroke in which the axonal degeneration and blood flow remodeling can be investigated.

Comprehensive optical and data management infrastructure for high-throughput light-sheet microscopy of whole mouse brains

Comprehensive mapping and quantification of neuronal projections in the central nervous system requires high-throughput imaging of large volumes with microscopic resolution. To this end, we have developed a confocal light-sheet microscope that has been optimized for three-dimensional (3-D) imaging of structurally intact clarified whole-mount mouse brains. We describe the optical and electromechanical arrangement of the microscope and give details on the organization of the microscope management software. The software orchestrates all components of the microscope, coordinates critical timing and synchronization, and has been written in a versatile and modular structure using the LabVIEW language. It can easily be adapted and integrated to other microscope systems and has been made freely available to the light-sheet community. The tremendous amount of data routinely generated by light-sheet microscopy further requires novel strategies for data handling and storage. To complete the full imaging pipeline of our high-throughput microscope, we further elaborate on big data management from streaming of raw images up to stitching of 3-D datasets. The mesoscale neuroanatomy imaged at micron-scale resolution in those datasets allows characterization and quantification of neuronal projections in unsectioned mouse brains.

MONSTER v1.1: a tool to extract and search for RNA non-branching structures

Background

Detection of RNA structure similarities is still one of the major computational problems in the discovery of RNA functions. A case in point is the study of the new appreciated long non-coding RNAs (lncRNAs), emerging as new players involved in many cellular processes and molecular interactions. Among several mechanisms of action, some lncRNAs show specific substructures that are likely to be instrumental for their functioning. For instance, it has been reported in literature that some lncRNAs have a guiding or scaffolding role by binding chromatin-modifying protein complexes. Thus, a functionally characterized lncRNA (reference) can be used to infer the function of others that are functionally unknown (target), based on shared structural motifs.

Methods

In our previous work we presented a tool, MONSTER v1.0, able to identify structural motifs shared between two full-length RNAs. Our procedure is mainly composed of two ad-hoc developed algorithms: nbRSSP_extractor for characterizing the folding of an RNA sequence by means of a sequence-structure descriptor (i.e., an array of non-overlapping substructures located on the RNA sequence and coded by dot-bracket notation); and SSD_finder, to enable an effective search engine for groups of matches (i.e., chains) common to the reference and target RNA based on a dynamic programming approach with a new score function. Here, we present an updated version of the previous one (MONSTER v1.1) accounting for the peculiar feature of lncRNAs that are not expected to have a unique fold, but appear to fluctuate among a large number of equally-stable folds. In particular, we improved our SSD_finder algorithm in order to take into account all the alternative equally-stable structures.

Results

We present an application of MONSTER v1.1 on lincRNAs, which are a specific class of lncRNAs located in genomic regions which do not overlap protein-coding genes. In particular, we provide reliable predictions of the shared chains between HOTAIR, ANRIL and COLDAIR. The latter are lincRNAs which interact with the same protein complexes of the Polycomb group and hence they are expected to share structural motifs.

Software availability: the software package is provided as additional file 1 ("archive_updated.zip").

Quantitative neuroanatomy of all Purkinje cells with light sheet microscopy and high-throughput image analysis

Characterizing the cytoarchitecture of mammalian central nervous system on a brain-wide scale is becoming a compelling need in neuroscience. For example, realistic modeling of brain activity requires the definition of quantitative features of large neuronal populations in the whole brain. Quantitative anatomical maps will also be crucial to classify the cytoarchtitectonic abnormalities associated with neuronal pathologies in a high reproducible and reliable manner. In this paper, we apply recent advances in optical microscopy and image analysis to characterize the spatial distribution of Purkinje cells (PCs) across the whole cerebellum. Light sheet microscopy was used to image with micron-scale resolution a fixed and cleared cerebellum of an L7-GFP transgenic mouse, in which all PCs are fluorescently labeled. A fast and scalable algorithm for fully automated cell identification was applied on the image to extract the position of all the fluorescent PCs. This vectorized representation of the cell population allows a thorough characterization of the complex three-dimensional distribution of the neurons, highlighting the presence of gaps inside the lamellar organization of PCs, whose density is believed to play a significant role in autism spectrum disorders. Furthermore, clustering analysis of the localized somata permits dividing the whole cerebellum in groups of PCs with high spatial correlation, suggesting new possibilities of anatomical partition. The quantitative approach presented here can be extended to study the distribution of different types of cell in many brain regions and across the whole encephalon, providing a robust base for building realistic computational models of the brain, and for unbiased morphological tissue screening in presence of pathologies and/or drug treatments.

A versatile clearing agent for multi-modal brain imaging

Extensive mapping of neuronal connections in the central nervous system requires high-throughput µm-scale imaging of large volumes. In recent years, different approaches have been developed to overcome the limitations due to tissue light scattering. These methods are generally developed to improve the performance of a specific imaging modality, thus limiting comprehensive neuroanatomical exploration by multi-modal optical techniques. Here, we introduce a versatile brain clearing agent (2,2′-thiodiethanol; TDE) suitable for various applications and imaging techniques. TDE is cost-efficient, water-soluble and low-viscous and, more importantly, it preserves fluorescence, is compatible with immunostaining and does not cause deformations at sub-cellular level. We demonstrate the effectiveness of this method in different applications: in fixed samples by imaging a whole mouse hippocampus with serial two-photon tomography; in combination with CLARITY by reconstructing an entire mouse brain with light sheet microscopy and in translational research by imaging immunostained human dysplastic brain tissue.

Neurophysiological features of motor cortex excitability and plasticity in Subcortical Ischemic Vascular Dementia: a TMS mapping study

OBJECTIVE

To evaluate neurophysiological features of M1 excitability and plasticity in Subcortical Ischemic Vascular Dementia (SIVD), by means of a TMS mapping study.

METHODS

Seven SIVD and nine AD patients, along with nine control subjects were tested. The M1 excitability was studied by resting thresholds, area and volume of active cortical sites for forearm and hand's examined muscles. For M1 plasticity, coordinates of the hot-spot and the center of gravity (CoG) were evaluated. The correlation between the degree of hyperexcitability and the amount of M1 plastic rearrangement was also calculated.

RESULTS

Multivariate analysis of excitability measures demonstrated similarly enhanced cortical excitability in AD and SIVD patients with respect to controls. SIVD patients showed a medial and frontal shift of CoG from the hot-spot, not statistically different from that observed in AD. A significant direct correlation was seen between parameters related to cortical excitability and those related to cortical plasticity.

CONCLUSIONS

The results suggest the existence of common compensatory mechanisms in different kind of dementing diseases supporting the idea that cortical hyperexcitability can promote cortical plasticity.

SIGNIFICANCE

This study characterizes neurophysiological features of motor cortex excitability and plasticity in SIVD, providing new insights on the correlation between cortical excitability and plasticity.

Large-scale automated identification of mouse brain cells in confocal light sheet microscopy images

MOTIVATION

Recently, confocal light sheet microscopy has enabled high-throughput acquisition of whole mouse brain 3D images at the micron scale resolution. This poses the unprecedented challenge of creating accurate digital maps of the whole set of cells in a brain.

RESULTS

We introduce a fast and scalable algorithm for fully automated cell identification. We obtained the whole digital map of Purkinje cells in mouse cerebellum consisting of a set of 3D cell center coordinates. The method is accurate and we estimated an F1 measure of 0.96 using 56 representative volumes, totaling 1.09 GVoxel and containing 4138 manually annotated soma centers.

AVAILABILITY AND IMPLEMENTATION

Source code and its documentation are available at http://bcfind.dinfo.unifi.it/. The whole pipeline of methods is implemented in Python and makes use of Pylearn2 and modified parts of Scikit-learn. Brain images are available on request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Virtual finger boosts three-dimensional imaging and microsurgery as well as terabyte volume image visualization and analysis

Three-dimensional (3D) bioimaging, visualization and data analysis are in strong need of powerful 3D exploration techniques. We develop virtual finger (VF) to generate 3D curves, points and regions-of-interest in the 3D space of a volumetric image with a single finger operation, such as a computer mouse stroke, or click or zoom from the 2D-projection plane of an image as visualized with a computer. VF provides efficient methods for acquisition, visualization and analysis of 3D images for roundworm, fruitfly, dragonfly, mouse, rat and human. Specifically, VF enables instant 3D optical zoom-in imaging, 3D free-form optical microsurgery, and 3D visualization and annotation of terabytes of whole-brain image volumes. VF also leads to orders of magnitude better efficiency of automated 3D reconstruction of neurons and similar biostructures over our previous systems. We use VF to generate from images of 1,107 Drosophila GAL4 lines a projectome of a Drosophila brain.

The classification of Crithidia luciliae immunofluorescence test (CLIFT) using a novel automated system

Introduction

In recent years, there has been an increased demand for computer-aided diagnosis (CAD) tools to support clinicians in the field of indirect immunofluorescence. To this aim, academic and industrial research is focusing on detecting antinuclear, anti-neutrophil, and anti-double-stranded (anti-dsDNA) antibodies. Within this framework, we present a CAD system for automatic analysis of dsDNA antibody images using a multi-step classification approach. The final classification of a well is based on the classification of all its images, and each image is classified on the basis of the labeling of its cells.

Methods

We populated a database of 342 images—74 positive (21.6%) and 268 negative (78.4%)— belonging to 63 consecutive sera: 15 positive (23.8%) and 48 negative (76.2%). We assessed system performance by using k-fold cross-validation. Furthermore, we successfully validated the recognition system on 83 consecutive sera, collected by using different equipment in a referral center, counting 279 images: 92 positive (33.0%) and 187 negative (67.0%).

Results

With respect to well classification, the system correctly classified 98.4% of wells (62 out of 63). Integrating information from multiple images of the same wells recovers the possible misclassifications that occurred at the previous steps (cell and image classification). This system, validated in a clinical routine fashion, provides recognition accuracy equal to 100%.

Conclusion

The data obtained show that automation is a viable alternative for Crithidia luciliae immunofluorescence test analysis.

Extensible visualization and analysis for multidimensional images using Vaa3D

Open-Source 3D Visualization-Assisted Analysis (Vaa3D) is a software platform for the visualization and analysis of large-scale multidimensional images. In this protocol we describe how to use several popular features of Vaa3D, including (i) multidimensional image visualization, (ii) 3D image object generation and quantitative measurement, (iii) 3D image comparison, fusion and management, (iv) visualization of heterogeneous images and respective surface objects and (v) extension of Vaa3D functions using its plug-in interface. We also briefly demonstrate how to integrate these functions for complicated applications of microscopic image visualization and quantitative analysis using three exemplar pipelines, including an automated pipeline for image filtering, segmentation and surface generation; an automated pipeline for 3D image stitching; and an automated pipeline for neuron morphology reconstruction, quantification and comparison. Once a user is familiar with Vaa3D, visualization usually runs in real time and analysis takes less than a few minutes for a simple data set.

A Kalman filter approach for denoising and deblurring 3-D microscopy images

This paper proposes a new method for removing noise and blurring from 3D microscopy images. The main contribution is the definition of a space-variant generating model of a 3-D signal, which is capable to stochastically describe a wide class of 3-D images. Unlike other approaches, the space-variant structure allows the model to consider the information on edge locations, if available. A suitable description of the image acquisition process, including blurring and noise, is then associated to the model. A state-space realization is finally derived, which is amenable to the application of standard Kalman filter as an image restoration algorithm. The so obtained method is able to remove, at each spatial step, both blur and noise, via a linear minimum variance recursive one-shot procedure, which does not require the simultaneous processing of the whole image. Numerical results on synthetic and real microscopy images confirm the merit of the approach.

Micron-scale resolution optical tomography of entire mouse brains with confocal light sheet microscopy

Understanding the architecture of mammalian brain at single-cell resolution is one of the key issues of neuroscience. However, mapping neuronal soma and projections throughout the whole brain is still challenging for imaging and data management technologies. Indeed, macroscopic volumes need to be reconstructed with high resolution and contrast in a reasonable time, producing datasets in the TeraByte range. We recently demonstrated an optical method (confocal light sheet microscopy, CLSM) capable of obtaining micron-scale reconstruction of entire mouse brains labeled with enhanced green fluorescent protein (EGFP). Combining light sheet illumination and confocal detection, CLSM allows deep imaging inside macroscopic cleared specimens with high contrast and speed. Here we describe the complete experimental pipeline to obtain comprehensive and human-readable images of entire mouse brains labeled with fluorescent proteins. The clearing and the mounting procedures are described, together with the steps to perform an optical tomography on its whole volume by acquiring many parallel adjacent stacks. We showed the usage of open-source custom-made software tools enabling stitching of the multiple stacks and multi-resolution data navigation. Finally, we illustrated some example of brain maps: the cerebellum from an L7-GFP transgenic mouse, in which all Purkinje cells are selectively labeled, and the whole brain from a thy1-GFP-M mouse, characterized by a random sparse neuronal labeling.

TeraStitcher - a tool for fast automatic 3D-stitching of teravoxel-sized microscopy images

Background

Further advances in modern microscopy are leading to teravoxel-sized tiled 3D images at high resolution, thus increasing the dimension of the stitching problem of at least two orders of magnitude. The existing software solutions do not seem adequate to address the additional requirements arising from these datasets, such as the minimization of memory usage and the need to process just a small portion of data.

Results

We propose a free and fully automated 3D Stitching tool designed to match the special requirements coming out of teravoxel-sized tiled microscopy images that is able to stitch them in a reasonable time even on workstations with limited resources. The tool was tested on teravoxel-sized whole mouse brain images with micrometer resolution and it was also compared with the state-of-the-art stitching tools on megavoxel-sized publicy available datasets. This comparison confirmed that the solutions we adopted are suited for stitching very large images and also perform well on datasets with different characteristics. Indeed, some of the algorithms embedded in other stitching tools could be easily integrated in our framework if they turned out to be more effective on other classes of images. To this purpose, we designed a software architecture which separates the strategies that use efficiently memory resources from the algorithms which may depend on the characteristics of the acquired images.

Conclusions

TeraStitcher is a free tool that enables the stitching of Teravoxel-sized tiled microscopy images even on workstations with relatively limited resources of memory (<8 GB) and processing power. It exploits the knowledge of approximate tile positions and uses ad-hoc strategies and algorithms designed for such very large datasets. The produced images can be saved into a multiresolution representation to be efficiently retrieved and processed. We provide TeraStitcher both as standalone application and as plugin of the free software Vaa3D.

Confocal light sheet microscopy: micron-scale neuroanatomy of the entire mouse brain

Elucidating the neural pathways that underlie brain function is one of the greatest challenges in neuroscience. Light sheet based microscopy is a cutting edge method to map cerebral circuitry through optical sectioning of cleared mouse brains. However, the image contrast provided by this method is not sufficient to resolve and reconstruct the entire neuronal network. Here we combined the advantages of light sheet illumination and confocal slit detection to increase the image contrast in real time, with a frame rate of 10 Hz. In fact, in confocal light sheet microscopy (CLSM), the out-of-focus and scattered light is filtered out before detection, without multiple acquisitions or any post-processing of the acquired data. The background rejection capabilities of CLSM were validated in cleared mouse brains by comparison with a structured illumination approach. We show that CLSM allows reconstructing macroscopic brain volumes with sub-cellular resolution. We obtained a comprehensive map of Purkinje cells in the cerebellum of L7-GFP transgenic mice. Further, we were able to trace neuronal projections across brain of thy1-GFP-M transgenic mice. The whole-brain high-resolution fluorescence imaging assured by CLSM may represent a powerful tool to navigate the brain through neuronal pathways. Although this work is focused on brain imaging, the macro-scale high-resolution tomographies affordable with CLSM are ideally suited to explore, at micron-scale resolution, the anatomy of different specimens like murine organs, embryos or flies.