Genome-wide sequencing, mRNA and miRNA expression, DNA copy number and methylation analyses were performed on 117 Wilms tumors, followed by targeted sequencing of 651 Wilms tumors. In addition to genes previously implicated in Wilms tumors (WT1, CTNNB1, FAM123B, DROSHA, DGCR8, XPO5, DICER1, SIX1, SIX2, MLLT1, MYCN, and TP53), mutations were identified in genes not previously recognized as recurrently involved in Wilms tumors, the most frequent being BCOR, BCORL1, NONO, MAX, COL6A3, ASXL1, MAP3K4, and ARID1A. DNA copy number changes resulted in recurrent 1q gain, MYCN amplification, LIN28B gain, and let-7a loss. Unexpected germline variants involved PALB2 and CHEK2. Integrated analyses support two major classes of genetic changes that preserve the progenitor state and/or interrupt normal development.

SUMMARY We report the most common single nucleotide substitution/deletion mutations in Favorable Histology Wilms Tumors (FHWT) to occur within SIX1/2 (7% of 534 tumors) and microRNA processing genes (miRNAPG) DGCR8 and DROSHA (15% of 534 tumors). Comprehensive analysis of 77 FHWTs indicates that tumors with SIX1/2 and/or miRNAPG mutations show a pre-induction metanephric mesenchyme gene expression pattern and are significantly associated with both perilobar nephrogenic rests and 11p15 imprinting aberrations. Significantly decreased expression of mature Let-7a and the miR-200 family (responsible for mesenchymal-to-epithelial transition) in miRNAPG-mutant tumors is associated with an undifferentiated blastemal histology. The combination of SIX and miRNAPG mutations in the same tumor is associated with evidence of RAS activation and a higher rate of relapse and death.

Wilms tumors (WT) have provided broad insights into the interface between development and tumorigenesis. Further understanding is confounded by their genetic, histologic, and clinical heterogeneity, the basis of which remains largely unknown. We evaluated 224 WT for global gene expression patterns; WT1, CTNNB1, and WTX mutation; and 11p15 copy number and methylation patterns. Five subsets were identified showing distinct differences in their pathologic and clinical features: these findings were validated in 100 additional WT. The gene expression pattern of each subset was compared with published gene expression profiles during normal renal development. A novel subset of epithelial WT in infants lacked WT1, CTNNB1, and WTX mutations and nephrogenic rests and displayed a gene expression pattern of the postinduction nephron, and none recurred. Three subsets were characterized by a low expression of WT1 and intralobar nephrogenic rests. These differed in their frequency of WT1 and CTNNB1 mutations, in their age, in their relapse rate, and in their expression similarities with the intermediate mesoderm versus the metanephric mesenchyme. The largest subset was characterized by biallelic methylation of the imprint control region 1, a gene expression profile of the metanephric mesenchyme, and both interlunar and perilobar nephrogenic rests. These data provide a biologic explanation for the clinical and pathologic heterogeneity seen within WT and enable the future development of subset-specific therapeutic strategies. Further, these data support a revision of the current model of WT ontogeny, which allows for an interplay between the type of initiating event and the developmental stage in which it occurs.