Identification of new biomarkers for breast cancer remains critical in order to enhance early detection of the disease and improve its prognosis. Towards this end, we performed an untargeted metabolomic analysis of breast ductal fluid using an ultra-performance liquid chromatography coupled with a quadrupole time-of-light (UPLC-QTOF) mass spectrometer. We investigated the metabolomic profiles of breast tumors using ductal fluid samples collected by ductal lavage (DL). We studied fluid from both the affected breasts and the unaffected contralateral breasts (as controls) from 43 women with confirmed unilateral breast cancer. Using this approach, we identified 1560 ions in the positive mode and 538 ions in the negative mode after preprocessing of the UPLC‑QTOF data. Paired t-tests applied on these data matrices identified 209 ions (positive and negative modes combined) with significant change in intensity level between affected and unaffected control breasts (adjusted p-values <0.05). Among these, 83 ions (39.7%) showed a fold change (FC) >1.2 and 66 ions (31.6%) were identified with putative compound names. The metabolites that we identified included endogenous metabolites such as amino acid derivatives (N-Acetyl-DL-tryptophan) or products of lipid metabolism such as N-linoleoyl taurine, trans-2-dodecenoylcarnitine, lysophosphatidylcholine LysoPC(18:2(9Z,12Z)), glycerophospholipids PG(18:0/0:0), and phosphatidylserine PS(20:4(5Z,8Z,11Z,14Z). Generalized LASSO regression further selected 21 metabolites when race, menopausal status, smoking, grade and TNM stage were adjusted for. A predictive conditional logistic regression model, using the LASSO selected 21 ions, provided diagnostic accuracy with the area under the curve of 0.956 (sensitivity/specificity of 0.907/0.884). This is the first study that shows the feasibility of conducting a comprehensive metabolomic profiling of breast tumors using breast ductal fluid to detect changes in the cellular microenvironment of the tumors and shows the potential for this approach to be used to improve detection of breast cancer.

Merlin/NF2 is a bona fide tumor suppressor whose mutations underlie inherited tumor syndrome neurofibromatosis type 2 (NF2), as well as various sporadic cancers including kidney cancer. Multiple Merlin/NF2 effector pathways including the Hippo-YAP/TAZ pathway have been identified. However, the molecular mechanisms underpinning the growth and survival of NF2-mutant tumors remain poorly understood. Using an inducible orthotopic kidney tumor model, we demonstrate that YAP/TAZ silencing is sufficient to induce regression of pre-established NF2-deficient tumors. Mechanistically, YAP/TAZ depletion diminishes glycolysis-dependent growth and increases mitochondrial respiration and reactive oxygen species (ROS) buildup, resulting in oxidative-stress-induced cell death when challenged by nutrient stress. Furthermore, we identify lysosome-mediated cAMP-PKA/EPAC-dependent activation of RAF-MEK-ERK signaling as a resistance mechanism to YAP/TAZ inhibition. Finally, unbiased analysis of TCGA primary kidney tumor transcriptomes confirms a positive correlation of a YAP/TAZ signature with glycolysis and inverse correlations with oxidative phosphorylation and lysosomal gene expression, supporting the clinical relevance of our findings.

Pituitary null cell adenoma is a challenging clinical condition, and its pathogenesis remains to be elucidated. We performed this study to determine the roles of C5orf66-AS1, NORAD, and TINCR in the pathogenesis and invasion of pituitary null cell adenomas. Expression of the three long non-coding RNAs in pituitary null cell adenoma tissues of 11 patients and normal pituitary tissues from four donors was examined by performing quantitative reverse transcription-polymerase chain reaction. We found that C5orf66-AS1 expression was lower in pituitary null cell adenoma tissues than in normal pituitary tissues. Moreover, C5orf66-AS1 expression level was significantly lower in invasive pituitary null cell adenomas than in non-invasive ones. After transfection of C5orf66-AS1 into pituitary adenoma cells, assessment of cell viability and invasion suggested that overexpressed C5orf66-AS1 inhibited cell viability and cell invasion. In silico algorithms predicted several cis- and trans-acting target genes of C5orf66-AS1, including PITX1 and SCGB3A1. In addition, expression of some of the predicted target genes was determined using microarray data of another cohort with pituitary null cell adenomas. It showed that some of these target genes were differentially expressed between pituitary null cell adenoma tissues and normal pituitary tissues as well as between invasive and non-invasive tumors. Co-expression analysis in RNA sequencing data showed that PAQR7 was the most correlated gene of C5orf66-AS1 and that several predicted trans-acting target genes, including SCGB3A1, were highly correlated with C5orf66-AS1. NORAD and TINCR expression was not statistically significant in the complete cohort; however, a negative correlation was observed between NORAD expression and maximum tumor diameter in some subgroups. These results indicate that C5orf66-AS1 suppresses the development and invasion of pituitary null cell adenomas. However, our results do not provide enough statistical evidence to support the roles of NORAD and TINCR in the development and invasion of pituitary null cell adenomas.

Sex is a key modifier of neurological disease outcomes. Microglia are implicated in neurological diseases and modulated by microRNAs, but it is unknown whether microglial microRNAs have sex-specific influences on disease. We show in mice that microglial microRNA expression differs in males and females and that loss of microRNAs leads to sex-specific changes in the microglial transcriptome and tau pathology. These findings suggest that microglial microRNAs influence tau pathogenesis in a sex-specific manner.

The identity and degree of heterogeneity of glial progenitors and their contributions to brain tumor malignancy remain elusive. By applying lineage-targeted single-cell transcriptomics, we uncover an unanticipated diversity of glial progenitor pools with unique molecular identities in developing brain. Our analysis identifies distinct transitional intermediate states and their divergent developmental trajectories in astroglial and oligodendroglial lineages. Moreover, intersectional analysis uncovers analogous intermediate progenitors during brain tumorigenesis, wherein oligodendrocyte-progenitor intermediates are abundant, hyper-proliferative, and progressively reprogrammed toward a stem-like state susceptible to further malignant transformation. Similar actively cycling intermediate progenitors are prominent components in human gliomas with distinct driver mutations. We further unveil lineage-driving networks underlying glial fate specification and identify Zfp36l1 as necessary for oligodendrocyte-astrocyte lineage transition and glioma growth. Together, our results resolve the dynamic repertoire of common and divergent glial progenitors during development and tumorigenesis and highlight Zfp36l1 as a molecular nexus for balancing glial cell-fate decision and controlling gliomagenesis.

Cortical state, defined by population-level neuronal activity patterns, determines sensory perception. While arousal-associated neuromodulators-including norepinephrine (NE)-reduce cortical synchrony, how the cortex resynchronizes remains unknown. Furthermore, general mechanisms regulating cortical synchrony in the wake state are poorly understood. Using in vivo imaging and electrophysiology in mouse visual cortex, we describe a critical role for cortical astrocytes in circuit resynchronization. We characterize astrocytes' calcium responses to changes in behavioral arousal and NE, and show that astrocytes signal when arousal-driven neuronal activity is reduced and bi-hemispheric cortical synchrony is increased. Using in vivo pharmacology, we uncover a paradoxical, synchronizing response to Adra1a receptor stimulation. We reconcile these results by demonstrating that astrocyte-specific deletion of Adra1a enhances arousal-driven neuronal activity, while impairing arousal-related cortical synchrony. Our findings demonstrate that astrocytic NE signaling acts as a distinct neuromodulatory pathway, regulating cortical state and linking arousal-associated desynchrony to cortical circuit resynchronization.

Intermittent hypoxemia (IH) may influence retinopathy of prematurity (ROP) development in preterm infants, however, previous studies had mixed results. This study aims to assess the influence and evaluate the predictive ability of IH measures on Type 1 ROP, a stage beyond which ROP treatment is indicated.

Precise regulation of the cell cycle is crucial to the growth and development of all organisms. Understanding the regulatory mechanism of the cell cycle is crucial to unraveling many complicated diseases, most notably cancer. Multiple sources of biological data are available to study the dynamic interactions among many genes that are related to the cancer cell cycle. Integrating these informative and complementary data sources can help to infer a mutually consistent gene transcriptional regulatory network with strong similarity to the underlying gene regulatory relationships in cancer cells.

Metabolomics is likely an ideal tool to assess tobacco smoke exposure and the impact of cigarette smoke on human exposure and health. To assess reproducibility and feasibility of this by UPLC-QTOF-MS, three experiments were designed for the assessment of smokers' blood. Experiment I was an analysis of 8 smokers with 8 replicates. Experiment II was an analysis of 62 pooled quality control (QC) samples from 7 nonsmokers' plasma placed as every tenth sample among a study of 613 samples from 160 smokers. Finally, to examine the feasibility of metabolomic study in assessing smoke exposure, Experiment III consisted of 9 smokers and 10 nonsmokers' serum to evaluate differences in their global metabolome. There was minimal measurement and sample preparation variation in all experiments, although some caution is needed when analyzing specific parts of the chromatogram. When assessing QC samples in the large scale study, QC clustering indicated high stability, reproducibility, and consistency. Finally, in addition to the identification of nicotine metabolites as expected, there was a characteristic profile distinguishing smokers from nonsmokers. Metabolites selected from putative identifications were verified by MS/MS, showing the potential to identify metabolic phenotypes and new metabolites relating to cigarette smoke exposure and toxicity.

Urine is increasingly being considered as a source of biomarker development in Duchenne Muscular Dystrophy (DMD), a severe, life-limiting disorder that affects approximately 1 in 4500 boys. In this study, we considered the mdx mice-a murine model of DMD-to discover biomarkers of disease, as well as pharmacodynamic biomarkers responsive to prednisolone, a corticosteroid commonly used to treat DMD. Longitudinal urine samples were analyzed from male age-matched mdx and wild-type mice randomized to prednisolone or vehicle control via liquid chromatography tandem mass spectrometry. A large number of metabolites (869 out of 6,334) were found to be significantly different between mdx and wild-type mice at baseline (Bonferroni-adjusted p-value < 0.05), thus being associated with disease status. These included a metabolite with m/z = 357 and creatine, which were also reported in a previous human study looking at serum. Novel observations in this study included peaks identified as biliverdin and hypusine. These four metabolites were significantly higher at baseline in the urine of mdx mice compared to wild-type, and significantly changed their levels over time after baseline. Creatine and biliverdin levels were also different between treated and control groups, but for creatine this may have been driven by an imbalance at baseline. In conclusion, our study reports a number of biomarkers, both known and novel, which may be related to either the mechanisms of muscle injury in DMD or prednisolone treatment.

Recent work examining astrocytic physiology centers on fluorescence imaging, due to development of sensitive fluorescent indicators and observation of spatiotemporally complex calcium activity. However, the field remains hindered in characterizing these dynamics, both within single cells and at the population level, because of the insufficiency of current region-of-interest-based approaches to describe activity that is often spatially unfixed, size-varying and propagative. Here we present an analytical framework that releases astrocyte biologists from region-of-interest-based tools. The Astrocyte Quantitative Analysis (AQuA) software takes an event-based perspective to model and accurately quantify complex calcium and neurotransmitter activity in fluorescence imaging datasets. We apply AQuA to a range of ex vivo and in vivo imaging data and use physiologically relevant parameters to comprehensively describe the data. Since AQuA is data-driven and based on machine learning principles, it can be applied across model organisms, fluorescent indicators, experimental modes, and imaging resolutions and speeds, enabling researchers to elucidate fundamental neural physiology.

Intramyocellular (IMCL), extramyocellular lipid (EMCL), and vitamin D deficiency are associated with muscle metabolic dysfunction. This study compared the change in [IMCL]:[EMCL] following the combined treatment of vitamin D and aerobic training (DAT) compared with vitamin D (D), aerobic training (AT), and control (CTL). Male and female subjects aged 60-80 years with a BMI ranging from 18.5-34.9 and vitamin D status of ≤32 ng/mL (25(OH)D) were recruited to randomized, prospective clinical trial double-blinded for supplement with a 2 × 2 factorial design. Cholecalciferol (Vitamin D3) (10,000 IU × 5 days/week) or placebo was provided for 13 weeks and treadmill aerobic training during week 13. Gastrocnemius IMCL and EMCL were measured with magnetic resonance spectroscopy (MRS) and MRI. Hybrid near-infrared diffuse correlation spectroscopy measured hemodynamics. Group differences in IMCL were observed when controlling for baseline IMCL (p = 0.049). DAT was the only group to reduce IMCL from baseline, while a mean increase was observed in all other groups combined (p = 0.008). IMCL reduction and the corresponding increase in rVO2 at study end (p = 0.011) were unique to DAT. Vitamin D, when combined with exercise, may potentiate the metabolic benefits of exercise by reducing IMCL and increasing tissue-level VO2 in healthy, older adults.

Microglia are increasingly recognized for their major contributions during brain development and neurodegenerative disease. It is currently unknown whether these functions are carried out by subsets of microglia during different stages of development and adulthood or within specific brain regions. Here, we performed deep single-cell RNA sequencing (scRNA-seq) of microglia and related myeloid cells sorted from various regions of embryonic, early postnatal, and adult mouse brains. We found that the majority of adult microglia expressing homeostatic genes are remarkably similar in transcriptomes, regardless of brain region. By contrast, early postnatal microglia are more heterogeneous. We discovered a proliferative-region-associated microglia (PAM) subset, mainly found in developing white matter, that shares a characteristic gene signature with degenerative disease-associated microglia (DAM). Such PAM have amoeboid morphology, are metabolically active, and phagocytose newly formed oligodendrocytes. This scRNA-seq atlas will be a valuable resource for dissecting innate immune functions in health and disease.

We investigated the influence of genetic variants on atherosclerosis using whole exome sequencing in cases and controls from the autopsy study "Pathobiological Determinants of Atherosclerosis in Youth (PDAY)". We identified a PDAY case group with the highest total amounts of raised lesions (n = 359) for comparisons with a control group with no detectable raised lesions (n = 626). In addition to the standard exome capture, we included genome-wide proximal promoter regions that contain sequences that regulate gene expression. Our statistical analyses included single variant analysis for common variants (MAF > 0.01) and rare variant analysis for low frequency and rare variants (MAF < 0.05). In addition, we investigated known CAD genes previously identified by meta-analysis of GWAS studies. We did not identify individual common variants that reached exome-wide significance using single variant analysis. In analysis limited to 60 CAD genes, we detected strong associations with COL4A2/COL4A1 that also previously showed associations with myocardial infarction and arterial stiffness, as well as coronary artery calcification. Likewise, rare variant analysis did not identify genes that reached exome-wide significance. Among the 60 CAD genes, the strongest association was with NBEAL1 that was also identified in gene-based analysis of whole exome sequencing for early onset myocardial infarction.

Continuous and longitudinal imaging of cerebral blood flow (CBF) variations provide vital information to investigate pathophysiology and interventions for a variety of neurological and cerebral diseases. An innovative noncontact speckle contrast diffuse correlation tomography (scDCT) system was downscaled and adapted for noninvasive imaging of CBF distributions in rat brain through intact scalp and skull. Algorithms for 2D mapping and 3D image reconstruction of CBF distributions were developed and optimized. The continuous imaging capability of the system was shown by imaging global CBF increases during CO2 inhalations and regional CBF decreases across two hemispheres during sequential unilateral and bilateral common carotid artery ligations. The longitudinal imaging capability was demonstrated by imaging CBF variations over a long recovery period of 14 days after an acute stroke. Compared to the 2D mapping method, the 3D imaging method reduces partial volume effects, but needs more computation time for image reconstruction. Results from this study generally agree with those reported in the literature using similar protocols to induce CBF changes in rats. The scDCT enables a relatively large penetration depth (up to ∼10 mm), which is sufficient for transcranial brain measurements in small animals and human neonates. Ultimately, we expect to provide a noninvasive noncontact cerebral imager for basic neuroscience research in small animal models and clinical applications in human neonates.

Recent discoveries that astrocytes exert proactive regulatory effects on neural information processing and that they are deeply involved in normal brain development and disease pathology have stimulated broad interest in understanding astrocyte functional roles in brain circuit. Measuring astrocyte functional status is now technically feasible, due to recent advances in modern microscopy and ultrasensitive cell-type specific genetically encoded Ca2+ indicators for chronic imaging. However, there is a big gap between the capability of generating large dataset via calcium imaging and the availability of sophisticated analytical tools for decoding the astrocyte function. Current practice is essentially manual, which not only limits analysis throughput but also risks introducing bias and missing important information latent in complex, dynamic big data. Here, we report a suite of computational tools, called Functional AStrocyte Phenotyping (FASP), for automatically quantifying the functional status of astrocytes. Considering the complex nature of Ca2+ signaling in astrocytes and low signal to noise ratio, FASP is designed with data-driven and probabilistic principles, to flexibly account for various patterns and to perform robustly with noisy data. In particular, FASP explicitly models signal propagation, which rules out the applicability of tools designed for other types of data. We demonstrate the effectiveness of FASP using extensive synthetic and real data sets. The findings by FASP were verified by manual inspection. FASP also detected signals that were missed by purely manual analysis but could be confirmed by more careful manual examination under the guidance of automatic analysis. All algorithms and the analysis pipeline are packaged into a plugin for Fiji (ImageJ), with the source code freely available online at https://github.com/VTcbil/FASP.

Conventional differential gene expression analysis by methods such as student's t-test, SAM, and Empirical Bayes often searches for statistically significant genes without considering the interactions among them. Network-based approaches provide a natural way to study these interactions and to investigate the rewiring interactions in disease versus control groups. In this paper, we apply weighted graphical LASSO (wgLASSO) algorithm to integrate a data-driven network model with prior biological knowledge (i.e., protein-protein interactions) for biological network inference. We propose a novel differentially weighted graphical LASSO (dwgLASSO) algorithm that builds group-specific networks and perform network-based differential gene expression analysis to select biomarker candidates by considering their topological differences between the groups.

The prefrontal cortex (PFC) is a hub for cognitive control, and dopamine profoundly influences its functions. In other brain regions, astrocytes sense diverse neurotransmitters and neuromodulators and, in turn, orchestrate regulation of neuroactive substances. However, basic physiology of PFC astrocytes, including which neuromodulatory signals they respond to and how they contribute to PFC function, is unclear. Here, we characterize divergent signaling signatures in mouse astrocytes of the PFC and primary sensory cortex, which show differential responsiveness to locomotion. We find that PFC astrocytes express receptors for dopamine but are unresponsive through the Gs/Gi-cAMP pathway. Instead, fast calcium signals in PFC astrocytes are time locked to dopamine release and are mediated by α1-adrenergic receptors both ex vivo and in vivo. Further, we describe dopamine-triggered regulation of extracellular ATP at PFC astrocyte territories. Thus, we identify astrocytes as active players in dopaminergic signaling in the PFC, contributing to PFC function though neuromodulator receptor crosstalk.

Pathologic alterations in epigenetic regulation have long been considered a hallmark of many cancers, including hepatocellular carcinoma (HCC). In a healthy individual, the relationship between DNA methylation and microRNA (miRNA) expression maintains a fine balance; however, disruptions in this harmony can aid in the genesis of cancer or the propagation of existing cancers. The balance between DNA methylation and microRNA expression and its potential disturbance in HCC can vary by race. There is emerging evidence linking epigenetic events including DNA methylation and miRNA expression to cancer disparities. In this paper, we evaluate the epigenetic mechanisms of racial heterogenity in HCC through an integrated analysis of DNA methylation, miRNA, and combined regulation of gene expression. Specifically, we generated DNA methylation, mRNA-seq, and miRNA-seq data through the analysis of tumor and adjacent non-tumor liver tissues from African Americans (AA) and European Americans (EA) with HCC. Using mixed ANOVA, we identified cytosine-phosphate-guanine (CpG) sites, mRNAs, and miRNAs that are significantly altered in HCC vs. adjacent non-tumor tissue in a race-specific manner. We observed that the methylome was drastically changed in EA with a significantly larger number of differentially methylated and differentially expressed genes than in AA. On the other hand, the miRNA expression was altered to a larger extent in AA than in EA. Pathway analysis functionally linked epigenetic regulation in EA to processes involved in immune cell maturation, inflammation, and vascular remodeling. In contrast, cellular proliferation, metabolism, and growth pathways are found to predominate in AA as a result of this epigenetic analysis. Furthermore, through integrative analysis, we identified significantly differentially expressed genes in HCC with disparate epigenetic regulation, associated with changes in miRNA expression for AA and DNA methylation for EA.

Recent single-cell RNA sequencing studies have revealed distinct microglial states in development and disease. These include proliferative region-associated microglia (PAM) in developing white matter and disease-associated microglia (DAM) prevalent in various neurodegenerative conditions. PAM and DAM share a similar core gene signature and other functional properties. However, the extent of the dynamism and plasticity of these microglial states, as well as their functional significance, remains elusive, partly due to the lack of specific tools. Here, we report the generation of an inducible Cre driver line, Clec7a-CreERT2, designed to target PAM and DAM in the brain parenchyma. Utilizing this tool, we profile labeled cells during development and in several disease models, uncovering convergence and context-dependent differences in PAM/DAM gene expression. Through long-term tracking, we demonstrate surprising levels of plasticity in these microglial states. Lastly, we specifically depleted DAM in cuprizone-induced demyelination, revealing their roles in disease progression and recovery.