Author Correction: Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors.

The original version of this Article omitted a declaration from the competing interests statement, which should have included the following: 'K.P.W. is President of Tempus Lab, Inc., Chicago, IL, USA'. This has now been corrected in both the PDF and HTML versions of the Article.

The Brussels Infant and Toddler Stool Scale: A Study on Interobserver Reliability.

OBJECTIVES: The Bristol Stool Form Scale (BSFS) is inadequate for non-toilet trained children. The Brussels Infant and Toddler Stool Scale (BITSS) was developed, consisting of 7 photographs of diapers containing stools of infants and toddlers. We aimed to evaluate interobserver reliability of stool consistency assessment among parents, nurses, and medical doctors (MDs) using the BITSS. METHODS: In this multicenter cross-sectional study (2016-2017), BITSS photographs were rated according to the BSFS. The reliability of the BITSS was evaluated using the overall proportion of perfect agreement and the linearly weighted κ statistic. RESULTS: A total of 2462 observers participated: 1181 parents (48.0%), 624 nurses (25.3%), and 657 MDs (26.7%). The best-performing BITSS photographs corresponded with BSFS type 7 (87.5%) and type 4 (87.6%), followed by the BITSS photographs representing BSFS type 6 (75.0%), BSFS type 5 (68.0%), BSFS type 1 (64.8%), and BSFS type 3 (64.6%). The weakest performing BITSS photograph corresponded with BSFS type 2 (49.7%). The overall weighted κ-value was 0.72 (95% CI 0.59-0.85; good agreement). Based on these results, photographs were categorized per stool group as hard (BSFS type 1-3), formed (BSFS type 4), loose (BSFS types 5 and 6), or watery (BSFS type 7) stools. According to this new categorization system, correct allocation for each photograph ranged from 83 to 96% (average: 90%). The overall proportion of correct allocations was 72.8%. CONCLUSIONS: BITSS showed good agreement with BSFS. Using the newly categorized BITSS photographs, the BITSS is reliable for the assessment of stools of non-toilet trained children in clinical practice and research. A multilanguage translated version of the BITSS can be downloaded at https://bitss-stoolscale.com/.

Author Correction: Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors.

The original version of this Article contained an error in the author affiliations. The affiliation of Kevin P. White with Tempus Labs, Inc., Chicago, IL, USA was inadvertently omitted.This has now been corrected in both the PDF and HTML versions of the Article.

Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors.

Homozygous deletions are rare in cancers and often target tumour suppressor genes. Here, we build a compendium of 2218 primary tumours across 12 human cancer types and systematically screen for homozygous deletions, aiming to identify rare tumour suppressors. Our analysis defines 96 genomic regions recurrently targeted by homozygous deletions. These recurrent homozygous deletions occur either over tumour suppressors or over fragile sites, regions of increased genomic instability. We construct a statistical model that separates fragile sites from regions showing signatures of positive selection for homozygous deletions and identify candidate tumour suppressors within those regions. We find 16 established tumour suppressors and propose 27 candidate tumour suppressors. Several of these genes (including MGMT, RAD17, and USP44) show prior evidence of a tumour suppressive function. Other candidate tumour suppressors, such as MAFTRR, KIAA1551, and IGF2BP2, are novel. Our study demonstrates how rare tumour suppressors can be identified through copy number meta-analysis.

Concurrent whole-genome haplotyping and copy-number profiling of single cells.

Methods for haplotyping and DNA copy-number typing of single cells are paramount for studying genomic heterogeneity and enabling genetic diagnosis. Before analyzing the DNA of a single cell by microarray or next-generation sequencing, a whole-genome amplification (WGA) process is required, but it substantially distorts the frequency and composition of the cell's alleles. As a consequence, haplotyping methods suffer from error-prone discrete SNP genotypes (AA, AB, BB) and DNA copy-number profiling remains difficult because true DNA copy-number aberrations have to be discriminated from WGA artifacts. Here, we developed a single-cell genome analysis method that reconstructs genome-wide haplotype architectures as well as the copy-number and segregational origin of those haplotypes by employing phased parental genotypes and deciphering WGA-distorted SNP B-allele fractions via a process we coin haplarithmisis. We demonstrate that the method can be applied as a generic method for preimplantation genetic diagnosis on single cells biopsied from human embryos, enabling diagnosis of disease alleles genome wide as well as numerical and structural chromosomal anomalies. Moreover, meiotic segregation errors can be distinguished from mitotic ones.

Single cell segmental aneuploidy detection is compromised by S phase.

BACKGROUND: Carriers of balanced translocations are at high risk for unbalanced gametes which can result in recurrent miscarriages or birth defects. Preimplantation genetic diagnosis (PGD) is often offered to select balanced embryos. This selection is currently mainly performed by array CGH on blastomeres. Current methodology does not take into account the phase of the cell cycle, despite the variable copy number status of different genomic regions in S phase. RESULTS: Cell lines derived from 3 patients with different chromosomal imbalances were used to evaluate the accuracy of single cell array CGH. The different cell cycle phases were sorted by flow cytometry and 10 single cells were picked per cell line per cell cycle phase, whole genome amplified and analyzed by BAC arrays, the most commonly used platform for PGD purposes. In contrast to G phase, where the imbalances were efficiently identified, less than half of the probes in the regions of interest indicated the presence of the aberration in 17 S-phase cells, resulting in reduced accuracy. CONCLUSIONS: The results demonstrate that the accuracy to detect segmental chromosomal imbalances is reduced in S-phase cells, which could be a source of misdiagnosis in PGD. Hence, the cell cycle phase of the analyzed cell is of great importance and should be taken into account during the analysis. This knowledge may guide future technological improvements.

Constitutional and somatic rearrangement of chromosome 21 in acute lymphoblastic leukaemia.

Changes in gene dosage are a major driver of cancer, known to be caused by a finite, but increasingly well annotated, repertoire of mutational mechanisms. This can potentially generate correlated copy-number alterations across hundreds of linked genes, as exemplified by the 2% of childhood acute lymphoblastic leukaemia (ALL) with recurrent amplification of megabase regions of chromosome 21 (iAMP21). We used genomic, cytogenetic and transcriptional analysis, coupled with novel bioinformatic approaches, to reconstruct the evolution of iAMP21 ALL. Here we show that individuals born with the rare constitutional Robertsonian translocation between chromosomes 15 and 21, rob(15;21)(q10;q10)c, have approximately 2,700-fold increased risk of developing iAMP21 ALL compared to the general population. In such cases, amplification is initiated by a chromothripsis event involving both sister chromatids of the Robertsonian chromosome, a novel mechanism for cancer predisposition. In sporadic iAMP21, breakage-fusion-bridge cycles are typically the initiating event, often followed by chromothripsis. In both sporadic and rob(15;21)c-associated iAMP21, the final stages frequently involve duplications of the entire abnormal chromosome. The end-product is a derivative of chromosome 21 or the rob(15;21)c chromosome with gene dosage optimized for leukaemic potential, showing constrained copy-number levels over multiple linked genes. Thus, dicentric chromosomes may be an important precipitant of chromothripsis, as we show rob(15;21)c to be constitutionally dicentric and breakage-fusion-bridge cycles generate dicentric chromosomes somatically. Furthermore, our data illustrate that several cancer-specific mutational processes, applied sequentially, can coordinate to fashion copy-number profiles over large genomic scales, incrementally refining the fitness benefits of aggregated gene dosage changes.

The genetic heterogeneity and mutational burden of engineered melanomas in zebrafish models.

BACKGROUND: Melanoma is the most deadly form of skin cancer. Expression of oncogenic BRAF or NRAS, which are frequently mutated in human melanomas, promote the formation of nevi but are not sufficient for tumorigenesis. Even with germline mutated p53, these engineered melanomas present with variable onset and pathology, implicating additional somatic mutations in a multi-hit tumorigenic process. RESULTS: To decipher the genetics of these melanomas, we sequence the protein coding exons of 53 primary melanomas generated from several BRAF(V600E) or NRAS(Q61K) driven transgenic zebrafish lines. We find that engineered zebrafish melanomas show an overall low mutation burden, which has a strong, inverse association with the number of initiating germline drivers. Although tumors reveal distinct mutation spectrums, they show mostly C > T transitions without UV light exposure, and enrichment of mutations in melanogenesis, p53 and MAPK signaling. Importantly, a recurrent amplification occurring with pre-configured drivers BRAF(V600E) and p53-/- suggests a novel path of BRAF cooperativity through the protein kinase A pathway. CONCLUSION: This is the first analysis of a melanoma mutational landscape in the absence of UV light, where tumors manifest with remarkably low mutation burden and high heterogeneity. Genotype specific amplification of protein kinase A in cooperation with BRAF and p53 mutation suggests the involvement of melanogenesis in these tumors. This work is important for defining the spectrum of events in BRAF or NRAS driven melanoma in the absence of UV light, and for informed exploitation of models such as transgenic zebrafish to better understand mechanisms leading to human melanoma formation.

Genome-wide copy number profiling of single cells in S-phase reveals DNA-replication domains.

Single-cell genomics is revolutionizing basic genome research and clinical genetic diagnosis. However, none of the current research or clinical methods for single-cell analysis distinguishes between the analysis of a cell in G1-, S- or G2/M-phase of the cell cycle. Here, we demonstrate by means of array comparative genomic hybridization that charting the DNA copy number landscape of a cell in S-phase requires conceptually different approaches to that of a cell in G1- or G2/M-phase. Remarkably, despite single-cell whole-genome amplification artifacts, the log2 intensity ratios of single S-phase cells oscillate according to early and late replication domains, which in turn leads to the detection of significantly more DNA imbalances when compared with a cell in G1- or G2/M-phase. Although these DNA imbalances may, on the one hand, be falsely interpreted as genuine structural aberrations in the S-phase cell's copy number profile and hence lead to misdiagnosis, on the other hand, the ability to detect replication domains genome wide in one cell has important applications in DNA-replication research. Genome-wide cell-type-specific early and late replicating domains have been identified by analyses of DNA from populations of cells, but cell-to-cell differences in DNA replication may be important in genome stability, disease aetiology and various other cellular processes.

Single-cell copy number variation detection.

Detection of chromosomal aberrations from a single cell by array comparative genomic hybridization (single-cell array CGH), instead of from a population of cells, is an emerging technique. However, such detection is challenging because of the genome artifacts and the DNA amplification process inherent to the single cell approach. Current normalization algorithms result in inaccurate aberration detection for single-cell data. We propose a normalization method based on channel, genome composition and recurrent genome artifact corrections. We demonstrate that the proposed channel clone normalization significantly improves the copy number variation detection in both simulated and real single-cell array CGH data.

Full sib pens of pigs are not suitable to identify variance component of associative effect: a simulation study using Gibbs Sampling.

BACKGROUND: Accounting for and quantifying the associative effect of each animal could improve both welfare of animals and response to selection. Because of the limitation of REML, Gibbs Sampling could be an alternative technique to estimate the variance component of the associative effect. The objective of this study was to investigate the estimation accuracy of the variance component of associative effect by using simulation via Gibbs Sampling. The simulated data comprised five generations of pigs. The breeding animals of each generation were selected randomly. In the simulation, variations were introduced for the methods of assigning pens (random, mixed sib and full sib), the number of pigs per pen (5 or 10), the number of breeding animals per generation (162 or 324) and the correlation between genetic direct effect and genetic associative effect (-0.5, 0.1 or +0.5). Each set of simulation was run for 30 replications. RESULTS: Random assignment or mixed sib assignment resulted in bias of estimated variance components in only 3 of 24 combinations. Furthermore, these 3 cases occurred with 162 breeding animals. With full sib assignment, 9 out of 12 groups of estimates significantly deviated from the true parameter value. The Root Mean Square Errors obtained with the full sib assignment were higher than with the other two methods of pen assignment in most of the cases. The Root Mean Square Errors obtained with datasets with 324 breeding animals were notably smaller than the datasets from 162 breeding animals. Within each method of pen assignment, the relative bias of the associative effect was significantly smaller with group size 10 than with group size 5. CONCLUSION: Full sib assignment caused difficulties to estimate variance components in most of the cases, due to a lack of identifiability. With random and mixed assignment, most data structures yielded unbiased results but increasing the number of breeding animals or group size improves the estimation. Thus to get identifiable and unbiased estimates of the genetic associative effect, it is recommended to avoid close genetic relationship between animals within one pen and to use sufficient numbers of breeding animals and sufficient group sizes.

Axillary lymph node status of operable breast cancers by combined steroid receptor and HER-2 status: triple positive tumours are more likely lymph node positive.

AIMS: To examine the frequency of axillary lymph node (ALN) invasion of operable breast cancers by their combined oestrogen receptor (ER), progesterone receptor (PR) and HER-2 status. METHODS: 2227 recently operated cases in one centre were retrieved from the Multidisciplinary Breast Centre database and stratified according to their combined immunohistochemical (IHC) expression of ER/PR/HER-2 status. An equivocal HER-2 status was further analysed by Fluorescence in situ Hybridisation (FISH). The following 6 groups were considered: ER(-)PR(-)HER-2(-) (NNN; triple negative), ER(-)PR(-)HER-2(+) (NNP), ER(+)PR(-)HER-2(-) (PNN), ER(+)PR(-)HER-2(+) (PNP), ER(+)PR(+)HER-2(- )(PPN), ER(+)PR(+)HER-2(+) (PPP; triple positive). For ALN, the following variables were tested in uni- and multivariate models: age at diagnosis (years), tumour size (mm), tumour grade, ER, PR, HER-2 and the combined steroid receptor and HER-2 status. Likelihood ratio chi(2)-tests were used for univariate analysis and logistic regression for multivariate analysis. RESULTS: Triple positive tumours had a higher likelihood of being ALN positive than others (56.2% versus 35.7%; P<0.0001). Univariate logistic regression also withheld age, size, grade and HER-2 as predictors of ALN involvement. Final multivariate logistic regression revealed age, size, grade and PPP versus non-PPP to be independent predictors of ALN involvement; the odds ratio (OR) and 95% CI for PPP versus non-PPP tumours was 2.169 (1.490-3.156). CONCLUSION: Our data provide insight into the natural history of triple positive breast carcinomas. Such tumours are more likely ALN positive than those with another steroid receptor and HER-2 status. How these findings correlate with breast cancer prognosis remains to be investigated.