The exosome and its nuclear specific subunit Rrp6 form a 3'-5' exonuclease complex that regulates diverse aspects of RNA biology including 3' end processing and degradation of a variety of noncoding RNAs (ncRNAs) and unstable transcripts. Known targets of the nuclear exosome include short (<1000 bp) RNAPII transcripts such as small noncoding RNAs (snRNAs), cryptic unstable transcripts (CUTs), and some stable unannotated transcripts (SUTs) that are terminated by an Nrd1, Nab3, and Sen1 (NNS) dependent mechanism. NNS-dependent termination is coupled to RNA 3' end processing and/or degradation by the Rrp6/exosome in yeast. Recent work suggests Nrd1 is necessary for transcriptome surveillance, regulating promoter directionality and suppressing antisense transcription independently of, or prior to, Rrp6 activity. It remains unclear whether Rrp6 is directly involved in termination; however, Rrp6 has been implicated in the 3' end processing and degradation of ncRNA transcripts including CUTs. To determine the role of Rrp6 in NNS termination globally, we performed RNA sequencing (RNA-Seq) on total RNA and perform ChIP-exo analysis of RNA Polymerase II (RNAPII) localization. Deletion of RRP6 promotes hyper-elongation of multiple NNS-dependent transcripts resulting from both improperly processed 3' RNA ends and faulty transcript termination at specific target genes. The defects in RNAPII termination cause transcriptome-wide changes in mRNA expression through transcription interference and/or antisense repression, similar to previously reported effects of depleting Nrd1 from the nucleus. Elongated transcripts were identified within all classes of known NNS targets with the largest changes in transcription termination occurring at CUTs. Interestingly, the extended transcripts that we have detected in our studies show remarkable similarity to Nrd1-unterminated transcripts at many locations, suggesting that Rrp6 acts with the NNS complex globally to promote transcription termination in addition to 3' end RNA processing and/or degradation at specific targets.

Epithelial to mesenchymal transition (EMT) can contribute to gastric cancer (GC) progression and recurrence following therapy. Tumor-associated neutrophils (TANs) are associated with poor outcomes in a variety of cancers. However, it is not clear whether TANs interact with the EMT process during GC development.

Genome-wide association studies (GWAS) of complex traits, such as alcohol use disorders (AUD), usually identify variants in non-coding regions and cannot by themselves distinguish whether the associated variants are functional or in linkage disequilibrium with the functional variants. Transcriptome studies can identify genes whose expression differs between alcoholics and controls. To test which variants associated with AUD may cause expression differences, we integrated data from deep RNA-seq and GWAS of four postmortem brain regions from 30 subjects with AUD and 30 controls to analyze allele-specific expression (ASE). We identified 88 genes with differential ASE in subjects with AUD compared to controls. Next, to test one potential mechanism contributing to the differential ASE, we analyzed single nucleotide polymorphisms (SNPs) in the 3' untranslated regions (3'UTR) of these genes. Of the 88 genes with differential ASE, 61 genes contained 437 SNPs in the 3'UTR with at least one heterozygote among the subjects studied. Using a modified PASSPORT-seq (parallel assessment of polymorphisms in miRNA target-sites by sequencing) assay, we identified 25 SNPs that affected RNA levels in a consistent manner in two neuroblastoma cell lines, SH-SY5Y and SK-N-BE(2). Many of these SNPs are in binding sites of miRNAs and RNA-binding proteins, indicating that these SNPs are likely causal variants of AUD-associated differential ASE. In sum, we demonstrate that a combination of computational and experimental approaches provides a powerful strategy to uncover functionally relevant variants associated with the risk for AUD.

Brain metastasis (br-met) develops in an immunologically unique br-met niche. Central nervous system-native myeloid cells (CNS-myeloids) and bone-marrow-derived myeloid cells (BMDMs) cooperatively regulate brain immunity. The phenotypic heterogeneity and specific roles of these myeloid subsets in shaping the br-met niche to regulate br-met outgrowth have not been fully revealed. Applying multimodal single-cell analyses, we elucidated a heterogeneous but spatially defined CNS-myeloid response during br-met outgrowth. We found Ccr2+ BMDMs minimally influenced br-met while CNS-myeloid promoted br-met outgrowth. Additionally, br-met-associated CNS-myeloid exhibited downregulation of Cx3cr1. Cx3cr1 knockout in CNS-myeloid increased br-met incidence, leading to an enriched interferon response signature and Cxcl10 upregulation. Significantly, neutralization of Cxcl10 reduced br-met, while rCxcl10 increased br-met and recruited VISTAHi PD-L1+ CNS-myeloid to br-met lesions. Inhibiting VISTA- and PD-L1-signaling relieved immune suppression and reduced br-met burden. Our results demonstrate that loss of Cx3cr1 in CNS-myeloid triggers a Cxcl10-mediated vicious cycle, cultivating a br-met-promoting, immune-suppressive niche.

Single nucleotide variants (SNVs) in intronic regions have yet to be systematically investigated for their disease-causing potential. Using known pathogenic and neutral intronic SNVs (iSNVs) as training data, we develop the RegSNPs-intron algorithm based on a random forest classifier that integrates RNA splicing, protein structure, and evolutionary conservation features. RegSNPs-intron showed excellent performance in evaluating the pathogenic impacts of iSNVs. Using a high-throughput functional reporter assay called ASSET-seq (ASsay for Splicing using ExonTrap and sequencing), we evaluate the impact of RegSNPs-intron predictions on splicing outcome. Together, RegSNPs-intron and ASSET-seq enable effective prioritization of iSNVs for disease pathogenesis.

While several multigene signatures are available for predicting breast cancer prognosis, particularly in early stage disease, effective molecular indicators are needed, especially for triple-negative carcinomas, to improve treatments and predict diagnostic outcomes. The objective of this study was to identify transcriptional regulatory networks to better understand mechanisms giving rise to breast cancer development and to incorporate this information into a model for predicting clinical outcomes.

Early human heart and brain development simultaneously occur during embryogenesis. Notably, in human newborns, congenital heart defects strongly associate with neurodevelopmental abnormalities, suggesting a common gene or complex underlying both cardiogenesis and neurogenesis. However, due to lack of in vivo studies, the molecular mechanisms that govern both early human heart and brain development remain elusive.

We report our discovery of an important player in the development of skin fibrosis, a hallmark of scleroderma. Scleroderma is a fibrotic disease, affecting 70,000 to 150,000 Americans. Fibrosis is a pathological wound healing process that produces an excessive extracellular matrix to interfere with normal organ function. Fibrosis contributes to nearly half of human mortality. Scleroderma has heterogeneous phenotypes, unpredictable outcomes, no validated biomarkers, and no effective treatment. Thus, strategies to slow down scleroderma progression represent an urgent medical need. While a pathological wound healing process like fibrosis leaves scars and weakens organ function, oral mucosa wound healing is a scarless process. After re-analyses of gene expression datasets from oral mucosa wound healing and skin fibrosis, we discovered that several pathways constitutively activated in skin fibrosis are transiently induced during oral mucosa wound healing process, particularly the amphiregulin (Areg) gene. Areg expression is upregulated ~ 10 folds 24hrs after oral mucosa wound but reduced to the basal level 3 days later. During bleomycin-induced skin fibrosis, a commonly used mouse model for skin fibrosis, Areg is up-regulated throughout the fibrogenesis and is associated with elevated cell proliferation in the dermis. To demonstrate the role of Areg for skin fibrosis, we used mice with Areg knockout, and found that Areg deficiency essentially prevents bleomycin-induced skin fibrosis. We further determined that bleomycin-induced cell proliferation in the dermis was not observed in the Areg null mice. Furthermore, we found that inhibiting MEK, a downstream signaling effector of Areg, by selumetinib also effectively blocked bleomycin-based skin fibrosis model. Based on these results, we concluded that the Areg-EGFR-MEK signaling axis is critical for skin fibrosis development. Blocking this signaling axis may be effective in treating scleroderma.

Single-cell RNA sequencing (scRNA-seq) is an evolving technology used to elucidate the cellular architecture of adult organs. Previous scRNA-seq on breast tissue utilized reduction mammoplasty samples, which are often histologically abnormal. We report a rapid tissue collection/processing protocol to perform scRNA-seq of breast biopsies of healthy women and identify 23 breast epithelial cell clusters. Putative cell-of-origin signatures derived from these clusters are applied to analyze transcriptomes of ~3,000 breast cancers. Gene signatures derived from mature luminal cell clusters are enriched in ~68% of breast cancers, whereas a signature from a luminal progenitor cluster is enriched in ~20% of breast cancers. Overexpression of luminal progenitor cluster-derived signatures in HER2+, but not in other subtypes, is associated with unfavorable outcome. We identify TBX3 and PDK4 as genes co-expressed with estrogen receptor (ER) in the normal breasts, and their expression analyses in >550 breast cancers enable prognostically relevant subclassification of ER+ breast cancers.

Kidney fibrosis is a hallmark of chronic kidney disease (CKD) progression that is caused by tubular injury and dysregulated lipid metabolism. Genetic abolition fatty acid-binding protein 4 (FABP4), a key lipid transporter, has been reported to suppress kidney interstitial fibrosis. However, the role and underlying mechanism of chemical inhibition of FABP4 in fibrotic kidney have not been well-documented. Here, we examined preemptive the effect of a FABP4 inhibitor, BMS309403, on lipid metabolism of tubular epithelial cells (TECs) and progression of kidney fibrosis. The expression of FABP4 was significantly elevated, concomitated with the accumulation of lipid droplets in TECs during kidney fibrosis. Treatment with BMS309403 alleviated lipid deposition of TECs, as well as interstitial fibrotic responses both in unilateral ureteral obstruction (UUO)-engaged mice and TGF-β-induced TECs. Moreover, BMS309403 administration enhanced fatty acid oxidation (FAO) in TECs by regulating peroxisome proliferator-activated receptor γ (PPARγ) and restoring FAO-related enzyme activities; In addition, BMS309403 markedly reduced cell lipotoxicity, such as endoplasmic reticulum (ER) stress and apoptosis in fibrotic kidney. Taken together, our results suggest that preemptive pharmacological inhibition of FABP4 by BMS309403 rebalances abnormal lipid metabolism in TECs and attenuates the progression of kidney fibrosis, thus may hold therapeutic potential for the treatment of fibrotic kidney diseases.

The Komen Tissue Bank is the only biorepository in the world for normal breast tissues from women. Below we report the acquisition and processing of breast tissue from volunteer donors and describe an experimental and analysis pipeline to generate a single-cell atlas. This atlas is based on single-cell RNA-seq and is useful to derive breast epithelial cell subcluster-specific gene expression signatures, which can be applied to breast cancer gene expression data to identify putative cell-of-origin. For complete details on the use and execution of this protocol, please refer to Bhat-Nakshatri et al. (2021).

Improved understanding of local breast biology that favors the development of estrogen receptor negative (ER-) breast cancer (BC) would foster better prevention strategies. We have previously shown that overexpression of specific lipid metabolism genes is associated with the development of ER- BC. We now report results of exposure of MCF-10A and MCF-12A cells, and mammary organoids to representative medium- and long-chain polyunsaturated fatty acids. This exposure caused a dynamic and profound change in gene expression, accompanied by changes in chromatin packing density, chromatin accessibility, and histone posttranslational modifications (PTMs). We identified 38 metabolic reactions that showed significantly increased activity, including reactions related to one-carbon metabolism. Among these reactions are those that produce S-adenosyl-L-methionine for histone PTMs. Utilizing both an in-vitro model and samples from women at high risk for ER- BC, we show that lipid exposure engenders gene expression, signaling pathway activation, and histone marks associated with the development of ER- BC.

Single-cell transcriptomics studies have begun to identify breast epithelial cell and stromal cell specific transcriptome differences between BRCA1/2 mutation carriers and non-carriers. We generated a single-cell transcriptome atlas of breast tissues from BRCA1, BRCA2 mutation carriers and compared this single-cell atlas of mutation carriers with our previously described single-cell breast atlas of healthy non-carriers. We observed that BRCA1 but not BRCA2 mutations altered the ratio between basal (basal-myoepithelial), luminal progenitor (luminal adaptive secretory precursor, LASP), and mature luminal (luminal hormone sensing) cells in breast tissues. A unique subcluster of cells within LASP cells is underrepresented in case of BRCA1 and BRCA2 mutation carriers compared with non-carriers. Both BRCA1 and BRCA2 mutations specifically altered transcriptomes in epithelial cells which are an integral part of NFκB, LARP1, and MYC signaling. Signaling pathway alterations in epithelial cells unique to BRCA1 mutations included STAT3, BRD4, SMARCA4, HIF2A/EPAS1, and Inhibin A signaling. BRCA2 mutations were associated with upregulation of IL6, PDK1, FOXO3, and TNFSF11 signaling. These signaling pathway alterations are sufficient to alter sensitivity of BRCA1/BRCA2-mutant breast epithelial cells to transformation as epithelial cells from BRCA1 mutation carriers overexpressing hTERT + PIK3CAH1047R generated adenocarcinomas, whereas similarly modified mutant BRCA2 cells generated basal carcinomas in NSG mice. Thus, our studies provide a high-resolution transcriptome atlas of breast epithelial cells of BRCA1 and BRCA2 mutation carriers and reveal their susceptibility to PIK3CA mutation-driven transformation.

Bempedoic acid, an innovative oral medication, has garnered significant interest in recent times due to its potential as a therapeutic intervention for hypercholesterolemia. Nonetheless, the outcomes of the initial investigations might have been more definitive and coherent. Our objective was to perform a quantitative meta-analysis in order to evaluate bempedoic acid's safety and effectiveness.

Cachexia affects the majority of cancer patients, with currently no effective treatments. Cachexia is defined by increased fatigue and loss of muscle function resulting from muscle and fat depletion. Previous studies suggest that chemotherapy may contribute to cachexia, although the causes responsible for this association are not clear. The purpose of this study was to investigate the mechanism(s) associated with chemotherapy-related effects on body composition and muscle function. Normal mice were administered chemotherapy regimens used for the treatment of colorectal cancer, such as Folfox (5-FU, leucovorin, oxaliplatin) or Folfiri (5-FU, leucovorin, irinotecan) for 5 weeks. The animals that received chemotherapy exhibited concurrent loss of muscle mass and muscle weakness. Consistently with previous findings, muscle wasting was associated with up-regulation of ERK1/2 and p38 MAPKs. No changes in ubiquitin-dependent proteolysis or in the expression of TGFβ-family members were detected. Further, marked decreases in mitochondrial content, associated with abnormalities at the sarcomeric level and with increase in the number of glycolytic fibers were observed in the muscle of mice receiving chemotherapy. Finally, ACVR2B/Fc or PD98059 prevented Folfiri-associated ERK1/2 activation and myofiber atrophy in C2C12 cultures. Our findings demonstrate that chemotherapy promotes MAPK-dependent muscle atrophy as well as mitochondrial depletion and alterations of the sarcomeric units. Therefore, these findings suggest that chemotherapy potentially plays a causative role in the occurrence of muscle loss and weakness. Moreover, the present observations provide a strong rationale for testing ACVR2B/Fc or MEK1 inhibitors in combination with anticancer drugs as novel strategies aimed at preventing chemotherapy-associated muscle atrophy.

Next-generation sequencing (NGS) studies have identified large numbers of genetic variants that are predicted to alter miRNA-mRNA interactions. We developed a novel high-throughput bioassay, PASSPORT-seq, that can functionally test in parallel 100s of these variants in miRNA binding sites (mirSNPs). The results are highly reproducible across both technical and biological replicates. The utility of the bioassay was demonstrated by testing 100 mirSNPs in HEK293, HepG2, and HeLa cells. The results of several of the variants were validated in all three cell lines using traditional individual luciferase assays. Fifty-five mirSNPs were functional in at least one of three cell lines (FDR ≤ 0.05); 11, 36, and 27 of them were functional in HEK293, HepG2, and HeLa cells, respectively. Only four of the variants were functional in all three cell lines, which demonstrates the cell-type specific effects of mirSNPs and the importance of testing the mirSNPs in multiple cell lines. Using PASSPORT-seq, we functionally tested 111 variants in the 3' UTR of 17 pharmacogenes that are predicted to alter miRNA regulation. Thirty-three of the variants tested were functional in at least one cell line.

Characterizing the biogeography of the microbiome of healthy humans is essential for understanding microbial associated diseases. Previous studies mainly focused on a single body habitat from a limited set of subjects. Here, we analyzed one of the largest microbiome datasets to date and generated a biogeographical map that annotates the biodiversity, spatial relationships, and temporal stability of 22 habitats from 279 healthy humans.

Although IL-9 has potent anti-tumor activity in adoptive cell transfer therapy, some models suggest that it can promote tumor growth. Here, we show that IL-9 signaling is associated with poor outcomes in patients with various forms of lung cancer, and is required for lung tumor growth in multiple mouse models. CD4+ T cell-derived IL-9 promotes the expansion of both CD11c+ and CD11c- interstitial macrophage populations in lung tumor models. Mechanistically, the IL-9/macrophage axis requires arginase 1 (Arg1) to mediate tumor growth. Indeed, adoptive transfer of Arg1+ but not Arg1- lung macrophages to Il9r-/- mice promotes tumor growth. Moreover, targeting IL-9 signaling using macrophage-specific nanoparticles restricts lung tumor growth in mice. Lastly, elevated expression of IL-9R and Arg1 in tumor lesions is associated with poor prognosis in lung cancer patients. Thus, our study suggests the IL-9/macrophage/Arg1 axis is a potential therapeutic target for lung cancer therapy.

Obesity is closely associated with adipose tissue inflammation and insulin resistance. Dysglycemia and type 2 diabetes results when islet β cells fail to maintain appropriate insulin secretion in the face of insulin resistance. To clarify the early transcriptional events leading to β-cell failure in the setting of obesity, we fed male C57BL/6J mice an obesogenic, high-fat diet (60% kcal from fat) or a control diet (10% kcal from fat) for one week, and islets from these mice (from four high-fat- and three control-fed mice) were subjected to single-cell RNA sequencing (sc-RNAseq) analysis. Islet endocrine cell types (α cells, β cells, δ cells, PP cells) and other resident cell types (macrophages, T cells) were annotated by transcript profiles and visualized using Uniform Manifold Approximation and Projection for Dimension Reduction (UMAP) plots. UMAP analysis revealed distinct cell clusters (11 for β cells, 5 for α cells, 3 for δ cells, PP cells, ductal cells, endothelial cells), emphasizing the heterogeneity of cell populations in the islet. Collectively, the clusters containing the majority of β cells showed the fewest gene expression changes, whereas clusters harboring the minority of β cells showed the most changes. We identified that distinct β-cell clusters downregulate genes associated with the endoplasmic reticulum stress response and upregulate genes associated with insulin secretion, whereas others upregulate genes that impair insulin secretion, cell proliferation, and cell survival. Moreover, all β-cell clusters negatively regulate genes associated with immune response activation. Glucagon-producing α cells exhibited patterns similar to β cells but, again, in clusters containing the minority of α cells. Our data indicate that an early transcriptional response in islets to an obesogenic diet reflects an attempt by distinct populations of β cells to augment or impair cellular function and/or reduce inflammatory responses as possible harbingers of ensuing insulin resistance.

Single-cell RNA sequencing reveals gene expression differences between individual cells and also identifies different cell populations that are present in the bulk starting material. To obtain an accurate assessment of patient samples, single-cell suspensions need to be generated as soon as possible once the tissue or sample has been collected. However, this requirement poses logistical challenges for experimental designs involving multiple samples from the same subject since these samples would ideally be processed at the same time to minimize technical variation in data analysis. Although cryopreservation has been shown to largely preserve the transcriptome, it is unclear whether the freeze-thaw process might alter gene expression profiles in a cell-type specific manner or whether changes in cell-type proportions might also occur. To address these questions in the context of multiple myeloma clinical studies, we performed single-cell RNA sequencing (scRNA-seq) to compare fresh and frozen cells isolated from bone marrow aspirates of six multiple myeloma patients, analyzing both myeloma cells (CD138+) and cells constituting the microenvironment (CD138-). We found that cryopreservation using 90% fetal calf serum and 10% dimethyl sulfoxide resulted in highly consistent gene expression profiles when comparing fresh and frozen samples from the same patient for both CD138+ myeloma cells (R ≥ 0.96) and for CD138- cells (R ≥ 0.9). We also demonstrate that CD138- cell-type proportions showed minimal alterations, which were mainly related to small differences in immune cell subtype sensitivity to the freeze-thaw procedures. Therefore, when processing fresh multiple myeloma samples is not feasible, cryopreservation is a useful option in single-cell profiling studies.