Methylome analysis identifies a Wilms tumor epigenetic biomarker detectable in blood

Background

Wilms tumor is the most common pediatric renal malignancy and there is a clinical need for a molecular biomarker to assess treatment response and predict relapse. The known mutated genes in this tumor type show low mutation frequencies, whereas aberrant methylation at 11p15 is by far the most common aberration. We therefore analyzed the epigenome, rather than the genome, to identify ubiquitous tumor-specific biomarkers.

Results

Methylome analysis of matched normal kidney and Wilms tumor identifies 309 preliminary methylation variable positions which we translate into three differentially methylated regions (DMRs) for use as tumor-specific biomarkers. Using two novel algorithms we show that these three DMRs are not confounded by cell type composition. We further show that these DMRs are not methylated in embryonic blastema but are intermediately methylated in Wilms tumor precursor lesions. We validate the biomarker DMRs using two independent sample sets of normal kidney and Wilms tumor and seven Wilms tumor histological subtypes, achieving 100% and 98% correct classification, respectively. As proof-of-principle for clinical utility, we successfully use biomarker DMR-2 in a pilot analysis of cell-free circulating DNA to monitor tumor response during treatment in ten patients.

Conclusions

These findings define the most common methylated regions in Wilms tumor known to date which are not associated with their embryonic origin or precursor stage. We show that this tumor-specific methylated DNA is released into the blood circulation where it can be detected non-invasively showing potential for clinical utility.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0434-y) contains supplementary material, which is available to authorized users.     

Methylation Imprinting of H19 and SNRPN Genes in Human Benign Ovarian Teratomas

In humans, studies of female germ cells are very limited by ethics. The current study investigated the usefulness of benign ovarian teratomas as a substitute for ova in analyses of imprinted genes. Twenty-five human benign ovarian teratomas were typed with 45 microsatellite DNA markers and classified according to their genotypic features. Two oppositely imprinted genes, H19 and SNRPN, were then chosen for analysis of their methylation states in these tumors. These analyses revealed that benign ovarian teratomas consist of a mixture of genetically and epigenetically heterogeneous cell populations. In contrast to previous reports, we could document only one case rising from germ cells by meiosis-II nondisjunction. H19 and SNRPN were methylated in individual teratomas to various degrees, ranging from normal somatic cell to expected ovum levels. The allele with residual methylation of H19 was consistent with that methylated in the patient's blood DNA, thus being of paternal origin. Degrees of H19 hypomethylation and SNRPN hypermethylation increased as the cellular origin of the tumors advanced in oogenesis and were closely correlated in individual teratomas. These results could be best explained by the assumption that the primary imprinting is a progressively organized process and suggest that the establishment of primary imprints on different genes might be mechanistically linked, even when those genes are oppositely imprinted.     

Oncofetal H19 RNA promotes tumor metastasis

The oncofetal H19 gene transcribes a long non-coding RNA(lncRNA) that is essential for tumor growth. Here we found that numerous established inducers of epithelial to mesenchymal transition(EMT) also induced H19/miR-675 expression. Both TGF-β and hypoxia concomitantly induced H19 and miR-675 with the induction of EMT markers. We identified the PI3K/AKT pathway mediating the inductions of Slug, H19 RNA and miR-675 in response to TGF-β treatment, while Slug induction depended on H19 RNA. In the EMT induced multidrug resistance model, H19 level was also induced. In a mouse breast cancer model, H19 expression was tightly correlated with metastatic potential. In patients, we detected high H19 expression in all common metastatic sites tested, regardless of tumor primary origin. H19 RNA suppressed the expression of E-cadherin protein. H19 up-regulated Slug expression concomitant with the suppression of E-cadherin protein through a mechanism that involved miR-675. Slug also up-regulated H19 expression and activated its promoter. Altogether, these results may support the existence of a positive feedback loop between Slug and H19/miR-675, that regulates E-cadherin expression. H19 RNA enhanced the invasive potential of cancer cells in vitro and enhanced tumor metastasis in vivo. Additionally, H19 knockdown attenuated the scattering and tumorigenic effects of HGF/SF. Our results present novel mechanistic insights into a critical role for H19 RNA in tumor progression and indicate a previously unknown link between H19/miR-675, Slug and E-cadherin in the regulation of cancer cell EMT programs.     

Methylation changes of H19 gene in sperms of X-irradiated mouse and maintenance in offspring

The nature of imprinting is just differential methylation of imprinted genes. Unlike the non-imprinted genes, the methylation pattern of imprinted genes established during the period of gametogenesis remains unchangeable after fertilization and during embryo development. It implies that gametogenesis is the key stage for methylation pattern of imprinted genes. The imprinting interfered by exogenous factors during this stage could be inherited to offspring and cause genetic effect. Now many studies have proved that ionizing irradiation could disturb DNA methylation. Here we choose BALB/c mice as a research model and X-ray as interfering source to further clarify it. We discovered that the whole-body irradiation of X-ray to male BALB/c mice could influence the methylation pattern of H19 gene in sperms, which resulted in some cytosines of partial CpG islands in the imprinting control region could not transform to methylated cytosines. Furthermore, by copulating the interfered male mice with normal female, we analyzed the promoter methylation pattern of H19 in offspring fetal liver and compared the same to the pattern of male parent in sperms. We found that the majority of methylation changes in offspring liver were related to the ones in their parent sperms. Our data proved that the changes of the H19 gene methylation pattern interfered by X-ray irradiation could be transmitted and maintained in the first-generation offspring.     

A tumor chromosome rearrangement further defines the 11p13 Wilms tumor locus.

A sporadic Wilms tumor, WT-21, with an (11;14)-(p13;q23) reciprocal translocation has been identified. The translocation is found in tumor cells, but not in the patients' circulating lymphocytes. Molecular analysis of somatic cell hybrids segregating the derivative translocation chromosomes reveals a submicroscopic interstitial deletion at the translocation breakpoint, as well as a cytologically undetectable interstitial deletion in the nontranslocation chromosome 11, resulting in a homozygous deletion in 11p13. Pulsed-field gel analysis of tumor DNA indicates that the two deletions are indistinguishable, and the homozygously deleted region is less than 875 kb. The homozygously deleted regions of three other sporadic Wilms tumors overlap with the deleted region in WT-21, and the candidate cDNA clone for the 11p13 Wilms tumor gene described by Call et al. (Cell 60, 509-520, 1990) is included in the deleted region. These findings strengthen previous conclusions regarding the obligate location for the 11p13 WT locus and support the suggestion that the Wilms tumor gene has been cloned.     

Smallest region of overlap in Wilms tumor deletions uniquely implicates an 11p13 zinc finger gene as the disease locus.

The development of Wilms tumor (WT) has been associated with the inactivation of a "tumor suppressor" locus in human chromosome 11 band p13. Several WTs that exhibit homozygous deletions of an 11p13 candidate WT gene in its entirety have been reported. We report here a partial deletion of the candidate gene which, upon comparison with other documented homozygous deletions, permitted a precise definition of the critical genomic target in Wilms tumor. The smallest region of overlap between these deletions is a 16-kb segment of DNA encompassing the 5' exon(s) of an 11p13 gene coding for a zinc finger protein, together with an associated CpG island. This finding supports the notion that the candidate gene in question corresponds to the 11p13 WT1 Wilms tumor locus.     

Combined analysis of fourteen nuclear genes refines the Ursidae phylogeny

Despite numerous studies, questions remain about the evolutionary history of Ursidae and additional independent genetic markers were needed to elucidate these ambiguities. For this purpose, we sequenced ten nuclear genes for all the eight extant bear species. By combining these new sequences with those of four other recently published nuclear markers, we provide new insights into the phylogenetic relationships of the Ursidae family members. The hypothesis that the giant panda was the first species to diverge among ursids is definitively confirmed and the precise branching order within the Ursus genus is clarified for the first time. Moreover, our analyses indicate that the American and the Asiatic black bears do not cluster as sister taxa, as had been previously hypothesised. Sun and sloth bears clearly appear as the most basal ursine species but uncertainties about their exact relationships remain. Since our larger dataset did not enable us to clarify this last question, identifying rare genomic changes in bear genomes could be a promising solution for further studies.     

Regulation of CRABP-II expression by MycN in Wilms tumor

Cellular retinoic acid binding protein II (CRABP-II) is overexpressed in a wide variety of cancers. Previously we have shown that CRABP-II expression levels are also elevated in neuroblastoma and Wilms tumors. To elucidate the molecular mechanisms underlying the abnormal expression of CRABP-II in Wilms tumor, we studied the expression of MycN and CRABP-II in these tumor samples. Our data revealed that CRABP-II is overexpressed in Wilms tumor compared to normal adjacent non-neoplastic tissue and its levels are even higher in late stage tumors. Its expression correlates with MycN expression in tumors. The tumors that do not express MycN have no CRABP-II expression. The expression of CRABP-II is also regulated by methylation and its promoter is unmethylated in tumors. Knockdown of MycN by small interfering RNA leads to downregulation of CRABP-II. Thus our results suggest that both MycN and DNA methylation are responsible for CRABP-II expression in pediatric tumors and demethylation of CRABP-II may be an early event in tumor development.     

Assembly and analysis of cosmid contigs in the CEA-gene family region of human chromosome 19.

The carcinoembryonic antigen (CEA)-like genes are members of a large gene family which is part of the immunoglobulin superfamily. The CEA family is divided into two major subgroups, the CEA-subgroup and the pregnancy-specific glycoprotein (PSG)-subgroup. In the course of an effort to develop a set of overlapping cosmids spanning human chromosome 19, we identified 245 cosmids in a human chromosome 19 cosmid library (6-7X redundant) by hybridization with an IgC-like domain fragment of the CEA gene. A fluorescence-based restriction enzyme digest fingerprinting strategy was used to assemble 212 probe-positive cosmids, along with 115 additional cosmids from a collection of approximately 8,000 randomly selected cosmids, into five contigs. Two of the contigs contain CEA-subgroup genes while the remaining three contigs contain PSG-subgroup genes. These five contigs range in size from 100 kb to over 300 kb and span an estimated 1 Mb. The CEA-like gene family was determined by fluorescence in situ hybridization to map in the q13.1-q13.2 region of human chromosome 19. Analysis of the two CEA-subgroup contigs provided verification of the contig assembly strategy and insight into the organization of 9 CEA-subgroup genes.     

Systematic sequencing of the Escherichia coli genome: analysis of the 0-2.4 min region.

A contiguous 111,402-nucleotide sequence corresponding to the 0 to 2.4 min region of the E. coli chromosome was determined as a first step to complete structural analysis of the genome. The resulting sequence was used to predict open reading frames and to search for sequence similarity against the PIR protein database. A number of novel genes were found whose predicted protein sequences showed significant homology with known proteins from various organisms, including several clusters of genes similar to those involved in fatty acid metabolism in bacteria (e.g., betT, baiF) and higher organisms, iron transport (sfuA, B, C) in Serratia marcescens, and symbiotic nitrogen fixation or electron transport (fixA, B, C, X) in Azorhizobium caulinodans. In addition, several genes and IS elements that had been mapped but not sequenced (e.g., leuA, B, C, D) were identified. We estimate that about 90 genes are represented in this region of the chromosome with little spacer.