Intrachromosomal Rearrangements Research Articles

Evolutionary patterns and functional effects of 3D chromatin structures in butterflies with extensive genome rearrangements

Chromosome rearrangements may distort 3D chromatin architectures and thus change gene regulation, yet how 3D chromatin structures evolve in insects is largely unknown. Here, we obtain chromosome-level genomes for four butterfly species, Graphium cloanthus, Graphium sarpedon, Graphium eurypylus with 2n = 30, 40, and 60, respectively, and Papilio bianor with 2n = 60. Together with large-scale Hi-C data, we find that inter-chromosome rearrangements very rarely disrupted the pre-existing 3D chromatin structure of ancestral chromosomes. However, some intra-chromosome rearrangements changed 3D chromatin structures compared to the ancestral configuration. We find that new TADs and subTADs have emerged across the rearrangement sites where their adjacent compartments exhibit uniform types. Two intra-chromosome rearrangements altered Rel and lft regulation, potentially contributing to wing patterning differentiation and host plant choice. Notably, butterflies exhibited chromatin loops between Hox gene cluster ANT-C and BX-C, unlike Drosophila. Our CRISPR-Cas9 experiments in butterflies confirm that knocking out the CTCF binding site of the loops in BX-C affected the phenotypes regulated by Antp in ANT-C, resulting in legless larva. Our results reveal evolutionary patterns of insect 3D chromatin structures and provide evidence that 3D chromatin structure changes can play important roles in the evolution of traits.

Abstract 4334: Insight into the structure of breakage fusion bridge events and chromothripsis in cancer cell lines using long-read sequencing

Abstract Cervical cancer is caused by human papillomavirus (HPV) infection and is the most common cause of cancer death in low-resource countries. We characterized 18 cervical and 4 head and neck cancer cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, and cancer driver mutations. This cell line panel represents the major molecular subtypes of cervical cancer. Structural variation analysis revealed recurrent chromosomal alteration and telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alteration in cancer that has not been studied using long-read sequencing. Analysis of the inversion sites revealed staggered ends consistent with exonuclease deletion of the DNA after breakage. Some BFB events are complex, involving interchromosomal or intrachromosomal insertions or rearrangements. In total, there were three different types of BFB events: Type I, II, and III. Type I BFB events have staggered ends but no additional insertion at the fusion site. Type II BFB events have an insertion of DNA from another chromosome at the junction, and Type III has an insertion of sequences from the same chromosome as the BFB. BFB cycles are thought to resolve with the addition of a telomere. However, a search for telomeric sequences did not identify any on chromosomes impacted by BFB. In many cancers, BFB events co-occur with chromothripsis; however, in our cervical cancers, only 1 out of 13 had chromosomal shattering near the BFB event, indicating that the mechanism of chromosome rearrangement may be different in cervical cancers. Five cell lines have a Chr11q BFB event, with YAP1/BIRC2/BIRC3 gene amplification. Analysis of 911 publicly available cervical cancer exomes showed that YAP1 amplification is associated with a 10-year earlier age of diagnosis of cervical cancer (x2, p<0.0001) and is three times more common in African Americans (x2=12, p=0.003). Therefore, cervical cancer patients with YAP1/BIRC2/BIRC3 amplification, especially those of African American ancestry, might benefit from anti-YAP1/BIRC2/BIRC3 therapies. Our panel provides models for all the major HPV types and subtypes of cervical cancer that will aid in developing new therapeutic approaches. In summary, we uncovered novel insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing. We are expanding long-read sequencing to other cancer cell lines to gain knowledge of events associated with BFB. Citation Format: Isabel Rodriguez, Nicole Rossi, Ayse Gokce Kesus, Yi Xie, Tanveer Ahmad, Asher Bryant, Hong Lou, Jesica Godinez Paredes, Rose Milano, Mikhail Kolmogorov, Sonam Tulsyan, Jazmyn Bess, Vera Mukhina, Daria Gaykalova, Laksh Malik, Kimberley J. Billingsley, Cornelis Blauwendraat, Lisa Mirabello, Michael Dean. Insight into the structure of breakage fusion bridge events and chromothripsis in cancer cell lines using long-read sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4334.

Understanding the chromosomal evolution in cuckoos (Aves, Cuculiformes): a journey through unusual rearrangements.

The Cuculiformes are a family of over 150 species that live in a range of habitats, such as forests, savannas, and deserts. Here, bacterial artificial chromosome (BAC) probes (75 from chicken and 14 from zebra finch macrochromosomes 1-10 +ZW and for microchromosomes 11-28 (except 16)) were used to investigate chromosome homologies between chicken and the squirrel cuckoo (Piaya cayana). In addition, repetitive DNA probes were applied to characterize the chromosome organization and to explore the role of these sequences in the karyotype evolution of P. cayana. We also applied BAC probes for chicken chromosome 17 and Z to the guira cuckoo (Guira guira) to test whether this species has an unusual Robertsonian translocation between a microchromosome and the Z chromosome, recently described in the smooth-billed ani (Crotophaga ani). Our results revealed extensive chromosome reorganization with inter- and intrachromosomal rearrangements in P. cayana, including a conspicuous chromosome size and heterochromatin polymorphism on chromosome pair 20. Furthermore, we confirmed that the Z-autosome Robertsonian translocation found in C. ani is also found in G. guira, not P. cayana. These findings suggest that this translocation occurred prior to the divergence between C. ani and G. guira, but after the divergence with P. cayana.

A Bird's-Eye View of Chromosomic Evolution in the Class Aves.

Birds (Aves) are the most speciose of terrestrial vertebrates, displaying Class-specific characteristics yet incredible external phenotypic diversity. Critical to agriculture and as model organisms, birds have adapted to many habitats. The only extant examples of dinosaurs, birds emerged ~150 mya and >10% are currently threatened with extinction. This review is a comprehensive overview of avian genome ("chromosomic") organization research based mostly on chromosome painting and BAC-based studies. We discuss traditional and contemporary tools for reliably generating chromosome-level assemblies and analyzing multiple species at a higher resolution and wider phylogenetic distance than previously possible. These results permit more detailed investigations into inter- and intrachromosomal rearrangements, providing unique insights into evolution and speciation mechanisms. The 'signature' avian karyotype likely arose ~250 mya and remained largely unchanged in most groups including extinct dinosaurs. Exceptions include Psittaciformes, Falconiformes, Caprimulgiformes, Cuculiformes, Suliformes, occasional Passeriformes, Ciconiiformes, and Pelecaniformes. The reasons for this remarkable conservation may be the greater diploid chromosome number generating variation (the driver of natural selection) through a greater possible combination of gametes and/or an increase in recombination rate. A deeper understanding of avian genomic structure permits the exploration of fundamental biological questions pertaining to the role of evolutionary breakpoint regions and homologous synteny blocks.

Insights into the mechanisms and structure of breakage-fusion-bridge cycles in cervical cancer using long-read sequencing

Cervical cancer is caused by human papillomavirus (HPV) infection, has few approved targeted therapeutics, and is the most common cause of cancer death in low-resource countries. We characterized 19 cervical and four head and neck cancer cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, chromosomal alterations, and cancer driver mutations. Structural variation analysis revealed telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alterations in cancer, and our study applies long-read sequencing to this important chromosomal rearrangement type. Analysis of the inversion sites revealed staggered ends consistent with exonuclease digestion of the DNA after breakage. Some BFB events are complex, involving inter- or intra-chromosomal insertions or rearrangements. None of the BFB breakpoints had telomere sequences added to resolve the dicentric chromosomes, and only one BFB breakpoint showed chromothripsis. Five cell lines have a chromosomal region 11q BFB event, with YAP1-BIRC3-BIRC2 amplification. Indeed, YAP1 amplification is associated with a 10-year-earlier age of diagnosis of cervical cancer and is three times more common in African American women. This suggests that individuals with cervical cancer and YAP1-BIRC3-BIRC2 amplification, especially those of African ancestry, might benefit from targeted therapy. In summary, we uncovered valuable insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing.

Monitoring Genomic Structural Rearrangements Resulting from Gene Editing.

The cytogenomics-based methodology of directional genomic hybridization (dGH) enables the detection and quantification of a more comprehensive spectrum of genomic structural variants than any other approach currently available, and importantly, does so on a single-cell basis. Thus, dGH is well-suited for testing and/or validating new advancements in CRISPR-Cas9 gene editing systems. In addition to aberrations detected by traditional cytogenetic approaches, the strand specificity of dGH facilitates detection of otherwise cryptic intra-chromosomal rearrangements, specifically small inversions. As such, dGH represents a powerful, high-resolution approach for the quantitative monitoring of potentially detrimental genomic structural rearrangements resulting from exposure to agents that induce DNA double-strand breaks (DSBs), including restriction endonucleases and ionizing radiations. For intentional genome editing strategies, it is critical that any undesired effects of DSBs induced either by the editing system itself or by mis-repair with other endogenous DSBs are recognized and minimized. In this paper, we discuss the application of dGH for assessing gene editing-associated structural variants and the potential heterogeneity of such rearrangements among cells within an edited population, highlighting its relevance to personalized medicine strategies.

The 3D Chromatin Landscape of B-Cell Precursor Childhood Acute Lymphoblastic Leukemia

Whilst the molecular pathogenesis of childhood B-cell precursor (BCP) acute lymphoblastic leukemia (ALL) has been studied extensively, its 3D chromatin landscape remains poorly explored. Genome-wide chromosome conformation capture methods have provided the tools to investigate the different units of chromatin organization, such as transcriptionally active (A) and inactive (B) compartments, topologically associating domains (TADs), and fine-scale chromatin loops and enhancer-promoter interactions. The aim of this study is to elucidate the chromatin architecture and topological gene regulation in childhood BCP ALL. To date, 29 primary patient samples were included, comprising the high hyperdiploid (HeH) (n=11), ETV6:: RUNX1-positive (n=8), BCR:: ABL1-positive (n=2), TCF3:: PBX1-positive (n=2), DUX4-rearranged (n=2), intrachromosomal amplification of chromosome 21 (iAMP21) (n=1), KMT2A-rearranged (n=1), near-haploid (n=1) and near-triploid (n=1) genetic subtypes. Leukemic blast cells obtained at diagnosis were analyzed using Micro-C, a high-resolution variation of Hi-C (average number of total reads = 1.4 billion, highest resolution = 5 kb) combined with pair-end sequencing. Chromatin contact heatmaps were generated for each case using Juicer and Cooler. A/B compartments were identified using FanC at 500 kb resolution and visualized in the software IGV, while TAD calling was performed by Juicer, Domaincaller and Insulation Score. Differential chromatin interaction loop calling was made using Pareidolia and Mustache, and structural variant (SV) calling was carried out by EagleC. Preliminary principal component analysis of the 29 cases, based on the first two eigenvectors of the contact matrix (500 kb resolution), showed that HeH, ETV6:: RUNX1-positive, TCF3:: PBX1-positive and DUX4-rearranged cases each clustered based on their chromatin 3D organization. Furthermore, the A/B compartments of 11 HeH and 8 ETV6:: RUNX1-positive cases were analyzed at 500 kb resolution and compartment shifts among the two subtypes were annotated. A total of 390 shifts were detected, where activating shifts (from B to A compartment) happened more often in HeH (263 shifts) than in ETV6:: RUNX1-positive cases (127 shifts). Analysis of TAD boundary strength at 25 kb resolution revealed that HeH cases displayed significantly weaker boundaries compared to ETV6:: RUNX1-positive cases. TAD boundary strength showed no bias towards the frequently gained or non-gained chromosomes in HeH ALL. By merging individual heatmaps of all HeH and ETV6:: RUNX1-positive cases using Cooler, we created subtype-specific profiles and compared the intensity of chromatin interactions between the two genetic subtypes. Chromatin interaction intensity analysis was then combined with previously published RNA-sequencing data to identify transcriptional dysregulation events that could be associated with chromatin interaction changes. Preliminary results show that there was a chromatin loop missing close to the well-known leukemia-related gene IKZF1 in HeH compared to ETV6:: RUNX1-positive cases; this gene also showed lower expression in the RNA-sequencing data. FLT3 was associated with weakened chromatin interactions and down-regulated in ETV6:: RUNX1-positive cases compared to HeH, in agreement with its known high expression in HeH. Finally, we performed screening of SVs using EagleC and Micro-C heatmaps in HeH and ETV6:: RUNX1-positive samples. Out of the 19 included cases, previous whole-genome sequencing (WGS) data were available for 16. We detected 75 SVs, of which 50 were intrachromosomal rearrangements and 25 were translocations. Micro-C heatmaps allowed visual detection of SVs smaller than 1 Mb and permitted identification of the type of SVs. WGS detected 61% of the SVs found in the HeH samples with Micro-C and 51% of those in the ETV6:: RUNX1-positive cases. In summary, we present the first high-resolution genome-wide map of chromatin 3D organization in pediatric ALL. Our results indicate that different subtypes of childhood BCP ALL have distinct 3D chromatin landscapes and that abnormal chromatin architectures affect the regulation of leukemia-related genes.

Unraveling the Unusual Subgenomic Organization in the Neopolyploid Free-Living Flatworm Macrostomum lignano.

Whole genome duplication (WGD) is an evolutionary event resulting in a redundancy of genetic material. Different mechanisms of WGD, allo- or autopolyploidization, lead to distinct evolutionary trajectories of newly formed polyploids. Genome studies on such species are important for understanding the early stages of genome evolution. However, assembling neopolyploid is a challenging task due to the presence of 2 homologous (or homeologous) chromosome sets and therefore the existence of the extended paralogous regions in its genome. Post-WGD evolution of polyploids includes cytogenetic diploidization leading to the formation of species, whose polyploid origin might be hidden by disomic inheritance. Earlier we uncovered the hidden polyploid origin of the free-living flatworms of the genus Macrostomum (Macrostomum lignano, M. janickei, and M. mirumnovem). Cytogenetic diploidization in these species is accompanied by intensive chromosomal rearrangements including chromosomes fusions. In this study, we unravel the M. lignano genome organization through generation and sequencing of 2 sublines of the commonly used inbred line of M. lignano (called DV1) differing only in a copy number of the largest chromosome (MLI1). Using nontrivial assembly free comparative analysis of their genomes, we deciphered DNA sequences belonging to MLI1 and validated them by sequencing the pool of microdissected MLI1. Here we presented the uncommon mechanism of genome rediplodization of M. lignano, which consists of (i) presence of 3 subgenomes, which emerged via formation of large fused chromosomes and its variants, and (ii) sustaining their heterozygosity through inter- and intrachromosomal rearrangements.

Chromosome-level genomes of two armyworms, Mythimna separata and Mythimna loreyi, provide insights into the biosynthesis and reception of sex pheromones.

Mythimna separata and Mythimna loreyi are global pests of gramineous cereals, heavily controlled with synthetic insecticides. Here, we generated two high-quality chromosome-level genome assemblies for M. separata (688 Mb) and M. loreyi (683 Mb). Our analysis identified Z and W chromosomes, with few genes and abundant transposable elements (TEs) found on the W chromosome. We also observed a recent explosion of long interspersed nuclear elements (LINEs), which contributed to the larger genomes of Mythimna. The two armyworms diverged ~10.5 MYA, with only three chromosomes have intrachromosomal rearrangements. Additionally, we observed a tandem repeat expansion of α-amylase genes in Mythimna, which may promote the digestion of carbohydrates and exacerbate their damage to crops. Furthermore, we inferred the sex pheromone biosynthesis pathway for M. separata, M. loreyi and Spodoptera frugiperda. We discovered that M. loreyi and S. frugiperda synthesized the same major constituents of sex pheromones through different pathways. Specifically, the double bonds in the dominant sex pheromone components of S. frugiperda were generated by Δ9- and Δ11-desaturase, while they were generated by Δ11-desaturase and chain-shortening reactions in M. loreyi. We also identified pheromone receptor (PR) genes and inferred their corresponding components. These findings provide a better understanding of sex pheromone communication and promote the development of a new pest control strategy involving pheromone traps, which are more effective and environmentally friendly than current strategies.

Abstract LB251: Blocking genomic instability delays acquired resistance to MAPK inhibitor therapy in melanoma

Abstract Background: In cutaneous melanoma, the burden of chromothripsis is high prior to targeted therapy, and additional chromothripsis appears to be a key evolutionary mechanism by which cancer rapidly generates and accumulates highly dynamic structural variants (SVs). Blocking cancer genomic instability may prevent tumor escape from targeted therapies. Methods: We assembled three cohorts of tissues for WGS-based analysis of SVs. The first cohort consisted of patient-matched normal tissues, BRAF V600MUT melanoma tumors before MAPKi therapy and at disease progression (n=10 normal tissues; n=10 pretreatment tumors; n=17 acquired-resistant tumors; n=10 patients). The second cohort consisted of rapid autopsy melanoma tissues (n=3 normal tissues; n=12 acquired-resistant tumors; n=6 metastatic organ sites). The third cohort consisted of cutaneous PDX tumors. To study acquired MAPKi-resistance at the whole-genome level, we subjected PDXs (n=6 models; 1 BRAF MUT and 5 NRAS MUT models) to MAPKi therapy in NSG mice at doses sufficient to elicit tumor regression, and then generated acquired MAPKi-resistant tumors. In total, we used vehicle-treated tumors (n=6), acquired-resistant tumors (n=12), and patient-matched normal tissues (n=6) to generate WGS data. Results: Analysis of genomic amplicons due to intrachromosomal complex genomic rearrangements (CGRs) and extrachromosomal circular DNAs (ecDNAs) uncovered a significant (unpaired Student’s t-test, p=0.0002) association between acquired-resistant tumors and CGRs and/or ecDNAs harboring bona fide MAPKi-resistance genes and revealed copy number amplification of BRAF (range 4.5-27), NRAS (range 5-13), HRAS (range 13-16), MYC (range 12-15) and EGFR (CN 4.6-5), known to drive acquired MAPKi-resistance. Moreover, we validated a recurrent ecDNA by direct isolation and high-depth sequencing using a new approach referred to as CRISPR-CATCH. This alternative technique confirmed the circularized junctions of a 890 kb, driver ecDNA within this acquired-resistant clinical tumor sample. Additionally, resistance-specific (versus sensitivity-specific) chromothriptic single-base substitutions (SBSs) enriched for signatures of defects in base excision repair (BER) and in DNA mismatch repair (MMR) (Wilcoxon rank sum test, p=0.04 and p=0.005 respectively) in 14 of 31 resistant tumors (10 of 16 patients). Moreover, breakpoint-sequence analysis inferred non-homologous end-joining (NHEJ) as critical, and homologous recombination repair (HRR) as adjunctive, to DNA double-stranded break repair underlying CGR and ecDNA formation harboring MAPK-reactivation or MAPKi resistance-driver genes. Inhibition of DNA-PKc or PARP1/2, even only initially during MAPKi treatment, suppressed acquired MAPKi-resistance in melanoma cell lines (BRAF V600MUT n=3 and NRAS Q61MUT n=3) and blunted the expansion of ecDNA + CGR genomic spans. In vivo, DNA-PKi in combination with MAPKi forestalled resistance, reduced ecDNA and CGRs size, and suppressed the contribution of NHEJ in 5 out of 5 cutaneous melanoma PDXs analyzed. Conclusions: Our findings advance the concepts that preventing—instead of reversing—acquired resistant phenotypes may be more clinically impactful and that targeting DNA-PKCS and NHEJ lies at the center of this approach in stabilizing cancer genomes during oncogene-targeted therapies. Citation Format: Prashanthi Dharanipragada, Xiao Zhang, Sixue Liu, Shirley H. Lomeli, Aayoung Hong, Yan Wang, Zhentao Yang, Agustin Vega-Crespo, Antoni Ribas, Stergios J. Moschos, Gatien Moriceau, Roger S. Lo. Blocking genomic instability delays acquired resistance to MAPK inhibitor therapy in melanoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 2 (Clinical Trials and Late-Breaking Research); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(8_Suppl):Abstract nr LB251.

Abstract 6542: Detection of androgen receptor splice variants from clinical sequencing

Abstract Androgen receptor (AR) splice variants (AR-Vs) have been widely studied on its important role in prostate cancer (PC) progression. In particular, AR-Vs lacking ligand binding domains, such as AR splice variant 7 (AR-v7) that links exon 3 and cryptic exon 3 or AR variant that losses exons 5 to 7 (ARv567es), have been demonstrated as one of the resistant mechanisms to androgen deprivation therapy in PC. Despite its clinical importance, the current algorithms to detect fusions such as TopHat-Fusion have limitation to detect intra-chromosomal rearrangements including AR-Vs. In this respect, we propose a novel approach to identify AR-Vs from targeted RNA sequencing (RNA-seq) dataset. In short, our two-step approach first gathers soft-clipped or divided reads that are adjacent to the splicing sites of AR-Vs and then yields the number of splitting reads that support the existence of AR-Vs. Next, the algorithm selects the paired-end reads whose one side and the other side are mapped to the preceding and following exons, respectively. The final number of spanning reads are calculated following to the removal of low-quality reads. We validated the proposed approach using two large-scale, independent RNA-seq datasets of PC samples: 34 samples from Seoul National University Hospital (SNUH) and 558 samples from The Cancer Genome Atlas (TCGA) dataset. Overall, our approach successfully identified samples with putative AR-Vs, including AR-v7 and ARv567es. In addition, our approach successfully detected AR-v7 with reliable number of supporting reads, from the RNA-seq dataset of AR-v7-expressing cell line. Finally, statistical analyses of 34 SNUH PC patients suggested potential detection threshold in clinical sequencing settings (at least 10 split reads and 1% proportion). In the SNUH dataset, The AR-v7 positive group, which included 6 patients (17.1%) with values above the threshold, had higher AR-v7 expression levels than the negative group (p=4.8E-04). The similar pattern was also found in the analysis of TCGA dataset. Owing to its flexibility, we also confirmed that the developed approach can be used to detection of other AR-Vs, including AR-v3 and ARv567es. We expect that the further in-depth analyses including larger samples and clinical outcomes can discover clinical applicability of AR-Vs. Citation Format: Suhyun Hwangbo, Sungyoung Lee, Sheehyun Kim, Hongseok Yun. Detection of androgen receptor splice variants from clinical sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6542.

Fold-back mechanism originating inv-dup-del rearrangements in chromosomes 13 and 15.

Intrachromosomal rearrangements involve a single chromosome and can be formed by several proposed mechanisms. We reported two patients with intrachromosomal duplications and deletions, whose rearrangements and breakpoints were characterized through karyotyping, chromosomal microarray, fluorescence in situ hybridization, whole-genome sequencing, and Sanger sequencing. Inverted duplications associated with terminal deletions, known as inv-dup-del rearrangements, were found in 13q and 15q in these patients. The presence of microhomology at the junction points led to the proposal of the Fold-back mechanism for their formation. The use of different high-resolution techniques allowed for a better characterization of the rearrangements, with Sanger sequencing of the junction points being essential to infer the mechanisms of formation as it revealed microhomologies that were missed by the previous techniques. A karyotype-phenotype correlation was also performed for the characterized rearrangements.

Generating in vitro models of NTRK-fusion mesenchymal neoplasia as tools for investigating kinase oncogenic activation and response to targeted therapy

The discovery of neurotrophic tyrosine receptor kinase (NTRK) gene fusions as pan-tumor oncogenic drivers has led to new personalized therapies in oncology. Recent studies investigating NTRK fusions among mesenchymal neoplasms have identified several emerging soft tissue tumor entities displaying various phenotypes and clinical behaviors. Among them, tumors resembling lipofibromatosis or malignant peripheral nerve sheath tumors often harbor intra-chromosomal NTRK1 rearrangements, while most infantile fibrosarcomas are characterized by canonical ETV6::NTRK3 fusions. However, appropriate cellular models to investigate mechanisms of how kinase oncogenic activation through gene fusions drives such a wide spectrum of morphology and malignancy are lacking. Progress in genome editing has facilitated the efficient generation of chromosomal translocations in isogenic cell lines. In this study we employ various strategies to model NTRK fusions, including LMNA::NTRK1 (interstitial deletion) and ETV6::NTRK3 (reciprocal translocation) in human embryonic stem (hES) cells and mesenchymal progenitors (hES-MP). Here, we undertake various methods to model non-reciprocal, intrachromosomal deletions/translocations by induction of DNA double strand breaks (DSBs) exploiting either the repair mechanisms of homology directed repair (HDR) or non-homologous end joining (NHEJ). Expression of LMNA::NTRK1 or ETV6::NTRK3 fusions in either hES cells or hES-MP did not affect cell proliferation. However, the level of mRNA expression of the fusion transcripts was significantly upregulated in hES-MP, and phosphorylation of the LMNA::NTRK1 fusion oncoprotein was noted only in hES-MP but not in hES cells. Similarly, an NTRK1-driven transcriptional profile related to neuronal and neuroectodermal lineage was upregulated mainly in hES-MP, supporting the importance of appropriate cellular context in modeling cancer relevant aberrations. As proof of concept of the validity of our in vitro models, phosphorylation was depleted by two TRK inhibitors, Entrectinib and Larotrectinib, currently used as targeted therapy for tumors with NTRK fusions.

Genome Evolution in Plants: Complex Thalloid Liverworts (Marchantiopsida).

Why do some genomes stay small and simple, while others become huge, and why are some genomes more stable? In contrast to angiosperms and gymnosperms, liverworts are characterized by small genomes with low variation in size and conserved chromosome numbers. We quantified genome evolution among five Marchantiophyta (liverworts), measuring gene characteristics, transposable element (TE) landscape, collinearity, and sex chromosome evolution that might explain the small size and limited variability of liverwort genomes. No genome duplications were identified among examined liverworts and levels of duplicated genes are low. Among the liverwort species, Lunularia cruciata stands out with a genome size almost twice that of the other liverwort species investigated here, and most of this increased size is due to bursts of Ty3/Gypsy retrotransposons. Intrachromosomal rearrangements between examined liverworts are abundant but occur at a slower rate compared with angiosperms. Most genes on L. cruciata scaffolds have their orthologs on homologous Marchantia polymorpha chromosomes, indicating a low degree of rearrangements between chromosomes. Still, translocation of a fragment of the female U chromosome to an autosome was predicted from our data, which might explain the uniquely small U chromosome in L. cruciata. Low levels of gene duplication, TE activity, and chromosomal rearrangements might contribute to the apparent slow rate of morphological evolution in liverworts.

Karyotype differentiation in the Nothobranchiusugandensis species group (Teleostei, Cyprinodontiformes), seasonal fishes from the east African inland plateau, in the context of phylogeny and biogeography.

The karyotype differentiation of the twelve known members of the Nothobranchiusugandensis Wildekamp, 1994 species group is reviewed and the karyotype composition of seven of its species is described herein for the first time using a conventional cytogenetic protocol. Changes in the architecture of eukaryotic genomes often have a major impact on processes underlying reproductive isolation, adaptation and diversification. African annual killifishes of the genus Nothobranchius Peters, 1868 (Teleostei: Nothobranchiidae), which are adapted to an extreme environment of ephemeral wetland pools in African savannahs, feature extensive karyotype evolution in small, isolated populations and thus are suitable models for studying the interplay between karyotype change and species evolution. The present investigation reveals a highly conserved diploid chromosome number (2n = 36) but a variable number of chromosomal arms (46-64) among members of the N.ugandensis species group, implying a significant role of pericentric inversions and/or other types of centromeric shift in the karyotype evolution of the group. When superimposed onto a phylogenetic tree based on molecular analyses of two mitochondrial genes the cytogenetic characteristics did not show any correlation with the phylogenetic relationships within the lineage. While karyotypes of many other Nothobranchius spp. studied to date diversified mainly via chromosome fusions and fissions, the N.ugandensis species group maintains stable 2n and the karyotype differentiation seems to be constrained to intrachromosomal rearrangements. Possible reasons for this difference in the trajectory of karyotype differentiation are discussed. While genetic drift seems to be a major factor in the fixation of chromosome rearrangements in Nothobranchius, future studies are needed to assess the impact of predicted multiple inversions on the genome evolution and species diversification within the N.ugandensis species group.

Chromosome-level genome assembly of the Pacific geoduck Panopea generosa reveals major inter- and intrachromosomal rearrangements and substantial expansion of the copine gene family.

The Pacific geoduck Panopea generosa (class Bivalvia, order Adapedonta, family Hiatellidae, genus Panopea) is the largest known burrowing bivalve with considerable commercial value. Pacific geoduck and other geoduck clams play important roles in maintaining ecosystem health for their filter feeding habit and coupling pelagic and benthic processes. Here, we report a high-quality chromosome-level genome assembly of P. generosa to characterize its phylogeny and molecular mechanisms of its life strategies. The assembled P. generosa genome consists of 19 chromosomes with a size of 1.47 Gb, a contig N50 length of 1.6 Mb, and a scaffold N50 length of 73.8 Mb. The BUSCO test of the genome assembly showed 93.0% completeness. Constructed chromosome synteny revealed many occurrences of inter- and intrachromosomal rearrangements between P. generosa and Sinonovacula constricta. Of the 35,034 predicted protein-coding genes, 30,700 (87.6%) could be functionally annotated in public databases, indicating the high quality of genome annotation. Comparison of gene copy numbers of gene families among P. generosa and 11 selected species identified 507 rapidly expanded P. generosa gene families that are functionally enriched in immune and gonad development and may be involved in its complex survival strategies. In particular, genes carrying the copine domains underwent additional duplications in P. generosa, which might be important for neuronal development and immune response. The availability of a fully annotated chromosome-level genome provides a valuable dataset for genetic breeding of P. generosa.

A treatment-refractory aggressive MDS-MLD with multiple highly complex chromosome 5 intrachromosomal rearrangements: a case report

BackgroundA patient with a myelodysplastic neoplasm exhibited a karyotype with multiple complex chromosome 5 rearrangements. This patient appeared to have a catastrophic cytogenetic event that manifested as a treatment-refractory aggressive form of disease, which lead to patient demise within one year. Both the clinical presentation and disease course were unusual based on the medical history and morphologic findings. Such cases of myelodysplastic syndrome with multilineage dysplasia (MDS-MLD) with complex abnormalities are not reported in the literature.Case presentationThe patient was a 62-year-old female who presented with pancytopenia and dyspnea. The morphologic appearance of the peripheral blood smear and bone marrow biopsy, along with flow cytometric findings, favored the diagnosis of MDS-MLD unclassifiable. Myelodysplastic syndrome (MDS) with multilineage dysplasia (MDS-MLD), is an MDS characterized by one or more cytopenias and dysplastic changes in two or more of the myeloid lineages (i.e., erythroid, granulocytic, and megakaryocytic). The bone marrow, in particular, showed prominent dysplasia, including the presence of atypical megakaryocytes with small hypolobated morphology reminiscent of those typically seen in MDS with isolated 5q deletion. Cytogenetic analysis, including interphase and metaphase FISH, karyotype and SNP chromosomal microarray were performed, as well as DNA sequencing studies. Cytogenetic analysis showed a very complex karyotype featuring multiple 5q intrachromosomal rearrangements including a pericentric inversion with multiple interspersed deletions and monosomy 7. FISH studies showed a partial deletion of the PDGFRβ gene, and SNP chromosomal microarray and targeted panel-based sequencing identified biallelic loss of function of the TP53 gene. Based on the pathologic findings, the patient was treated for MDS but did not respond to either lenalidomide or azacitidine.ConclusionThe genetic changes described, in particular, the complex intrachromosomal rearrangements of chromosome 5, suggest the occurrence of a sudden catastrophic event that led to an aggressive course in the patient’s disease. Conventional karyotyping, metaphase and interphase FISH, SNP chromosomal microarray and NGS helped to identify the complex genetic changes seen in this case. This highlights the importance of utilizing a multimodality approach to fully characterize complex chromosomal events that may significantly impact disease progression, treatment and survival.

Comparative Studies of Karyotypes in the Cervidae Family

The family Cervidae is the second most diverse family in the infraorder Pecora and is characterized by a striking variability in the diploid chromosome numbers among species, ranging from 6 to 70. Chromosomal rearrangements in Cervidae have been studied in detail by chromosome painting. There are many comparative cytogenetic data for both subfamilies (Cervinae and Capreolinae) based on homologies with chromosomes of cattle and Chinese muntjac. Previously it was found that interchromosomal rearrangements are the major type of rearrangements occurring in the Cervidae family. Here, we build a detailed chromosome map of a female reindeer (Rangifer tarandus, 2n = 70, Capreolinae) and a female black muntjac (Muntiacus crinifrons, 2n = 8, Cervinae) with dromedary homologies to find out what other types of rearrangements may have underlined the variability of Cervidae karyotypes. To track chromosomal rearrangements and the distribution of nucleolus organizer regions not only during Cervidae but also Pecora evolution, we summarized new data and compared them with chromosomal maps of other already studied species. We discuss changes in the pecoran ancestral karyotype in the light of new painting data. We show that intrachromosomal rearrangements in autosomes of Cervidae are more frequent than previously thought: at least 13 inversions in evolutionary breakpoint regions were detected.

Functional Profiling Identifies Selective Enhancer Dependencies in MYC-Intact and MYC-Rearranged Diffuse Large B-Cell Lymphoma

Strong expression of MYC and MYC-regulated genes defines the molecular high-grade / double-hit signature subgroup of germinal center B-cell diffuse large B-cell lymphoma (GCB-DLBCL) which shows poor clinical outcomes with current therapy. About half of such cases show large-scale chromosomal rearrangements between the MYC locus and heterologous enhancer-bearing loci, but there is little understanding of the specific regulatory elements that activate MYC in absence of a rearrangement, or in the setting of non-immunoglobulin MYC rearrangement partners (RPs). We used high-throughput CRISPR-interference functional profiling (CRISPRi) to pinpoint cis-regulatory element dependencies in MYC-intact and MYC-rearranged lymphoma cell lines at the level of individual nucleosome-free regions (NFRs). Custom sgRNA libraries targeted 233 individual NFRs (median 42 sgRNA / NFR after stringent quality filtering) across the MYC topological domain (TAD) and MYC RP loci in six GCB-DLBCL cell lines and two comparator mantle cell lymphoma (MCL) cell lines. Promoter-targeting sgRNAs against 587 sequence-specific transcription factor (TF) genes and 65 essential gene controls allowed for interrogation of trans-factor dependencies in the same assay. CRISPRi profiling in one GCB-DLBCL and one MCL cell line with intra-chromosomal rearrangements of chromosome 8 showed strong, unique dependencies on enhancers within the cell-line-specific RPs, supporting an "enhancer-hijacking" function of these rearrangements. Three profiled cell lines showed structurally distinct variants of a recurrent t(3;8)(q27;q24) rearrangement that links the MYC gene to the distal super-enhancer region 150-240 kB upstream of BCL6. Although the distal BCL6 super-enhancers have been shown to sustain the essential BCL6 oncogene itself, we surprisingly found only weak dependency on the distal BCL6 enhancers for the three GCB-DLBCL cell lines lacking t(3;8) rearrangements. In contrast, all three GCB-DLBCL cell lines with t(3;8) rearrangements showed strong dependency on specific NFRs within the distal BCL6 super-enhancer. Validation assays showed that concurrent targeting of two BCL6 enhancer NFRs lead to even greater growth suppression and loss of MYC expression in t(3;8)-bearing cell lines. In both GCB-DLBCL cell lines lacking MYC rearrangements, we identified a strong, novel dependency on a specific enhancer ("GCB-ME1") in the 3' portion of the MYC TAD. This enhancer was acetylated in MYC-intact GCB-DLBCL patient samples and normal centroblasts, but not in MCL or MYC-rearranged DLBCL samples. We confirmed that sustained MYC expression in MYC-intact cell lines requires activity of GCB-ME1. GCB-ME1 showed little or no essentiality in the four MYC-rearranged DLBCL cell lines or in the MYC-intact MCL cell line SP-49, which is dependent on previously identified Notch/RBPJ-regulated 5' MYC enhancers (Ryan et al., 2017). ChIP-Seq showed that both GCB-ME1 and all of the most strongly essential BCL6 enhancers in t(3;8) cell lines were strongly bound by a triad of TFs, MEF2B, POU2F2, and POU2AF1, that were previously implicated in BCL6 activation. These genes were required to sustain MYC expression in both MYC-intact and t(3;8) GCB-DLBCL, and were highly essential in those cell lines, but POU2F2 and MEF2B were not essential in some cell lines with alternate MYC rearrangements. Our data show that DLBCL is highly sensitive to selective disruption of discrete MYC-activating enhancers. Transcriptional regulatory factor dependencies in DLBCL may be shaped by the requirement to sustain MYC activation via the native GCB-ME1 enhancer or alternately via enhancers in specific RP loci, a finding with implications for therapeutic targeting of transcriptional regulatory mechanisms in DLBCL. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal

EPCO-11. GATEKEEPER INACTIVATION DRIVES TUMOR PROGRESSION TO GRADE IV ASTROCYTOMA

Abstract IDH-mutant low-grade diffuse astrocytomas frequently progress to grade IV astrocytomas with implications for patient prognosis. To better understand this process, we applied whole-genome and transcriptome sequencing to matched tumor samples collected before and after progression to grade IV astrocytomas from five patients. All tumors carried an IDH1 mutation. The number of chromosomal rearrangements was increased between 1.3 and 3.5-fold in the tumors upon progression, with the exception of one case, in which the increase was only 1.03-fold. This case exhibited a hypermutation signature caused by homozygous deletion of the MSH2 gene, which encodes a member of the DNA mismatch repair complex. The most common genomic alterations acquired at progression were homozygous deletions in the CDKN2A/ RB1 -pathway or hemizygous deletion of PTEN. Additionally, PDGFRA was amplified in two grade IV tumors, with concordantly increased expression. For one of these cases, a PDGFRA-amplified subclone is likely to be present already in the low-grade astrocytoma. We further detected intrachromosomal rearrangements closeby the genes NRG3 in the progressed tumors as well as in the The Cancer Genome Atlas (TCGA) cohort. The expression of NRG3 decreased with increasing grade in the TCGA cohort and the gene was frequently deleted. Lower NRG3 expression was associated with shorter survival in the TCGA cohort. Several miRNAs showed differential expression upon progression. For two miRNAs the predicted targets were associated with cell cycle regulation and we detected inverse correlation between miRNA and target mRNA expression. While progression seems to occur via different pathways, the predicted outcome for many of the alterations was the inactivation of tumor suppressor genes and further dysregulation of cell proliferation.