Erika Thompson

Landscape of DNA Virus Associations across Human Malignant Cancers: Analysis of 3,775 Cases Using RNA-Seq

ABSTRACT Elucidation of tumor-DNA virus associations in many cancer types has enhanced our knowledge of fundamental oncogenesis mechanisms and provided a basis for cancer prevention initiatives. RNA-Seq is a novel tool to comprehensively assess such associations. We interrogated RNA-Seq data from 3,775 malignant neoplasms in The Cancer Genome Atlas database for the presence of viral sequences. Viral integration sites were also detected in expressed transcripts using a novel approach. The detection capacity of RNA-Seq was compared to available clinical laboratory data. Human papillomavirus (HPV) transcripts were detected using RNA-Seq analysis in head-and-neck squamous cell carcinoma, uterine endometrioid carcinoma, and squamous cell carcinoma of the lung. Detection of HPV by RNA-Seq correlated with detection by in situ hybridization and immunohistochemistry in squamous cell carcinoma tumors of the head and neck. Hepatitis B virus and Epstein-Barr virus (EBV) were detected using RNA-Seq in hepatocellular carcinoma and gastric carcinoma tumors, respectively. Integration sites of viral genes and oncogenes were detected in cancers harboring HPV or hepatitis B virus but not in EBV-positive gastric carcinoma. Integration sites of expressed viral transcripts frequently involved known coding areas of the host genome. No DNA virus transcripts were detected in acute myeloid leukemia, cutaneous melanoma, low- and high-grade gliomas of the brain, and adenocarcinomas of the breast, colon and rectum, lung, prostate, ovary, kidney, and thyroid. In conclusion, this study provides a large-scale overview of the landscape of DNA viruses in human malignant cancers. While further validation is necessary for specific cancer types, our findings highlight the utility of RNA-Seq in detecting tumor-associated DNA viruses and identifying viral integration sites that may unravel novel mechanisms of cancer pathogenesis.

Gene program for cardiac cell survival induced by transient ischemia in conscious pigs

Therapy for ischemic heart disease has been directed traditionally at limiting cell necrosis. We determined by genome profiling whether ischemic myocardium can trigger a genetic program promoting cardiac cell survival, which would be a novel and potentially equally important mechanism of salvage. Although cardiac genomics is usually performed in rodents, we used a swine model of ischemia/reperfusion followed by ventricular dysfunction (stunning), which more closely resembles clinical conditions. Gene expression profiles were compared by subtractive hybridization between ischemic and normal tissue of the same hearts. About one-third (23/74) of the nuclear-encoded genes that were up-regulated in ischemic myocardium participate in survival mechanisms (inhibition of apoptosis, cytoprotection, cell growth, and stimulation of translation). The specificity of this response was confirmed by Northern blot and quantitative PCR. Unexpectedly, this program also included genes not previously described in cardiomyocytes. Up-regulation of survival genes was more profound in subendocardium over subepicardium, reflecting that this response in stunned myocardium was proportional to the severity of the ischemic insult. Thus, in a swine model that recapitulates human heart disease, nonlethal ischemia activates a genomic program of cell survival that relates to the time course of myocardial stunning and differs transmurally in relation to ischemic stress, which induced the stunning. Understanding the genes up-regulated during myocardial stunning, including those not previously described in the heart, and developing strategies that activate this program may open new avenues for therapy in ischemic heart disease.

VirusSeq: software to identify viruses and their integration sites using next-generation sequencing of human cancer tissue

SUMMARY We developed a new algorithmic method, VirusSeq, for detecting known viruses and their integration sites in the human genome using next-generation sequencing data. We evaluated VirusSeq on whole-transcriptome sequencing (RNA-Seq) data of 256 human cancer samples from The Cancer Genome Atlas. Using these data, we showed that VirusSeq accurately detects the known viruses and their integration sites with high sensitivity and specificity. VirusSeq can also perform this function using whole-genome sequencing data of human tissue. AVAILABILITY VirusSeq has been implemented in PERL and is available at http://odin.mdacc.tmc.edu/∼xsu1/VirusSeq.html. CONTACT xsu1@mdanderson.org SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Comparative Sequence Analysis of Mycobacterium leprae and the New Leprosy-Causing Mycobacterium lepromatosis

Mycobacterium lepromatosis is a newly discovered leprosy-causing organism. Preliminary phylogenetic analysis of its 16S rRNA gene and a few other gene segments revealed significant divergence from Mycobacterium leprae, a well-known cause of leprosy, that justifies the status of M. lepromatosis as a new species. In this study we analyzed the sequences of 20 genes and pseudogenes (22,814 nucleotides). Overall, the level of matching of these sequences with M. leprae sequences was 90.9%, which substantiated the species-level difference; the levels of matching for the 16S rRNA genes and 14 protein-encoding genes were 98.0% and 93.1%, respectively, but the level of matching for five pseudogenes was only 79.1%. Five conserved protein-encoding genes were selected to construct phylogenetic trees and to calculate the numbers of synonymous substitutions (dS values) and nonsynonymous substitutions (dN values) in the two species. Robust phylogenetic trees constructed using concatenated alignment of these genes placed M. lepromatosis and M. leprae in a tight cluster with long terminal branches, implying that the divergence occurred long ago. The dS and dN values were also much higher than those for other closest pairs of mycobacteria. The dS values were 14 to 28% of the dS values for M. leprae and Mycobacterium tuberculosis, a more divergent pair of species. These results thus indicate that M. lepromatosis and M. leprae diverged approximately 10 million years ago. The M. lepromatosis pseudogenes analyzed that were also pseudogenes in M. leprae showed nearly neutral evolution, and their relative ages were similar to those of M. leprae pseudogenes, suggesting that they were pseudogenes before divergence. Taken together, the results described above indicate that M. lepromatosis and M. leprae diverged from a common ancestor after the massive gene inactivation event described previously for M. leprae.

Preleukaemic clonal haemopoiesis and risk of therapy-related myeloid neoplasms: a case-control study

Background Therapy-related myeloid neoplasms (t-MNs) are often fatal secondary malignancies. Risk factors for t-MNs are not well understood. Recent studies suggested that individuals with clonal hematopoiesis have higher risk of developing hematological malignancies. We hypothesized that cancer patients with clonal hematopoiesis have increased risk of developing t-MNs. Methods We conducted a retrospective case-control study to compare the prevalence of clonal hematopoiesis between patients who developed t-MNs (cases) and who did not develop t-MNs (control). For cases, we studied14 patients with various types of cancers who developed t-MNs and whose paired samples of t-MN bone marrow (BM) and peripheral blood (PB) that were previously obtained at the time of primary cancer diagnosis were available. Fifty four patients with lymphoma who received combination chemotherapy and did not develop t-MNs after at least 5 years of follow up were studied as a control. We performed molecular barcode sequencing of 32 genes on the pre-treatment PB samples to detect clonal hematopoiesis. For the t-MN cases, we also performed targeted gene sequencing on t-MN BM samples and investigated clonal evolution from clonal hematopoiesis to t-MNs. To confirm association between clonal hematopoiesis and t-MN development, we also analyzed prevalence of clonal hematopoiesis in a separate cohort of 74 patients with lymphoma. All of these patients were treated under the prospective randomized trial of frontline chemotherapy with cyclophosphamide, doxorubicin, vincristine, and prednisone (CHOP) with or without melatonin and 5 (7%) of them had developed t-MNs. Findings In 14 patients with t-MNs, we detected pre-leukemic mutations in 10 of their prior PB samples (71%). In control, clonal hematopoiesis was detected in 17 patients (31%), and the cumulative incidence of t-MNs at 5 years was significantly higher in patients with clonal hematopoiesis (30% [95% CI: 16% – 51%] vs. 7% [95% CI: 2% – 21%], P = 0.016). In the separate cohort, 5 patients (7%) developed t-MNs and 4 (80%) of them had clonal hematopoiesis. The cumulative incidence of t-MNs at 10 years was significantly higher in patients with clonal hematopoiesis (29% [95% CI: 8%–53%] vs. 0% [95% CI: 0%–0%], P = 0.0009). Multivariate Fine and Gray model showed that having clonal hematopoiesis significantly increased the risk of t-MN development (HR = 13.7, P = 0.013). Interpretation Pre-leukemic clonal hematopoiesis is frequently detected in patients with t-MNs at the time of their primary cancer diagnosis and before patients were exposed to chemotherapy/radiation therapy. Detection of clonal hematopoiesis significantly increased the risk of t-MN development in patients with lymphoma. These data suggest potential approaches of screening clonal hematopoiesis in cancer patients to identify patients at risk of t-MNs and warrants a validation in prospective trial investigating a role of clonal hematopoiesis as a predictive marker for t-MNs.

Streptococcus mitis strains causing severe clinical disease in cancer patients.

The genetically diverse viridans group streptococci (VGS) are increasingly recognized as the cause of a variety of human diseases. We used a recently developed multilocus sequence analysis scheme to define the species of 118 unique VGS strains causing bacteremia in patients with cancer; Streptococcus mitis (68 patients) and S. oralis (22 patients) were the most frequently identified strains. Compared with patients infected with non–S. mitis strains, patients infected with S. mitis strains were more likely to have moderate or severe clinical disease (e.g., VGS shock syndrome). Combined with the sequence data, whole-genome analyses showed that S. mitis strains may more precisely be considered as >2 species. Furthermore, we found that multiple S. mitis strains induced disease in neutropenic mice in a dose-dependent fashion. Our data define the prominent clinical effect of the group of organisms currently classified as S. mitis and lay the groundwork for increased understanding of this understudied pathogen.

Next-generation sequencing of translocation renal cell carcinoma reveals novel RNA splicing partners and frequent mutations of chromatin-remodeling genes.

Purpose: MITF/TFE translocation renal cell carcinoma (TRCC) is a rare subtype of kidney cancer. Its incidence and the genome-wide characterization of its genetic origin have not been fully elucidated. Experimental Design: We performed RNA and exome sequencing on an exploratory set of TRCC (n = 7), and validated our findings using The Cancer Genome Atlas (TCGA) clear-cell RCC (ccRCC) dataset (n = 460). Results: Using the TCGA dataset, we identified seven TRCC (1.5%) cases and determined their genomic profile. We discovered three novel partners of MITF/TFE (LUC7L3, KHSRP, and KHDRBS2) that are involved in RNA splicing. TRCC displayed a unique gene expression signature as compared with other RCC types, and showed activation of MITF, the transforming growth factor β1 and the PI3K complex targets. Genes differentially spliced between TRCC and other RCC types were enriched for MITF and ID2 targets. Exome sequencing of TRCC revealed a distinct mutational spectrum as compared with ccRCC, with frequent mutations in chromatin-remodeling genes (six of eight cases, three of which were from the TCGA). In two cases, we identified mutations in INO80D, an ATP-dependent chromatin-remodeling gene, previously shown to control the amplitude of the S phase. Knockdown of INO80D decreased cell proliferation in a novel cell line bearing LUC7L3–TFE3 translocation. Conclusions: This genome-wide study defines the incidence of TRCC within a ccRCC-directed project and expands the genomic spectrum of TRCC by identifying novel MITF/TFE partners involved in RNA splicing and frequent mutations in chromatin-remodeling genes. Clin Cancer Res; 20(15); 4129–40. ©2014 AACR.

Sequence type 1 group B Streptococcus, an emerging cause of invasive disease in adults, evolves by small genetic changes

Significance Serotype V group B Streptococcus (GBS) infection rates in humans have steadily increased during the past several decades. We determined that 92% of bloodstream infections caused by serotype V GBS in Houston and Toronto are caused by genetically related strains called sequence type (ST) 1. Whole-genome analysis of 202 serotype V ST-1 strains revealed the molecular relationship among these strains and that they are closely related to a bovine strain. Moreover, we found that a subset of GBS genes is under selective evolutionary pressure, indicating that proteins produced by these genes likely contribute to GBS host–pathogen interaction. These data will assist in understanding how bacteria adapt to cause disease in humans, thereby potentially informing new preventive and therapeutic strategies. The molecular mechanisms underlying pathogen emergence in humans is a critical but poorly understood area of microbiologic investigation. Serotype V group B Streptococcus (GBS) was first isolated from humans in 1975, and rates of invasive serotype V GBS disease significantly increased starting in the early 1990s. We found that 210 of 229 serotype V GBS strains (92%) isolated from the bloodstream of nonpregnant adults in the United States and Canada between 1992 and 2013 were multilocus sequence type (ST) 1. Elucidation of the complete genome of a 1992 ST-1 strain revealed that this strain had the highest homology with a GBS strain causing cow mastitis and that the 1992 ST-1 strain differed from serotype V strains isolated in the late 1970s by acquisition of cell surface proteins and antimicrobial resistance determinants. Whole-genome comparison of 202 invasive ST-1 strains detected significant recombination in only eight strains. The remaining 194 strains differed by an average of 97 SNPs. Phylogenetic analysis revealed a temporally dependent mode of genetic diversification consistent with the emergence in the 1990s of ST-1 GBS as major agents of human disease. Thirty-one loci were identified as being under positive selective pressure, and mutations at loci encoding polysaccharide capsule production proteins, regulators of pilus expression, and two-component gene regulatory systems were shown to affect the bacterial phenotype. These data reveal that phenotypic diversity among ST-1 GBS is mainly driven by small genetic changes rather than extensive recombination, thereby extending knowledge into how pathogens adapt to humans.

Dual-Site Phosphorylation of the Control of Virulence Regulator Impacts Group A Streptococcal Global Gene Expression and Pathogenesis

Phosphorylation relays are a major mechanism by which bacteria alter transcription in response to environmental signals, but understanding of the functional consequences of bacterial response regulator phosphorylation is limited. We sought to characterize how phosphorylation of the control of virulence regulator (CovR) protein from the major human pathogen group A Streptococcus (GAS) influences GAS global gene expression and pathogenesis. CovR mainly serves to repress GAS virulence factor-encoding genes and has been shown to homodimerize following phosphorylation on aspartate-53 (D53) in vitro. We discovered that CovR is phosphorylated in vivo and that such phosphorylation is partially heat-stable, suggesting additional phosphorylation at non-aspartate residues. Using mass spectroscopy along with targeted mutagenesis, we identified threonine-65 (T65) as an additional CovR phosphorylation site under control of the serine/threonine kinase (Stk). Phosphorylation on T65, as mimicked by the recombinant CovR T65E variant, abolished in vitro CovR D53 phosphorylation. Similarly, isoallelic GAS strains that were either unable to be phosphorylated at D53 (CovR-D53A) or had functional constitutive phosphorylation at T65 (CovR-T65E) had essentially an identical gene repression profile to each other and to a CovR-inactivated strain. However, the CovR-D53A and CovR-T65E isoallelic strains retained the ability to positively influence gene expression that was abolished in the CovR-inactivated strain. Consistent with these observations, the CovR-D53A and CovR-T65E strains were hypervirulent compared to the CovR-inactivated strain in a mouse model of invasive GAS disease. Surprisingly, an isoalleic strain unable to be phosphorylated at CovR T65 (CovR-T65A) was hypervirulent compared to the wild-type strain, as auto-regulation of covR gene expression resulted in lower covR gene transcript and CovR protein levels in the CovR-T65A strain. Taken together, these data establish that CovR is phosphorylated in vivo and elucidate how the complex interplay between CovR D53 activating phosphorylation, T65 inhibiting phosphorylation, and auto-regulation impacts streptococcal host-pathogen interaction.

Characterization of the Effect of the Histidine Kinase CovS on Response Regulator Phosphorylation in Group A Streptococcus

ABSTRACT Two-component gene regulatory systems (TCSs) are a major mechanism by which bacteria respond to environmental stimuli and thus are critical to infectivity. For example, the control of virulence regulator/sensor kinase (CovRS) TCS is central to the virulence of the major human pathogen group A Streptococcus (GAS). Here, we used a combination of quantitative in vivo phosphorylation assays, isoallelic strains that varied by only a single amino acid in CovS, and transcriptome analyses to characterize the impact of CovS on CovR phosphorylation and GAS global gene expression. We discovered that CovS primarily serves to phosphorylate CovR, thereby resulting in the repression of virulence factor-encoding genes. However, a GAS strain selectively deficient in CovS phosphatase activity had a distinct transcriptome relative to that of its parental strain, indicating that both CovS kinase and phosphatase activities influence the CovR phosphorylation status. Surprisingly, compared to a serotype M3 strain, serotype M1 GAS strains had high levels of phosphorylated CovR, low transcript levels of CovR-repressed genes, and strikingly different responses to environmental cues. Moreover, the inactivation of CovS in the serotype M1 background resulted in a greater decrease in phosphorylated CovR levels and a greater increase in the transcript levels of CovR-repressed genes than did CovS inactivation in a serotype M3 strain. These data clarify the influence of CovS on the CovR phosphorylation status and provide insight into why serotype M1 GAS strains have high rates of spontaneous mutations in covS during invasive GAS infection, thus providing a link between TCS molecular function and the epidemiology of deadly bacterial infections.