Identification of novel genetic alterations in pediatric cytogenetically normal acute myeloid leukemia by next-generation sequencing

Pediatric acute myeloid leukemia (AML) is a molecularly heterogeneous disease that arises from genetic alterations in pathways that regulate self-renewal and myeloid differentiation. While the majority of patients carry recurrent chromosomal translocations, almost 20% of childhood AML do not show any recognizable cytogenetic alteration and are defined as cytogenetically normal (CN)-AML. CN-AML patients have always showed a great variability in response to therapy and overall outcome, underlining the presence of unknown genetic changes, not detectable by conventional analyses, but relevant for pathogenesis, and outcome of AML. The development of novel genome-wide techniques such as next-generation sequencing, have tremendously improved our ability to interrogate the cancer genome. Based on this background, the aim of this research study was to investigate the mutational landscape of pediatric CN-AML patients negative for all the currently known somatic mutations reported in AML through whole-transcriptome sequencing (RNA-seq). RNA-seq performed on diagnostic leukemic blasts from 19 pediatric CN-AML cases revealed a considerable incidence of cryptic chromosomal rearrangements, with the identification of 21 putative fusion genes. Several of the fusion genes that were identified in this study are recurrent and might have a prognostic and/or therapeutic relevance. A paradigm of that is the CBFA2T3-GLIS2 fusion, which has been demonstrated to be a common alteration in pediatric CN-AML, predicting poor outcome. Important findings have been also obtained in the identification of novel therapeutic targets. On one side, the identification of NUP98-JARID1A fusion suggests the use of disulfiram; on the other, here we describe alteration-activating tyrosine kinases, providing functional data supporting the use of tyrosine kinase inhibitors to specifically inhibit leukemia cells. This study provides new insights in the knowledge of genetic alterations underlying pediatric AML, defines novel prognostic markers and putative therapeutic targets, and prospectively ensures a correct risk stratification and risk-adapted therapy also for the “all-neg” AML subgroup.

Defining the regions of homology between human chromosome 16p12--13 and mouse chromosomes

In this study, the evolutionary relationship between human chromosome 16p12-p13 and mouse chromosomes was investigated by determining the order of marker loci in the region and then identifying the chromosomal locations of the homologous loci in mice. Eighteen genes from human 16 were mapped to fifteen subchromosomal regions by a variety of mapping approaches.^ Thirteen of the genes were mapped in the mouse. Linkage analysis with backcross mice and segregation analysis in a mouse - Chinese Hamster Ovary (CHO) somatic cell hybrid panel informative for different regions of mouse genome were used. The results assigned the thirteen genes to three different mouse chromosomes.^ A group of six genes on mouse 16 was found to be closely linked to Scid. The order of Myh11 and Mrp remains ambiguous since no recombination was detected in backcross analysis. Their relative position in human is also uncertain since they were shown to be very close to each other. For the other mouse loci, an unambiguous gene order could be determined and was found to be identical to that in human. Therefore, they comprise a new conserved linkage group between the two species. The orientation of the group was inverted relative to the centromeres, i.e. the proximal loci in one species become distal in another. The size of the group was estimated to be from 4.4 to 8 Mb and 10 to 32 cM in human. In mouse, it was about 21 cM in the backcross analysis. The two boundaries of the conserved linkage were defined within a 1 Mb range. It is now possible to predict the locations of mouse homologs for some human disease genes based on their locations on human 16p.^ The six human 16p genes that map to MMU7 showed a different gene order in mouse than in human. No recombination was found between Crym and Umod while Crym was distal to D16S79A and proximal to D16S92. The location of Stp and Cdr2 with respect to the above four loci was not determined since they were not mapped in the same set of backcross mice. These genes greatly expanded an existing conserved synteny group between the human 16p12-p13 region and the MMU7. It now consists of eleven loci that span a region of probably more than 10 Mb in human. The gene order derived from this study provided further evidence for chromosomal rearrangements within the conserved synteny. (Abstract shortened by UMI.) ^

GCN5 maintains genome stability during development

The histone acetyltransferase, GCN5, is essential for survival of mice during embryogenesis. GCN5 null embryos die early during development due to increased apoptosis. We have demonstrated that the increased apoptosis in associated with increased p53 protein levels. Loss of p53 rescues the embryonic apoptosis in the GCN5 null embryos. These results raised the question of what molecular trigger leads to p53 stabilization and cell death in the absence of GCN5. p53 is generally referred to as the gatekeeper of the cell, monitoring cellular responses to DNA damage, genotoxic stress, and other unfavorable conditions in the cell. Therefore, we examined individual cells in wild type and mutant embryos for gross chromosomal aberrations that might trigger a genome integrity checkpoint. Karyotype analysis indicates that approximately 30% of the cells in an E8.5 GCN5 null embryo display chromosomal aberrations, predominantly chromosomal end adhesions and associations. In wild type E8.5 embryos, only 6% of the cells have chromosomal aberrations. Recent data using telomeric FISH demonstrates that cells from GCN5 null embryos have a decreased telomeric signal. Telomere maintenance is essential for maintaining genome integrity. Telomeric defects are associated with loss of chromosomes and chromosomal rearrangements that can lead to detrimental gene fusions involved in many types of cancers. Little is known about the chromatin structures present near the telomeric ends, or whether any of the telomere-associated proteins are subject to post-translational modification such as acetylation. Our results are the first data to demonstrate the involvement of a histone acetyltransferase, GCN5, in maintaining genome integrity through telomere maintenance and/or capping. ^

Whole genome analysis of Fusobacterium nucleatum, Rickettsia typhi and Francisella tularensis

The genomes of Fusobacterium nucleatum subspecies polymorphum strain ATCC 10953, Rickettsia typhi strain Wilmington, and Francisella tularensis subspecies holarctica strain OSU18 were sequenced, annotated, and analyzed. Each genome was then compared to the sequenced genomes of closely related bacteria. The genome of F. nucleatum ATCC 10953 was compared to two additional F. nucleatum subspecies, subspecies nucleatum and subspecies vincentii. This analysis revealed substantial evidence of horizontal gene transfer along with considerable genetic diversity within the species of F. nucleatum. R. typhi was compared to R. prowazekii and R. conorii. This analysis uncovered a hotspot for chromosomal rearrangements in the Spotted Fever Group but not the Typhus Group Rickettsia and revealed the close genetic relationship between the Typhus Group rickettsial species. F. tularensis OSU18 was compared to two additional F. tularensis strains. These comparisons uncovered significant chromosomal rearrangements between F. tularensis subspecies due to recombination between insertion sequence elements. ^

Modification of EWS/WT1 functional properties by phosphorylation

In many human cancers, tumor-specific chromosomal rearrangements are known to create chimeric products with the ability to transform cells. The EWS/WT1 protein is such a fusion product, resulting from a t(11;22) chromosomal translocation in desmoplastic small round cell tumors, where 265 aa from the EWS amino terminus are fused to the DNA binding domain of the WT1 tumor suppressor gene. Herein, we find that EWS/WT1 is phosphorylated in vivo on serine and tyrosine residues and that this affects DNA binding and homodimerization. We also show that EWS/WT1 can interact with, and is a substrate for, modification on tyrosine residues by c-Abl. Tyrosine phosphorylation of EWS/WT1 by c-Abl negatively regulates its DNA binding properties. These results indicate that the biological activity of EWS/WT1 is closely linked to its phosphorylation status.

Posttranslational modification of TEL and TEL/AML1 by SUMO-1 and cell-cycle-dependent assembly into nuclear bodies

The E-26 transforming specific (ETS)-related gene, TEL, also known as ETV6, encodes a strong transcription repressor that is rearranged in several recurring chromosomal rearrangements associated with leukemia and congenital fibrosarcoma. TEL is a nuclear phosphoprotein that is widely expressed in all normal tissues. TEL contains a DNA-binding domain at the C terminus and a helix–loop–helix domain (also called a pointed domain) at the N terminus. The pointed domain is necessary for homotypic dimerization and for interaction with the ubiquitin-conjugating enzyme UBC9. Here we show that the interaction with UBC9 leads to modification of TEL by conjugating it to SUMO-1. The SUMO-1-modified TEL localizes to cell-cycle-specific nuclear speckles that we named TEL bodies. We also show that the leukemia-associated fusion protein TEL/AML1 is modified by SUMO-1 and found in the TEL bodies, in a pattern quite different from what we observe and report for AML1. Therefore, SUMO-1 modification of TEL could be a critical signal necessary for normal functioning of the protein. In addition, the modification by SUMO-1 of TEL/AML1 could lead to abnormal localization of the fusion protein, which could have consequences that include contribution to neoplastic transformation.

Normal cardiovascular development in mice deficient for 16 genes in 550 kb of the velocardiofacial/ DiGeorge syndrome region

Hemizygous interstitial deletions in human chromosome 22q11 are associated with velocardiofacial syndrome and DiGeorge syndrome and lead to multiple congenital abnormalities, including cardiovascular defects. The gene(s) responsible for these disorders is thought to reside in a 1.5-Mb region of 22q11 in which 27 genes have been identified. We have used Cre-mediated recombination of LoxP sites in embryonic stem cells and mice to generate a 550-kb deletion encompassing 16 of these genes in the corresponding region on mouse chromosome 16. Mice heterozygous for this deletion are normal and do not exhibit cardiovascular abnormalities. Because mice with a larger deletion on mouse chromosome 16 do have heart defects, the results allow us to exclude these 16 genes as being solely, or in combination among themselves, responsible for the cardiovascular abnormalities in velocardiofacial/DiGeorge syndrome. We also generated mice with a duplication of the 16 genes that may help dissect the genetic basis of “cat eye” and derivative 22 syndromes that are characterized by extra copies of portions of 22q11, including these 16 genes. We also describe a strategy for selecting cell lines with defined chromosomal rearrangements. The method is based on reconstitution of a dominant selection marker after Cre-mediated recombination of LoxP sites. Therefore it should be widely applicable to many cell lines.

Deletion of long-range regulatory elements upstream of SOX9 causes campomelic dysplasia

Campomelic dysplasia (CD) is a rare, neonatal human chondrodysplasia characterized by bowing of the long bones and often associated with male-to-female sex-reversal. Patients present with either heterozygous mutations in the SOX9 gene or chromosome rearrangements mapping at least 50 kb upstream of SOX9. Whereas mutations in SOX9 ORF cause haploinsufficiency, the effects of translocations 5′ to SOX9 are unclear. To test whether these rearrangements also cause haploinsufficiency by altering spatial and temporal expression of SOX9, we generated mice transgenic for human SOX9-lacZ yeast artificial chromosomes containing variable amounts of DNA sequences upstream of SOX9. We show that elements necessary for SOX9 expression during skeletal development are highly conserved between mouse and human and reveal that a rearrangement upstream of SOX9, similar to those observed in CD patients, leads to a substantial reduction of SOX9 expression, particularly in chondrogenic tissues. These data demonstrate that important regulatory elements are scattered over a large region upstream of SOX9 and explain how particular aspects of the CD phenotype are caused by chromosomal rearrangements 5′ to SOX9.

The chimeric genes AML1/MDS1 and AML1/EAP inhibit AML1B activation at the CSF1R promoter, but only AML1/MDS1 has tumor-promoter properties.

The (3;21)(q26;q22) translocation associated with treatment-related myelodysplastic syndrome, treatment-related acute myeloid leukemia, and blast crisis of chronic myeloid leukemia results in the expression of the chimeric genes AML1/EAP, AML1/MDS1, and AML1/EVI1. AML1 (CBFA2), which codes for the alpha subunit of the heterodimeric transcription factor CBF, is also involved in the t(8;21), and the gene coding for the beta subunit (CBFB) is involved in the inv(16). These are two of the most common recurring chromosomal rearrangements in acute myeloid leukemia. CBF corresponds to the murine Pebp2 factor, and CBF binding sites are found in a number of eukaryotic and viral enhancers and promoters. We studied the effects of AML1/EAP and AML1/MDS1 at the AML1 binding site of the CSF1R (macrophage-colony-stimulating factor receptor gene) promoter by using reporter gene assays, and we analyzed the consequences of the expression of both chimeric proteins in an embryonic rat fibroblast cell line (Rat1A) in culture and after injection into athymic nude mice. Unlike AML1, which is an activator of the CSF1R promoter, the chimeric proteins did not transactivate the CSF1R promoter site but acted as inhibitors of AML1 (CBFA2). AML1/EAP and AML1/MDS1 expressed in adherent Rat1A cells decreased contact inhibition of growth, and expression of AML1/MDS1 was associated with acquisition of the ability to grow in suspension culture. Expression of AML1/MDS1 increased the tumorigenicity of Rat1A cells injected into athymic nude mice, whereas AML1/EAP expression prevented tumor growth. These results suggest that expression of AML1/EAP and AML1/MDS1 can interfere with normal AML1 function, and that AML1/MDS1 has tumor-promoting properties in an embryonic rat fibroblast cell line.

Frequent somatic hypermutation of the 5' noncoding region of the BCL6 gene in B-cell lymphoma.

The BCL6 gene encodes a zinc-finger transcription factor and is altered by chromosomal arrangements in its 5' noncoding region in approximately 30% of diffuse large-cell lymphoma (DLCL). We report here that, in 22/30 (73%) DLCL and 7/15 (47%) follicular lymphoma (FL), but not in other tumor types, the BCL6 gene is also altered by multiple (1.4 x 10(-3) -1.6 x 10(-2) per bp), often biallelic, mutations clustering in its 5' noncoding region. These mutations are of somatic origin and are found in cases displaying either normal or rearranged BLC6 alleles indicating their independence from chromosomal rearrangements and linkage to immunoglobulin genes. These alterations identify a mechanism of genetic instability in malignant B cells and may have been selected during lymphomagenesis for their role in altering BCL6 expression.

Fusion of the TEL gene on 12p13 to the AML1 gene on 21q22 in acute lymphoblastic leukemia.

Chromosomal rearrangements involving band 12p13 are found in a wide variety of human leukemias but are particularly common in childhood acute lymphoblastic leukemia. The genes involved in these rearrangements, however, have not been identified. We now report the cloning of a t(12;21) translocation breakpoint involving 12p13 and 21q22 in two cases of childhood pre-B acute lymphoblastic leukemia, in which t(12;21) rearrangements were not initially apparent. The consequence of the translocation is fusion of the helix-loop-helix domain of TEL, an ETS-like putative transcription factor, to the DNA-binding and transactivation domains of the transcription factor AML1. These data show that TEL, previously shown to be fused to the platelet-derived growth factor receptor beta in chronic myelomonocytic leukemia, can be implicated in the pathogenesis of leukemia through its fusion to either a receptor tyrosine kinase or a transcription factor. The TEL-AML1 fusion also indicates that translocations affecting the AML1 gene can be associated with lymphoid, as well as myeloid, malignancy.

The su(Hw) protein bound to gypsy sequences in one chromosome can repress enhancer-promoter interactions in the paired gene located in the other homolog.

The suppressor of Hairy-wing [su(Hw)] protein exerts a polar effect on gene expression by repressing the function of transcriptional enhancers located distally from the promoter with respect to the location of su(Hw) binding sequences. The directionality of this effect suggests that the su(Hw) protein specifically interferes with the basic mechanism of enhancer action. Moreover, mutations in modifier of mdg4 [mod(mdg4)] result in the repression of expression of a gene when the su(Hw) protein is bound to sequences in the copy of this gene located in the homologous chromosome. This effect is dependent on the presence of the su(Hw) binding region from the gypsy retrotransposon in at least one of the chromosomes and is enhanced by the presence of additional gypsy sequences in the other homology. This phenomenon is inhibited by chromosomal rearrangements that disrupt pairing, suggesting that close apposition between the two copies of the affected gene is important for trans repression of transcription. These results indicate that, in the absence of mod-(mdg4) product, the su(Hw) protein present in one chromosome can act in trans and inactivate enhancers located in the other homolog.

Caracterização de rearranjos cromossômicos citogeneticamente equilibrados associados a quadros clínicos

Este estudo teve como objetivos (a) identificar mecanismos pelos quais rearranjos cromossômicos citogeneticamente equilibrados possam estar associados de maneira causal a determinados quadros clínicos e (b) contribuir para a compreensão dos mecanismos de formação desses rearranjos. Para isso, foram estudados 45 rearranjos cromossômicos citogeneticamente equilibrados (29 translocações, 10 inversões e seis rearranjos complexos), detectados em pacientes que apresentavam malformações congênitas, comprometimento do desenvolvimento neuropsicomotor ou déficit intelectual. Foram 31 rearranjos cromossômicos esporádicos, três familiais que segregavam com o quadro clínico e mais 11 rearranjos cromossômicos herdados de genitores fenotipicamente normais. Inicialmente os pontos de quebra desses rearranjos foram mapeados por hibridação in situ fluorescente (FISH). A busca por microdeleções e duplicações genômicas foi realizada por a-CGH. A investigação dos pontos de quebra prosseguiu com a aplicação da técnica de Mate-Pair Sequencing (MPS), que permite localizar as quebras em segmentos de 100 pb - 1 kb, na maioria dos casos. Para obter os segmentos de junção das quebras no nível de pares de bases, os segmentos delimitados por MPS foram sequenciados pelo método de Sanger. A análise por aCGH revelou microdeleções ou microduplicações localizadas nos cromossomos rearranjados, em 12 dos 45 pacientes investigados (27%). A análise de 27 rearranjos por MPS permitiu a caracterização dos pontos de junção das quebras. MPS expandiu o número de pontos de quebra, detectados por análise do cariótipo ou aCGH, de 114 para 156 (em resolução < 2kb, na maioria dos casos). O número de pontos de quebra/rearranjo variou de 2 a 20. Os 156 pontos de quebra resultaram em 86 variantes estruturais equilibradas e outras 32 variantes não equilibradas. Perdas e ganhos de segmentos submiscroscópicos nos cromossomos rearranjados constituíram a principal causa ou, provavelmente, contribuíram para o quadro clínico de 12 dos 45 pacientes. Em cinco desses 12 rearranjos foram detectadas por MPS a interrupção de genes já relacionados à doença, ou provável alteração de sua região reguladora, contribundo para o quadro clínico. Em quatro dos 33 rearranjos não associados a perdas ou ganhos de segmentos, a análise por MPS revelou a interrupção de genes que já foram anteriormente relacionados a doenças, explicando-se, assim, as características clínicas dos portadores; outro rearranjo pode ter levando alteração da expressão gênica de gene sensível a dosagem e ao quadro clínico. Um rearranjo cromossômico familial, identificado na análise após bandamento G como uma translocação equilibrada, t(2;22)(p14;q12), segregava com quadro de atraso do desenvolvimento neuropsicomotor e dificuldade de aprendizado associados a dismorfismos. A combinação das análises por FISH, aCGH e MPS revelou que se tratava, na verdade, de rearranjo complexo entre os cromossomos 2, 5 e 22, incluindo 10 quebras. A segregação de diferentes desequilíbrios submicroscópicos em indivíduos afetados e clinicamente normais permitiu a compreensão da variabilidade clínica observada na família. Rearranjos equilibrados detectados em indivíduos afetados, mas herdados de genitores clinicamente normais, são, em geral, considerados como não tendo relação com o quadro clínico, apesar da possibilidade de desequilíbrios cromossômicos gerados por permuta desigual na meiose do genitor portador do rearranjo. Neste trabalho, a investigação de 11 desses rearranjos por aCGH não revelou perdas ou ganhos de segmentos nos cromossomos rearranjados. No entanto, a análise por aCGH da portadora de um desses rearranjos - inv(12)mat - revelou deleção de 8,7 Mb no cromossomo 8, como causa de seu fenótipo clínico. Essa deleção estava relacionada com outro rearranjo equilibrado também presente em sua mãe, independente da inversão. Para compreender os mecanismos de formação de rearranjos citogeneticamente equilibrados, investigamos os segmentos de junção no nível de pares de base. A análise por MPS que levou, na maioria dos casos, ao mapeamento dos pontos de quebras em segmentos <1kb permitiu o sequenciamento pelo método de Sanger de 51 segmentos de junções de 17 rearranjos. A ocorrência de blunt fusions ou inserções e deleções <10 pb, e a ausência de homologia ou a presença de micro homologia de 2 pb a 4 pb de extensão indicaram o mecanismo de junção de extremidades não homólogas (non-homologous end joinging; NHEJ), na maioria das 51 junções caracterizadas. As características de três dos quatro rearranjos mais complexos, com 17-20 quebras, indicaram sua formação pelo mecanismo de chromothripsis. Este estudo mostra a importância da análise genômica de variações de número de cópias por microarray, juntamente com o mapeamento dos pontos de quebra por MPS, para determinar a estrutura de rearranjos cromossômicos citogeneticamente equilibrados e seu impacto clínico. O mapeamento dos segmentos de junção por MPS, permitindo o sequenciamento pelo método de Sanger, foi essencial para a compreensão de mecanismos de formação desses rearranjos

Padrões temporais e grau de diversificação cariotípica em espécies atlânticas da família Acanthuridae (Perciformes)

The Acanthuridae family is a representative group from the marine fish that plays a key role in ecological dynamics of coral reefs. Three species are common along coastal reefs of Western Atlantic: Acanthurus coeruleus, Acanthurus bahianus and Acanthurus chirurgus. In the present study, cytogenetic data are presented for these three species Acanthurus based on classical cytogenetic methods and mapping of repetitive sequences such as ribosomal 18S and 5S rDNA and telomeric repeats to improve their karyotype evolutionary analyses. The cytogenetic pattern of these species indicated sequential steps of chromosomal rearrangements dating back 19 to 5 millions of years ago (M.a.) that accounted for their interspecific differences. A. coeruleus (2n=48; 2sm+4st+42a), A. bahianus (2n=36; 12m+2sm+4st+18a) and A. chirurgus (2n=34; 12m+2sm+4st+16a) share an older set of three chromosomal pairs that were originated through pericentric inversions. A set of six large metacentric pairs formed by Robertsonian (Rb) translocations found in A. bahianus and A. chirurgus and a putative in tandem fusion found in A. chirurgus are more recent events. The lack of interstitial telomeric sequences (ITS) in spite of several centric fusions in A. bahianus and A. chirurgus might be related to the long period of time after their occurrence (estimated in 5 M.a.). Furthermore, the homeologies among the chromosome pairs bearing ribosomal genes, in addition to other structural features, highlight large conserved chromosomal regions in the three species. Our findings indicate that macrostructural changes occurred during the cladogenesis of these species were not followed by conspicuous microstructural rearrangements in the karyotypes.

How does DNA break during chromosomal translocations?

Chromosomal translocations are one of the most common types of genetic rearrangements and are molecular signatures for many types of cancers. They are considered as primary causes for cancers, especially lymphoma and leukemia. Although many translocations have been reported in the last four decades, the mechanism by which chromosomes break during a translocation remains largely unknown. In this review, we summarize recent advances made in understanding the molecular mechanism of chromosomal translocations.