Chromophobe renal cell carcinoma (chRCC) is a renal tumor subtype with a good prognosis, characterized by multiple chromosomal copy number variations (CNV). The World Health Organization (WHO) chRCC classification guidelines define a classic and an eosinophilic variant. Large cells with reticular cytoplasm and prominent cell membranes (pale cells) are characteristic for classic chRCC. Classic and eosinophilic variants were defined in 42 Swiss chRCCs, 119 Japanese chRCCs and in whole-slide digital images of 66 chRCCs from the Cancer Genome Atlas (TCGA) kidney chromophobe (KICH) dataset. 32 of 42 (76.2%) Swiss chRCCs, 90 of 119 (75.6%) Japanese chRCCs and 53 of 66 (80.3%) TCGA-KICH were classic chRCCs. There was no survival difference between eosinophilic and classic chRCC in all three cohorts. To identify a genotype/phenotype correlation, we performed a genome-wide CNV analysis using Affymetrix OncoScan® CNV Assay (Affymetrix/Thermo Fisher Scientific, Waltham, MA, USA) in 33 Swiss chRCCs. TCGA-KICH subtypes were compared with TCGA CNV data. In the combined Swiss and TCGA-KICH cohorts, losses of chromosome 1, 2, 6, 10, 13, and 17 were significantly more frequent in classic chRCC (p < 0.05, each), suggesting that classic chRCC are characterized by higher chromosomal instability. This molecular difference justifies the definition of two chRCC variants. Absence of pale cells could be used as main histological criterion to define the eosinophilic variant of chRCC.
Approximately 70% of clear cell renal cell carcinoma (ccRCC) is characterized by the biallelic inactivation of von Hippel-Lindau (VHL) on chromosome 3p. ELOC-mutated (Elongin C-mutated) renal cell carcinoma containing biallelic ELOC inactivations with chromosome 8q deletions is considered a novel subtype of renal cancer possessing a morphologic overlap with ccRCC, renal cell carcinoma (RCC) with fibromyomatous stroma exhibiting Tuberous Sclerosis Complex (TSC)/mammalian Target of Rapamycin (mTOR) mutations, and clear cell papillary tumor. However, the frequency and consequences of ELOC alterations in wild-type VHL and mutated VHL RCC are unclear. In this study, we characterize 123 renal tumors with clear cell morphology and known VHL mutation status to assess the morphologic and molecular consequences of ELOC inactivation. Using OncoScan and whole-exome sequencing, we identify 18 ELOC-deleted RCCs, 3 of which contain ELOC mutations resulting in the biallelic inactivation of ELOC. Biallelic ELOC and biallelic VHL aberrations were mutually exclusive; however, 2 ELOC-mutated RCCs showed monoallelic VHL alterations. Furthermore, no mutations in TSC1, TSC2, or mTOR were identified in ELOC-mutated RCC with biallelic ELOC inactivation. Using High Ambiguity Driven biomolecular DOCKing, we report a novel ELOC variant containing a duplication event disrupting ELOC-VHL interaction alongside the frequently seen Y79C alteration. Using hyper reaction monitoring mass spectrometry, we show RCCs with biallelic ELOC alterations have significantly reduced ELOC expression but similar carbonic anhydrase 9 and vascular endothelial growth factor A expression compared with classical ccRCC with biallelic VHL inactivation. The absence of biallelic VHL and TSC1, TSC2, or mTOR inactivation in RCC with biallelic ELOC inactivation (ELOC mutation in combination with ELOC deletions on chromosome 8q) supports the notion of a novel, molecularly defined tumor entity.
Chromophobe renal cell carcinoma (chRCC) patients have good prognosis. Only 5%–10% patients die of metastatic disease after tumorectomy, but tumor progression cannot be predicted by histopathological parameters alone. chRCC are characterized by losses of many chromosomes, whereas gene mutations are rare. In this study, we aim at identifying genes indicating chRCC progression. A bioinformatic approach was used to correlate chromosomal loss and mRNA expression from 15287 genes from The Cancer Genome Atlas (TCGA) database. All genes in TCGA chromophobe renal cancer dataset (KICH) for which a significant correlation between chromosomal loss and mRNA expression was shown, were identified and their associations with outcome was assessed. Genome-wide DNA copy-number alterations were analyzed by Affymetrix OncoScan® CNV FFPE Microarrays in a second cohort of Swiss chRCC. In both cohorts, tumors with loss of chromosomes 2, 6, 10, 13, 17 and 21 had signs of tumor progression. There were 4654 genes located on these chromosomes, and 13 of these genes had reduced mRNA levels, which was associated with poor outcome in chRCC. Decreased CDKN1A expression at mRNA (p = 0.02) and protein levels (p = 0.02) were associated with short overall survival and were independent predictors of prognosis (p < 0.01 and <0.05 respectively). CDKN1A expression status is a prognostic biomarker independent of tumor stage. CDKN1A immunohistochemistry may be used to identify chRCC patients at greater risk of disease progression.
The large amount and diversity of viral genomic datasets generated by next-generation sequencing technologies poses a set of challenges for computational data analysis workflows, including rigorous quality control, scaling to large sample sizes, and tailored steps for specific applications. Here, we present V-pipe 3.0, a computational pipeline designed for analyzing next-generation sequencing data of short viral genomes. It is developed to enable reproducible, scalable, adaptable, and transparent inference of genetic diversity of viral samples. By presenting 2 large-scale data analysis projects, we demonstrate the effectiveness of V-pipe 3.0 in supporting sustainable viral genomic data science.
SARS-CoV-2, the virus responsible for the current COVID-19 pandemic, is evolving into different genetic variants by accumulating mutations as it spreads globally. In addition to this diversity of consensus genomes across patients, RNA viruses can also display genetic diversity within individual hosts, and co-existing viral variants may affect disease progression and the success of medical interventions. To systematically examine the intra-patient genetic diversity of SARS-CoV-2, we processed a large cohort of 3939 publicly-available deeply sequenced genomes with specialised bioinformatics software, along with 749 recently sequenced samples from Switzerland. We found that the distribution of diversity across patients and across genomic loci is very unbalanced with a minority of hosts and positions accounting for much of the diversity. For example, the D614G variant in the Spike gene, which is present in the consensus sequences of 67.4% of patients, is also highly diverse within hosts, with 29.7% of the public cohort being affected by this coexistence and exhibiting different variants. We also investigated the impact of several technical and epidemiological parameters on genetic heterogeneity and found that age, which is known to be correlated with poor disease outcomes, is a significant predictor of viral genetic diversity. Author Summary: Since it arose in late 2019, the new coronavirus (SARS-CoV-2) behind the COVID-19 pandemic has mutated and evolved during its global spread. Individual patients may host different versions, or variants, of the virus, hallmarked by different mutations. We examine the diversity of genetic variants coexisting within patients across a cohort of 3939 publicly accessible samples and 749 recently sequenced samples from Switzerland. We find that a small number of patients carry most of the diversity, and that patients with more diversity tend to be older. We also find that most of the diversity is concentrated in certain regions and positions of the virus genome. In particular, we find that a variant reported to increase infectivity is among the most diverse positions. Our study provides a large-scale survey of within-patient diversity of the SARS-CoV-2 genome.