UROSCAN and UROSCANSEQ: a large-scale multicenter effort towards translation of molecular bladder cancer subtypes into clinical practice - from biobank to RNA-sequencing in real time.

BACKGROUND: Bladder cancer is molecularly one of the most heterogenous malignancies characterized by equally heterogenous clinical outcomes. Standard morphological assessment with pathology and added immunohistochemical analyses is unable to fully address the heterogeneity, but up to now treatment decisions have been made based on such information only. Bladder cancer molecular subtypes will likely provide means for a more personalized bladder cancer care. METHODS: To facilitate further development of bladder cancer molecular subtypes and clinical translation, the UROSCAN-biobank was initiated in 2013 to achieve systematic biobanking of preoperative blood and fresh frozen tumor tissue in a population-based setting. In a second phase, we established in 2018 a parallel logistic pipeline for molecular profiling by RNA-sequencing, to develop and validate clinical implementation of molecular subtyping and actionable molecular target identification in real-time. RESULTS: Until June 2021, 1825 individuals were included in the UROSCAN-biobank, of which 1650 (90%) had primary bladder cancer, 127 (7%) recurrent tumors, and 48 (3%) unknown tumor status. In 159 patients, multiple tumors were sampled, and metachronous tumors were collected in 83 patients. Between 2016 and 2020 the UROSCAN-biobanking included 1122/2999 (37%) of all primary bladder cancer patients in the Southern Healthcare Region. Until June 2021, the corresponding numbers subjected to RNA-sequencing and molecular subtyping was 605 (UROSCANSEQ), of which 52 (9%) samples were not sequenced due to inadequate RNA-quality (n = 47) or technical failure/lost sample (n = 5). CONCLUSIONS: The UROSCAN-biobanking and UROSCANSEQ-infrastructure for molecular subtyping by real-time RNA-sequencing represents, to our knowledge, the largest effort of evaluating population-wide molecular classification of bladder cancer.

Single-cell sequencing in translational cancer research and challenges to meet clinical diagnostic needs.

The ability to capture alterations in the genome or transcriptome by next-generation sequencing has provided critical insight into molecular changes and programs underlying cancer biology. With the rapid technological development in single-cell sequencing, it has become possible to study individual cells at the transcriptional, genetic, epigenetic, and protein level. Using single-cell analysis, an increased resolution of fundamental processes underlying cancer development is obtained, providing comprehensive insights otherwise lost by sequencing of entire (bulk) samples, in which molecular signatures of individual cells are averaged across the entire cell population. Here, we provide a concise overview on the application of single-cell analysis of different modalities within cancer research by highlighting key articles of their respective fields. We furthermore examine the potential of existing technologies to meet clinical diagnostic needs and discuss current challenges associated with this translation.

Silencing of microRNA families by seed-targeting tiny LNAs.

The challenge of understanding the widespread biological roles of animal microRNAs (miRNAs) has prompted the development of genetic and functional genomics technologies for miRNA loss-of-function studies. However, tools for exploring the functions of entire miRNA families are still limited. We developed a method that enables antagonism of miRNA function using seed-targeting 8-mer locked nucleic acid (LNA) oligonucleotides, termed tiny LNAs. Transfection of tiny LNAs into cells resulted in simultaneous inhibition of miRNAs within families sharing the same seed with concomitant upregulation of direct targets. In addition, systemically delivered, unconjugated tiny LNAs showed uptake in many normal tissues and in breast tumors in mice, coinciding with long-term miRNA silencing. Transcriptional and proteomic profiling suggested that tiny LNAs have negligible off-target effects, not significantly altering the output from mRNAs with perfect tiny LNA complementary sites. Considered together, these data support the utility of tiny LNAs in elucidating the functions of miRNA families in vivo.

Two genetic pathways, t(1;10) and amplification of 3p11-12, in myxoinflammatory fibroblastic sarcoma, haemosiderotic fibrolipomatous tumour, and morphologically similar lesions.

Myxoinflammatory fibroblastic sarcoma (MIFS) is a low-grade malignant neoplasm for which limited genetic information, including a t(1;10)(p22;q24) and amplification of chromosome 3 material, is available. To further characterize these aberrations, we have investigated eight soft tissue sarcomas diagnosed as MIFS, haemosiderotic fibrolipomatous tumour (HFT), myxoid spindle cell/pleomorphic sarcoma with MIFS features, and inflammatory malignant fibrous histiocytoma/undifferentiated pleomorphic sarcoma with prominent inflammation (IMFH) harbouring a t(1;10) or variants thereof and/or ring chromosomes with possible involvement of chromosome 3. Using chromosome banding, fluorescence in situ hybridization, array-based comparative genomic hybridization, global gene expression, and real-time quantitative PCR analyses, we identified the breakpoint regions on chromosomes 1 and 10, demonstrated and delineated the commonly amplified region on chromosome 3, and assessed the consequences of these alterations for gene expression. The breakpoints in the t(1;10) mapped to TGFBR3 in 1p22 and in or near MGEA5 in 10q24, resulting in transcriptional up-regulation of NPM3 and particularly FGF8, two consecutive genes located close to MGEA5. The ring chromosomes contained a commonly amplified 1.44 Mb region in 3p11-12, which was associated with increased expression of VGLL3 and CHMP2B. The identified genetic aberrations were not confined to MIFS; an identical t(1;10) was also found in a case of HFT and the amplicon in 3p was seen in an IMFH.

Tiling resolution array CGH and high density expression profiling of urothelial carcinomas delineate genomic amplicons and candidate target genes specific for advanced tumors.

BACKGROUND: Urothelial carcinoma (UC) is characterized by nonrandom chromosomal aberrations, varying from one or a few changes in early-stage and low-grade tumors, to highly rearranged karyotypes in muscle-invasive lesions. Recent array-CGH analyses have shed further light on the genomic changes underlying the neoplastic development of UC, and have facilitated the molecular delineation amplified and deleted regions to the level of specific candidate genes. In the present investigation we combine detailed genomic information with expression information to identify putative target genes for genomic amplifications. METHODS: We analyzed 38 urothelial carcinomas by whole-genome tiling resolution array-CGH and high density expression profiling to identify putative target genes in common genomic amplifications. When necessary expression profiling was complemented with Q-PCR of individual genes. RESULTS: Three genomic segments were frequently and exclusively amplified in high grade tumors; 1q23, 6p22 and 8q22, respectively. Detailed mapping of the 1q23 segment showed a heterogeneous amplification pattern and no obvious commonly amplified region. The 6p22 amplicon was defined by a 1.8 Mb core region present in all amplifications, flanked both distally and proximally by segments amplified to a lesser extent. By combining genomic profiles with expression profiles we could show that amplification of E2F3, CDKAL1, SOX4, and MBOAT1 as well as NUP153, AOF1, FAM8A1 and DEK in 6p22 was associated with increased gene expression. Amplification of the 8q22 segment was primarily associated with YWHAZ (14-3-3-zeta) and POLR2K over expression. The possible importance of the YWHA genes in the development of urothelial carcinomas was supported by another recurrent amplicon paralogous to 8q22, in 2p25, where increased copy numbers lead to enhanced expression of YWHAQ (14-3-3-theta). Homozygous deletions were identified at 10 different genomic locations, most frequently affecting CDKN2A/CDKN2B in 9p21 (32%). Notably, the latter occurred mutually exclusive with 6p22 amplifications. CONCLUSION: The presented data indicates 6p22 as a composite amplicon with more than one possible target gene. The data also suggests that amplification of 6p22 and homozygous deletions of 9p21 may have complementary roles. Furthermore, the analysis of paralogous regions that showed genomic amplification indicated altered expression of YWHA (14-3-3) genes as important events in the development of UC.

Tiling resolution array comparative genomic hybridization, expression and methylation analyses of dup(1q) in Burkitt lymphomas and pediatric high hyperdiploid acute lymphoblastic leukemias reveal clustered near-centromeric breakpoints and overexpression of genes in 1q22-32.3.

Although gain of 1q occurs in 25% of Burkitt lymphomas (BLs) and 10% of pediatric high hyperdiploid acute lymphoblastic leukemias (ALLs), little is known about the origin, molecular genetic characteristics and functional outcome of dup(1q) in these disorders. Ten dup(1q)-positive BLs/ALLs were investigated by tiling resolution (32k) array CGH analysis, which revealed that the proximal breakpoints in all cases were near-centromeric, in eight of them clustering within a 1.4 Mb segment in 1q12-21.1. The 1q distal breakpoints were heterogeneous, being more distal in the ALLs than in the BLs. The minimally gained segments in the ALLs and BLs were 57.4 Mb [dup(1)(q22q32.3)] and 35 Mb [dup(1)(q12q25.2)], respectively. Satellite II DNA on 1q was not hypomethylated, as ascertained by Southern blot analyses of 15 BLs/ALLs with and without gain of 1q, indicating that aberrant methylation was not involved in the origin of dup(1q), as previously suggested for other neoplasms with 1q rearrangements. Global gene expression analyses revealed that five genes in the minimally 57.4 Mb gained region--B4GALT3, DAP3, RGS16, TMEM183A and UCK2--were significantly overexpressed in dup(1q)-positive ALLs compared with high hyperdiploid ALLs without dup(1q). The DAP3 and UCK2 genes were among the most overexpressed genes in the BL case with gain of 1q investigated. The DAP3 protein has been reported to be highly expressed in invasive glioblastoma multiforme cells, whereas expression of the UCK2 protein has been correlated with sensitivity to anticancer drugs. However, involvement of these genes in dup(1q)-positive ALLs and BLs has previously not been reported.

High-resolution genomic profiles of breast cancer cell lines assessed by tiling BAC array comparative genomic hybridization.

A BAC-array platform for comparative genomic hybridization was constructed from a library of 32,433 clones providing complete genome coverage, and evaluated by screening for DNA copy number changes in 10 breast cancer cell lines (BT474, MCF7, HCC1937, SK-BR-3, L56Br-C1, ZR-75-1, JIMT1, MDA-MB-231, MDA-MB-361, and HCC2218) and one cell line derived from fibrocystic disease of the breast (MCF10A). These were also characterized by gene expression analysis and found to represent all five recently described breast cancer subtypes using the "intrinsic gene set" and centroid correlation. Three cell lines, HCC1937 and L56BrC1 derived from BRCA1 mutation carriers and MDA-MB-231, were of basal-like subtype and characterized by a high frequency of low-level gains and losses of typical pattern, including limited deletions on 5q. Four estrogen receptor positive cell lines were of luminal A subtype and characterized by a different pattern of aberrations and high-level amplifications, including ERBB2 and other 17q amplicons in BT474 and MDA-MB-361. SK-BR-3 cells, characterized by a complex genome including ERBB2 amplification, massive high-level amplifications on 8q and a homozygous deletion of CDH1 at 16q22, had an expression signature closest to luminal B subtype. The effects of gene amplifications were verified by gene expression analysis to distinguish targeted genes from silent amplicon passengers. JIMT1, derived from an ERBB2 amplified trastuzumab resistant tumor, was of the ERBB2 subtype. Homozygous deletions included other known targets such as PTEN (HCC1937) and CDKN2A (MDA-MB-231, MCF10A), but also new candidate suppressor genes such as FUSSEL18 (HCC1937) and WDR11 (L56Br-C1) as well as regions without known genes. The tiling BAC-arrays constitute a powerful tool for high-resolution genomic profiling suitable for cancer research and clinical diagnostics.

Tiling resolution array comparative genomic hybridization analysis of a fibrosarcoma of bone.

Fibrosarcoma of bone is a rare malignant tumor accounting for less than 5% of all primary malignant bone neoplasms. There is very limited knowledge regarding the molecular genetics of this tumor, and there are no cytogenetic data available. In the present study, a fibrosarcoma deriving from the left iliac bone of a 10-year-old girl was characterized using cytogenetics, fluorescence in situ hybridization (FISH), and whole genome tiling resolution array comparative genomic hybridization (CGH). Cytogenetic and FISH analyses revealed a ring chromosome 6 as the sole acquired aberration, a finding corroborated by array CGH. The ring formation, however, did not result in any gain of genetic material. Nor did the breakpoints in 6p25 and 6q14 seem to affect any known gene loci in such a way that the ring formation could have resulted in the creation of a fusion gene or in the exchange of regulatory sequences. Thus, a reasonable interpretation of the pathogenetic significance of the ring formation would be that it resulted in the loss of one or more putative tumor suppressor gene loci distal to the two breakpoints.

Microarray analyses reveal strong influence of DNA copy number alterations on the transcriptional patterns in pancreatic cancer: implications for the interpretation of genomic amplifications.

DNA copy number alterations are believed to play a major role in the development and progression of human neoplasms. Although most of these genomic imbalances have been associated with dysregulation of individual genes, their large-scale transcriptional consequences remain unclear. Pancreatic carcinomas frequently display gene copy number variation of entire chromosomes as well as of chromosomal subregions. These changes range from homozygous deletions to high-level amplifications and are believed to constitute key genetic alterations in the cellular transformation of this tumor type. To investigate the transcriptional consequences of the most drastic genomic changes, that is, genomic amplifications, and to analyse the genome-wide transcriptional effects of DNA copy number changes, we performed expression profiling of 29 pancreatic carcinoma cell lines and compared the results with matching genomic profiling data. We show that a strong association between DNA copy numbers and mRNA expression levels is present in pancreatic cancer, and demonstrate that as much as 60% of the genes within highly amplified genomic regions display associated overexpression. Consequently, we identified 67 recurrently overexpressed genes located in seven precisely mapped commonly amplified regions. The presented findings indicate that more than one putative target gene may be of importance in most pancreatic cancer amplicons.

Genome-wide array-based comparative genomic hybridization reveals multiple amplification targets and novel homozygous deletions in pancreatic carcinoma cell lines.

Pancreatic carcinomas display highly complex chromosomal abnormalities, including many structural and numerical aberrations. There is ample evidence indicating that some of these abnormalities, such as recurrent amplifications and homozygous deletions, contribute to tumorigenesis by altering expression levels of critical oncogenes and tumor suppressor genes. To increase the understanding of gene copy number changes in pancreatic carcinomas and to identify key amplification/deletion targets, we applied genome-wide array-based comparative genomic hybridization to 31 pancreatic carcinoma cell lines. Two different microarrays were used, one containing 3,565 fluorescence in situ hybridization-verified bacterial artificial chromosome clones and one containing 25,468 cDNA clones representing 17,494 UniGene clusters. Overall, the analyses revealed a high genomic complexity, with several copy number changes detected in each case. Specifically, 60 amplicons at 32 different locations were identified, most frequently located within 8q (8 cases), 12p (7 cases), 7q (5 cases), 18q (5 cases), 19q (5 cases), 6p (4 cases), and 8p (4 cases). Amplifications of 8q and 12p were mainly clustered at 8q23-24 and 12p11-12, respectively, whereas amplifications on other chromosome arms were more dispersed. Furthermore, our analyses identified several novel homozygously deleted segments located to 9p24, 9p21, 9q32, 10p12, 10q22, 12q24, and 18q23. The individual complexity and aberration patterns varied substantially among cases, i.e., some cell lines were characterized mainly by high-level amplifications, whereas others showed primarily whole-arm imbalances and homozygous deletions. The described amplification and deletion targets are likely to contain genes important in pancreatic tumorigenesis.

Pancreatic carcinoma cell lines with SMAD4 inactivation show distinct expression responses to TGFB1.

Transforming growth factor beta-1 (TGFB1)-induced gene expression was studied in five pancreatic carcinoma cell lines and one known TGFB1-sensitive cell line (HaCaT) by use of high-density filter-based cDNA microarrays representing over 4,000 human genes. The results indicate a complex cellular response to TGFB1 with 10% of the investigated genes showing altered expression after 3 or 48 hr of TGFB1 exposure. The tumor cell lines displayed a gradually inversed gene expression pattern, which correlated with reduced sensitivity to TGFB1, as compared to the HaCaT cell line. In the HaCaT cells, several proapoptotic genes showed increased expression in response to TGFB1, whereas the expression of antiapoptotic genes was decreased. In contrast, two pancreatic carcinoma cell lines, previously found to be growth stimulated by TGFB1, displayed an expression pattern opposite to that of these genes. Similarly, the expression of other functional groups of genes, such as cell cycle and transcription factor related genes, was almost completely reversed in these two tumor cell lines. Importantly, three of the five investigated pancreatic carcinoma cell lines responded to TGFB1, although they had SMAD4 inactivations, suggesting that the observed gene expression changes in these cell lines must be accomplished by SMAD-independent pathways.

Detailed genomic mapping and expression analyses of 12p amplifications in pancreatic carcinomas reveal a 3.5-Mb target region for amplification.

Previous cytogenetic and comparative genomic hybridization (CGH) analyses have shown that the gain of chromosome arm 12p is frequent in pancreatic carcinomas. We investigated 15 pancreatic carcinoma cell lines using CGH, fluorescence in situ hybridization (FISH), and semiquantitative polymerase chain reaction (PCR) to characterize 12p amplifications in detail. The CGH analysis revealed gains of 12p in four of the cell lines and local amplification within 12p11-12 in six cell lines. By FISH analysis, using precisely mapped YAC clones, the commonly amplified region was found to be approximately 5 Mb. The amplified segment extended from YAC 753f12, covering the KRAS2 locus, to YAC 891f1, close to the centromere. A semiquantitative PCR methodology was used to estimate genomic copy numbers of 14 precisely mapped expressed sequence tags (ESTs) and sequence-tagged sites, located within this interval. The level of amplification ranged from two- to 12-fold. The produced gene copy profiles revealed a 3.5-Mb segment with various local amplifications. This region includes KRAS2 and ranges from D12S1617 to sts-N38796. Two of the cell lines (primary and metastatic tumor from the same patient) showed amplification peaks within the distal region of this segment, two had peaks within the proximal region, one showed subpeaks in both regions, and one displayed amplification of the entire region. Chromosome segment-specific cDNA array analysis of 29 expressed sequences within the whole interval between D12S1617 and sts-N38796 indicated overexpression of four ESTs, two corresponding to DEC2 and PPFIBP1, and two to ESTs with unknown function. Expression analysis of these and of KRAS2 showed specific overexpression in the six cell lines with local 12p amplifications. These findings indicate two target regions within the 3.5-Mb segment in 12p11-12, one proximal including PPFIBP1, and one distal including KRAS2.