Integrative oncogene-dependency mapping identifies RIT1 vulnerabilities and synergies in lung cancer.

CRISPR-based cancer dependency maps are accelerating advances in cancer precision medicine, but adequate functional maps are limited to the most common oncogenes. To identify opportunities for therapeutic intervention in other rarer subsets of cancer, we investigate the oncogene-specific dependencies conferred by the lung cancer oncogene, RIT1. Here, genome-wide CRISPR screening in KRAS, EGFR, and RIT1-mutant isogenic lung cancer cells identifies shared and unique vulnerabilities of each oncogene. Combining this genetic data with small-molecule sensitivity profiling, we identify a unique vulnerability of RIT1-mutant cells to loss of spindle assembly checkpoint regulators. Oncogenic RIT1M90I weakens the spindle assembly checkpoint and perturbs mitotic timing, resulting in sensitivity to Aurora A inhibition. In addition, we observe synergy between mutant RIT1 and activation of YAP1 in multiple models and frequent nuclear overexpression of YAP1 in human primary RIT1-mutant lung tumors. These results provide a genome-wide atlas of oncogenic RIT1 functional interactions and identify components of the RAS pathway, spindle assembly checkpoint, and Hippo/YAP1 network as candidate therapeutic targets in RIT1-mutant lung cancer.

Deep learning using tumor HLA peptide mass spectrometry datasets improves neoantigen identification.

Neoantigens, which are expressed on tumor cells, are one of the main targets of an effective antitumor T-cell response. Cancer immunotherapies to target neoantigens are of growing interest and are in early human trials, but methods to identify neoantigens either require invasive or difficult-to-obtain clinical specimens, require the screening of hundreds to thousands of synthetic peptides or tandem minigenes, or are only relevant to specific human leukocyte antigen (HLA) alleles. We apply deep learning to a large (N = 74 patients) HLA peptide and genomic dataset from various human tumors to create a computational model of antigen presentation for neoantigen prediction. We show that our model, named EDGE, increases the positive predictive value of HLA antigen prediction by up to ninefold. We apply EDGE to enable identification of neoantigens and neoantigen-reactive T cells using routine clinical specimens and small numbers of synthetic peptides for most common HLA alleles. EDGE could enable an improved ability to develop neoantigen-targeted immunotherapies for cancer patients.

Burly1 is a mouse QTL for lean body mass that maps to a 0.8-Mb region of chromosome 2.

To fine map a mouse QTL for lean body mass (Burly1), we used information from intercross, backcross, consomic, and congenic mice derived from the C57BL/6ByJ (host) and 129P3/J (donor) strains. The results from these mapping populations were concordant and showed that Burly1 is located between 151.9 and 152.7 Mb (rs33197365 to rs3700604) on mouse chromosome 2. The congenic region harboring Burly1 contains 26 protein-coding genes, 11 noncoding RNA elements (e.g., lncRNA), and 4 pseudogenes, with 1949 predicted functional variants. Of the protein-coding genes, 7 have missense variants, including genes that may contribute to lean body weight, such as Angpt41, Slc52c3, and Rem1. Lean body mass was increased by the B6-derived variant relative to the 129-derived allele. Burly1 influenced lean body weight at all ages but not food intake or locomotor activity. However, congenic mice with the B6 allele produced more heat per kilogram of lean body weight than did controls, pointing to a genotype effect on lean mass metabolism. These results show the value of integrating information from several mapping populations to refine the map location of body composition QTLs and to identify a short list of candidate genes.

Adiposity QTL Adip20 decomposes into at least four loci when dissected using congenic strains.

An average mouse in midlife weighs between 25 and 30 g, with about a gram of tissue in the largest adipose depot (gonadal), and the weight of this depot differs between inbred strains. Specifically, C57BL/6ByJ mice have heavier gonadal depots on average than do 129P3/J mice. To understand the genetic contributions to this trait, we mapped several quantitative trait loci (QTLs) for gonadal depot weight in an F2 intercross population. Our goal here was to fine-map one of these QTLs, Adip20 (formerly Adip5), on mouse chromosome 9. To that end, we analyzed the weight of the gonadal adipose depot from newly created congenic strains. Results from the sequential comparison method indicated at least four rather than one QTL; two of the QTLs were less than 0.5 Mb apart, with opposing directions of allelic effect. Different types of evidence (missense and regulatory genetic variation, human adiposity/body mass index orthologues, and differential gene expression) implicated numerous candidate genes from the four QTL regions. These results highlight the value of mouse congenic strains and the value of this sequential method to dissect challenging genetic architecture.

Genomic Activation of PPARG Reveals a Candidate Therapeutic Axis in Bladder Cancer.

The PPARG gene encoding the nuclear receptor PPARγ is activated in bladder cancer, either directly by gene amplification or mutation, or indirectly by mutation of the RXRA gene, which encodes the heterodimeric partner of PPARγ. Here, we show that activating alterations of PPARG or RXRA lead to a specific gene expression signature in bladder cancers. Reducing PPARG activity, whether by pharmacologic inhibition or genetic ablation, inhibited proliferation of PPARG-activated bladder cancer cells. Our results offer a preclinical proof of concept for PPARG as a candidate therapeutic target in bladder cancer. Cancer Res; 77(24); 6987-98. ©2017 AACR.

Genetic analysis of impaired trimethylamine metabolism using whole exome sequencing.

BACKGROUND: Trimethylaminuria (TMAU) is a genetic disorder whereby people cannot convert trimethylamine (TMA) to its oxidized form (TMAO), a process that requires the liver enzyme FMO3. Loss-of-function variants in the FMO3 gene are a known cause of TMAU. In addition to the inability to metabolize TMA precursors like choline, patients often emit a characteristic odor because while TMAO is odorless, TMA has a fishy smell. The Monell Chemical Senses Center is a research institute with a program to evaluate people with odor complaints for TMAU. METHODS: Here we evaluated ten subjects by (1) odor evaluation by a trained sensory panel, (2) analysis of their urine concentration of TMA relative to TMAO before and after choline ingestion, and (3) whole exome sequencing as well as subsequent variant analysis of all ten samples to investigate the genetics of TMAU. RESULTS: While all subjects reported they often emitted a fish-like odor, none had this malodor during sensory evaluation. However, all were impaired in their ability to produce >90% TMAO/TMA in their urine and thus met the criteria for TMAU. To probe for genetic causes, the exome of each subject was sequenced, and variants were filtered by genes with a known (FMO3) or expected effect on TMA metabolism function (other oxidoreductases). We filtered the remaining variants by allele frequency and predicated functional effects. We identified one subject that had a rare loss-of-function FMO3 variant and six with more common decreased-function variants. In other oxidoreductases genes, five subjects had four novel rare single-nucleotide polymorphisms as well as one rare insertion/deletion. Novel in this context means no investigators have previously linked these variants to TMAU although they are in dbSNP. CONCLUSIONS: Thus, variants in genes other than FMO3 may cause TMAU and the genetic variants identified here serve as a starting point for future studies of impaired TMA metabolism.

Metagenomic Characterization of Microbial Communities In Situ Within the Deeper Layers of the Ileum in Crohn's Disease.

BACKGROUND & AIMS: Microbial dysbiosis and aberrant host-microbe interactions in the gut are believed to contribute to the development and progression of Crohn's disease (CD). Microbiome studies in CD typically have focused on microbiota in feces or superficial mucosal layers of the colon because accessing DNA from deeper layers of the bowel is challenging. In this study, we analyzed the deep tissue microbiome in patients who underwent surgical resection of the small intestine. METHODS: Paraffin blocks were obtained from 12 CD patients undergoing ileocecal resection, and healthy ileum samples (inflammatory bowel disease-free controls) were obtained from 12 patients undergoing surgery for right-sided colon cancer. Diseased and healthy-appearing ileum was identified using microscopy, and paraffin blocks were macrodissected using a core needle to specifically isolate DNA. Illumina Whole Genome Sequencing was used for microbial sequence identification and subsequent taxonomic classification using the PathSeq tool. RESULTS: We observed significant differences between the microbiome of CD samples vs inflammatory bowel disease-free controls, including depletion of Bacteroidetes and Clostridia. Notably, microbial composition at the phyla level did not differ markedly between healthy and diseased areas of CD patients. However, we observed enrichment of potentially pathogenic organisms at the species level. CONCLUSIONS: Our study showed dysbiosis within deeper layers of the ileum of CD patients, specifically enrichment of enterotoxigenic Staphylococcus aureus and an environmental Mycobacterium species not described previously. Future studies with larger cohort sizes are warranted to confirm these findings. Studies would benefit from effective microbial DNA extraction methods from paraffin sections and host nucleic acid depletion approaches to increase microbial read coverage.

Body Composition QTLs Identified in Intercross Populations Are Reproducible in Consomic Mouse Strains.

Genetic variation contributes to individual differences in obesity, but defining the exact relationships between naturally occurring genotypes and their effects on fatness remains elusive. As a step toward positional cloning of previously identified body composition quantitative trait loci (QTLs) from F2 crosses of mice from the C57BL/6ByJ and 129P3/J inbred strains, we sought to recapture them on a homogenous genetic background of consomic (chromosome substitution) strains. Male and female mice from reciprocal consomic strains originating from the C57BL/6ByJ and 129P3/J strains were bred and measured for body weight, length, and adiposity. Chromosomes 2, 7, and 9 were selected for substitution because previous F2 intercross studies revealed body composition QTLs on these chromosomes. We considered a QTL confirmed if one or both sexes of one or both reciprocal consomic strains differed significantly from the host strain in the expected direction after correction for multiple testing. Using these criteria, we confirmed two of two QTLs for body weight (Bwq5-6), three of three QTLs for body length (Bdln3-5), and three of three QTLs for adiposity (Adip20, Adip26 and Adip27). Overall, this study shows that despite the biological complexity of body size and composition, most QTLs for these traits are preserved when transferred to consomic strains; in addition, studying reciprocal consomic strains of both sexes is useful in assessing the robustness of a particular QTL.

DNA copy number analysis of metastatic urothelial carcinoma with comparison to primary tumors.

BACKGROUND: To date, there have been no reports characterizing the genome-wide somatic DNA chromosomal copy-number alteration landscape in metastatic urothelial carcinoma. We sought to characterize the DNA copy-number profile in a cohort of metastatic samples and compare them to a cohort of primary urothelial carcinoma samples in order to identify changes that are associated with progression from primary to metastatic disease. METHODS: Using molecular inversion probe array analysis we compared genome-wide chromosomal copy-number alterations between 30 metastatic and 29 primary UC samples. Whole transcriptome RNA-Seq analysis was also performed in primary and matched metastatic samples which was available for 9 patients. RESULTS: Based on a focused analysis of 32 genes in which alterations may be clinically actionable, there were significantly more amplifications/deletions in metastases (8.6% vs 4.5%, p < 0.001). In particular, there was a higher frequency of E2F3 amplification in metastases (30% vs 7%, p = 0.046). Paired primary and metastatic tissue was available for 11 patients and 3 of these had amplifications of potential clinical relevance in metastases that were not in the primary tumor including ERBB2, CDK4, CCND1, E2F3, and AKT1. The transcriptional activity of these amplifications was supported by RNA expression data. CONCLUSIONS: The discordance in alterations between primary and metastatic tissue may be of clinical relevance in the era of genomically directed precision cancer medicine.

The tumor virus landscape of AIDS-related lymphomas.

Immunodeficiency dramatically increases susceptibility to cancer as a result of reduced immune surveillance and enhanced opportunities for virus-mediated oncogenesis. Although AIDS-related lymphomas (ARLs) are frequently associated with known oncogenic viruses, many cases contain no known transforming virus. To discover novel transforming viruses, we profiled a set of ARL samples using whole transcriptome sequencing. We determined that Epstein-Barr virus (EBV) was the only virus detected in the tumor samples of this cohort, suggesting that if unidentified pathogens exist in this disease, they are present in <10% of cases or undetectable by our methods. To evaluate the role of EBV in ARL pathogenesis, we analyzed viral gene expression and found highly heterogeneous patterns of viral transcription across samples. We also found significant heterogeneity of viral antigen expression across a large cohort, with many patient samples presenting with restricted type I viral latency, indicating that EBV latency proteins are under increased immunosurveillance in the post-combined antiretroviral therapies era. Furthermore, EBV infection of lymphoma cells in HIV-positive individuals was associated with a distinct host gene expression program. These findings provide insight into the joint host-virus regulatory network of primary ARL tumor samples and expand our understanding of virus-associated oncogenesis. Our findings may also have therapeutic implications, as treatment may be personalized to target specific viral and virus-associated host processes that are only present in a subset of patients.

Malawi polyomavirus is a prevalent human virus that interacts with known tumor suppressors.

Malawi polyomavirus (MWPyV) is a recently identified human polyomavirus. Serology for MWPyV VP1 indicates that infection frequently occurs in childhood and reaches a prevalence of 75% in adults. The MWPyV small T antigen (ST) binds protein phosphatase 2A (PP2A), and the large T antigen (LT) binds pRb, p107, p130, and p53. However, the MWPyV LT was less stable than the simian virus 40 (SV40) LT and was unable to promote the growth of normal cells. This report confirms that MWPyV is a widespread human virus expressing T antigens with low transforming potential.

In search of a candidate pathogen for giant cell arteritis: sequencing-based characterization of the giant cell arteritis microbiome.

OBJECTIVE: To characterize the microbiome of the temporal artery in patients with giant cell arteritis (GCA), and to apply an unbiased and comprehensive shotgun sequencing-based approach to determine whether there is an enrichment of candidate pathogens in the affected tissue. METHODS: Temporal artery biopsy specimens were collected from patients at a single institution over a period of 4 years, and unbiased DNA sequencing was performed on 17 formalin-fixed, paraffin-embedded specimens. Twelve of the 17 patients fulfilled the clinical and histopathologic criteria for GCA, and the other 5 patients served as controls. Using PathSeq software, human DNA sequences were computationally subtracted, and the remaining non-human DNA sequences were taxonomically classified using a comprehensive microbial sequence database. The relative abundance of microbes was inferred based on read counts assigned to each organism. Comparison of the microbial diversity between GCA cases and controls was carried out using hierarchical clustering and linear discriminant analysis of effect size. RESULTS: Propionibacterium acnes and Escherichia coli were the most abundant microorganisms in 16 of the 17 samples, and Moraxella catarrhalis was the most abundant organism in 1 control sample. Pathogens previously described to be correlated with GCA were not differentially abundant in cases compared to controls. There was not a significant burden of likely pathogenic viruses. CONCLUSION: DNA sequencing of temporal artery biopsy specimens from GCA cases, in comparison with non-GCA controls, showed no evidence of previously identified candidate GCA pathogens. A single pathogen was not clearly and consistently associated with GCA in this case series.

Landscape of genomic alterations in cervical carcinomas.

Cervical cancer is responsible for 10-15% of cancer-related deaths in women worldwide. The aetiological role of infection with high-risk human papilloma viruses (HPVs) in cervical carcinomas is well established. Previous studies have also implicated somatic mutations in PIK3CA, PTEN, TP53, STK11 and KRAS as well as several copy-number alterations in the pathogenesis of cervical carcinomas. Here we report whole-exome sequencing analysis of 115 cervical carcinoma-normal paired samples, transcriptome sequencing of 79 cases and whole-genome sequencing of 14 tumour-normal pairs. Previously unknown somatic mutations in 79 primary squamous cell carcinomas include recurrent E322K substitutions in the MAPK1 gene (8%), inactivating mutations in the HLA-B gene (9%), and mutations in EP300 (16%), FBXW7 (15%), NFE2L2 (4%), TP53 (5%) and ERBB2 (6%). We also observe somatic ELF3 (13%) and CBFB (8%) mutations in 24 adenocarcinomas. Squamous cell carcinomas have higher frequencies of somatic nucleotide substitutions occurring at cytosines preceded by thymines (Tp*C sites) than adenocarcinomas. Gene expression levels at HPV integration sites were statistically significantly higher in tumours with HPV integration compared with expression of the same genes in tumours without viral integration at the same site. These data demonstrate several recurrent genomic alterations in cervical carcinomas that suggest new strategies to combat this disease.

Sequence-based discovery of Bradyrhizobium enterica in cord colitis syndrome.

BACKGROUND: Immunosuppression is associated with a variety of idiopathic clinical syndromes that may have infectious causes. It has been hypothesized that the cord colitis syndrome, a complication of umbilical-cord hematopoietic stem-cell transplantation, is infectious in origin. METHODS: We performed shotgun DNA sequencing on four archived, paraffin-embedded endoscopic colon-biopsy specimens obtained from two patients with cord colitis. Computational subtraction of human and known microbial sequences and assembly of residual sequences into a bacterial draft genome were performed. We used polymerase-chain-reaction (PCR) assays and fluorescence in situ hybridization to determine whether the corresponding bacterium was present in additional patients and controls. RESULTS: DNA sequencing of the biopsy specimens revealed more than 2.5 million sequencing reads that did not match known organisms. These sequences were computationally assembled into a 7.65-Mb draft genome showing a high degree of homology with genomes of bacteria in the bradyrhizobium genus. The corresponding newly discovered bacterium was provisionally named Bradyrhizobium enterica. PCR identified B. enterica nucleotide sequences in biopsy specimens from all three additional patients with cord colitis whose samples were tested, whereas B. enterica sequences were absent in samples obtained from healthy controls and patients with colon cancer or graft-versus-host disease. CONCLUSIONS: We assembled a novel bacterial draft genome from the direct sequencing of tissue specimens from patients with cord colitis. Association of these sequences with cord colitis suggests that B. enterica may be an opportunistic human pathogen. (Funded by the National Cancer Institute and others.)

Genetic analysis of chemosensory traits in human twins.

We explored genetic influences on the perception of taste and smell stimuli. Adult twins rated the chemosensory aspects of water, sucrose, sodium chloride, citric acid, ethanol, quinine hydrochloride, phenylthiocarbamide (PTC), potassium chloride, calcium chloride, cinnamon, androstenone, Galaxolide™, cilantro, and basil. For most traits, individual differences were stable over time and some traits were heritable (h(2) from 0.41 to 0.71). Subjects were genotyped for 44 single nucleotide polymorphisms within and near genes related to taste and smell. The results of these association analyses confirmed previous genotype-phenotype results for PTC, quinine, and androstenone. New associations were detected for ratings of basil and a bitter taste receptor gene, TAS2R60, and between cilantro and variants in three genes (TRPA1, GNAT3, and TAS2R50). The flavor of ethanol was related to variation within an olfactory receptor gene (OR7D4) and a gene encoding a subunit of the epithelial sodium channel (SCNN1D). Our study demonstrates that person-to-person differences in the taste and smell perception of simple foods and drinks are partially accounted for by genetic variation within chemosensory pathways.

Sequence analysis of mutations and translocations across breast cancer subtypes.

Breast carcinoma is the leading cause of cancer-related mortality in women worldwide, with an estimated 1.38 million new cases and 458,000 deaths in 2008 alone. This malignancy represents a heterogeneous group of tumours with characteristic molecular features, prognosis and responses to available therapy. Recurrent somatic alterations in breast cancer have been described, including mutations and copy number alterations, notably ERBB2 amplifications, the first successful therapy target defined by a genomic aberration. Previous DNA sequencing studies of breast cancer genomes have revealed additional candidate mutations and gene rearrangements. Here we report the whole-exome sequences of DNA from 103 human breast cancers of diverse subtypes from patients in Mexico and Vietnam compared to matched-normal DNA, together with whole-genome sequences of 22 breast cancer/normal pairs. Beyond confirming recurrent somatic mutations in PIK3CA, TP53, AKT1, GATA3 and MAP3K1, we discovered recurrent mutations in the CBFB transcription factor gene and deletions of its partner RUNX1. Furthermore, we have identified a recurrent MAGI3-AKT3 fusion enriched in triple-negative breast cancer lacking oestrogen and progesterone receptors and ERBB2 expression. The MAGI3-AKT3 fusion leads to constitutive activation of AKT kinase, which is abolished by treatment with an ATP-competitive AKT small-molecule inhibitor.

Genomic analysis identifies association of Fusobacterium with colorectal carcinoma.

The tumor microenvironment of colorectal carcinoma is a complex community of genomically altered cancer cells, nonneoplastic cells, and a diverse collection of microorganisms. Each of these components may contribute to carcinogenesis; however, the role of the microbiota is the least well understood. We have characterized the composition of the microbiota in colorectal carcinoma using whole genome sequences from nine tumor/normal pairs. Fusobacterium sequences were enriched in carcinomas, confirmed by quantitative PCR and 16S rDNA sequence analysis of 95 carcinoma/normal DNA pairs, while the Bacteroidetes and Firmicutes phyla were depleted in tumors. Fusobacteria were also visualized within colorectal tumors using FISH. These findings reveal alterations in the colorectal cancer microbiota; however, the precise role of Fusobacteria in colorectal carcinoma pathogenesis requires further investigation.

Body fat distribution and organ weights of 14 common strains and a 22-strain consomic panel of rats.

The goal of this study was to determine the adiposity of a range of rat strains, including a panel of consomics, to estimate heritability. To that end, we assessed the body fat distribution and organ weights of groups of adult male rats from 3 outbred strains, 11 inbred strains and 22 consomic strains. We measured the weights of the gonadal, retroperitoneal, mesenteric, femoral, subscapular and pericardial white fat depots, the subscapular brown fat depot, the kidneys, liver, heart, spleen, and brain. Strains were compared for the measured weight of each of these adipose depots and organs, and also for these weights adjusted statistically for body size. All individual adipose depot and organ weights were highly heritable, in most cases h(2)>0.50. The fourteen inbred and outbred rat strains were not very different in body length but there was a three-fold difference in body weight, and up to a twenty-fold difference in the weight of some adipose depots. Comparison of the FHH-Chr n(BN) consomic strains with the FHH host strain revealed 98 quantitative trait loci (QTLs) for body composition and organ weight, with the introgressed chromosome reducing weight or adiposity in most cases. These results can be used to guide the choice of appropriate rat strains for future studies of the genetic architecture of obesity and body size.

Psychophysical dissection of genotype effects on human bitter perception.

The purpose of this study was to define the effects of individual polymorphisms within the haplotypes of the TAS2R38 taste receptor gene on human bitter taste perception. A racially and ethnically diverse sample of children and adults (N = 980) was phenotyped for thresholds of 6-n-propylthiouracil (PROP) and genotyped for 3 polymorphisms of the TAS2R38 gene (A49P, V262A, I296V). Subjects were grouped according to their diplotype (i.e., specific combinations of haplotypes) and compared for PROP thresholds. By contrasting subjects with particular diplotypes, we found that in addition to A49P, V262A and I296V were related to the ability of the subjects to detect PROP. The V262A variant site affected the ability of subjects to detect mid-range concentrations of PROP, whereas the I296V variant site affected the ability of subjects to perceive PROP at the lowest concentration. These data agree with results from previous studies using cell-based assays for 2 variant sites (A49P and V262A) but not those for the I296V variant site. The reason for the discordant results is not known but it highlights the need for psychophysical as well as cell-based methods to understand the genotype-phenotype relationship for taste receptors. Human PROP sensitivity is determined by the combination of each of these 3 polymorphisms within the TAS2R38 gene.