Single-cell analyses reveal evolution mimicry during the specification of breast cancer subtype.

Background: The stem or progenitor antecedents confer developmental plasticity and unique cell identities to cancer cells via genetic and epigenetic programs. A comprehensive characterization and mapping of the cell-of-origin of breast cancer using novel technologies to unveil novel subtype-specific therapeutic targets is still absent. Methods: We integrated 195,144 high-quality cells from normal breast tissues and 406,501 high-quality cells from primary breast cancer samples to create a large-scale single-cell atlas of human normal and cancerous breasts. Potential heterogeneous origin of malignant cells was explored by contrasting cancer cells against reference normal epithelial cells. Multi-omics analyses and both in vitro and in vivo experiments were performed to screen and validate potential subtype-specific treatment targets. Novel biomarkers of identified immune and stromal cell subpopulations were validated by immunohistochemistry in our cohort. Results: Tumor stratification based on cancer cell-of-origin patterns correlated with clinical outcomes, genomic aberrations and diverse microenvironment constitutions. We found that the luminal progenitor (LP) subtype was robustly associated with poor prognosis, genomic instability and dysfunctional immune microenvironment. However, the LP subtype patients were sensitive to neoadjuvant chemotherapy (NAC), PARP inhibitors (PARPi) and immunotherapy. The LP subtype-specific target PLK1 was investigated by both in vitro and in vivo experiments. Besides, large-scale single-cell profiling of breast cancer inspired us to identify a range of clinically relevant immune and stromal cell subpopulations, including subsets of innate lymphoid cells (ILCs), macrophages and endothelial cells. Conclusion: The present single-cell study revealed the cellular repertoire and cell-of-origin patterns of breast cancer. Combining single-cell and bulk transcriptome data, we elucidated the evolution mimicry from normal to malignant subtypes and expounded the LP subtype with vital clinical implications. Novel immune and stromal cell subpopulations of breast cancer identified in our study could be potential therapeutic targets. Taken together, Our findings lay the foundation for the precise prognostic and therapeutic stratification of breast cancer.

Single-cell analysis of white adipose tissue reveals the tumor-promoting adipocyte subtypes.

BACKGROUND: The tumor-adipose microenvironment (TAME) is characterized by the enrichment of adipocytes, and is considered a special ecosystem that supports cancer progression. However, the heterogeneity and diversity of adipocytes in TAME remains poorly understood. METHODS: We conducted a single-cell RNA sequencing analysis of adipocytes in mouse and human white adipose tissue (WAT). We analyzed several adipocyte subtypes to evaluate their relationship and potential as prognostic factors for overall survival (OS). The potential drugs are screened by using bioinformatics methods. The tumor-promoting effects of a typical adipocyte subtype in breast cancer are validated by performing in vitro functional assays and immunohistochemistry (IHC) in clinical samples. RESULTS: We profiled a comprehensive single-cell atlas of adipocyte in mouse and human WAT and described their characteristics, origins, development, functions and interactions with immune cells. Several cancer-associated adipocyte subtypes, namely DPP4+ adipocytes in visceral adipose and ADIPOQ+ adipocytes in subcutaneous adipose, are identified. We found that high levels of these subtypes are associated with unfavorable outcomes in four typical adipose-associated cancers. Some potential drugs including Trametinib, Selumetinib and Ulixertinib are discovered. Emphatically, knockdown of adiponectin receptor 1 (AdipoR1) and AdipoR2 impaired the proliferation and invasion of breast cancer cells. Patients with AdipoR2-high breast cancer display significantly shorter relapse-free survival (RFS) than those with AdipoR2-low breast cancer. CONCLUSION: Our results provide a novel understanding of TAME at the single-cell level. Based on our findings, several adipocyte subtypes have negative impact on prognosis. These cancer-associated adipocytes may serve as key prognostic predictor and potential targets for treatment in the future.

The ALOX5 inhibitor Zileuton regulates tumor-associated macrophage M2 polarization by JAK/STAT and inhibits pancreatic cancer invasion and metastasis.

5-lipoxygenase (encoded by ALOX5) plays an important role in immune regulation. Zileuton is currently the only approved ALOX5 inhibitor. However, the mechanisms of ALOX5 and Zileuton in progression of pancreatic cancer remain unclear. Therefore, we investigated the effects of Zileuton on tumor-associated macrophage M2 polarization and pancreatic cancer invasion and metastasis, both in vivo and in vitro. In bulk RNA sequencing (RNA-seq) and single-cell RNA sequencing (scRNA-seq) analyses, we found a significant association between elevated levels of ALOX5 and poor survival, adverse stages, M2 macrophage infiltration, and the activation of JAK/STAT pathways in macrophages. In clinical samples, immunofluorescence, quantitative real-time PCR and immunohistochemical results verified the high expression of ALOX5 in pancreatic cancer, primarily in macrophages. We constructed PANC-1 human pancreatic cancer cells and macrophages overexpressing ALOX5 using lentivirus. In PANC-1 pancreatic cancer cells, low-dose Zileuton inhibited PANC-1 cell invasion and migration by blocking ALOX5. In macrophages, ALOX5 induced the M2-like phenotype through the JAK/STAT pathway and promoted the chemotaxis of macrophages towards PANC-1 cells, while Zileuton can inhibit these effects. We constructed the nude mouse model of in situ transplantation tumor of pancreatic cancer. After treatment with Zileuton, the mice showed increased survival rates and reduced liver metastasis. These findings indicate that ALOX5 regulates tumor-associated macrophage M2 polarization via the JAK/STAT pathway and promotes invasion and metastasis in pancreatic cancer. Zileuton can inhibit these effects by inhibiting ALOX5. These results provide a theoretical basis for the potential use of Zileuton in the treatment of pancreatic cancer.

Neuroendocrine regulations in tissue-specific immunity: From mechanism to applications in tumor.

Immune responses in nonlymphoid tissues play a vital role in the maintenance of homeostasis. Lots of evidence supports that tissue-specific immune cells provide defense against tumor through the localization in different tissue throughout the body, and can be regulated by diverse factors. Accordingly, the distribution of nervous tissue is also tissue-specific which is essential in the growth of corresponding organs, and the occurrence and development of tumor. Although there have been many mature perspectives on the neuroendocrine regulation in tumor microenvironment, the neuroendocrine regulation of tissue-specific immune cells has not yet been summarized. In this review, we focus on how tissue immune responses are influenced by autonomic nervous system, sensory nerves, and various neuroendocrine factors and reversely how tissue-specific immune cells communicate with neuroendocrine system through releasing different factors. Furthermore, we pay attention to the potential mechanisms of neuroendocrine-tissue specific immunity axis involved in tumors. This may provide new insights for the immunotherapy of tumors in the future.

Single-cell and spatially resolved analysis uncovers cell heterogeneity of breast cancer.

The heterogeneity and the complex cellular architecture have a crucial effect on breast cancer progression and response to treatment. However, deciphering the neoplastic subtypes and their spatial organization is still challenging. Here, we combine single-nucleus RNA sequencing (snRNA-seq) with a microarray-based spatial transcriptomics (ST) to identify cell populations and their spatial distribution in breast cancer tissues. Malignant cells are clustered into distinct subpopulations. These cell clusters not only have diverse features, origins and functions, but also emerge to the crosstalk within subtypes. Furthermore, we find that these subclusters are mapped in distinct tissue regions, where discrepant enrichment of stromal cell types are observed. We also inferred the abundance of these tumorous subpopulations by deconvolution of large breast cancer RNA-seq cohorts, revealing differential association with patient survival and therapeutic response. Our study provides a novel insight for the cellular architecture of breast cancer and potential therapeutic strategies.

Sequence duplication in 3' UTR modulates virus replication and virulence of Japanese encephalitis virus.

Japanese encephalitis virus (JEV), an important neurotropic pathogen, belongs to the genus Flavivirus of the family Flaviviridae and has caused huge threat to public health. It is still obscure regarding the functions of stem-loop (SL) and dumbbell (DB) domains of JEV 3' UTR in viral replication and virulence. In the current study, using the infectious clone of JEV SA14 strain as a backbone, we constructed a series of deletion mutants of 3' UTR to investigate their effects on virus replication. The results showed that partial deletions within SL or DB domain had no apparent effects on virus replication in both mammalian (BHK-21) and mosquito (C6/36) cells, suggesting that they were not involved in viral host-specific replication. However, the entire SL domain deletion (ΔVR) significantly reduced virus replication in both cell lines, indicating the important role of the complete SL domain in virus replication. The revertant of ΔVR mutant virus was obtained by serial passage in BHK-21 cells that acquired a duplication of DB domain (DB-dup) in the 3' UTR, which greatly restored virus replication as well as the capability to produce the subgenomic flavivirus RNAs (sfRNAs). Interestingly, the DB-dup mutant virus was highly attenuated in C57BL/6 mice despite replicating similar to WT JEV. These findings demonstrate the significant roles of the duplicated structures in 3' UTR in JEV replication and provide a novel strategy for the design of live-attenuated vaccines.

Effects of a contusion load on spinal cord with different curvatures.

The cervical spine injury is a complicated procedure in the combination of different injury loads and postures. The aim of this study is to investigate the injury mechanism considering different types of cervical curvatures subjected to contusion loads. A finite element model of a cervical spinal cord was constructed. Gray matter, white matter and pia matter were modeled and hyperelastic material properties were assigned. Convergence analysis and validation analysis were carried out. The model was simulated in 3 different spinal curvatures and loaded by 2 directions with 4 compression degrees. The maximum von Mises stress in the whole model was concentrated in the pia matter in all loading cases. When investigating spinal cord injury, the pia matter must be considered. For all three curvatures, the stress in the gray matter and white matter was higher in front-to-back loading condition than that in back-to-front loading condition. The front-to-back impact may cause a larger damage. A back-to-front load damaged the structure around the central canal and a front-to-back contusion load damaged the anterior horn of the spinal cord at most time. From the view of the maximum stress, the lordotic curvature did not show significant buffering effect. However, the pathological curvature had large areas affected and the lordotic curvature showed some benefits to some degree from the view of stress distribution.

Effects of contusion load on cervical spinal cord: A finite element study.

Injury of cervical spine is a common injury of locomotor system usually accompanied by spinal cord injury, however the injury mechanism of contusion load to the spinal cord is not clear. This study aims to investigate its injury mechanism associated with the contusion load, with different extents of spinal cord compression. A finite element model of cervical spinal cord was established and two scenarios of contusion injury loading conditions, i.e. back-to-front and front-to-back loads, were adopted. Four different compression displacements were applied to the middle section of the cervical spinal cord. The distributions of von Mises stress in middle transverse cross section were obtained from the finite element analysis. For the back-to-front loading scenario, the stress concentration was found in the area at and near the central canal and the damage may lead to the central canal syndrome from biomechanical point of view. With the front-to-back load, the maximum von Mises stress located in central canal area of gray matter when subject to 10% compression, whilst it appeared at the anterior horn when the compression increased. For the white matter, the maximum von Mises stress appeared in the area of the anterior funiculus. This leads to complicated symptoms given rise by damage to multiple locations in the cervical spinal cord. The illustrative results demonstrated the need of considering different loading scenarios in understanding the damage mechanisms of the cervical spinal cord, particularly when the loading conditions were given rise by different pathophysiological causes.

Development and evaluation of one-step multiplex real-time RT-PCR assay for simultaneous detection of Zika virus and Chikungunya virus.

Zika virus (ZIKV) and chikungunya virus (CHIKV) are important human pathogens and mosquito-borne arboviruses, which have resembling history, common vectors, circulating regions, and indistinguishable clinical symptoms. Wide geographical range that is suitable for ZIKV and CHIKV transmission underlines the concern about the impact of epidemic and endemic infections on burden of public health. In the present study, a highly sensitive and specific one-step multiplex real-time RT-PCR assay was developed and evaluated for simultaneous detection and quantification of ZIKV and CHIKV. The single reaction assay employs two pairs of primers and two TaqMan probes that differentiate ZIKV and CHIKV infections. The entire viral genomic RNA in vitro transcribed from full-length infectious clones were used to generate the standard curves for absolute quantification in subsequent tests. The detection limit of the one-step multiplex assay was 1 and 0.5 PFU for infectious ZIKV and CHIKV, respectively. The assessment of specificity indicated this assay is highly specific to targeted viruses showing no amplification of a variety of other flaviviruses. Our assay was able to detect geographically separated and phylogenetically diverse strains of ZIKV and CHIKV. On the applicability of monitoring viral multiplication in cells and testing clinical samples, the one-step multiplex assay provided efficient and accurate determination. The one-step multiplex real-time RT-PCR assay offers a valuable tool for detection of ZIKV and CHIKV and potentially contributes to general surveillance and clinical treatment.

Detection, isolation, and characterization of chikungunya viruses associated with the Pakistan outbreak of 2016-2017.

The chikungunya virus (CHIKV) is a mosquito-transmitted alphavirus, which has infected millions of people in Africa, Asia, Americas, and Europe since it reemerged in India and Indian Ocean regions in 2005-2006. Starting in the middle of November 2016, CHIKV has been widely spread, and more than 4,000 cases of infections in humans were confirmed in Pakistan. Here, we report the first isolation and characterization of CHIKV from the Pakistan outbreak. Eight CHIKV strains were newly isolated from human serum samples using a cell culture procedure. A full-length genome sequence and eight complete envelope (E1) sequences of CHIKV from Pakistan were obtained in this study. Alignment of the CHIKV E1 sequences revealed that the eight new CHIKV isolates were highly homogeneous, with only two nonsynonymous substitutions found at generally conserved sites (E99 and Q235). Based on the comparison of 342 E1 sequences, the two nonsynonymous mutations were located in well-recognized domains associated with viral functions such as the cell fusion and vector specificity, suggesting their potential functional importance. Phylogenetic analysis indicated that the CHIKV strains from Pakistan originated from CHIKV circulating in the Indian region. This study helps elucidate the epidemics of CHIKV in Pakistan and also provides a foundation for studies of evolution and expansion of CHIKV in South Asia.

Development of a stable Japanese encephalitis virus replicon cell line for antiviral screening.

Japanese encephalitis virus (JEV), an important pathogen in Eastern and Southern Asia and the Pacific, has spread to Australia and other territories in recent years. Although the vaccine for JEV has been used in some countries, development of efficient antiviral drugs is still an urgent requirement. Replicon systems have been widely used in the research of viral replication and antiviral screening for West Nile virus (WNV), yellow fever virus (YFV) and dengue virus (DENV). Here, a novel JEV replicon harboring the Rluc and Pac gene (JEV-Pac-Rluc-Rep) was constructed. Furthermore, we established a BHK-21 cell line harboring JEV-Pac-Rluc-Rep (BHK-21/PAC/Rluc cell line) through continuous puromycin selection. Characterization of cell line stability showed that the replicon RNA could persistently replicate in this cell line for at least up to 10 rounds of passage. Using a known flavivirus inhibitor, the JEV replicon cell line was validated for antiviral screening. The JEV replicon cell line will be a valuable tool for both compound screening and viral replication studies.

Recovery of the Zika virus through an in vitro ligation approach.

In this study, an in vitro ligation method was developed to assemble a full-length infectious cDNA clone of the Zika virus (ZIKV). Four contiguous cDNA subclones covering the complete ZIKV genome were constructed with unique BglI restriction sites at the ends of each fragment. The BglI restriction sites only allow in vitro ligation to happen between interconnecting restriction sites from adjacent cDNA fragments, resulting in an intact full-length cDNA of ZIKV. RNA transcripts derived from the full-length cDNA were infectious. The recombinant virus replicated as efficiently as the wild-type virus with similar growth kinetics and plaque morphologies in Vero and C6/36 cells. Both viruses were inhibited by NITD008 treatment. This in vitro ligation method will facilitate manipulation of the viral genome through genetic modifications of four separated subclones of ZIKV for the rapid and rational development of candidate vaccines and viral replication study.

West Nile Virus NS1 Antagonizes Interferon Beta Production by Targeting RIG-I and MDA5.

West Nile virus (WNV) is a mosquito-borne flavivirus that causes epidemics of encephalitis and viscerotropic disease worldwide. This virus has spread rapidly and has posed a significant public health threat since the outbreak in New York City in 1999. The interferon (IFN)-mediated antiviral response represents an important component of virus-host interactions and plays an essential role in regulating viral replication. Previous studies have suggested that multifunctional nonstructural proteins encoded by flaviviruses antagonize the host IFN response via various means in order to establish efficient viral replication. In this study, we demonstrated that the nonstructural protein 1 (NS1) of WNV antagonizes IFN-ß production, most likely through suppression of retinoic acid-inducible gene I (RIG-I)-like receptor (RLR) activation. In a dual-luciferase reporter assay, WNV NS1 significantly inhibited the activation of the IFN-ß promoter after Sendai virus infection or poly(I·C) treatment. NS1 also suppressed the activation of the IFN-ß promoter when it was stimulated by interferon regulatory factor 3 (IRF3)/5D or its upstream molecules in the RLR signaling pathway. Furthermore, NS1 blocked the phosphorylation and nuclear translocation of IRF3 upon stimulation by various inducers. Mechanistically, WNV NS1 targets RIG-I and melanoma differentiation-associated gene 5 (MDA5) by interacting with them and subsequently causing their degradation by the proteasome. Furthermore, WNV NS1 inhibits the K63-linked polyubiquitination of RIG-I, thereby inhibiting the activation of downstream sensors in the RLR signaling pathway. Taken together, our results reveal a novel mechanism by which WNV NS1 interferes with the host antiviral response.IMPORTANCE WNV Nile virus (WNV) has received increased attention since its introduction to the United States. However, the pathogenesis of this virus is poorly understood. This study demonstrated that the nonstructural protein 1 (NS1) of WNV antagonizes the induction of interferon beta (IFN-ß) by interacting with and degrading retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDA5), which are crucial viral sensors in the host innate immune system. Further experiments suggested that NS1-mediated inhibition of the RIG-I-like receptor (RLR) signaling pathway involves inhibition of RIG-I K63-linked polyubiquitination and that the proteasome plays a role in RIG-I degradation. This study provides new insights into the regulation of WNV NS1 in the RLR signaling pathway and reveals a novel mechanism by which WNV evades the host innate immune response. The novel findings may guide us to discover new therapeutic targets and develop effective vaccines for WNV infections.

Generation and characterization of West Nile pseudo-infectious reporter virus for antiviral screening.

West Nile virus (WNV), a mosquito-borne flavivirus, is an important neurotropic human pathogen. As a biosafety level-3 (BSL-3) agent, WNV is strictly to BSL-3 laboratories for experimentations, thus greatly hindering the development of vaccine and antiviral drug. Here, we developed a novel pseudo-infectious WNV reporter virus expressing the Gaussia luciferase (Gluc). A stable 293TNS1 cell line expressing NS1 was selected for trans-supplying NS1 protein to support the replication of WNV-ΔNS1 virus and WNV-ΔNS1-Gluc reporter virus with large-fragment deletion of NS1. WNV-ΔNS1 virus and WNV-Gluc-ΔNS1 reporter virus were confined to complete their replication cycle in this 293TNS1 cell line, displaying nearly identical growth kinetics to WT WNV although the viral titers were lower than those of WT WNV. The reporter gene was stably maintained in virus genome at least within three rounds of passage in 293TNS1 cell line. Using a known flaviviruses inhibitor, NITD008, we demonstrated that the pseudo-infectious WNV-Gluc-ΔNS1 could be used for antiviral screening. Furthermore, a high-throughput screening (HTS) assay in a 96-well format was optimized and validated using several known WNV inhibitors, indicating that the optimized HTS assay was suitable for high-throughput screening WNV inhibitors. Our work provides a stable and safe tool to handle WNV outside of a BSL-3 facility and facilitates high throughput screening for anti-WNV drugs.

Detection of Zika virus by SYBR green one-step real-time RT-PCR.

The ongoing Zika virus (ZIKV) outbreak has rapidly spread to new areas of Americas, which were the first transmissions outside its traditional endemic areas in Africa and Asia. Due to the link with newborn defects and neurological disorder, numerous infected cases throughout the world and various mosquito vectors, the virus has been considered to be an international public health emergency. In the present study, we developed a SYBR Green based one-step real-time RT-PCR assay for rapid detection of ZIKV. Our results revealed that the real-time assay is highly specific and sensitive in detection of ZIKV in cell samples. Importantly, the replication of ZIKV at different time points in infected cells could be rapidly monitored by the real-time RT-PCR assay. Specifically, the real-time RT-PCR showed acceptable performance in measurement of infectious ZIKV RNA. This assay could detect ZIKV at a titer as low as 1PFU/mL. The real-time RT-PCR assay could be a useful tool for further virology surveillance and diagnosis of ZIKV.

Development of Neutralization Assay Using an eGFP Chikungunya Virus.

Chikungunya virus (CHIKV), a member of the Alphavirus genus, is an important human emerging/re-emerging pathogen. Currently, there are no effective antiviral drugs or vaccines against CHIKV infection. Herein, we construct an infectious clone of CHIKV and an eGFP reporter CHIKV (eGFP-CHIKV) with an isolated strain (assigned to Asian lineage) from CHIKV-infected patients. The eGFP-CHIKV reporter virus allows for direct visualization of viral replication through the levels of eGFP expression. Using a known CHIKV inhibitor, ribavirin, we confirmed that the eGFP-CHIKV reporter virus could be used to identify inhibitors against CHIKV. Importantly, we developed a novel and reliable eGFP-CHIKV reporter virus-based neutralization assay that could be used for rapid screening neutralizing antibodies against CHIKV.

Generation of a recombinant West Nile virus stably expressing the Gaussia luciferase for neutralization assay.

West Nile virus (WNV) is a neurotropic human pathogen that has caused increasing infected cases over recent years. There is currently no licensed vaccine or effective drug for prevention and treatment of WNV infection in humans. To facilitate antiviral drug discovery and neutralizing antibody detection, a WNV cDNA clone containing a luciferase reporter gene was constructed through incorporating Gaussia luciferase (Gluc) gene within the capsid-coding region of WNV genome. Transfection of BHK-21 cells with the cDNA clone-derived RNA generated luciferase reporter WNV (WNV-Gluc) and the stable WNV-Gluc with high titers (>10(7)PFU/ml) was obtained through plaque purification. Luciferase activity was used to effectively quantify the viral production of WNV-Gluc. Using the reporter virus WNV-Gluc, we developed a luciferase based assay in a 12-well format for evaluating neutralizing antibodies. The reporter virus could be a powerful tool for epidemiological investigation of WNV, vaccine evaluation, antiviral drug screening, and the study of WNV replication and pathogenesis.

Investigation of the genotype III to genotype I shift in Japanese encephalitis virus and the impact on human cases.

Japanese encephalitis is a mosquito borne disease and is the leading cause of viral encephalitis in the Asia-Pacific area. The causative agent, Japanese encephalitis virus (JEV) can be phylogenetically classified into five genotypes based on nucleotide sequence. In recent years, genotype I (GI) has displaced genotype III (GIII) as the dominant lineage, but the mechanisms behind this displacement event requires elucidation. In an earlier study, we compared host variation over time between the two genotypes and observed that GI appears to have evolved to achieve more efficient infection in hosts in the replication cycle, with the tradeoff of reduced infectivity in secondary hosts such as humans. To further investigate this phenomenon, we collected JEV surveillance data on human cases and, together with sequence data, and generated genotype/case profiles from seven Asia-Pacific countries and regions to characterize the GI/GIII displacement event. We found that, when comprehensive and consistent vaccination and surveillance data was available, and the GIII to GI shift occurred within a well-defined time period, there was a statistically significant drop in JEV human cases. Our findings provide further support for the argument that GI is less effective in infecting humans, who represent a dead end host. However, experimental investigation is necessary to confirm this hypothesis. The study highlights the value of alternative approaches to investigation of epidemics, as well as the importance of effective data collection for disease surveillance and control.