Tumor DNA Methylation Profiling Is a Robust and Independent Prognostic Marker for Pediatric Patients With Adrenocortical Tumors

Pediatric adrenocortical tumors (pACTs) display complex genomic background and lack robust prognostic biomarkers. However, pACT methylation profiling was recently associated with patient prognosis. In order to evaluate whether tumor DNA methylation is associated with patient prognosis, we studied a Brazilian cohort of pACT in which most of the patients were carriers of the germline P53 p.R337H mutation (86%). We analyzed the DNA methylation profile of 57 tumor samples (MethylationEPIC BeadChip Array, Illumina) and the respective patients’ clinicopathological features associated with outcome and overall survival. Median age at diagnosis of the 40 girls and 17 boys was 2.1 years (range: 0.2–16.6). Reduced overall survival was observed in patients diagnosed after 4 years of age (n=43, HR=26.5, p<0.0001, 5-year survival (5ys)=43%), in those presenting with Cushing’s syndrome (n=7, HR=12.5, p<0.0001, 5ys=28%), advanced/non localized disease (IPACTR stages III and IV; n=11, HR=12.9, p<0.0001, 5ys=33%) or carcinomas (n=10, HR=4.9, p=0.04 5ys=56%). Unsupervised analysis of methylation data (Euclidean distance and Ward’s method; Pvclust package, R, version 4.0.3) included all 766,031 probes filtered in by the quality control step, and revealed two groups with distinct methylation signatures with a bootstrap support of 99%, namely pACT1 and pACT2. The overall variance of differentially methylated probes (DMPs; FDR p-values <10–6 and, simultaneously, median M-value difference >2) between the groups was 0.084% (n=642 probes). Tumors from pACT1 were enriched for hypomethylation in genes’ body region and for hypermethylation in promoter regions (p<0.0001). Gene set analysis showed that DMPs were mostly enriched in gene sets involved in nervous system differentiation and development, mRNA binding, and peptide secretion and related to pathways in cancer, cell adhesion, and calcium signaling. Methylation was not associated with germline P53 p.R337H mutation (p=0.57). Compared to pACT2, pACT1 group comprised patients older at diagnosis (p=0.0004) and higher frequencies of patients diagnosed after 4 years of age (80 vs 13%; p<0.0001), presenting with Cushing syndrome (40 vs. 4%; p=0.01), non-localized/advanced disease (70 vs 8%; p=0.0001), carcinomas (71 vs 17%. p=0.009), who needed adjuvant chemotherapy (90 vs. 9%; p<0.0001), and who experienced post-surgery recurrence or metastasis (90 vs 6%; p<0.0001). pACT1 methylation signature was associated with reduced overall survival from the disease (HR=51, p=0.0001, 5ys: 20 vs 98%), and remained associated after adjusting for prognostic features. Interestingly, different tumor methylation signatures were associated with pACT clinical presentation and evolution, regardless the presence of the P53 p.R337H mutation. Our data supports the analysis of methylation profiling as a robust and independent prognostic biomarker for pACT.

HIF1A is Overexpressed in Medulloblastoma and its Inhibition Reduces Proliferation and Increases EPAS1 and ATG16L1 Methylation

BACKGROUND

Genetic and epigenetic modifications are closely related to tumor initiation and progression and can provide guidance for understanding tumor functioning, potentially leading to the discovery of new therapies. Studies have associated hypoxia-related genes to tumor progression and chemo/radioresistance in brain tumors. Information on the expression profile of hypoxiarelated genes in pediatric medulloblastoma, although scarce, may reveal relevant information that could support treatment decisions.

OBJECTIVE

Our study focused on evaluation the of CA9, CA12, HIF1A, EPAS1, SCL2A1 and VEGF genes in 41 pediatric fresh-frozen medulloblastoma sample. Additionally, we analyzed the effect of hypoxia and normoxia in the pediatric medulloblastoma cell-line UW402. Furthermore, we assessed the effects of HIF1A knockdown in cell-proliferation and methylation levels of genes related to hypoxia, apoptosis and autophagy.

METHOD

qPCR was performed to evaluate mRNA levels, and Western blot to confirm HIF1A silencing in both patient samples and cell line. Pyrosequencing was performed to asses the methylation levels after HIF1A knockdown in the UW402 cell line.

RESULTS

A higher HIF1A mRNA level was observed in MB patients when compared to the cerebellum (non-tumor match). In UW402 MB cell-line, chemically induced hypoxic resulted in an increase of mRNA levels of HIF1A, VEGF, SCL2A1 and CA9 genes. Additionally, HIF1A knockdown induced a decrease in the expression of hypoxia related genes and a decrease of 30% in cell proliferation was also observed. Also, a significant increase in the methylation of ATG16L1 promoter and decrease in the methylation of EPAS1 promoter were observed after HIF1A knockdown.

CONCLUSION

HIF1A knockdown in medulloblastoma cells lead to decreased cellular proliferation, suggesting that HIF1A can be a potential therapeutic target to be explored in the medulloblastoma. However, the mechanisms behind HIF1A protein stabilization and function are very complex and more data need to be generated to potentially use HIF1A as a therapeutical target.

Downregulation of 14q32 microRNAs in Primary Human Desmoplastic Medulloblastoma

Medulloblastoma (MB) is one of the most common pediatric cancers, likely originating from abnormal development of cerebellar progenitor neurons. MicroRNA (miRNA) has been shown to play an important role in the development of the central nervous system. Microarray analysis was used to investigate miRNA expression in desmoplastic MB from patients diagnosed at a young age (1 or 2 years old). Normal fetal or newborn cerebellum was used as control. A total of 84 differentially expressed miRNAs (64 downregulated and 20 upregulated) were found. Most downregulated miRNAs (32/64) were found to belong to the cluster of miRNAs at the 14q32 locus, suggesting that this miRNA locus is regulated as a module in MB. Possible mechanisms of 14q32 miRNAs downregulation were investigated by the analysis of publicly available gene expression data sets. First, expression of estrogen-related receptor-γ (ESRRG), a reported positive transcriptional regulator of some 14q32 miRNAs, was found downregulated in desmoplastic MB. Second, expression of the parentally imprinted gene MEG3 was lower in MB in comparison to normal cerebellum, suggesting a possible epigenetic silencing of the 14q32 locus. miR-129-5p (11p11.2/7q32.1), miR-206 (6p12.2), and miR-323-3p (14q32.2), were chosen for functional studies in DAOY cells. Overexpression of miR-129-5p using mimics decreased DAOY proliferation. No effect was found with miR-206 or miR-323 mimics.

Assessment of BCL2/J(H) translocation in healthy individuals exposed to low-level radiation of 137CsCl in Goiânia, Goiás, Brazil

Healthy radio-exposed individuals who received low levels of Cesium-137 radiation during the accident that occurred in Goiânia in 1987, their families and controls were tested for the detection of t(14;18)-rearranged B cells in peripheral blood by using a highly sensitive, real-time quantitative PCR method. The chromosomal translocation t(14;18)(q32;q21) is characteristic of follicular lymphoma and is a frequent abnormality observed in other types of non-Hodgkin's lymphoma. This translocation leads to constitutive activation of the BCL2 oncogene by the enhancers of the immunoglobulin heavy-chain locus. In healthy individuals, the same translocation may also be found in a small fraction of peripheral blood lymphocytes, and positive cells might serve as an indicator for environmental exposure to carcinogens and possibly correlate with the cumulative risk of developing t(14;18)- positive non-Hodgkin's lymphoma. Twenty healthy radio-exposed individuals, 10 relatives and 10 non-exposed healthy individuals were tested for the detection of this translocation. Only 1 non-exposed individual was positive for the chromosomal translocation, and healthy radio-exposed individuals presented lower levels of cells bearing the BCL2/J(H) rearrangement when compared to the levels of the patients with follicular lymphoma before treatment. However, evaluation of more cells would be required to confirm the total absence of circulating cells bearing BCL2/J(H) rearrangement.

TAL1/SCL is downregulated upon histone deacetylase inhibition in T-cell acute lymphoblastic leukemia cells

The transcription factor T-cell acute lymphocytic leukemia (TAL)-1 is a major T-cell oncogene associated with poor prognosis in T-cell acute lymphoblastic leukemia (T-ALL). TAL1 binds histone deacetylase 1 and incubation with histone deacetylase inhibitors (HDACis) promotes apoptosis of leukemia cells obtained from TAL1 transgenic mice. Here, we show for the first time that TAL1 protein expression is strikingly downregulated upon histone deacetylase inhibition in T-ALL cells. This is due to decreased TAL1 gene transcription in cells with native TAL1 promoter, and due to impaired TAL1 mRNA translation in cells that harbor the TAL1(d) microdeletion and consequently express TAL1 under the control of the SCL/TAL1 interrupting locus (SIL) promoter. Notably, HDACi-triggered apoptosis of T-ALL cells is significantly reversed by TAL1 forced overexpression. Our results indicate that the HDACi-mediated apoptotic program in T-ALL cells is partially dependent on their capacity to downregulate TAL1 and provide support for the therapeutic use of HDACi in T-ALL.

Conservation of glycolytic oscillations in Saccharomyces cerevisiae and human pancreatic beta-cells: a study of metabolic robustness

The present study compares two computer models of the first part of glucose catabolism in different organisms in search of evolutionarily conserved characteristics of the glycolysis cycle and proposes the main parameters that define the stable steady-state or oscillatory behavior of the glycolytic system. It is suggested that in both human pancreatic beta-cells and Saccharomyces cerevisiae there are oscillations that, despite differences in wave form and period of oscillation, share the same robustness strategy: the oscillation is not controlled by only one but by at least two parameters that will have more or less control over the pathway flux depending on the initial state of the system as well as on extra-cellular conditions. This observation leads to two important interpretations: the first is that in both S. cerevisiae and human beta-cells, despite differences in enzyme kinetics and mechanism of feedback control, evolution seems to have kept an oscillatory behavior coupled to the glucose concentration outside the cytoplasm, and the second is that the development of drugs to regulate metabolic dysfunctions in more complex systems may require further study, not only determining which enzyme is controlling the flux of the system but also under which conditions and how its control is maintained by the enzyme or transferred to other enzymes in the pathway as the drug starts acting.

Rapid, reliable and inexpensive quality assessment of biotinylated cRNA

The interpretation of oligonucleotide array experiments depends on the quality of the target cRNA used. cRNA target quality is assessed by quantitative analysis of the representation of 5' and 3' sequences of control genes using commercially available Test arrays. The Test array provides an economically priced means of determining the quality of labeled target prior to analysis on whole genome expression arrays. This manuscript validates the use of a duplex RT-PCR assay as a faster (6 h) and less expensive (<10 US dollars) but equally accurate alternative to the Test arrays in determining biotinylated cRNA quality. Forty-one different cRNA samples were hybridized to HG-U133A microarrays from Affymetrix. Ten cRNA samples with a beta-actin 3'/5' ratio >6 were chosen and classified as degraded cRNAs, and 31 samples with a beta-actin 3'/5' ratio <6 were selected as good quality cRNAs. Blinded samples were then used for the RT-PCR assay. After gel electrophoresis, optical densities of the amplified 3' and 5' fragments of beta-actin were measured and the 3'/5' ratio was calculated. There was a strong correlation (r(2) = 0.6802) between the array and the RT-PCR beta-actin 3'/5' ratios. Moreover, the RT-PCR 3'/5' ratio was significantly different (P < 0.0001) between undegraded (mean +/- SD, 0.34 +/- 0.09) and degraded (1.71 +/- 0.83) samples. None of the other parameters analyzed, such as i) the starting amount of RNA, ii) RNA quality assessed using the Bioanalyzer Chip technology, or iii) the concentration and OD260/OD280 ratio of the purified biotinylated cRNA, correlated with cRNA quality.

Detection of oncogenic human papillomavirus in sporadic retinoblastoma

PURPOSE

To investigate the presence of human papillomavirus (HPV) DNA in tumour tissue from patients with unilateral retinoblastoma.

METHODS

Samples of paraffin-embedded tumour tissue from 43 children with unilateral retinoblastoma were collected to investigate the presence of HPV DNA using polymerase chain reaction (PCR) and dot blot hybridization.

RESULTS

Oncogenic HPV DNA types 16 and 35 were detected in 12 (27.9%) of 43 tumour specimens. A higher frequency of differentiated tumours (63.3%) was observed among the HPV-positive tumours.

CONCLUSIONS

Future studies are necessary to demonstrate an association between HPV and sporadic retinoblastoma.

The DapA gene encoding the lysine biosynthetic enzyme dihydrodipicolinate synthase from Coix lacryma-jobi: cloning, characterization, and expression analysis

Dihydrodipicolinate synthase (DHPS) is the main enzyme of a specific branch of the aspartate pathway leading to lysine biosynthesis in higher plants. We have cloned and characterized the DHPS-encoding Dap)A gene from the maize-related grass Coix lacryiana-jobi. The DapA open reading frame is interrupted by two introns and encodes the 326 amino acid-long Coix DHPS protein, which is 95% identical to the maize DHPS protein. Coix DNA gel blot analysis with maize DHPS cDNA as a probe showed a single strongly hybridizing band along with faint bands. RNA gel blot analysis showed that DHPS transcripts are present in coleoptiles, embryos, endosperms, and roots but are almost undetectable in blades of young leaves of both Coix and maize. The 5'-flanking region of the DapA gene contains a TGACTC GCN4-like element located 372 bp upstream the putative translation start codon. Steady-state levels of DHPS mRNA were slightly reduced in the endosperms and embryos of the maize lysine-rich opaque2 mutants when compared with those in normal kernels. Selective binding assay with the maize Opaque2 protein (O2) showed that the GCN4-like element is not an O2 binding site, suggesting that the DHPS gene is not under the control of O2.

Identification of a DNA-binding factor that recognizes an alpha-coixin promoter and interacts with a Coix Opaque-2 like protein

Transient expression and electrophoretic mobility shift assay were used to investigate the cis elements and the DNA-binding proteins involved in the regulation of expression of a 22 kDa zein-like alpha-coixin gene. A set of unidirectional deletions was generated in a 962 bp fragment of the alpha-coixin promoter that had been previously fused to the reporter gene GUS. The constructs were assayed by transient expression in immature maize endosperm. There was no significant decrease in GUS activity as deletions progressed from -1084 to -238. However, deletion from -238 to -158, which partially deleted the O2c box, resulted in a dramatic decrease in GUS activity emphasizing the importance of the O2 box in the quantitative expression of the gene. The -238 promoter fragment interacted with Coix endosperm nuclear proteins to form 5 DNA-protein complexes, C1-C5, as detected by EMSA. The same retarded complexes were observed when the -158 promoter fragment was used in the binding reactions. Reactions with nuclear extracts isolated from Coix endosperms harvested from 6 to 35 days after pollination revealed that the 5 DNA-protein complexes that interact with the alpha-coixin promoter are differentially assembled during seed development. Deletion analysis carried out on the -238/ATG promoter fragment showed that a 35 bp region from -86 to -51 is essential for the formation of the complexes observed. When nuclear extracts were incubated with an antiserum raised against the maize Opaque-2 protein, the formation of 4 complexes, C1, C3, C4 and C5, was prevented indicating that an Opaque-2 like protein participates in the formation of those complexes. Complex C2 was not affected by the addition of the O2 antibody, suggesting the existence of a novel nuclear factor, CBF1, that binds to the promoter and makes protein-protein associations with other proteins present in Coix endosperm nuclei.

Cooperative DNA binding and sequence discrimination by the Opaque2 bZIP factor

The maize Opaque2 (O2) protein is a basic leucine zipper transcription factor that controls the expression of distinct classes of endosperm genes through the recognition of different cis-acting elements in their promoters. The O2 target region in the promoter of the alpha-coixin gene was analyzed in detail and shown to comprise two closely adjacent binding sites, named O2u and O2d, which are related in sequence to the GCN4 binding site. Quantitative DNase footprint analysis indicated that O2 binding to alpha-coixin target sites is best described by a cooperative model. Transient expression assays showed that the two adjacent sites act synergistically. This synergy is mediated in part by cooperative DNA binding. In tobacco protoplasts, O2 binding at the O2u site is more important for enhancer activity than is binding at the O2d site, suggesting that the architecture of the O2-DNA complex is important for interaction with the transcriptional machinery.

The molecular and functional characterization of an Opaque2 homologue gene from Coix and a new classification of plant bZIP proteins

The seed storage proteins of Coix, sorghum and maize are codified by homologous genes which are coordinately expressed in the endosperm in a temporal-specific fashion. Opaque2 (O2), a bZIP protein originally isolated from maize, has been described as a transcription activator of alpha- and beta-prolamin genes. The isolation and characterization of cDNA and genomic clones encoding the Opaque2 homologue from Coix are reported here. The coding region of the Coix O2 gene is interrupted by five introns and codifies a polypeptide of 408 amino acids. Comparison of the deduced amino acid sequence with two different sequences of maize O2 protein showed that the Coix O2 protein is similar to the maize O2 isolated from W22 maize inbred line. The Coix O2 protein has the same binding specificity and expression pattern of the maize O2. The O2 proteins together with OHP1, OsBZIPPA, SPA, CPRF2 and RITA1 were assigned to one of the five bZIP plant families in an updated classification of plant bZIP according to bZIP domain similarity.

The involvement of Opaque 2 on beta-prolamin gene regulation in maize and Coix suggests a more general role for this transcriptional activator

The maize opaque 2 (o2) mutation is known to have numerous pleiotropic effects. Some polypeptides have their expression depressed while others are enhanced. The best characterized effects of the o2 mutation are those exerted on endosperm genes encoding the storage protein class of the 22 kDa alpha-zeins and the ribosome inactivating protein b-32. The Opaque 2 (O2) locus encodes a basic domain-leucine zipper DNA-binding factor, O2, which transcriptionally regulates these genes. In the maize-related grass Coix lacryma-jobi, an O2-homologous protein regulates the 25 kDa alpha-coixin family. We show in this paper that O2 transcriptionally regulates the structurally and developmentally different class of the beta-prolamins. A new O2-binding box was identified in beta-prolamin genes from maize and Coix that, together with the boxes previously identified in other endosperm expressed genes, forms a curious collection of O2 cis elements. This may have regulatory implications on the role of O2 in the mechanism that controls coordinated gene expression in the developing endosperm. Considering that the O2 locus controls at least three distinct classes of genes in maize endosperm, we propose that the O2 protein may play a more general role in maize endosperm development than previously conceived.

Structural characterization and promoter activity analysis of the gamma-kafirin gene from sorghum

A genomic clone encoding the gamma-kafirin gene from sorghum was isolated and sequenced. A 2938 bp sequenced fragment includes an intronless open reading frame of 636 nucleotides encoding a putative polypeptide of 212 amino acids. Comparison of the deduced amino acid sequence of gamma-kafirin with the published sequences of gamma-prolamins of maize, and Coix revealed highly conserved domains. The N-terminal region of these proteins contains the conserved hexapeptide PPPVHL, which is repeated eight times in gamma-zein, four times in gamma-kafirin and three times in gamma-coixin. The number of PPPVHL repeats accounts predominantly for the differences in the molecular weights of gamma-prolamins. Several putative regulatory sequences common to the gamma-kafirin and gamma-zein genes were identified in both the 5' and the 3' flanking regions. Putative GCN4-like regulatory sequences were found at positions -192 and -476 in the 5' flanking region of gamma-kafirin. In the 3' noncoding region, three putative polyadenylation signals, two AATAAT and one AATGAA, were found at positions +658, +716, and +785, respectively. In order to investigate the role of the putative GCN4-like motifs and other possible cis-acting element(s) of the gamma-kafirin promoter, a series of deleted and chimeric promoter constructs were introduced into maize, Coix and sorghum tissues by particle bombardment. Histochemical analysis of beta-glucuronidase (GUS) activity in different tissues indicated that the element(s) responsible for tissue specificity is probably located in the 285-bp proximal region of the promoter, while the remaining promoter sequence seems to carry the element(s) responsible for the quantitative response.

The transcriptional activator Opaque2 recognizes two different target sequences in the 22-kD-like alpha-prolamin genes

The maize Opaque2 (O2) protein is a "leucine zipper" DNA binding factor that interacts with the sequence TCCACGTAGA in the promoters of the 22-kD alpha-zein genes and activates its transcription. A completely different consensus sequence (GATGAPyPuTGPu) identified in b-32, a gene that encodes an abundant albumin that is also under control of the O2 locus, can also be bound by the O2 protein. We showed that the gene encoding the 22-kD-like alpha-coixin, the alpha-prolamin of the maize-related grass Coix, can also be transactivated by the O2 protein. A binding assay in vitro and footprint analysis demonstrated that the GACATGTC sequence of the alpha-coixin promoter can be recognized and protected by the maize O2 protein. Employing transient expression experiments in immature maize endosperm and tobacco mesophyll protoplasts, we demonstrated that the O2 protein can activate expression of the beta-glucuronidase reporter gene placed under the control of the 22-kD-like alpha-coixin promoter. We also demonstrated that a 22-kD-like alpha-coixin pseudogene promoter is transactivated by the maize O2 protein.

Sequence analysis of 22 kDa-like alpha-coixin genes and their comparison with homologous zein and kafirin genes reveals highly conserved protein structure and regulatory elements

Several genomic and cDNA clones encoding the 22 kDa-like alpha-coixin, the alpha-prolamin of Coix seeds, were isolated and sequenced. Three contiguous 22 kDa-like alpha-coixin genes designated alpha-3A, alpha-3B and alpha-3C were found in the 15 kb alpha-3 genomic clone. The alpha-3A and alpha-3C genes presented in-frame stop codons at position +652. The two genes with truncated ORFs are flanking the alpha-3B gene, suggesting that the three alpha-coixin genes may have arisen by tandem duplication and that the stop codon was introduced before the duplication. Comparison of the deduced amino acid sequences of alpha-coixin clones with the published sequences of 22 kDa alpha-zein and 22 kDa-like alpha-kafirin revealed a highly conserved protein structure. The protein consists of an N-terminus, containing the signal peptide, followed by ten highly conserved tandem repeats of 15-20 amino acids flanked by polyglutamines, and a short C-terminus. The difference between the 22 kDa-like alpha-prolamins and the 19 kDa alpha-zein lies in the fact that the 19 kDa protein is exactly one repeat motif shorter than the 22 kDa proteins. Several putative regulatory sequences common to the zein and kafirin genes were identified within both the 5' and 3' flanking regions of alpha-3B. Nucleotide sequences that match the consensus TATA, CATC and the ca. -300 prolamin box are present at conserved positions in alpha-3B relative to zein and kafirin genes. Two putative Opaque-2 boxes are present in alpha-3B that occupies approximately the same positions as those identified for the 22 kDa alpha-zein and alpha-kafirin genes. Southern hybridization, using a fragment of a maize Opaque-2 cDNA clone as a probe, confirmed the presence of Opaque-2 homologous sequences in the Coix and sorghum genomes. The overall results suggest that the structural and regulatory genes involved in the expression of the 22 kDa-like alpha-prolamin genes of Coix, sorghum and maize, originated from a common ancestor, and that variations were introduced in the structural and regulatory sequences after species separation.