Genomic, transcriptomic and epigenomic sequencing data of the B-cell leukemia cell line REH.

OBJECTIVES: The aim of this data paper is to describe a collection of 33 genomic, transcriptomic and epigenomic sequencing datasets of the B-cell acute lymphoblastic leukemia (ALL) cell line REH. REH is one of the most frequently used cell lines for functional studies of pediatric ALL, and these data provide a multi-faceted characterization of its molecular features. The datasets described herein, generated with short- and long-read sequencing technologies, can both provide insights into the complex aberrant karyotype of REH, and be used as reference datasets for sequencing data quality assessment or for methods development. DATA DESCRIPTION: This paper describes 33 datasets corresponding to 867 gigabases of raw sequencing data generated from the REH cell line. These datasets include five different approaches for whole genome sequencing (WGS) on four sequencing platforms, two RNA sequencing (RNA-seq) techniques on two different sequencing platforms, DNA methylation sequencing, and single-cell ATAC-sequencing.

Functional analysis of the HD-Zip transcription factor genes Oshox12 and Oshox14 in rice.

The homeodomain-leucine zipper (HD-Zip) transcription factor family plays vital roles in plant development and morphogenesis as well as responses to biotic and abiotic stresses. In barley, a recessive mutation in Vrs1 (HvHox1) changes two-rowed barley to six-rowed barley, which improves yield considerably. The Vrs1 gene encodes an HD-Zip subfamily I transcription factor. Phylogenetic analysis has shown that the rice HD-Zip I genes Oshox12 and Oshox14 are the closest homologues of Vrs1. Here, we show that Oshox12 and Oshox14 are ubiquitously expressed with higher levels in developing panicles. Trans-activation assays in yeast and rice protoplasts demonstrated that Oshox12 and Oshox14 can bind to a specific DNA sequence, AH1 (CAAT(A/T)ATTG), and activate reporter gene expression. Overexpression of Oshox12 and Oshox14 in rice resulted in reduced panicle length and a dwarf phenotype. In addition, Oshox14 overexpression lines showed a deficiency in panicle exsertion. Our findings suggest that Oshox12 and Oshox14 may be involved in the regulation of panicle development. This study provides a significant advancement in understanding the functions of HD-Zip transcription factors in rice.

Deep targeted sequencing in pediatric acute lymphoblastic leukemia unveils distinct mutational patterns between genetic subtypes and novel relapse-associated genes.

To characterize the mutational patterns of acute lymphoblastic leukemia (ALL) we performed deep next generation sequencing of 872 cancer genes in 172 diagnostic and 24 relapse samples from 172 pediatric ALL patients. We found an overall greater mutational burden and more driver mutations in T-cell ALL (T-ALL) patients compared to B-cell precursor ALL (BCP-ALL) patients. In addition, the majority of the mutations in T-ALL had occurred in the original leukemic clone, while most of the mutations in BCP-ALL were subclonal. BCP-ALL patients carrying any of the recurrent translocations ETV6-RUNX1, BCR-ABL or TCF3-PBX1 harbored few mutations in driver genes compared to other BCP-ALL patients. Specifically in BCP-ALL, we identified ATRX as a novel putative driver gene and uncovered an association between somatic mutations in the Notch signaling pathway at ALL diagnosis and increased risk of relapse. Furthermore, we identified EP300, ARID1A and SH2B3 as relapse-associated genes. The genes highlighted in our study were frequently involved in epigenetic regulation, associated with germline susceptibility to ALL, and present in minor subclones at diagnosis that became dominant at relapse. We observed a high degree of clonal heterogeneity and evolution between diagnosis and relapse in both BCP-ALL and T-ALL, which could have implications for the treatment efficiency.

DNA methylation-based subtype prediction for pediatric acute lymphoblastic leukemia.

BACKGROUND: We present a method that utilizes DNA methylation profiling for prediction of the cytogenetic subtypes of acute lymphoblastic leukemia (ALL) cells from pediatric ALL patients. The primary aim of our study was to improve risk stratification of ALL patients into treatment groups using DNA methylation as a complement to current diagnostic methods. A secondary aim was to gain insight into the functional role of DNA methylation in ALL. RESULTS: We used the methylation status of ~450,000 CpG sites in 546 well-characterized patients with T-ALL or seven recurrent B-cell precursor ALL subtypes to design and validate sensitive and accurate DNA methylation classifiers. After repeated cross-validation, a final classifier was derived that consisted of only 246 CpG sites. The mean sensitivity and specificity of the classifier across the known subtypes was 0.90 and 0.99, respectively. We then used DNA methylation classification to screen for subtype membership of 210 patients with undefined karyotype (normal or no result) or non-recurrent cytogenetic aberrations ('other' subtype). Nearly half (n = 106) of the patients lacking cytogenetic subgrouping displayed highly similar methylation profiles as the patients in the known recurrent groups. We verified the subtype of 20% of the newly classified patients by examination of diagnostic karyotypes, array-based copy number analysis, and detection of fusion genes by quantitative polymerase chain reaction (PCR) and RNA-sequencing (RNA-seq). Using RNA-seq data from ALL patients where cytogenetic subtype and DNA methylation classification did not agree, we discovered several novel fusion genes involving ETV6, RUNX1, and PAX5. CONCLUSIONS: Our findings indicate that DNA methylation profiling contributes to the clarification of the heterogeneity in cytogenetically undefined ALL patient groups and could be implemented as a complementary method for diagnosis of ALL. The results of our study provide clues to the origin and development of leukemic transformation. The methylation status of the CpG sites constituting the classifiers also highlight relevant biological characteristics in otherwise unclassified ALL patients.

The mutational landscape in pediatric acute lymphoblastic leukemia deciphered by whole genome sequencing.

Genomic characterization of pediatric acute lymphoblastic leukemia (ALL) has identified distinct patterns of genes and pathways altered in patients with well-defined genetic aberrations. To extend the spectrum of known somatic variants in ALL, we performed whole genome and transcriptome sequencing of three B-cell precursor patients, of which one carried the t(12;21)ETV6-RUNX1 translocation and two lacked a known primary genetic aberration, and one T-ALL patient. We found that each patient had a unique genome, with a combination of well-known and previously undetected genomic aberrations. By targeted sequencing in 168 patients, we identified KMT2D and KIF1B as novel putative driver genes. We also identified a putative regulatory non-coding variant that coincided with overexpression of the growth factor MDK. Our results contribute to an increased understanding of the biological mechanisms that lead to ALL and suggest that regulatory variants may be more important for cancer development than recognized to date. The heterogeneity of the genetic aberrations in ALL renders whole genome sequencing particularly well suited for analysis of somatic variants in both research and diagnostic applications.

Interaction between the GROWTH-REGULATING FACTOR and KNOTTED1-LIKE HOMEOBOX families of transcription factors.

KNOTTED1-LIKE HOMEOBOX (KNOX) genes are important regulators of meristem function, and a complex network of transcription factors ensures tight control of their expression. Here, we show that members of the GROWTH-REGULATING FACTOR (GRF) family act as players in this network. A yeast (Saccharomyces cerevisiae) one-hybrid screen with the upstream sequence of the KNOX gene Oskn2 from rice (Oryza sativa) resulted in isolation of OsGRF3 and OsGRF10. Specific binding to a region in the untranslated leader sequence of Oskn2 was confirmed by yeast and in vitro binding assays. ProOskn2:ß-glucuronidase reporter expression was down-regulated by OsGRF3 and OsGRF10 in vivo, suggesting that these proteins function as transcriptional repressors. Likewise, we found that the GRF protein BGRF1 from barley (Hordeum vulgare) could act as a repressor on an intron sequence in the KNOX gene Hooded/Barley Knotted3 (Bkn3) and that AtGRF4, AtGRF5, and AtGRF6 from Arabidopsis (Arabidopsis thaliana) could repress KNOTTED-LIKE FROM ARABIDOPSIS THALIANA2 (KNAT2) promoter activity. OsGRF overexpression phenotypes in rice were consistent with aberrant meristematic activity, showing reduced formation of tillers and internodes and extensive adventitious root/shoot formation on nodes. These effects were associated with down-regulation of endogenous Oskn2 expression by OsGRF3. Conversely, RNA interference silencing of OsGRF3, OsGRF4, and OsGRF5 resulted in dwarfism, delayed growth and inflorescence formation, and up-regulation of Oskn2. These data demonstrate conserved interactions between the GRF and KNOX families of transcription factors in both monocot and dicot plants.

Accurate detection of subclonal single nucleotide variants in whole genome amplified and pooled cancer samples using HaloPlex target enrichment.

BACKGROUND: Target enrichment and resequencing is a widely used approach for identification of cancer genes and genetic variants associated with diseases. Although cost effective compared to whole genome sequencing, analysis of many samples constitutes a significant cost, which could be reduced by pooling samples before capture. Another limitation to the number of cancer samples that can be analyzed is often the amount of available tumor DNA. We evaluated the performance of whole genome amplified DNA and the power to detect subclonal somatic single nucleotide variants in non-indexed pools of cancer samples using the HaloPlex technology for target enrichment and next generation sequencing. RESULTS: We captured a set of 1528 putative somatic single nucleotide variants and germline SNPs, which were identified by whole genome sequencing, with the HaloPlex technology and sequenced to a depth of 792-1752. We found that the allele fractions of the analyzed variants are well preserved during whole genome amplification and that capture specificity or variant calling is not affected. We detected a large majority of the known single nucleotide variants present uniquely in one sample with allele fractions as low as 0.1 in non-indexed pools of up to ten samples. We also identified and experimentally validated six novel variants in the samples included in the pools. CONCLUSION: Our work demonstrates that whole genome amplified DNA can be used for target enrichment equally well as genomic DNA and that accurate variant detection is possible in non-indexed pools of cancer samples. These findings show that analysis of a large number of samples is feasible at low cost, even when only small amounts of DNA is available, and thereby significantly increases the chances of indentifying recurrent mutations in cancer samples.

The homeodomain-leucine zipper (HD-Zip) class I transcription factors ATHB7 and ATHB12 modulate abscisic acid signalling by regulating protein phosphatase 2C and abscisic acid receptor gene activities.

Plants perceiving drought activate multiple responses to improve survival, including large-scale alterations in gene expression. This article reports on the roles in the drought response of two Arabidopsis thaliana homeodomain-leucine zipper class I genes; ATHB7 and ATHB12, both strongly induced by water-deficit and abscisic acid (ABA). ABA-mediated transcriptional regulation of both genes is shown to depend on the activity of protein phosphatases type 2C (PP2C). ATHB7 and ATHB12 are, thus, targets of the ABA signalling mechanism defined by the PP2Cs and the PYR/PYL family of ABA receptors, with which the PP2C proteins interact. Our results from chromatin immunoprecipitation and gene expression analyses demonstrate that ATHB7 and ATHB12 act as positive transcriptional regulators of PP2C genes, and thereby as negative regulators of abscisic acid signalling. In support of this notion, our results also show that ATHB7 and ATHB12 act to repress the transcription of genes encoding the ABA receptors PYL5 and PYL8 in response to an ABA stimulus. In summary, we demonstrate that ATHB7 and ATHB12 have essential functions in the primary response to drought, as mediators of a negative feedback effect on ABA signalling in the plant response to water deficit.

A genome-wide survey of HD-Zip genes in rice and analysis of drought-responsive family members.

The homeodomain leucine zipper (HD-Zip) genes encode transcription factors that have diverse functions in plant development and have often been implicated in stress adaptation. The HD-Zip genes are the most abundant group of homeobox (HB) genes in plants and do not occur in other eukaryotes. This paper describes the complete annotation of the HD-Zip families I, II and III from rice and compares these gene families with Arabidopsis in a phylogeny reconstruction. Orthologous pairs of rice and Arabidopsis HD-Zip genes were predicted based on neighbour joining and maximum parsimony (MP) trees with support of conserved intron-exon organization. Additionally, a number of HD-Zip genes appeared to be unique to rice. Searching of EST and cDNA databases and expression analysis using RT-PCR showed that 30 out of 31 predicted rice HD-Zip genes are expressed. Most HD-Zip genes were broadly expressed in mature plants and seedlings, but others showed more organ specific patterns. Like in Arabidopsis and other dicots, a subset of the rice HD-Zip I and II genes was found to be regulated by drought stress. We identified both drought-induced and drought-repressed HD-Zip genes and demonstrate that these genes are differentially regulated in drought-sensitive versus drought-tolerant rice cultivars. The drought-repressed HD-Zip family I gene, Oshox4, was selected for promoter-GUS analysis, showing that drought-responsiveness of Oshox4 is controlled by the promoter and that Oshox4 expression is predominantly vascular-specific. Loss-of-function analysis of Oshox4 revealed no specific phenotype, but overexpression analysis suggested a role for Oshox4 in elongation and maturation processes.