Immune Activation following Irbesartan Treatment in a Colorectal Cancer Patient: A Case Study.

Colorectal cancers are one of the most prevalent tumour types worldwide and, despite the emergence of targeted and biologic therapies, have among the highest mortality rates. The Personalized OncoGenomics (POG) program at BC Cancer performs whole genome and transcriptome analysis (WGTA) to identify specific alterations in an individual's cancer that may be most effectively targeted. Informed using WGTA, a patient with advanced mismatch repair-deficient colorectal cancer was treated with the antihypertensive drug irbesartan and experienced a profound and durable response. We describe the subsequent relapse of this patient and potential mechanisms of response using WGTA and multiplex immunohistochemistry (m-IHC) profiling of biopsies before and after treatment from the same metastatic site of the L3 spine. We did not observe marked differences in the genomic landscape before and after treatment. Analyses revealed an increase in immune signalling and infiltrating immune cells, particularly CD8+ T cells, in the relapsed tumour. These results indicate that the observed anti-tumour response to irbesartan may have been due to an activated immune response. Determining whether there may be other cancer contexts in which irbesartan may be similarly valuable will require additional studies.

Whole-genome and transcriptome analysis enhances precision cancer treatment options.

BACKGROUND: Recent advances are enabling delivery of precision genomic medicine to cancer clinics. While the majority of approaches profile panels of selected genes or hotspot regions, comprehensive data provided by whole-genome and transcriptome sequencing and analysis (WGTA) present an opportunity to align a much larger proportion of patients to therapies. PATIENTS AND METHODS: Samples from 570 patients with advanced or metastatic cancer of diverse types enrolled in the Personalized OncoGenomics (POG) program underwent WGTA. DNA-based data, including mutations, copy number and mutation signatures, were combined with RNA-based data, including gene expression and fusions, to generate comprehensive WGTA profiles. A multidisciplinary molecular tumour board used WGTA profiles to identify and prioritize clinically actionable alterations and inform therapy. Patient responses to WGTA-informed therapies were collected. RESULTS: Clinically actionable targets were identified for 83% of patients, of which 37% of patients received WGTA-informed treatments. RNA expression data were particularly informative, contributing to 67% of WGTA-informed treatments; 25% of treatments were informed by RNA expression alone. Of a total 248 WGTA-informed treatments, 46% resulted in clinical benefit. RNA expression data were comparable to DNA-based mutation and copy number data in aligning to clinically beneficial treatments. Genome signatures also guided therapeutics including platinum, poly-ADP ribose polymerase inhibitors and immunotherapies. Patients accessed WGTA-informed treatments through clinical trials (19%), off-label use (35%) and as standard therapies (46%) including those which would not otherwise have been the next choice of therapy, demonstrating the utility of genomic information to direct use of chemotherapies as well as targeted therapies. CONCLUSIONS: Integrating RNA expression and genome data illuminated treatment options that resulted in 46% of treated patients experiencing positive clinical benefit, supporting the use of comprehensive WGTA profiling in clinical cancer care.

Successful targeting of the NRG1 pathway indicates novel treatment strategy for metastatic cancer.

BACKGROUND: NRG1 fusion-positive lung cancers have emerged as potentially actionable events in lung cancer, but clinical support is currently limited and no evidence of efficacy of this approach in cancers beyond lung has been shown. PATIENTS AND METHODS: Here, we describe two patients with advanced cancers refractory to standard therapies. Patient 1 had lung adenocarcinoma and patient 2 cholangiocarcinoma. Whole-genome and transcriptome sequencing were carried out for these cases with select findings validated by fluorescence in situ hybridization. RESULTS: Both tumors were found to be positive for NRG1 gene fusions. In patient 1, an SDC4-NRG1 gene fusion was detected, similar gene fusions having been described in lung cancers previously. In patient 2, a novel ATP1B1-NRG1 gene fusion was detected. Cholangiocarcinoma is not a disease type in which NRG1 fusions had been described previously. Integrative genome analysis was used to assess the potential functional significance of the detected genomic events including the gene fusions, prioritizing therapeutic strategies targeting the HER-family of growth factor receptors. Both patients were treated with the pan HER-family kinase inhibitor afatinib and both displayed significant and durable response to treatment. Upon progression sites of disease were sequenced. The lack of obvious genomic events to describe the disease progression indicated that broad transcriptomic or epigenetic mechanisms could be attributed to the lack of prolonged response to afatinib. CONCLUSION: These observations lend further support to the use of pan HER-tyrosine kinase inhibitors for the treatment of NRG1 fusion-positive in both cancers of lung and hepatocellular origin and indicate more broadly that cancers found to be NRG1 fusion-positive may benefit from such a clinical approach regardless of their site of origin. CLINICAL TRIAL INFORMATION: Personalized Oncogenomics (POG) Program of British Columbia: Utilization of Genomic Analysis to Better Understand Tumour Heterogeneity and Evolution (NCT02155621).

Response to angiotensin blockade with irbesartan in a patient with metastatic colorectal cancer.

BACKGROUND: A patient suffering from metastatic colorectal cancer, treatment-related toxicity and resistance to standard chemotherapy and radiation was assessed as part of a personalized oncogenomics initiative to derive potential alternative therapeutic strategies. PATIENTS AND METHODS: Whole-genome and transcriptome sequencing was used to interrogate a metastatic tumor refractory to standard treatments of a patient with mismatch repair-deficient metastatic colorectal cancer. RESULTS: Integrative genomic analysis indicated overexpression of the AP-1 transcriptional complex suggesting experimental therapeutic rationales, including blockade of the renin-angiotensin system. This led to the repurposing of the angiotensin II receptor antagonist, irbesartan, as an anticancer therapy, resulting in the patient experiencing a dramatic and durable response. CONCLUSIONS: This case highlights the utility of comprehensive integrative genomic profiling and bioinformatics analysis to provide hypothetical rationales for personalized treatment options.

RPS6KA2, a putative tumour suppressor gene at 6q27 in sporadic epithelial ovarian cancer.

We had previously defined by allele loss studies a minimal region at 6q27 (between D6S264 and D6S297) to contain a putative tumour suppressor gene. The p90 ribosomal S6 kinase-3 gene (p90 Rsk-3, RPS6KA2) maps in this interval. It is a serine-threonine kinase that signals downstream of the mitogen-activated protein kinase pathway. It is expressed in normal ovarian epithelium, whereas reduced or absent in tumours or cell lines. We show that RPS6KA2 is monoallelically expressed in the ovary suggesting that loss of a single expressed allele is sufficient to cause complete loss of expression in cancer cells. Further, we have identified two new isoforms of RPS6KA2 with an alternative start codon. Homozygous deletions were identified within the RPS6KA2 gene in two cell lines. Re-expression of RPS6KA2 in ovarian cancer cell lines suppressed colony formation. In UCI101 cells, the expression of RPS6KA2 reduced proliferation, caused G1 arrest, increased apoptosis, reduced levels of phosphorylated extracellular signal-regulated kinase and altered other cell cycle proteins. In contrast, small interfering RNA against RPS6KA2 showed the opposite effect in 41M cells. The above results suggest that RPS6KA2 is a putative tumour suppressor gene to explain allele loss at 6q27.

Small regions of overlapping deletions on 6q26 in human astrocytic tumours identified using chromosome 6 tile path array-CGH.

Deletions of chromosome 6 are a common abnormality in diverse human malignancies including astrocytic tumours, suggesting the presence of tumour suppressor genes (TSG). In order to help identify candidate TSGs, we have constructed a chromosome 6 tile path microarray. The array contains 1,780 clones (778 P1-derived artificial chromosome and 1,002 bacterial artificial chromosome) that cover 98.3% of the published chromosome 6 sequences. A total of 104 adult astrocytic tumours (10 diffuse astrocytomas, 30 anaplastic astrocytomas (AA), 64 glioblastomas (GB)) were analysed using this array. Single copy number change was successfully detected and the result was in general concordant with a microsatellite analysis. The pattern of copy number change was complex with multiple interstitial deletions/gains. However, a predominance of telomeric 6q deletions was seen. Two small common and overlapping regions of deletion at 6q26 were identified. One was 1,002 kb in size and contained PACRG and QKI, while the second was 199 kb and harbours a single gene, ARID1B. The data show that the chromosome 6 tile path array is useful in mapping copy number changes with high resolution and accuracy. We confirmed the high frequency of chromosome 6 deletions in AA and GB, and identified two novel commonly deleted regions that may harbour TSGs.

The DNA sequence and analysis of human chromosome 6.

Chromosome 6 is a metacentric chromosome that constitutes about 6% of the human genome. The finished sequence comprises 166,880,988 base pairs, representing the largest chromosome sequenced so far. The entire sequence has been subjected to high-quality manual annotation, resulting in the evidence-supported identification of 1,557 genes and 633 pseudogenes. Here we report that at least 96% of the protein-coding genes have been identified, as assessed by multi-species comparative sequence analysis, and provide evidence for the presence of further, otherwise unsupported exons/genes. Among these are genes directly implicated in cancer, schizophrenia, autoimmunity and many other diseases. Chromosome 6 harbours the largest transfer RNA gene cluster in the genome; we show that this cluster co-localizes with a region of high transcriptional activity. Within the essential immune loci of the major histocompatibility complex, we find HLA-B to be the most polymorphic gene on chromosome 6 and in the human genome.

beta-1,3-Glucuronyltransferase-1 gene implicated as a candidate for a schizophrenia-like psychosis through molecular analysis of a balanced translocation.

We have mapped and sequenced both chromosome breakpoints of a balanced t(6;11)(q14.2;q25) chromosome translocation that segregates with a schizophrenia-like psychosis. Bioinformatics analysis of the regions revealed a number of confirmed and predicted transcripts. No confirmed transcripts are disrupted by either breakpoint. The chromosome 6 breakpoint region is gene poor, the closest transcript being the serotonin receptor 1E (HTR1E) at 625 kb telomeric to the breakpoint. The chromosome 11 breakpoint is situated close to the telomere. The closest gene, beta-1,3-glucuronyltransferase (B3GAT1 or GlcAT-P), is 299 kb centromeric to the breakpoint. B3GAT1 is the key enzyme during the biosynthesis of the carbohydrate epitope HNK-1, which is present on a number of cell adhesion molecules important in neurodevelopment. Mice deleted for the B3GAT1 gene show defects in hippocampal long-term potentiation and in spatial memory formation. We propose that the translocation causes a positional effect on B3GAT1, affecting expression levels and making it a plausible candidate for the psychosis found in this family. More generally, regions close to telomeres are highly polymorphic in both sequence and length in the general population and several studies have implicated subtelomeric deletions as a common cause of idiopathic mental retardation. This leads us to the hypothesis that polymorphic or other variation of the 11q telomere may affect the activity of B3GAT1 and be a risk factor for schizophrenia and related psychoses in the general population.

Cyclin D3 is a target gene of t(6;14)(p21.1;q32.3) of mature B-cell malignancies.

Chromosomal translocation t(6;14)(p21.1;q32.3) has been reported as a rare but recurrent event not only in myeloma and plasma cell leukemia but also in diffuse large B-cell non-Hodgkin lymphoma (B-NHL) (diffuse large B-cell lymphoma [DLBCL]) and splenic lymphoma with villous lymphocytes (SLVL); however, the nature of the target gene(s) has not been determined. This study identified t(6;14)(p21.1;q32.3) in 3 cases of transformed extranodal marginal zone B-NHL, in 1 case of SLVL, and in 1 case of a low-grade B-cell lymphoproliferative disorder. In a sixth case, a CD5(+) DLBCL, the translocation was identified by molecular cloning in the absence of cytogenetically detectable change. Two chromosomal translocation breakpoints were cloned by using long-distance inverse polymerase chain reaction methods. Comparison with the genomic sequence for chromosome 6p21.1 showed breakpoints approximately 59 and 73.5 kilobases 5' of the cyclin D3 (CCND3) gene with no other identifiable transcribed sequences in the intervening region. Although Southern blotting with derived genomic 6p21.1 probes failed to detect other rearrangements, fluorescent in situ hybridization assays, using BAC (bacterial artificial chromosome) clones spanning and flanking the CCND3 locus, along with probes for IGH confirmed localization of 6p21.1 breakpoints within the same region, as well as fusion of the CCND3 and IGH loci. Furthermore, in all cases, high-level expression of CCND3 was demonstrated at RNA and/or protein levels by Northern and Western blotting and by immunohistochemistry. These data implicate CCND3 as a dominant oncogene in the pathogenesis and transformation in several histologic subtypes of mature B-cell malignancies with t(6;14)(p21.1;q32.3) and suggest that CCND3 overexpression seen in about 10% of DLBCL cases may have a genetic basis.

The genomic structure and promoter region of the human parkin gene.

Parkin has been identified as a causative gene of the autosomal recessive juvenile parkinsonism (AR-JP). In this study, we determined the genomic structure of the Parkin gene and identified a core promoter region based on the DNA sequence of 1.4 Mb. The 5'-flanking region contained no apparent TATA or CAAT box elements but several putative cis-elements for various transcription factors. The GC- and CpG-rich regions were observed not only in the 5'-flanking sequence but also in the 5'-part of the first intron of Parkin. We identified an exact starting point of Parkin transcription. A core promoter region was determined by transfecting a series of deletion constructs with a dual luciferase reporter system into human neuroblastoma cells. Furthermore, we located a neighboring novel gene in a head-to-head direction with Parkin with only a 198-bp interval.

Vanin genes are clustered (human 6q22-24 and mouse 10A2B1) and encode isoforms of pantetheinase ectoenzymes.

The mouse Vanin-1 molecule plays a role in thymic reconstitution following damage by irradiation. We recently demonstrated that it is a membrane pantetheinase (EC 3.56.1.-). This molecule is the prototypic member of a larger Vanin family encoded by at least two mouse (Vanin-1 and Vanin-3) and three human (VNN1, VNN2, VNN3) orthologous genes. We now report (1) the structural characterization of the human and mouse Vanin genes and their organization in clusters on the 6q22-24 and 10A2B1 chromosomes, respectively; (2) identification of the human VNN3 gene and the demonstration that the mouse Vanin-3 molecule is secreted by cells, and (3) that the Vanin genes encode different isoforms of the mammalian pantetheinase activity. Thus, the Vanin family represents a novel class of secreted or membrane-associated ectoenzymes. We discuss here their possible role in processes pertaining to tissue repair in the context of oxidative stress.

Characterization of clustered MHC-linked olfactory receptor genes in human and mouse.

Olfactory receptor (OR) loci frequently cluster and are present on most human chromosomes. They are members of the seven transmembrane receptor (7-TM) superfamily and, as such, are part of one of the largest mammalian multigene families, with an estimated copy number of up to 1000 ORs per haploid genome. As their name implies, ORs are known to be involved in the perception of odors and possibly also in other, nonolfaction-related, functions. Here, we report the characterization of ORs that are part of the MHC-linked OR clusters in human and mouse (partial sequence only). These clusters are of particular interest because of their possible involvement in olfaction-driven mate selection. In total, we describe 50 novel OR loci (36 human, 14 murine), making the human MHC-linked cluster the largest sequenced OR cluster in any organism so far. Comparative and phylogenetic analyses confirm the cluster to be MHC-linked but divergent in both species and allow the identification of at least one ortholog that will be useful for future regulatory and functional studies. Quantitative feature analysis shows clear evidence of duplications of blocks of OR genes and reveals the entire cluster to have a genomic environment that is very different from its neighboring regions. Based on in silico transcript analysis, we also present evidence of extensive long-distance splicing in the 5'-untranslated regions and, for the first time, of alternative splicing within the single coding exon of ORs. Taken together with our previous finding that ORs are also polymorphic, the presented data indicate that the expression, function, and evolution of these interesting genes might be more complex than previously thought.

The physical maps for sequencing human chromosomes 1, 6, 9, 10, 13, 20 and X.

We constructed maps for eight chromosomes (1, 6, 9, 10, 13, 20, X and (previously) 22), representing one-third of the genome, by building landmark maps, isolating bacterial clones and assembling contigs. By this approach, we could establish the long-range organization of the maps early in the project, and all contig extension, gap closure and problem-solving was simplified by containment within local regions. The maps currently represent more than 94% of the euchromatic (gene-containing) regions of these chromosomes in 176 contigs, and contain 96% of the chromosome-specific markers in the human gene map. By measuring the remaining gaps, we can assess chromosome length and coverage in sequenced clones.

Comparison of human genetic and sequence-based physical maps.

Recombination is the exchange of information between two homologous chromosomes during meiosis. The rate of recombination per nucleotide, which profoundly affects the evolution of chromosomal segments, is calculated by comparing genetic and physical maps. Human physical maps have been constructed using cytogenetics, overlapping DNA clones and radiation hybrids; but the ultimate and by far the most accurate physical map is the actual nucleotide sequence. The completion of the draft human genomic sequence provides us with the best opportunity yet to compare the genetic and physical maps. Here we describe our estimates of female, male and sex-average recombination rates for about 60% of the genome. Recombination rates varied greatly along each chromosome, from 0 to at least 9 centiMorgans per megabase (cM Mb(-1)). Among several sequence and marker parameters tested, only relative marker position along the metacentric chromosomes in males correlated strongly with recombination rate. We identified several chromosomal regions up to 6 Mb in length with particularly low (deserts) or high (jungles) recombination rates. Linkage disequilibrium was much more common and extended for greater distances in the deserts than in the jungles.

Cloning and characterization of two overlapping genes in a subregion at 6q21 involved in replicative senescence and schizophrenia.

Two new genes were cloned from region 6q21 and characterized. One gene, C6orf4-6, expresses three mRNA isoforms diverging at the 5' and 3' ends, and encodes two protein isoforms that differ by nine amino acids at their amino terminus. The second gene, C6UAS, is transcribed in the antisense orientation from the complementary strand of C6orf4-6. C6UAS overlaps the second exon of C6orf4, where the start codon of protein isoform 1 is located. C6UAS has no apparent ORF and most likely represents a structural RNA gene that is transcribed but not translated. This feature and the antisense polarity of transcription suggest that C6UAS could play a regulatory role on the expression of C6orf4, as indicated by a significant decrease of endogenous C6orf4 expression after transfection of C6UAS cDNA in human fibroblasts. Neither C6UAS nor C6orf4-6 genes show any homology with known human genes. The two genes were cloned from a subregion at 6q21 containing a replicative senescence gene, a tumor suppressor gene and a gene involved in hereditary schizophrenia. In addition, the common fragile site FRA6F was mapped in the same region. Cloning and characterization of C6orf4-6 and C6UAS may help to clarify the structure and the functional role of this important region.

An imprinted locus associated with transient neonatal diabetes mellitus.

Recently, we reported the localization of a gene for transient neonatal diabetes mellitus (TNDM), a rare form of childhood diabetes, to an approximately 5.4 Mb region of chromosome 6q24. We have also shown that TNDM is associated with both paternal uniparental disomy (UPD) of chromosome 6 and paternal duplications of the critical region. The sequencing of P1-derived artificial chromosome clones from within the region of interest has allowed us to further localize the gene and to investigate the methylation status of the region. The gene is now known to reside in a 300-400 kb region of 6q24 which contains several CpG islands. At one island we have demonstrated differential DNA methylation between patients with paternal UPD of chromosome 6 and normal controls. In addition, two patients with TNDM, in whom neither paternal UPD of chromosome 6 nor duplication of 6q24 have been found, show a DNA methylation pattern identical to that of patients with paternal UPD of chromosome 6. Control individuals show a hemizygous methylation pattern. These results show that TNDM can be associated with a methylation change and identify a novel methylation imprint on chromosome 6 associated with TNDM.

Localisation of a gene for transient neonatal diabetes mellitus to an 18.72 cR3000 (approximately 5.4 Mb) interval on chromosome 6q.

Transient neonatal diabetes mellitus (TNDM) is a rare condition which presents with intrauterine growth retardation, dehydration, and failure to thrive. The condition spontaneously resolves before 1 year of age but predisposes patients to type 2 diabetes later in life. We have previously shown that, in some cases, TNDM is associated with paternal uniparental disomy (UPD) of chromosome 6 and suggested that an imprinted gene responsible for TNDM lies within a region of chromosome 6q. By analysing three families, two with duplications (family A and patient C) and one with several affected subjects with normal karyotypes (family B), we have further defined the TNDM critical region. In patient A, polymorphic microsatellite repeat analysis identified a duplicated region of chromosome 6, flanked by markers D6S472 and D6S311. This region was identified on the Sanger Centre's chromosome 6 radiation hybrid map (http://www.sanger.ac.uk/HGP/Chr6) and spanned approximately 60 cR3000. Using markers within the region, 418 unique P1 derived artificial chromosomes (PACs) have been isolated and used to localise the distal breakpoints of the two duplications. Linkage analysis of the familial case with a normal karyotype identified a recombination within the critical region. This recombination has been identified on the radiation hybrid map and defines the proximal end of the region of interest. We therefore propose that an imprinted gene for TNDM lies within an 18.72 cR3000 (approximately 5.4 Mb) interval on chromosome 6q24.1-q24.3 between markers D6S1699 and D6S1010.

Mutations in a gene encoding a novel protein tyrosine phosphatase cause progressive myoclonus epilepsy.

Lafora's disease (LD; OMIM 254780) is an autosomal recessive form of progressive myoclonus epilepsy characterized by seizures and cumulative neurological deterioration. Onset occurs during late childhood and usually results in death within ten years of the first symptoms. With few exceptions, patients follow a homogeneous clinical course despite the existence of genetic heterogeneity. Biopsy of various tissues, including brain, revealed characteristic polyglucosan inclusions called Lafora bodies, which suggested LD might be a generalized storage disease. Using a positional cloning approach, we have identified at chromosome 6q24 a novel gene, EPM2A, that encodes a protein with consensus amino acid sequence indicative of a protein tyrosine phosphatase (PTP). mRNA transcripts representing alternatively spliced forms of EPM2A were found in every tissue examined, including brain. Six distinct DNA sequence variations in EPM2A in nine families, and one homozygous microdeletion in another family, have been found to cosegregate with LD. These mutations are predicted to cause deleterious effects in the putative protein product, named laforin, resulting in LD.