Large-scale Rearrangements Research Articles

BIOM-59. INDIVIDUALIZED TUMOR-SPECIFIC AMPLIFIED DNA JUNCTIONS FOR PLASMA-BASED MONITORING IN PATIENTS WITH HIGH-GRADE GLIOMAS

Abstract Accurate response assessment in patients with high-grade gliomas (HGGs) is complicated by post-treatment radiographic changes that are equivocal for treatment effect versus progression. Plasma-based liquid biopsy assays for circulating tumor DNA (ctDNA) have had limited sensitivity to date, due in part to the presence of the blood-brain barrier that limits the diffusion of tumor-derived analytes out of the central nervous system. To address this challenge, we developed a ctDNA plasma-based assay to detect individualized tumor-derived amplified junctions in patients with HGGs. These junctions occur due to large-scale chromosomal rearrangements, such as those associated with epidermal growth factor receptor (EGFR) or cyclin-dependent kinase 4 (CDK4) amplification. To identify these amplified junctions, whole genome sequencing was performed on resected tissue. Primers were designed to flank junction breakpoints and were then utilized on each patient’s plasma samples to detect the amplified junction via qPCR. Tumor-specific amplified junctions were robustly detected in plasma from 14/14 (100%) patients with IDH-mutant or wild-type grade 4 astrocytomas whose tumor genomes demonstrated amplification-related junctions, located at chromosomes 7, 11, or 12. In comparison, ctDNA could not be detected in three control cases that did not have an amplification or did not have high copy number junctions. Importantly, serial plasma samples obtained prior to and at disease recurrence in 4 patients demonstrated an increasing abundance of amplified junctions that accurately correlated with disease progression in 4/4 patients (100%). In two of these patients, increasing abundance of amplified junctions preceded radiographic changes compatible with disease progression by approximately 2 or 4 months, respectively. In conclusion, tumor-specific amplified junctions could reliably be detected in plasma from patients with IDH-mutant or wild-type grade 4 astrocytomas and correlated with changes in disease burden. Further work is currently ongoing to evaluate these promising results in a larger cohort of patients with serial plasma samples.

High frequency of mitochondrial DNA rearrangements in the peripheral blood of adults with intellectual disability.

Mitochondrial DNA (mtDNA) rearrangements are recognised factors in mitochondrial disorders and ageing, but their involvement in neurodevelopmental disorders, particularly intellectual disability (ID) and autism spectrum disorder (ASD), remains poorly understood. Previous studies have reported mitochondrial dysfunction in individuals with both ID and ASD. The aim of this study was to investigate the prevalence of large-scale mtDNA rearrangements in ID and ID with comorbid ASD (ID-ASD). We used mtDNA-targeted next-generation sequencing and the MitoSAlt high-throughput computational pipeline in peripheral blood samples from 76 patients with ID (mean age 52.5years, 37% female), 59 patients with ID-ASD (mean age 41.3years, 46% female) and 32 healthy controls (mean age 42.4years, 47% female) from Catalonia. The study revealed a high frequency of mtDNA rearrangements in patients with ID, with 10/76 (13.2%) affected individuals. However, the prevalence was significantly lower in patients with ID-ASD 1/59 (1.7%) and in HC 1/32 (3.1%). Among the mtDNA rearrangements, six were identified as deletions (median size 6937bp and median heteroplasmy level 2.3%) and six as duplications (median size 10 455bp and median heteroplasmy level 1.9%). One of the duplications, MT-ATP6 m.8765-8793dup (29bp), was present in four individuals with ID with a median heteroplasmy level of 3.9%. Our results show that mtDNA rearrangements are frequent in patients with ID, but not in those with ID-ASD, when compared to HC. Additionally, MitoSAlt has demonstrated high sensitivity and accuracy in detecting mtDNA rearrangements, even at very low heteroplasmy levels in blood samples. While the high frequency of mtDNA rearrangements in ID is noteworthy, the role of these rearrangements is currently unclear and needs to be confirmed with further data, particularly in post-mitotic tissues and through age-matched control studies.

Structure of Extracellular Vesicles and Their Effect on Bacterial Biofilms Change with the Development of Antibiotic Resistance in the Probiotic Strain Lactiplantibacillus plantarum 8p-a3.

Significant changes in lactobacillus-derived extracellular membrane vesicles, which transfer lipids, polysaccharides, proteins, and nucleic acids, were for the first time observed to accompany the development of resistance to antibiotics (amoxicillin and clarithromycin) in vitro in the probiotic strain Lactiplantibacillus plantarum 8p-a3. The changes occur together with large-scale genome rearrangements, changes in the profile of phenotypic sensitivity to antimicrobials of various groups, and evolution of virulence. Changes in vesicles affected their structure, composition, and activity against biofilms of opportunistic bacteria. The data provide for a better understanding of the molecular mechanisms of microbial survival under selective pressure of antimicrobials, the functional potential of probiotic vesicles, and probiotic safety assessments.

Mutating a flexible region of the RSV F protein can stabilize the prefusion conformation.

The respiratory syncytial virus (RSV) fusion (F) glycoprotein is highly immunogenic in its prefusion (pre-F) conformation. However, the protein is unstable, and its conformation must be stabilized for it to function effectively as an immunogen in vaccines. We present a mutagenesis strategy to arrest the RSV F protein in its pre-F state by blocking localized changes in protein structure that accompany large-scale conformational rearrangements. We generated a series of mutants and screened them in vitro to assess their potential for forming a stable pre-F. In animals, the immunogenicity of a representative mutant F protein, with a conformation confirmed by cryo-electron microscopy, elicited levels of neutralizing antibodies and protection against RSV-induced lung damage that were comparable to those of DS-Cav1, a pre-F used in a licensed vaccine.

Centromeric transposable elements and epigenetic status drive karyotypic variation in the eastern hoolock gibbon.

Great apes have maintained a stable karyotype with few large-scale rearrangements; in contrast, gibbons have undergone a high rate of chromosomal rearrangements coincident with rapid centromere turnover. Here we characterize assembled centromeres in the Eastern hoolock gibbon, Hoolock leuconedys (HLE), finding a diverse group of transposable elements (TEs) that differ from the canonical alpha satellites found across centromeres of other apes. We find that HLE centromeres contain a CpG methylation centromere dip region, providing evidence this epigenetic feature is conserved in the absence of satellite arrays; nevertheless, we report a variety of atypical centromeric features, including protein-coding genes and mismatched replication timing. Further, large structural variations define HLE centromeres and distinguish them from other gibbons. Combined with differentially methylated TEs, topologically associated domain boundaries, and segmental duplications at chromosomal breakpoints, we propose that a "perfect storm" of multiple genomic attributes with propensities for chromosome instability shaped gibbon centromere evolution.

Independent expansion, selection, and hypervariability of the TBC1D3 gene family in humans.

TBC1D3 is a primate-specific gene family that has expanded in the human lineage and has been implicated in neuronal progenitor proliferation and expansion of the frontal cortex. The gene family and its expression have been challenging to investigate because it is embedded in high-identity and highly variable segmental duplications. We sequenced and assembled the gene family using long-read sequencing data from 34 humans and 11 nonhuman primate species. Our analysis shows that this particular gene family has independently duplicated in at least five primate lineages, and the duplicated loci are enriched at sites of large-scale chromosomal rearrangements on Chromosome 17. We find that all human copy-number variation maps to two distinct clusters located at Chromosome 17q12 and that humans are highly structurally variable at this locus, differing by as many as 20 copies and ∼1 Mbp in length depending on haplotypes. We also show evidence of positive selection, as well as a significant change in the predicted human TBC1D3 protein sequence. Last, we find that, despite multiple duplications, human TBC1D3 expression is limited to a subset of copies and, most notably, from a single paralog group: TBC1D3-CDKL These observations may help explain why a gene potentially important in cortical development can be so variable in the human population.

Genomic diversity, antibiotic resistance, and virulence in South African Enterococcus faecalis and Enterococcus lactis isolates

This study presents the empirical findings of an in-depth genomic analysis of Enterococcus faecalis and Enterococcus lactis isolates from South Africa. It offers valuable insights into their genetic characteristics and their significant implications for public health. The study uncovers nuanced variations in the gene content of these isolates, despite their similar GC contents, providing a comprehensive view of the evolutionary diversity within the species. Genomic islands are identified, particularly in E. faecalis, emphasizing its propensity for horizontal gene transfer and genetic diversity, especially in terms of antibiotic resistance genes. Pangenome analysis reveals the existence of a core genome, accounting for a modest proportion of the total genes, with 2157 core genes, 1164 shell genes, and 4638 cloud genes out of 7959 genes in 52 South African E. faecalis genomes (2 from this study, 49 south Africa genomes downloaded from NCBI, and E. faecalis reference genome). Detecting large-scale genomic rearrangements, including chromosomal inversions, underscores the dynamic nature of bacterial genomes and their role in generating genetic diversity. The study uncovers an array of antibiotic resistance genes, with trimethoprim, tetracycline, glycopeptide, and multidrug resistance genes prevalent, raising concerns about the effectiveness of antibiotic treatment. Virulence gene profiling unveils a diverse repertoire of factors contributing to pathogenicity, encompassing adhesion, biofilm formation, stress resistance, and tissue damage. These empirical findings provide indispensable insights into these bacteria’s genomic dynamics, antibiotic resistance mechanisms, and virulence potential, underlining the pressing need to address antibiotic resistance and implement robust control measures.

A multiomic characterization of the leukemia cell line REH using short- and long-read sequencing.

The B-cell acute lymphoblastic leukemia (ALL) cell line REH, with the t(12;21) ETV6::RUNX1 translocation, is known to have a complex karyotype defined by a series of large-scale chromosomal rearrangements. Taken from a 15-yr-old at relapse, the cell line offers a practical model for the study of pediatric B-ALL. In recent years, short- and long-read DNA and RNA sequencing have emerged as a complement to karyotyping techniques in the resolution of structural variants in an oncological context. Here, we explore the integration of long-read PacBio and Oxford Nanopore whole-genome sequencing, IsoSeq RNA sequencing, and short-read Illumina sequencing to create a detailed genomic and transcriptomic characterization of the REH cell line. Whole-genome sequencing clarified the molecular traits of disrupted ALL-associated genes including CDKN2A, PAX5, BTG1, VPREB1, and TBL1XR1, as well as the glucocorticoid receptor NR3C1 Meanwhile, transcriptome sequencing identified seven fusion genes within the genomic breakpoints. Together, our extensive whole-genome investigation makes high-quality open-source data available to the leukemia genomics community.

Mechanisms and pathologies of human mitochondrial DNA replication and deletion formation.

Human mitochondria possess a multi-copy circular genome, mitochondrial DNA (mtDNA), that is essential for cellular energy metabolism. The number of copies of mtDNA per cell, and their integrity, are maintained by nuclear-encoded mtDNA replication and repair machineries. Aberrant mtDNA replication and mtDNA breakage are believed to cause deletions within mtDNA. The genomic location and breakpoint sequences of these deletions show similar patterns across various inherited and acquired diseases, and are also observed during normal ageing, suggesting a common mechanism of deletion formation. However, an ongoing debate over the mechanism by which mtDNA replicates has made it difficult to develop clear and testable models for how mtDNA rearrangements arise and propagate at a molecular and cellular level. These deletions may impair energy metabolism if present in a high proportion of the mtDNA copies within the cell, and can be seen in primary mitochondrial diseases, either in sporadic cases or caused by autosomal variants in nuclear-encoded mtDNA maintenance genes. These mitochondrial diseases have diverse genetic causes and multiple modes of inheritance, and show notoriously broad clinical heterogeneity with complex tissue specificities, which further makes establishing genotype-phenotype relationships challenging. In this review, we aim to cover our current understanding of how the human mitochondrial genome is replicated, the mechanisms by which mtDNA replication and repair can lead to mtDNA instability in the form of large-scale rearrangements, how rearranged mtDNAs subsequently accumulate within cells, and the pathological consequences when this occurs.

Major chromosome rearrangements in intergeneric wheat × rye hybrids in compatible and incompatible crosses detected by GBS read coverage analysis

The presence of incompatibility alleles in primary amphidiploids constitutes a reproductive barrier in newly synthesized wheat-rye hybrids. To overcome this barrier, the genome stabilization process includes large-scale chromosome rearrangements. In incompatible crosses resulting in fertile amphidiploids, the elimination of one of the incompatible alleles Eml-A1 or Eml-R1b can occur already in the somatic tissue of the wheat × rye hybrid embryo. We observed that the interaction of incompatible loci Eml-A1 of wheat and Eml-R1b of rye after overcoming embryo lethality leads to hybrid sterility in primary triticale. During subsequent seed reproductions (R1, R2 or R3) most of the chromosomes of A, B, D and R subgenomes undergo rearrangement or eliminations to increase the fertility of the amphidiploid by natural selection. Genotyping-by-sequencing (GBS) coverage analysis showed that improved fertility is associated with the elimination of entire and partial chromosomes carrying factors that either cause the disruption of plant development in hybrid plants or lead to the restoration of the euploid number of chromosomes (2n = 56) in the absence of one of the incompatible alleles. Highly fertile offspring obtained in compatible and incompatible crosses can be successfully adapted for the production of triticale pre-breeding stocks.

The chromosome-scale reference genome for the pinfish (Lagodon rhomboides) provides insights into their evolutionary and demographic history.

The pinfish (Lagodon rhomboides) is an ecologically, economically, and culturally relevant member of the family Sparidae, playing crucial roles in the marine food webs of the western Atlantic Ocean and Gulf of Mexico. Despite their high abundance and ecological importance, there is a scarcity of genomic resources for this species. We assembled and annotated a chromosome-scale genome for the pinfish, resulting in a highly contiguous 785 Mb assembly of 24 scaffolded chromosomes. The high-quality assembly contains 98.9% complete BUSCOs and shows strong synteny to other chromosome-scale genomes of fish in the family Sparidae, with a limited number of large-scale genomic rearrangements. Leveraging this new genomic resource, we found evidence of significant expansions of dietary gene families over the evolutionary history of the pinfish, which may be associated with an ontogenetic shift from carnivory to herbivory seen in this species. Estimates of historical patterns of population demography using this new reference genome identified several periods of population growth and contraction which were associated with ancient climatic shifts and sea level changes. This genome serves as a valuable reference for future studies of population genomics and differentiation and provides a much-needed genomic resource for this western Atlantic sparid.

CAM evolution is associated with gene family expansion in an explosive bromeliad radiation.

The subgenus Tillandsia (Bromeliaceae) belongs to one of the fastest radiating clades in the plant kingdom and is characterized by the repeated evolution of Crassulacean acid metabolism (CAM). Despite its complex genetic basis, this water-conserving trait has evolved independently across many plant families and is regarded as a key innovation trait and driver of ecological diversification in Bromeliaceae. By producing high-quality genome assemblies of a Tillandsia species pair displaying divergent photosynthetic phenotypes, and combining genome-wide investigations of synteny, transposable element (TE) dynamics, sequence evolution, gene family evolution, and temporal differential expression, we were able to pinpoint the genomic drivers of CAM evolution in Tillandsia. Several large-scale rearrangements associated with karyotype changes between the 2 genomes and a highly dynamic TE landscape shaped the genomes of Tillandsia. However, our analyses show that rewiring of photosynthetic metabolism is mainly obtained through regulatory evolution rather than coding sequence evolution, as CAM-related genes are differentially expressed across a 24-h cycle between the 2 species but are not candidates of positive selection. Gene orthology analyses reveal that CAM-related gene families manifesting differential expression underwent accelerated gene family expansion in the constitutive CAM species, further supporting the view of gene family evolution as a driver of CAM evolution.

The Number and Pattern of Viral Genomic Reassortments are not Necessarily Identifiable from Segment Trees.

Reassortment is an evolutionary process common in viruses with segmented genomes. These viruses can swap whole genomic segments during cellular co-infection, giving rise to novel progeny formed from the mixture of parental segments. Since large-scale genome rearrangements have the potential to generate new phenotypes, reassortment is important to both evolutionary biology and public health research. However, statistical inference of the pattern of reassortment events from phylogenetic data is exceptionally difficult, potentially involving inference of general graphs in which individual segment trees are embedded. In this paper, we argue that, in general, the number and pattern of reassortment events are not identifiable from segment trees alone, even with theoretically ideal data. We call this fact the fundamental problem of reassortment, which we illustrate using the concept of the "first-infection tree," a potentially counterfactual genealogy that would have been observed in the segment trees had no reassortment occurred. Further, we illustrate four additional problems that can arise logically in the inference of reassortment events and show, using simulated data, that these problems are not rare and can potentially distort our observation of reassortment even in small data sets. Finally, we discuss how existing methods can be augmented or adapted to account for not only the fundamental problem of reassortment, but also the four additional situations that can complicate the inference of reassortment.

A macroevolutionary role for chromosomal fusion and fission in Erebia butterflies.

The impact of large-scale chromosomal rearrangements, such as fusions and fissions, on speciation is a long-standing conundrum. We assessed whether bursts of change in chromosome numbers resulting from chromosomal fusion or fission are related to increased speciation rates in Erebia, one of the most species-rich and karyotypically variable butterfly groups. We established a genome-based phylogeny and used state-dependent birth-death models to infer trajectories of karyotype evolution. We demonstrated that rates of anagenetic chromosomal changes (i.e., along phylogenetic branches) exceed cladogenetic changes (i.e., at speciation events), but, when cladogenetic changes occur, they are mostly associated with chromosomal fissions rather than fusions. We found that the relative importance of fusion and fission differs among Erebia clades of different ages and that especially in younger, more karyotypically diverse clades, speciation is more frequently associated with cladogenetic chromosomal changes. Overall, our results imply that chromosomal fusions and fissions have contrasting macroevolutionary roles and that large-scale chromosomal rearrangements are associated with bursts of species diversification.

The rotavirus VP5*/VP8* conformational transition permeabilizes membranes to Ca2.

Rotaviruses infect cells by delivering into the cytosol a transcriptionally active inner capsid particle (a "double-layer particle": DLP). Delivery is the function of a third, outer layer, which drives uptake from the cell surface into small vesicles from which the DLPs escape. In published work, we followed stages of rhesus rotavirus (RRV) entry by live-cell imaging and correlated them with structures from cryogenic electron microscopy and tomography (cryo-EM and cryo-ET). The virus appears to wrap itself in membrane, leading to complete engulfment and loss of Ca2+ from the vesicle produced by the wrapping. One of the outer-layer proteins, VP7, is a Ca2+-stabilized trimer; loss of Ca2+ releases both VP7 and the other outer-layer protein, VP4, from the particle. VP4, activated by cleavage into VP8* and VP5*, is a trimer that undergoes a large-scale conformational rearrangement, reminiscent of the transition that viral fusion proteins undergo to penetrate a membrane. The rearrangement of VP5* thrusts a 250-residue, C-terminal segment of each of the three subunits outward, while allowing the protein to remain attached to the virus particle and to the cell being infected. We proposed that this segment inserts into the membrane of the target cell, enabling Ca2+ to cross. In the work reported here, we show the validity of key aspects of this proposed sequence. By cryo-EM studies of liposome-attached virions ("triple-layer particles": TLPs) and single-particle fluorescence imaging of liposome-attached TLPs, we confirm insertion of the VP4 C-terminal segment into the membrane and ensuing generation of a Ca2+ "leak". The results allow us to formulate a molecular description of early events in entry. We also discuss our observations in the context of other work on double-strand RNA virus entry.

Abstract 2936: Nanopore adaptive sampling detects nucleotide variants and improves large scale rearrangement characterization for diagnosis of cancer predisposition

Abstract Molecular diagnosis has become highly significant for patient management. Facing new challenges to generate more comprehensive data in reduced time and cost, we tested Nanopore (ONT) adaptive sampling to analyze a gene panel for germline predisposition to cancer. We tested 30 clinical germline samples carrying known alterations on BRCA1, BRCA2, PALB2 and MLH1 genes to detect large scale rearrangements (LSR) and single nucleotide variants (SNV). We used adaptive sampling to enrich the full sequence (from promoters to 3’UTR) of 152 cancer predisposition genes on ONT R10.4.1 MinION flowcells and the Q20 chemistry. The analysis was performed using the long-read bioinformatics pipeline developed by SeqOne Genomics. We observed a mean enrichment of 7.85 (3.91-9.85) and a mean coverage on target genes of 14.73X (4.71-35). In our cohort, 11 samples had LSR that were all detected with ONT sequencing. In addition to perfectly detect the locus of the LSR (starts and stops always in intronic regions), we discovered that a known MLPA amplification of exon 13 in the BRCA1 gene was in fact an amplification of both exons 12 and 13. This 2-exon duplication was also confirmed on a patient’s relative, originally diagnosed with an exon 13 duplication as well. Another LSR was a total deletion of the BRCA1 gene. ONT sequencing highlighted this complete deletion of BRCA1 was the consequence of a large deletion of almost 140 kbp carrying 5 different genes (BRCA1, IFI35, NBR2, RND2, and VAT1). ONT sequencing was also able to detect all class-5 SNVs (point mutations and small indel) present in 16 patient samples at low coverage (between 6 and 27X). We faced bioinformatics detection issues for 3 mutations located in homopolymer regions using standard pipelines. These were all confirmed correctly when visualizing the raw data. As we analyzed complete genes and more genes than with short read sequencing, we detected novel unknown variants. To check the reliability of those novel variant by ONT, we wondered whether we could confirm results obtained at a low coverage (10X) with average quality (Phread score >/= 15). For that we randomly selected 8 new variants we analyzed by Sanger sequencing, and all of these were confirmed, suggesting that variants detected with ONT, even at a 10X coverage and a Phread score higher than 15, could be considered as real variants. In conclusion, we showed that ONT adaptive sampling sequencing is suitable for the analysis of germline alterations. It improves the characterization of LSR and is able to detect SNV even at low coverage. Nevertheless, an improvement of bioinformatics detection of variants located in homopolymer regions should be considered. A larger comparison study between ONT adaptive sampling and short read sequencings should be performed on hundreds of samples to adapt international guidelines to third generation sequencing for germline diagnosis. Citation Format: Sandy Chevrier, Corentin Richard, Marie Mille, Denis Bertrand, Michael Blum, Nicolas Philippe, Romain Boidot. Nanopore adaptive sampling detects nucleotide variants and improves large scale rearrangement characterization for diagnosis of cancer predisposition [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 2936.

Abstract 409: Deciphering chromosomal instability in consensus molecular subtypes (CMS) in CRC: Insights from an integrative multi-omics approach

Abstract CRC is a heterogeneous disease with different molecular subtypes, which can have a profound impact on treatment response and patient outcomes. CMS classification based on gene expression profiles has been shown to differentiate patients into clinically relevant sub-groups. However, deep multi-omics characterization of the CMS classes based on high quality patient data is lacking. Through the integration of multi-omics data, we aim to enhance our understanding of CRC’s molecular heterogeneity and thereby facilitate the development of personalized therapeutic strategies. We performed a multi-omics analysis of CRC patients to identify potential new features of CMS groups. 890 fresh-frozen, surgically resected CRC tumor and adjacent normal samples were used to analyze whole genome and RNA sequencing data. DNA and RNA was prepared using Qiagen AllPrep Universal and KAPA Hyper Prep kit (DNA), as well as TruSeq Stranded Total mRNA kits and sequenced on a NovaSeq6000 system. We identified CMS subtypes based on RNA-Seq data and performed all statistical analyses using R. Our analysis revealed distinct molecular characteristics across the four CMS subtypes, including unique driver mutations, signaling cascades, and immune cell infiltration profiles. The poor prognosis associated with CMS4 patients was corroborated by our clinical data. Notably, we observed high global chromosomal instability (CIN) profiles for both CMS2 and CMS4 and identified recurrent global aneuploidy patterns within these subtypes that link alterations of known cancer genes to these large-scale structural events. We found higher levels of CIN in CMS2 compared to CMS4, which appears to be predominantly driven by numerical, rather than structural CIN. Specifically, in CMS2, SMAD2/4 is commonly impacted by deletion on chromosome arm 18q, while PLCG1 is more influenced by amplifications on 20q respectively, suggesting potential values as CRC CMS biomarkers. In conclusion, our integrative multi-omics analysis offers a comprehensive understanding of CRC molecular subtypes, which can inform personalized treatment strategies. Our findings underscore the importance of considering not only specific driver mutations but also large-scale structural rearrangements in the context of CMS subtypes. Citation Format: Julia Bischof, Jonathan Woodsmith, David Church. Deciphering chromosomal instability in consensus molecular subtypes (CMS) in CRC: Insights from an integrative multi-omics approach [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 409.

Abstract 405: Translocation detection in cancer using low-pass pore-c sequencing

Abstract Complex chromosomal rearrangements, including translocations, play a critical role in oncogenesis and are often identified as recurrent genetic aberrations in hematologic malignancies and solid tumors[1]. Translocation detection by karyotyping of an individual’s metaphase chromosomes remains challenging for detection of balanced rearrangements, which have no gain or loss of genetic material. Rearrangement detection using whole-genome sequencing is a viable alternative, but it relies on breakpoint-spanning reads and requires high sequencing depth (30-60x) with long reads to achieve high sensitivity, especially in repetitive regions of the genome[2]. Here, we describe a workflow for translocation detection at low sequencing depth with Oxford Nanopore’s long read chromatin conformation capture technique, Pore-C[3]. Importantly, for translocation calling, Pore-C does not rely on breakpoint-spanning reads but rather the high intrachromosomal interaction frequency of genomic regions around the breakpoint, thus largely lowering sequencing depth requirements and reducing mapping issues in low complexity regions. We prepared, barcoded, and sequenced Pore-C libraries of 3 cancer cell lines including lung cancer (A549), Acute Monocytic Leukemia (THP-1), melanoma (COLO 829) and its matching normal (COLO 829 BL) on a single Q20+ MinION flow cell (FLO-MIN114). Genome-wide contact maps of low-pass data (<1.5 Gbps per sample) were compared against higher depth (>30 Gbps per sample) Hi-C maps [4], revealing that low-pass Pore-C successfully captured large-scale genomic rearrangement. For example, chromosomal pairs chr8-chr11 and chr15-chr19 in A549, chr9-chr11 and chr1-chr20 in THP1, and chr7-chr15 in Colo829 were detected at <0.5X depth of sequencing coverage. These translocation events were further validated by breakpoint analysis using adaptive sampling, where targeted regions were previously confirmed by independent studies[4],[5]. Low-pass Pore-C detects translocations in an unbiased manner and does not require prior knowledge of the translocation structure. Combined with its simple sample-preparation workflow and the capability to provide genome-wide copy-number information in a single experiment, it can serve as a cost-effective and comprehensive tool for cancer genomic studies. [1] Chromosomal translocations in human cancer, Nature, 372, 143 (1994) [2] Hi-C as a tool for precise detection and characterisation of chromosomal rearrangements and copy number variation in human tumors, Genome Biology, 18, 125 (2017) [3] Identifying synergistic high-order 3D chromatin conformations from genome-scale nanopore concatemer sequencing, Nature Biotechnology, 40, 1488(2022) [4] Chromosomal translocations detection in cancer cells using chromosomal conformation capture data, Genes (Basel), 13, 7 (2022) [5] A multi-platform reference for somatic structural variation detection, Cell Genomics, 2, 6 (2022) Citation Format: Scott Hickey, Xiaoguang Dai, Sergey Aganezov, John Beaulaurier, Eoghan Harrington, Sissel Juul. Translocation detection in cancer using low-pass pore-c 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 405.

Computational Models of Claudin Assembly in Tight Junctions and Strand Properties.

Claudins are one of the major components of tight junctions (TJs) that polymerize within the cell membrane and form interactions between cells. Some claudins seal the paracellular space, limiting paracellular flux, while others form selectively permeable ion channels that control the paracellular permeability of small ions. Claudin strands are known to be dynamic and reshape within TJs to accommodate large-scale movements and rearrangements of epithelial tissues. Here, we summarize the recent computational and modeling studies on claudin assembly into tetrameric ion channels and their polymerization into μm long strands within the membrane. Computational studies ranging from all-atom molecular dynamics, coarse-grained simulations, and hybrid-resolution simulations elucidate the molecular nature of claudin assembly and function and provide a framework that describes the lateral flexibility of claudin strands.

Comparative Mitogenomics Analysis Revealed Evolutionary Divergence among Neopestalotiopsis Species Complex (Fungi: Xylariales).

The genus Neopestalotiopsis consists of obligate parasites that cause ring spot, scab, and leaf blight diseases in higher plant species. We assembled the three complete mitogenomes for the guava fruit ring spot pathogen, Neopestalotiopsis cubana. The mitogenomes are circular, with sizes of 38,666 bp, 33,846 bp, and 32,593 bp. The comparative analyses with Pestalotiopsis fici showed that N. cubana differs greatly from it in the length of the mitogenomes and the number of introns. Moreover, they showed significant differences in the gene content and tRNAs. The two genera showed little difference in gene skewness and codon preference for core protein-coding genes (PCGs). We compared gene sequencing in the mitogenomes of the order Xylariales and found large-scale gene rearrangement events, such as gene translocations and the duplication of tRNAs. N. cubana shows a unique evolutionary position in the phylum Ascomycota constructed in phylogenetic analyses. We also found a more concentrated distribution of evolutionary pressures on the PCGs of Neopestalotiopsis in the phylum Ascomycota and that they are under little selective pressure compared to other species and are subjected to purifying selection. This study explores the evolutionary dynamics of the mitogenomes of Neopestalotiopsis and provides important support for genetic and taxonomic studies.