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.

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.

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.

Mycophenolic acid (MPA) is an immunosuppressant commonly used to prevent renal transplant rejection and treat glomerulonephritis. MPA inhibits IMPDH2 within stimulated lymphocytes, reducing guanosine synthesis. Despite the widespread use of MPA, interindividual variability in response remains with rates of allograft rejection up to 15% and approximately half of individuals fail to achieve complete remission to lupus nephritis. We sought to identify contributors to interindividual variability in MPA response, hypothesizing that the HPRT1 salvage guanosine synthesis contributes to variability. MPA sensitivity was measured in 40 healthy individuals using an ex vivo lymphocyte viability assay. Measurement of candidate gene expression (n ± 40) and single-cell RNA-sequencing (n ± 6) in lymphocytes was performed at baseline, poststimulation, and post-MPA treatment. After stimulation, HPRT1 expression was 2.1-fold higher in resistant individuals compared with sensitive individuals (P ± 0.049). Knockdown of HPRT1 increased MPA sensitivity (12%; P ± 0.003), consistent with higher expression levels in resistant individuals. Sensitive individuals had higher IMPDH2 expression and 132% greater stimulation. In lymphocyte subpopulations, differentially expressed genes between sensitive and resistant individuals included KLF2 and LTB. Knockdown of KLF2 and LTB aligned with the predicted direction of effect on proliferation. In sensitive individuals, more frequent receptor-ligand interactions were observed after stimulation (P ± 0.0004), but fewer interactions remained after MPA treatment (P ± 0.0014). These data identify a polygenic transcriptomic signature in lymphocyte subpopulations predictive of MPA response. The degree of lymphocyte stimulation, HPRT1, KLF2, and LTB expression may serve as markers of MPA efficacy.

Myasthenia gravis (MG) is an autoimmune disorder resulting from antibodies against the proteins at the neuromuscular junction. Emerging evidence indicates that long non-coding RNAs (lncRNAs), acting as competing endogenous RNAs (ceRNAs), are involved in various diseases. However, the regulatory mechanisms of ceRNAs underlying MG remain largely unknown. In this study, we constructed a lncRNA-mediated ceRNA network involved in MG using a multi-step computational strategy. Functional annotation analysis suggests that these lncRNAs may play crucial roles in the immunological mechanism underlying MG. Importantly, through manual literature mining, we found that lncRNA SNHG16 (small nucleolar RNA host gene 16), acting as a ceRNA, plays important roles in the immune processes. Further experiments showed that SNHG16 expression was upregulated in peripheral blood mononuclear cells (PBMCs) from MG patients compared to healthy controls. Luciferase reporter assays confirmed that SNHG16 is a target of the microRNA (miRNA) let-7c-5p. Subsequent experiments indicated that SNHG16 regulates the expression of the key MG gene interleukin (IL)-10 by sponging let-7c-5p in a ceRNA manner. Furthermore, functional assays showed that SNHG16 inhibits Jurkat cell apoptosis and promotes cell proliferation by sponging let-7c-5p. Our study will contribute to a deeper understanding of the regulatory mechanism of MG and will potentially provide new therapeutic targets for MG patients.

Phytophthora sojae is an oomycete that causes stem and root rot disease in soybean. P. sojae delivers many RxLR effector proteins, including Avr1b, into host cells to promote infection. We show here that Avr1b interacts with the soybean U-box protein, GmPUB1-1, in yeast two-hybrid, pull down, and bimolecular fluorescence complementation (BIFC) assays. GmPUB1-1, and a homeologous copy GmPUB1-2, are induced by infection and encode 403 amino acid proteins with U-Box domains at their N-termini. Non-synonymous mutations in the Avr1b C-terminus that abolish suppression of cell death also abolished the interaction of Avr1b with GmPUB1-1, while deletion of the GmPUB1-1 C-terminus, but not the U box, abolished the interaction. BIFC experiments suggested that the GmPUB1-1-Avr1b complex is targeted to the nucleus. In vitro ubiquitination assays demonstrated that GmPUB1-1 possesses E3 ligase activity. Silencing of the GmPUB1 genes in soybean cotyledons resulted in loss of recognition of Avr1b by gene products encoded by Rps1-b and Rps1-k. The recognition of Avr1k (which did not interact with GmPUB1-1) by Rps1-k plants was not, however, affected following GmPUB1-1 silencing. Furthermore, over-expression of GmPUB1-1 in particle bombardment experiments triggered cell death suggesting that GmPUB1 may be a positive regulator of effector-triggered immunity. In a yeast two-hybrid system, GmPUB1-1 also interacted with a number of other RxLR effectors including Avr1d, while Avr1b and Avr1d interacted with a number of other infection-induced GmPUB proteins, suggesting that the pathogen uses a multiplex of interactions of RxLR effectors with GmPUB proteins to modulate host immunity.

Myasthenia gravis (MG) is a complex neurological autoimmune disease with a pathogenetic mechanism that has yet to be elucidated. Emerging evidence has revealed that genes, non-coding RNAs and genetic variants play significant roles in the pathogenesis of MG. However, the molecular mechanisms of single nucleotide polymorphisms (SNPs) located on lncRNAs could disturb lncRNA-mediated ceRNA regulatory functions still unclear in MG. In this study, we collated 276 experimentally confirmed MG risk genes and 192 MG risk miRNAs. We then constructed a lncRNA-mediated ceRNA network for MG based on multi-step computational strategies. Next, we systematically integrated risk pathways and identified candidate SNPs in lncRNAs for MG based on data acquired from public databases. In addition, we constructed a pathway-based lncRNA-SNP mediated network (LSPN) that contained 128 lncRNAs targeting 8 MG risk pathways. By analyzing network, we propose a latent mechanism for how the "lncRNA-SNP-mRNA-pathway" axis affects the pathogenesis of MG. Moreover, 25 lncRNAs and 51 SNPs on lncRNAs were extracted from the "lncRNA-SNP-mRNA-pathway" axis. Finally, functional analyses demonstrated lncRNA-SNPs mediated ceRNA regulation pairs associated with MG participated in the MAPK signaling pathway. In summary, we constructed MG-specific lncRNA-SNPs mediated ceRNA regulatory networks based on pathway in the present study, which was helpful to elucidate the roles of lncRNA-SNPs in the pathogenesis of MG and provide novel insights into mechanism of lncRNA-SNPs as potential genetic risk biomarkers of MG.

The efficacy and longevity of medical implants and devices is largely determined by the host immune response, which extends along a continuum from pro-inflammatory/pro-fibrotic to anti-inflammatory/pro-regenerative. Using a rat subcutaneous implantation model, along with histological and transcriptomics analyses, we characterized the tissue response to a collagen polymeric scaffold fabricated from polymerizable type I oligomeric collagen (Oligomer) in comparison to commercial synthetic and collagen-based products. In contrast to commercial biomaterials, no evidence of an immune-mediated foreign body reaction, fibrosis, or bioresorption was observed with Oligomer scaffolds for beyond 60 days. Oligomer scaffolds were noninflammatory, eliciting minimal innate inflammation and immune cell accumulation similar to sham surgical controls. Genes associated with Th2 and regulatory T cells were instead upregulated, implying a novel pathway to immune tolerance and regenerative remodeling for biomaterials.

Populations that are born and raised at high altitude develop under conditions of chronic developmental hypoxia (CDH), which results in pulmonary adaptations of increased lung volume and diffusion capacity to increase gas exchange. It is not clear how CDH may alter allergic inflammation in the lung. In this study, we sought to characterize the impact of CDH on immune cell populations in the rat lung during a murine model of asthma. Rats were bred and raised in either hypoxic (15% oxygen, CDH) or normobaric room air (20% oxygen). At 3-weeks of age, animals were sensitized to ovalbumin (OVA) or physiologic saline (phosphate-buffered saline [PBS]) as a control, followed by three consecutive days of intra-nasal OVA or PBS at 6-weeks of age. We then assessed airway reactivity and allergic-associated cytokine levels. This was followed by single-cell transcriptomic profiling of lung cell populations. In scRNA-seq analysis, we assessed differentially expressed genes, differentially enriched functional pathways, immune cell exhaustion/activation markers, and immune cell secretory products. Our results show that while OVA heightened airway reactivity, CDH suppressed airway reactivity in OVA-challenged and control animals. Through scRNA-seq analysis, we further demonstrate that CDH alters the transcriptional landscape in the lung and alters transcriptional programs in immune cells. These data define CDH-dependent changes in the lung that impact airway reactivity.