The Canadian VirusSeq Data Portal & Duotang: open resources for SARS-CoV-2 viral sequences and genomic epidemiology.

The COVID-19 pandemic led to a large global effort to sequence SARS-CoV-2 genomes from patient samples to track viral evolution and inform public health response. Millions of SARS-CoV-2 genome sequences have been deposited in global public repositories. The Canadian COVID-19 Genomics Network (CanCOGeN - VirusSeq), a consortium tasked with coordinating expanded sequencing of SARS-CoV-2 genomes across Canada early in the pandemic, created the Canadian VirusSeq Data Portal, with associated data pipelines and procedures, to support these efforts. The goal of VirusSeq was to allow open access to Canadian SARS-CoV-2 genomic sequences and enhanced, standardized contextual data that were unavailable in other repositories and that meet FAIR standards (Findable, Accessible, Interoperable and Reusable). The Portal data submission pipeline contains data quality checking procedures and appropriate acknowledgement of data generators that encourages collaboration. Here we also highlight Duotang, a web platform that presents genomic epidemiology and modeling analyses on circulating and emerging SARS-CoV-2 variants in Canada. Duotang presents dynamic changes in variant composition of SARS-CoV-2 in Canada and by province, estimates variant growth, and displays complementary interactive visualizations, with a text overview of the current situation. The VirusSeq Data Portal and Duotang resources, alongside additional analyses and resources computed from the Portal (COVID-MVP, CoVizu), are all open-source and freely available. Together, they provide an updated picture of SARS-CoV-2 evolution to spur scientific discussions, inform public discourse, and support communication with and within public health authorities. They also serve as a framework for other jurisdictions interested in open, collaborative sequence data sharing and analyses.

Evolution of chromosome-arm aberrations in breast cancer through genetic network rewiring.

The basal breast cancer subtype is enriched for triple-negative breast cancer (TNBC) and displays consistent large chromosomal deletions. Here, we characterize evolution and maintenance of chromosome 4p (chr4p) loss in basal breast cancer. Analysis of The Cancer Genome Atlas data shows recurrent deletion of chr4p in basal breast cancer. Phylogenetic analysis of a panel of 23 primary tumor/patient-derived xenograft basal breast cancers reveals early evolution of chr4p deletion. Mechanistically we show that chr4p loss is associated with enhanced proliferation. Gene function studies identify an unknown gene, C4orf19, within chr4p, which suppresses proliferation when overexpressed-a member of the PDCD10-GCKIII kinase module we name PGCKA1. Genome-wide pooled overexpression screens using a barcoded library of human open reading frames identify chromosomal regions, including chr4p, that suppress proliferation when overexpressed in a context-dependent manner, implicating network interactions. Together, these results shed light on the early emergence of complex aneuploid karyotypes involving chr4p and adaptive landscapes shaping breast cancer genomes.

A cost-effective sequencing method for genetic studies combining high-depth whole exome and low-depth whole genome.

Whole genome sequencing (WGS) at high-depth (30X) allows the accurate discovery of variants in the coding and non-coding DNA regions and helps elucidate the genetic underpinnings of human health and diseases. Yet, due to the prohibitive cost of high-depth WGS, most large-scale genetic association studies use genotyping arrays or high-depth whole exome sequencing (WES). Here we propose a cost-effective method which we call "Whole Exome Genome Sequencing" (WEGS), that combines low-depth WGS and high-depth WES with up to 8 samples pooled and sequenced simultaneously (multiplexed). We experimentally assess the performance of WEGS with four different depth of coverage and sample multiplexing configurations. We show that the optimal WEGS configurations are 1.7-2.0 times cheaper than standard WES (no-plexing), 1.8-2.1 times cheaper than high-depth WGS, reach similar recall and precision rates in detecting coding variants as WES, and capture more population-specific variants in the rest of the genome that are difficult to recover when using genotype imputation methods. We apply WEGS to 862 patients with peripheral artery disease and show that it directly assesses more known disease-associated variants than a typical genotyping array and thousands of non-imputable variants per disease-associated locus.

Genome-wide enhancer-associated tandem repeats are expanded in cardiomyopathy.

BACKGROUND: Cardiomyopathy is a clinically and genetically heterogeneous heart condition that can lead to heart failure and sudden cardiac death in childhood. While it has a strong genetic basis, the genetic aetiology for over 50% of cardiomyopathy cases remains unknown. METHODS: In this study, we analyse the characteristics of tandem repeats from genome sequence data of unrelated individuals diagnosed with cardiomyopathy from Canada and the United Kingdom (n = 1216) and compare them to those found in the general population. We perform burden analysis to identify genomic and epigenomic features that are impacted by rare tandem repeat expansions (TREs), and enrichment analysis to identify functional pathways that are involved in the TRE-associated genes in cardiomyopathy. We use Oxford Nanopore targeted long-read sequencing to validate repeat size and methylation status of one of the most recurrent TREs. We also compare the TRE-associated genes to those that are dysregulated in the heart tissues of individuals with cardiomyopathy. FINDINGS: We demonstrate that tandem repeats that are rarely expanded in the general population are predominantly expanded in cardiomyopathy. We find that rare TREs are disproportionately present in constrained genes near transcriptional start sites, have high GC content, and frequently overlap active enhancer H3K27ac marks, where expansion-related DNA methylation may reduce gene expression. We demonstrate the gene silencing effect of expanded CGG tandem repeats in DIP2B through promoter hypermethylation. We show that the enhancer-associated loci are found in genes that are highly expressed in human cardiomyocytes and are differentially expressed in the left ventricle of the heart in individuals with cardiomyopathy. INTERPRETATION: Our findings highlight the underrecognized contribution of rare tandem repeat expansions to the risk of cardiomyopathy and suggest that rare TREs contribute to ∼4% of cardiomyopathy risk. FUNDING: Government of Ontario (RKCY), The Canadian Institutes of Health Research PJT 175329 (RKCY), The Azrieli Foundation (RKCY), SickKids Catalyst Scholar in Genetics (RKCY), The University of Toronto McLaughlin Centre (RKCY, SM), Ted Rogers Centre for Heart Research (SM), Data Sciences Institute at the University of Toronto (SM), The Canadian Institutes of Health Research PJT 175034 (SM), The Canadian Institutes of Health Research ENP 161429 under the frame of ERA PerMed (SM, RL), Heart and Stroke Foundation of Ontario & Robert M Freedom Chair in Cardiovascular Science (SM), Bitove Family Professorship of Adult Congenital Heart Disease (EO), Canada Foundation for Innovation (SWS, JR), Canada Research Chair (PS), Genome Canada (PS, JR), The Canadian Institutes of Health Research (PS).

SARS-CoV-2 rapidly evolves lineage-specific phenotypic differences when passaged repeatedly in immune-naïve mice.

The persistence of SARS-CoV-2 despite the development of vaccines and a degree of herd immunity is partly due to viral evolution reducing vaccine and treatment efficacy. Serial infections of wild-type (WT) SARS-CoV-2 in Balb/c mice yield mouse-adapted strains with greater infectivity and mortality. We investigate if passaging unmodified B.1.351 (Beta) and B.1.617.2 (Delta) 20 times in K18-ACE2 mice, expressing the human ACE2 receptor, in a BSL-3 laboratory without selective pressures, drives human health-relevant evolution and if evolution is lineage-dependent. Late-passage virus causes more severe disease, at organism and lung tissue scales, with late-passage Delta demonstrating antibody resistance and interferon suppression. This resistance co-occurs with a de novo spike S371F mutation, linked with both traits. S371F, an Omicron-characteristic mutation, is co-inherited at times with spike E1182G per Nanopore sequencing, existing in different within-sample viral variants at others. Both S371F and E1182G are linked to mammalian GOLGA7 and ZDHHC5 interactions, which mediate viral-cell entry and antiviral response. This study demonstrates SARS-CoV-2's tendency to evolve with phenotypic consequences, its evolution varying by lineage, and suggests non-dominant quasi-species contribution.

Panmixia in the American eel extends to its tropical range of distribution: Biological implications and policymaking challenges.

The American eel (Anguilla rostrata) has long been regarded as a panmictic fish and has been confirmed as such in the northern part of its range. In this paper, we tested for the first time whether panmixia extends to the tropical range of the species. To do so, we first assembled a reference genome (975 Mbp, 19 chromosomes) combining long (PacBio and Nanopore and short (Illumina paired-end) reads technologies to support both this study and future research. To test for population structure, we estimated genotype likelihoods from low-coverage whole-genome sequencing of 460 American eels, collected at 21 sampling sites (in seven geographic regions) ranging from Canada to Trinidad and Tobago. We estimated genetic distance between regions, performed ADMIXTURE-like clustering analysis and multivariate analysis, and found no evidence of population structure, thus confirming that panmixia extends to the tropical range of the species. In addition, two genomic regions with putative inversions were observed, both geographically widespread and present at similar frequencies in all regions. We discuss the implications of lack of genetic population structure for the species. Our results are key for the future genomic research in the American eel and the implementation of conservation measures throughout its geographic range. Additionally, our results can be applied to fisheries management and aquaculture of the species.

Linked-read based analysis of the medulloblastoma genome.

Introduction: Medulloblastoma is the most common type of malignant pediatric brain tumor with group 4 medulloblastomas (G4 MBs) accounting for 40% of cases. However, the molecular mechanisms that underlie this subgroup are still poorly understood. Point mutations are detected in a large number of genes at low incidence per gene while the detection of complex structural variants in recurrently affected genes typically requires the application of long-read technologies. Methods: Here, we applied linked-read sequencing, which combines the long-range genome information of long-read sequencing with the high base pair accuracy of short read sequencing and very low sample input requirements. Results: We demonstrate the detection of complex structural variants and point mutations in these tumors, and, for the first time, the detection of extrachromosomal DNA (ecDNA) with linked-reads. We provide further evidence for the high heterogeneity of somatic mutations in G4 MBs and add new complex events associated with it. Discussion: We detected several enhancer-hijacking events, an ecDNA containing the MYCN gene, and rare structural rearrangements, such a chromothripsis in a G4 medulloblastoma, chromoplexy involving 8 different chromosomes, a TERT gene rearrangement, and a PRDM6 duplication.

Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population.

BACKGROUND: Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource. RESULTS: We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue. CONCLUSIONS: By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

PTHrP Regulates Fatty Acid Metabolism via Novel lncRNA in Breast Cancer Initiation and Progression Models.

Parathyroid hormone-related peptide (PTHrP) is the primary cause of malignancy-associated hypercalcemia (MAH). We previously showed that PTHrP ablation, in the MMTV-PyMT murine model of breast cancer (BC) progression, can dramatically prolong tumor latency, slow tumor growth, and prevent metastatic spread. However, the signaling mechanisms using lineage tracing have not yet been carefully analyzed. Here, we generated Pthrpflox/flox; Cre+ mT/mG mice (KO) and Pthrpwt/wt; Cre+ mT/mG tumor mice (WT) to examine the signaling pathways under the control of PTHrP from the early to late stages of tumorigenesis. GFP+ mammary epithelial cells were further enriched for subsequent RNA sequencing (RNAseq) analyses. We observed significant upregulation of cell cycle signaling and fatty acid metabolism in PTHrP WT tumors, which are linked to tumor initiation and progression. Next, we observed that the expression levels of a novel lncRNA, GM50337, along with stearoyl-Coenzyme A desaturase 1 (Scd1) are significantly upregulated in PTHrP WT but not in KO tumors. We further validated a potential human orthologue lncRNA, OLMALINC, together with SCD1 that can be regulated via PTHrP in human BC cell lines. In conclusion, these novel findings could be used to develop targeted strategies for the treatment of BC and its metastatic complications.

Apposition of Fibroblasts With Metaplastic Gastric Cells Promotes Dysplastic Transition.

BACKGROUND & AIMS: Elements of field cancerization, including atrophic gastritis, metaplasia, and dysplasia, promote gastric cancer development in association with chronic inflammation. However, it remains unclear how stroma changes during carcinogenesis and how the stroma contributes to progression of gastric preneoplasia. Here we investigated heterogeneity of fibroblasts, one of the most important elements in the stroma, and their roles in neoplastic transformation of metaplasia. METHODS: We used single-cell transcriptomics to evaluate the cellular heterogeneity of mucosal cells from patients with gastric cancer. Tissue sections from the same cohort and tissue microarrays were used to identify the geographical distribution of distinct fibroblast subsets. We further evaluated the role of fibroblasts from pathologic mucosa in dysplastic progression of metaplastic cells using patient-derived metaplastic gastroids and fibroblasts. RESULTS: We identified 4 subsets of fibroblasts within stromal cells defined by the differential expression of PDGFRA, FBLN2, ACTA2, or PDGFRB. Each subset was distributed distinctively throughout stomach tissues with different proportions at each pathologic stage. The PDGFRα+ subset expanded in metaplasia and cancer compared with normal, maintaining a close proximity with the epithelial compartment. Co-culture of metaplasia- or cancer-derived fibroblasts with gastroids showing the characteristics of spasmolytic polypeptide-expressing metaplasia-induced disordered growth, loss of metaplastic markers, and increases in markers of dysplasia. Culture of metaplastic gastroids with conditioned media from metaplasia- or cancer-derived fibroblasts also promoted dysplastic transition. CONCLUSIONS: These findings indicate that fibroblast associations with metaplastic epithelial cells can facilitate direct transition of metaplastic spasmolytic polypeptide-expressing metaplasia cell lineages into dysplastic lineages.

Genetic analyses of DNA repair pathway associated genes implicate new candidate cancer predisposing genes in ancestrally defined ovarian cancer cases.

Not all familial ovarian cancer (OC) cases are explained by pathogenic germline variants in known risk genes. A candidate gene approach involving DNA repair pathway genes was applied to identify rare recurring pathogenic variants in familial OC cases not associated with known OC risk genes from a population exhibiting genetic drift. Whole exome sequencing (WES) data of 15 OC cases from 13 families tested negative for pathogenic variants in known OC risk genes were investigated for candidate variants in 468 DNA repair pathway genes. Filtering and prioritization criteria were applied to WES data to select top candidates for further analyses. Candidates were genotyped in ancestry defined study groups of 214 familial and 998 sporadic OC or breast cancer (BC) cases and 1025 population-matched controls and screened for additional carriers in 605 population-matched OC cases. The candidate genes were also analyzed in WES data from 937 familial or sporadic OC cases of diverse ancestries. Top candidate variants in ERCC5, EXO1, FANCC, NEIL1 and NTHL1 were identified in 5/13 (39%) OC families. Collectively, candidate variants were identified in 7/435 (1.6%) sporadic OC cases and 1/566 (0.2%) sporadic BC cases versus 1/1025 (0.1%) controls. Additional carriers were identified in 6/605 (0.9%) OC cases. Tumour DNA from ERCC5, NEIL1 and NTHL1 variant carriers exhibited loss of the wild-type allele. Carriers of various candidate variants in these genes were identified in 31/937 (3.3%) OC cases of diverse ancestries versus 0-0.004% in cancer-free controls. The strategy of applying a candidate gene approach in a population exhibiting genetic drift identified new candidate OC predisposition variants in DNA repair pathway genes.

Genomic and cytogenetic analysis of the Ceratitis capitata temperature-sensitive lethal region.

Genetic sexing strains (GSS) are an important tool in support of sterile insect technique (SIT) applications against insect pests and disease vectors. The yet unknown temperature-sensitive lethal (tsl) gene and the recently identified white pupae (wp) gene have been used as selectable markers in the most successful GSS developed so far, the Ceratitis capitata (medfly) VIENNA 8 GSS. The molecular identification of the tsl gene may open the way for its use as a marker for the development of GSS in other insect pests and disease vectors of SIT importance. Prior studies have already shown that the tsl gene is located on the right arm of chromosome 5, between the wp and Zw loci (tsl genomic region). In the present study, we used genomic, transcriptomic, bioinformatic, and cytogenetic approaches to characterize and analyze this genomic region in wild-type and tsl mutant medfly strains. Our results suggested the presence of 561 genes, with 322 of them carrying SNPs and/or insertion-deletion (indel) mutations in the tsl genomic region. Furthermore, comparative transcriptomic analysis indicated the presence of 32 differentially expressed genes, and bioinformatic analysis revealed the presence of 33 orthologs with a described heat-sensitive phenotype of Drosophila melanogaster in this region. These data can be used in functional genetic studies to identify the tsl gene(s) and the causal mutation(s) responsible for the temperature-sensitive lethal phenotype in medfly, and potentially additional genes causing a similar phenotype.

Improved platform for breast cancer circulating tumor cell enrichment and characterization with next-generation sequencing technology.

Circulating tumor cells (CTCs) represent cells shed from the primary tumor or metastatic sites and can be used to monitor treatment response and tumor recurrence. However, CTCs circulate in extremely low numbers making in-depth analysis beyond simple enumeration challenging when collected from peripheral blood. Furthermore, tumor heterogeneity, a hallmark of many tumors, especially breast cancer, further complicates CTC characterization. To overcome this limitation, we developed a platform based on the large-scale isolation of CTCs by apheresis, allowing us to collect CTCs in large numbers, which were preserved live in liquid nitrogen for further characterization. Flow cytometry followed by cell sorting (FACS) was performed using a combination of antibodies directed against cell surface markers of white blood cells (CD45) and epithelial tumor cells (CK8). Analysis of subpopulations CD45+/- and CK8+/- by bulk RNA sequencing (RNAseq) and the CD45-/CK8 positive population by single-cell RNAseq was performed. The CD45- population was enriched using CD45 magnetic beads separation and examined by IHC for pan-cytokeratin and immunofluorescence (IF) for specific markers, including the elusive circulating cancer stem cells (CSCs). CSC-rich mammospheres were grown in vitro for further analysis and treated to examine their response to chemotherapeutic agents. Finally, mammospheres were transplanted into the mammary fat pad and bone of immunodeficient mice to examine tumor growth in vivo. This platform enables the detection and collection of CTCs in early and late-stage breast cancer patients of every subtype. Markers including CD44/24, ALDH1 and CXCR4 were identified by IF and showed high expression following mammosphere culture, which responded predictably to chemotherapeutic agents. Mammospheres were also transplanted into nude mice and induced tumors in the mammary fat pad and bone following intra-tibial transplantation. Finally, bulk RNA analysis of the FACS isolated CD45+/- and CK8+/- cells showed a clear separation of CD45- away from CD45+ populations. Single-cell RNAseq of the FACS isolated CD45-/CK8+ cells showed the presence of 4-5 clusters, confirming the high degree of heterogeneity of CTCs. Our platform for large-scale isolation of CTCs using apheresis is suitable for an in-depth analysis of the cancer phenotype and may eventually allow evaluation in real-time of the disease process to optimize cancer regimens.

Molecular Genetic Characteristics of FANCI, a Proposed New Ovarian Cancer Predisposing Gene.

FANCI was recently identified as a new candidate ovarian cancer (OC)-predisposing gene from the genetic analysis of carriers of FANCI c.1813C>T; p.L605F in OC families. Here, we aimed to investigate the molecular genetic characteristics of FANCI, as they have not been described in the context of cancer. We first investigated the germline genetic landscape of two sisters with OC from the discovery FANCI c.1813C>T; p.L605F family (F1528) to re-affirm the plausibility of this candidate. As we did not find other conclusive candidates, we then performed a candidate gene approach to identify other candidate variants in genes involved in the FANCI protein interactome in OC families negative for pathogenic variants in BRCA1, BRCA2, BRIP1, RAD51C, RAD51D, and FANCI, which identified four candidate variants. We then investigated FANCI in high-grade serous ovarian carcinoma (HGSC) from FANCI c.1813C>T carriers and found evidence of loss of the wild-type allele in tumour DNA from some of these cases. The somatic genetic landscape of OC tumours from FANCI c.1813C>T carriers was investigated for mutations in selected genes, copy number alterations, and mutational signatures, which determined that the profiles of tumours from carriers were characteristic of features exhibited by HGSC cases. As other OC-predisposing genes such as BRCA1 and BRCA2 are known to increase the risk of other cancers including breast cancer, we investigated the carrier frequency of germline FANCI c.1813C>T in various cancer types and found overall more carriers among cancer cases compared to cancer-free controls (p = 0.007). In these different tumour types, we also identified a spectrum of somatic variants in FANCI that were not restricted to any specific region within the gene. Collectively, these findings expand on the characteristics described for OC cases carrying FANCI c.1813C>T; p.L605F and suggest the possible involvement of FANCI in other cancer types at the germline and/or somatic level.

Transcriptional reprogramming of skeletal muscle stem cells by the niche environment.

Adult stem cells are indispensable for tissue regeneration, but their function declines with age. The niche environment in which the stem cells reside plays a critical role in their function. However, quantification of the niche effect on stem cell function is lacking. Using muscle stem cells (MuSC) as a model, we show that aging leads to a significant transcriptomic shift in their subpopulations accompanied by locus-specific gain and loss of chromatin accessibility and DNA methylation. By combining in vivo MuSC transplantation and computational methods, we show that the expression of approximately half of all age-altered genes in MuSCs from aged male mice can be restored by exposure to a young niche environment. While there is a correlation between gene reversibility and epigenetic alterations, restoration of gene expression occurs primarily at the level of transcription. The stem cell niche environment therefore represents an important therapeutic target to enhance tissue regeneration in aging.

Single-Cell Transcriptomics Identifies Brain Endothelium Inflammatory Networks in Experimental Autoimmune Encephalomyelitis.

BACKGROUND AND OBJECTIVES: Multiple sclerosis (MS) is a neuroinflammatory and neurodegenerative disease characterized by infiltration of immune cells in multifocal areas of the CNS. The specific molecular processes allowing autoreactive immune cells to enter the CNS compartment through the blood-brain barrier remain elusive. METHODS: Using endothelial cell (EC) enrichment and single-cell RNA sequencing, we characterized the cells implicated in the neuroinflammatory processes in experimental autoimmune encephalomyelitis, an animal model of MS. Validations on human MS brain sections of the most differentially expressed genes in venous ECs were performed using immunohistochemistry and confocal microscopy. RESULTS: We found an upregulation of genes associated with antigen presentation and interferon in most populations of CNS-resident cells, including ECs. Interestingly, instead of transcriptionally distinct profiles, a continuous gradient of gene expression separated the arteriovenous zonation of the brain vasculature. However, differential gene expression analysis presented more transcriptomic alterations on the venous side of the axis, suggesting a prominent role of venous ECs in neuroinflammation. Furthermore, analysis of ligand-receptor interactions identified important potential molecular communications between venous ECs and infiltrated immune populations. To confirm the relevance of our observation in the context of human disease, we validated the protein expression of the most upregulated genes (Ackr1 and Lcn2) in MS lesions. DISCUSSION: In this study, we provide a landscape of the cellular heterogeneity associated with neuroinflammation. We also present important molecular insights for further exploration of specific cell processes that promote infiltration of immune cells inside the brain of experimental autoimmune encephalomyelitis mice.

Deep Intronic FGF14 GAA Repeat Expansion in Late-Onset Cerebellar Ataxia.

BACKGROUND: The late-onset cerebellar ataxias (LOCAs) have largely resisted molecular diagnosis. METHODS: We sequenced the genomes of six persons with autosomal dominant LOCA who were members of three French Canadian families and identified a candidate pathogenic repeat expansion. We then tested for association between the repeat expansion and disease in two independent case-control series - one French Canadian (66 patients and 209 controls) and the other German (228 patients and 199 controls). We also genotyped the repeat in 20 Australian and 31 Indian index patients. We assayed gene and protein expression in two postmortem cerebellum specimens and two induced pluripotent stem-cell (iPSC)-derived motor-neuron cell lines. RESULTS: In the six French Canadian patients, we identified a GAA repeat expansion deep in the first intron of FGF14, which encodes fibroblast growth factor 14. Cosegregation of the repeat expansion with disease in the families supported a pathogenic threshold of at least 250 GAA repeats ([GAA]≥250). There was significant association between FGF14 (GAA)≥250 expansions and LOCA in the French Canadian series (odds ratio, 105.60; 95% confidence interval [CI], 31.09 to 334.20; P<0.001) and in the German series (odds ratio, 8.76; 95% CI, 3.45 to 20.84; P<0.001). The repeat expansion was present in 61%, 18%, 15%, and 10% of French Canadian, German, Australian, and Indian index patients, respectively. In total, we identified 128 patients with LOCA who carried an FGF14 (GAA)≥250 expansion. Postmortem cerebellum specimens and iPSC-derived motor neurons from patients showed reduced expression of FGF14 RNA and protein. CONCLUSIONS: A dominantly inherited deep intronic GAA repeat expansion in FGF14 was found to be associated with LOCA. (Funded by Fondation Groupe Monaco and others.).

Detection and genomic analysis of BRAF fusions in Juvenile Pilocytic Astrocytoma through the combination and integration of multi-omic data.

BACKGROUND: Juvenile Pilocytic Astrocytomas (JPAs) are one of the most common pediatric brain tumors, and they are driven by aberrant activation of the mitogen-activated protein kinase (MAPK) signaling pathway. RAF-fusions are the most common genetic alterations identified in JPAs, with the prototypical KIAA1549-BRAF fusion leading to loss of BRAF's auto-inhibitory domain and subsequent constitutive kinase activation. JPAs are highly vascular and show pervasive immune infiltration, which can lead to low tumor cell purity in clinical samples. This can result in gene fusions that are difficult to detect with conventional omics approaches including RNA-Seq. METHODS: To this effect, we applied RNA-Seq as well as linked-read whole-genome sequencing and in situ Hi-C as new approaches to detect and characterize low-frequency gene fusions at the genomic, transcriptomic and spatial level. RESULTS: Integration of these datasets allowed the identification and detailed characterization of two novel BRAF fusion partners, PTPRZ1 and TOP2B, in addition to the canonical fusion with partner KIAA1549. Additionally, our Hi-C datasets enabled investigations of 3D genome architecture in JPAs which showed a high level of correlation in 3D compartment annotations between JPAs compared to other pediatric tumors, and high similarity to normal adult astrocytes. We detected interactions between BRAF and its fusion partners exclusively in tumor samples containing BRAF fusions. CONCLUSIONS: We demonstrate the power of integrating multi-omic datasets to identify low frequency fusions and characterize the JPA genome at high resolution. We suggest that linked-reads and Hi-C could be used in clinic for the detection and characterization of JPAs.