A variational algorithm to detect the clonal copy number substructure of tumors from scRNA-seq data.

Single-cell RNA sequencing is the reference technology to characterize the composition of the tumor microenvironment and to study tumor heterogeneity at high resolution. Here we report Single CEll Variational ANeuploidy analysis (SCEVAN), a fast variational algorithm for the deconvolution of the clonal substructure of tumors from single-cell RNA-seq data. It uses a multichannel segmentation algorithm exploiting the assumption that all the cells in a given copy number clone share the same breakpoints. Thus, the smoothed expression profile of every individual cell constitutes part of the evidence of the copy number profile in each subclone. SCEVAN can automatically and accurately discriminate between malignant and non-malignant cells, resulting in a practical framework to analyze tumors and their microenvironment. We apply SCEVAN to datasets encompassing 106 samples and 93,322 cells from different tumor types and technologies. We demonstrate its application to characterize the intratumor heterogeneity and geographic evolution of malignant brain tumors.

Arabidopsis transcription factor TCP4 represses chlorophyll biosynthesis to prevent petal greening.

Green petals pose a challenge for pollinators to distinguish flowers from leaves, but they are valuable as a specialty flower trait. However, little is understood about the molecular mechanisms that underlie the development of green petals. Here, we report that CINCINNATA (CIN)-like TEOSINTE BRANCHED 1/CYCLOIDEA/PCF (TCP) proteins play key roles in the control of petal color. The septuple tcp2/3/4/5/10/13/17 mutant produced flowers with green petals due to chlorophyll accumulation. Expression of TCP4 complemented the petal phenotype of tcp2/3/4/5/10/13/17. We found that chloroplasts were converted into leucoplasts in the distal parts of wild-type petals but not in the proximal parts during flower development, whereas plastid conversion was compromised in the distal parts of tcp2/3/4/5/10/13/17 petals. TCP4 and most CIN-like TCPs were predominantly expressed in distal petal regions, consistent with the green-white pattern in wild-type petals and the petal greening observed in the distal parts of tcp2/3/4/5/10/13/17 petals. RNA-sequencing data revealed that most chlorophyll biosynthesis genes were downregulated in the white distal parts of wild-type petals, but these genes had elevated expression in the distal green parts of tcp2/3/4/5/10/13/17 petals and the green proximal parts of wild-type petals. We revealed that TCP4 repressed chlorophyll biosynthesis by directly binding to the promoters of PROTOCHLOROPHYLLIDE REDUCTASE (PORB), DIVINYL REDUCTASE (DVR), and SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), which are known to promote petal greening. We found that the conversion of chloroplasts to leucoplasts and the green coloration in the proximal parts of petals appeared to be conserved among plant species. Our findings uncover a major molecular mechanism that underpins the formation of petal color patterns and provide a foundation for the breeding of plants with green flowers.

Revealing atomic-scale molecular diffusion of a plant-transcription factor WRKY domain protein along DNA.

Transcription factor (TF) target search on genome is highly essential for gene expression and regulation. High-resolution determination of TF diffusion along DNA remains technically challenging. Here, we constructed a TF model system using the plant WRKY domain protein in complex with DNA from crystallography and demonstrated microsecond diffusion dynamics of WRKY on DNA by employing all-atom molecular-dynamics (MD) simulations. Notably, we found that WRKY preferentially binds to one strand of DNA with significant energetic bias compared with the other, or nonpreferred strand. The preferential DNA-strand binding becomes most prominent in the static process, from nonspecific to specific DNA binding, but less distinct during diffusive movements of the domain protein on the DNA. Remarkably, without employing acceleration forces or bias, we captured a complete one-base-pair stepping cycle of the protein tracking along major groove of DNA with a homogeneous poly-adenosine sequence, as individual hydrogen bonds break and reform at the protein-DNA binding interface. Further DNA-groove tracking motions of the protein forward or backward, with occasional sliding as well as strand crossing to minor groove of DNA, were also captured. The processive diffusion of WRKY along DNA has been further sampled via coarse-grained MD simulations. The study thus provides structural dynamics details on diffusion of a small TF domain protein, suggests how the protein approaches a specific recognition site on DNA, and supports further high-precision experimental detection. The stochastic movements revealed in the TF diffusion also provide general clues about how other protein walkers step and slide along DNA.

Pathway-based classification of glioblastoma uncovers a mitochondrial subtype with therapeutic vulnerabilities.

The transcriptomic classification of glioblastoma (GBM) has failed to predict survival and therapeutic vulnerabilities. A computational approach for unbiased identification of core biological traits of single cells and bulk tumors uncovered four tumor cell states and GBM subtypes distributed along neurodevelopmental and metabolic axes, classified as proliferative/progenitor, neuronal, mitochondrial and glycolytic/plurimetabolic. Each subtype was enriched with biologically coherent multiomic features. Mitochondrial GBM was associated with the most favorable clinical outcome. It relied exclusively on oxidative phosphorylation for energy production, whereas the glycolytic/plurimetabolic subtype was sustained by aerobic glycolysis and amino acid and lipid metabolism. Deletion of the glucose-proton symporter SLC45A1 was the truncal alteration most significantly associated with mitochondrial GBM, and the reintroduction of SLC45A1 in mitochondrial glioma cells induced acidification and loss of fitness. Mitochondrial, but not glycolytic/plurimetabolic, GBM exhibited marked vulnerability to inhibitors of oxidative phosphorylation. The pathway-based classification of GBM informs survival and enables precision targeting of cancer metabolism.

A map of tumor-host interactions in glioma at single-cell resolution.

BACKGROUND: Single-cell RNA sequencing is the reference technique for characterizing the heterogeneity of the tumor microenvironment. The composition of the various cell types making up the microenvironment can significantly affect the way in which the immune system activates cancer rejection mechanisms. Understanding the cross-talk signals between immune cells and cancer cells is of fundamental importance for the identification of immuno-oncology therapeutic targets. RESULTS: We present a novel method, single-cell Tumor-Host Interaction tool (scTHI), to identify significantly activated ligand-receptor interactions across clusters of cells from single-cell RNA sequencing data. We apply our approach to uncover the ligand-receptor interactions in glioma using 6 publicly available human glioma datasets encompassing 57,060 gene expression profiles from 71 patients. By leveraging this large-scale collection we show that unexpected cross-talk partners are highly conserved across different datasets in the majority of the tumor samples. This suggests that shared cross-talk mechanisms exist in glioma. CONCLUSIONS: Our results provide a complete map of the active tumor-host interaction pairs in glioma that can be therapeutically exploited to reduce the immunosuppressive action of the microenvironment in brain tumor.

Mn2+ Directly Activates cGAS and Structural Analysis Suggests Mn2+ Induces a Noncanonical Catalytic Synthesis of 2'3'-cGAMP.

DNA binding allosterically activates the cytosolic DNA sensor cGAS (cyclic GMP-AMP [cGAMP] synthase) to synthesize 2'3'-cGAMP, using Mg2+ as the metal cofactor that catalyzes two nucleotidyl-transferring reactions. We previously found that Mn2+ potentiates cGAS activation, but the underlying mechanism remains unclear. Here, we report that Mn2+ directly activates cGAS. Structural analysis reveals that Mn2+-activated cGAS undergoes globally similar conformational changes to DNA-activated cGAS but forms a unique η1 helix to widen the catalytic pocket, allowing substrate entry and cGAMP synthesis. Strikingly, in Mn2+-activated cGAS, the linear intermediates pppGpG and pGpA take an inverted orientation in the active pocket, suggesting a noncanonical but accelerated cGAMP cyclization without substrate flip-over. Moreover, unlike the octahedral coordination around Mg2+, the two catalytic Mn2+ are coordinated by triphosphate moiety of the inverted substrate, independent of the catalytic triad residues. Our findings thus uncover Mn2+ as a cGAS activator that initiates noncanonical 2'3'-cGAMP synthesis.

Structural insights into immunoglobulin M.

Immunoglobulin M (IgM) plays a pivotal role in both humoral and mucosal immunity. Its assembly and transport depend on the joining chain (J-chain) and the polymeric immunoglobulin receptor (pIgR), but the underlying molecular mechanisms of these processes are unclear. We report a cryo-electron microscopy structure of the Fc region of human IgM in complex with the J-chain and pIgR ectodomain. The IgM-Fc pentamer is formed asymmetrically, resembling a hexagon with a missing triangle. The tailpieces of IgM-Fc pack into an amyloid-like structure to stabilize the pentamer. The J-chain caps the tailpiece assembly and bridges the interaction between IgM-Fc and the polymeric immunoglobulin receptor, which undergoes a large conformational change to engage the IgM-J complex. These results provide a structural basis for the function of IgM.

Surveying brain tumor heterogeneity by single-cell RNA-sequencing of multi-sector biopsies.

Brain tumors are among the most challenging human tumors for which the mechanisms driving progression and heterogeneity remain poorly understood. We combined single-cell RNA-seq with multi-sector biopsies to sample and analyze single-cell expression profiles of gliomas from 13 Chinese patients. After classifying individual cells, we generated a spatial and temporal landscape of glioma that revealed the patterns of invasion between the different sub-regions of gliomas. We also used single-cell inferred copy number variations and pseudotime trajectories to inform on the crucial branches that dominate tumor progression. The dynamic cell components of the multi-region biopsy analysis allowed us to spatially deconvolute with unprecedented accuracy the transcriptomic features of the core and those of the periphery of glioma at single-cell level. Through this rich and geographically detailed dataset, we were also able to characterize and construct the chemokine and chemokine receptor interactions that exist among different tumor and non-tumor cells. This study provides the first spatial-level analysis of the cellular states that characterize human gliomas. It also presents an initial molecular map of the cross-talks between glioma cells and the surrounding microenvironment with single-cell resolution.

HeLa-CCL2 cell heterogeneity studied by single-cell DNA and RNA sequencing.

The HeLa cells are the earliest and mostly used laboratory human cells for biomedical particularly cancer research. They were derived from a patient's cervical cancerous tissue, and are known for their heterogeneous cellular origin and variable genomic landscapes. Single-cell sequencing techniques with faithful linear and uniformly amplified genomes (DNA) and transcriptomes (RNA) may facilitate the study of cellular differences at the individual cell level. In this work, we have performed single-cell DNA and RNA sequencing with HeLa-CCL2 cells to study their heterogeneity. We have studied the complexity of copy number variations (CNVs) of HeLa-CCL2 genome at the single cell level, and revealed the transcriptomic heterogeneity of HeLa-CCL2. We also analyzed the relationship between genome and transcriptome at the single-cell level, and found overall correlation between CNV and transcriptome expression patterns. Finally, we concluded that although single-cell sequencing techniques are applicable to study heterogeneous cells such as HeLa-CCL2, the data analyses need to be more careful and well controlled.

Tumor Volume Is Better Than Diameter for Predicting the Prognosis of Patients with Early-Stage Non-small Cell Lung Cancer.

BACKGROUND: This study aimed to investigate whether tumor volume (TV) is better than diameter for predicting the prognosis of patients with early-stage non-small cell lung cancer (NSCLC) after complete resection. METHODS: This study retrospectively reviewed the clinicopathologic characteristics of 274 patients with early-stage NSCLC who had received pretreatment computed tomography (CT) scans and complete resection. TV was semi-automatically measured from CT scans using an imaging software program. The optimal cutoff of TV was determined by X-tile software. Disease-free survival (DFS) and overall survival (OS) were assessed by the Kaplan-Meier method. The prognostic significance of TV and other variables was assessed by Cox proportional hazards regression analysis. RESULTS: Using 3.046 cm3 and 8.078 cm3 as optimal cutoff values of TV, the patients were separated into three groups. A larger TV was significantly associated with poor DFS and OS in the multivariable analysis. Kaplan-Meier curves of DFS and OS showed significant differences on the basis of TV among patients with stage 1a disease, greatest tumor diameter (GTD) of 2 cm or smaller, and GTD of 2-3 cm, respectively. Using two TV cutoff points, three categories of TV were created. In 54 cases (19.7%), patients migrated from the GTD categories of 2 cm or smaller, 2-3 cm, and larger than 3 cm into the TV categories of 3.046 cm3 or smaller, 3.046-8.078 cm3, and larger than 8.078 cm3. CONCLUSION: TV is an independent prognostic factor of DFS and OS for early-stage NSCLC. The findings show that TV is better than GTD for predicting the prognosis of patients with early-stage NSCLC.

Mutation Profile of Resected EGFR-Mutated Lung Adenocarcinoma by Next-Generation Sequencing.

BACKGROUND: The efficacy of adjuvant targeted therapy for operable lung cancer is still under debate. Comprehensive genetic profiling is needed for detecting co-mutations in resected epidermal growth factor receptor (EGFR)-mutated lung adenocarcinoma (ADC), which may interfere the efficacy of adjuvant tyrosine kinase inhibitor (TKI) treatment. MATERIALS AND METHODS: Mutation profiling of 416 cancer-relevant genes was conducted for 139 resected stage I-IIIa lung ADCs with EGFR mutations using targeted next-generation sequencing. Co-mutation profiles were systematically analyzed. RESULTS: Rare EGFR alterations other than exon 19 deletion and L858R, such as L861Q (∼3%) and G719A (∼2%), were identified at low frequencies. Approximately 10% of patients had mutations in EGFR exon 20 that could confer resistance to first-generation TKIs. Ninety-one percent of patients harbored at least one co-mutation in addition to the major EGFR mutation. TP53 was the top mutated gene and was found more frequently mutated at later stage. Markedly, NF1 mutations were found only in stage II-III ADCs. Conversely, RB1 mutations were more frequent in stage I ADCs, whereas APC mutations were observed exclusively in this group. Thirty-four percent of patients with EGFR TKI-sensitizing mutations had genetic alterations involving EGFR downstream effectors or bypass pathways that could affect the response to EGFR TKIs, such as PIK3CA, BRCA1, and NOTCH1. CONCLUSION: Operable lung ADCs with EGFR TKI-sensitizing mutations are associated with a high proportion of co-mutations. Mutation profiling of these resected tumors could facilitate in determining the applicability and efficacy of adjuvant EGFR TKI therapeutic strategy. IMPLICATIONS FOR PRACTICE: The efficacy of adjuvant epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor (TKI) therapy for lung cancer harboring EGFR mutation after surgical resection is still under debate. Next-generation sequencing of 416 cancer-relevant genes in 139 resected lung cancers revealed the co-mutational landscape with background EGFR mutation. Notably, the study identified potential EGFR TKI-resistant mutations in 34.71% of patients with a drug-sensitizing EGFR mutation and who were naive in terms of targeted therapy. A comprehensive mutation profiling of these resected tumors could facilitate in determining the applicability and efficacy of adjuvant EGFR TKI therapeutic strategy for these patients.

Novel Mechanism for Cyclic Dinucleotide Degradation Revealed by Structural Studies of Vibrio Phosphodiesterase V-cGAP3.

3'3'-cyclic GMP-AMP (3'3'-cGAMP) belongs to a family of the bacterial secondary messenger cyclic dinucleotides. It was first discovered in the Vibrio cholerae seventh pandemic strains and is involved in efficient intestinal colonization and chemotaxis regulation. Phosphodiesterases (PDEs) that degrade 3'3'-cGAMP play important regulatory roles in the relevant signaling pathways, and a previous study has identified three PDEs in V. cholerae, namely, V-cGAP1, V-cGAP2, and V-cGAP3, functioning in 3'3'-cGAMP degradation. We report the crystal structure, biochemical, and structural analyses of V-cGAP3, providing a foundation for understanding the mechanism of 3'3'-cGAMP degradation and regulation in general. Our crystal and molecular dynamic (MD)-simulated structures revealed that V-cGAP3 contains tandem HD-GYP domains within its N- and C-terminal domains, with similar three-dimensional topologies despite their low-sequence identity. Biochemical and structural analyses showed that the N-terminal domain plays a mechanism of positive regulation for the catalytic C-terminal domain. We also demonstrated that the other homologous Vibrio PDEs, V-cGAP1/2, likely function via a similar mechanism.

Loss of MED12 Induces Tumor Dormancy in Human Epithelial Ovarian Cancer via Downregulation of EGFR.

A high rate of disease relapse makes epithelial ovarian cancer (EOC) the leading cause of death among all gynecologic malignancies. These relapses are often due to tumor dormancy. Here we identify the RNA polymerase II transcriptional mediator subunit 12 (MED12) as an important molecular regulator of tumor dormancy. MED12 knockout (KO) induced dormancy of EOC cells in vitro and in vivo, and microarray analysis showed that MED12 KO decreased expression of EGFR. Restoration of EGFR expression in MED12 KO cells restored proliferation. Additionally, MED12 bound to the promoter of EGFR, and correlation studies showed that MED12 expression positively correlated with EGFR expression in EOC patient samples. Clinical data demonstrated that chemotherapy-resistant patients expressed lower levels of MED12 compared with responsive patients. Overall, our data show that MED12 plays an important role in regulating dormancy of EOC through regulation of EGFR.Significance: MED12 is identified as a novel, important regulator of tumor dormancy in human ovarian cancer. Cancer Res; 78(13); 3532-43. ©2018 AACR.

Polo-like kinase 1 (PLK1)-dependent phosphorylation of methylenetetrahydrofolate reductase (MTHFR) regulates replication via histone methylation.

Methylenetetrahydrofolate reductase (MTHFR) is a key enzyme regulating the folate cycle and its genetic variations have been associated with various human diseases. Previously we identified that MTHFR is phosphorylated by cyclin-dependent kinase 1 (CDK1) at T34 and MTHFR underlies heterochromatin maintenance marked by H3K9me3 levels. Herein we demonstrate that pT34 creates a binding motif that docks MTHFR to the polo-binding domain (PBD) of polo-like kinase 1 (PLK1), a fundamental kinase that orchestrates many cell cycle events. We show that PLK1 phosphorylates MTHFR at T549 in vitro and in vivo. Further, we uncovered a role of MTHFR in replication. First, MTHFR depletion increased the fraction of cells in S phase. This defect could not be rescued by siRNA resistant plasmids harboring T549A, but could be restored by overproduction of Suv4-20H2, the H4K20 methyltransferase. Moreover, siMTHFR attenuated H4K20me3 levels, which could be rescued by Suv4-20H2 overproduction. More importantly, we also investigated MTHFR-E429A, the protein product of an MTHFR single nucleotide variant. MTHFR-E429A overexpression also increased S phase cells and decreased H4K20me3 levels, and it is linked to a poor glioma prognosis in the Chinese population. Collectively, we have unveiled a vital role of PLK1-dependent phosphorylation of MTHFR in replication via histone methylation, and implicate folate metabolism with glioma.

The prognostic impact of tumor volume on stage I non-small cell lung cancer.

OBJECTIVES: The purpose of this study was to investigate the prognostic impact of tumor volume (TV) on patients with stage I non-small cell lung cancer (NSCLC) after complete resection. MATERIALS AND METHODS: We retrospectively reviewed the clinicopathological characteristics of 274 patients with stage I NSCLC who had received preoperative chest computed tomography (CT) scans and complete resection. TV was semi-automatically measured from chest CT scans by using an imaging software program. The optimal cutoff values of TV were determined by X-tile software. Disease-free survival (DFS) and overall survival (OS) were compared using Kaplan-Meier analysis. Univariate and multivariate analyses were performed to identify risk factors for DFS and OS. RESULTS: By using 3.046cm3 and 8.078cm3 as two optimal cutoff values of TV, the patients were separated into three groups. The 5-year DFS and OS for patients with TV≤3.046cm3, 3.046-8.078cm3, and>8.078cm3 were 88.0%, 73.6%, and 62.1%, respectively (P<0.001), and 91.4%, 84.5%, and 73.3%, respectively (p<0.001). Multivariate analysis showed that age and TV were independent factors associated with DFS. Sex, age, histology, visceral pleural invasion, and TV were independent factors associated with OS. Stage Ia patients might be separated into three groups on the basis of TV with significantly different DFS and OS. Patients with tumor diameter≤2cm and 2-3cm were also stratified into two groups with significantly different DFS and OS on the basis of TV, respectively. CONCLUSION: TV is an independent risk factor for DFS and OS for stage I NSCLC after complete resection. TV might provide additional prognostic information over tumor diameter in patients with stage I NSCLC.

Effect of body mass index on survival of patients with stage I non-small cell lung cancer.

BACKGROUND: Body mass index (BMI) has a U-shaped association with lung cancer risk. However, the effect of BMI on prognosis is controversial. This retrospective study aimed to investigate the effect of BMI on the survival of patients with stage I non-small cell lung cancer (NSCLC) after surgical resection. METHODS: In total, 624 consecutive stage I NSCLC patients who underwent radical resection were classified into four groups according to their BMI: underweight (BMI < 18.5 kg/m2), normal weight (BMI = 18.5-22.4 kg/m2), overweight (BMI = 22.5-28.0 kg/m2), and obese (BMI > 28.0 kg/m2). The effect of BMI on progression-free survival (PFS) and overall survival (OS) was estimated using the Kaplan-Meier method and Cox proportional hazards model. Postoperative complications in each group were analyzed using the Chi square test or Fisher's exact test. RESULTS: A univariate analysis showed that PFS and OS were longer in the overweight group than in other groups (both P < 0.05). A multivariate analysis showed that OS was longer in the overweight group than in other groups (compared with the other three groups in combination: hazard ratio [HR] = 1.87, 95% confidence interval [CI] 1.30-2.68, P = 0.003; compared with the underweight group: HR = 2.24, 95% CI 1.18-4.25, P = 0.013; compared with the normal weight group: HR = 1.58, 95% CI 1.07-2.33, P = 0.022; compared with the obese group: HR = 2.87, 95% CI 1.48-5.59, P = 0.002), but PFS was similar among the groups (HR = 1.28, 95% CI 0.97-1.68, P = 0.080). A subgroup analysis showed an association between being overweight and prolonged OS in patients at stage T1a (P = 0.024), T1b (P = 0.051), and T2a (P = 0.02), as well as in patients with a non-smoking history (P = 0.001). Overweight patients had lower rates of postoperative complications, such as respiratory failure (compared with the underweight and obese groups: P = 0.014), myocardial infarction (compared with the obese group: P = 0.033), and perioperative death (compared with the other three groups: P = 0.016). CONCLUSIONS: Preoperative BMI is an independent prognostic factor for stage I NSCLC patients after resection, with overweight patients having a favorable prognosis.

The effects of cytosine methylation on general transcription factors.

DNA methylation on CpG sites is the most common epigenetic modification. Recently, methylation in a non-CpG context was found to occur widely on genomic DNA. Moreover, methylation of non-CpG sites is a highly controlled process, and its level may vary during cellular development. To study non-CpG methylation effects on DNA/protein interactions, we have chosen three human transcription factors (TFs): glucocorticoid receptor (GR), brain and muscle ARNT-like 1 (BMAL1) - circadian locomotor output cycles kaput (CLOCK) and estrogen receptor (ER) with methylated or unmethylated DNA binding sequences, using single-molecule and isothermal titration calorimetry assays. The results demonstrated that these TFs interact with methylated DNA with different effects compared with their cognate DNA sequences. The effects of non-CpG methylation on transcriptional regulation were validated by cell-based luciferase assay at protein level. The mechanisms of non-CpG methylation influencing DNA-protein interactions were investigated by crystallographic analyses and molecular dynamics simulation. With BisChIP-seq assays in HEK-293T cells, we found that GR can recognize highly methylated sites within chromatin in cells. Therefore, we conclude that non-CpG methylation of DNA can provide a mechanism for regulating gene expression through directly affecting the binding of TFs.

Younger women have a better prognosis among patients with esophageal squamous cell carcinoma after esophagectomy.

BACKGROUND: Epidemiological studies have suggested a potential role for sex hormones in esophageal squamous cell carcinoma (ESCC) etiology. However, the significance of gender as an independent prognostic factor remains uncertain. Our retrospective study was designed to investigate the prognostic role of gender for survival after esophagectomy. METHODS: Data from 674 patients with ESCC who underwent surgical resection were retrospectively analysed. Age 55 years was selected as a surrogate for menopause. Patients were assigned to age-gender groups: A (female younger than age 55 years), B (female age 55 years and older), C (male younger than age 55 years) and D (male age 55 years and older). Univariate and multivariate analyses were performed to identify prognostic factors for survival. RESULTS: There were 520 males and 154 females with median age was 58 years. The 1-, 3-, 5-year survival rates of group A, B, C and D were 93.6%, 70.2%, 61.7% vs. 86.9%, 47.7%, 40.2% vs. 77.8%, 43.9%, 37.0% and 80.3%, 47.9%, 36.6%, respectively (P=0.003). Multivariate analysis suggested that age-gender groups and pTNM staging were independent prognostic factors. CONCLUSIONS: Among patients with ESCC after esophagectomy, women younger than age 55 years attained a favorable prognosis.

Rat and human STINGs profile similarly towards anticancer/antiviral compounds.

Cyclic dinucleotides (CDNs) and antitumor/antiviral agents (DMXAA and CMA) trigger STING-dependent innate immunity activation. Accumulative evidences have showed that DMXAA and CMA selectively activate mouse, but not human STING signaling. The mechanism underlying this species selectivity remains poorly understood. In this report, we have shown that human and rat STINGs display more similar signaling profiles toward DMXAA and CMA than that of human and mouse STINGs, suggesting that rat is more suitable for preclinical testing of STING-targeted drugs. We have also determined the crystal structures of both apo rat STING and its complex with cyclic GMP-AMP with 2'5' and 3'5' phosphodiester linkage (2'3'-cGAMP), a human endogenous CDN. Structure-guided biochemical analysis also revealed the functional importance of the connecting loop (A140-N152) between membrane and cytosolic domains in STING activation. Taken together, these findings reveal that rat STING is more closely related to human STING in terms of substrate preference, serving as a foundation for the development of STING-targeted drugs.