Beyond multilayer perceptrons: Investigating complex topologies in neural networks.

This study delves into the crucial aspect of network topology in artificial neural networks (NNs) and its impact on model performance. Addressing the need to comprehend how network structures influence learning capabilities, the research contrasts traditional multilayer perceptrons (MLPs) with models built on various complex topologies using novel network generation techniques. Drawing insights from synthetic datasets, the study reveals the remarkable accuracy of complex NNs, particularly in high-difficulty scenarios, outperforming MLPs. Our exploration extends to real-world datasets, highlighting the task-specific nature of optimal network topologies and unveiling trade-offs, including increased computational demands and reduced robustness to graph damage in complex NNs compared to MLPs. This research underscores the pivotal role of complex topologies in addressing challenging learning tasks. However, it also signals the necessity for deeper insights into the complex interplay among topological attributes influencing NN performance. By shedding light on the advantages and limitations of complex topologies, this study provides valuable guidance for practitioners and paves the way for future endeavors to design more efficient and adaptable neural architectures across various applications.

Causal influence of brainstem response to transcutaneous vagus nerve stimulation on cardiovagal outflow.

BACKGROUND: The autonomic response to transcutaneous auricular vagus nerve stimulation (taVNS) has been linked to the engagement of brainstem circuitry modulating autonomic outflow. However, the physiological mechanisms supporting such efferent vagal responses are not well understood, particularly in humans. HYPOTHESIS: We present a paradigm for estimating directional brain-heart interactions in response to taVNS. We propose that our approach is able to identify causal links between the activity of brainstem nuclei involved in autonomic control and cardiovagal outflow. METHODS: We adopt an approach based on a recent reformulation of Granger causality that includes permutation-based, nonparametric statistics. The method is applied to ultrahigh field (7T) functional magnetic resonance imaging (fMRI) data collected on healthy subjects during taVNS. RESULTS: Our framework identified taVNS-evoked functional brainstem responses with superior sensitivity compared to prior conventional approaches, confirming causal links between taVNS stimulation and fMRI response in the nucleus tractus solitarii (NTS). Furthermore, our causal approach elucidated potential mechanisms by which information is relayed between brainstem nuclei and cardiovagal, i.e., high-frequency heart rate variability, in response to taVNS. Our findings revealed that key brainstem nuclei, known from animal models to be involved in cardiovascular control, exert a causal influence on taVNS-induced cardiovagal outflow in humans. CONCLUSION: Our causal approach allowed us to noninvasively evaluate directional interactions between fMRI BOLD signals from brainstem nuclei and cardiovagal outflow.

Heritability of human "directed" functional connectome.

INTRODUCTION: The functional connectivity patterns in the brain are highly heritable; however, it is unclear how genetic factors influence the directionality of such "information flows." Studying the "directionality" of the brain functional connectivity and assessing how heritability modulates it can improve our understanding of the human connectome. METHODS: Here, we investigated the heritability of "directed" functional connections using a state-space formulation of Granger causality (GC), in conjunction with blind deconvolution methods accounting for local variability in the hemodynamic response function. Such GC implementation is ideal to explore the directionality of functional interactions across a large number of networks. Resting-state functional magnetic resonance imaging data were drawn from the Human Connectome Project (total n = 898 participants). To add robustness to our findings, the dataset was randomly split into a "discovery" and a "replication" sample (each with n = 449 participants). The two cohorts were carefully matched in terms of demographic variables and other confounding factors (e.g., education). The effect of shared environment was also modeled. RESULTS: The parieto- and prefronto-cerebellar, parieto-prefrontal, and posterior-cingulate to hippocampus connections showed the highest and most replicable heritability effects with little influence by shared environment. In contrast, shared environmental factors significantly affected the visuo-parietal and sensory-motor directed connectivity. CONCLUSION: We suggest a robust role of heritability in influencing the directed connectivity of some cortico-subcortical circuits implicated in cognition. Further studies, for example using task-based fMRI and GC, are warranted to confirm the asymmetric effects of genetic factors on the functional connectivity within cognitive networks and their role in supporting executive functions and learning.

Can machine learning of post-procedural cone-beam CT images in acute ischemic stroke improve the detection of 24-h hemorrhagic transformation? A preliminary study.

PURPOSE: Hemorrhagic transformation (HT) is an independent predictor of unfavorable outcome in acute ischemic stroke (AIS) patients undergoing endovascular thrombectomy (EVT). Its early identification could help tailor AIS management. We hypothesize that machine learning (ML) applied to cone-beam computed tomography (CB-CT), immediately after EVT, improves performance in 24-h HT prediction. METHODS: We prospectively enrolled AIS patients undergoing EVT, post-procedural CB-CT, and 24-h non-contrast CT (NCCT). Three raters independently analyzed imaging at four anatomic levels qualitatively and quantitatively selecting a region of interest (ROI) < 5 mm2. Each ROI was labeled as "hemorrhagic" or "non-hemorrhagic" depending on 24-h NCCT. For each level of CB-CT, Mean Hounsfield Unit (HU), minimum HU, maximum HU, and signal- and contrast-to-noise ratios were calculated, and the differential HU-ROI value was compared between both hemispheres. The number of anatomic levels affected was computed for lesion volume estimation. ML with the best validation performance for 24-h HT prediction was selected. RESULTS: One hundred seventy-two ROIs from affected hemispheres of 43 patients were extracted. Ninety-two ROIs were classified as unremarkable, whereas 5 as parenchymal contrast staining, 29 as ischemia, 7 as subarachnoid hemorrhages, and 39 as HT. The Bernoulli Naïve Bayes was the best ML classifier with a good performance for 24-h HT prediction (sensitivity = 1.00; specificity = 0.75; accuracy = 0.82), though precision was 0.60. CONCLUSION: ML demonstrates high-sensitivity but low-accuracy 24-h HT prediction in AIS. The automated CB-CT imaging evaluation resizes sensitivity, specificity, and accuracy rates of visual interpretation reported in the literature so far. A standardized quantitative interpretation of CB-CT may be warranted to overcome the inter-operator variability.

Physically constrained neural networks for inferring physiological system models.

Systems biology and systems neurophysiology in particular have recently emerged as powerful tools for a number of key applications in the biomedical sciences. Nevertheless, such models are often based on complex combinations of multiscale (and possibly multiphysics) strategies that require ad hoc computational strategies and pose extremely high computational demands. Recent developments in the field of deep neural networks have demonstrated the possibility of formulating nonlinear, universal approximators to estimate solutions to highly nonlinear and complex problems with significant speed and accuracy advantages in comparison with traditional models. After synthetic data validation, we use so-called physically constrained neural networks (PINN) to simultaneously solve the biologically plausible Hodgkin-Huxley model and infer its parameters and hidden time-courses from real data under both variable and constant current stimulation, demonstrating extremely low variability across spikes and faithful signal reconstruction. The parameter ranges we obtain are also compatible with prior knowledge. We demonstrate that detailed biological knowledge can be provided to a neural network, making it able to fit complex dynamics over both simulated and real data.

Spatiotemporal learning of dynamic positron emission tomography data improves diagnostic accuracy in breast cancer.

Positron emission tomography (PET) can reveal metabolic activity in a voxelwise manner. PET analysis is commonly performed in a static manner by analyzing the standardized uptake value (SUV) obtained from the plateau region of PET acquisitions. A dynamic PET acquisition can provide a map of the spatiotemporal concentration of the tracer in vivo, hence conveying information about radiotracer delivery to tissue, its interaction with the target and washout. Therefore, tissue-specific biochemical properties are embedded in the shape of time activity curves (TACs), which are generally used for kinetic analysis. Conventionally, TACs are employed along with information about blood plasma activity concentration, i.e., the arterial input function (AIF), and specific compartmental models to obtain a full quantitative analysis of PET data. The main drawback of this approach is the need for invasive procedures requiring arterial blood sample collection during the whole PET scan. In this paper, we address the challenge of improving PET diagnostic accuracy through an alternative approach based on the analysis of time signal intensity patterns. Specifically, we demonstrate the diagnostic potential of tissue TACs provided by dynamic PET acquisition using various deep learning models. Our framework is shown to outperform the discriminative potential of classical SUV analysis, hence paving the way for more accurate PET-based lesion discrimination without additional acquisition time or invasive procedures. Clinical Relevance- The diagnostic accuracy of dynamic PET data exploited by deep-learning based time signal intensity pattern analysis is superior to that of static SUV imaging.

Patient-clinician brain concordance underlies causal dynamics in nonverbal communication and negative affective expressivity.

Patient-clinician concordance in behavior and brain activity has been proposed as a potential key mediator of mutual empathy and clinical rapport in the therapeutic encounter. However, the specific elements of patient-clinician communication that may support brain-to-brain concordance and therapeutic alliance are unknown. Here, we investigated how pain-related, directional facial communication between patients and clinicians is associated with brain-to-brain concordance. Patient-clinician dyads interacted in a pain-treatment context, during synchronous assessment of brain activity (fMRI hyperscanning) and online video transfer, enabling face-to-face social interaction. In-scanner videos were used for automated individual facial action unit (AU) time-series extraction. First, an interpretable machine-learning classifier of patients' facial expressions, from an independent fMRI experiment, significantly distinguished moderately painful leg pressure from innocuous pressure stimuli. Next, we estimated neural-network causality of patient-to-clinician directional information flow of facial expressions during clinician-initiated treatment of patients' evoked pain. We identified a leader-follower relationship in which patients predominantly led the facial communication while clinicians responded to patients' expressions. Finally, analyses of dynamic brain-to-brain concordance showed that patients' mid/posterior insular concordance with the clinicians' anterior insula cortex, a region identified in previously published data from this study1, was associated with therapeutic alliance, and self-reported and objective (patient-to-clinician-directed causal influence) markers of negative-affect expressivity. These results suggest a role of patient-clinician concordance of the insula, a social-mirroring and salience-processing brain node, in mediating directional dynamics of pain-directed facial communication during therapeutic encounters.

Classification of real-world pathological phonocardiograms through multi-instance learning.

Heart auscultation is an inexpensive and fundamental technique to effectively to diagnose cardiovascular disease. However, due to relatively high human error rates even when auscultation is performed by an experienced physician, and due to the not universal availability of qualified personnel e.g. in developing countries, a large body of research is attempting to develop automated, computational tools for detecting abnormalities in heart sounds. The large heterogeneity of achievable data quality and devices, the variety o possible heart pathologies, and a generally poor signal-to-noise ratio make this problem extremely challenging. We present an accurate classification strategy for diagnosing heart sounds based on 1) automatic heart phase segmentation, 2) state-of-the art filters drawn from the filed of speech synthesis (mel-frequency cepstral representation), and 3) an ad-hoc multi-branch, multi-instance artificial neural network based on convolutional layers and fully connected neuronal ensembles which separately learns from each heart phase, hence leveraging their different physiological significance. We demonstrate that it is possible to train our architecture to reach very high performances, e.g. an AUC of 0.87 or a sensitivity of 0.97. Our machine-learning-based tool could be employed for heart sound classification, especially as a screening tool in a variety of situations including telemedicine applications.

Echo state network models for nonlinear Granger causality.

While Granger causality (GC) has been often employed in network neuroscience, most GC applications are based on linear multivariate autoregressive (MVAR) models. However, real-life systems like biological networks exhibit notable nonlinear behaviour, hence undermining the validity of MVAR-based GC (MVAR-GC). Most nonlinear GC estimators only cater for additive nonlinearities or, alternatively, are based on recurrent neural networks or long short-term memory networks, which present considerable training difficulties and tailoring needs. We reformulate the GC framework in terms of echo-state networks-based models for arbitrarily complex networks, and characterize its ability to capture nonlinear causal relations in a network of noisy Duffing oscillators, showing a net advantage of echo state GC (ES-GC) in detecting nonlinear, causal links. We then explore the structure of ES-GC networks in the human brain employing functional MRI data from 1003 healthy subjects drawn from the human connectome project, demonstrating the existence of previously unknown directed within-brain interactions. In addition, we examine joint brain-heart signals in 15 subjects where we explore directed interaction between brain networks and central vagal cardiac control in order to investigate the so-called central autonomic network in a causal manner. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

A novel multi-branch architecture for state of the art robust detection of pathological phonocardiograms.

Heart auscultation is an inexpensive and fundamental technique to effectively diagnose cardiovascular disease. However, due to relatively high human error rates even when auscultation is performed by an experienced physician, and due to the not universal availability of qualified personnel, e.g. in developing countries, many efforts are made worldwide to propose computational tools for detecting abnormalities in heart sounds. The large heterogeneity of achievable data quality and devices, the variety of possible heart pathologies, and a generally poor signal-to-noise ratio make this problem very challenging. We present an accurate classification strategy for diagnosing heart sounds based on (1) automatic heart phase segmentation, (2) state-of-the art filters drawn from the field of speech synthesis (mel-frequency cepstral representation) and (3) an ad hoc multi-branch, multi-instance artificial neural network based on convolutional layers and fully connected neuronal ensembles which separately learns from each heart phase hence implicitly leveraging their different physiological significance. We demonstrate that it is possible to train our architecture to reach very high performances, e.g. an area under the curve of 0.87 or a sensitivity of 0.97. Our machine-learning-based tool could be employed for heartsound classification, especially as a screening tool in a variety of situations including telemedicine applications. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Deep computational pathology in breast cancer.

Deep Learning (DL) algorithms are a set of techniques that exploit large and/or complex real-world datasets for cross-domain and cross-discipline prediction and classification tasks. DL architectures excel in computer vision tasks, and in particular image processing and interpretation. This has prompted a wave of disruptingly innovative applications in medical imaging, where DL strategies have the potential to vastly outperform human experts. This is particularly relevant in the context of histopathology, where whole slide imaging (WSI) of stained tissue in conjuction with DL algorithms for their interpretation, selection and cancer staging are beginning to play an ever increasing role in supporting human operators in visual assessments. This has the potential to reduce everyday workload as well as to increase precision and reproducibility across observers, centers, staining techniques and even pathologies. In this paper we introduce the most common DL architectures used in image analysis, with a focus on histopathological image analysis in general and in breast histology in particular. We briefly review how, state-of-art DL architectures compare to human performance on across a number of critical tasks such as mitotic count, tubules analysis and nuclear pleomorphism analysis. Also, the development of DL algorithms specialized to pathology images have been enormously fueled by a number of world-wide challenges based on large, multicentric image databases which are now publicly available. In turn, this has allowed most recent efforts to shift more and more towards semi-supervised learning methods, which provide greater flexibility and applicability. We also review all major repositories of manually labelled pathology images in breast cancer and provide an in-depth discussion of the challenges specific to training DL architectures to interpret WSI data, as well as a review of the state-of-the-art methods for interpretation of images generated from immunohistochemical analysis of breast lesions. We finally discuss the future challenges and opportunities which the adoption of DL paradigms is most likely to pose in the field of pathology for breast cancer detection, diagnosis, staging and prognosis. This review is intended as a comprehensive stepping stone into the field of modern computational pathology for a transdisciplinary readership across technical and medical disciplines.

Thalamic neuroinflammation as a reproducible and discriminating signature for chronic low back pain.

ABSTRACT: Using positron emission tomography, we recently demonstrated elevated brain levels of the 18 kDa translocator protein (TSPO), a glial activation marker, in chronic low back pain (cLBP) patients, compared to healthy controls (HCs). Here, we first sought to replicate the original findings in an independent cohort (15 cLBP, 37.8 ± 12.5 y/o; 18 HC, 48.2 ± 12.8 y/o). We then trained random forest machine learning algorithms based on TSPO imaging features combining discovery and replication cohorts (totaling 25 cLBP, 42.4 ± 13.2 y/o; 27 HC, 48.9 ± 12.6 y/o), to explore whether image features other than the mean contain meaningful information that might contribute to the discrimination of cLBP patients and HC. Feature importance was ranked using SHapley Additive exPlanations values, and the classification performance (in terms of area under the curve values) of classifiers containing only the mean, other features, or all features was compared using the DeLong test. Both region-of-interest and voxelwise analyses replicated the original observation of thalamic TSPO signal elevations in cLBP patients compared to HC (P < 0.05). The random forest-based analyses revealed that although the mean is a discriminating feature, other features demonstrate similar level of importance, including the maximum, kurtosis, and entropy. Our observations suggest that thalamic neuroinflammatory signal is a reproducible and discriminating feature for cLBP, further supporting a role for glial activation in human cLBP, and the exploration of neuroinflammation as a therapeutic target for chronic pain. This work further shows that TSPO signal contains a richness of information that the simple mean might fail to capture completely.

Radiomics in breast cancer classification and prediction.

Breast Cancer (BC) is the common form of cancer in women. Its diagnosis and screening are usually performed through different imaging modalities such as mammography, magnetic resonance imaging and ultrasound. However, mammography and ultrasound-imaging techniques have limited sensitivity and specificity both in identifying lesions and in differentiating malign from benign lesions, especially in presence of dense breast parenchyma. Due to the higher resolution of magnetic resonance images, MRI represents the method with the higher specificity and sensitivity among all the available tools, in both lesions' identification and diagnosis. However, especially for diagnosis, even MRI has limitations that are only partially solved if combined with mammography. Unfortunately, due to the limits of all these imaging tools, in order to have a certain diagnosis, patients often receive painful and costly bioptics procedures. In this context, several computational approaches have been developed to increase sensitivity, while maintaining the same specificity, in BC diagnosis and screening. Amongst these, radiomics has been increasingly gaining ground in oncology to improve cancer diagnosis, prognosis and treatment. Radiomics derives multiple quantitative features from single or multiple medical imaging modalities, highlighting image traits which are not visible to the naked eye and hence significantly augmenting the discriminatory and predictive potential of medical imaging. This review article aims to summarize the state of the art in radiomics-based BC research. The dominating evidence extracted from the literature points towards a high potential of radiomics in disentangling malignant from benign breast lesions, classifying BC types and grades and also in predicting treatment response and recurrence risk. In the era of personalized medicine, radiomics has the potential to improve diagnosis, prognosis, prediction, monitoring, image-based intervention, and assessment of therapeutic response in BC.

Serological determinants of COVID-19.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection spreaded rapidly worldwide, as far as it has become a global pandemic. Therefore, the introduction of serological tests for determination of IgM and IgG antibodies has become the main diagnostic tool, useful for tracking the spread of the virus and for consequently allowing its containment. In our study we compared point of care test (POCT) lateral flow immunoassay (FIA) vs automated chemiluminescent immunoassay (CLIA), in order to assess their specificity and sensibility for COVID-19 antibodies detection. RESULTS: We find that different specificities and sensitivities for IgM and IgG tests. Notably IgM POCT FIA method vs CLIA method (gold standard) has a low sensitivity (0.526), while IgG POCT FIA method vs CLIA method (gold standard) test has a much higher sensitivity (0.937); further, with respect of IgG, FIA and CLIA could arguably provide equivalent information. CONCLUSIONS: FIA method could be helpful in assessing in short time, the possible contagiousness of subjects that for work reasons cannot guarantee "social distancing".

Unsupervised stratification in neuroimaging through deep latent embeddings.

There is growing evidence that the use of stringent and dichotomic diagnostic categories in many medical disciplines (particularly 'brain sciences' as neurology and psychiatry) is an oversimplification. Although clear diagnostic boundaries remain useful for patients, families, and their access to dedicated NHS and health care services, the traditional dichotomic categories are not helpful to describe the complexity and large heterogeneity of symptoms across many and overlapping clinical phenotypes. With the advent of 'big' multimodal neuroimaging databases, data-driven stratification of the wide spectrum of healthy human physiology or disease based on neuroimages is theoretically become possible. However, this conceptual framework is hampered by severe computational constraints. In this paper we present a novel, deep learning based encode-decode architecture which leverages several parameter efficiency techniques generate latent deep embedding which compress the information contained in a full 3D neuroimaging volume by a factor 1000 while still retaining anatomical detail and hence rendering the subsequent stratification problem tractable. We train our architecture on 1003 brain scan derived from the human connectome project and demonstrate the faithfulness of the obtained reconstructions. Further, we employ a data driven clustering technique driven by a grid search in hyperparameter space to identify six different strata within the 1003 healthy community dwelling individuals which turn out to correspond to highly significant group differences in both physiological and cognitive data. Indicating that the well-known relationships between such variables and brain structure can be probed in an unsupervised manner through our novel architecture and pipeline. This opens the door to a variety of previously inaccessible applications in the realm of data driven stratification of large cohorts based on neuroimaging data.Clinical Relevance -With our approach, each person can be described and classified within a multi-dimensional space of data, where they are uniquely classified according to their individual anatomy, physiology and disease-related anatomical and physiological alterations.

Frequency dependent functional brain reorganization in anesthesia is specific to drug concentration.

The differential effects of general anesthesia on brain activity in terms of drug selection, concentration and combination remain to be elucidated. Using fMRI, it has been shown that increasing doses of sevoflurane is associated with progressive breakdown in brain functional connectivity, while EEG studies have shown that higher activity in the delta band is associated with unconsciousness. Despite these promising results, the band- specific neural substrates of brain changes which occur during sevoflurane anesthesia have not yet been investigated. To this end, we employ high-density EEG-based brain connectivity estimates and graph theoretical analysis in a protocol of progressive sevoflurane administration (conditions: baseline, 1.1%, 2.1%, 2.8%, recovery), both at a global (whole-brain) and at a local (sensor-specific) level in 12 healthy subjects (7 males, mean age 25 ± 4.7 years). We show a statistically significant dependence of global strength, clustering coefficient and efficiency on sevoflurane concentration in the slow delta, beta 1 and beta 2 bands. Interestingly, high and low-frequency bands behaved in an opposite manner as a function of condition. We also found significant band*condition interactive effects in clustering coefficient, efficiency and strength both on local and global scales.

Uncovering complex central autonomic networks at rest: a functional magnetic resonance imaging study on complex cardiovascular oscillations.

This study aims to uncover brain areas that are functionally linked to complex cardiovascular oscillations in resting-state conditions. Multi-session functional magnetic resonance imaging (fMRI) and cardiovascular data were gathered from 34 healthy volunteers recruited within the human connectome project (the '100-unrelated subjects' release). Group-wise multi-level fMRI analyses in conjunction with complex instantaneous heartbeat correlates (entropy and Lyapunov exponent) revealed the existence of a specialized brain network, i.e. a complex central autonomic network (CCAN), reflecting what we refer to as complex autonomic control of the heart. Our results reveal CCAN areas comprised the paracingulate and cingulate gyri, temporal gyrus, frontal orbital cortex, planum temporale, temporal fusiform, superior and middle frontal gyri, lateral occipital cortex, angular gyrus, precuneous cortex, frontal pole, intracalcarine and supracalcarine cortices, parahippocampal gyrus and left hippocampus. The CCAN visible at rest does not include the insular cortex, thalamus, putamen, amygdala and right caudate, which are classical CAN regions peculiar to sympatho-vagal control. Our results also suggest that the CCAN is mainly involved in complex vagal control mechanisms, with possible links with emotional processing networks.

A novel approach for geographical risk mapping of morbidity and mortality rates: the case of Val D'Agri, Italy.

While associations between exposure to air pollutants and increased morbidity and mortality are well established, few rigorous studies on this issue are available. The aim of the current study is to implement a new approach to the spatial analysis of mortality and morbidity, based on testing for the presence of the same association in other areas of similar size. Additionally, we perform a case study in Val d'Agri (VA), an area of Basilicata Region, Southern Italy, where oil and natural gas extraction began in 1998. In order to examine the spatial distribution of morbidity and mortality in the region of interest, Hospital discharge (2001-2013) and mortality (2003-2014) rates for the main environment-related diseases were calculated. In addition, a comparison between the period 1980-1998 and the period 1999-2014 was performed for cardiovascular disease mortality. For the period under study, a neutral scenario emerged for cancer and respiratory diseases, where we found no differences in morbidity and mortality as compared to the national benchmark. In some cases significantly lower values (as compared to the nation-wide benchmark) were found. Conversely, a slight excess in morbidity and mortality (as compared to the nation-wide benchmark) emerged for cardiovascular diseases. Still, this excess was common to a number of municipalities in the surroundings of VA, and appeared to be already present in 1980. Higher rates of cardiovascular diseases, lower rates of neoplastic disorders no differences in mortality for respiratory causes (as compared to the nation-wide benchmark) were found in multiple areas of the region, and were therefore not specific to VA. In summary, our data do not support the hypothesis of a role of industrial activities related to oil extraction in VA in determining mortality and morbidity patterns and trends.

Novel insights into breast cancer progression and metastasis: A multidisciplinary opportunity to transition from biology to clinical oncology.

According to the most recent epidemiological studies, breast cancer shows the highest incidence and the second leading cause of death in women. Cancer progression and metastasis are the main events related to poor survival of breast cancer patients. This can be explained by the presence of highly resistant to chemo- and radiotherapy stem cells in many breast tumor tissues. In this context, numerous studies highlighted the possible involvement of epithelial to mesenchymal transition phenomenon as biological program to generate cancer stem cells, and thus participate to both metastatic and drug resistance process. Therefore, the comprehension of mechanisms (both cellular and molecular) involved in breast cancer occurrence and progression can lay the foundation for the development of new diagnostic and therapeutical protocols. In this review, we reported the most important findings in the field of breast cancer highlighting the most recent data concerning breast tumor biology, diagnosis and therapy.