Genomic analysis across 53 canine cancer types reveals novel mutations and high clinical actionability potential.

A genomic understanding of the oncogenic processes and individual variability of human cancer has steadily fueled improvement in patient outcomes over the past 20 years. Mutations within tumour tissues are routinely assessed through clinical genomic diagnostic assays by academic and commercial laboratories to facilitate diagnosis, prognosis and effective treatment stratification. The application of genomics has unveiled a wealth of mutation-based biomarkers in canine cancers, suggesting that the transformative principles that have revolutionized human cancer medicine can be brought to bear in veterinary oncology. To advance clinical genomics and genomics-guided medicine in canine oncology, we have developed and validated a canine cancer next-generation sequencing gene panel for the identification of multiple mutation types in clinical specimens. With this panel, we examined the genomic landscapes of 828 tumours from 813 dogs, spanning 53 cancer types. We identified 7856 alterations, encompassing copy number variants, single nucleotide variants, indels and internal tandem duplications. Additionally, we evaluated the clinical utility of these alterations by incorporating a biomarker framework from comprehensive curation of primary canine literature and inferences from human cancer genomic biomarker literature and clinical diagnostics. Remarkably, nearly 90% of the cases exhibited mutations with diagnostic, prognostic or therapeutic implications. Our work represents a thorough assessment of genomic landscapes in a large cohort of canine cancers, the first of its kind for its comprehensive inclusion of multiple mutation types and structured annotation of biomarkers, demonstrating the clinical potential of leveraging mutation-based biomarkers in veterinary oncology.

Genomic tumor analysis provides clinical guidance for the management of diagnostically challenging cancers in dogs.

OBJECTIVE: To evaluate the diagnostic, prognostic, and therapeutic utility of a cancer genomic diagnostic assay (SearchLight DNA; Vidium Animal Health) for diagnostically ambiguous cancer cases. ANIMALS: 69 privately owned dogs with ambiguous cancer diagnoses and for which the genomic assay was performed. PROCEDURES: Genomic assay reports generated between September 28, 2020, and July 31, 2022, for dogs with malignancy or suspected malignancy were reviewed to determine the assay's clinical utility defined as providing diagnostic clarity, prognostic information, and/or therapeutic options. RESULTS: Genomic analysis provided diagnostic clarity in 37 of 69 cases (54%; group 1) and therapeutic and/or prognostic information in 22 of the remaining 32 cases (69%; group 2) for which the diagnosis remained elusive. Overall, the genomic assay was clinically useful in 86% (59/69) of cases. CLINICAL RELEVANCE: To our knowledge, this was the first study to evaluate the multifaceted clinical utility of a single cancer genomic test in veterinary medicine. Study findings supported the use of tumor genomic testing for dogs with cancer, particularly those that are diagnostically ambiguous and therefore inherently challenging to manage. This evidence-driven genomic assay provided diagnostic guidance, prognostic support, and therapeutic options for most patients with an unclear cancer diagnosis that would otherwise have an unsubstantiated clinical plan. Furthermore, 38% (26/69) of samples were easily obtained aspirates. Sample factors (sample type, percentage of tumor cells, and number of mutations) did not influence diagnostic yield. Our study demonstrated the value of genomic testing for the management of canine cancer.

Identifying treatment options for BRAFV600 wild-type metastatic melanoma: A SU2C/MRA genomics-enabled clinical trial.

Although combination BRAF and MEK inhibitors are highly effective for the 40-50% of cutaneous metastatic melanomas harboring BRAFV600 mutations, targeted agents have been ineffective for BRAFV600wild-type (wt) metastatic melanomas. The SU2C Genomics-Enabled Medicine for Melanoma Trial utilized a Simon two-stage optimal design to assess whether comprehensive genomic profiling improves selection of molecular-based therapies for BRAFV600wt metastatic melanoma patients who had progressed on standard-of-care therapy, which may include immunotherapy. Of the response-evaluable patients, binimetinib was selected for 20 patients randomized to the genomics-enabled arm, and nine were treated on the alternate treatment arm. Response rates for 27 patients treated with targeted recommendations included one (4%) partial response, 18 (67%) with stable disease, and eight (30%) with progressive disease. Post-trial genomic and protein pathway activation mapping identified additional drug classes that may be considered for future studies. Our results highlight the complexity and heterogeneity of metastatic melanomas, as well as how the lack of response in this trial may be associated with limitations including monotherapy drug selection and the dearth of available single and combination molecularly-driven therapies to treat BRAFV600wt metastatic melanomas.

ACValidator: A novel assembly-based approach for in silico verification of circular RNAs.

Circular RNAs (circRNAs) are evolutionarily conserved RNA species that are formed when exons "back-splice" to each other. Current computational algorithms to detect these back-splicing junctions produce divergent results, and hence there is a need for a method to distinguish true-positive circRNAs. To this end, we developed Assembly based CircRNA Validator (ACValidator) for in silico verification of circRNAs. ACValidator extracts reads from a user-defined window on either side of a circRNA junction and assembles them to generate contigs. These contigs are aligned against the circRNA sequence to find contigs spanning the back-spliced junction. When evaluated on simulated datasets, ACValidator achieved over ∼80% sensitivity on datasets with an average of 10 circRNA-supporting reads and with read lengths of at least 100 bp. In experimental datasets, ACValidator produced higher verification percentages for samples treated with ribonuclease R compared to nontreated samples. Our workflow is applicable to non-polyA-selected RNAseq datasets and can also be used as a candidate selection strategy for prioritizing experimental validations. All workflow scripts are freely accessible on our GitHub page https://github.com/tgen/ACValidator along with detailed instructions to set up and run ACValidator.

Multiscale, multimodal analysis of tumor heterogeneity in IDH1 mutant vs wild-type diffuse gliomas.

Glioma is recognized to be a highly heterogeneous CNS malignancy, whose diverse cellular composition and cellular interactions have not been well characterized. To gain new clinical- and biological-insights into the genetically-bifurcated IDH1 mutant (mt) vs wildtype (wt) forms of glioma, we integrated data from protein, genomic and MR imaging from 20 treatment-naïve glioma cases and 16 recurrent GBM cases. Multiplexed immunofluorescence (MxIF) was used to generate single cell data for 43 protein markers representing all cancer hallmarks, Genomic sequencing (exome and RNA (normal and tumor) and magnetic resonance imaging (MRI) quantitative features (protocols were T1-post, FLAIR and ADC) from whole tumor, peritumoral edema and enhancing core vs equivalent normal region were also collected from patients. Based on MxIF analysis, 85,767 cells (glioma cases) and 56,304 cells (GBM cases) were used to generate cell-level data for 24 biomarkers. K-means clustering was used to generate 7 distinct groups of cells with divergent biomarker profiles and deconvolution was used to assign RNA data into three classes. Spatial and molecular heterogeneity metrics were generated for the cell data. All features were compared between IDH mt and IDHwt patients and were finally combined to provide a holistic/integrated comparison. Protein expression by hallmark was generally lower in the IDHmt vs wt patients. Molecular and spatial heterogeneity scores for angiogenesis and cell invasion also differed between IDHmt and wt gliomas irrespective of prior treatment and tumor grade; these differences also persisted in the MR imaging features of peritumoral edema and contrast enhancement volumes. A coherent picture of enhanced angiogenesis in IDHwt tumors was derived from multiple platforms (genomic, proteomic and imaging) and scales from individual proteins to cell clusters and heterogeneity, as well as bulk tumor RNA and imaging features. Longer overall survival for IDH1mt glioma patients may reflect mutation-driven alterations in cellular, molecular, and spatial heterogeneity which manifest in discernable radiological manifestations.

Leveraging Spatial Variation in Tumor Purity for Improved Somatic Variant Calling of Archival Tumor Only Samples.

Archival tumor samples represent a rich resource of annotated specimens for translational genomics research. However, standard variant calling approaches require a matched normal sample from the same individual, which is often not available in the retrospective setting, making it difficult to distinguish between true somatic variants and individual-specific germline variants. Archival sections often contain adjacent normal tissue, but this tissue can include infiltrating tumor cells. As existing comparative somatic variant callers are designed to exclude variants present in the normal sample, a novel approach is required to leverage adjacent normal tissue with infiltrating tumor cells for somatic variant calling. Here we present lumosVar 2.0, a software package designed to jointly analyze multiple samples from the same patient, built upon our previous single sample tumor only variant caller lumosVar 1.0. The approach assumes that the allelic fraction of somatic variants and germline variants follow different patterns as tumor content and copy number state change. lumosVar 2.0 estimates allele specific copy number and tumor sample fractions from the data, and uses a to model to determine expected allelic fractions for somatic and germline variants and to classify variants accordingly. To evaluate the utility of lumosVar 2.0 to jointly call somatic variants with tumor and adjacent normal samples, we used a glioblastoma dataset with matched high and low tumor content and germline whole exome sequencing data (for true somatic variants) available for each patient. Both sensitivity and positive predictive value were improved when analyzing the high tumor and low tumor samples jointly compared to analyzing the samples individually or in-silico pooling of the two samples. Finally, we applied this approach to a set of breast and prostate archival tumor samples for which tumor blocks containing adjacent normal tissue were available for sequencing. Joint analysis using lumosVar 2.0 detected several variants, including known cancer hotspot mutations that were not detected by standard somatic variant calling tools using the adjacent tissue as presumed normal reference. Together, these results demonstrate the utility of leveraging paired tissue samples to improve somatic variant calling when a constitutional sample is not available.

E6201, an intravenous MEK1 inhibitor, achieves an exceptional response in BRAF V600E-mutated metastatic malignant melanoma with brain metastases.

Malignant melanoma (MM) exhibits a high propensity for central nervous system dissemination with ~50% of metastatic MM patients developing brain metastases (BM). Targeted therapies and immune checkpoint inhibitors have improved overall survival for MM patients with BM. However, responses are usually of short duration and new agents that effectively penetrate the blood brain barrier (BBB) are needed. Here, we report a MM patient with BM who experienced an exceptional response to E6201, an ATP-competitive MEK1 inhibitor, on a Phase 1 study, with ongoing near-complete response and overall survival extending beyond 8 years. Whole exome and transcriptome sequencing revealed a high mutational burden tumor (22 mutations/Megabase) with homozygous BRAF V600E mutation. Correlative preclinical studies demonstrated broad activity for E6201 across BRAF V600E mutant melanoma cell lines and effective BBB penetration in vivo. Together, these results suggest that E6201 may represent a potential new treatment option for BRAF-mutant MM patients with BM.

Identification of therapeutic targets in chordoma through comprehensive genomic and transcriptomic analyses.

Chordoma is a rare, orphan cancer arising from embryonal precursors of bone. Surgery and radiotherapy (RT) provide excellent local control, often at the price of significant morbidity because of the structures involved and the need for relatively high doses of RT; however, recurrence remains high. Although our understanding of the genetic changes that occur in chordoma is evolving rapidly, this knowledge has yet to translate into treatments. We performed comprehensive DNA (paired tumor/normal whole-exome and shallow whole-genome) and RNA (tumor whole-transcriptome) next-generation sequencing analyses of archival sacral and clivus chordoma specimens. Incorporation of transcriptomic data enabled the identification of gene overexpression and expressed DNA alterations, thus providing additional support for potential therapeutic targets. In three patients, we identified alterations that may be amenable to off-label FDA-approved treatments for other tumor types. These alterations include FGFR1 overexpression (ponatinib, pazopanib) and copy-number duplication of CDK4 (palbociclib) and ERBB3 (gefitinib). In a third patient, germline DNA demonstrated predicted pathogenic changes in CHEK2 and ATM, which may have predisposed the patient to developing chordoma at a young age and may also be associated with potential sensitivity to PARP inhibitors because of homologous recombination repair deficiency. Last, in the fourth patient, a missense mutation in IGF1R was identified, suggesting potential activity for investigational anti-IGF1R strategies. Our findings demonstrate that chordoma patients present with aberrations in overlapping pathways. These results provide support for targeting the IGF1R/FGFR/EGFR and CDK4/6 pathways as treatment strategies for chordoma patients. This study underscores the value of comprehensive genomic and transcriptomic analysis in the development of rational, individualized treatment plans for chordoma.

Circular RNA expression and regulatory network prediction in posterior cingulate astrocytes in elderly subjects.

BACKGROUND: Circular RNAs (circRNAs) are a novel class of endogenous, non-coding RNAs that form covalently closed continuous loops and that are both highly conserved and abundant in the mammalian brain. A role for circRNAs in sponging microRNAs (miRNAs) has been proposed, but the circRNA-miRNA-mRNA interaction networks in human brain cells have not been defined. Therefore, we identified circRNAs in RNA sequencing data previously generated from astrocytes microdissected from the posterior cingulate (PC) of Alzheimer's disease (AD) patients (N = 10) and healthy elderly controls (N = 10) using four circRNA prediction algorithms - CIRI, CIRCexplorer, find_circ and KNIFE. RESULTS: Overall, utilizing these four tools, we identified a union of 4438 unique circRNAs across all samples, of which 70.3% were derived from exonic regions. Notably, the widely reported CDR1as circRNA was detected in all samples across both groups by find_circ. Given the putative miRNA regulatory function of circRNAs, we identified potential miRNA targets of circRNAs, and further, delineated circRNA-miRNA-mRNA networks using in silico methods. Pathway analysis of the genes regulated by these miRNAs identified significantly enriched immune response pathways, which is consistent with the known function of astrocytes as immune sensors in the brain. CONCLUSIONS: In this study, we performed circRNA detection on cell-specific transcriptomic data and identified potential circRNA-miRNA-mRNA regulatory networks in PC astrocytes. Given the known function of astrocytes in cerebral innate immunity and our identification of significantly enriched immune response pathways, the circRNAs we identified may be associated with such key functions. While we did not detect recurrent differentially expressed circRNAs in the context of healthy controls or AD, we report for the first time circRNAs and their potential regulatory impact in a cell-specific and region-specific manner in aged subjects. These predicted regulatory network and pathway analyses may help provide new insights into transcriptional regulation in the brain.

Whole Genome Library Construction for Next Generation Sequencing.

With the rapid evolution of genomics technologies over the past decade, whole genome sequencing (WGS) has become an increasingly accessible tool in biomedical research. WGS applications include analysis of genomic DNA from single individuals, multiple related family members, and tumor/normal samples from the same patient in the context of oncology. A number of different modalities are available for performing WGS; this chapter focuses on wet lab library construction procedures for complex short insert WGS libraries using the KAPA Hyper Prep Kit (Kapa Biosystems), and includes a discussion of appropriate quality control measures for sequencing on the Illumina HiSeq2000 platform. Additional modifications to the protocol for long insert WGS library construction, to assess structural alterations and copy number changes, are also described.

Whole Exome Library Construction for Next Generation Sequencing.

Whole exome sequencing (WES) is a DNA sequencing strategy that provides a survey of base substitutions across coding genomic locations and other regions of interest. As the coding portion of the genome encompasses only 1-2% of the entire genome, this approach represents a more cost-effective strategy to detect DNA alterations that may alter protein function, compared to whole genome sequencing. Although the research community has and is currently delineating the functional implications of sequence changes in noncoding regions of the genome, WES is a currently available assay that provides valuable information for both discovery research and precision medicine applications. In this chapter, we present a WES library preparation protocol using the KAPA Hyper Prep Kit with Agilent SureSelect Human All Exon V5+UTR probes that demonstrates high DNA-to-library conversion efficiency for sequencing on the Illumina HiSeq platform.

Integrated genomic analyses reveal frequent TERT aberrations in acral melanoma.

Genomic analyses of cutaneous melanoma (CM) have yielded biological and therapeutic insights, but understanding of non-ultraviolet (UV)-derived CMs remains limited. Deeper analysis of acral lentiginous melanoma (ALM), a rare sun-shielded melanoma subtype associated with worse survival than CM, is needed to delineate non-UV oncogenic mechanisms. We thus performed comprehensive genomic and transcriptomic analysis of 34 ALM patients. Unlike CM, somatic alterations were dominated by structural variation and absence of UV-derived mutation signatures. Only 38% of patients demonstrated driver BRAF/NRAS/NF1 mutations. In contrast with CM, we observed PAK1 copy gains in 15% of patients, and somatic TERT translocations, copy gains, and missense and promoter mutations, or germline events, in 41% of patients. We further show that in vitro TERT inhibition has cytotoxic effects on primary ALM cells. These findings provide insight into the role of TERT in ALM tumorigenesis and reveal preliminary evidence that TERT inhibition represents a potential therapeutic strategy in ALM.

Clinical Implementation of Integrated Genomic Profiling in Patients with Advanced Cancers.

DNA focused panel sequencing has been rapidly adopted to assess therapeutic targets in advanced/refractory cancer. Integrated Genomic Profiling (IGP) utilising DNA/RNA with tumour/normal comparisons in a Clinical Laboratory Improvement Amendments (CLIA) compliant setting enables a single assay to provide: therapeutic target prioritisation, novel target discovery/application and comprehensive germline assessment. A prospective study in 35 advanced/refractory cancer patients was conducted using CLIA-compliant IGP. Feasibility was assessed by estimating time to results (TTR), prioritising/assigning putative therapeutic targets, assessing drug access, ascertaining germline alterations, and assessing patient preferences/perspectives on data use/reporting. Therapeutic targets were identified using biointelligence/pathway analyses and interpreted by a Genomic Tumour Board. Seventy-five percent of cases harboured 1-3 therapeutically targetable mutations/case (median 79 mutations of potential functional significance/case). Median time to CLIA-validated results was 116 days with CLIA-validation of targets achieved in 21/22 patients. IGP directed treatment was instituted in 13 patients utilising on/off label FDA approved drugs (n = 9), clinical trials (n = 3) and single patient IND (n = 1). Preliminary clinical efficacy was noted in five patients (two partial response, three stable disease). Although barriers to broader application exist, including the need for wider availability of therapies, IGP in a CLIA-framework is feasible and valuable in selection/prioritisation of anti-cancer therapeutic targets.

A somatic reference standard for cancer genome sequencing.

Large-scale multiplexed identification of somatic alterations in cancer has become feasible with next generation sequencing (NGS). However, calibration of NGS somatic analysis tools has been hampered by a lack of tumor/normal reference standards. We thus performed paired PCR-free whole genome sequencing of a matched metastatic melanoma cell line (COLO829) and normal across three lineages and across separate institutions, with independent library preparations, sequencing, and analysis. We generated mean mapped coverages of 99X for COLO829 and 103X for the paired normal across three institutions. Results were combined with previously generated data allowing for comparison to a fourth lineage on earlier NGS technology. Aggregate variant detection led to the identification of consensus variants, including key events that represent hallmark mutation types including amplified BRAF V600E, a CDK2NA small deletion, a 12 kb PTEN deletion, and a dinucleotide TERT promoter substitution. Overall, common events include >35,000 point mutations, 446 small insertion/deletions, and >6,000 genes affected by copy number changes. We present this reference to the community as an initial standard for enabling quantitative evaluation of somatic mutation pipelines across institutions.

MKS1 regulates ciliary INPP5E levels in Joubert syndrome.

BACKGROUND: Joubert syndrome (JS) is a recessive ciliopathy characterised by a distinctive brain malformation 'the molar tooth sign'. Mutations in >27 genes cause JS, and mutations in 12 of these genes also cause Meckel-Gruber syndrome (MKS). The goals of this work are to describe the clinical features of MKS1-related JS and determine whether disease causing MKS1 mutations affect cellular phenotypes such as cilium number, length and protein content as potential mechanisms underlying JS. METHODS: We measured cilium number, length and protein content (ARL13B and INPP5E) by immunofluorescence in fibroblasts from individuals with MKS1-related JS and in a three-dimensional (3D) spheroid rescue assay to test the effects of disease-related MKS1 mutations. RESULTS: We report MKS1 mutations (eight of them previously unreported) in nine individuals with JS. A minority of the individuals with MKS1-related JS have MKS features. In contrast to the truncating mutations associated with MKS, all of the individuals with MKS1-related JS carry ≥ 1 non-truncating mutation. Fibroblasts from individuals with MKS1-related JS make normal or fewer cilia than control fibroblasts, their cilia are more variable in length than controls, and show decreased ciliary ARL13B and INPP5E. Additionally, MKS1 mutant alleles have similar effects in 3D spheroids. CONCLUSIONS: MKS1 functions in the transition zone at the base of the cilium to regulate ciliary INPP5E content, through an ARL13B-dependent mechanism. Mutations in INPP5E also cause JS, so our findings in patient fibroblasts support the notion that loss of INPP5E function, due to either mutation or mislocalisation, is a key mechanism underlying JS, downstream of MKS1 and ARL13B.

Evidence Suggesting That Discontinuous Dosing of ALK Kinase Inhibitors May Prolong Control of ALK+ Tumors.

The anaplastic lymphoma kinase (ALK) is chromosomally rearranged in a subset of certain cancers, including 2% to 7% of non-small cell lung cancers (NSCLC) and ∼70% of anaplastic large cell lymphomas (ALCL). The ALK kinase inhibitors crizotinib and ceritinib are approved for relapsed ALK(+) NSCLC, but acquired resistance to these drugs limits median progression-free survival on average to ∼10 months. Kinase domain mutations are detectable in 25% to 37% of resistant NSCLC samples, with activation of bypass signaling pathways detected frequently with or without concurrent ALK mutations. Here we report that, in contrast to NSCLC cells, drug-resistant ALCL cells show no evidence of bypassing ALK by activating alternate signaling pathways. Instead, drug resistance selected in this setting reflects upregulation of ALK itself. Notably, in the absence of crizotinib or ceritinib, we found that increased ALK signaling rapidly arrested or killed cells, allowing a prolonged control of drug-resistant tumors in vivo with the administration of discontinuous rather than continuous regimens of drug dosing. Furthermore, even when drug resistance mutations were detected in the kinase domain, overexpression of the mutant ALK was toxic to tumor cells. We confirmed these findings derived from human ALCL cells in murine pro-B cells that were transformed to cytokine independence by ectopic expression of an activated NPM-ALK fusion oncoprotein. In summary, our results show how ALK activation functions as a double-edged sword for tumor cell viability, with potential therapeutic implications.

Alzheimer's disease is associated with altered expression of genes involved in immune response and mitochondrial processes in astrocytes.

Alzheimer's disease (AD) is characterized by deficits in cerebral metabolic rates of glucose in the posterior cingulate (PC) and precuneus in AD subjects, and in APOEε4 carriers, decades before the onset of measureable cognitive deficits. However, the cellular and molecular basis of this phenotype remains to be clarified. Given the roles of astrocytes in energy storage and brain immunity, we sought to characterize the transcriptome of AD PC astrocytes. Cells were laser capture microdissected from AD (n = 10) and healthy elderly control (n = 10) subjects for RNA sequencing. We generated >5.22 billion reads and compared sequencing data between controls and AD patients. We identified differentially expressed mitochondria-related genes including TRMT61B, FASTKD2, and NDUFA4L2, and using pathway and weighted gene coexpression analyses, we identified differentially expressed immune response genes. A number of these genes, including CLU, C3, and CD74, have been implicated in beta amyloid generation or clearance. These data provide key insights into astrocyte-specific contributions to AD, and we present this data set as a publicly available resource.

Integrated genomic characterization reveals novel, therapeutically relevant drug targets in FGFR and EGFR pathways in sporadic intrahepatic cholangiocarcinoma.

Advanced cholangiocarcinoma continues to harbor a difficult prognosis and therapeutic options have been limited. During the course of a clinical trial of whole genomic sequencing seeking druggable targets, we examined six patients with advanced cholangiocarcinoma. Integrated genome-wide and whole transcriptome sequence analyses were performed on tumors from six patients with advanced, sporadic intrahepatic cholangiocarcinoma (SIC) to identify potential therapeutically actionable events. Among the somatic events captured in our analysis, we uncovered two novel therapeutically relevant genomic contexts that when acted upon, resulted in preliminary evidence of anti-tumor activity. Genome-wide structural analysis of sequence data revealed recurrent translocation events involving the FGFR2 locus in three of six assessed patients. These observations and supporting evidence triggered the use of FGFR inhibitors in these patients. In one example, preliminary anti-tumor activity of pazopanib (in vitro FGFR2 IC50≈350 nM) was noted in a patient with an FGFR2-TACC3 fusion. After progression on pazopanib, the same patient also had stable disease on ponatinib, a pan-FGFR inhibitor (in vitro, FGFR2 IC50≈8 nM). In an independent non-FGFR2 translocation patient, exome and transcriptome analysis revealed an allele specific somatic nonsense mutation (E384X) in ERRFI1, a direct negative regulator of EGFR activation. Rapid and robust disease regression was noted in this ERRFI1 inactivated tumor when treated with erlotinib, an EGFR kinase inhibitor. FGFR2 fusions and ERRFI mutations may represent novel targets in sporadic intrahepatic cholangiocarcinoma and trials should be characterized in larger cohorts of patients with these aberrations.

Joubert syndrome: brain and spinal cord malformations in genotyped cases and implications for neurodevelopmental functions of primary cilia.

Joubert syndrome (JS) is an autosomal recessive ciliopathy characterized by hypotonia, ataxia, abnormal eye movements, and intellectual disability. The brain is malformed, with severe vermian hypoplasia, fourth ventriculomegaly, and "molar tooth" appearance of the cerebral and superior cerebellar peduncles visible as consistent features on neuroimaging. Neuropathological studies, though few, suggest that several other brain and spinal cord structures, such as the dorsal cervicomedullary junction, may also be affected in at least some patients. Genetically, JS is heterogeneous, with mutations in 13 genes accounting for approximately 50% of patients. Here, we compare neuropathologic findings in five subjects with JS, including four with defined mutations in OFD1 (2 siblings), RPGRIP1L, or TCTN2. Characteristic findings in all JS genotypes included vermian hypoplasia, fragmented dentate and spinal trigeminal nuclei, hypoplastic pontine and inferior olivary nuclei, and nondecussation of corticospinal tracts. Other common findings, seen in multiple genotypes but not all subjects, were dorsal cervicomedullary heterotopia, nondecussation of superior cerebellar peduncles, enlarged arcuate nuclei, hypoplastic reticular formation, hypoplastic medial lemnisci, and dorsal spinal cord disorganization. Thus, while JS exhibits significant neuropathologic as well as genetic heterogeneity, no genotype-phenotype correlations are apparent as yet. Our findings suggest that primary cilia are important for neural patterning, progenitor proliferation, cell migration, and axon guidance in the developing human brain and spinal cord.