Search

ABCA4 is the most frequently mutated gene leading to inherited retinal disease (IRD) with over 2200 pathogenic variants reported to date. Of these, ~1% are copy number variants (CNVs) involving the deletion or duplication of genomic regions, typically >50 nucleotides in length. An in-depth assessment of the current literature based on the public database LOVD, regarding the presence of known CNVs and structural variants in ABCA4 , and additional sequencing analysis of ABCA4 using single-molecule Molecular Inversion Probes (smMIPs) for 148 probands highlighted recurrent and novel CNVs associated with ABCA4 -associated retinopathies. An analysis of the coverage depth in the sequencing data led to the identification of eleven deletions (six novel and five recurrent), three duplications (one novel and two recurrent) and one complex CNV. Of particular interest was the identification of a complex defect, i.e., a 15.3 kb duplicated segment encompassing exon 31 through intron 41 that was inserted at the junction of a downstream 2.7 kb deletion encompassing intron 44 through intron 47. In addition, we identified a 7.0 kb tandem duplication of intron 1 in three cases. The identification of CNVs in ABCA4 can provide patients and their families with a genetic diagnosis whilst expanding our understanding of the complexity of diseases caused by ABCA4 variants.

Importance: Previous studies indicated that female sex might be a modifier in Stargardt disease, which is an ABCA4-associated retinopathy., Objective: To investigate whether women are overrepresented among individuals with ABCA4-associated retinopathy who are carrying at least 1 mild allele or carrying nonmild alleles., Data Sources: Literature data, data from 2 European centers, and a new study. Data from a Radboudumc database and from the Rotterdam Eye Hospital were used for exploratory hypothesis testing., Study Selection: Studies investigating the sex ratio in individuals with ABCA4-AR and data from centers that collected ABCA4 variant and sex data. The literature search was performed on February 1, 2023; data from the centers were from before 2023., Data Extraction and Synthesis: Random-effects meta-analyses were conducted to test whether the proportions of women among individuals with ABCA4-associated retinopathy with mild and nonmild variants differed from 0.5, including subgroup analyses for mild alleles. Sensitivity analyses were performed excluding data with possibly incomplete variant identification. χ2 Tests were conducted to compare the proportions of women in adult-onset autosomal non-ABCA4-associated retinopathy and adult-onset ABCA4-associated retinopathy and to investigate if women with suspected ABCA4-associated retinopathy are more likely to obtain a genetic diagnosis. Data analyses were performed from March to October 2023., Main Outcomes and Measures: Proportion of women per ABCA4-associated retinopathy group. The exploratory testing included sex ratio comparisons for individuals with ABCA4-associated retinopathy vs those with other autosomal retinopathies and for individuals with ABCA4-associated retinopathy who underwent genetic testing vs those who did not., Results: Women were significantly overrepresented in the mild variant group (proportion, 0.59; 95% CI, 0.56-0.62; P < .001) but not in the nonmild variant group (proportion, 0.50; 95% CI, 0.46-0.54; P = .89). Sensitivity analyses confirmed these results. Subgroup analyses on mild variants showed differences in the proportions of women. Furthermore, in the Radboudumc database, the proportion of adult women among individuals with ABCA4-associated retinopathy (652/1154 = 0.56) was 0.10 (95% CI, 0.05-0.15) higher than among individuals with other retinopathies (280/602 = 0.47)., Conclusions and Relevance: This meta-analysis supports the likelihood that sex is a modifier in developing ABCA4-associated retinopathy for individuals with a mild ABCA4 allele. This finding may be relevant for prognosis predictions and recurrence risks for individuals with ABCA4-associated retinopathy. Future studies should further investigate whether the overrepresentation of women is caused by differences in the disease mechanism, by differences in health care-seeking behavior, or by health care discrimination between women and men with ABCA4-AR.

Precision medicine is rapidly gaining recognition in the field of (ultra)rare conditions, where only a few individuals in the world are affected. Clinical trial design for a small number of patients is extremely challenging, and for this reason, the development of N-of-1 strategies is explored to accelerate customized therapy design for rare cases. A strong candidate for this approach is Stargardt disease (STGD1), an autosomal recessive macular degeneration characterized by high genetic and phenotypic heterogeneity. STGD1 is caused by pathogenic variants in ABCA4 , and amongst them, several deep-intronic variants alter the pre-mRNA splicing process, generally resulting in the insertion of pseudoexons (PEs) into the final transcript. In this study, we describe a 10-year-old girl harboring the unique deep-intronic ABCA4 variant c.6817-713A>G. Clinically, she presents with typical early-onset STGD1 with a high disease symmetry between her two eyes. Molecularly, we designed antisense oligonucleotides (AONs) to block the produced PE insertion. Splicing rescue was assessed in three different in vitro models: HEK293T cells, fibroblasts, and photoreceptor precursor cells, the last two being derived from the patient. Overall, our research is intended to serve as the basis for a personalized N-of-1 AON-based treatment to stop early vision loss in this patient.

Inherited macular dystrophies (iMDs) are a group of genetic disorders, which affect the central region of the retina. To investigate the genetic basis of iMDs, we used single-molecule Molecular Inversion Probes to sequence 105 maculopathy-associated genes in 1352 patients diagnosed with iMDs. Within this cohort, 39.8% of patients were considered genetically explained by 460 different variants in 49 distinct genes of which 73 were novel variants, with some affecting splicing. The top five most frequent causative genes were ABCA4 (37.2%), PRPH2 (6.7%), CDHR1 (6.1%), PROM1 (4.3%) and RP1L1 (3.1%). Interestingly, variants with incomplete penetrance were revealed in almost one-third of patients considered solved (28.1%), and therefore, a proportion of patients may not be explained solely by the variants reported. This includes eight previously reported variants with incomplete penetrance in addition to CDHR1 :c.783G>A and CNGB3 :c.1208G>A. Notably, segregation analysis was not routinely performed for variant phasing-a limitation, which may also impact the overall diagnostic yield. The relatively high proportion of probands without any putative causal variant (60.2%) highlights the need to explore variants with incomplete penetrance, the potential modifiers of disease and the genetic overlap between iMDs and age-related macular degeneration. Our results provide valuable insights into the genetic landscape of iMDs and warrant future exploration to determine the involvement of other maculopathy genes.

Purpose: Late-onset Stargardt disease is a subtype of Stargardt disease type 1 (STGD1), defined by an age of onset of 45 years or older. We describe the disease characteristics, underlying genetics, and disease progression of late-onset STGD1 and highlight the differences from geographic atrophy., Design: Retrospective cohort study., Participants: Seventy-one patients with late-onset STGD1., Methods: Medical files were reviewed for clinical data including age at onset, initial symptoms, and best-corrected visual acuity. A quantitative and qualitative assessment of retinal pigment epithelium (RPE) atrophy was performed on fundus autofluorescence images and OCT scans., Main Outcome Measures: Age at onset, genotype, visual acuity, atrophy growth rates, and loss of external limiting membrane, ellipsoid zone, and RPE., Results: Median age at onset was 55.0 years (range, 45-82 years). A combination of a mild and severe variant in ATP-binding cassette subfamily A member 4 (ABCA4) was the most common genotype (n = 49 [69.0%]). The most frequent allele, c.5603A→T (p.Asn1868Ile), was present in 43 of 71 patients (60.6%). No combination of 2 severe variants was found. At first presentation, all patients have flecks. Foveal-sparing atrophy was present in 33.3% of eyes, whereas 21.1% had atrophy with foveal involvement. Extrafoveal atrophy was present in 38.9% of eyes, and no atrophy was evident in 6.7% of eyes. Time-to-event curves showed a median duration of 15.4 years (95% confidence interval, 11.1-19.6 years) from onset to foveal involvement. The median visual acuity decline was -0.03 Snellen decimal per year (interquartile range [IQR], -0.07 to 0.00 Snellen decimal; 0.03 logarithm of the minimum angle of resolution). Median atrophy growth was 0.590 mm 2 /year (IQR, 0.046-1.641 mm 2 /year) for definitely decreased autofluorescence and 0.650 mm 2 /year (IQR, 0.299-1.729 mm 2 /year) for total decreased autofluorescence., Conclusions: Late-onset STGD1 is a subtype of STGD1 with most commonly 1 severe and 1 mild ABCA4 variant. The general patient presents with typical fundus flecks and retinal atrophy in a foveal-sparing pattern with preserved central vision. Misdiagnosis as age-related macular degeneration should be avoided to prevent futile invasive treatments with potential complications. In addition, correct diagnosis lends patients with late-onset STGD1 the opportunity to participate in potentially beneficial therapeutic trials for STGD1., Financial Disclosure(s): The author(s) have no proprietary or commercial interest in any materials discussed in this article., (Copyright © 2023 American Academy of Ophthalmology. Published by Elsevier Inc. All rights reserved.)

The ABCA4 gene, involved in Stargardt disease, has a high percentage of splice-altering pathogenic variants, some of which cause complex RNA defects. Although antisense oligonucleotides (AONs) have shown promising results in splicing modulation, they have not yet been used to target complex splicing defects. Here, we performed AON-based rescue studies on ABCA4 complex splicing defects. Intron 13 variants c.1938-724A>G, c.1938-621G>A, c.1938-619A>G, and c.1938-514A>G all lead to the inclusion of different pseudo-exons (PEs) with and without an upstream PE (PE1). Intron 44 variant c.6148-84A>T results in multiple PE inclusions and/or exon skipping events. Five novel AONs were designed to target these defects. AON efficacy was assessed by in vitro splice assays using midigenes containing the variants of interest. All screened complex splicing defects were effectively rescued by the AONs. Although varying levels of efficacy were observed between AONs targeting the same PEs, for all variants at least one AON restored splicing to levels comparable or better than wildtype. In conclusion, AONs are a promising approach to target complex splicing defects in ABCA4 .

Missense variants in ABCA4 constitute ~50% of causal variants in Stargardt disease (STGD1). Their pathogenicity is attributed to their direct effect on protein function, whilst their potential impact on pre-mRNA splicing disruption remains poorly understood. Interestingly, synonymous ABCA4 variants have previously been classified as 'severe' variants based on in silico analyses. Here, we systemically investigated the role of synonymous and missense variants in ABCA4 splicing by combining computational predictions and experimental assays. To identify variants of interest, we used SpliceAI to ascribe defective splice predictions on a dataset of 5579 biallelic STGD1 probands. We selected those variants with predicted delta scores for acceptor/donor gain > 0.20, and no previous reports on their effect on splicing. Fifteen ABCA4 variants were selected, 4 of which were predicted to create a new splice acceptor site and 11 to create a new splice donor site. In addition, three variants of interest with delta scores < 0.20 were included. The variants were introduced in wild-type midigenes that contained 4-12 kb of ABCA4 genomic sequence, which were subsequently expressed in HEK293T cells. By using RT-PCR and Sanger sequencing, we identified splice aberrations for 16 of 18 analyzed variants. SpliceAI correctly predicted the outcomes for 15 out of 18 variants, illustrating its reliability in predicting the impact of coding ABCA4 variants on splicing. Our findings highlight a causal role for coding ABCA4 variants in splicing aberrations, improving the severity assessment of missense and synonymous ABCA4 variants, and guiding to new treatment strategies for STGD1., (© The Author(s) 2023. Published by Oxford University Press.)

The ABCA4 gene is the most frequently mutated Mendelian retinopathy-associated gene. Biallelic variants lead to a variety of phenotypes, however, for thousands of cases the underlying variants remain unknown. Here, we aim to shed further light on the missing heritability of ABCA4-associated retinopathy by analyzing a large cohort of macular dystrophy probands. A total of 858 probands were collected from 26 centers, of whom 722 carried no or one pathogenic ABCA4 variant, while 136 cases carried two ABCA4 alleles, one of which was a frequent mild variant, suggesting that deep-intronic variants (DIVs) or other cis-modifiers might have been missed. After single molecule molecular inversion probes (smMIPs)-based sequencing of the complete 128-kb ABCA4 locus, the effect of putative splice variants was assessed in vitro by midigene splice assays in HEK293T cells. The breakpoints of copy number variants (CNVs) were determined by junction PCR and Sanger sequencing. ABCA4 sequence analysis solved 207 of 520 (39.8%) naive or unsolved cases and 70 of 202 (34.7%) monoallelic cases, while additional causal variants were identified in 54 of 136 (39.7%) probands carrying two variants. Seven novel DIVs and six novel non-canonical splice site variants were detected in a total of 35 alleles and characterized, including the c.6283-321C>G variant leading to a complex splicing defect. Additionally, four novel CNVs were identified and characterized in five alleles. These results confirm that smMIPs-based sequencing of the complete ABCA4 gene provides a cost-effective method to genetically solve retinopathy cases and that several rare structural and splice altering defects remain undiscovered in Stargardt disease cases., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Introduction: Inherited retinal dystrophies (IRDs) can be caused by variants in more than 280 genes. The ATP-binding cassette transporter type A4 ( ABCA4 ) gene is one of these genes and has been linked to Stargardt disease type 1 (STGD1), fundus flavimaculatus, cone-rod dystrophy (CRD), and pan-retinal CRD. Approximately 25% of the reported ABCA4 variants affect RNA splicing. In most cases, it is necessary to perform a functional assay to determine the effect of these variants. Methods: Whole genome sequencing (WGS) was performed in one Spanish proband with Stargardt disease. The putative pathogenicity of c.6480-35A>G on splicing was investigated both in silico and in vitro . The in silico approach was based on the deep-learning tool SpliceAI. For the in vitro approach we used a midigene splice assay in HEK293T cells, based on a previously established wild-type midigene (BA29) containing ABCA4 exons 46 to 48. Results: Through the analysis of WGS data, we identified two candidate variants in ABCA4 in one proband: a previously described deletion, c.699&#95;768+342del (p.(Gln234Phefs*5)), and a novel branchpoint variant, c.6480-35A>G. Segregation analysis confirmed that the variants were in trans . For the branchpoint variant, SpliceAI predicted an acceptor gain with a high score (0.47) at position c.6480-47. A midigene splice assay in HEK293T cells revealed the inclusion of the last 47 nucleotides of intron 47 creating a premature stop codon and allowed to categorize the variant as moderately severe. Subsequent analysis revealed the presence of this variant as a second allele besides c.1958G>A p.(Arg653His) in an additional Spanish proband in a large cohort of IRD cases. Conclusion: A splice-altering effect of the branchpoint variant, confirmed by the midigene splice assay, along with the identification of this variant in a second unrelated individual affected with STGD, provides sufficient evidence to classify the variant as likely pathogenic. In addition, this research highlights the importance of studying non-coding regions and performing functional assays to provide a conclusive molecular diagnosis., Competing Interests: Author AL-L was employed by Miramoon Pharma S.L. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Rodríguez-Hidalgo, de Bruijn, Corradi, Rodenburg, Lara-López, Valverde-Megías, Ávila-Fernández, Fernandez-Caballero, Del Pozo-Valero, Corominas, Gilissen, Irigoyen, Cremers, Ayuso, Ruiz-Ederra and Roosing.)

Purpose: To determine the disease pathogenesis associated with the frequent ABCA4 variant c.5714+5G>A (p.[=,Glu1863Leufs*33])., Methods: Patient-derived photoreceptor precursor cells were generated to analyze the effect of c.5714+5G>A on splicing and perform a quantitative analysis of c.5714+5G>A products. Patients with c.5714+5G>A in trans with a null allele (i.e., c.5714+5G>A patients; n = 7) were compared with patients with two null alleles (i.e., double null patients; n = 11); with a special attention to the degree of RPE atrophy (area of definitely decreased autofluorescence and the degree of photoreceptor impairment (outer nuclear layer thickness and pattern electroretinography amplitude)., Results: RT-PCR of mRNA from patient-derived photoreceptor precursor cells showed exon 40 and exon 39/40 deletion products, as well as the normal transcript. Quantification of products showed 52.4% normal and 47.6% mutant ABCA4 mRNA. Clinically, c.5714+5G>A patients displayed significantly better structural and functional preservation of photoreceptors (thicker outer nuclear layer, presence of tubulations, higher pattern electroretinography amplitude) than double null patients with similar degrees of RPE loss, whereas double null patients exhibited signs of extensive photoreceptor ,damage even in the areas with preserved RPE., Conclusions: The prototypical STGD1 sequence of events of primary RPE and secondary photoreceptor damage is congruous with c.5714+5G>A, but not the double null genotype, which implies different and genotype-dependent disease mechanisms. We hypothesize that the relative photoreceptor sparing in c.5714+5G>A patients results from the remaining function of the ABCA4 transporter originating from the normally spliced product, possibly by decreasing the direct bisretinoid toxicity on photoreceptor membranes.

Over 15% of probands in a large cohort of more than 1500 inherited retinal degeneration patients present with a clinical diagnosis of Stargardt disease (STGD1), a recessive form of macular dystrophy caused by biallelic variants in the ABCA4 gene. Participants were clinically examined and underwent either target capture sequencing of the exons and some pathogenic intronic regions of ABCA4, sequencing of the entire ABCA4 gene or whole genome sequencing. ABCA4 c.4539 + 2028C > T, p.[= ,Arg1514Leufs*36] is a pathogenic deep intronic variant that results in a retina-specific 345-nucleotide pseudoexon inclusion. Through analysis of the Irish STGD1 cohort, 25 individuals across 18 pedigrees harbour ABCA4 c.4539 + 2028C > T and another pathogenic variant. This includes, to the best of our knowledge, the only two homozygous patients identified to date. This provides important evidence of variant pathogenicity for this deep intronic variant, highlighting the value of homozygotes for variant interpretation. 15 other heterozygous incidents of this variant in patients have been reported globally, indicating significant enrichment in the Irish population. We provide detailed genetic and clinical characterization of these patients, illustrating that ABCA4 c.4539 + 2028C > T is a variant of mild to intermediate severity. These results have important implications for unresolved STGD1 patients globally with approximately 10% of the population in some western countries claiming Irish heritage. This study exemplifies that detection and characterization of founder variants is a diagnostic imperative., (© 2023. The Author(s).)

Macular dystrophies are a group of individually rare but collectively common inherited retinal dystrophies characterised by central vision loss and loss of visual acuity. Single molecule Molecular Inversion Probes (smMIPs) have proved effective in identifying genetic variants causing macular dystrophy. Here, a previously established smMIPs panel tailored for genes associated with macular diseases has been used to examine 57 UK macular dystrophy cases, achieving a high solve rate of 63.2% (36/57). Among 27 bi-allelic STGD1 cases, only three novel ABCA4 variants were identified, illustrating that the majority of ABCA4 variants in Caucasian STGD1 cases are currently known. We examined cases with ABCA4 -associated disease in detail, comparing our results with a previously reported variant grading system, and found this model to be accurate and clinically useful. In this study, we showed that ABCA4 -associated disease could be distinguished from other forms of macular dystrophy based on clinical evaluation in the majority of cases (34/36).

Macular degenerations (MDs) are a subgroup of retinal disorders characterized by central vision loss. Knowledge is still lacking on the extent of genetic and nongenetic factors influencing inherited MD (iMD) and age-related MD (AMD) expression. Single molecule Molecular Inversion Probes (smMIPs) have proven effective in sequencing the ABCA4 gene in patients with Stargardt disease to identify associated coding and noncoding variation, however many MD patients still remain genetically unexplained. We hypothesized that the missing heritability of MDs may be revealed by smMIPs-based sequencing of all MD-associated genes and risk factors. Using 17,394 smMIPs, we sequenced the coding regions of 105 iMD and AMD-associated genes and noncoding or regulatory loci, known pseudo-exons, and the mitochondrial genome in two test cohorts that were previously screened for variants in ABCA4. Following detailed sequencing analysis of 110 probands, a diagnostic yield of 38% was observed. This established an ''MD-smMIPs panel," enabling a genotype-first approach in a high-throughput and cost-effective manner, whilst achieving uniform and high coverage across targets. Further analysis will identify known and novel variants in MD-associated genes to offer an accurate clinical diagnosis to patients. Furthermore, this will reveal new genetic associations for MD and potential genetic overlaps between iMD and AMD., (© 2022 The Authors. Human Mutation published by Wiley Periodicals LLC.)

Purpose: The effect of noncoding variants is often unknown in the absence of functional assays. Here, we characterized an ABCA4 intron 7 variant, c.859-25A>G, identified in Palestinian probands with Stargardt disease (STGD) or cone-rod dystrophy (CRD). We investigated the effect of this variant on the ABCA4 mRNA and retinal phenotype, and its prevalence in Palestine., Methods: The ABCA4 gene was sequenced completely or partially in 1998 cases with STGD or CRD. The effect of c.859-25A>G on splicing was investigated in silico using SpliceAI and in vitro using splice assays. Homozygosity mapping was performed for 16 affected individuals homozygous for c.859-25A>G. The clinical phenotype was assessed using functional and structural analyses including visual acuity, full-field electroretinography, and multimodal imaging., Results: The smMIPs-based ABCA4 sequencing revealed c.859-25A>G in 10 Palestinian probands from Hebron and Jerusalem. SpliceAI predicted a significant effect of this putative branchpoint-inactivating variant on the nearby intron 7 splice acceptor site. Splice assays revealed exon 8 skipping and two partial inclusions of intron 7, each having a deleterious effect. Additional genotyping revealed another 46 affected homozygous or compound heterozygous individuals carrying variant c.859-25A>G. Homozygotes shared a genomic segment of 59.6 to 87.9 kb and showed severe retinal defects on ophthalmoscopic evaluation., Conclusions: The ABCA4 variant c.859-25A>G disrupts a predicted branchpoint, resulting in protein truncation because of different splice defects, and is associated with early-onset STGD1 when present in homozygosity. This variant was found in 25/525 Palestinian inherited retinal dystrophy probands, representing one of the most frequent inherited retinal disease-causing variants in West-Bank Palestine.

Recurrence risk calculations in autosomal recessive diseases are complicated when the effect of genetic variants and their population frequencies and penetrances are unknown. An example of this is Stargardt disease (STGD1), a frequent recessive retinal disease caused by bi-allelic pathogenic variants in ABCA4. In this cross-sectional study, 1,619 ABCA4 variants from 5,579 individuals with STGD1 were collected and categorized by (1) severity based on statistical comparisons of their frequencies in STGD1-affected individuals versus the general population, (2) their observed versus expected homozygous occurrence in STGD1-affected individuals, (3) their occurrence in combination with established mild alleles in STGD1-affected individuals, and (4) previous functional and clinical studies. We used the sum allele frequencies of these severity categories to estimate recurrence risks for offspring of STGD1-affected individuals and carriers of pathogenic ABCA4 variants. The risk for offspring of an STGD1-affected individual with the "severe|severe" genotype or a "severe|mild with complete penetrance" genotype to develop STGD1 at some moment in life was estimated at 2.8%-3.1% (1 in 36-32 individuals) and 1.6%-1.8% (1 in 62-57 individuals), respectively. The risk to develop STGD1 in childhood was estimated to be 2- to 4-fold lower: 0.68%-0.79% (1 in 148-126) and 0.34%-0.39% (1 in 296-252), respectively. In conclusion, we established personalized recurrence risk calculations for STGD1-affected individuals with different combinations of variants. We thus propose an expanded genotype-based personalized counseling to appreciate the variable recurrence risks for STGD1-affected individuals. This represents a conceptual breakthrough because risk calculations for STGD1 may be exemplary for many other inherited diseases., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Inherited retinal diseases (IRDs) are a major cause of visual impairment. These clinically heterogeneous disorders are caused by pathogenic variants in more than 270 genes. As 30-40% of cases remain genetically unexplained following conventional genetic testing, we aimed to obtain a genetic diagnosis in an IRD cohort in which the genetic cause was not found using whole-exome sequencing or targeted capture sequencing. We performed whole-genome sequencing (WGS) to identify causative variants in 100 unresolved cases. After initial prioritization, we performed an in-depth interrogation of all noncoding and structural variants in genes when one candidate variant was detected. In addition, functional analysis of putative splice-altering variants was performed using in vitro splice assays. We identified the genetic cause of the disease in 24 patients. Causative coding variants were observed in genes such as ATXN7, CEP78, EYS, FAM161A, and HGSNAT. Gene disrupting structural variants were also detected in ATXN7, PRPF31, and RPGRIP1. In 14 monoallelic cases, we prioritized candidate noncanonical splice sites or deep-intronic variants that were predicted to disrupt the splicing process based on in silico analyses. Of these, seven cases were resolved as they carried pathogenic splice defects. WGS is a powerful tool to identify causative variants residing outside coding regions or heterozygous structural variants. This approach was most efficient in cases with a distinct clinical diagnosis. In addition, in vitro splice assays provide important evidence of the pathogenicity of rare variants., (© 2021. The Author(s).)

Autism spectrum disorders (ASD) are characterized by a high degree of genetic heterogeneity. Genomic studies identified common pathological processes underlying the heterogeneous clinical manifestations of ASD, and transcriptome analyses revealed that gene networks involved in synapse development, neuronal activity, and immune function are deregulated in ASD. Mouse models provide unique tools to investigate the neurobiological basis of ASD; however, a comprehensive approach to identify transcriptional abnormalities in different ASD models has never been performed. Here we used two well-recognized ASD mouse models, BTBR T(+) Itpr3 (tf) /J (BTBR) and Engrailed-2 knockout (En2 (-/-)), to identify conserved ASD-related molecular signatures. En2 (-/-) mice bear a mutation within the EN2 transcription factor homeobox, while BTBR is an inbred strain with unknown genetic defects. Hippocampal RNA samples from BTBR, En2 (-/-) and respective control (C57Bl/6J and En2 (+/+)) adult mice were assessed for differential gene expression using microarrays. A total of 153 genes were similarly deregulated in the BTBR and En2 (-/-) hippocampus. Mouse phenotype and gene ontology enrichment analyses were performed on BTBR and En2 (-/-) hippocampal differentially expressed genes (DEGs). Pathways represented in both BTBR and En2 (-/-) hippocampal DEGs included abnormal behavioral response and chemokine/MAP kinase signaling. Genes involved in abnormal function of the immune system and abnormal synaptic transmission/seizures were significantly represented among BTBR and En2 (-/-) DEGs, respectively. Interestingly, both BTBR and En2 (-/-) hippocampal DEGs showed a significant enrichment of ASD and schizophrenia (SCZ)-associated genes. Specific gene sets were enriched in the two models: microglial genes were significantly enriched among BTBR DEGs, whereas GABAergic/glutamatergic postsynaptic genes, FMRP-interacting genes and epilepsy-related genes were significantly enriched among En2 (-/-) DEGs. Weighted correlation network analysis (WGCNA) performed on BTBR and En2 (-/-) hippocampal transcriptomes together identified six modules significantly enriched in ASD-related genes. Each of these modules showed a specific enrichment profile in neuronal and glial genes, as well as in genes associated to ASD comorbidities such as epilepsy and SCZ. Our data reveal significant transcriptional similarities and differences between the BTBR and En2 (-/-) hippocampus, indicating that transcriptome analysis of ASD mouse models may contribute to identify novel molecular targets for pharmacological studies.