SummaryIncidence and mortality for sex-unspecific cancers is higher among men and is largely unexplained1,2. Furthermore, age-related loss of chromosome Y (LOY) is frequent in normal haematopoietic cells3,4, but the phenotypic consequences of LOY have been elusive5–10. From analysis of 1153 elderly men, we report that LOY was associated with risks of all-cause mortality (HR=1.91, 95% CI=1.17-3.13, events=637) and non-haematological cancer mortality (HR=3.62, CI=1.56-8.41, events=132). LOY affected at least 8.2% of subjects in this cohort and median survival among men with LOY was 5.5 years shorter. Risk of all-cause mortality and LOY was validated in an independent cohort (HR=3.66), in which 20.5% of subjects displayed LOY. These results illustrate the impact of post-zygotic mosaicism on disease risk, could explain why males are more frequently affected by cancer and suggest that chromosome Y is important in processes beyond sex determination. LOY in blood could become a predictive biomarker of male carcinogenesis.

The exploration of copy-number variation (CNV), notably of somatic cells, is an understudied aspect of genome biology. Any differences in the genetic makeup between twins derived from the same zygote represent an irrefutable example of somatic mosaicism. We studied 19 pairs of monozygotic twins with either concordant or discordant phenotype by using two platforms for genome-wide CNV analyses and showed that CNVs exist within pairs in both groups. These findings have an impact on our views of genotypic and phenotypic diversity in monozygotic twins and suggest that CNV analysis in phenotypically discordant monozygotic twins may provide a powerful tool for identifying disease-predisposition loci. Our results also imply that caution should be exercised when interpreting disease causality of de novo CNVs found in patients based on analysis of a single tissue in routine disease-related DNA diagnostics.

Two major types of genetic variation are known: single nucleotide polymorphisms (SNPs), and a more recently discovered structural variation, involving changes in copy number (CNVs) of kilobase- to megabase-sized chromosomal segments. It is unknown whether CNVs arise in somatic cells, but it is, however, generally assumed that normal cells are genetically identical. We tested 34 tissue samples from three subjects and, having analyzed for each tissue < or =10(-6) of all cells expected in an adult human, we observed at least six CNVs, affecting a single organ or one or more tissues of the same subject. The CNVs ranged from 82 to 176 kb, often encompassing known genes, potentially affecting gene function. Our results indicate that humans are commonly affected by somatic mosaicism for stochastic CNVs, which occur in a substantial fraction of cells. The majority of described CNVs were previously shown to be polymorphic between unrelated subjects, suggesting that some CNVs previously reported as germline might represent somatic events, since in most studies of this kind, only one tissue is typically examined and analysis of parents for the studied subjects is not routinely performed. A considerable number of human phenotypes are a consequence of a somatic process. Thus, our conclusions will be important for the delineation of genetic factors behind these phenotypes. Consequently, biobanks should consider sampling multiple tissues to better address mosaicism in the studies of somatic disorders.

Whole-genome sequencing (WGS) is a fundamental technology for research to advance precision medicine, but the limited availability of portable and user-friendly workflows for WGS analyses poses a major challenge for many research groups and hampers scientific progress. Here we present Sarek, an open-source workflow to detect germline variants and somatic mutations based on sequencing data from WGS, whole-exome sequencing (WES), or gene panels. Sarek features (i) easy installation, (ii) robust portability across different computer environments, (iii) comprehensive documentation, (iv) transparent and easy-to-read code, and (v) extensive quality metrics reporting. Sarek is implemented in the Nextflow workflow language and supports both Docker and Singularity containers as well as Conda environments, making it ideal for easy deployment on any POSIX-compatible computers and cloud compute environments. Sarek follows the GATK best-practice recommendations for read alignment and pre-processing, and includes a wide range of software for the identification and annotation of germline and somatic single-nucleotide variants, insertion and deletion variants, structural variants, tumour sample purity, and variations in ploidy and copy number. Sarek offers easy, efficient, and reproducible WGS analyses, and can readily be used both as a production workflow at sequencing facilities and as a powerful stand-alone tool for individual research groups. The Sarek source code, documentation and installation instructions are freely available at https://github.com/nf-core/sarek and at https://nf-co.re/sarek/.

MicroRNAs (miRNAs) are small RNAs (w22 bp) that post-transcriptionally regulate protein expression and are found to be differentially expressed in a number of human cancers. There is increasing evidence to suggest that miRNAs could be useful in cancer diagnosis, prognosis, and therapy. We performed miRNA microarray expression profiling on a cohort of 12 benign and 12 malignant pheochromocytomas and identified a number of differentially expressed miRNAs. These results were validated in a separate cohort of ten benign and ten malignant samples using real-time quantitative RT-PCR; benign samples had a minimum follow-up of at least 2 years. It was found that IGF2 as well as its intronic miR-483-5p was over-expressed, while miR-15a and miR-16 were under-expressed in malignant tumours compared with benign tumours. These miRNAs were found to be diagnostic and prognostic markers for malignant pheochromocytoma. The functional role of miR-15a and miR-16 was investigated in vitro in the rat PC12 pheochromocytoma cell line, and these miRNAs were found to regulate cell proliferation via their effect on cyclin D1 and apoptosis. These data indicate that miRNAs play a pivotal role in the biology of malignant pheochromocytoma, and represent an important class of diagnostic and prognostic biomarkers and therapeutic targets warranting further investigation.

To further explore the extent of structural large-scale variation in the human genome, we assessed copy number variations (CNVs) in a series of 71 healthy subjects from three ethnic groups. CNVs were analyzed using comparative genomic hybridization (CGH) to a BAC array covering the human genome, using DNA extracted from peripheral blood, thus avoiding any culture-induced rearrangements. By applying a newly developed computational algorithm based on Hidden Markov modeling, we identified 1,078 autosomal CNVs, including at least two neighboring/overlapping BACs, which represent 315 distinct regions. The average size of the sequence polymorphisms was approximately 350 kb and involved in total approximately 117 Mb or approximately 3.5% of the genome. Gains were about four times more common than deletions, and segmental duplications (SDs) were overrepresented, especially in larger deletion variants. This strengthens the notion that SDs often define hotspots of chromosomal rearrangements. Over 60% of the identified autosomal rearrangements match previously reported CNVs, recognized with various platforms. However, results from chromosome X do not agree well with the previously annotated CNVs. Furthermore, data from single BACs deviating in copy number suggest that our above estimate of total variation is conservative. This report contributes to the establishment of the common baseline for CNV, which is an important resource in human genetics.

Pheochromocytomas and abdominal paragangliomas are adrenal and extra-adrenal catecholamine-producing tumours. They arise due to heritable cancer syndromes, or more frequently occur sporadically due to an unknown genetic cause. The majority of cases are benign, but malignant tumours are observed. Previous comparative genomic hybridization (CGH) and loss of heterozygosity studies have shown frequent deletions of chromosome arms 1p, 3q and 22q in pheochromocytomas. We applied high-resolution whole-genome array CGH on 53 benign and malignant pheochromocytomas and paragangliomas to narrow down candidate regions as well as to identify chromosomal alterations more specific to malignant tumours. Minimal overlapping regions (MORs) were identified on 16 chromosomes, with the most frequent MORs of deletion (R32%) occurring on chromosome arms 1p, 3q, 11p/q, 17p and 22q, while the chromosome arms 1q, 7p, 12q and 19p harboured the most common MORs of gain (R14%). The most frequent MORs (61-75%) in the pheochromocytomas were identified at 1p, and the four regions of common losses encompassed 1p36, 1p32-31, 1p22-21 and 1p13. Tumours that did not show 1p loss generally demonstrated aberrations on chromosome 11. Gain of chromosomal material was significantly more frequent among the malignant cases. Moreover, gain at 19q, trisomy 12 and loss at 11q were positively associated with malignant pheochromocytomas, while 1q gain was commonly observed in the malignant paragangliomas. Our study revealed novel and narrow recurrent chromosomal regions of loss and gain at several autosomes, a prerequisite for identifying candidate tumour suppressor genes and oncogenes involved in the development of adrenal and extra-adrenal catecholamine-producing tumours.

Breast cancer is a major cause of morbidity and mortality in women and its metastatic spread is the principal reason behind the fatal outcome. Metastasis-related research of breast cancer is however underdeveloped when compared with the abundant literature on primary tumors. We applied an unexplored approach comparing at high resolution the genomic profiles of primary tumors and synchronous axillary lymph node metastases from 13 patients with breast cancer. Overall, primary tumors displayed 20% higher number of aberrations than metastases. In all but two patients, we detected in total 157 statistically significant differences between primary lesions and matched metastases. We further observed differences that can be linked to metastatic disease and there was also an overlapping pattern of changes between different patients. Many of the differences described here have been previously linked to poor patient survival, suggesting that this is a viable approach toward finding biomarkers for disease progression and definition of new targets useful for development of anticancer drugs. Frequent genetic differences between primary tumors and metastases in breast cancer also question, at least to some extent, the role of primary tumors as a surrogate subject of study for the systemic disease.