A rare KMT2A::CBL transcript in an acute monoblastic leukemia patient with an unfavorable outcome

BACKGROUND

Lysine [K] methyltransferase 2A (KMT2A, previously known as MLL) gene rearrangements are common in acute leukemias of various lineages and are associated with features such as chemotherapy resistance and rapid relapse. KMT2A::CBL is a rare fusion of unknown pathogenesis generated by a unique interstitial deletion of chromosome 11 that has been reported across a wide age range in both acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL) patients. The leukemogenic effect of the KMT2A::CBL rearrangement and its association with clinical prognosis have not been well clarified.

METHODS AND RESULTS

We report the case of a 64-year-old female who was diagnosed with acute monoblastic leukemia (M5a) and who acquired the rare KMT2A::CBL fusion. The patient received multiple cycles of therapy but did not achieve remission and eventually succumbed to severe infection and disease progression. Additionally, we characterized the predicted KMT2A-CBL protein structure in this case to reveal the underlying leukemogenic mechanisms and summarized reported cases of hematological malignancies with KMT2A::CBL fusion to investigate the correlation of gene rearrangements with clinical outcomes.

CONCLUSIONS

This report provides novel insights into the leukemogenic potential of the KMT2A::CBL rearrangement and the correlation between gene rearrangements and clinical outcomes.

Large-scale genomic rearrangements boost SCRaMbLE in Saccharomyces cerevisiae

Synthetic Chromosome Rearrangement and Modification by LoxP-mediated Evolution (SCRaMbLE) is a promising tool to study genomic rearrangements. However, the potential of SCRaMbLE to study genomic rearrangements is currently hindered, because a strain containing all 16 synthetic chromosomes is not yet available. Here, we construct SparLox83R, a yeast strain containing 83 loxPsym sites distributed across all 16 chromosomes. SCRaMbLE of SparLox83R produces versatile genome-wide genomic rearrangements, including inter-chromosomal events. Moreover, when combined with synthetic chromosomes, SCRaMbLE of hetero-diploids with SparLox83R leads to increased diversity of genomic rearrangements and relatively faster evolution of traits compared to hetero-diploids only with wild-type chromosomes. Analysis of the SCRaMbLEd strain with increased tolerance to nocodazole demonstrates that genomic rearrangements can perturb the transcriptome and 3D genome structure and consequently impact phenotypes. In summary, a genome with sparsely distributed loxPsym sites can serve as a powerful tool for studying the consequence of genomic rearrangements and accelerating strain engineering in Saccharomyces cerevisiae.

Evaluation of genetic risk of apparently balanced chromosomal rearrangement carriers by breakpoint characterization

PURPOSE

To report genetic characteristics and associated risk of chromosomal breaks due to chromosomal rearrangements in large samples.

METHODS

MicroSeq, a technique that combines chromosome microdissection and next-generation sequencing, was used to identify chromosomal breakpoints. Long-range PCR and Sanger sequencing were used to precisely characterize 100 breakpoints in 50 ABCR carriers.

RESULTS

In addition to the recurrent regions of balanced rearrangement breaks in 8q24.13, 11q11.23, and 22q11.21 that had been documented, we have discovered a 10-Mb region of 12q24.13-q24.3 that could potentially be a sparse region of balanced rearrangement breaks. We found that 898 breakpoints caused gene disruption and a total of 188 breakpoints interrupted genes recorded in OMIM. The percentage of breakpoints that disrupted autosomal dominant genes recorded in OMIM was 25.53% (48/188). Fifty-four of the precisely characterized breakpoints had 1-8-bp microhomologous sequences.

CONCLUSION

Our findings provide a reference for the evaluation of the pathogenicity of mutations in related genes that cause protein truncation in clinical practice. According to the characteristics of breakpoints, non-homologous end joining and microhomology-mediated break-induced replication may be the main mechanism for ABCRs formation.

The analyses of the complete mitochondrial genomes of three crabs revealed novel gene rearrangements and phylogenetic relationships of Brachyura

BACKGROUND

Brachyura crab is the largest branch of Decapoda crustacean. Phylogenetic relationships within Brachyura remain controversial to be investigated. The mitochondrial genome (mitogenome) is an important molecular marker for studying the phylogenetic relationships of Brachyura.

METHODS AND RESULTS

To understand the phylogeny of Brachyura, the three complete mitogenomes from Charybdis annulata, Leptodius exaratus, and Spider crab were sequenced and annotated. Their full length was 15,747, 15,716, and 16,608 bp long, respectively. The first two crabs both contained 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. However, Spider crab contained 13 PCGs, two rRNA genes, 25 tRNA genes and a control region. The mitogenomes of each of the three crabs exhibited high AT content (67.8%, 69.1%, and 70.8%), with negative AT skews (-0.014, - 0.028, and - 0.017) and GC skews (-0.269, - 0.286, and - 0.341). The gene order of C. annulata was identical to the ancestor of Brachyura. Compared with the ancestor of Brachyura, L. exaratus exhibited the gene rearrangements of Val (V)-rrnS-control region, and Spider crab had the four copies of Lys (K). Phylogenetic analyses indicated that C. annulata belonged to Portunidae family, Portunoidea superfamilies, L. exaratus belonged to Xanthidae family, Xanthoidea superfamilies, and Spider crab belonged to Mithracidae family, Majoidea superfamilies. Phylogenetic analyses showed that the two species (Somanniathelphusa boyangensis and Huananpotamon lichuanense) belonging to the Potamoidea were sister groups to the Thoracotremata, thus supporting the conclusion that Heterotremata is polyphyletic.

CONCLUSION

The results of this study enriched the crab mitogenome database and enabled us to better understand the phylogenetic relationships of Brachyura.

Enhancers within the Ig V Gene Region Orchestrate Chromatin Topology and Regulate V Gene Rearrangement Frequency to Shape the B Cell Receptor Repertoire Specificities

Effective Ab-mediated responses depend on a highly diverse Ab repertoire with the ability to bind a wide range of epitopes in disease-causing agents. The generation of this repertoire depends on the somatic recombination of the variable (V), diversity (D), and joining (J) genes in the Ig loci of developing B cells. It has been known for some time that individual V, D, and J gene segments rearrange at different frequencies, but the mechanisms behind this unequal V gene usage have not been well understood. However, recent work has revealed that newly described enhancers scattered throughout the V gene-containing portion of the Ig loci regulate the V gene recombination frequency in a regional manner. Deletion of three of these enhancers revealed that these elements exert many layers of control during V(D)J recombination, including long-range chromatin interactions, epigenetic milieu, chromatin accessibility, and compartmentalization.

Comprehensive profiling of pathogenic germline large genomic rearrangements in a pan-cancer analysis

The presence of large genomic rearrangements (LGRs) has been heavily investigated in breast and ovarian cancer. However, correlations between LGRs and cancer types beyond these two have not been extensively profiled, likely due to the highly inefficient methods of detecting these types of alterations. This study utilized next-generation sequencing (NGS) to analyze and classify the germline LGR profile in 17 025 cancer patients across 22 cancer types. We characterized newly identified LGRs based on predicted pathogenicity and took a closer look at genes that acquire both germline and somatic mutations within our samples. The detection method for LGRs was validated using droplet digital polymerase chain reaction (ddPCR) assay of commonly investigated LGR genes. In total, 15 659 samples from across 22 cancer types were retained for analysis after filtering. We observed that, in our cohort, the cancer types with the highest proportion of germline LGRs were ovarian cancer (4.7%), renal cell carcinoma (2.5%), breast cancer (2%), glioma (1.8%) and thyroid carcinoma (1.8%). Annotation of detected germline variants revealed several genes-MSH2, FANCA and PMS2-that contain novel LGRs. We observed co-occurrences between germline LGRs in MSH2 and somatic single nucleotide variants/insertion and deletions (SNVs/InDels) in BRCA2, KTM2B, KDM5A, CHD8, and HNF1A. Furthermore, our analysis showed that samples with pathogenic and likely pathogenic germline LGRs tended to also have higher mutational burden, chromosomal instability, and microsatellite instability ratio compared to samples with pathogenic germline SNVs/InDels. In this study, we demonstrated the prevalence of pathogenic germline LGRs beyond breast and ovarian cancer. The profiles of these pathogenic or likely pathogenic alterations will fuel further investigations and highlight new understanding of LGRs across multiple cancer types.

B-lymphoblastic leukemia/lymphoma with MYC and BCL2 gene rearrangements shows evidence for clonal evolution and mitotic recombination

BACKGROUND

B-lymphoblastic leukemia/lymphomas (B-ALL/LBL) are uncommon neoplasms that may be associated with a variety of cytogenetic and molecular changes. The mechanisms by which these changes arise have not been fully described.

AIMS/PURPOSE

This report describes an unusual case of B-ALL/LBL with complex clonal evolution that includes BCL2 and MYC gene rearrangements.

METHODS

Immunophenotyping was performed by immunohistochemistry and flow cytometry. Traditional G-band karyotyping was accompanied by fluorescence in-situ hybridization (FISH) using break-apart and dual fusion probes. Single nucleotide polymorphisms were assessed using a high-density DNA microarray.

RESULTS

The karyotype of the blasts showed reciprocal translocation of chromosomes 4 and 18, reciprocal translocation of chromosomes 8 and 14 with two copies of the oncogenic translocation derivative(14)t(8;14), and no normal chromosome 14. FISH studies showed complex IGH-BCL2 and IGH-MYC fusion signals.

CONCLUSIONS

A clonal evolution model involving multiple chromosomal translocations and mitotic recombination is postulated to account for the karyotype, FISH, and microarray results but leaves unresolved the exact order of the evolutionary changes.

A Cre-LoxP-based approach for combinatorial chromosome rearrangements in human HAP1 cells

Alterations of human karyotype caused by chromosomal rearrangements are often associated with considerable phenotypic effects. Studying molecular mechanisms underlying these effects requires an efficient and scalable experimental model. Here, we propose a Cre-LoxP-based approach for the generation of combinatorial diversity of chromosomal rearrangements. We demonstrate that using the developed system, both intra- and inter-chromosomal rearrangements can be induced in the human haploid HAP1 cells, although the latter is significantly less effective. The obtained genetically modified HAP1 cell line can be used to dissect genomic effects associated with intra-chromosomal structural variations.

Small RNA-mediated genome rearrangement pathways in ciliates

Most eukaryotes employ a combination of transcriptional and post-transcriptional silencing mechanisms to suppress transposons, yet ciliates employ a more extreme approach. They separate germline and somatic functions into distinct nuclei, enabling the elimination of transposons from the active somatic genome through diverse small RNA-mediated genome rearrangement pathways during sexual processes.

Complex genomic rearrangements: an underestimated cause of rare diseases

Complex genomic rearrangements (CGRs) are known contributors to disease but are often missed during routine genetic screening. Identifying CGRs requires (i) identifying copy number variants (CNVs) concurrently with inversions, (ii) phasing multiple breakpoint junctions incis, as well as (iii) detecting and resolving structural variants (SVs) within repeats. We demonstrate how combining cytogenetics and new sequencing methodologies is being successfully applied to gain insights into the genomic architecture of CGRs. In addition, we review CGR patterns and molecular features revealed by studying constitutional genomic disorders. These data offer invaluable lessons to individuals interested in investigating CGRs, evaluating their clinical relevance and frequency, as well as assessing their impact(s) on rare genetic diseases.

Two complete mitochondrial genomes in Scolopendra and a comparative analysis of tRNA rearrangements in centipedes

BACKGROUND

Centipedes are one of the oldest terrestrial arthropods belonging to the sub phylum Myriapoda. With the expansion of our understanding of the application of the two centipedes Scolopendra morsitans and Scolopendra hainanum, belonging to the order Scolopendromorpha, an exhaustive classification was required. Although consensus has been reached on the phylogeny of Chilopoda based on morphological traits, recent analyses based on molecular data exhibited differences in results.

METHODS AND RESULTS

The mitochondrial genome sequences of S. morsitans and S. hainanum were obtained by next-generation sequencing. S. morsitans contains 13 PCGs, two rRNAs, 11 tRNAs, and one CR. whereas S. hainanum contains 12 PCGs, of which ATP8 remains unpredicted, two rRNAs, 14 tRNAs, and one CR. An obvious tRNA rearrangement was found in the genus Scolopendra. S. morsitans exhibited a loss of trnW, trnC, trnI, trnK, trnD, trnA, trnN, trnQ, trnF, trnT, trnS, trnL, and trnV, and a repeat of trnR and trnL. S. hainanum exhibited a loss of trnQ, trnC, trnW, trnI, trnD, trnQ, trnP, and trnV. Phylogenetic analyses of centipedes based on 12 PCGs supported the sister relationship between the orders Geophilomorpha and Lithobiomorpha and a close relationship between Scolopendra dehaani and S. hainanum.

CONCLUSIONS

The new mitogenomes determined in this study provide new genomic resources for gene rearrangements and contribute to the understanding of the evolution of gene rearrangement in Chilopoda.

Next-Generation Sequencing-Based Clonality Detection of Immunoglobulin Gene Rearrangements in B-Cell Lymphoma

Immunoglobulin (IG) clonality assessment is a widely used supplementary test for the diagnosis of suspected lymphoid malignancies. The specific rearrangements of the immunoglobulin (IG) heavy and light chain genes act as a unique hallmark of a B-cell lymphoma, a feature that is used in clonality assessment. The widely used BIOMED-2/EuroClonality IG clonality assay, visualized by GeneScanning or heteroduplex analysis, has an unprecedented high detection rate because of the complementarity of this approach. However, the BIOMED-2/EuroClonality clonality assays have been developed for the assessment of specimens with optimal DNA quality. Further improvements for the assessment of samples with suboptimal DNA quality, such as from formalin-fixed paraffin-embedded (FFPE) specimens or specimens with a limited tumor burden, are required. The EuroClonality-NGS Working Group recently developed a next-generation sequencing (NGS)-based clonality assay for the detection of the IG heavy and kappa light chain rearrangements, using the same complementary approach as in the conventional assay. By employing next-generation sequencing, both the sensitivity and specificity of the clonality assay have increased, which not only is very useful for diagnostic clonality testing but also allows robust comparison of clonality patterns in a patient with multiple lymphoma's that have suboptimal DNA quality. Here, we describe the protocols for IG-NGS clonality assessment that are compatible for Ion Torrent and Illumina sequencing platforms including pre-analytical DNA isolation, the analytical phase, and the post-analytical data analysis.

The Use of Fluorescence in situ Hybridization to Confirm PRKACA Gene Rearrangement in Fibrolamellar Hepatocellular Carcinoma: A Validation Study

OBJECTIVE

The objectives of this study are to define the specificity of the DNAJB1-PRKACA fusion transcript for the fibrolamellar subtype of hepatocellular carcinoma (FL-HCC) by testing a targeted sampling of other hepatic neoplasms/proliferations and extrahepatic neoplasms seen in children and young adults and to develop a FISH assay using a commercially available PRKACA break apart probe for use in a CLIA-certified clinical laboratory.

METHODS

Formalin fixed paraffin embedded tissue sections from 12 FL-HCC cases, 142 cases of other hepatic neoplasms/proliferations (conventional HCC, focal nodular hyperplasia (FNH), hepatocellular adenoma (HA) and hepatoblastoma (HB)) and extrahepatic neoplasms (neuroblastoma (NB), Wilms tumor (WT) and Gastrointestinal neuroendocrine tumor (GNET)) and 60 matched background normal control tissues underwent fluorescence in situ hybridization (FISH) testing using a break apart probe targeting the PRKACA gene locus on chromosome 19 using standard techniques.

RESULTS

The PRKACA gene rearrangement was detected in 11/12 (92%) FL-HCC cases and 1/94 (1%) of conventional HCC cases. All other cases and background control tissues were negative for the PRKACA gene rearrangement. These findings establish a test sensitivity of 91.7% and specificity of 99.5%.

CONCLUSION

This study shows that, using standard techniques, FISH testing with a commercially available break apart probe targeting the PRKACA gene can be used as a surrogate for the DNAJB1-PRKACA fusion commonly found in FL-HCC. Also, the PRKACA gene rearrangement is not expressed in other hepatic neo-plasms/proliferations or extrahepatic neoplasms seen in children and young adults. Finally, FISH testing can be used as a diagnostic tool to confirm the diagnosis of FL-HCC, in the appropriate clinical setting.

Single-cell transcriptomics reveals a distinct developmental state of KMT2A-rearranged infant B-cell acute lymphoblastic leukemia

KMT2A-rearranged infant ALL is an aggressive childhood leukemia with poor prognosis. Here, we investigated the developmental state of KMT2A-rearranged infant B-cell acute lymphoblastic leukemia (B-ALL) using bulk messenger RNA (mRNA) meta-analysis and examination of single lymphoblast transcriptomes against a developing bone marrow reference. KMT2A-rearranged infant B-ALL was uniquely dominated by an early lymphocyte precursor (ELP) state, whereas less adverse NUTM1-rearranged infant ALL demonstrated signals of later developing B cells, in line with most other childhood B-ALLs. We compared infant lymphoblasts with ELP cells and revealed that the cancer harbored hybrid myeloid-lymphoid features, including nonphysiological antigen combinations potentially targetable to achieve cancer specificity. We validated surface coexpression of exemplar combinations by flow cytometry. Through analysis of shared mutations in separate leukemias from a child with infant KMT2A-rearranged B-ALL relapsing as AML, we established that KMT2A rearrangement occurred in very early development, before hematopoietic specification, emphasizing that cell of origin cannot be inferred from the transcriptional state.

Comparison of myeloid neoplasms with nonclassic 3q26.2/MECOM versus classic inv(3)/t(3;3) rearrangements reveals diverse clinicopathologic features, genetic profiles, and molecular mechanisms of MECOM activation

MECOM rearrangements are recurrent in myeloid neoplasms and associated with poor prognosis. However, only inv(3)(q21q26.2) and t(3;3)(q21;q26.2), the classic MECOM rearrangements resulting in RPN1-MECOM rearrangement with Mecom overexpression and GATA2 haploinsufficiency, define the distinct subtype of acute myeloid leukemia (AML), and serve as presumptive evidence for myelodysplastic syndrome based on the current World Health Organization classification. Myeloid neoplasms with nonclassic 3q26.2/MECOM rearrangements have been found to be clinically aggressive, but comparative analysis of clinicopathologic and genomic features is limited. We retrospectively studied cohorts of myeloid neoplasms with classic and nonclassic MECOM rearrangements. Cases with classic rearrangements consisted predominantly of AML, often with inv(3) or t(3;3) as the sole chromosome abnormality, whereas the group of nonclassic rearrangements included a variety of myeloid neoplasms, often with complex karyotype without TP53 mutations and similarly dismal overall survival. Immunohistochemistry revealed Mecom protein overexpression in both groups, but overexpression in cases with nonclassic rearrangements was mediated through a mechanism other than GATA2 distal enhancer involvement typical for classic rearrangement. Our results demonstrated that myeloid neoplasms with nonclassic 3q26.2/MECOM rearrangements encompass a diverse group of diseases with poor clinical outcome, overexpression of Mecom protein as a result of the nonclassic mechanism of MECOM activation.

Comparative Analysis of Mitochondrial Genomes in Two Subspecies of the Sunwatcher Toad-Headed Agama (Phrynocephalus helioscopus): Prevalent Intraspecific Gene Rearrangements in Phrynocephalus

Intraspecific rearrangements of mitochondrial genomes are rarely reported in reptiles, even in vertebrates. The sunwatcher toad-headed agama, Phryncoephalus helioscopus, can serve as an excellent model for investigating the dynamic mitogenome structure at intraspecific level. To date, seven subspecies of P. helioscopus are well recognized, but little is known about the mitogenomic evolution among different subspecies. In this study, complete mitogenomes of subspecies P. helioscopus varius II and P. helioscopus cameranoi were determined by next-generation sequencing, and another P. helioscopus varius I retrieved from GenBank was compiled for comparative analysis. The nucleotide composition and the codon usage are similar to those previously published from toad-headed agamas. P. helioscopus varius II and P. helioscopus cameranoi have 23 tRNA genes, including standard 22 tRNA genes and one extra tRNA-Phe (tRNA-Phe duplication). Gene order and phylogenetic analyses in the genus Phrynocephalus support prevalent intraspecific gene rearrangement in P. helioscopus and other congener species including P. erythrurus, P. vlangalii, and P. forsythii. Six different mitochondrial gene arrangements are observed in Phrynocephalus. Overall, the occurrence of rearrangements may result from multiple independent structural dynamic events. The split of the two subspecies in P. helioscopus was dated at approximately 2.34 million years ago (Ma). Two types of gene rearrangements are found in the three mitogenomes of P. helioscopus, and this intraspecific rearrangement phenomenon can be explained by the tandem duplication/random loss (TDRL) model. Post duplication, the alternative loss types can occur in 0.23–0.72 Ma, suggesting that the duplication and fixation of these rearrangements can occur quite quickly. These findings highlight the need for more mitogenomes at the population level in order to better understand the potentially rampant intraspecific mitogenomic reorganization in Phrynocephalus.

Distribution of copy number variations and rearrangement endpoints in human cancers with a review of literature

Copy number variations (CNVs) which include deletions, duplications, inversions, translocations, and other forms of chromosomal re-arrangements are common to human cancers. In this report we investigated the pattern of these variations with the goal of understanding whether there exist specific cancer signatures. We used re-arrangement endpoint data deposited on the Catalogue of Somatic Mutations in Cancers (COSMIC) for our analysis. Indeed, we find that human cancers are characterized by specific patterns of chromosome rearrangements endpoints which in turn result in cancer specific CNVs. A review of the literature reveals tissue specific mutations which either drive these CNVs or appear as a consequence of CNVs because they confer an advantage to the cancer cell. We also identify several rearrangement endpoints hotspots that were not previously reported. Our analysis suggests that in addition to local chromosomal architecture, CNVs are driven by the internal cellular or nuclear physiology of each cancer tissue.

Orbital TFE3-Rearranged Perivascular Epithelioid Cell Tumor: A Case Report and Review of the Literature

ABSTRACT

Perivascular epithelioid cell tumors (PEComas) are infrequent mesenchymal neoplasms that have particular histological and immunohistochemical features. Only a few cases have been described in the eye and orbit. This report presents a 28-year-old man who consulted for a painless left orbital mass. With the presumptive diagnosis of cavernous hemangioma, a surgical excisional biopsy was performed. Histopathological examination showed a well-delimited tumor composed of epithelioid cells with an eosinophilic cytoplasm and oval nucleus. The tumor cells stained diffusely for HMB-45 and transcription factor E3 (TFE3) and were focally positive for actin. There was no reactivity to S100 or desmin. Genetic testing revealed a TFE3 rearrangement. Based on these results, an extremely rare orbital TFE3-rearranged PEComa was diagnosed. Although no recurrence was seen at last follow-up, a review of the literature shows experience is limited regarding orbital PEComas and their malignant potential. Further research is needed to establish management guidelines, their association with the tuberous sclerosis complex, and the role of genetic mutations such as TFE3 rearrangement.

Assessing acute myeloid leukemia susceptibility in rearrangement-driven patients by DNA breakage at topoisomerase II and CCCTC-binding factor/cohesin binding sites

An initiating DNA double strand break (DSB) event precedes the formation of cancer-driven chromosomal abnormalities, such as gene rearrangements. Therefore, measuring DNA breaks at rearrangement-participating regions can provide a unique tool to identify and characterize susceptible individuals. Here, we developed a highly sensitive and low-input DNA break mapping method, the first of its kind for patient samples. We then measured genome-wide DNA breakage in normal cells of acute myeloid leukemia (AML) patients with KMT2A (previously MLL) rearrangements, compared to that of nonfusion AML individuals, as a means to evaluate individual susceptibility to gene rearrangements. DNA breakage at the KMT2A gene region was significantly greater in fusion-driven remission individuals, as compared to nonfusion individuals. Moreover, we identified select topoisomerase II (TOP2)-sensitive and CCCTC-binding factor (CTCF)/cohesin-binding sites with preferential DNA breakage in fusion-driven patients. Importantly, measuring DSBs at these sites, in addition to the KMT2A gene region, provided greater predictive power when assessing individual break susceptibility. We also demonstrated that low-dose etoposide exposure further elevated DNA breakage at these regions in fusion-driven AML patients, but not in nonfusion patients, indicating that these sites are preferentially sensitive to TOP2 activity in fusion-driven AML patients. These results support that mapping of DSBs in patients enables discovery of novel break-prone regions and monitoring of individuals susceptible to chromosomal abnormalities, and thus cancer. This will build the foundation for early detection of cancer-susceptible individuals, as well as those preferentially susceptible to therapy-related malignancies caused by treatment with TOP2 poisons.

Rearrangement-mediated cis-regulatory alterations in advanced patient tumors reveal interactions with therapy

The global impact of somatic structural variants (SVs) on gene regulation in advanced tumors with complex treatment histories has been mostly uncharacterized. Here, using whole-genome and RNA sequencing from 570 recurrent or metastatic tumors, we report the altered expression of hundreds of genes in association with nearby SV breakpoints, including oncogenes and G-protein-coupled receptor-related genes such as PLEKHG2. A significant fraction of genes with SV-expression associations correlate with worse patient survival in primary and advanced cancers, including SRD5A1. In many instances, SV-expression associations involve retrotransposons being translocated near genes. High overall SV burden is associated with treatment with DNA alkylating agents or taxanes and altered expression of metabolism-associated genes. SV-expression associations within tumors from topoisomerase I inhibitor-treated patients include chromatin-related genes. Within anthracycline-treated tumors, SV breakpoints near chromosome 1p genes include PDE4B. Patient treatment and history can help understand the widespread SV-mediated cis-regulatory alterations found in cancer.

Heuristics for Genome Rearrangement Distance With Replicated Genes

In comparative genomics, one goal is to find similarities between genomes of different organisms. Comparisons using genome features like genes, gene order, and regulatory sequences are carried out with this purpose in mind. Genome rearrangements are mutational events that affect large extensions of the genome. They are responsible for creating extant species with conserved genes in different positions across genomes. Close species - from an evolutionary point of view - tend to have the same set of genes or share most of them. When we consider gene order to compare two genomes, it is possible to use a parsimony criterion to estimate how close the species are. We are interested in the shortest sequence of genome rearrangements capable of transforming one genome into the other, which is named rearrangement distance. Reversal is one of the most studied genome rearrangements events. This event acts in a segment of the genome, inverting the position and the orientation of genes in it. Transposition is another widely studied event. This event swaps the position of two consecutive segments of the genome. When the genome has no gene repetition, a common approach is to map it as a permutation such that each element represents a conserved block. When genomes have replicated genes, this mapping is usually performed using strings. The number of replicas depends on the organisms being compared, but in many scenarios, it tends to be small. In this work, we study the rearrangement distance between genomes with replicated genes considering that the orientation of genes is unknown. We present four heuristics for the problem of genome rearrangement distance with replicated genes. We carry out experiments considering the exclusive use of the reversals or transpositions events, as well as the version in which both events are allowed. We developed a database of simulated genomes and compared our results with other algorithms from the literature. The experiments showed that our heuristics with more sophisticated rules presented a better performance than the known algorithms to estimate the evolutionary distance between genomes with replicated genes. In order to validate the application of our algorithms in real data, we construct a phylogenetic tree based on the distance provided by our algorithm and compare it with a know tree from the literature.

Sorting Signed Permutations by Inverse Tandem Duplication Random Losses

Gene order evolution of unichromosomal genomes, for example mitochondrial genomes, has been modelled mostly by four major types of genome rearrangements: inversions, transpositions, inverse transpositions, and tandem duplication random losses. Generalizing models that include all those rearrangements while admitting computational tractability are rare. In this paper, we study such a rearrangement model, namely the inverse tandem duplication random loss (iTDRL) model, where an iTDRL duplicates and inverts a continuous segment of a gene order followed by the random loss of one of the redundant copies of each gene. The iTDRL rearrangement has currently been proposed by several authors suggesting it to be a possible mechanisms of mitochondrial gene order evolution. We initiate the algorithmic study of this new model of genome rearrangement by proving that a shortest rearrangement scenario that transforms one given gene order into another given gene order can be obtained in quasilinear time. Furthermore, we show that the length of such a scenario, i.e., the minimum number of iTDRLs in the transformation, can be computed in linear time.

Sorting Permutations by Intergenic Operations

Genome Rearrangements are events that affect large stretches of genomes during evolution. Many mathematical models have been used to estimate the evolutionary distance between two genomes based on genome rearrangements. However, most of them focused on the (order of the) genes of a genome, disregarding other important elements in it. Recently, researchers have shown that considering regions between each pair of genes, called intergenic regions, can enhance distance estimation in realistic data. Two of the most studied genome rearrangements are the reversal, which inverts a sequence of genes, and the transposition, which occurs when two adjacent gene sequences swap their positions inside the genome. In this work, we study the transposition distance between two genomes, but we also consider intergenic regions, a problem we name Sorting by Intergenic Transpositions. We show that this problem is NP-hard and propose two approximation algorithms, with factors 3.5 and 2.5, considering two distinct definitions for the problem. We also investigate the signed reversal and transposition distance between two genomes considering their intergenic regions. This second problem is called Sorting by Signed Intergenic Reversals and Intergenic Transpositions. We show that this problem is NP-hard and develop two approximation algorithms, with factors 3 and 2.5. We check how these algorithms behave when assigning weights for genome rearrangements. Finally, we implemented all these algorithms and tested them on real and simulated data.

Computing the Inversion-Indel Distance

The inversion distance, that is the distance between two unichromosomal genomes with the same content allowing only inversions of DNA segments, can be exactly computed thanks to a pioneering approach of Hannenhalli and Pevzner from 1995. In 2000, El-Mabrouk extended the inversion model to perform the comparison of unichromosomal genomes with unequal contents, combining inversions with insertions and deletions (indels) of DNA segments, giving rise to the inversion-indel distance. However, only a heuristic was provided for its computation. In 2005, Yancopoulos, Attie and Friedberg started a new branch of research by introducing the generic double cut and join (DCJ) operation, that can represent several genome rearrangements (including inversions). In 2006, Bergeron, Mixtacki and Stoye showed that the DCJ distance can be computed in linear time with a very simple procedure. As a consequence, in 2010 we gave a linear-time algorithm to compute the DCJ-indel distance. This result allowed the inversion-indel model to be revisited from another angle. In 2013, we could show that, when the diagram that represents the relation between the two compared genomes has no bad components, the inversion-indel distance is equal to the DCJ-indel distance. In the present work we complete the study of the inversion-indel distance by giving the first algorithm to compute it exactly even in the presence of bad components.

Sorting Signed Permutations by Intergenic Reversals

Genome rearrangements are mutations affecting large portions of a genome, and a reversal is one of the most studied genome rearrangements in the literature through the Sorting by Reversals (SbR) problem. SbR is solvable in polynomial time on signed permutations (i.e., the gene orientation is known), and it is NP-hard on unsigned permutations. This problem (and many others considering genome rearrangements) models genome as a list of its genes in the order they appear, ignoring all other information present in the genome. Recent works claimed that the incorporation of the size of intergenic regions, i.e., sequences of nucleotides between genes, may result in better estimators for the real distance between genomes. Here we introduce the Sorting Signed Permutations by Intergenic Reversals problem, that sorts a signed permutation using reversals both on gene order and intergenic sizes. We show that this problem is NP-hard by a reduction from the 3-partition problem. Then, we propose a 2-approximation algorithm for it. Finally, we also incorporate intergenic indels (i.e., insertions or deletions of intergenic regions) to overcome a limitation of sorting by conservative events (such as reversals) and propose two approximation algorithms.

Radiation-induced gliomas represent H3-/IDH-wild type pediatric gliomas with recurrent PDGFRA amplification and loss of CDKN2A/B

Long-term complications such as radiation-induced second malignancies occur in a subset of patients following radiation-therapy, particularly relevant in pediatric patients due to the long follow-up period in case of survival. Radiation-induced gliomas (RIGs) have been reported in patients after treatment with cranial irradiation for various primary malignancies such as acute lymphoblastic leukemia (ALL) and medulloblastoma (MB). We perform comprehensive (epi-) genetic and expression profiling of RIGs arising after cranial irradiation for MB (n = 23) and ALL (n = 9). Our study reveals a unifying molecular signature for the majority of RIGs, with recurrent PDGFRA amplification and loss of CDKN2A/B and an absence of somatic hotspot mutations in genes encoding histone 3 variants or IDH1/2, uncovering diagnostic markers and potentially actionable targets.

Extensive genomic rearrangements mediated by repetitive sequences in plastomes of Medicago and its relatives

BACKGROUND

Although plastomes are highly conserved with respect to gene content and order in most photosynthetic angiosperms, extensive genomic rearrangements have been reported in Fabaceae, particularly within the inverted repeat lacking clade (IRLC) of Papilionoideae. Two hypotheses, i.e., the absence of the IR and the increased repeat content, have been proposed to affect the stability of plastomes. However, this is still unclear for the IRLC species. Here, we aimed to investigate the relationships between repeat content and the degree of genomic rearrangements in plastomes of Medicago and its relatives Trigonella and Melilotus, which are nested firmly within the IRLC.

RESULTS

We detected abundant repetitive elements and extensive genomic rearrangements in the 75 newly assembled plastomes of 20 species, including gene loss, intron loss and gain, pseudogenization, tRNA duplication, inversion, and a second independent IR gain (IR ~ 15 kb in Melilotus dentata) in addition to the previous first reported cases in Medicago minima. We also conducted comparative genomic analysis to evaluate plastome evolution. Our results indicated that the overall repeat content is positively correlated with the degree of genomic rearrangements. Some of the genomic rearrangements were found to be directly linked with repetitive sequences. Tandem repeated sequences have been detected in the three genes with accelerated substitution rates (i.e., accD, clpP, and ycf1) and their length variation could be explained by the insertions of tandem repeats. The repeat contents of the three localized hypermutation regions around these three genes with accelerated substitution rates are also significantly higher than that of the remaining plastome sequences.

CONCLUSIONS

Our results suggest that IR reemergence in the IRLC species does not ensure their plastome stability. Instead, repeat-mediated illegitimate recombination is the major mechanism leading to genome instability, a pattern in agreement with recent findings in other angiosperm lineages. The plastome data generated herein provide valuable genomic resources for further investigating the plastome evolution in legumes.

Chromoanagenesis from radiation-induced genome damage in Populus

Chromoanagenesis is a genomic catastrophe that results in chromosomal shattering and reassembly. These extreme single chromosome events were first identified in cancer, and have since been observed in other systems, but have so far only been formally documented in plants in the context of haploid induction crosses. The frequency, origins, consequences, and evolutionary impact of such major chromosomal remodeling in other situations remain obscure. Here, we demonstrate the occurrence of chromoanagenesis in poplar (Populus sp.) trees produced from gamma-irradiated pollen. Specifically, in this population of siblings carrying indel mutations, two individuals exhibited highly frequent copy number variation (CNV) clustered on a single chromosome, one of the hallmarks of chromoanagenesis. Using short-read sequencing, we confirmed the presence of clustered segmental rearrangement. Independently, we identified and validated novel DNA junctions and confirmed that they were clustered and corresponded to these rearrangements. Our reconstruction of the novel sequences suggests that the chromosomal segments have reorganized randomly to produce a novel rearranged chromosome but that two different mechanisms might be at play. Our results indicate that gamma irradiation can trigger chromoanagenesis, suggesting that this may also occur when natural or induced mutagens cause DNA breaks. We further demonstrate that such events can be tolerated in poplar, and even replicated clonally, providing an attractive system for more in-depth investigations of their consequences.

Plant breeders often use radiation treatment to produce variation, with the goal of identifying new varieties with superior traits. We studied a population of poplar trees produced by gamma irradiation of pollen, and asked what kind of DNA changes were associated with this variation. We found many changes, most often in the form of added (insertions) or removed (deletions) pieces of DNA. We also found two lines with much more drastic changes. In those lines, we observed massive reorganization. We characterized these two lines in detail and found that catastrophic pulverization and random reassembly only occurred on a single chromosome. Looking closely at how the pieces were put back together suggest that the rearrangements in these two lines may have resulted from two slightly different mechanisms. This type of rearrangement is commonly observed in human cancer cells, but has rarely been observed in plants. We demonstrated here that they can be induced by gamma irradiation, indicating this type of event might be more widespread than we expected. Characterizing such genome restructuring instances helps to understand how genome instability can remodel chromosomes and affect genome function.

R-loops and regulatory changes in chronologically ageing fission yeast cells drive non-random patterns of genome rearrangements

Aberrant repair of DNA double-strand breaks can recombine distant chromosomal breakpoints. Chromosomal rearrangements compromise genome function and are a hallmark of ageing. Rearrangements are challenging to detect in non-dividing cell populations, because they reflect individually rare, heterogeneous events. The genomic distribution of de novo rearrangements in non-dividing cells, and their dynamics during ageing, remain therefore poorly characterized. Studies of genomic instability during ageing have focussed on mitochondrial DNA, small genetic variants, or proliferating cells. To characterize genome rearrangements during cellular ageing in non-dividing cells, we interrogated a single diagnostic measure, DNA breakpoint junctions, using Schizosaccharomyces pombe as a model system. Aberrant DNA junctions that accumulated with age were associated with microhomology sequences and R-loops. Global hotspots for age-associated breakpoint formation were evident near telomeric genes and linked to remote breakpoints elsewhere in the genome, including the mitochondrial chromosome. Formation of breakpoint junctions at global hotspots was inhibited by the Sir2 histone deacetylase and might be triggered by an age-dependent de-repression of chromatin silencing. An unexpected mechanism of genomic instability may cause more local hotspots: age-associated reduction in an RNA-binding protein triggering R-loops at target loci. This result suggests that biological processes other than transcription or replication can drive genome rearrangements. Notably, we detected similar signatures of genome rearrangements that accumulated in old brain cells of humans. These findings provide insights into the unique patterns and possible mechanisms of genome rearrangements in non-dividing cells, which can be promoted by ageing-related changes in gene-regulatory proteins.

Concomitance of a novel RMDN2-ALK fusion and an EML4-ALK fusion in a lung adenocarcinoma

The anaplastic lymphoma kinase (ALK) fusions/rearrangements in non-small cell lung cancer (NSCLC) act as oncogenic driver mutations. ALK tyrosine kinase inhibitors have anti-tumor activities in ALK-positive NSCLC. Although the EML4-ALK fusion is common in NSCLC, concomitance of an additional ALK fusion together with an EML4-ALK fusion is not common. Here, we present a lung adenocarcinoma with two ALK fusions, a novel RMDN2-ALK fusion accompanied by an EML4-ALK fusion, detected by a targeted next generation sequencing assay. The genomic translocation breakpoints of the RMDN2-ALK fusion were mapped to intron 2 for RMDN2 and exon 15 for ALK, and EML4-ALK breakpoints were mapped to intron 13 for EML4 and intron 19 for ALK. ALK break-apart FISH detected multiple ALK rearrangements, a gene fusion panel (NanoString) test confirmed the EML4-ALK fusion, and RNA-sequencing revealed two ALK fusions. The RMDN2 gene locates at the short arm of chromosome 2 between ALK and EML4 genes. The intact ALK kinase domain fused to RMDN2. Genome-wide copy number variants were found in multiple chromosome arms and the short arm of chromosome 2, suggestive of complex rearrangements. Further detailed analyses of breakpoints and copy number variants may shed light on mechanisms of their formation and pathogenesis in lung malignancies.

Multiple endocrine neoplasia 2A with RET mutation p.Cys611Tyr: A case report

RATIONALE

Multiple endocrine neoplasia 2A (MEN2A) is a rare autosomal-dominant genetic syndrome, frequently misdiagnosed or neglected clinically, resulting in delayed therapy to patients.

PATIENT CONCERNS

A 47-year-old Chinese male patient underwent laparoscopic right adrenal tumorectomy, and postoperative pathology confirmed the tumor as pheochromocytoma (PHEO). He was readmitted to the department of endocrinology and metabolism due to constant increase in carcinoembryonic antigen (CEA) at 5 months after the operation.

DIAGNOSIS

The patient was confirmed with medullary thyroid carcinoma (MTC), multiple neck lymph node metastasis, and pituitary microadenoma. The p.Cys611Tyr (c.1832G>A, C611Y) mutation was detected. Therefore, he was diagnosed with MEN2A.

INTERVENTIONS

He underwent total thyroidectomy. The gene-sequencing analysis of his family was conducted, and the C611Y mutation was detected in his daughter.

OUTCOMES

The level of carcinoembryonic antigen decreased significantly after thyroidectomy in this patient. Long-term follow-up management was conducted. Elevated serum calcitonin and bilateral thyroid nodules were found in his 13-year-old daughter. Thus, MEN2A was highly suspected and she was suggested to undergo total thyroidectomy.

CONCLUSION

Patients with MEN2A should be screened regularly and managed by a multidisciplinary team.

Relaxed selection underlies genome erosion in socially parasitic ant species

Inquiline ants are highly specialized and obligate social parasites that infiltrate and exploit colonies of closely related species. They have evolved many times convergently, are often evolutionarily young lineages, and are almost invariably rare. Focusing on the leaf-cutting ant genus Acromyrmex, we compared genomes of three inquiline social parasites with their free-living, closely-related hosts. The social parasite genomes show distinct signatures of erosion compared to the host lineages, as a consequence of relaxed selective constraints on traits associated with cooperative ant colony life and of inquilines having very small effective population sizes. We find parallel gene losses, particularly in olfactory receptors, consistent with inquiline species having highly reduced social behavioral repertoires. Many of the genomic changes that we uncover resemble those observed in the genomes of obligate non-social parasites and intracellular endosymbionts that branched off into highly specialized, host-dependent niches.

Clinical utility of next generation sequencing to detect IGH/IL3 rearrangements [t(5;14)(q31.1;q32.1)] in B-lymphoblastic leukemia/lymphoma

The t(5;14)(q31.1;q32.1) associated with B-lymphoblastic leukemia/lymphoma (B-ALL/LBL) is a rare, recurrent genetic abnormality recognized as a distinct entity by the 2017 World Health Organization (WHO) classification. In these cases, the IGH enhancer region (14q32.1) is juxtaposed to the vicinity of the IL3 gene (5q31.1), resulting in increased production of interleukin-3 (IL3) and subsequently a characteristic reactive eosinophilia. B-ALL with t(5;14)(q31.1;q32.1) may have a low lymphoblast count that can complicate detection of t(5;14)(q31.1;q32.1) by conventional chromosome studies. We have identified four patients with IGH/IL3 rearrangements despite normal conventional chromosome studies in each case [one patient had a non-clonal t(5;14)(q31;q32) finding]. Fluorescence in situ hybridization utilizing a laboratory-developed IGH break-apart probe set identified IGH rearrangements in three of four cases, and a next generation sequencing (NGS) based assay, mate-pair sequencing (MPseq), was required to characterize the IGH/IL3 rearrangements in each case. Three patients demonstrated a balanced t(5;14)(q31.1;q32.1) while one patient had a cryptic insertion of the IL3 gene into the IGH region. These results demonstrate that NGS-based assays, such as MPseq, confer an advantage in the detection of IGH/IL3 rearrangements that are otherwise challenging to characterize by traditional cytogenetic methodologies.

Interspecies Chromosome Mapping in Caprimulgiformes, Piciformes, Suliformes, and Trogoniformes (Aves): Cytogenomic Insight into Microchromosome Organization and Karyotype Evolution in Birds

Interchromosomal rearrangements involving microchromosomes are rare events in birds. To date, they have been found mostly in Psittaciformes, Falconiformes, and Cuculiformes, although only a few orders have been analyzed. Hence, cytogenomic studies focusing on microchromosomes in species belonging to different bird orders are essential to shed more light on the avian chromosome and karyotype evolution. Based on this, we performed a comparative chromosome mapping for chicken microchromosomes 10 to 28 using interspecies BAC-based FISH hybridization in five species, representing four Neoaves orders (Caprimulgiformes, Piciformes, Suliformes, and Trogoniformes). Our results suggest that the ancestral microchromosomal syntenies are conserved in Pteroglossus inscriptus (Piciformes), Ramphastos tucanus tucanus (Piciformes), and Trogon surrucura surrucura (Trogoniformes). On the other hand, chromosome reorganization in Phalacrocorax brasilianus (Suliformes) and Hydropsalis torquata (Caprimulgiformes) included fusions involving both macro- and microchromosomes. Fissions in macrochromosomes were observed in P. brasilianus and H. torquata. Relevant hypothetical Neognathae and Neoaves ancestral karyotypes were reconstructed to trace these rearrangements. We found no interchromosomal rearrangement involving microchromosomes to be shared between avian orders where rearrangements were detected. Our findings suggest that convergent evolution involving microchromosomal change is a rare event in birds and may be appropriate in cytotaxonomic inferences in orders where these rearrangements occurred.

Computing the Rearrangement Distance of Natural Genomes

The computation of genomic distances has been a very active field of computational comparative genomics over the past 25 years. Substantial results include the polynomial-time computability of the inversion distance by Hannenhalli and Pevzner in 1995 and the introduction of the double cut and join distance by Yancopoulos et al. in 2005. Both results, however, rely on the assumption that the genomes under comparison contain the same set of unique markers (syntenic genomic regions, sometimes also referred to as genes). In 2015, Shao et al. relax this condition by allowing for duplicate markers in the analysis. This generalized version of the genomic distance problem is NP-hard, and they give an integer linear programming (ILP) solution that is efficient enough to be applied to real-world datasets. A restriction of their approach is that it can be applied only to balanced genomes that have equal numbers of duplicates of any marker. Therefore, it still needs a delicate preprocessing of the input data in which excessive copies of unbalanced markers have to be removed. In this article, we present an algorithm solving the genomic distance problem for natural genomes, in which any marker may occur an arbitrary number of times. Our method is based on a new graph data structure, the multi-relational diagram, that allows an elegant extension of the ILP by Shao et al. to count runs of markers that are under- or over-represented in one genome with respect to the other and need to be inserted or deleted, respectively. With this extension, previous restrictions on the genome configurations are lifted, for the first time enabling an uncompromising rearrangement analysis. Any marker sequence can directly be used for the distance calculation. The evaluation of our approach shows that it can be used to analyze genomes with up to a few 10,000 markers, which we demonstrate on simulated and real data.

Comprehensive analysis of cancer breakpoints reveals signatures of genetic and epigenetic contribution to cancer genome rearrangements

Understanding mechanisms of cancer breakpoint mutagenesis is a difficult task and predictive models of cancer breakpoint formation have to this time failed to achieve even moderate predictive power. Here we take advantage of a machine learning approach that can gather important features from big data and quantify contribution of different factors. We performed comprehensive analysis of almost 630,000 cancer breakpoints and quantified the contribution of genomic and epigenomic features-non-B DNA structures, chromatin organization, transcription factor binding sites and epigenetic markers. The results showed that transcription and formation of non-B DNA structures are two major processes responsible for cancer genome fragility. Epigenetic factors, such as chromatin organization in TADs, open/closed regions, DNA methylation, histone marks are less informative but do make their contribution. As a general trend, individual features inside the groups show a relatively high contribution of G-quadruplexes and repeats and CTCF, GABPA, RXRA, SP1, MAX and NR2F2 transcription factors. Overall, the cancer breakpoint landscape can be represented by well-predicted hotspots and poorly predicted individual breakpoints scattered across genomes. We demonstrated that hotspot mutagenesis has genomic and epigenomic factors, and not all individual cancer breakpoints are just random noise but have a definite mutation signature. Besides we found a long-range action of some features on breakpoint mutagenesis. Combining omics data, cancer-specific individual feature importance and adding the distant to local features, predictive models for cancer breakpoint formation achieved 70-90% ROC AUC for different cancer types; however precision remained low at 2% and the recall did not exceed 50%. On the one hand, the power of models strongly correlates with the size of available cancer breakpoint and epigenomic data, and on the other hand finding strong determinants of cancer breakpoint formation still remains a challenge. The strength of predictive signals of each group and of each feature inside a group can be converted into cancer-specific breakpoint mutation signatures. Overall our results add to the understanding of cancer genome rearrangement processes.

Challenges and Opportunities in NUT Carcinoma Research

NUT carcinoma (NC) is a type of aggressive cancer driven by chromosome translocations. Fusion genes between a DNA-binding protein, such as bromodomain and extraterminal domain (BET) proteins, and the testis-specific protein NUTM1 generated by these translocations drive the formation of NC. NC can develop in very young children without significant accumulation of somatic mutations, presenting a relatively clean model to study the genetic etiology of oncogenesis. However, after 20 years of research, a few challenging questions still remain for understanding the mechanism and developing therapeutics for NC. In this short review, we first briefly summarize the current knowledge regarding the molecular mechanism and targeted therapy development of NC. We then raise three challenging questions: (1) What is the cell of origin of NC? (2) How does the germline analogous epigenetic reprogramming process driven by the BET-NUTM1 fusion proteins cause NC? and (3) How will BET-NUTM1 targeted therapies be developed? We propose that with the unprecedented technological advancements in genome editing, animal models, stem cell biology, organoids, and chemical biology, we have unique opportunities to address these challenges.

An efficient Oligo-FISH painting system for revealing chromosome rearrangements and polyploidization in Triticeae

A chromosome-specific painting technique has been developed which combines the most recent approaches of the companion disciplines of molecular cytogenetics and genome research. We developed seven oligonucleotide (oligo) pools derivd from single-copy sequences on chromosomes 1 to 7 of barley (Hordeum vulgare L.) and corresponding collinear regions of wheat (Triticum aestivum L.). The seven groups of pooled oligos comprised between 10 986 and 12 496 45-bp monomers, and these then produced stable fluorescence in situ hybridization (FISH) signals on chromosomes of each linkage group of wheat and barley. The pooled oligo probes were applied to high-throughput karyotyping of the chromosomes of other Triticeae species in the genera Secale, Aegilops, Thinopyrum, and Dasypyrum, and the study also extended to some wheat-alien amphiploids and derived lines. We demonstrated that a complete set of whole-chromosome oligo painting probes facilitated the study of inter-species chromosome homologous relationships and visualized non-homologous chromosomal rearrangements in Triticeae species and some wheat-alien species derivatives. When combined with other non-denaturing FISH procedures using tandem-repeat oligos, the newly developed oligo painting techniques provide an efficient tool for the study of chromosome structure, organization, and evolution among any wild Triticeae species with non-sequenced genomes.

Burkitt lymphoma in the modern era: real-world outcomes and prognostication across 30 US cancer centers

We examined adults with untreated Burkitt lymphoma (BL) from 2009 to 2018 across 30 US cancer centers. Factors associated with progression-free survival (PFS) and overall survival (OS) were evaluated in univariate and multivariate Cox models. Among 641 BL patients, baseline features included the following: median age, 47 years; HIV+, 22%; Eastern Cooperative Oncology Group (ECOG) performance status (PS) 2 to 4, 23%; >1 extranodal site, 43%; advanced stage, 78%; and central nervous system (CNS) involvement, 19%. Treatment-related mortality was 10%, with most common causes being sepsis, gastrointestinal bleed/perforation, and respiratory failure. With 45-month median follow-up, 3-year PFS and OS rates were 64% and 70%, respectively, without differences by HIV status. Survival was better for patients who received rituximab vs not (3-year PFS, 67% vs 38%; OS, 72% vs 44%; P < .001) and without difference based on setting of administration (ie, inpatient vs outpatient). Outcomes were also improved at an academic vs community cancer center (3-year PFS, 67% vs 46%, P = .006; OS, 72% vs 53%, P = .01). In multivariate models, age ≥ 40 years (PFS, hazard ratio [HR] = 1.70, P = .001; OS, HR = 2.09, P < .001), ECOG PS 2 to 4 (PFS, HR = 1.60, P < .001; OS, HR = 1.74, P = .003), lactate dehydrogenase > 3× normal (PFS, HR = 1.83, P < .001; OS, HR = 1.63, P = .009), and CNS involvement (PFS, HR = 1.52, P = .017; OS, HR = 1.67, P = .014) predicted inferior survival. Furthermore, survival varied based on number of factors present (0, 1, 2 to 4 factors) yielding 3-year PFS rates of 91%, 73%, and 50%, respectively; and 3-year OS rates of 95%, 77%, and 56%, respectively. Collectively, outcomes for adult BL in this real-world analysis appeared more modest compared with results of clinical trials and smaller series. In addition, clinical prognostic factors at diagnosis identified patients with divergent survival rates.

Comprehensive Overview of Gene Rearrangements in Childhood T-Cell Acute Lymphoblastic Leukaemia

Acute lymphoblastic leukaemia (ALL) is a relevant form of childhood neoplasm, as it accounts for over 80% of all leukaemia cases. T-cell ALL constitutes a genetically heterogeneous cancer derived from T-lymphoid progenitors. The diagnosis of T-ALL is based on morphologic, immunophenotypic, cytogenetic, and molecular features, thus the results are used for patient stratification. Due to the expression of surface and intracellular antigens, several subtypes of T-ALL can be distinguished. Although the aetiology of T-ALL remains unclear, a wide spectrum of rearrangements and mutations affecting crucial signalling pathways has been described so far. Due to intensive chemotherapy regimens and supportive care, overall cure rates of more than 80% in paediatric T-ALL patients have been accomplished. However, improved knowledge of the mechanisms of relapse, drug resistance, and determination of risk factors are crucial for patients in the high-risk group. Even though some residual disease studies have allowed the optimization of therapy, the identification of novel diagnostic and prognostic markers is required to individualize therapy. The following review summarizes our current knowledge about genetic abnormalities in paediatric patients with T-ALL. As molecular biology techniques provide insights into the biology of cancer, our study focuses on new potential therapeutic targets and predictive factors which may improve the outcome of young patients with T-ALL.

Chromosome 2q14.3 microdeletion encompassing CNTNAP5 gene in a patient carrying a complex chromosomal rearrangement

We report a patient with loss of chromosome region 2q14.3 encompassing exon 1 of the gene CNTNAP5. The deletion occurred in association with a de novo complex chromosomal rearrangement, characterized by routine G-banding, fluorescence in situ hybridization and microarray analysis. The presented patient's phenotype is dominated by severe early childhood weight gain, severe speech delay and behavioural problems. To our knowledge, a few similar patients have been reported previously. CNTNAP5 is a member of the neurexin gene family and is associated with autism spectrum disorder and potentially other behavioural and neurodevelopmental disorders. Recent data point to its possible role in obesity and/or metabolism. The phenotype of the herein presented pediatric patient corroborates CNTNAP5's pathogenic role in human disease.

Distinctive evolutionary pattern of organelle genomes linked to the nuclear genome in Selaginellaceae

Plastids and mitochondria are endosymbiotic organelles that store genetic information. The genomes of these organelles generally exhibit contrasting patterns regarding genome architecture and genetic content. However, they have similar genetic features in Selaginellaceae, and little is known about what causes parallel evolution. Here, we document the multipartite plastid genomes (plastomes) and the highly divergent mitochondrial genomes (mitogenomes) from spikemoss obtained by combining short- and long-reads. The 188-kb multipartite plastome has three ribosomal operon copies in the master genomic conformation, creating the alternative subgenomic conformation composed of 110- and 78-kb subgenomes. The long-read data indicated that the two different genomic conformations were present in almost equal proportions in the plastomes of Selaginella nipponica. The mitogenome of S. nipponica was assembled into 27 contigs with a total size of 110 kb. All contigs contained directly arranged repeats at both ends, which introduced multiple conformations. Our results showed that plastomes and mitogenomes share high tRNA losses, GC-biased nucleotides, elevated substitution rates and complicated organization. The exploration of nuclear-encoded organelle DNA replication, recombination and repair proteins indicated that, several single-targeted proteins, particularly plastid-targeted recombinase A1, have been lost in Selaginellaceae; conversely, the dual-targeted proteins remain intact. According to the reported function of recombinase A1, we propose that the plastomes of spikemoss often fail to pair homologous sequences during recombination, and the dual-targeted proteins play a key role in the convergent genetic features of plastomes and mitogenomes. Our results provide a distinctive evolutionary pattern of the organelle genomes in Selaginellaceae and evidence of their convergent evolution.

Clinical activity of brigatinib in ROS1-rearranged non-small cell lung cancer

BACKGROUND

Brigatinib is a potent ROS1 inhibitor. The existing data on its clinical activity in ROS1-rearranged non-small cell lung cancer (NSCLC) are limited to four cases.

METHODS

Six patients with ROS1-rearranged advanced NSCLC treated with brigatinib were identified through search of the internal databases of four participating cancer centers. Four additional patients were selected by PubMed and Google Scholar search. The objective response rate (ORR), progression-free survival (PFS) (RECIST v.1.1), duration of treatment (DOT), and safety were assessed.

RESULTS

Of eight patients evaluable for response assessment (crizotinib naive-1, crizotinib resistant -7), three patients demonstrated a partial response (ORR-37%). One crizotinib-naive patient had an ongoing response at 21.6 months. Of seven crizotinib-resistant patients, two patients demonstrated a partial response (ORR-29%), and one patient (14%) had stable disease. PFS, available in four crizotinib-resistant patients, was 7.6 + , 2.9, 2.0, and 0.4 months. In crizotinib-resistant patients, DOT was 9.7 + , 7.7 + , 7.6 + , 4.0, 2.0, 1.1, 0.4 months, and was not reported in two patients. Genomic profiling in one responder revealed no ROS1 alteration, suggesting that the response was attributable to "off-target" brigatinib activity. In two patients with progressive disease, genomic profiling demonstrated a cMET exon 14 mutation + KRAS G12A mutation in one case, and a persisting ROS1-CD74 fusion + TP53 K139N, FGFR2 E250G, ATM G2695D, and NF1 R2258Q mutations in the other. No grade 3-5 toxicity was observed.

CONCLUSION

Brigatinib demonstrated modest activity in crizotinib-resistant ROS1-rearranged NSCLC. Its intracranial and systemic activity should be assessed in correlation with the underlying molecular mechanism of crizotinib resistance.

Resolving Breakpoints of Chromosomal Rearrangements at the Nucleotide Level Using Sanger Sequencing

Novel cytogenetic tools are increasingly based on genome sequencing for detecting chromosomal abnormalities. Different sequence-based techniques optimized for diagnosis of structural variants can be useful for narrowing down the localization of breakpoints of chromosomal abnormalities, but do not offer nucleotide resolution of breakpoints for proper interpretation of gene disruption. This protocol presents the characterization of structural variants at nucleotide resolution using Sanger sequencing after low-pass large-insert genome sequencing or other long-molecule methods. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Primer design for junction amplification at translocations and inversions Basic Protocol 2: Amplification of derivative chromosomes using a long-range polymerase Alternate Protocol: Amplification of derivative chromosomes using a hot-start polymerase Basic Protocol 3: Preparation of DNA for Sanger sequencing Basic Protocol 4: Interpretation and reporting of breakpoints based on Sanger sequencing.

Evaluation of Single-Molecule Sequencing Technologies for Structural Variant Detection in Two Swedish Human Genomes

Long-read single molecule sequencing is increasingly used in human genomics research, as it allows to accurately detect large-scale DNA rearrangements such as structural variations (SVs) at high resolution. However, few studies have evaluated the performance of different single molecule sequencing platforms for SV detection in human samples. Here we performed Oxford Nanopore Technologies (ONT) whole-genome sequencing of two Swedish human samples (average 32× coverage) and compared the results to previously generated Pacific Biosciences (PacBio) data for the same individuals (average 66× coverage). Our analysis inferred an average of 17k and 23k SVs from the ONT and PacBio data, respectively, with a majority of them overlapping with an available multi-platform SV dataset. When comparing the SV calls in the two Swedish individuals, we find a higher concordance between ONT and PacBio SVs detected in the same individual as compared to SVs detected by the same technology in different individuals. Downsampling of PacBio reads, performed to obtain similar coverage levels for all datasets, resulted in 17k SVs per individual and improved overlap with the ONT SVs. Our results suggest that ONT and PacBio have a similar performance for SV detection in human whole genome sequencing data, and that both technologies are feasible for population-scale studies.

Genome-wide mapping of spontaneous genetic alterations in diploid yeast cells

Genomic alterations including single-base mutations, deletions and duplications, translocations, mitotic recombination events, and chromosome aneuploidy generate genetic diversity. We examined the rates of all of these genetic changes in a diploid strain of Saccharomyces cerevisiae by whole-genome sequencing of many independent isolates (n = 93) subcloned about 100 times in unstressed growth conditions. The most common alterations were point mutations and small (<100 bp) insertion/deletions (n = 1,337) and mitotic recombination events (n = 1,215). The diploid cells of most eukaryotes are heterozygous for many single-nucleotide polymorphisms (SNPs). During mitotic cell divisions, recombination can produce derivatives of these cells that have become homozygous for the polymorphisms, termed loss-of-heterozygosity (LOH) events. LOH events can change the phenotype of the cells and contribute to tumor formation in humans. We observed two types of LOH events: interstitial events (conversions) resulting in a short LOH tract (usually less than 15 kb) and terminal events (mostly cross-overs) in which the LOH tract extends to the end of the chromosome. These two types of LOH events had different distributions, suggesting that they may have initiated by different mechanisms. Based on our results, we present a method of calculating the probability of an LOH event for individual SNPs located throughout the genome. We also identified several hotspots for chromosomal rearrangements (large deletions and duplications). Our results provide insights into the relative importance of different types of genetic alterations produced during vegetative growth.

Allele-Specific Chromosome Removal after Cas9 Cleavage in Human Embryos

Correction of disease-causing mutations in human embryos holds the potential to reduce the burden of inherited genetic disorders and improve fertility treatments for couples with disease-causing mutations in lieu of embryo selection. Here, we evaluate repair outcomes of a Cas9-induced double-strand break (DSB) introduced on the paternal chromosome at the EYS locus, which carries a frameshift mutation causing blindness. We show that the most common repair outcome is microhomology-mediated end joining, which occurs during the first cell cycle in the zygote, leading to embryos with non-mosaic restoration of the reading frame. Notably, about half of the breaks remain unrepaired, resulting in an undetectable paternal allele and, after mitosis, loss of one or both chromosomal arms. Correspondingly, Cas9 off-target cleavage results in chromosomal losses and hemizygous indels because of cleavage of both alleles. These results demonstrate the ability to manipulate chromosome content and reveal significant challenges for mutation correction in human embryos.

The complete mitochondrial genome of medicinal fungus Taiwanofungus camphoratus reveals gene rearrangements and intron dynamics of Polyporales

Taiwanofungus camphoratus is a highly valued medicinal mushroom that is endemic to Taiwan, China. In the present study, the mitogenome of T. camphoratus was assembled and compared with other published Polyporales mitogenomes. The T. camphoratus mitogenome was composed of circular DNA molecules, with a total size of 114,922 bp. Genome collinearity analysis revealed large-scale gene rearrangements between the mitogenomes of Polyporales, and T. camphoratus contained a unique gene order. The number and classes of introns were highly variable in 12 Polyporales species we examined, which proved that numerous intron loss or gain events occurred in the evolution of Polyporales. The Ka/Ks values for most core protein coding genes in Polyporales species were less than 1, indicating that these genes were subject to purifying selection. However, the rps3 gene was found under positive or relaxed selection between some Polyporales species. Phylogenetic analysis based on the combined mitochondrial gene set obtained a well-supported topology, and T. camphoratus was identified as a sister species to Laetiporus sulphureus. This study served as the first report on the mitogenome in the Taiwanofungus genus, which will provide a basis for understanding the phylogeny and evolution of this important fungus.

Extended TAQing system for large-scale plant genome reorganization

We previously developed a large-scale genome restructuring technology called the TAQing system. It can induce genomic rearrangements by introducing transient and conditional formation of DNA double-strand breaks (DSBs) via heat activation of a restriction enzyme TaqI, which can cleave DNA at 5'-TCGA-3' sequences in the genome at higher temperatures (37-42°C). Such heat treatment sometimes confers lethal damage in certain plant species and TaqI cannot induce rearrangements in AT-rich regions. To overcome such problems we developed an extended TAQing (Ex-TAQing) system, which enables the use of a wider range of restriction enzymes active at standard plant-growing temperatures. We established the Ex-TAQing system using MseI that can efficiently cleave DNA at room temperature (at temperatures ranging from 22 to 25°C) and the 5'-TTAA-3' sequence which is highly abundant in the Arabidopsis genome. A synthetic intron-spanning MseI gene, which was placed downstream of a heat-shock-inducible promoter, was conditionally expressed upon milder heat treatment (33°C) to generate DSBs in Arabidopsis chromosomes. Genome resequencing revealed various types of genomic rearrangements, including copy number variations, translocation and direct end-joining at MseI cleavage sites. The Ex-TAQing system could induce large-scale rearrangements in diploids more frequently (17.4%, n = 23) than the standard TAQing system. The application of this system to tetraploids generated several strains with chromosomal rearrangements associated with beneficial phenotypes, such as high salinity stress tolerance and hypersensitivity to abscisic acid. We have developed the Ex-TAQing system, allowing more diverse patterns of genomic rearrangements, by employing various types of endonucleases and have opened a way to expand the capacity for artificial genome reorganization.