Inferring the impact of climate upon the transmission of SARS-CoV-2 has been confounded by variability in testing, unknown disease introduction rates, and changing weather. Here we present a data ...model that accounts for dynamic testing rates and variations in disease introduction rates. We apply this model to data from Colombia, whose varied and seasonless climate, central port of entry, and swift, centralized response to the COVID-19 pandemic present an opportune environment for assessing the impact of climate factors on the spread of COVID-19. We observe strong attenuation of transmission in climates with sustained daily temperatures above 30 degrees Celsius and simultaneous mean relative humidity below 78%, with outbreaks occurring at high humidity even where the temperature is high. We hypothesize that temperature and relative humidity comodulate the infectivity of SARS-CoV-2 within respiratory droplets.
Antibodies that target the immune checkpoint receptor programmed cell death protein 1 (PD-1) have resulted in prolonged and beneficial responses toward a variety of human cancers. However, anti-PD-1 ...therapy in some patients provides no benefit and/or results in adverse side effects. The factors that determine whether patients will be drug sensitive or resistant are not fully understood; therefore, genomic assessment of exceptional responders can provide important insight into patient response. Here, we identified a patient with endometrial cancer who had an exceptional response to the anti-PD-1 antibody pembrolizumab. Clinical grade targeted genomic profiling of a pretreatment tumor sample from this individual identified a mutation in DNA polymerase epsilon (POLE) that associated with an ultramutator phenotype. Analysis of The Cancer Genome Atlas (TCGA) revealed that the presence of POLE mutation associates with high mutational burden and elevated expression of several immune checkpoint genes. Together, these data suggest that cancers harboring POLE mutations are good candidates for immune checkpoint inhibitor therapy.
Glutaminase (GLS) isoenzymes GLS1 and GLS2 are key enzymes for glutamine metabolism. Interestingly, GLS1 and GLS2 display contrasting functions in tumorigenesis with elusive mechanism; GLS1 promotes ...tumorigenesis, whereas GLS2 exhibits a tumor-suppressive function. In this study, we found that GLS2 but not GLS1 binds to small GTPase Rac1 and inhibits its interaction with Rac1 activators guanine-nucleotide exchange factors, which in turn inhibits Rac1 to suppress cancer metastasis. This function of GLS2 is independent of GLS2 glutaminase activity. Furthermore, decreased GLS2 expression is associated with enhanced metastasis in human cancer. As a p53 target, GLS2 mediates p53's function in metastasis suppression through inhibiting Rac1. In summary, our results reveal that GLS2 is a novel negative regulator of Rac1, and uncover a novel function and mechanism whereby GLS2 suppresses metastasis. Our results also elucidate a novel mechanism that contributes to the contrasting functions of GLS1 and GLS2 in tumorigenesis.
Response to immune checkpoint therapy can be associated with a high mutation burden, but other mechanisms are also likely to be important. We identified a patient with metastatic gastric cancer with ...meaningful clinical benefit from treatment with the anti-programmed death-ligand 1 (PD-L1) antibody avelumab. This tumor showed no evidence of high mutation burden or mismatch repair defect but was strongly positive for presence of Epstein-Barr virus (EBV) encoded RNA. Analysis of The Cancer Genome Atlas gastric cancer data (25 EBV+, 80 microsatellite-instable MSI, 310 microsatellite-stable MSS) showed that EBV-positive tumors were MSS. Two-sided Wilcoxon rank-sum tests showed that: 1) EBV-positive tumors had low mutation burden (median = 2.07 vs 3.13 in log10 scale, P < 10-12) but stronger evidence of immune infiltration (median ImmuneScore 2212 vs 1295, P < 10-4; log2 fold-change of CD8A = 1.85, P < 10-6) compared with MSI tumors, and 2) EBV-positive tumors had higher expression of immune checkpoint pathway (PD-1, CTLA-4 pathway) genes in RNA-seq data (log2 fold-changes: PD-1 = 1.85, PD-L1 = 1.93, PD-L2 = 1.50, CTLA-4 = 1.31, CD80 = 0.89, CD86 = 1.31, P < 10-4 each), and higher lymphocytic infiltration by histology (median tumor-infiltrating lymphocyte score = 3 vs 2, P < .001) compared with MSS tumors. These data suggest that EBV-positive low-mutation burden gastric cancers are a subset of MSS gastric cancers that may respond to immune checkpoint therapy.
MicroRNA regulate mRNA levels in a tissue specific way, either by inducing degradation of the transcript or by inhibiting translation or transcription. Putative mRNA targets of microRNA identified ...from seed sequence matches are available in many databases. However, such matches have a high false positive rate and cannot identify tissue specificity of regulation.
We describe a simple method to identify direct mRNA targets of microRNA dysregulated in cancers from expression level measurements in patient matched tumor/normal samples. The word "direct" is used here in a strict sense to: a) represent mRNA which have an exact seed sequence match to the microRNA in their 3'UTR, b) the seed sequence match is strictly conserved across mouse, human, rat and dog genomes, c) the mRNA and microRNA expression levels can distinguish tumor from normal with high significance and d) the microRNA/mRNA expression levels are strongly and significantly anti-correlated in tumor and/or normal samples. We apply and validate the method using clear cell Renal Cell Carcinoma (ccRCC) and matched normal kidney samples, limiting our analysis to mRNA targets which undergo degradation of the mRNA transcript because of a perfect seed sequence match. Dysregulated microRNA and mRNA are first identified by comparing their expression levels in tumor vs normal samples. Putative dysregulated microRNA/mRNA pairs are identified from these using seed sequence matches, requiring that the seed sequence be conserved in human/dog/rat/mouse genomes. These are further pruned by requiring a strong anti-correlation signature in tumor and/or normal samples. The method revealed many new regulations in ccRCC. For instance, loss of miR-149, miR-200c and mir-141 causes gain of function of oncogenes (KCNMA1, LOX), VEGFA and SEMA6A respectively and increased levels of miR-142-3p, miR-185, mir-34a, miR-224, miR-21 cause loss of function of tumor suppressors LRRC2, PTPN13, SFRP1, ERBB4, and (SLC12A1, TCF21) respectively. We also found strong anti-correlation between VEGFA and the miR-200 family of microRNA: miR-200a*, 200b, 200c and miR-141. Several identified microRNA/mRNA pairs were validated on an independent set of matched ccRCC/normal samples. The regulation of SEMA6A by miR-141 was verified by a transfection assay.
We describe a simple and reliable method to identify direct gene targets of microRNA in any cancer. The constraints we impose (strong dysregulation signature for microRNA and mRNA levels between tumor/normal samples, evolutionary conservation of seed sequence and strong anti-correlation of expression levels) remove spurious matches and identify a subset of robust, tissue specific, functional mRNA targets of dysregulated microRNA.
Objectives To identify a robust panel of microRNA signatures that can classify tumor from normal kidney using microRNA expression levels. Mounting evidence suggests that microRNAs are key players in ...essential cellular processes and that their expression pattern can serve as diagnostic biomarkers for cancerous tissues. Methods We selected 28 clear-cell type human renal cell carcinoma (ccRCC), samples from patient-matched specimens to perform high-throughput, quantitative real-time polymerase chain reaction analysis of microRNA expression levels. The data were subjected to rigorous statistical analyses and hierarchical clustering to produce a discrete set of microRNAs that can robustly distinguish ccRCC from their patient-matched normal kidney tissue samples with high confidence. Results Thirty-five microRNAs were found that can robustly distinguish ccRCC from their patient-matched normal kidney tissue samples with high confidence. Among this set of 35 signature microRNAs, 26 were found to be consistently downregulated and 9 consistently upregulated in ccRCC relative to normal kidney samples. Two microRNAs, namely, MiR - 155 and miR - 21 , commonly found to be upregulated in other cancers, and miR - 210 , induced by hypoxia, were also identified as overexpressed in ccRCC in our study. MicroRNAs identified as downregulated in our study can be correlated to common chromosome deletions in ccRCC. Conclusions Our analysis is a comprehensive, statistically relevant study that identifies the microRNAs dysregulated in ccRCC, which can serve as the basis of molecular markers for diagnosis.
To improve cancer therapy, it is critical to target metastasizing cells. Circulating tumor cells (CTCs) are rare cells found in the blood of patients with solid tumors and may play a key role in ...cancer dissemination. Uncovering CTC phenotypes offers a potential avenue to inform treatment. However, CTC transcriptional profiling is limited by leukocyte contamination; an approach to surmount this problem is single cell analysis. Here we demonstrate feasibility of performing high dimensional single CTC profiling, providing early insight into CTC heterogeneity and allowing comparisons to breast cancer cell lines widely used for drug discovery.
We purified CTCs using the MagSweeper, an immunomagnetic enrichment device that isolates live tumor cells from unfractionated blood. CTCs that met stringent criteria for further analysis were obtained from 70% (14/20) of primary and 70% (21/30) of metastatic breast cancer patients; none were captured from patients with non-epithelial cancer (n = 20) or healthy subjects (n = 25). Microfluidic-based single cell transcriptional profiling of 87 cancer-associated and reference genes showed heterogeneity among individual CTCs, separating them into two major subgroups, based on 31 highly expressed genes. In contrast, single cells from seven breast cancer cell lines were tightly clustered together by sample ID and ER status. CTC profiles were distinct from those of cancer cell lines, questioning the suitability of such lines for drug discovery efforts for late stage cancer therapy.
For the first time, we directly measured high dimensional gene expression in individual CTCs without the common practice of pooling such cells. Elevated transcript levels of genes associated with metastasis NPTN, S100A4, S100A9, and with epithelial mesenchymal transition: VIM, TGFß1, ZEB2, FOXC1, CXCR4, were striking compared to cell lines. Our findings demonstrate that profiling CTCs on a cell-by-cell basis is possible and may facilitate the application of 'liquid biopsies' to better model drug discovery.
Esophageal adenocarcinoma (EAC) is strongly associated with Barrett’s esophagus (BE), a pre-malignant condition resulting from gastric reflux. Esophageal squamous cell carcinoma (ESCC), the other ...major subtype of esophageal cancer, shows strong association with smoking and alcohol intake and no association with gastric reflux. In this study, we constructed and validated gene expression signatures of EAC vs. ESCC tumors using publicly available datasets, and subsequently assessed the enrichment levels of these signatures in commonly used EAC and ESCC cell lines, normal esophageal tissues and normal esophageal cell lines, and primary BE tissues and BE cell lines. We found that unlike ESCC cell lines which were quite similar to primary ESCC tumors, EAC cell lines were considerably different from primary EAC tumors but still more similar to EAC tumors than ESCC tumors, as the genes up in EAC vs. ESCC (EAChi) had considerably lower expression in EAC cell lines than EAC tumors. However, more surprisingly, unlike various normal cell lines (EPC2, Het-1A) which were very similar to various tissues from normal esophagus, BE cell lines (BAR-T, CP-A) were extremely different from primary BE tissues, as BE cell lines had substantially lower levels of EAChi and substantially higher levels of ESCChi gene expression. This ESCC-like profile of the BAR-T remained unaltered even after prolonged exposure to an acidic bile mixture in vitro resulting in malignant transformation of this cell line. However, primary BE tissues had EAC-like gene expression profiles as expected. Only one EAC case from the Cancer Genome Atlas resembled BE cell lines, and while it had the clinical profile and some mutational features of EAC, it had some mutational features, the copy number alteration profile, and the gene expression profile of ESCC instead. These incomprehensible changes in gene expression patterns may result in ambiguous changes in the phenotype and warrants careful evaluation to inform selection of appropriate in vitro tools for future studies on esophageal adenocarcinoma.
It is well known that the dinucleotide CpG is under-represented in the genomic DNA of many vertebrates. This is commonly thought to be due to the methylation of cytosine residues in this dinucleotide ...and the corresponding high rate of deamination of 5-methycytosine, which lowers the frequency of this dinucleotide in DNA. Surprisingly, many single-stranded RNA viruses that replicate in these vertebrate hosts also have a very low presence of CpG dinucleotides in their genomes. Viruses are obligate intracellular parasites and the evolution of a virus is inexorably linked to the nature and fate of its host. One therefore expects that virus and host genomes should have common features. In this work, we compare evolutionary patterns in the genomes of ssRNA viruses and their hosts. In particular, we have analyzed dinucleotide patterns and found that the same patterns are pervasively over- or under-represented in many RNA viruses and their hosts suggesting that many RNA viruses evolve by mimicking some of the features of their host's genes (DNA) and likely also their corresponding mRNAs. When a virus crosses a species barrier into a different host, the pressure to replicate, survive and adapt, leaves a footprint in dinucleotide frequencies. For instance, since human genes seem to be under higher pressure to eliminate CpG dinucleotide motifs than avian genes, this pressure might be reflected in the genomes of human viruses (DNA and RNA viruses) when compared to those of the same viruses replicating in avian hosts. To test this idea we have analyzed the evolution of the influenza virus since 1918. We find that the influenza A virus, which originated from an avian reservoir and has been replicating in humans over many generations, evolves in a direction strongly selected to reduce the frequency of CpG dinucleotides in its genome. Consistent with this observation, we find that the influenza B virus, which has spent much more time in the human population, has adapted to its human host and exhibits an extremely low CpG dinucleotide content. We believe that these observations directly show that the evolution of RNA viral genomes can be shaped by pressures observed in the host genome. As a possible explanation, we suggest that the strong selection pressures acting on these RNA viruses are most likely related to the innate immune response and to nucleotide motifs in the host DNA and RNAs.
Immune checkpoint blockade leads to unprecedented responses in many cancers. Although currently available agents mostly target the PD-1 and CTLA-4 pathways, agents targeting the immune checkpoint ...protein LAG-3 are under active clinical development, and early clinical data show that
expression is a biomarker of response to LAG-3 blockade. To determine which cancers may benefit most from LAG-3 blockade, we performed a pan-cancer analysis of The Cancer Genome Atlas dataset to identify genomic and immunologic correlates of
expression. High mutation burden, and expression of exogenous virus (EBV, HPV) or endogenous retrovirus (
), were associated with overexpression of
in multiple cancers. Although CD8
T-cell marker (
) and
were strongly co-expressed with each other and with
in most cancers, there were three notable exceptions: HPV+ head-neck squamous cell cancer, renal cell cancer, and glioblastoma. These results may have important implications for guiding development clinical trials of LAG-3 blockade.